Mantle convection modeling with viscoelastic/brittle lithosphere: Numerical methodology and plate tectonic modeling

The earth's tectonic plates are strong, viscoelastic shells which make up the outermost part of a thermally convecting, predominantly viscous layer. Brittle failure of the lithosphere occurs when stresses are high. In order to build a realistic simulation of the planet's evolution, the complete viscoelastic/brittle convection system needs to be considered. A particle-in-cell finite element method is demonstrated which can simulate very large deformation viscoelasticity with a strain-dependent yield stress. This is applied to a plate-deformation problem. Numerical accuracy is demonstrated relative to analytic benchmarks, and the characteristics of the method are discussed.

On the theory of mixing-frequency electron spin-echo envelope modulation spectroscopy

It has recently been stated that the parametrization of the time variables in the one-dimensional (I-D) mixing-frequency electron spin-echo envelope modulation (MIF-ESEEM) experiment is incorrect and hence the wrong frequencies for correlated nuclear transitions are predicted. This paper is a direct response to such a claim, its purpose being to show that the parametrization in land 2-D MIF-ESEEM experiments possesses the same form as that used in other 4-pulse incrementation schemes and predicts the same correlation frequencies. We show that the parametrization represents a shearing transformation of the 2-D time-domain and relate the resulting frequency domain spectrum to the HYSCORE spectrum in terms of a skew-projection. It is emphasized that the parametrization of the time-domain variables may be chosen arbitrarily and affects neither the computation of the correct nuclear frequencies nor the resulting resolution. The usefulness or otherwise of the MIF parameters \gamma\ > 1 is addressed, together with the validity of the original claims of the authors with respect to resolution enhancement in cases of purely homogeneous and inhomogeneous broadening. Numerical simulations are provided to illustrate the main points.

Finite-difference time-domain-based studies of MRI pulsed field gradient-induced eddy currents inside the human body

In modern magnetic resonance imaging (MRI), patients are exposed to strong, rapidly switching magnetic gradient fields that, in extreme cases, may be able to elicit nerve stimulation. This paper presents theoretical investigations into the spatial distribution of induced current inside human tissues caused by pulsed z-gradient fields. A variety of gradient waveforms have been studied. The simulations are based on a new, high-definition, finite-difference time-domain method and a realistic inhomogeneous 10-mm resolution human body model with appropriate tissue parameters. it was found that the eddy current densities are affected not only by the pulse sequences but by many parameters such as the position of the body inside the gradient set, the local biological material properties and the geometry of the body. The discussion contains a comparison of these results with previous results found in the literature. This study and the new methods presented herein will help to further investigate the biological effects caused by the switched gradient fields in a MRI scan. (C) 2002 Wiley Periodicals, Inc.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

Apparent strain localization and shear wave dispersion in elastic fault gouge with microrotations

Shear deformation of fault gouge or other particulate materials often results in observed strain localization, or more precisely, the localization of measured deformation gradients. In conventional elastic materials the strain localization cannot take place therefore this phenomenon is attributed to special types of non-elastic constitutive behaviour. For particulate materials however the Cosserat continuum which takes care of microrotations independent of displacements is a more appropriate model. In elastic Cosserat continuum the localization in displacement gradients is possible under some combinations of the generalized Cosserat elastic moduli. The same combinations of parameters also correspond to a considerable dispersion in shear wave propagation which can be used for independent experimental verification of the proposed mechanism of apparent strain localization in fault gouge.

Surface characterization of 6H-SiC (0001) substrates in indentation and abrasive machining

Surface characterization of 6H-SiC (0001) substrates in indentation and abrasive machining was carried out to investigate microfracture, residual damage, and surface roughness associated with material removal and surface generation. Brittle versus plastic deformation was studied using Vickers indention and nano-indentation. To characterize the abrasive machining response, the 6H-SiC (0001) substrates were ground using diamond wheels with grit sizes of 25, 15 and 7 mum, and then polished with diamond suspensions of 3 and 0.05 mum. It is found that in indentation, there was a scale effect for brittle versus plastic deformation in 6H-SiC substrates. Also, in grinding, the scales of fracture and surface roughness of the substrates decreased with a decrease in diamond grit size. However, in polishing, a reduction in grit size of diamond suspensions gave no significant improvement in surface roughness. Furthermore, the results showed that fracture-free 6H-SiC (0001) surfaces were generated in polishing with the existence of the residual crystal defects, which were associated with the origin of defects in single crystal growth. (C) 2003 Elsevier Ltd. All rights reserved.

A study of high-speed milling characteristics of nitinol

High-speed milling (HSM) has many advantages over conventional machining. Among these advantages, the lower cutting force associated with the machining process is of particular significance for Nitinol alloys because their machined surfaces show less strain hardening. In this article, a systematic study has been carried out to investigate the machining characteristics of a Ni50.6Ti49.4 alloy in HSM. The effects of cutting speed, feed rate, and depth of cut on machined surface characteristics and tool wear are studied. It is found that an increase in cutting speed has resulted in a better surface finish and less work hardening. This is attributed to the reduction of chip cross-sectional area or chip thickness, which thus leads to a lower cutting force or load.

The Western Alps from the Jurassic to Oligocene: spatio-temporal constraints and evolutionary reconstructions

Despite extensive research in the last 150 years, the regional tectonic reconstruction of the Western Alps has remained controversial. The curved orogenic belt consists of several ribbon-like continental terranes (Sesia/Austroalpine, Internal Crystalline Massifs, Brianconnais), which are separated by two or more ophiolitic sutures (Piemonte, Valais, Antrona?, Lanzo/ Canavese?). High-pressure (HP) metamorphism of each terrane occurred during distinct orogenic episodes: at similar to65 Ma in the Sesia/Austroalpine, at similar to45 Ma in the Piemonte zone and at similar to35 Ma in the Internal Crystalline Massifs. It is suggested that these events reflect individual accretionary episodes, which together with kinematic indicators and the speed and direction of plate motions, provide constraints for the discussed reconstruction model. The model involves a prolonged orogenic history that took place during relative convergence of Europe and Adria (here considered as a promontory of the African plate). The first accretionary event involved the Sesia/Austroalpine terrane. Final closure of the Piemonte Ocean occurred during the Eocene (similar to45 Ma) and involved ultra-high-pressure (UHP) metamorphism of the Piemonte oceanic crust. Incorporation of the Brianconnais terrane in the accretionary wedge occurred thereafter, possibly during or after subduction of the Valais Ocean in the late Eocene (45-35 Ma). This subduction was terminated at ca. 35 Ma, when the Internal Crystalline Massifs (i.e. the assumed internal parts of the Brianconnais terrane) were buried into great depths and underwent HP and UHP metamorphism. (C) 2004 Elsevier B.V. All rights reserved.

Continental extension: From core complexes to rigid block faulting

Extension of overthickened continental crust is commonly characterized by an early core complex stage of extension followed by a later stage of crustal-scale rigid block faulting. These two stages are clearly recognized during the extensional destruction of the Alpine orogen in northeast Corsica, where rigid block faulting overprinting core complex formation eventually led to crustal separation and the formation of a new oceanic backarc basin (the Ligurian Sea). Here we investigate the geodynamic evolution of continental extension by using a novel, fully coupled thermomechanical numerical model of the continental crust. We consider that the dynamic evolution is governed by fault weakening, which is generated by the evolution of the natural-state variables (i.e., pressure, deviatoric stress, temperature, and strain rate) and their associated energy fluxes. Our results show the appearance of a detachment layer that controls the initial separation of the brittle crust on characteristic listric faults, and a core complex formation that is exhuming strongly deformed rocks of the detachment zone and relatively undeformed crustal cores. This process is followed by a transitional period, characterized by an apparent tectonic quiescence, in which deformation is not localized and energy stored in the upper crust is transferred downward and causes self-organized mobilization of the lower crust. Eventually, the entire crust ruptures on major crosscutting faults, shifting the tectonic regime from core complex formation to wholesale rigid block faulting.

Effect of sodium dodecylbenzene sulfonate on the motion of three-phase contact lines on the Wilhelmy plate surface

The combined approach of the molecular-kinetic and hydrodynamic theories for description of the motion of three-phase gas-liquid-solid contact lines has been examined using the Wilhelmy plate method. The whole dynamic meniscus has been divided into molecular, hydrodynamic, and static-like regions. The Young-Laplace equation and the molecular-kinetic and hydrodynamic dewetting theories have been applied to describe the meniscus profiles and contact angle. The dissipative forces accompanying the dynamic dewetting have also been investigated. The experiments with a Wilhelmy plate made from an acrylic polymer sheet were carried out using a computerized apparatus for contact angle analysis (OCA 20, DataPhysics, Germany). The extrapolated dynamic contact angle versus velocity of the three-phase contact line for Milli-Q water and 5 x 10(-4) M SDBS solution was experimentally obtained and compared with the combined MHD models with low and moderate Reynolds numbers. The models predict similar results for the extrapolated contact angle. SDBS decreases the equilibrium contact angle and increases the molecular jumping length but does not affect the molecular frequency significantly. The hydrodynamic deformation of the meniscus, viscous dissipation, and friction were also influenced by the SDBS surfactant. (c) 2005 Elsevier Inc. All rights reserved.

Effects of Si content on defect band formation in hypoeutectic Al-Si die castings

Al-3-11% Si alloys have been high-pressure die-cast and characterized microstructurally. Alstruc was used to calculate the solidification characteristics and fraction of eutectic. Defect bands were observed at all Si contents, although their constitution, position and distinctiveness were a function of Si content. The defect bands contain a higher fraction Al-Si eutectic than the surroundings in all alloys, and porosity was additionally found in the band in AlSi3. With decreasing Si content, the defect bands formed closer to the casting surface, became more prevalent and also the width of the bands decreased. These differences are discussed by considering the effect of Si content on the distribution of solid in the mushy wall layers and on the feeding potentials of the alloys. The observations are consistent with the mechanism proposed by Gourlay et al. in which bands form due to deformation within the solidifying mushy wall layers. (c) 2005 Elsevier B.V. All rights reserved.

Quantification of orientation of pore patterns in X-ray images of deformed clay

The use of X-ray imaging to assess variability in ceramic fabrication is common in archaeological studies aimed at examining ancient pottery technologies. In this paper, a method based on the measurement of individual pores orientation is presented. This method is successfully applied to ceramic specimens of known origin whose structure signified different deformation histories.

A numerical study of effect of grain boundaries on elastic and plastic properties in nanocomposite materials

Nanocomposite materials have received considerable attention in recent years due to their novel properties. Grain boundaries are considered to play an important role in nanostructured materials. This work focuses on the finite element analysis of the effect of grain boundaries on the overall mechanical properties of aluminium/alumina composites. A grain boundary is incorporated into the commonly used unit cell model to investigate its effect on material properties. By combining the unit cell model with an indentation model, coupled with experimental indentation measurements, the ''effective'' plastic property of the grain boundary is estimated. In addition, the strengthening mechanism is also discussed based on the Estrin-Mecking model.

Modelling of cold-formed purlin-sheeting systems .2. Simplified model

A number of theoretical and experimental investigations have been made into the nature of purlin-sheeting systems over the past 30 years. These systems commonly consist of cold-formed zed or channel section purlins, connected to corrugated sheeting. They have proven difficult to model due to the complexity of both the purlin deformation and the restraint provided to the purlin by the sheeting. Part 1 of this paper presented a non-linear elasto plastic finite element model which, by incorporating both the purlin and the sheeting in the analysis, allowed the interaction between the two components of the system to be modelled. This paper presents a simplified version of the first model which has considerably decreased requirements in terms of computer memory, running time and data preparation. The Simplified Model includes only the purlin but allows for the sheeting's shear and rotational restraints by modelling these effects as springs located at the purlin-sheeting connections. Two accompanying programs determine the stiffness of these springs numerically. As in the Full Model, the Simplified Model is able to account for the cross-sectional distortion of the purlin, the shear and rotational restraining effects of the sheeting, and failure of the purlin by local buckling or yielding. The model requires no experimental or empirical input and its validity is shown by its goon con elation with experimental results. (C) 1997 Elsevier Science Ltd.

A continuum model for periodic two-dimensional block structures

A continuum model for regular block structures is derived by replacing the difference quotients of the discrete equations by corresponding differential quotients. The homogenization procedure leads to an anisotropic Cosserat Continuum. For elastic block interactions the dispersion relations of the discrete and the continuous models are derived and compared. Yield criteria for block tilting and sliding are formulated. An extension of the theory for large deformation is proposed. (C) 1997 by John Wiley & Sons, Ltd.