Using three-dimensional microfluidic networks for solving computationally hard problems

This paper describes the design of a parallel algorithm that uses moving fluids in a three-dimensional microfluidic system to solve a nondeterministically polynomial complete problem (the maximal clique problem) in polynomial time. This algorithm relies on (i) parallel fabrication of the microfluidic system, (ii) parallel searching of all potential solutions by using fluid flow, and (iii) parallel optical readout of all solutions. This algorithm was implemented to solve the maximal clique problem for a simple graph with six vertices. The successful implementation of this algorithm to compute solutions for small-size graphs with fluids in microchannels is not useful, per se, but does suggest broader application for microfluidics in computation and control.

Illite and hydrocarbon exploration

Illite is a general term for the dioctahedral mica-like clay mineral common in sedimentary rocks, especially shales. Illite is of interest to the petroleum industry because it can provide a K-Ar isotope date that constrains the timing of basin heating events. It is critical to establish that hydrocarbon formation and migration occurred after the formation of the trap (anticline, etc.) that is to hold the oil. Illite also may precipitate in the pores of sandstone reservoirs, impeding fluid flow. Illite in shales is a mixture of detrital mica and its weathering products with diagenetic illite formed by reaction with pore fluids during burial. K-Ar ages are apparent ages of mixtures of detrital and diagenetic end members, and what we need are the ages of the end members themselves. This paper describes a methodology, based on mineralogy and crystallography, for interpreting the K-Ar ages from illites in sedimentary rocks and for estimating the ages of the end members.

What electrical measurements can say about changes in fault systems.

Earthquake zones in the upper crust are usually more conductive than the surrounding rocks, and electrical geophysical measurements can be used to map these zones. Magnetotelluric (MT) measurements across fault zones that are parallel to the coast and not too far away can also give some important information about the lower crustal zone. This is because the long-period electric currents coming from the ocean gradually leak into the mantle, but the lower crust is usually very resistive and very little leakage takes place. If a lower crustal zone is less resistive it will be a leakage zone, and this can be seen because the MT phase will change as the ocean currents leave the upper crust. The San Andreas Fault is parallel to the ocean boundary and close enough to have a lot of extra ocean currents crossing the zone. The Loma Prieta zone, after the earthquake, showed a lot of ocean electric current leakage, suggesting that the lower crust under the fault zone was much more conductive than normal. It is hard to believe that water, which is responsible for the conductivity, had time to get into the lower crustal zone, so it was probably always there, but not well connected. If this is true, then the poorly connected water would be at a pressure close to the rock pressure, and it may play a role in modifying the fluid pressure in the upper crust fault zone. We also have telluric measurements across the San Andreas Fault near Palmdale from 1979 to 1990, and beginning in 1985 we saw changes in the telluric signals on the fault zone and east of the fault zone compared with the signals west of the fault zone. These measurements were probably seeing a better connection of the lower crust fluids taking place, and this may result in a fluid flow from the lower crust to the upper crust. This could be a factor in changing the strength of the upper crust fault zone.

Modelo computacional para avaliação do desempenho hidrodinâmico de embarcações de planeio em águas calmas.

Em geral, uma embarcação de planeio é projetada para atingir elevados níveis de velocidade. Esse atributo de desempenho está diretamente relacionado ao porte da embarcação e à potência instalada em sua planta propulsiva. Tradicionalmente, durante o projeto de uma embarcação, as análises de desempenho são realizadas através de resultados de embarcações já existentes, retirados de séries sistemáticas ou de embarcações já desenvolvidas pelo estaleiro e/ou projetista. Além disso, a determinação dos atributos de desempenho pode ser feita através de métodos empíricos e/ou estatísticos, onde a embarcação é representada através de seus parâmetros geométricos principais; ou a partir de testes em modelos em escala reduzida ou protótipos. No caso específico de embarcações de planeio, o custo dos testes em escala reduzida é muito elevado em relação ao custo de projeto. Isso faz com que a maioria dos projetistas não opte por ensaios experimentais das novas embarcações em desenvolvimento. Ao longo dos últimos anos, o método de Savitsky foi largamente utilizado para se realizar estimativas de potência instalada de uma embarcação de planeio. Esse método utiliza um conjunto de equações semi-empíricas para determinar os esforços atuantes na embarcação, a partir dos quais é possível determinar a posição de equilíbrio de operação e a força propulsora necessária para navegar em uma dada velocidade. O método de Savitsky é muito utilizado nas fases iniciais de projeto, onde a geometria do casco ainda não foi totalmente definida, pois utiliza apenas as características geométricas principais da embarcação para realização das estimativas de esforços. À medida que se avança nas etapas de projeto, aumenta o detalhamento necessário das estimativas de desempenho. Para a realização, por exemplo, do projeto estrutural é necessária uma estimativa do campo de pressão atuante no fundo do casco, o qual não pode ser determinado pelo método de Savitsky. O método computacional implementado nesta dissertação, tem o objetivo de determinar as características do escoamento e o campo de pressão atuante no casco de uma embarcação de planeio navegando em águas calmas. O escoamento é determinado através de um problema de valor de contorno, no qual a superfície molhada no casco é considerada um corpo esbelto. Devido ao uso da teoria de corpo esbelto o problema pode ser tratado, separadamente, em cada seção, onde as condições de contorno são forçadamente respeitadas através de uma distribuição de vórtices.

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.