Comparison of surfactant lipids between pleural and pulmonary lining fluids

Saturated phospholipids (PCs), particularly dipalmitoylphosphatidylcholine (DPPC), predominate in surfactant lining the alveoli, although little is known about the relationship between saturated and unsaturated PCs on the outer surface of the lung, the pleura. Seven healthy cats were anesthetized and a bronchoalveolar lavage (BAL) was performed, immediately followed by a pleural lavage (PL). Lipid was extracted from lavage fluid and then analyzed for saturated, primarily dipalmitoylphosphatidylcholine (DPPC), and unsaturated PC species using high-performance liquid chromatography (HPLC) with combined fluorescence and ultraviolet detection. Dilution of epithelial lining fluid (ELF) in lavage fluids was corrected for using the urea method. The concentration of DPPC in BAL fluid (85.3 +/- 15.7 mu g/mL) was significantly higher (P=0.021) than unsaturated PCs (similar to 40 mu g/mL). However, unsaturated PCs (similar to 34 mu g/mL), particularly stearoyl-linoleoyl-phosphatidylcholine (SLPC; 17.4 +/- 6.8), were significantly higher (P = 0.021) than DPPC (4.3 +/- 1.8 mu g/mL) in PL fluid. These results show that unsaturated PCs appear functionally more important in the pleural cavity, which may have implications for surfactant replenishment following pleural disease or thoracic surgery. (c) 2005 Published by Elsevier Ltd.

Drop formation and breakup of low viscosity elastic fluids: Effects of molecular weight and concentration

The dynamics of drop formation and pinch-off have been investigated for a series of low viscosity elastic fluids possessing similar shear viscosities, but differing substantially in elastic properties. On initial approach to the pinch region, the viscoelastic fluids all exhibit the same global necking behavior that is observed for a Newtonian fluid of equivalent shear viscosity. For these low viscosity dilute polymer solutions, inertial and capillary forces form the dominant balance in this potential flow regime, with the viscous force being negligible. The approach to the pinch point, which corresponds to the point of rupture for a Newtonian fluid, is extremely rapid in such solutions, with the sudden increase in curvature producing very large extension rates at this location. In this region the polymer molecules are significantly extended, causing a localized increase in the elastic stresses, which grow to balance the capillary pressure. This prevents the necked fluid from breaking off, as would occur in the equivalent Newtonian fluid. Alternatively, a cylindrical filament forms in which elastic stresses and capillary pressure balance, and the radius decreases exponentially with time. A (0+1)-dimensional finitely extensible nonlinear elastic dumbbell theory incorporating inertial, capillary, and elastic stresses is able to capture the basic features of the experimental observations. Before the critical "pinch time" t(p), an inertial-capillary balance leads to the expected 2/3-power scaling of the minimum radius with time: R-min similar to(t(p)-t)(2/3). However, the diverging deformation rate results in large molecular deformations and rapid crossover to an elastocapillary balance for times t>t(p). In this region, the filament radius decreases exponentially with time R-min similar to exp[(t(p)-t)/lambda(1)], where lambda(1) is the characteristic time constant of the polymer molecules. Measurements of the relaxation times of polyethylene oxide solutions of varying concentrations and molecular weights obtained from high speed imaging of the rate of change of filament radius are significantly higher than the relaxation times estimated from Rouse-Zimm theory, even though the solutions are within the dilute concentration region as determined using intrinsic viscosity measurements. The effective relaxation times exhibit the expected scaling with molecular weight but with an additional dependence on the concentration of the polymer in solution. This is consistent with the expectation that the polymer molecules are in fact highly extended during the approach to the pinch region (i.e., prior to the elastocapillary filament thinning regime) and subsequently as the filament is formed they are further extended by filament stretching at a constant rate until full extension of the polymer coil is achieved. In this highly extended state, intermolecular interactions become significant, producing relaxation times far above theoretical predictions for dilute polymer solutions under equilibrium conditions. (C) 2006 American Institute of Physics

Elasticity, yielding and episodicity in simple models of mantle convection

We explore the implications of refinements in the mechanical description of planetary constituents on the convection modes predicted by finite-element simulations. The refinements consist in the inclusion of incremental elasticity, plasticity (yielding) and multiple simultaneous creep mechanisms in addition to the usual visco-plastic models employed in the context of unified plate-mantle models. The main emphasis of this paper rests on the constitutive and computational formulation of the model. We apply a consistent incremental formulation of the non-linear governing equations avoiding the computationally expensive iterations that are otherwise necessary to handle the onset of plastic yield. In connection with episodic convection simulations, we point out the strong dependency of the results on the choice of the initial temperature distribution. Our results also indicate that the inclusion of elasticity in the constitutive relationships lowers the mechanical energy associated with subduction events.

Stirring in 3-d spherical models of convection in the Earth's mantle

On a global scale basalts from mid-ocean ridges are strikingly more homogeneous than basalts from intraplate volcanism. The observed geochemical heterogeneity argues strongly for the existence of distinct reservoirs in the Earth's mantle. It is an unresolved problem of Geodynamics as to how these findings can be reconciled with large-scale convection. We review observational constraints, and investigate stirring properties of numerical models of mantle convection. Conditions in the early Earth may have supported layered convection with rapid stirring in the upper layers. Material that has been altered near the surface is transported downwards by small-scale convection. Thereby a layer of homogeneous depleted material develops above pristine mantle. As the mantle cools over Earth history, the effects leading to layering become reduced and models show the large-scale convection favoured for the Earth today. Laterally averaged, the upper mantle below the lithosphere is least affected by material that has experienced near-surface differentiation. The geochemical signature obtained during the previous episode of small-scale convection may be preserved there for the longest time. Additionally, stirring is less effective in the high viscosity layer of the central lower mantle [1, 2], supporting the survival of medium-scale heterogeneities there. These models are the first, using 3-d spherical geometry and mostly Earth-like parameters, to address the suggested change of convective style. Although the models are still far from reproducing our planet, we find that proposal might be helpful towards reconciling geochemical and geophysical constraints.

An anisotropic viscous representation of Mohr-Coulomb failure for use is modelling coupled mantle-continent dynamics

In mantle convection models it has become common to make use of a modified (pressure sensitive, Boussinesq) von Mises yield criterion to limit the maximum stress the lithosphere can support. This approach allows the viscous, cool thermal boundary layer to deform in a relatively plate-like mode even in a fully Eulerian representation. In large-scale models with embedded continental crust where the mobile boundary layer represents the oceanic lithosphere, the von Mises yield criterion for the oceans ensures that the continents experience a realistic broad-scale stress regime. In detailed models of crustal deformation it is, however, more appropriate to choose a Mohr-Coulomb yield criterion based upon the idea that frictional slip occurs on whichever one of many randomly oriented planes happens to be favorably oriented with respect to the stress field. As coupled crust/mantle models become more sophisticated it is important to be able to use whichever failure model is appropriate to a given part of the system. We have therefore developed a way to represent Mohr-Coulomb failure within a code which is suited to mantle convection problems coupled to large-scale crustal deformation. Our approach uses an orthotropic viscous rheology (a different viscosity for pure shear to that for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation - neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane which is formally equivalent in any symmetric formulation. It is also important to allow strain-weakening of the material. The material should remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-Integration-point finite element code and will show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion, and comparisons between mantle convection models using the von Mises versus the Mohr-Coulomb yield criteria. The formulation itself is general and applies to 2D and 3D problems, although it is somewhat more complicated to identify the slip plane in 3D.

Anisotropy model for mantle convection

The paper presents a new theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including folding and kinking in multi-layered visco-elastic rock. The orientation of slip planes in the context of crystallographic slip is determined by the normal vector, the so-called director of these surfaces. The model is applied to simulate anisotropic natural mantle convection. We compare the evolution of the director and approximately steady states of isotropic and anisotropic convection. The isotropic case has a simple steady state solution, whereas the orthotropic convection model produces a continuously evolving patterning in tile core of the convection cell which makes only a near-steady condition possible, in which the thermal boundary layer appears to be well aligned with the flow and hence as observed in seismic tomomgraphy strong anistropic.

Geochemical and Sr, Nd, Pb isotope investigation of the New Caledonia ophiolite

The New Caledonia ophiolite hosts one of the largest obducted mantle section in the world, hence providing a unique insight for the study of upper mantle processes. These mantle rocks belong to an “atypical” ophiolitic sequence, which is dominated by refractory harzburgites but it also includes minor spinel and plagioclase lherzolites. Upper crust is notably absent in the ophiolite, with the exception of some mafic-ultramafic cumulates cropping out in the southern part of the island. Although the New Caledonia ophiolite has been under investigation for decades, its ultra-depleted nature has made its characterization an analytical challenge, so that few trace element data are available, while isotopic data are completely missing. In this thesis a comprehensive geochemical study (major, trace element and Sr-Nd-Pb isotopes) of the peridotites and the associated intrusive mafic rocks from the New Caledonia ophiolite has been carried out. The peridotites are low-strain tectonites showing porphyroclastic textures. Spinel lherzolites are undepleted lithotypes, as attested by the presence of 7-8 vol% of Na2O and Al2O3-rich clinopyroxene (up to 0.5 wt% Na2O; 6.5 wt% Al2O3), Fo content of olivine (88.5-90.0 mol%) and low Cr# of spinel (13-17). Conversely, harzburgites display a refractory nature, proven by the remarkable absence of primary clinopyroxene, very high Fo content in olivine (90.9-92.9 mol%), high Mg# in orthopyroxene (89.8-94.2) and Cr# in spinel (39-71). REE contents show abyssal-type patterns for spinel lherzolites, while harzburgites display U-shaped patterns, typical of fore-arc settings. Spinel lherzolites REE compositions are consistent with relatively low degree (8-9%) of fractional melting of a DMM source, starting in the garnet stability field. Conversely, REE models for harzburgites indicate high melting degrees (20-25%) of a DMM mantle source under spinel faies conditions, consistent with hydrous melting in forearc setting. Plagioclase lherzolites exhibit melt impregnation microtextures, Cr- and TiO2-enriched spinels and REE, Ti, Y, Zr progressive increase with respect to spinel lherzolites. Impregnation models indicate that plagioclase lherzolites may derive from spinel lherzolites by entrapment of highly depleted MORB melts in the shallow oceanic lithosphere. Mafic intrusives are olivine gabbronorites with a very refractory composition, as attested by high Fo content of olivine (87.3-88.9 mol.%), very high Mg# of clinopyroxene (87.7-92.2) and extreme anorthitic content of plagioclase (An = 90-96 mol%). The high Mg#, low TiO2 concentrations in pyroxenes and the anorthitic composition of plagioclase point out an origin from ultra-depleted primitive magmas in a convergent setting. Geochemical trace element models show that the parental melts of gabbronorites are primitive magmas with striking depleted compositions, bearing only in part similarities with the primitive boninitic melts of Bonin Islands. The first Sr, Nd and Pb isotope data obtained for the New Caledonia ophiolite highlight the presence of DM mantle source variably modified by different processes. Nd-Sr-Pb isotopic ratios for the lherzolites (+6.98≤epsilon Ndi≤+10.97) indicate a DM source that suffered low-temperature hydrothermal reactions. Harzburgites are characterized by a wide variation of Sr, Nd and Pb isotopic values, extending from DM-type to EM2 compositions (-0.82≤ epsilon Ndi≤+17.55), suggesting that harzburgite source was strongly affected by subduction-related processes. Conversely, combined trace element and Sr-Nd-Pb isotopic data for gabbronorites indicate a derivation from a source with composition similar to Indian-type mantle, but affected by fluid input in subduction environment. These geochemical features point out an evolution in a pre-Eocenic marginal basin setting, possibly in the proximity of a transform fault, for the lherzolites. Conversely, the harzburgites acquired their main geochemical and isotopic fingerprint in subduction zone setting.

Laminar heat transfer to non-Newtonian fluids in pipes with abrupt changes in diameter

In this thesis the results of experimental work performed to determine local heat transfer coefficients for non-Newtonian fluids in laminar flow through pipes with abrupt discontinuities are reported. The fluids investigated were water-based polymeric solutiorrs of time-indpendent, pseudoplastic materials, with flow indices "n" ranging from 0.39 to 0.9.The tube configurations were a 3.3 :1 sudden convergence, and a 1: 3.3 sudden divergence.The condition of a prescribed uniform wall heat flux was considered, with both upstream and downstream tube sections heated. Radial temperature traverses were also under­ taken primarily to justify the procedures used in estimating the tube wall and bulk fluid temperatures and secondly to give further insight into the mechanism of heat transfer beyond a sudden tube expansion. A theoretical assessment of the influence of viscous dissipation on a non-Newtonian pseudoplastic fluid of' arbitrary index "n" was carried out. The effects of other secondary factors such as free convection and temperature-dependent consistency were evaluated empirically. In the present investigations, the test conditions were chosen to minimise the effects of natural convection and the estimates of viscous heat generation showed the effect to be insignificant with the polymeric concentrations tested here. The final results have been presented as the relationships between local heat transfer coef'ficient and axial distance downstream of the discontinuities and relationships between dimensionless wall temperature and reduced radius. The influence of Reynolds number, Prandtl number, non-Newtonian index and heat flux have been indicated.