Multiscale coupling and multiphysics approaches in earth sciences : applications

Geoscientists are confronted with the challenge of assessing nonlinear phenomena that result from multiphysics coupling across multiple scales from the quantum level to the scale of the earth and from femtoseconds to the 4.5 Ga of history of our planet. We neglect in this review electromagnetic modelling of the processes in the Earth’s core, and focus on four types of couplings that underpin fundamental instabilities in the Earth. These are thermal (T), hydraulic (H), mechanical (M) and chemical (C) processes which are driven and controlled by the transfer of heat to the Earth’s surface. Instabilities appear as faults, folds, compaction bands, shear/fault zones, plate boundaries and convective patterns. Convective patterns emerge from buoyancy overcoming viscous drag at a critical Rayleigh number. All other processes emerge from non-conservative thermodynamic forces with a critical critical dissipative source term, which can be characterised by the modified Gruntfest number Gr. These dissipative processes reach a quasi-steady state when, at maximum dissipation, THMC diffusion (Fourier, Darcy, Biot, Fick) balance the source term. The emerging steady state dissipative patterns are defined by the respective diffusion length scales. These length scales provide a fundamental thermodynamic yardstick for measuring instabilities in the Earth. The implementation of a fully coupled THMC multiscale theoretical framework into an applied workflow is still in its early stages. This is largely owing to the four fundamentally different lengths of the THMC diffusion yardsticks spanning micro-metre to tens of kilometres compounded by the additional necessity to consider microstructure information in the formulation of enriched continua for THMC feedback simulations (i.e., micro-structure enriched continuum formulation). Another challenge is to consider the important factor time which implies that the geomaterial often is very far away from initial yield and flowing on a time scale that cannot be accessed in the laboratory. This leads to the requirement of adopting a thermodynamic framework in conjunction with flow theories of plasticity. This framework allows, unlike consistency plasticity, the description of both solid mechanical and fluid dynamic instabilities. In the applications we show the similarity of THMC feedback patterns across scales such as brittle and ductile folds and faults. A particular interesting case is discussed in detail, where out of the fluid dynamic solution, ductile compaction bands appear which are akin and can be confused with their brittle siblings. The main difference is that they require the factor time and also a much lower driving forces to emerge. These low stress solutions cannot be obtained on short laboratory time scales and they are therefore much more likely to appear in nature than in the laboratory. We finish with a multiscale description of a seminal structure in the Swiss Alps, the Glarus thrust, which puzzled geologists for more than 100 years. Along the Glarus thrust, a km-scale package of rocks (nappe) has been pushed 40 km over its footwall as a solid rock body. The thrust itself is a m-wide ductile shear zone, while in turn the centre of the thrust shows a mm-cm wide central slip zone experiencing periodic extreme deformation akin to a stick-slip event. The m-wide creeping zone is consistent with the THM feedback length scale of solid mechanics, while the ultralocalised central slip zones is most likely a fluid dynamic instability.

Magnetic convection heat transfer in an open ended enclosure filled with paramagnetic fluids

Numerical simulations of thermomagnetic convection of paramagnetic fluids placed in a micro-gravity condition (g nearly 0) and under a uniform vertical gradient magnetic field in an open ended square enclosure with ramp heating temperature condition applied on a vertical wall is investigated in this study. In presence of the strong magnetic gradient field thermal convection of the paramagnetic fluid might take place even in a zero-gravity environment as a direct consequence of temperature differences occurring within the fluid. The thermal boundary layer develops adjacent to the hot wall as soon as the ramp temperature condition is applied on it. There are two scenario that can be observed based on the ramp heating time. The steady state of the thermal boundary layer can be reached before the ramp time is finished or vice versa. If the ramp time is larger than the quasi-steady time then the thermal boundary layer is in a quasi-steady mode with convection balancing conduction after the quasi-steady time. Further increase of the heat input simply accelerates the flow to maintain the proper thermal balance. Finally, the boundary layer becomes completely steady state when the ramp time is finished. Effects of magnetic Rayleigh number, Prandtl number and paramagnetic fluid parameter on the flow pattern and heat transfer are presented.

Crystalline Si nanoparticles below crystallization threshold : effects of collisional heating in non-thermal atmospheric-pressure microplasmas

Nucleation and growth of highly crystalline silicon nanoparticles in atmospheric-pressure low-temperature microplasmas at gas temperatures well below the Si crystallization threshold and within a short (100 μs) period of time are demonstrated and explained. The modeling reveals that collision-enhanced ion fluxes can effectively increase the heat flux on the nanoparticle surface and this heating is controlled by the ion density. It is shown that nanoparticles can be heated to temperatures above the crystallization threshold. These combined experimental and theoretical results confirm the effective heating and structure control of Si nanoparticles at atmospheric pressure and low gas temperatures.

Lanthanum oxide nanostructured films synthesized using hot dense and extremely non-equilibrium plasma for nanoelectronic device applications

Lanthanum oxide (La2O3) nanostructured films are synthesized on a p-type silicon wafer by ablation of La2O3 pellet due to interaction with hot dense argon plasmas in a modified dense plasma focus (DPF) device. The nanostructured films are investigated using scanning electron microscopy (SEM), atomic force microscopy (AFM), and X-ray diffraction (XRD) spectra. SEM study shows the formation of nano-films having nano-size structures with the average nanostructures size ~25, ~53, and ~45 nm for one, two, and three DPF shots, respectively. The nanostructures sizes and morphology of nano-films are consistent between the AFM and SEM analyses. XRD spectra confirms nano-sized La2O3 with an average grain size ~34, ~51, and ~42 nm for one, two, and three DPF shots, respectively. The electrical properties such as current-voltage and capacitance-voltage (C-V) characteristics of the Al-La2O3-Si metal-oxide- semiconductor (MOS) capacitor structure are measured. The current conduction mechanism of the MOS capacitors is also demonstrated. The C-V characteristics are further used to obtain the electrical parameters such as the dielectric constant, oxide thickness, flat-band capacitance, and flat-band voltage of the MOS capacitors. These measurements demonstrate significantly lower leakage currents without any commonly used annealing or doping, thereby revealing a significant improvement of the MOS nanoelectronic device performance due to the incorporation of the DPF-produced La2O3 nano-films.

Si quantum dots embedded in an amorphous SiC matrix : nanophase control by non-equilibrium plasma hydrogenation

Nanophase nc-Si/a-SiC films that contain Si quantum dots (QDs) embedded in an amorphous SiC matrix were deposited on single-crystal silicon substrates using inductively coupled plasma-assisted chemical vapor deposition from the reactive silane and methane precursor gases diluted with hydrogen at a substrate temperature of 200 °C. The effect of the hydrogen dilution ratio X (X is defined as the flow rate ratio of hydrogen-to-silane plus methane gases), ranging from 0 to 10.0, on the morphological, structural, and compositional properties of the deposited films, is extensively and systematically studied by scanning electron microscopy, high-resolution transmission electron microscopy, X-ray diffraction, Raman spectroscopy, Fourier-transform infrared absorption spectroscopy, and X-ray photoelectron spectroscopy. Effective nanophase segregation at a low hydrogen dilution ratio of 4.0 leads to the formation of highly uniform Si QDs embedded in the amorphous SiC matrix. It is also shown that with the increase of X, the crystallinity degree and the crystallite size increase while the carbon content and the growth rate decrease. The obtained experimental results are explained in terms of the effect of hydrogen dilution on the nucleation and growth processes of the Si QDs in the high-density plasmas. These results are highly relevant to the development of next-generation photovoltaic solar cells, light-emitting diodes, thin-film transistors, and other applications.

Rapid, low-temperature synthesis of nc-Si in high-density, non-equilibrium plasmas : enabling nanocrystallinity at very low hydrogen dilution

Nanocrystalline silicon thin films were deposited on single-crystal silicon and glass substrates simultaneously by inductively coupled plasma-assisted chemical vapor deposition from the reactive silane reactant gas diluted with hydrogen at a substrate temperature of 200 °C. The effect of hydrogen dilution ratio X (X is defined as the flow rate ratio of hydrogen to silane gas), ranging from 1 to 20, on the structural and optical properties of the deposited films, is extensively investigated by Raman spectroscopy, X-ray diffraction, Fourier transform infrared absorption spectroscopy, UV/VIS spectroscopy, and scanning electron microscopy. Our experimental results reveal that, with the increase of the hydrogen dilution ratio X, the deposition rate Rd and hydrogen content CH are reduced while the crystalline fraction Fc, mean grain size δ and optical bandgap ETauc are increased. In comparison with other plasma enhanced chemical vapor deposition methods of nanocrystalline silicon films where a very high hydrogen dilution ratio X is routinely required (e.g. X > 16), we have achieved nanocrystalline silicon films at a very low hydrogen dilution ratio of 1, featuring a high deposition rate of 1.57 nm/s, a high crystalline fraction of 67.1%, a very low hydrogen content of 4.4 at.%, an optical bandgap of 1.89 eV, and an almost vertically aligned columnar structure with a mean grain size of approximately 19 nm. We have also shown that a sufficient amount of atomic hydrogen on the growth surface essential for the formation of nanocrystalline silicon is obtained through highly-effective dissociation of silane and hydrogen molecules in the high-density inductively coupled plasmas. © 2009 The Royal Society of Chemistry.

Assembly and annotation of a non-model gastropod (Nerita melanotragus) transcriptome: a comparison of De novo assemblers

Background The sequencing, de novo assembly and annotation of transcriptome datasets generated with next generation sequencing (NGS) has enabled biologists to answer genomic questions in non-model species with unprecedented ease. Reliable and accurate de novo assembly and annotation of transcriptomes, however, is a critically important step for transcriptome assemblies generated from short read sequences. Typical benchmarks for assembly and annotation reliability have been performed with model species. To address the reliability and accuracy of de novo transcriptome assembly in non-model species, we generated an RNAseq dataset for an intertidal gastropod mollusc species, Nerita melanotragus, and compared the assembly produced by four different de novo transcriptome assemblers; Velvet, Oases, Geneious and Trinity, for a number of quality metrics and redundancy. Results Transcriptome sequencing on the Ion Torrent PGM™ produced 1,883,624 raw reads with a mean length of 133 base pairs (bp). Both the Trinity and Oases de novo assemblers produced the best assemblies based on all quality metrics including fewer contigs, increased N50 and average contig length and contigs of greater length. Overall the BLAST and annotation success of our assemblies was not high with only 15-19% of contigs assigned a putative function. Conclusions We believe that any improvement in annotation success of gastropod species will require more gastropod genome sequences, but in particular an increase in mollusc protein sequences in public databases. Overall, this paper demonstrates that reliable and accurate de novo transcriptome assemblies can be generated from short read sequencers with the right assembly algorithms. Keywords: Nerita melanotragus; De novo assembly; Transcriptome; Heat shock protein; Ion torrent

Non-invasive techniques for imaging in-vivo human corneal sub basal nerve mapping and migration rate in diabetic neuropathy

Purpose Corneal confocal microscopy (CCM) is a rapid non-invasive ophthalmic technique, which has been shown to diagnose and stratify the severity of diabetic neuropathy. Current morphometric techniques assess individual static images of the subbasal nerve plexus; this work explores the potential for non-invasive assessment of the wide-field morphology and dynamic changes of this plexus in vivo. Methods In this pilot study, laser scanning CCM was used to acquire maps (using a dynamic fixation target and semi-automated tiling software) of the central corneal sub-basal nerve plexus in 4 diabetic patients with and 6 without neuropathy and in 2 control subjects. Nerve migration was measured in an additional 7 diabetic patients with neuropathy, 4 without neuropathy and in 2 control subjects by repeating a modified version of the mapping procedure within 2-8 weeks, thus facilitating re-identification of distinctive nerve landmarks in the 2 montages. The rate of nerve movement was determined from these data and normalised to a weekly rate (µm/week), using customised software. Results Wide-field corneal nerve fibre length correlated significantly with the Neuropathy Disability Score (r = -0.58, p < 0.05), vibration perception (r = -0.66, p < 0.05) and peroneal conduction velocity (r = 0.67, p < 0.05). Central corneal nerve fibre length did not correlate with any of these measures of neuropathy (p > 0.05 for all). The rate of corneal nerve migration was 14.3 ± 1.1 µm/week in diabetic patients with neuropathy, 19.7 ± 13.3µm/week in diabetic patients without neuropathy, and 24.4 ± 9.8µm/week in control subjects; however, these differences were not significantly different (p = 0.543). Conclusions Our data demonstrate that it is possible to capture wide-field images of the corneal nerve plexus, and to quantify the rate of corneal nerve migration by repeating this procedure over a number of weeks. Further studies on larger sample sizes are required to determine the utility of this approach for the diagnosis and monitoring of diabetic neuropathy.

Heat strain evaluation of overt and covert body armour in a hot and humid environment

The aim of this study was to elucidate the thermophysiological effects of wearing lightweight non-military overt and covert personal body armour (PBA) in a hot and humid environment. Eight healthy males walked on a treadmill for 120 min at 22% of their heart rate reserve in a climate chamber simulating 31 °C (60%RH) wearing either no armour (control), overt or covert PBA in addition to a security guard uniform, in a randomised controlled crossover design. No significant difference between conditions at the end of each trial was observed in core temperature, heart rate or skin temperature (P > 0.05). Covert PBA produced a significantly greater amount of body mass change (−1.81 ± 0.44%) compared to control (−1.07 ± 0.38%, P = 0.009) and overt conditions (−1.27 ± 0.44%, P = 0.025). Although a greater change in body mass was observed after the covert PBA trial; based on the physiological outcome measures recorded, the heat strain encountered while wearing lightweight, non-military overt or covert PBA was negligible compared to no PBA. Practitioner summary The wearing of bullet proof vests or body armour is a requirement of personnel engaged in a wide range of occupations including police, security, customs and even journalists in theatres of war. This randomised controlled crossover study is the first to examine the thermophysiological effects of wearing lightweight non-military overt and covert personal body armour (PBA) in a hot and humid environment. We conclude that the heat strain encountered while wearing both overt and covert lightweight, non-military PBA was negligible compared to no PBA.

Fouling of waste heat recovery : numerical and experimental results

Numerical results are presented to investigate the performance of a partly-filled porous heat exchanger for waste heat recovery units. A parametric study was conducted to investigate the effects of inlet velocity and porous block height on the pressure drop of the heat exchanger. The focus of this work is on modelling the interface of a porous and non-porous region. As such, numerical simulation of the problem is conducted along with hot-wire measurements to better understand the physics of the problem. Results from the two sources are then compared to existing theoretical predictions available in the literature which are unable to predict the existence of two separation regions before and after the porous block. More interestingly, a non-uniform interface velocity was observed along the streamwise direction based on both numerical and experimental data.

Natural convection in a triangular enclosure heated from below and non-uniformly cooled from top

This study is concerned with transient natural convection in an isosceles triangular enclosure subject to non-uniformly cooling at the inclined surfaces and uniformly heating at the base. The numerical simulations of the unsteady flows over a range of Rayleigh numbers and aspect ratios are carried out using Finite Volume Method. Since the upper inclined surfaces are linearly cooled and the bottom surface is heated, the flow is potentially unstable. It is revealed from the numerical simulations that the transient flow development in the enclosure can be classified into three distinct stages; an early stage, a transitional stage, and a steady stage. The flow inside the enclosure depends significantly on the governing parameters, Rayleigh number and aspect ratio. The effect of Rayleigh number and aspect ratio on the flow development and heat transfer rate are discussed. The key finding for this study is to analyze the pitchfork bifurcation of the flow about the geometric center line. The heat transfer through the roof and the ceiling as a form of Nusselt number is reported in this study.

Multiphase models for simulating hot and slightly miscible DNAPL (dense non-aqueous phase fluids) in a saturated rock fracture under deformation

In the present study a two dimensional model is first developed to show the behaviour of dense non-aqueous phase liquids (DNAPL) within a rough fracture. To consider the rough fracture, the fracture is imposed with variable apertures along its plane. It is found that DNAPL follows preferential pathways. In next part of the study the above model is further extended for non-isothermal DNAPL flow and DNAPL-water interphase mass transfer phenomenon. These two models are then coupled with joint deformation due to normal stresses. The primary focus of these models is specifically to elucidate the influence of joint alteration due to external stress and fluid pressures on flow driven energy transport and interphase mass transfer. For this, it is assumed that the critical value for joint alteration is associated with external stress and average of water and DNAPL pressures in multiphase system and the temporal and spatial evolution of joint alteration are determined for its further influence on energy transport and miscible phase transfer. The developed model has been studied to show the influence of deformation on DNAPL flow. Further this preliminary study demonstrates the influence of joint deformation on heat transport and phase miscibility via multiphase flow velocities. It is seen that the temperature profile changes and shows higher diffusivity due to deformation and although the interphase miscibility value decreases but the lateral dispersion increases to a considerably higher extent.

Role of Melt Convection on the Thermal Performance of Heat Sinks with Phase Change Material

The role of melt convection oil the performance of beat sinks with Phase Change Material (PCM) is presented in this paper. The beat sink consists of aluminum plate fins embedded in PCM and heat flux is supplied from the bottom. The design of such a heat sink requires optimization with respect to its geometrical parameters. The objective of the optimization is to maximize the heat sink operation time for the prescribed heat flux and the critical chip temperature. The parameters considered for optimization are fin number and fill thickness. The height and base plate thickness of heat sink are kept constant in the present analysis. An enthalpy based CFD model is developed, which is capable Of Simulating phase change and associated melt convection. The CFD model is Coupled with Genetic Algorithm (GA) for carrying out the optimization. Two cases are considered, one without melt convection (conduction regime) and the other with convection. It is found that the geometrical optimizations of heat sinks are different for the two cases, indicating the importance of inch convection in the design of heat sinks with PCMs.

Boundary-layer convective heat transfer in a convergent channel

In this paper the numerical solution of the heat transfer problem in a convergent channel with uniform and non-uniform wall temperatures under boundary-layer approximations has been presented. Also, a semi-analytical solution for uniform wall temperature has been obtained.

Constitutive modeling of shape memory alloy wire with non-local rate kinetics

A constitutive modeling approach for shape memory alloy (SMA) wire by taking into account the microstructural phase inhomogeneity and the associated solid-solid phase transformation kinetics is reported in this paper. The approach is applicable to general thermomechanical loading. Characterization of various scales in the non-local rate sensitive kinetics is the main focus of this paper. Design of SMA materials and actuators not only involve an optimal exploitation of the hysteresis loops during loading-unloading, but also accounts for fatigue and training cycle identifications. For a successful design of SMA integrated actuator systems, it is essential to include the microstructural inhomogeneity effects and the loading rate dependence of the martensitic evolution, since these factors play predominant role in fatigue. In the proposed formulation, the evolution of new phase is assumed according to Weibull distribution. Fourier transformation and finite difference methods are applied to arrive at the analytical form of two important scaling parameters. The ratio of these scaling parameters is of the order of 10(6) for stress-free temperature-induced transformation and 10(4) for stress-induced transformation. These scaling parameters are used in order to study the effect of microstructural variation on the thermo-mechanical force and interface driving force. It is observed that the interface driving force is significant during the evolution. Increase in the slopes of the transformation start and end regions in the stress-strain hysteresis loop is observed for mechanical loading with higher rates.