Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Heat flow for Yang-Mills-Higgs fields, part I

The Yang-Mills-Higgs field generalizes the Yang-Mills field. The authors establish the local existence and uniqueness of the weak solution to the heat flow for the Yang-Mills-Higgs field in a vector bundle over a compact Riemannian 4-manifold, and show that the weak solution is gauge-equivalent to a smooth solution and there are at most finite singularities at the maximum existing time.

Heat flow for the Yang-Mills-Higgs field and the Hermitian Yang-Mills-Higgs metric

For a parameter lambda > 0, we study a type of vortex equations, which generalize the well-known Hermitian-Einstein equation, for a connection A and a section phi of a holomorphic vector bundle E over a Kahler manifold X. We establish a global existence of smooth solutions to heat flow for a self-dual Yang-Mills-Higgs field on E. Assuming the lambda -stability of (E, phi), we prove the existence of the Hermitian Yang-Mills-Higgs metric on the holomorphic bundle E by studying the limiting behaviour of the gauge flow.

Similitude applied to centrifugal scaling of unsaturated flow

Centrifuge experiments modeling single-phase flow in prototype porous media typically use the same porous medium and permeant. Then, well-known scaling laws are used to transfer the results to the prototype. More general scaling laws that relax these restrictions are presented. For permeants that are immiscible with an accompanying gas phase, model-prototype (i.e., centrifuge model experiment-target system) scaling is demonstrated. Scaling is shown to be feasible for Miller-similar (or geometrically similar) media. Scalings are presented for a more, general class, Lisle-similar media, based on the equivalence mapping of Richards' equation onto itself. Whereas model-prototype scaling of Miller-similar media can be realized easily for arbitrary boundary conditions, Lisle-similarity in a finite length medium generally, but not always, involves a mapping to a moving boundary problem. An exception occurs for redistribution in Lisle-similar porous media, which is shown to map to spatially fixed boundary conditions. Complete model-prototype scalings for this example are derived.

Model for gas-liquid flow through wet porous medium

A method involving bubbling of air through a fibrous filter immersed in water has recently been investigated (Agranovski et al. [1]). Experimental results showed that the removal efficiency for ultra-fine aerosols by such filters was greatly increased compared to dry filters. Nuclear Magnetic Resonance (NMR) imaging was used to examine the wet filter and to determine the nature of the gas flow inside the filter (Agranovski et al. [2]). It was found that tortuous preferential pathways (or flow tubes) develop within the filter through which the air flows and the distribution of air and water inside the porous medium has been investigated. The aim of this paper is to investigate the geometry of the pathways and to make estimates of the flow velocities and particle removal efficiency in such pathways. A mathematical model of the flow of air along the preferred pathways has been developed and verified experimentally. Even for the highest realistic gas velocity the flow field was essentially laminar (Re approximate to 250). We solved Laplace's equation for stream function to map trajectories of particles and gas molecules to investigate the possibility of their removal from the carrier.

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.

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.

Determination of fluoxetine and norfluoxetine enantiomers in human plasma by polypyrrole-coated capillary in-tube solid-phase microextraction coupled with liquid chromatography-fluorescence detection

A sensitive, selective, and reproducible in-tube polypyrrole-coated capillary (PPY) solid-phase microextraction and liquid chromatographic method for fluoxetine and norfluoxetine enantiomers analysis in plasma samples has been developed, validated, and further applied to the analysis of plasma samples from elderly patients undergoing therapy with antidepressants. Important factors in the optimization of in-tube SPME efficiency are discussed, including the sample draw/eject volume, draw/eject cycle number, draw/eject flow-rate, sample pH, and influence of plasma proteins. Separation of the analytes was achieved with a Chiralcel OD-R column and a mobile phase consisting of potassium hexafluorophosphate 7.5 mM and sodium phosphate 0.25 M solution, pH 3.0, and acetonitrile (75:25, v/v) in the isocratic mode, at a flow rate of 1.0 mL/min. Detection was carried out by fluorescence absorbance at Ex/Em 230/290 nm. The multifunctional porous surface structure of the PPY-coated film provided high precision and accuracy for enantiomers. Compared with other commercial capillaries, PPY-coated capillary showed better extraction efficiency for all the analytes. The quantification limits of the proposed method were 10 ng/mL for R- and S-fluoxetine, and 15 ng/mL for R- and S-norfluoxetine, with a coefficient of variation lower than 13%. The response of the method for enantiomers is linear over a dynamic range, from the limit of quantification to 700ng/mL, with correlation coefficients higher than 0.9940. The in-tube SPME/LC method can therefore be successfully used to analyze plasma samples from ageing patients undergoing therapy with fluoxetine. (C) 2009 Elsevier B.V. All rights reserved.

Stabilization of native protein fold by intein-mediated covalent cyclization

A mutant version of the N-terminal domain of Escherichia coli DnaB helicase was used as a model system to assess the stabilization against unfolding gained by covalent cyclization. Cyclization was achieved in vivo by formation of an amide bond between the N and C termini with the help of a split mini-intein. Linear and circular proteins were constructed to be identical in amino acid sequence. Mutagenesis of Phe102 to Glu rendered the protein monomeric even at high concentration. A difference in free energy of unfolding, DeltaDeltaG, between circular and linear protein of 2.3(+/-0.5) kcal mol(-1) was measured at 10degreesC by circular dichroism. A theoretical estimate of the difference in conformational entropy of linear and circular random chains in a three-dimensional cubic lattice model predicted DeltaDeltaG = 2.3 kcal mol(-1), suggesting that stabilization by protein cyclization is driven by the reduced conformational entropy of the unfolded state. Amide-proton exchange rates measured by NMR spectroscopy and mass spectrometry showed a uniform, approximately tenfold decrease of the exchange rates of the most slowly exchanging amide protons, demonstrating that cyclization globally decreases the unfolding rate of the protein. The amide proton exchange was found to follow EX1 kinetics at near-neutral pH, in agreement with an unusually slow refolding I measured by stopped-flow circular dichroism. rate of less than 4 min(-1) The linear and circular proteins differed more in their unfolding than in their folding rates. Global unfolding of the N-terminal domain of E. coli DnaB is thus promoted strongly by spatial separation of the N and C termini, whereas their proximity is much less important for folding. (C) 2005 Elsevier Ltd. All rights reserved.

Social separation and diazepam withdrawal increase anxiety in the elevated plus-maze and serotonin turnover in the median raphe and hippocampus

The present work aimed to evaluate the effects of social separation for 14 days (chronic stress) and of withdrawal from a 14-day treatment with diazepam (acute stress) on the exploratory behaviour of male rats in the elevated plus-maze and on serotonin (5-hydroxytryptamine) turnover in different brain structures. Social separation had an anxiogenic effect, evidenced by fewer entries into, and less time spent on the open arms of the elevated plus-maze. Separation also selectively increased 5-hydroxytryptamine turnover in the hippocampus and median raphe nucleus. Diazepam withdrawal had a similar anxiogenic effect in grouped animals and increased 5-hydroxytryptamine turnover in the same brain structures. Chronic treatment with imipramine during the 14 days of separation prevented the behavioural and neurochemical changes caused by social separation. It is suggested that the increase in anxiety determined by both acute and chronic stress is mediated by the activation of the median raphe nucleus-hippocampal 5-hydroxytryptamine pathway.

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.