Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.

Shear deformation at 29% solid during solidification of magnesium alloy AZ91 and aluminium alloy A356

Partially solid commercial Al-Si and Mg-Al alloys have been deformed in shear during solidification using vane rheometry. The dendritic mush was deformed for a short period at 29% solid and allowed to cool naturally after deformation. Both alloys exhibited yield point behaviour and deformation was highly localised at the surface of maximum shear stress. The short period of deformation was found to have a distinct impact on the as-cast microstructure leading to fragmented dendrites in the deformation region of both alloys. In the case of the Mg-Al alloy, a concentrated region of interdendritic porosity was also observed in the deformation region. Concentrated porosity was not observed in the Al-Si alloy. (c) 2005 Elsevier B.V. All rights reserved.

On the overabundance of light rare earth elements in terrestrial zircons and its implication for Earth's earliest magmatic differentiation

We present whole-rock and zircon rare earth element (REE) data from two early Archaean gneisses (3.81 Ga and 3.64 Ga) from the Itsaq gneiss complex, south-west Greenland. Both gneisses represent extremely rare examples of unaltered, fresh and relatively undeformed igneous rocks of such antiquity. Cathodoluminescence imaging of their zircons indicates a single crystallisation episode with no evidence for either later metamorphic and/or anatectic reworking or inheritance of earlier grains. Uniform, single-population U/Pb age data confirm the structural simplicity of these zircons. One sample, a 3.64 Ga granodioritic gneiss from the Gothabsfjord, yields a chondrite-normalised REE pattern with a positive slope from La to Lu as well as substantial positive Ce and slight negative Eu anomalies, features generally considered to be typical of igneous zircon. In contrast, the second sample, a 3.81 Ga tonalite from south of the Isua Greenstone Belt, has variable but generally much higher light REE abundances, with similar middle to heavy REE. Calculation of zircon/melt distribution coefficients (D-REE(zircon/melt)) from each sample yields markedly different values for the trivalent REE (i.e. Ce and Eu omitted) and simple application of one set of D-REE(zircon/melt) to model the melt composition for the other sample yields concentrations that are in error by up to two orders of magnitude for the light REE (La-Nd). The observed light REE overabundance in the 3.81 Ga tonalite is a commonly observed feature in terrestrial zircons for which a number of explanations ranging from lattice strain to disequilibrium crystallisation have been proposed and are further investigated herein. Regardless of the cause of light REE overabundance, our study shows that simple application of zircon/melt distribution coefficients is not an unambiguous method for ascertaining original melt composition. In this context, recent studies that use REE data to claim that > 4.3 Ga Hadean detrital zircons originally crystallised from an evolved magma, in turn suggesting the operation of geological processes in the early Earth analogous to those of the present day (e.g. subduction and melting of hydrated oceanic crust), must be regarded with caution. Indeed, comparison of terrestrial Hadean and > 3.9 Ga lunar highland zircons shows remarkable similarities in the light REE, even though subduction processes that have been used to explain the terrestrial zircons have never operated on the Moon. (C) 2002 Elsevier Science B.V. All rights reserved.

Dynamic stability of laminated FGM plates based on higher-order shear deformation theory

This paper conducts a dynamic stability analysis of symmetrically laminated FGM rectangular plates with general out-of-plane supporting conditions, subjected to a uniaxial periodic in-plane load and undergoing uniform temperature change. Theoretical formulations are based on Reddy's third-order shear deformation plate theory, and account for the temperature dependence of material properties. A semi-analytical Galerkin-differential quadrature approach is employed to convert the governing equations into a linear system of Mathieu-Hill equations from which the boundary points on the unstable regions are determined by Bolotin's method. Free vibration and bifurcation buckling are also discussed as subset problems. Numerical results are presented in both dimensionless tabular and graphical forms for laminated plates with FGM layers made of silicon nitride and stainless steel. The influences of various parameters such as material composition, layer thickness ratio, temperature change, static load level, boundary constraints on the dynamic stability, buckling and vibration frequencies are examined in detail through parametric studies.

Martensite Crystallography - the apparent controversy between the infinitesimal deformation approach and the phenomenological theory of martensitic transformations

The aim of this article is to demonstrate that the apparent controversy between the infinitesimal deformation (ID) approach and the phenomenological theory of martensitic transformations (PTMTs) in predicting the crystallographic characteristics of a martensitic transformation is entirely based on unjustified approximations associated with the way in which the ID calculations are performed. When applied correctly, the ID approach is shown to be absolutely identical to the PTMT. Nevertheless, there may be some advantages in using the ID approach. In particular, it is somewhat simpler than the PTMT; it is based on a physical concept that is easier to understand and, most important, it may provide a tool for investigating some of the features of martensitic transformations that have eluded explanation via the PTMT.

Segmented polyurethane nanocomposites: Impact of controlled particle size nanofillers on the morphological response to uniaxial deformation

A series of TPU nanocomposites were prepared by incorporating organically modified layered silicates with controlled particle size. To our knowledge, this is the first study into the effects of layered silicate diameter in polymer nanocomposites utilizing the same mineral for each size fraction. The tensile properties of these materials were found to be highly dependent upon the size of the layered silicates. A decrease in disk diameter was associated with a sharp upturn in the stress-strain curve and a pronounced increase in tensile strength. Results from SAXS/SANS experiments showed that the layered silicates did not affect the bulk TPU microphase structure and the morphological response of the host TPU to deformation or promote/hinder strain-induced soft segment crystallization. The improved tensile properties of the nanocomposites containing the smaller nanofillers resulted from the layered silicates aligning in the direction of strain and interacting with the TPU sequences via secondary bonding. This phenomenon contributes predominantly above 400% strain once the microdomain architecture has largely been disassembled. Large tactoids that are unable to align in the strain direction lead to concentrated tensile stresses between the polymer and filler, instead of desirable shear stresses, resulting in void formation and reduced tensile properties. In severe cases, such as that observed for the composite containing the largest silicate, these voids manifest visually as stress whitening.

Bajo de la Alumbrera copper-gold deposit: Stable isotope evidence for a porphyry-related hydrothermal system dominated by magmatic aqueous fluids

Alteration zones at the gold-rich Bajo de la Alumbrera porphyry copper deposit in northwestern Argentina are centered on several porphyritic intrusions. They are zoned from a central copper-iron sulfide and gold-mineralized potassic (biotite-K-feldspar +/- quartz) core outward to propylitic (chlorite-illite-epidote-calcite) assemblages. A mineralized intermediate argillic alteration assemblage (chlorite-illite +/- pyrite) has overprinted the potassic alteration zone across the top and sides of the deposit and is itself zoned outward into phyllic (quartzinuscovite-illite +/- pyrite) alteration. This study contributes new data to previously reported delta(18)O and delta D compositions of fluids responsible for the alteration at Bajo de la Alumbrera, and the data are used to infer likely ore-forming processes. Measured and calculated delta(18)O and delta D values of fluids (+8.3 to +10.2 and -33 to -81 parts per thousand, respectively) confirm a primary magmatic origin for the earliest potassic alteration phase. Lower temperature potassic alteration formed from magmatic fluids with lower delta D values (down to -123 parts per thousand). These depleted compositions are distinct from meteoric water and consistent with degassing and volatile exsolution of magmatic fluids derived from an underlying magma. Variability in the calculated composition of fluid associated with potassic alteration is explained in terms of phase separation (or boiling). if copper-iron sulfide deposition occurred during cooling (as proposed elsewhere), this cooling was largely a result of phase separation. Magmatic water was directly involved in the formation of overprinting intermediate argillic alteration assemblages at Bajo de la Alumbrera. Calculated delta(18)O and delta D values of fluids associated with this alteration range from +4.8 to +8.1 and -31 to -71 per mil, respectively Compositions determined for fluids associated with phyllic alteration (-0.8 to +10.2 and -31 to -119 parts per thousand) overlap with the values determined for the intermediate argillic alteration. We infer that phyllic alteration assemblages developed during two stages; the first was a high-temperature (400 degrees-300 degrees C) stage with D-depleted water (delta D = -66 to -119 parts per thousand). This compositional range may have resulted from magma degassing and/or the injection of new magmatic water into a compositionally evolved hydrothermal system. The isotopic variations also can be explained by increased fluid-rock interaction. The second stage of phyllic alteration occurred at a lower temperature (similar to 200 degrees C), and variations in the modeled isotopic compositions imply mixing of magmatic and meteoric waters. Ore deposition that occurred late in the evolution of the hydrothermal system was probably associated with further cooling of the magmatic fluid, in part caused by fluid-rock interaction and phase separation. Changing pH and/or oxygen fuoracity may have caused additional ore deposition. The ingress of meteoric water appears to postdate the bulk of mineralization and occurred as the system at Bajo de la Alumbrera waned.

Contact angle effects on microdroplet deformation using CFD

In this paper we use computational fluid dynamics (CFD) to study the effect of contact angle on droplet shape as it moves through a contraction. A new non-dimensional number is proposed in order to predict situations where the deformed droplet will form a slug in the contraction and thus have the opportunity to interact with the channel wall. It is proposed that droplet flow into a contraction is a useful method to ensure that a droplet will wet a channel surface without a trapped lubrication film, and thus help ensure that a slug will remain attached to the wall downstream of the contraction. We demonstrate that when a droplet is larger than a contraction, capillary and Reynolds numbers, and fluid properties may not be sufficient to fully describe the droplet dynamics through a contraction. We show that, with everything else constant, droplet shape and breakup can be controlled simply by changing the wetting properties of the channel wall. CFD simulations with contact angles ranging from 30 degrees to 150 degrees show that lower contact angles can induce droplet breakup while higher contact angles can form slugs with contact angle dependent shape. Crown Copyright (c) 2005 Published by Elsevier Inc. All rights reserved.

A parametric study of droplet deformation through a microfluidic contraction: Low viscosity Newtonian droplets

Numerical simulations are conducted to investigate how a droplet of Newtonian liquid. entrained in a higher viscosity Newtonian liquid, behaves when passing through an axisymmetric microfluidic contraction. Simulations are performed using a transient Volume of Fluid finite volume algorithm, and cover ranges of Reynolds and Weber numbers relevant to microfluidic flows. Results are presented for a droplet to surrounding fluid viscosity ratio of 0.001. In contrast to behaviour at higher viscosity ratios obtained previously by the authors, shear and interfacial tension driven instabilities often develop along the droplet Surface. leading to complex shape development, and in some instances, droplet breakup. (c) 2006 Elsevier Ltd. All rights reserved.

Effect of adsorption deformation on thermodynamic characteristics of a fluid in slit pores at sub-critical conditions

Solvation. pressure due to adsorption of fluids in porous materials is the cause of elastic deformation of an adsorbent, which is accessible to direct experimental measurements. Such a deformation contributes to the Helmholtz free energy of the whole adsorbent-adsorbate system due to accumulation of compression or tension energy by the solid. It means that in the general case the solid has to be considered as not solely a source of the external potential field for the fluid confined in the pore volume, but also as thermodynamically nonmert component of the solid-fluid system. We present analysis of nitrogen adsorption isotherms and heat of adsorption in slit graphitic pores accounting for the adsorption deformation by means of nonlocal density functional theory. (c) 2006 Elsevier Ltd. All rights reserved.