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.

Optical coupling of Au nanoparticles on vertical graphenes to maximize SERS response

Gold particle interaction with few-layer graphenes is of interest for the development of numerous optical nanodevices. The results of numerical studies of the coupling of gold nanoparticles with few-layer vertical graphene sheets are presented. The field strengths are computed and the optimum nanoparticle configurations for the formation of SERS hotpots are obtained. The nanoparticles are modeled as 8 nm diameter spheres atop 1.5 nm (5 layers) graphene sheet. The vertical orientation is of particular interest as it is possible to use both sides of the graphene structure and potentially double the number of particles in the system. Our results show that with the addition of an opposing particle a much stronger signal can be obtained as well as the particle separation can be controlled by the number of atomic carbon layers. These results provide further insights and contribute to the development of next-generation plasmonic devices based on nanostructures with hybrid dimensionality.

Biaxial fiber Bragg grating accelerometer using axial and transverse forces

We demonstrate the first biaxial fiber Bragg grating (FBG) accelerometer using axial and transverse forces. An inertial object is fixed at the middle of two FBGs inscribed in one fiber. The difference between the resonant wavelengths of the two FBGs can distinguish the acceleration in the axial direction, while being insensitive in the transverse direction. The average of the resonant wavelengths of the two FBGs can distinguish the acceleration in the transverse direction, while being insensitive in the axial direction. In the experiments, when the transverse direction was vertical, the crest-to-trough sensitivity at 5 Hz and resonant frequency of the average were 0.545 nm/g and 34.42 Hz, respectively. When the axial direction was vertical, those of the difference were 0.0454 nm/g and 900 Hz, respectively. For each FBG, the crest-to-trough sensitivity at 5 Hz and resonant frequency in the transverse/vertical direction were 24 and 1/26 times those in the axial/vertical direction, respectively.

Improvement of the mode II interface fracture toughness of glass fiber reinforced plastics/aluminum laminates through vapor grown carbon fiber interleaves

The effects of acid treatment, vapor grown carbon fiber (VGCF) interlayer and the angle, i.e., 0° and 90°, between the rolling stripes of an aluminum (Al) plate and the fiber direction of glass fiber reinforced plastics (GFRP) on the mode II interlaminar mechanical properties of GFRP/Al laminates were investigated. The experimental results of an end notched flexure test demonstrate that the acid treatment and the proper addition of VGCF can effectively improve the critical load and mode II fracture toughness of GFRP/Al laminates. The specimens with acid treatment and 10 g m−2 VGCF addition possess the highest mode II fracture toughness, i.e., 269% and 385% increases in the 0° and 90° specimens, respectively compared to those corresponding pristine ones. Due to the induced anisotropy by the rolling stripes on the aluminum plate, the 90° specimens possess 15.3%–73.6% higher mode II fracture toughness compared to the 0° specimens. The improvement mechanisms were explored by the observation of crack propagation path and fracture surface with optical, laser scanning and scanning electron microscopies. Moreover, finite element analyses were carried out based on the cohesive zone model to verify the experimental fracture toughness and to predict the interface shear strength between the aluminum plates and GFRP laminates.

A novel method of Fiber Bragg Grating (FBG) strain modulation and its application in FBG accelerometers

This thesis has systemically investigated the possibility of improving one type of optical fiber sensors by using a novel mechanism. Many parameters of the sensor have been improved, and one outcome of this innovation is that civil structures, such as bridges and high-rise buildings, may be operated more safely and used longer.

The pattern of coupling dynamics between postural motion, isotonic hand movements and physiological tremor

This study was designed to examine differences in the coupling dynamics between upper limb motion, physiological tremor and whole body postural sway in young healthy adults. Acceleration of the hand and fingers, forearm EMG activity and postural sway data were recorded. Estimation of the degree of bilateral and limb motion-postural sway coupling was determined by cross correlation, coherence and Cross-ApEn analyses. The results of the analysis revealed that, under postural tremor conditions, there was no significant coupling between limbs, muscles or sway across all metrics of coupling. In contrast, performing a rapid alternating flexion/extension movement about the wrist joint (with one or both limbs) resulted in stronger coupling between limb motion and postural sway. These results support the view that, for physiological tremor responses, the control of postural sway is maintained independent to tremor in the upper limb. However, increasing the level of movement about a distal segment of one arm (or both) leads to increased coupling throughout the body. The basis for this increased coupling would appear to be related to the enhanced neural drive to task-specific muscles within the upper limb.

Compensation for higher order mode coupling between inclined coupling slots and radiating slots in planar slotted waveguide arrays

Summary form only given. Geometric simplicity, efficiency and polarization purity make slot antenna arrays ideal solutions for many radar, communications and navigation applications, especially when high power, light weight and limited scan volume are priorities. Resonant arrays of longitudinal slots have a slot spacing of one-half guide wavelength at the design frequency, so that the slots are located at the standing wave peaks. Planar arrays are implemented using a number of rectangular waveguides (branch line guides), arranged side-by-side, while waveguides main lines located behind and at right angles to the branch lines excite the radiating waveguides via centered-inclined coupling slots. Planar slotted waveguide arrays radiate broadside beams and all radiators are designed to be in phase.

Excitation of bound plasmons along nanoscale stripe waveguides: A comparison of end and grating coupling techniques

In this paper we excite bound long range stripe plasmon modes with a highly focused laser beam. We demonstrate highly confined plasmons propagating along a 50 μm long silver stripe 750 nm wide and 30 nm thick. Two excitation techniques were studied: focusing the laser spot onto the waveguide end and focusing the laser spot onto a silver grating. By comparing the intensity of the out-coupling photons at the end of the stripe for both grating and end excitation we are able to show that gratings provide an increase of a factor of two in the output intensity and thus out-coupling of plasmons excited by this technique are easier to detect. Authors expect that the outcome of this paper will prove beneficial for the development of passive nano-optical devices based on stripe waveguides, by providing insight into the different excitation techniques available and the advantages of each technique.

Interlaminar mechanical properties of carbon fiber reinforced plastic laminates modified with graphene oxide interleaf

By taking the advantage of the excellent mechanical properties and high specific surface area of graphene oxide (GO) sheets, we develop a simple and effective strategy to improve the interlaminar mechanical properties of carbon fiber reinforced plastic (CFRP) laminates. With the incorporation of graphene oxide reinforced epoxy interleaf into the interface of CFRP laminates, the Mode-I fracture toughness and resistance were greatly increased. The experimental results of double cantilever beam (DCB) tests demonstrated that, with 2 g/m2 addition of GO, the Mode-I fracture toughness and resistance of the specimen increase by 170.8% and 108.0%, respectively, compared to those of the plain specimen. The improvement mechanisms were investigated by the observation of fracture surface with scanning electron microscopies. Moreover, finite element analyses were performed based on the cohesive zone model to verify the experimental fracture toughness and to predict the interfacial tensile strength of CFRP laminates.

Extending genetic linkage analysis to diffusion tensor images to map single gene effects on brain fiber architecture

We extended genetic linkage analysis - an analysis widely used in quantitative genetics - to 3D images to analyze single gene effects on brain fiber architecture. We collected 4 Tesla diffusion tensor images (DTI) and genotype data from 258 healthy adult twins and their non-twin siblings. After high-dimensional fluid registration, at each voxel we estimated the genetic linkage between the single nucleotide polymorphism (SNP), Val66Met (dbSNP number rs6265), of the BDNF gene (brain-derived neurotrophic factor) with fractional anisotropy (FA) derived from each subject's DTI scan, by fitting structural equation models (SEM) from quantitative genetics. We also examined how image filtering affects the effect sizes for genetic linkage by examining how the overall significance of voxelwise effects varied with respect to full width at half maximum (FWHM) of the Gaussian smoothing applied to the FA images. Raw FA maps with no smoothing yielded the greatest sensitivity to detect gene effects, when corrected for multiple comparisons using the false discovery rate (FDR) procedure. The BDNF polymorphism significantly contributed to the variation in FA in the posterior cingulate gyrus, where it accounted for around 90-95% of the total variance in FA. Our study generated the first maps to visualize the effect of the BDNF gene on brain fiber integrity, suggesting that common genetic variants may strongly determine white matter integrity.

Brain fiber architecture, genetics, and intelligence: a high angular resolution diffusion imaging (HARDI) study

We developed an analysis pipeline enabling population studies of HARDI data, and applied it to map genetic influences on fiber architecture in 90 twin subjects. We applied tensor-driven 3D fluid registration to HARDI, resampling the spherical fiber orientation distribution functions (ODFs) in appropriate Riemannian manifolds, after ODF regularization and sharpening. Fitting structural equation models (SEM) from quantitative genetics, we evaluated genetic influences on the Jensen-Shannon divergence (JSD), a novel measure of fiber spatial coherence, and on the generalized fiber anisotropy (GFA) a measure of fiber integrity. With random-effects regression, we mapped regions where diffusion profiles were highly correlated with subjects' intelligence quotient (IQ). Fiber complexity was predominantly under genetic control, and higher in more highly anisotropic regions; the proportion of genetic versus environmental control varied spatially. Our methods show promise for discovering genes affecting fiber connectivity in the brain.

Mapping genetic influences on brain fiber architecture with high angular resolution diffusion imaging (HARDI)

We report the first 3D maps of genetic effects on brain fiber complexity. We analyzed HARDI brain imaging data from 90 young adult twins using an information-theoretic measure, the Jensen-Shannon divergence (JSD), to gauge the regional complexity of the white matter fiber orientation distribution functions (ODF). HARDI data were fluidly registered using Karcher means and ODF square-roots for interpol ation; each subject's JSD map was computed from the spatial coherence of the ODFs in each voxel's neighborhood. We evaluated the genetic influences on generalized fiber anisotropy (GFA) and complexity (JSD) using structural equation models (SEM). At each voxel, genetic and environmental components of data variation were estimated, and their goodness of fit tested by permutation. Color-coded maps revealed that the optimal models varied for different brain regions. Fiber complexity was predominantly under genetic control, and was higher in more highly anisotropic regions. These methods show promise for discovering factors affecting fiber connectivity in the brain.