Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Double diffusion-driven convective instability of three-dimensional fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the double diffusion-driven convective instability of three-dimensional fluid-saturated geological fault zones when they are heated uniformly from below. The fault zone is assumed to be more permeable than its surrounding rocks. In particular, we have derived exact analytical solutions to the total critical Rayleigh numbers of the double diffusion-driven convective flow. Using the corresponding total critical Rayleigh numbers, the double diffusion-driven convective instability of a fluid-saturated three-dimensional geological fault zone system has been investigated. The related theoretical analysis demonstrates that: (1) The relative higher concentration of the chemical species at the top of the three-dimensional geological fault zone system can destabilize the convective flow of the system, while the relative lower concentration of the chemical species at the top of the three-dimensional geological fault zone system can stabilize the convective flow of the system. (2) The double diffusion-driven convective flow modes of the three-dimensional geological fault zone system are very close each other and therefore, the system may have the similar chance to pick up different double diffusion-driven convective flow modes, especially in the case of the fault thickness to height ratio approaching 0. (3) The significant influence of the chemical species diffusion on the convective instability of the three-dimensional geological fault zone system implies that the seawater intrusion into the surface of the Earth is a potential mechanism to trigger the convective flow in the shallow three-dimensional geological fault zone system.

Characterization of optically driven fluid stress fields with optical tweezers

We present a controlled stress microviscometer with applications to complex fluids. It generates and measures microscopic fluid velocity fields, based on dual beam optical tweezers. This allows an investigation of bulk viscous properties and local inhomogeneities at the probe particle surface. The accuracy of the method is demonstrated in water. In a complex fluid model (hyaluronic acid), we observe a strong deviation of the flow field from classical behavior. Knowledge of the deviation together with an optical torque measurement is used to determine the bulk viscosity. Furthermore, we model the observed deviation and derive microscopic parameters.

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

Arterial pulse pressure variation predicting fluid responsiveness in critically ill patients (Retracted article. See vol. 31, pg. 542, 2009)

In critically ill patients, it is important to predict which patients will have their systemic blood flow increased in response to volume expansion to avoid undesired hypovolemia and fluid overloading. Static parameters such as the central venous pressure, the pulmonary arterial occlusion pressure, and the left ventricular end-diastolic dimension cannot accurately discriminate between responders and nonresponders to a fluid challenge. In this regard, respiratory-induced changes in arterial pulse pressure have been demonstrated to accurately predict preload responsiveness in mechanically ventilated patients. Some experimental and clinical studies confirm the usefulness of arterial pulse pressure as a useful tool to guide fluid therapy in critically ill patients.

Fluid Resuscitation With Isotonic or Hypertonic Saline Solution Avoids Intraneural Calcium Influx After Traumatic Brain Injury Associated With Hemorrhagic Shock

Background: Calcium is one of the triggers involved in ischemic neuronal death. Because hypotension is a strong predictor of outcome in traumatic brain injury (TBI), we tested the hypothesis that early fluid resuscitation blunts calcium influx in hemorrhagic shock associated to TBI. Methods: Fifteen ketamine-halothane anesthetized mongrel dogs (18.7 kg +/- 1.4 kg) underwent unilateral cryogenic brain injury. Blood was shed in 5 minutes to a target mean arterial pressure of 40 mm Hg to 45 mm Hg and maintained at these levels for 20 minutes (shed blood volume = 26 mL/kg +/- 7 mL/kg). Animals were then randomized into three groups: CT (controls, no fluid resuscitation), HS (7.5% NaCl, 4 mL/kg, in 5 minutes), and LR (lactate Ringer`s, 33 mL/kg, in 15 minutes). Twenty minutes later, a craniotomy was performed and cerebral biopsies were obtained next to the lesion (""clinical penumbra"") and from the corresponding contralateral side (""lesion`s mirror"") to determine intracellular calcium by fluorescence signals of Fura-2-loaded cells. Results: Controls remained hypotensive and in a low-flow state, whereas fluid resuscitation improved hemodynamic profile. There was a significant increase in intracellular calcium in the injured hemisphere in CT (1035 nM +/- 782 nM), compared with both HS (457 nM +/- 149 nM, p = 0.028) and LR (392 nM +/- 178 nM, p = 0.017), with no differences between HS and LR (p = 0.38). Intracellular calcium at the contralateral, uninjured hemisphere was 438 nM +/- 192 nM in CT, 510 nM +/- 196 nM in HS, and 311 nM +/- 51 nM in LR, with no significant differences between them. Conclusion: Both small volume hypertonic saline and large volume lactated Ringer`s blunts calcium influx in early stages of TBI associated to hemorrhagic shock. No fluid resuscitation strategy promotes calcium influx and further neural damage.

Simulation of mixing in unsteady flow through a periodically obstructed channel

Fluid mixing in steady and unsteady Bow through a channel containing periodic square obstructions has been studied using a finite-difference simulation to determine fluid velocities, followed by the use of passive marker particle advection to look at fluid transport out of the cavities formed between each of the obstructions. The geometry and Bow conditions were chosen from the work by Perkins (1989, M.S. Thesis, Lehigh University; 1992, Ph.D. Thesis, Lehigh University); who investigated heat transfer enhancement due to unsteady flow through such an obstructed channel. Particle advection shows that Bow regimes which are predicted to give good mixing based on snapshots of instantaneous streamline contour plots were not necessarily able to efficiently mix fluid which started in the cavity regions throughout the channel. The use of Poincare sections shows regular regions existing under these conditions which inhibit efficient fluid transport. These regular regions are found to disappear when the unsteady Bow velocity is increased. (C) 1997 Elsevier Science Ltd.

Effects of Fluid Resuscitation on Cardiovascular Performance After Posttraumatic Pneumonectomy

Background: Several factors have been implicated in the high-mortality rate of posttraumatic pneumonectomy. In this study, we evaluated the hemodynamic and echocardiographic changes induced by pneumonectomy and fluid resuscitation after hemorrhagic shock. Methods: Fourteen dogs were bled to a target mean arterial pressure of 40 mmHg. The animals were assigned to two groups: control (no fluid resuscitation) and lactated Ringer`s (3 x shed blood volume). The left pulmonary hilum was cross clamped, and the animals were observed for 60 minutes. Systemic hemodynamics was evaluated using Swan-Ganz, arterial catheter, and ultrasonic flow probe. Systemic O(2)-derived variables were calculated. Ejection fraction was determined by two-dimensional echocardiography. Results: Fluid resuscitation improved the mean arterial pressure and systemic oxygen delivery. After pneumonectomy, no significant increase in right ventricular pressure was observed in the LR group. No signs of major ventricular dilation or changes in arterial oxygenation were observed. Conclusion: Our data suggest that pneumonectomy is not associated with early pulmonary hypertension; gentle fluid resuscitation improves cardiovascular performance and is not associated with an increase in right ventricular pressure.

Fluid Replacement With Hypertonic or Isotonic Solutions Guided by Mixed Venous Oxygen Saturation in Experimental Hypodynamic Sepsis

Background: Splanchnic perfusion is prone to early injury and persists despite normalization of global hemodynamic variables in sepsis. Volume replacement guided by oxygen derived variables has been recommended in the management of septic patients. Our hypothesis was that a hypertonic isoneotic solution Would improve the benefits of crystalloids replacement guided by mixed venous oxygen saturation. Methods: Seventeen anesthetized and mechanically ventilated mongrel dogs received an intravenous infusion of live E. coli in 30 minutes. They were then randomized into three groups: control group (n = 3) bacterial infusion without treatment; normal saline (n = 7), initial fluid replacement with 32 mL/kg of normal saline during 20 minutes; hypertonic solution (n = 7), initial fluid replacement with 4 mL/kg of hypertonic solution during 5 minutes. After 30 and 60 Minutes, additional boluses of normal saline were administered when mixed venous oxygen saturation remained below 70%. Mean arterial pressure, cardiac output; regional blood flows, systemic and regional oxygen-derived variables, and lactate levels were assessed. Animals were observed for 90 minutes and then killed. Hystopathological analysis including apoptosis detection using terminal deoxynucleotidil transferase mediated dUTP-biotin nick end labeling was performed. Results: A hypodynamic septic shock was observed after bacterial infusion. Both the fluid-treated groups presented similar transient benefits in systemic and regional variables. A greater degree of gut epithelial cells apoptosis was observed in normal saline-treated animals. Conclusions: Although normalization of mixed venous oxygen saturation was not associated with restoration of markers of splanchnic or other systemic perfusion variables, the initial fluid savings with hypertonic saline and its latter effect on gut apoptosis may be of interest in sepsis management.

Can eccentric arterial plaques alone cause flow stagnation points and favour thrombus incorporation?

We have used an experimental model of aorta stenosis, with a Plexiglas plug, simulating a stable atheromatous plaque that promotes local turbulence and thrombosis. With animal survival of more than 24 h, we followed the partial fibrinolysis of the thrombus as well as its posterior organization and incorporation to the arterial wall as a neointima for up to 30 days. The mushroom plug form permitted the development of recirculation and stasis areas around it, favouring this evolution. Despite noted limitations, this study demonstrates that thrombus incorporation can contribute to plaque extension, as it can promote recirculation and stasis areas.

A small disturbance in the strong vortex flow

A small disturbance in the axisymmetric, bathtub-like flow with strong vorticity is considered and the asymptotic representation of the solution is found. It is shown that if the disturbance is smaller than a certain critical scale, the conventional strong vortex approximation cannot describe the field generated by the disturbance not only in the vicinity of the disturbance but also at the distances much larger than the critical scale. (C) 2001 American Institute of Physics.

Moderately strong vorticity in a bathtub-type flow

Theoretical and numerical analysis is performed for an inviscid axisymmetric vortical bathtub-type flow. The level of vorticity is kept high so that the image of the flow on the radial-axial plane (r-z plane) is not potential. The most significant findings are: (1) the region of validity of the strong vortex approximation is separated from the drain by a buffer region, (2) the power-law asymptote of the stream function, specified by Delta Psi similar to r(4/3) Deltaz, appears near the axis when vorticity in the flow is sufficiently strong and (3) the local Rossby number in the region of the 4/3 power-law the initial vorticity level in the flow and the global Rossby number.