Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

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.

Effect of material anisotropy on the onset of convective flow in three-dimensional fluid-saturated faults

In order to investigate the effect of material anisotropy on convective instability of three-dimensional fluid-saturated faults, an exact analytical solution for the critical Rayleigh number of three-dimensional convective flow has been obtained. Using this critical Rayleigh number, effects of different permeability ratios and thermal conductivity ratios on convective instability of a vertically oriented three-dimensional fault have been examined in detail. It has been recognized that (1) if the fault material is isotropic in the horizontal direction, the horizontal to vertical permeability ratio has a significant effect on the critical Rayleigh number of the three-dimensional fault system, but the horizontal to vertical thermal conductivity ratio has little influence on the convective instability of the system, and (2) if the fault material is isotropic in the fault plane, the thermal conductivity ratio of the fault normal to plane has a considerable effect on the critical Rayleigh number of the three-dimensional fault system, but the effect of the permeability ratio of the fault normal to plane on the critical Rayleigh number of three-dimensional convective flow is negligible.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

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.

Thermal Characteristics of Polybia scutellaris Nests (Hymenoptera: Vespidae) Using Computational Fluid Dynamics: A Possible Adaptation to Tropical Climates

Polybia scutellaris constructs huge nests characterized by numerous spinal projections on the surface. We investigated the thermal characteristics of P scutellaris nests in order to determine whether their nest temperature is homeothermically maintained and whether the spines play a role in the thermoregulation of the nests. In order to examine these hypotheses, we measured the nest temperature in a active nest and in an abandoned nest. The temperature in the active nest was almost stable at 27 degrees C, whereas that of the abandoned nest varied with changes in the ambient temperature, suggesting that nest temperature was maintained by the thermogenesis of colony individuals. In order to predict the thermal properties of the spines, a numerical simulation was employed. To construct a 3D-model of a P scutellaris nest, the nest architecture was simplified into an outer envelope and the surface spines, for both of which the initial temperature was set at 27 degrees C. The physical properties of the simulated nest were regarded to be those of wood since the nest of this species is constructed from plant materials. When the model was exposed to cool air (12 degrees C), the temperature was lower in the models with more spines. On the other hand, when the nest was heated (42 degrees C), the temperature increase was smaller in models with more spines. It is suggested that the spines act as a heat radiator, not as an insulator, against the changes in ambient temperature.

Fast-convergent iterative scheme for filtering velocity signals and finding Kolmogorov scales

The present fundamental knowledge of fluid turbulence has been established primarily from hot- and cold-wire measurements. Unfortunately, however, these measurements necessarily suffer from contamination by noise since no certain method has previously been available to optimally filter noise from the measured signals. This limitation has impeded our progress of understanding turbulence profoundly. We address this limitation by presenting a simple, fast-convergent iterative scheme to digitally filter signals optimally and find Kolmogorov scales definitely. The great efficacy of the scheme is demonstrated by its application to the instantaneous velocity measured in a turbulent jet.

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.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Online monitoring of pulse pressure variation to guide fluid therapy after cardiac surgery

BACKGROUND: The arterial pulse pressure variation induced by mechanical ventilation (Delta PP) has been shown to be a predictor of fluid responsiveness. Until now, Delta PP has had to be calculated offline (from a computer recording or a paper printing of the arterial pressure curve), or to be derived from specific cardiac output monitors, limiting the widespread use of this parameter. Recently, a method has been developed for the automatic calculation and real-time monitoring of Delta PP using standard bedside monitors. Whether this method is to predict reliable predictor of fluid responsiveness remains to be determined. METHODS: We conducted a prospective clinical study in 59 mechanically ventilated patients in the postoperative period of cardiac surgery. Patients studied were considered at low risk for complications related to fluid administration (pulmonary artery occlusion pressure <20 mm Hg, left ventricular ejection fraction >= 40%). All patients were instrumented with an arterial line and a pulmonary artery catheter. Cardiac filling pressures and cardiac output were measured before and after intravascular fluid administration (20 mL/kg of lactated Ringer`s solution over 20 min), whereas Delta PP was automatically calculated and continuously monitored. RESULTS: Fluid administration increased cardiac output by at least 15% in 39 patients (66% = responders). Before fluid administration, responders and nonresponders were comparable with regard to right atrial and pulmonary artery occlusion pressures. In contrast, Delta PP was significantly greater in responders than in nonresponders, (17% +/- 3% vs 9% +/- 2%, P < 0.001). The Delta PP cut-off value of 12% allowed identification of responders with a sensitivity of 97% and a specificity of 95%. CONCLUSION: Automatic real-time monitoring of Delta PP is possible using a standard bedside rnonitor and was found to be a reliable method to predict fluid responsiveness after cardiac surgery. Additional studies are needed to determine if this technique can be used to avoid the complications of fluid administration in high-risk patients.

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.

Changes in purines concentration in the cerebrospinal fluid of patients experiencing pain: A case-control study

This study analyzes the relationship between extracellular purines and pain perception in humans. Cerebrospinal fluid (CSF) levels of purines and their metabolites were compared between patients displaying acute and/or chronic pain syndromes and control subjects. The CSF levels of IMP, inosine, guanosine and uric acid were significantly increased in the chronic pain group and correlated with pain severity (P<0.05). Patients displaying both chronic and acute pain presented similar changes in the CSF purines concentration (P<0.05). However, in the acute pain group, only CSF inosine and uric acid levels were significantly increased (P<0.05). These findings suggest that purines, in special inosine, guanosine and uric acid, are associated with the spinal mechanisms underlying nociception. (C) 2010 Elsevier Ireland Ltd. All rights reserved.

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.