Analytical solution to the zero-inertia problem for surge flow phenomena in nonprismatic channels

Surge flow phenomena. e.g.. as a consequence of a dam failure or a flash flood, represent free boundary problems. ne extending computational domain together with the discontinuities involved renders their numerical solution a cumbersome procedure. This contribution proposes an analytical solution to the problem, It is based on the slightly modified zero-inertia (ZI) differential equations for nonprismatic channels and uses exclusively physical parameters. Employing the concept of a momentum-representative cross section of the moving water body together with a specific relationship for describing the cross sectional geometry leads, after considerable mathematical calculus. to the analytical solution. The hydrodynamic analytical model is free of numerical troubles, easy to run, computationally efficient. and fully satisfies the law of volume conservation. In a first test series, the hydrodynamic analytical ZI model compares very favorably with a full hydrodynamic numerical model in respect to published results of surge flow simulations in different types of prismatic channels. In order to extend these considerations to natural rivers, the accuracy of the analytical model in describing an irregular cross section is investigated and tested successfully. A sensitivity and error analysis reveals the important impact of the hydraulic radius on the velocity of the surge, and this underlines the importance of an adequate description of the topography, The new approach is finally applied to simulate a surge propagating down the irregularly shaped Isar Valley in the Bavarian Alps after a hypothetical dam failure. The straightforward and fully stable computation of the flood hydrograph along the Isar Valley clearly reflects the impact of the strongly varying topographic characteristics on the How phenomenon. Apart from treating surge flow phenomena as a whole, the analytical solution also offers a rigorous alternative to both (a) the approximate Whitham solution, for generating initial values, and (b) the rough volume balance techniques used to model the wave tip in numerical surge flow computations.

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.

Combined in-fermenter extraction and cross-flow microfiltration for improved inclusion body processing

In this study we demonstrate a new in-fermenter chemical extraction procedure that degrades the cell wall of Escherichia coli and releases inclusion bodies (IBs) into the fermentation medium. We then prove that cross-flow microfiltration can be used to remove 91% of soluble contaminants from the released IBs. The extraction protocol, based on a combination of Triton X-100, EDTA, and intracellular T7 lysozyme, effectively released most of the intracellular soluble content without solubilising the IBs. Cross-flow microfiltration using a 0.2 mum ceramic membrane successfully recovered the granulocyte macrophagecolony stimulating factor (GM-CSF) IBs with removal of 91% of the soluble contaminants and virtually no loss of IBs to the permeate. The filtration efficiency, in terms of both flux and transmission, was significantly enhanced by infermenter Benzonase(R) digestion of nucleic acids following chemical extraction. Both the extraction and filtration methods exerted their efficacy directly on a crude fermentation broth, eliminating the need for cell recovery and re-suspension in buffer. The processes demonstrated here can all be performed using just a fermenter and a single cross-flow filtration unit, demonstrating a high level of process intensification. Furthermore, there is considerable scope to also use the microfiltration system to subsequently solubilise the IBs, to separate the denatured protein from cell debris, and to refold the protein using diafiltration. In this way refolded protein can potentially be obtained, in a relatively pure state, using only two unit operations. (C) 2004 Wiley Periodicals Inc.

Continental extension: From core complexes to rigid block faulting

Extension of overthickened continental crust is commonly characterized by an early core complex stage of extension followed by a later stage of crustal-scale rigid block faulting. These two stages are clearly recognized during the extensional destruction of the Alpine orogen in northeast Corsica, where rigid block faulting overprinting core complex formation eventually led to crustal separation and the formation of a new oceanic backarc basin (the Ligurian Sea). Here we investigate the geodynamic evolution of continental extension by using a novel, fully coupled thermomechanical numerical model of the continental crust. We consider that the dynamic evolution is governed by fault weakening, which is generated by the evolution of the natural-state variables (i.e., pressure, deviatoric stress, temperature, and strain rate) and their associated energy fluxes. Our results show the appearance of a detachment layer that controls the initial separation of the brittle crust on characteristic listric faults, and a core complex formation that is exhuming strongly deformed rocks of the detachment zone and relatively undeformed crustal cores. This process is followed by a transitional period, characterized by an apparent tectonic quiescence, in which deformation is not localized and energy stored in the upper crust is transferred downward and causes self-organized mobilization of the lower crust. Eventually, the entire crust ruptures on major crosscutting faults, shifting the tectonic regime from core complex formation to wholesale rigid block faulting.

Driving mechanisms for groundwater flow and salt transport in a subterranean estuary

This paper presents field measurements and numerical simulations of groundwater dynamics in the intertidal zone of a sandy meso-tidal beach. The study, focusing on vertical hydraulic gradients and pore water salinities, reveals that tides and waves provide important forcing mechanisms for flow and salt transport in the nearshore aquifer. Such forcing, interacting with the beach morphology, enhances the exchange between the aquifer and ocean. The spatial and temporal variations of vertical hydraulic gradients demonstrate the complexity and dynamic nature of the processes and the extent of mixing between fresh groundwater and seawater in a subterranean estuary''. These results provide evidence of a potentially important reaction zone in the nearshore aquifer driven by oceanic oscillations. Land-derived contaminants may undergo important biogeochemical transformations in this zone prior to discharge.

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.

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.

Turbulent spots in boundary layers in a free-piston shock-tunnel flow

In this paper, experiments to detect turbulent spots in the transitional boundary layers, formed on a flat plate in a free-piston shock tunnel how, are reported. Experiments indicate that thin-film heat-transfer gauges are suitable for identifying turbulent-spot activity and can be used to identify parameters such as the convection rate of spots and the intermittency of turbulence.

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.

Low birth weight in response to salt restriction during pregnancy is not due to alterations in uterine-placental blood flow or the placental and peripheral renin-angiotensin system

A number of studies conducted in humans and in animals have observed that events occurring early in life are associated with the development of diseases in adulthood. Salt overload and restriction during pregnancy and lactation are responsible for functional (hemodynamic and hormonal) and structural alterations in adult offspring. Our group observed that lower birth weight and insulin resistance in adulthood is associated with salt restriction during pregnancy On the other hand, perinatal salt overload is associated with higher blood pressure and higher renal angiotensin II content in adult offspring. Therefore, we hypothesised that renin-angiotensin system (RAS) function is altered by changes in sodium intake during pregnancy. Such changes may influence fetoplacental blood flow and thereby fetal nutrient supply, with effects on growth in utero and, consequently, on birth weight. Female Wistar rats were fed low-salt (LS), normal-salt (NS), or high-salt (HS) diet, starting before conception and continuing until day 19 of pregnancy, Blood pressure, heart rate, fetuses and dams` body weight, placentae weight and litter size were measured on day 19 of pregnancy. Cardiac output, uterine and placental blood flow were also determined on day 19. Expressions of renin-angiotensin system components and of the TNF-alpha gene were evaluated in the placentae. Plasma renin activity (PRA) and plasma and tissue angiotensin-converting enzyme (ACE) activity, as well as plasma and placental levels of angiotensins I, II, and 1-7 were measured. Body weight and kidney mass were greater in HS than in NS and LS dams. Food intake did not differ among the maternal groups. Placental weight was lower in LS dams than in NS and HS dams. Fetal weight was lower in the US group than in the NS and HS groups. The PRA was greater in IS dams than in NS and HS dams, although ACE activity (serum, cardiac, renal, and placental) was unaffected by the level of sodium intake. Placental levels of angiotensins I and II were lower in the HS group than in the ISIS and IS groups. Placental angiotensin receptor type 1 (AT(1)) gene expression and levels of thiobarbituric acid reactive substances (TBARS) were higher in HS dams, as were uterine blood flow and cardiac output. The degree of salt intake did not influence plasma sodium, potassium or creatinine. Although fractional sodium excretion was higher in HS dams than in NS and LS dams, fractional potassium excretion was unchanged. In conclusion, findings from this study indicate that the reduction in fetal weight in response to salt restriction during pregnancy does not involve alterations in uterine-placental perfusion or the RAS. Moreover, no change in fetal weight is observed in response to salt overload during pregnancy. However, salt overload did lead to an increase in placental weight and uterine blood flow associated with alterations in maternal plasma and placental RAS. Therefore, these findings indicate that changes in salt intake during pregnancy lead to alterations in uterine-placental perfusion and fetal growth. (C) 2008 Elsevier Inc. All rights reserved.