Thermo-Osmosis in Charged Nanochannels: Effects of Surface Charge and Ionic Strength.

Thermo-osmosis refers to fluid migration due to the temperature gradient. The mechanistic understanding of thermo-osmosis in charged nano-porous media is still incomplete, while it is important for several environmental and energy applications, such as low-grade waste heat recovery, wastewater recovery, fuel cells, and nuclear waste storage. This paper presents results from a series of molecular dynamics simulations of thermo-osmosis in charged silica nanochannels that advance the understanding of the phenomenon. Simulations with pure water and water with dissolved NaCl are considered. First, the effect of surface charge on the sign and magnitude of the thermo-osmotic coefficient is quantified. This effect was found to be mainly linked to the structural modifications of an aqueous electrical double layer (EDL) caused by the nanoconfinement and surface charges. In addition, the results illustrate that the surface charges reduce the self-diffusivity and thermo-osmosis of interfacial liquid. The thermo-osmosis was found to change direction when the surface charge density exceeds -0.03C · m-2. It was found that the thermo-osmotic flow and self-diffusivity increase with the concentration of NaCl. The fluxes of solvent and solute are decoupled by considering the Ludwig-Soret effect of NaCl ions to identify the main mechanisms controlling the behavior. In addition to the advance in microscopic quantification and mechanistic understanding of thermo-osmosis, the work provides approaches to investigate a broader category of coupled heat and mass transfer problems in nanoscale space.

Thermal diffusion of ionic species in charged nanochannels.

Diffusion of ions due to temperature gradients (known as thermal diffusion) in charged nanochannels is of interest in several engineering fields, including energy recovery and environmental protection. This paper presents a fundamental investigation of the thermal diffusion of sodium chloride in charged silica nanochannels performed by molecular dynamics (MD). The results reveal the effects of nanoconfinement and surface charges on the sign and magnitude of the Soret coefficient. It is shown that the sign and magnitude of the Soret coefficient are controlled by the structural modifications of the interfacial solutions. These modifications include the ionic solvation and hydrogen bond structure induced by the nanoconfinement and surface charges. The results show that both nanoconfinement and surface charges can make the solutions more thermophilic. Furthermore, the thermal diffusion of solutions in boundary layers is significantly different from that of solutions in bulk fluid, contributing to the overall difference between the thermal diffusivity of pore fluid and that associated with bulk fluid. The findings provide further understanding of thermal diffusion in nano-porous systems. The proposed MD simulation methodology is applicable to a wider category of coupled heat and mass transfer problems in nanoscale spaces.

Image-Based vs. Parametric Modelling of Concrete Meso-Structures.

Damage initiation and crack propagation in concrete are associated with localisation of energy dissipation by the concrete meso-structure. Meso-scale models are, therefore, required for realistic analysis of concrete non-linear behaviour. Such models are constructed either from X-ray Computed Tomography images (image-based modelling) or by in silico meso-structure generation (parametric modelling), while both approaches are widely used and their advantages and disadvantages are recognised, little work is done on comparing their performance in predicting measured macroscopic behaviour with equivalent constitutive relations for meso-structural features. This work uses microstructure characterisation and mechanical behaviour data to construct, validate and compare the two modelling approaches. The macroscopic behaviour obtained with both meso-structural models is found to be in good agreement with experimental data. Differences are observed only between the predicted distributions of damage within specimens. These outcomes suggest that the computationally simpler parametric meso-structures are sufficient to derive stress-strain behaviour for engineering-scale models in the absence of other environmental factors. The observed differences in damage distribution could be important for analysis of coupled behaviour, e.g., mass transport and chemical reactions affecting local mechanical properties and being affected by local damage. Establishing the importance of damage distribution is such cases requires further research.

Incorporation of Obstacle Hardening into Local Approach to Cleavage Fracture to Predict Temperature Effects in the Ductile to Brittle Transition Regime.

Ductile-to-brittle-transition refers to observable change in fracture mode with decreasing temperature-from slow ductile crack growth to rapid cleavage. It is exhibited by body-centred cubic metals and presents a challenge for integrity assessment of structural components made of such metals. Local approaches to cleavage fracture, based on Weibull stress as a cleavage crack-driving force, have been shown to predict fracture toughness at very low temperatures. However, they are ineffective in the transition regime without the recalibration of Weibull stress parameters, which requires further testing and thus diminishes their predictive capability. We propose new Weibull stress formulation with thinning function based on obstacle hardening model, which modifies the number of cleavage-initiating features with temperature. Our model is implemented as a post-processor of finite element analysis results. It is applied to analyses of standard compact tension specimens of typical reactor pressure vessel steel, for which deformation and fracture toughness properties in the transition regime are available. It is shown that the new Weibull stress is independent of temperature, and of Weibull shape parameter, within the experimental error. It accurately predicts the fracture toughness at any temperature in the transition regime without relying upon empirical fits for the first time.

Thermo-osmosis in hydrophilic nanochannels: mechanism and size effect.

Understanding thermo-osmosis in nanoscale channels and pores is essential for both theoretical advances of thermally induced mass flow and a wide range of emerging industrial applications. We present a new mechanistic understanding and quantification of thermo-osmosis at nanometric/sub-nanometric length scales and link the outcomes with the non-equilibrium thermodynamics of the phenomenon. The work is focused on thermo-osmosis of water in quartz slit nanochannels, which is analysed by molecular dynamics (MD) simulations of mechano-caloric and thermo-osmotic systems. We investigate the applicability of Onsager reciprocal relation, irreversible thermodynamics, and continuum fluid mechanics at the nanoscale. Further, we analyse the effects of channel size on the thermo-osmosis coefficient, and show, for the first time, that these arise from specific liquid structures dictated by the channel size. The mechanical conditions of the interfacial water under different temperatures are quantified using a continuum approach (pressure tensor distribution) and a discrete approach (body force per molecule) to elucidate the underlying mechanism of thermo-osmosis. The results show that the fluid molecules located in the boundary layers adjacent to the solid surfaces experience a driving force which generates the thermo-osmotic flow. While the findings provide a fundamental understanding of thermo-osmosis, the methods developed provide a route for analysis of the entire class of coupled heat and mass transport phenomena in nanoscale structures.

Effect of colloids on non-Fickian transport of strontium in sediments elucidated by continuous-time random walk analysis.

Understanding the influence of colloids on radionuclide migration is of significance to evaluate environmental risks for radioactive waste disposals. In order to formulate an appropriate modelling framework that can quantify and interpret the anomalous transport of Strontium (Sr) in the absence and presence of colloids, the continuous time random walk (CTRW) approach is implemented in this work using available experimental information. The results show that the transport of Sr and its recovery are enhanced in the presence of colloids. The causes can be largely attributed to the trap-release processes, e.g. electrostatic interactions of Sr, colloids and natural sediments, and differences in pore structures, which gave rise to the varying interstitial velocities of dissolved and, if any, colloid-associated Sr. Good agreement between the CTRW simulations and the column-scale observations is demonstrated. Regardless of the presence of colloids, the CTRW modelling captures the characteristics of non-Fickian anomalous transport (0 < ß < 2) of Sr. In particular, a range of 0 < ß < 1, corresponding to the cases with greater recoveries, reveal strongly non-Fickian transport with distinctive earlier arrivals and tailing effects, likely due to the physicochemical heterogeneities, i.e. the repulsive interactions and/or the macro-pores originating from local heterogeneities. The results imply that colloids can increase the Sr transport as a barrier of Sr sorption onto sediments herein, apart from often being carriers of sored radionuclides in aqueous phase. From a modelling perspective, the findings show that the established CTRW model is valid for quantifying the non-Fickian and promoted transport of Sr with colloids.

Structural modelling of the cardiovascular system.

Computational modelling of the cardiovascular system offers much promise, but represents a truly interdisciplinary challenge, requiring knowledge of physiology, mechanics of materials, fluid dynamics and biochemistry. This paper aims to provide a summary of the recent advances in cardiovascular structural modelling, including the numerical methods, main constitutive models and modelling procedures developed to represent cardiovascular structures and pathologies across a broad range of length and timescales; serving as an accessible point of reference to newcomers to the field. The class of so-called hyperelastic materials provides the theoretical foundation for the modelling of how these materials deform under load, and so an overview of these models is provided; comparing classical to application-specific phenomenological models. The physiology is split into components and pathologies of the cardiovascular system and linked back to constitutive modelling developments, identifying current state of the art in modelling procedures from both clinical and engineering sources. Models which have originally been derived for one application and scale are shown to be used for an increasing range and for similar applications. The trend for such approaches is discussed in the context of increasing availability of high performance computing resources, where in some cases computer hardware can impact the choice of modelling approach used.

Review of pore network modelling of porous media: Experimental characterisations, network constructions and applications to reactive transport.

Pore network models have been applied widely for simulating a variety of different physical and chemical processes, including phase exchange, non-Newtonian displacement, non-Darcy flow, reactive transport and thermodynamically consistent oil layers. The realism of such modelling, i.e. the credibility of their predictions, depends to a large extent on the quality of the correspondence between the pore space of a given medium and the pore network constructed as its representation. The main experimental techniques for pore space characterisation, including direct imaging, mercury intrusion porosimetry and gas adsorption, are firstly summarised. A review of the main pore network construction techniques is then presented. Particular focus is given on how such constructions are adapted to the data from experimentally characterised pore systems. Current applications of pore network models are considered, with special emphasis on the effects of adsorption, dissolution and precipitation, as well as biomass growth, on transport coefficients. Pore network models are found to be a valuable tool for understanding and predicting meso-scale phenomena, linking single pore processes, where other techniques are more accurate, and the homogenised continuum porous media, used by engineering community.

Measurement and modelling of reactive transport in geological barriers for nuclear waste containment.

Compacted clays are considered as excellent candidates for barriers to radionuclide transport in future repositories for nuclear waste due to their very low hydraulic permeability. Diffusion is the dominant transport mechanism, controlled by a nano-scale pore system. Assessment of the clays' long-term containment function requires adequate modelling of such pore systems and their evolution. Existing characterisation techniques do not provide complete pore space information for effective modelling, such as pore and throat size distributions and connectivity. Special network models for reactive transport are proposed here using the complimentary character of the pore space and the solid phase. This balances the insufficient characterisation information and provides the means for future mechanical-physical-chemical coupling. The anisotropy and heterogeneity of clays is represented using different length parameters and percentage of pores in different directions. Resulting networks are described as mathematical graphs with efficient discrete calculus formulation of transport. Opalinus Clay (OPA) is chosen as an example. Experimental data for the tritiated water (HTO) and U(vi) diffusion through OPA are presented. Calculated diffusion coefficients of HTO and uranium species are within the ranges of the experimentally determined data in different clay directions. This verifies the proposed pore network model and validates that uranium complexes are diffusing as neutral species in OPA. In the case of U(vi) diffusion the method is extended to account for sorption and convection. Rather than changing pore radii by coarse grained mathematical formula, physical sorption is simulated in each pore, which is more accurate and realistic.