On the development of advanced techniques for mixed-elastohydrodynamic lubrification modelling of journal and sliding bearing systems.

The present thesis is focused on the development of a thorough mathematical modelling and computational solution framework aimed at the numerical simulation of journal and sliding bearing systems operating under a wide range of lubrication regimes (mixed, elastohydrodynamic and full film lubrication regimes) and working conditions (static, quasi-static and transient conditions). The fluid flow effects have been considered in terms of the Isothermal Generalized Equation of the Mechanics of the Viscous Thin Films (Reynolds equation), along with the massconserving p-Ø Elrod-Adams cavitation model that accordingly ensures the so-called JFO complementary boundary conditions for fluid film rupture. The variation of the lubricant rheological properties due to the viscous-pressure (Barus and Roelands equations), viscous-shear-thinning (Eyring and Carreau-Yasuda equations) and density-pressure (Dowson-Higginson equation) relationships have also been taken into account in the overall modelling. Generic models have been derived for the aforementioned bearing components in order to enable their applications in general multibody dynamic systems (MDS), and by including the effects of angular misalignments, superficial geometric defects (form/waviness deviations, EHL deformations, etc.) and axial motion. The bearing exibility (conformal EHL) has been incorporated by means of FEM model reduction (or condensation) techniques. The macroscopic in fluence of the mixedlubrication phenomena have been included into the modelling by the stochastic Patir and Cheng average ow model and the Greenwood-Williamson/Greenwood-Tripp formulations for rough contacts. Furthermore, a deterministic mixed-lubrication model with inter-asperity cavitation has also been proposed for full-scale simulations in the microscopic (roughness) level. According to the extensive mathematical modelling background established, three significant contributions have been accomplished. Firstly, a general numerical solution for the Reynolds lubrication equation with the mass-conserving p - Ø cavitation model has been developed based on the hybridtype Element-Based Finite Volume Method (EbFVM). This new solution scheme allows solving lubrication problems with complex geometries to be discretized by unstructured grids. The numerical method was validated in agreement with several example cases from the literature, and further used in numerical experiments to explore its exibility in coping with irregular meshes for reducing the number of nodes required in the solution of textured sliding bearings. Secondly, novel robust partitioned techniques, namely: Fixed Point Gauss-Seidel Method (PGMF), Point Gauss-Seidel Method with Aitken Acceleration (PGMA) and Interface Quasi-Newton Method with Inverse Jacobian from Least-Squares approximation (IQN-ILS), commonly adopted for solving uid-structure interaction problems have been introduced in the context of tribological simulations, particularly for the coupled calculation of dynamic conformal EHL contacts. The performance of such partitioned methods was evaluated according to simulations of dynamically loaded connecting-rod big-end bearings of both heavy-duty and high-speed engines. Finally, the proposed deterministic mixed-lubrication modelling was applied to investigate the in fluence of the cylinder liner wear after a 100h dynamometer engine test on the hydrodynamic pressure generation and friction of Twin-Land Oil Control Rings.

Interplay between cardiac function and heart development

Mechanotransduction refers to the conversion of mechanical forces into biochemical or electrical signals that initiate structural and functional remodeling in cells and tissues. The heart is a kinetic organ whose form changes considerably during development and disease. This requires cardiomyocytes to be mechanically durable and able to mount coordinated responses to a variety of environmental signals on different time scales, including cardiac pressure loading and electrical and hemodynamic forces. During physiological growth, myocytes, endocardial and epicardial cells have to adaptively remodel to these mechanical forces. Here we review some of the recent advances in the understanding of how mechanical forces influence cardiac development, with a focus on fluid flow forces. This article is part of a Special Issue entitled: Cardiomyocyte Biology: Integration of Develomental and Environmental Cues in the Heart edited by Marcus Schaub and Hughes Abriel.

Stable isotope ratios of interstitial fluids from ODP Leg 110 sites

Delta18O values of pore waters from the northern Barbados accretionary prism range from -0.3 to -3.6? and reflect pervasive reaction of volcanic ash to form smectite within the sedimentary sequence and continued low temperature alteration of basalt in the underlying ocean crust with the overprint of diffusive exchange between water in the sediment pores and the open ocean. Delta D values of pore waters in sediments sampled seaward of the deformation front drop from +5? at the sediment surface to -6? at the deepest levels sampled. These changes may also be related to alteration processes but remain largely enigmatic. Sediment deformation caused by impingement of the Caribbean plate on the Atlantic plate has instigated migration of chemically and isotopically distinct fluid along faults and coarse-grained sedimentary beds; delta18O values of pore waters are also locally affected by thrust stacking which increases diffusive pathlengths and possibly modifies diagenetic reaction rates in Pleistocene sediments. Migrating fluids are distinguished by anomalous delta18O values that are as much as 1? higher than those of surrounding fluids. Uncertainties in hydrogen isotope fractionation resulting from processes occurring under these conditions hinder identification of the hydrogen isotope composition of expelled fluid. Stable isotope analyses of pore waters help constrain the fluid migration history of the accretionary prism by limiting the source of fluids, the paths along which fluid flows, and the timing of faulting and subsequent fluid flow.

(Table T1) Physical properties of minicores from ODP Leg 193 sites, PACMANUS field

Permeability of the ocean crust is one of the most crucial parameters for constraining submarine fluid flow systems. Active hydrothermal fields are dynamic areas where fluid flow strongly affects the geochemistry and biology of the surrounding environment. There have been few permeability measurements in these regions, especially in felsic-hosted hydrothermal systems. We present a data set of 38 permeability and porosity measurements from the PACMANUS hydrothermal field, an actively venting, felsic hydrothermal field in the eastern Manus Basin. Permeability was measured using a complex transient method on 2.54-cm minicores. Permeability varies greatly between the samples, spanning over five orders of magnitude. Permeability decreases with both depth and decreasing porosity. When the alteration intensity of individual samples is considered, relationships between depth and porosity and permeability become more clearly defined. For incompletely altered samples (defined as >5% fresh rock), permeability and porosity are constant with depth. For completely altered samples (defined as <5% fresh rock), permeability and porosity decrease with depth. On average, the permeability values from the PACMANUS hydrothermal field are greater than those in other submarine environments using similar core-scale laboratory measurements; the average permeability, 4.5 x 10-16 m**2, is two to four orders of magnitude greater than in other areas. Although the core-scale permeability is higher than in other seafloor environments, it is still too low to obtain the fluid velocities observed in the PACMANUS hydrothermal field based on simplified analytical calculations. It is likely that core-scale permeability measurements are not representative of bulk rock permeability of the hydrothermal system overall, and that the latter is predominantly fracture controlled.

(Tables 1, 3) Component composition of organic particles in samples from ODP Site 131-808

Analysis of the palynofacies and miospore thermal alteration indices (TAI) of sediments from ODP Site 808 in the Nankai Trough was undertaken to determine (1) the source, depositional environment, and diagenesis of organic matter in the accreted sediments, and (2) the thermal structure and history of the prism and its relationship to fluid flow. Using the Hartax classification system, two palynofacies were recognized in the sedimentary sequence. Facies 1 occurs within the upper 600 m of trench-wedge turbidites (sedimentation rate > 1 km/m.y.) and contains >50% inertite particles. The rest of the assemblage is dominated by well-preserved phytoclasts and contains small amounts of poorly preserved phytoclasts and well-preserved scleratoclasts. Facies 2 occurs within the Shikoku Basin hemipelagites (600-1300 m below seafloor; sedimentation rate <150 m/m.y.) and contains over two-thirds inertite particles. The rest of the assemblage is dominated by poorly preserved phytoclasts. Miospores and marine phytoplankton compose only a small percentage of both palynofacies. Degraded organic matter is most noticeable in Facies 2, whereas its presence in Facies 1 is overshadowed by the high influx of well-preserved primary organic matter. Most of the degraded organic matter and inertite is interpreted to be reworked. Some of the degraded organic matter may be primary, and may have experienced more biodegradation and thermal alteration in Facies 2 than in Facies 1. TAI values indicate an immature stage of organic maturation (< 2) down to about 900 mbsf. Below this, samples show an increase with depth to a mature stage, reaching peak levels of about 3 just above basement. Samples from within the thrust fault and decollement zones do not show levels of maturity significantly greater than those of surrounding samples, leaving uncertain whether hot fluids have migrated along these fault boundaries in the past.

Isotopic composition of carbonates from the Logachev Hydrothermal Field, MAR 14°45'N

We investigate the Logatchev Hydrothermal Field at the Mid-Atlantic Ridge, 14°45' N to constrain the calcium isotope hydrothermal flux into the ocean. During the transformation of seawater to a hydrothermal solution, the Ca concentration of pristine seawater ([Ca]_SW) increases from about 10 mM to about 32 mM in the hydrothermal fluid endmember ([Ca]_HydEnd) and thereby adopts a d44/40Ca_HydEnd of -0.95+/-0.07 per mil relative to seawater (SW) and a 87Sr/86Sr isotope ratio of 0.7034(4). We demonstrate that d44/40Ca_HydEnd is higher than that of the bedrock at the Logatchev field. From mass balance calculations, we deduce a d44/40Ca of -1.17+/-0.04 per mil (SW) for the host-rocks in the reaction zone and -1.45+/-0.05 per mil (SW) for the isotopic composition of the entire hydrothermal cell of the Logatchev field. The values are isotopically lighter than the currently assumed d44/40Ca for Bulk Earth of -0.92+/-0.18 per mil (SW) [Skulan J., DePaolo D. J. and Owens T. L. (1997) Biological control of calcium isotopic abundances in the global calcium cycle. Geochim. Cosmochim. Acta 61,(12) 2505-2510] and challenge previous assumptions of no Ca isotope fractionation between hydrothermal fluid and the oceanic crust [Zhu P. and Macdougall J. D. (1998) Calcium isotopes in the marine environment and the oceanic calcium cycle. Geochim. Cosmochim. Acta 62,(10) 1691-1698; Schmitt A. -D., Chabeaux F. and Stille P. (2003) The calcium riverine and hydrothermal isotopic fluxes and the oceanic calcium mass balance. Earth Planet. Sci. Lett. 6731, 1-16]. Here we propose that Ca isotope fractionation along the fluid flow pathway of the Logatchev field occurs during the precipitation of anhydrite. Two anhydrite samples from the Logatchev Hydrothermal Field show an average fractionation of about D44/40Ca = -0.5 per mil relative to their assumed parental solutions. Ca isotope ratios in aragonites from carbonate veins from ODP drill cores indicate aragonite precipitation directly from seawater at low temperatures with an average d44/40Ca of -1.54+/-0.08 per mil (SW). The relatively large fractionation between the aragonite precipitates and seawater in combination with their frequent abundance in weathered mafic and ultramafic rocks suggest a reconsideration of the marine Ca isotope budget, in particular with regard to ocean crust alteration.

Depth profiles of dissolved sulfide and sulfate in the pore water of hydrocarbon seep sediments offshore Pakistan

The interaction between fluid seepage, bottom water redox, and chemosynthetic communities was studied at cold seeps across one of the world's largest oxygen minimum zones (OMZ) located at the Makran convergent continental margin. Push cores were obtained from seeps within and below the core-OMZ with a remotely operated vehicle. Extracted sediment pore water was analyzed for sulfide and sulfate concentrations. Depending on oxygen availability in the bottom water, seeps were either colonized by microbial mats or by mats and macrofauna. The latter, including ampharetid polychaetes and vesicomyid clams, occurred in distinct benthic habitats, which were arranged in a concentric fashion around gas orifices. At most sites colonized by microbial mats, hydrogen sulfide was exported into the bottom water. Where macrofauna was widely abundant, hydrogen sulfide was retained within the sediment. Numerical modeling of pore water profiles was performed in order to assess rates of fluid advection and bioirrigation. While the magnitude of upward fluid flow decreased from 11 cm yr**-1 to <1 cm yr**-1 and the sulfate/methane transition (SMT) deepened with increasing distance from the central gas orifice, the fluxes of sulfate into the SMT did not significantly differ (6.6-9.3 mol m**-2 yr**-1). Depth-integrated rates of bioirrigation increased from 120 cm yr**-1 in the central habitat, characterized by microbial mats and sparse macrofauna, to 297 cm yr**-1 in the habitat of large and few small vesicomyid clams. These results reveal that chemosynthetic macrofauna inhabiting the outer seep habitats below the core-OMZ efficiently bioirrigate and thus transport sulfate down into the upper 10 to 15 cm of the sediment. In this way the animals deal with the lower upward flux of methane in outer habitats by stimulating rates of anaerobic oxidation of methane (AOM) with sulfate high enough to provide hydrogen sulfide for chemosynthesis. Through bioirrigation, macrofauna engineer their geochemical environment and fuel upward sulfide flux via AOM. Furthermore, due to the introduction of oxygenated bottom water into the sediment via bioirrigation, the depth of the sulfide sink gradually deepens towards outer habitats. We therefore suggest that - in addition to the oxygen levels in the water column, which determine whether macrofaunal communities can develop or not - it is the depth of the SMT and thus of sulfide production that determines which chemosynthetic communities are able to exploit the sulfide at depth. We hypothesize that large vesicomyid clams, by efficiently expanding the sulfate zone down into the sediment, could cut off smaller or less mobile organisms, as e.g. small clams and sulfur bacteria, from the sulfide source.

(Table 1) Sr-isotope values of representative sediment samples from ODP sites at the Costa Rica margin

Two sealed borehole hydrologic observatories (CORKs) were installed in two active hydrogeochemical systems at the Costa Rica subduction zone to investigate the relationship between tectonics, fluid flow, and fluid composition. The observatories were deployed during Ocean Drilling Program (ODP) Leg 205 at Site 1253, ~ 0.2 km seaward of the trench, in the upper igneous basement, and at Site 1255, ~ 0.5 km landward of the trench, in the décollement. Downhole instrumentation was designed to monitor formation fluid flow rates, composition, pressure, and temperature. The two-year records collected by this interdisciplinary effort constitute the first co-registered hydrological, chemical, and physical dataset from a subduction zone, providing critical information on the average and transient state of the subduction thrust and upper igneous basement. The continuous records at ODP Site 1253 show that the uppermost igneous basement is highly permeable hosting an average fluid flow rate of 0.3 m/yr, and indicate that the fluid sampled in the basement is a mixture between seawater (~ 50%) and a subduction zone fluid originating within the forearc (~ 50%). These results suggest that the uppermost basement serves as an efficient pathway for fluid expelled from the forearc that should be considered in models of subduction zone hydrogeology and deformation. Three transients in fluid flow rates were observed along the décollement at ODP Site 1255, two of which coincided with stepwise increases in formation pressure. These two transients are the result of aseismic slip dislocations that propagated up-dip from the seismogenic zone over the course of ~ 2 weeks terminating before reaching ODP Site 1255 and the trench. The nature and temporal behavior of strain and the associated hydrological response during these slow slip events may be an analog for the response of the seaward part of the subduction prism during or soon after large subduction zone earthquakes.

(Table 1) Deuterium ratios, chloride content and hydrate volume of pore water and hydrate meltwater from ODP Site 164-996

Site 996 is located above the Blake Diapir where numerous indications of vertical fluid migration and the presence of hydrate existed prior to Ocean Drilling Program (ODP) Leg 164. Direct sampling of hydrates and visual observations of hydrate-filled veins that could be traced 30-40 cm along cores suggest a connection between fluid migration and hydrate formation. The composition of pore water squeezed from sediment cores showed large variations due to melting of hydrate during core recovery and influence of saline water from the evaporitic diapir below. Analysis of water released during hydrate decomposition experiments showed that the recovered hydrates contained significant amounts of pore water. Solutions of the transport equations for deuterium (d2H) and chloride (Cl-) were used to determine maximum (d2H) and minimum (Cl-) in situ concentrations of these species. Minimum in situ concentrations of hydrate were estimated by combining these results with Cl- and d2H values measured on hydrate meltwaters and pore waters obtained by squeezing of sediments, by the means of a method based on analysis of distances in the two-dimensional Cl- d2H space. The computed Cl- and d2H distribution indicates that the minimum hydrate amount solutions are representative of the actual hydrate amount. The highest and mean hydrate concentrations estimates from our model are 31% and 10% of the pore space, respectively. These concentrations agree well with visual core observations, supporting the validity of the model assumptions. The minimum in situ Cl- concentrations were used to constrain the rates of upward fluid migration. Simulation of all available data gave a mean flow rate of 0.35 m/k.y. (range: 0.125-0.5 m/k.y.).

Interferometric measurement in a shock tube of dissociation rates for air and its component gases.

At head of title: Cornell University, Graduate School of Aeronautical Engineering.

Chronic shallow seawater circulation driven by subsurface gas dynamics, at the Southern Hydrate Ridge seep system, Cascadia

Senior thesis written for Oceanography 445

Rare earth element fractionation in authigenic illite–smectite from Late Permian clastic rocks, Bowen Basin, Australia: implications for physico-chemical environments of fluids during illitization

REE analyses were performed on authigenic illitic clay. minerals from Late Permian mudrocks, sandstones and bentonites from the Bowen Basin (Australia). The mixed-layer illite-smectite exhibit REE patterns with an obvious fractionation of the HREE from the LREE and MREE, which is an apparent function of degree of illitization reaction. The highly illitic (R greater than or equal to 3) illite-smectite from the northern Bowen Basin show a depletion of LREE relative to the less illitic (R=0 and 1) clays. In contrast, an enrichment of HREE for the illite-rich clays relative to less. illitic clays is evident for the southern Bowen Basin samples. The North American Shale Composite-normalized (La/Lu)(sn) ratios show negative correlations with the illite content in illite-smectite and positive correlations with the delta(18)O values of the clays for both the northern and southern Bowen Basin samples. These correlations indicate that the increasing depletion of LREE in hydrothermal fluids is a function of increasing water/rock ratios in the northern Bowen Basin. Good negative correlations between (La/Lu)(sn) ratios and illite content in illite-smectite from the southern Bowen Basin suggest the involvement of fluids with higher alkalinity and higher pH in low water/ rock ratio conditions. Increasing HREE enrichment with delta(18)O decrease indicates the effect of increasing temperature at low water/rock ratios in the southern Bowen Basin. Results of the present study confirm the conclusions of some earlier studies suggesting that REE in illitic clay minerals are mobile and fractionated during illitization and that this fact should be considered in studies of sedimentary processes and in identifying provenance. Moreover, our results show that REE systematic of illitic clay minerals can be applied as an useful technique to gain information about physico-chemical conditions during thermal and fluid flow events in certain sedimentary basins. (C) 2003 Elsevier Science B.V. All rights reserved.

Three-dimensional simulation of the rolling of a saturated fibro-porous media

This paper describes recent advances made in computational modelling of the sugar cane liquid extraction process. The saturated fibro-porous material is rolled between circumferentially grooved rolls, which enhance frictional grip and provide a low-resistance path for liquid flow during the extraction process. Previously reported two-dimensional (2D) computational models, account for the large deformation of the porous material by solving the fully coupled governing fibre stress and fluid-flow equations using finite element techniques. While the 2D simulations provide much insight into the overarching cause-effect relationships, predictions of mechanical quantities such as roll separating force and particularly torque as a function of roll speed and degree of compression are not satisfactory for industrial use. It is considered that the unsatisfactory response in roll torque prediction may be due to the stress levels that exist between the groove tips and roots which have been largely neglected in the geometrically simplified 2D model. This paper gives results for both two- and three-dimensional finite element models and highlights their strengths and weaknesses in predicting key milling parameters. (c) 2005 Elsevier B.V. All rights reserved.