The unsteady flow of a weakly compressible fluid in a thin porous layer. I: Two-dimensional theory

We consider the problem of determining the pressure and velocity fields for a weakly compressible fluid flowing in a two-dimensional reservoir in an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting or extracting fluid. Numerical solution of this problem may be expensive, particularly in the case that the depth scale of the layer h is small compared to the horizontal length scale l. This is a situation which occurs frequently in the application to oil reservoir recovery. Under the assumption that epsilon=h/l<<1, we show that the pressure field varies only in the horizontal direction away from the wells (the outer region). We construct two-term asymptotic expansions in epsilon in both the inner (near the wells) and outer regions and use the asymptotic matching principle to derive analytical expressions for all significant process quantities. This approach, via the method of matched asymptotic expansions, takes advantage of the small aspect ratio of the reservoir, epsilon, at precisely the stage where full numerical computations become stiff, and also reveals the detailed structure of the dynamics of the flow, both in the neighborhood of wells and away from wells.

The unsteady flow of a weakly compressible fluid in a thin porous layer II: three-dimensional theory

We consider the problem of determining the pressure and velocity fields for a weakly compressible fluid flowing in a three-dimensional layer, composed of an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting and/or extracting fluid. Numerical solution of this three-dimensional evolution problem may be expensive, particularly in the case that the depth scale of the layer h is small compared to the horizontal length scale l, a situation which occurs frequently in the application to oil and gas reservoir recovery and which leads to significant stiffness in the numerical problem. Under the assumption that $\epsilon\propto h/l\ll 1$, we show that, to leading order in $\epsilon$, the pressure field varies only in the horizontal directions away from the wells (the outer region). We construct asymptotic expansions in $\epsilon$ in both the inner (near the wells) and outer regions and use the asymptotic matching principle to derive expressions for all significant process quantities. The only computations required are for the solution of non-stiff linear, elliptic, two-dimensional boundary-value, and eigenvalue problems. This approach, via the method of matched asymptotic expansions, takes advantage of the small aspect ratio of the layer, $\epsilon$, at precisely the stage where full numerical computations become stiff, and also reveals the detailed structure of the dynamics of the flow, both in the neighbourhood of wells and away from wells.

The unsteady flow of a weakly compressible fluid in a thin porous layer III: three-dimensional computations

We describe a novel method for determining the pressure and velocity fields for a weakly compressible fluid flowing in a thin three-dimensional layer composed of an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting and/or extracting fluid. Our approach uses the method of matched asymptotic expansions to derive expressions for all significant process quantities, the computation of which requires only the solution of linear, elliptic, two-dimensional boundary value and eigenvalue problems. In this article, we provide full implementation details and present numerical results demonstrating the efficiency and accuracy of our scheme.

Amniotic fluid iodine concentrations do not vary in pregnant women with varying iodine intake.

Iodine deficiency is an important clinical and public health problem. Its prevention begins with an adequate intake of iodine during pregnancy. International agencies recommend at least 200 microg iodine per d for pregnant women. We assessed whether iodine concentrations in the amniotic fluid of healthy pregnant women are independent of iodine intake. This cross-sectional, non-interventional study included 365 consecutive women who underwent amniocentesis to determine the fetal karyotype. The amniocentesis was performed with abdominal antisepsis using chlorhexidine. The iodine concentration was measured in urine and amniotic fluid. The study variables were the intake of iodized salt and multivitamin supplements or the prescription of a KI supplement. The mean level of urinary iodine was 139.0 (SD 94.5) microg/l and of amniotic fluid 15.81 (SD 7.09) microg/l. The women who consumed iodized salt and those who took a KI supplement had significantly higher levels of urinary iodine than those who did not (P = 0.01 and P = 0.004, respectively). The urinary iodine levels were not significantly different in the women who took a multivitamin supplement compared with those who did not take this supplement, independently of iodine concentration or multivitamin supplement. The concentrations of iodine in the amniotic fluid were similar, independent of the dietary iodine intake. Urine and amniotic fluid iodine concentrations were weakly correlated, although the amniotic fluid values were no higher in those women taking a KI supplement. KI prescription at recommended doses increases the iodine levels in the mother without influencing the iodine levels in the amniotic fluid.

Systematic weakly nonlinear analysis of interfacial instabilities in Hele-Shaw flows

We develop a systematic method to derive all orders of mode couplings in a weakly nonlinear approach to the dynamics of the interface between two immiscible viscous fluids in a Hele-Shaw cell. The method is completely general: it applies to arbitrary geometry and driving. Here we apply it to the channel geometry driven by gravity and pressure. The finite radius of convergence of the mode-coupling expansion is found. Calculation up to third-order couplings is done, which is necessary to account for the time-dependent Saffman-Taylor finger solution and the case of zero viscosity contrast. The explicit results provide relevant analytical information about the role that the viscosity contrast and the surface tension play in the dynamics of the system. We finally check the quantitative validity of different orders of approximation and a resummation scheme against a physically relevant, exact time-dependent solution. The agreement between the low-order approximations and the exact solution is excellent within the radius of convergence, and is even reasonably good beyond this radius.

Systematic weakly nonlinear analysis of radial viscous fingering

We present a weakly nonlinear analysis of the interface dynamics in a radial Hele-Shaw cell driven by both injection and rotation. We extend the systematic expansion introduced in [E. Alvarez-Lacalle et al., Phys. Rev. E 64, 016302 (2001)] to the radial geometry, and compute explicitly the first nonlinear contributions. We also find the necessary and sufficient condition for the uniform convergence of the nonlinear expansion. Within this region of convergence, the analytical predictions at low orders are compared satisfactorily to exact solutions and numerical integration of the problem. This is particularly remarkable in configurations (with no counterpart in the channel geometry) for which the interplay between injection and rotation allows that condition to be satisfied at all times. In the case of the purely centrifugal forcing we demonstrate that nonlinear couplings make the interface more unstable for lower viscosity contrast between the fluids.

Magmatic-dominated fluid evolution in the Jurassic Nambija gold skarn deposits (southeastern Ecuador)

The Jurassic (approximately 145 Ma) Nambija oxidized gold skarns are hosted by the Triassic volcanosedimentary Piuntza unit in the sub-Andean zone of southeastern Ecuador. The skarns consist dominantly of granditic garnet (Ad(20-98)) with subordinate pyroxene (Di(46-92)Hd(17-42)Jo(0-19)) and epidote and are spatially associated with porphyritic quartz-diorite to granodiorite intrusions. Endoskarn is developed at the intrusion margins and grades inwards into a potassic alteration zone. Exoskarn has an outer K- and Na-enriched zone in the volcanosedimentary unit. Gold mineralization is associated with the weakly developed retrograde alteration of the exoskarn and occurs mainly in sulfide-poor vugs and milky quartz veins and veinlets in association with hematite. Fluid inclusion data for the main part of the prograde stage indicate the coexistence of high-temperature (500A degrees C to > 600A degrees C), high-salinity (up to 65 wt.% eq. NaCl), and moderate- to low-salinity aqueous-carbonic fluids interpreted to have been trapped at pressures around 100-120 MPa, corresponding to about 4-km depth. Lower-temperature (510-300A degrees C) and moderate- to low-salinity (23-2 wt.% eq. NaCl) aqueous fluids are recorded in garnet and epidote of the end of the prograde stage. The microthermometric data (Th from 513A degrees C to 318A degrees C and salinity from 1.0 to 23 wt.% eq. NaCl) and delta(18)O values between 6.2aEuro degrees and 11.5aEuro degrees for gold-bearing milky quartz from the retrograde stage suggest that the ore-forming fluid was dominantly magmatic. Pressures during the early retrograde stage were in the range of 50-100 MPa, in line with the evidence for CO(2) effervescence and probable local boiling. The dominance of magmatic low-saline to moderately saline oxidizing fluids during the retrograde stage is consistent with the depth of the skarn system, which could have delayed the ingression of external fluids until relatively low temperatures were reached. The resulting low water-to-rock ratios explain the weak retrograde alteration and the compositional variability of chlorite, essentially controlled by host rock compositions. Gold was precipitated at this stage as a result of cooling and pH increase related to CO(2) effervescence, which both result in destabilization of gold-bearing chloride complexes. Significant ingression of external fluids took place after gold deposition only, as recorded by delta(18)O values of 0.4aEuro degrees to 6.2aEuro degrees for fluids depositing quartz (below 350A degrees C) in sulfide-rich barren veins. Low-temperature (< 300A degrees C) meteoric fluids (delta(18)O(water) between -10.0aEuro degrees and -2.0aEuro degrees) are responsible for the precipitation of late comb quartz and calcite in cavities and veins and indicate mixing with cooler fluids of higher salinities (about 100A degrees C and 25 wt.% eq. NaCl). The latter are similar to low-temperature fluids (202-74.5A degrees C) with delta(18)O values of -0.5aEuro degrees to 3.1aEuro degrees and salinities in the range of 21.1 to 17.3 wt.% eq. CaCl(2), trapped in calcite of late veins and interpreted as basinal brines. Nambija represents a deep equivalent of the oxidized gold skarn class, the presence of CO(2) in the fluids being partly a consequence of the relatively deep setting at about 4-km depth. As in other Au-bearing skarn deposits, not only the prograde stage but also the gold-precipitating retrograde stage is dominated by fluids of magmatic origin.

Chemoconvection patterns in the Methylene-blue-glucose system: Weakly nonlinear analysis

The oxidation of solutions of glucose with methylene-blue as a catalyst in basic media can induce hydrodynamic overturning instabilities, termed chemoconvection in recognition of their similarity to convective instabilities. The phenomenon is due to gluconic acid, the marginally dense product of the reaction, which gradually builds an unstable density profile. Experiments indicate that dominant pattern wavenumbers initially increase before gradually decreasing or can even oscillate for long times. Here, we perform a weakly nonlinear analysis for an established model of the system with simple kinetics, and show that the resulting amplitude equation is analogous to that obtained in convection with insulating walls. We show that the amplitude description predicts that dominant pattern wavenumbers should decrease in the long term, but does not reproduce the aforementioned increasing wavenumber behavior in the initial stages of pattern development. We hypothesize that this is due to horizontally homogeneous steady states not being attained before pattern onset. We show that the behavior can be explained using a combination of pseudo-steady-state linear and steady-state weakly nonlinear theories. The results obtained are in qualitative agreement with the analysis of experiments.

Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.

Limits of validity for a semiclassical mean-field two-fluid model for Bose-Einstein condensation thermodynamics

We reinvestigate the Bose-Einstein condensation (BEC) thermodynamics of a weakly interacting dilute Bose gas under the action of a trap using a semi-classical two-fluid mean-field model in order to find the domain of applicability of the model. Such a model is expected to break down once the condition of diluteness and weak interaction is violated. We find that this breakdown happens for values of coupling and density near the present experimental scenario of BEG. With the increase of the interaction coupling and density the model may lead to unphysical results for thermodynamic observables. (C) 2000 Published by Elsevier B.V. B.V, All rights reserved.

Shear viscosity from thermal fluctuations in relativistic conformal fluid dynamics

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The PSE-3D instability analysis methodology for flows depending strongly on two and weakly on the third spatial dimension

The present contribution discusses the development of a PSE-3D instability analysis algorithm, in which a matrix forming and storing approach is followed. Alternatively to the typically used in stability calculations spectral methods, new stable high-order finitedifference-based numerical schemes for spatial discretization 1 are employed. Attention is paid to the issue of efficiency, which is critical for the success of the overall algorithm. To this end, use is made of a parallelizable sparse matrix linear algebra package which takes advantage of the sparsity offered by the finite-difference scheme and, as expected, is shown to perform substantially more efficiently than when spectral collocation methods are used. The building blocks of the algorithm have been implemented and extensively validated, focusing on classic PSE analysis of instability on the flow-plate boundary layer, temporal and spatial BiGlobal EVP solutions (the latter necessary for the initialization of the PSE-3D), as well as standard PSE in a cylindrical coordinates using the nonparallel Batchelor vortex basic flow model, such that comparisons between PSE and PSE-3D be possible; excellent agreement is shown in all aforementioned comparisons. Finally, the linear PSE-3D instability analysis is applied to a fully three-dimensional flow composed of a counter-rotating pair of nonparallel Batchelor vortices.

Range of validity of weakly non-linear theory in Rayleigh-Bénard problem

In this paper we examine the equilibrium states of finite amplitude flow in a horizontal fluid layer with differential heating between the two rigid boundaries. The solutions to the Navier-Stokes equations are obtained by means of a perturbation method for evaluating the Landau constants and through a Newton-Raphson iterative method that results from the Fourier expansion of the solutions that bifurcate above the linear stability threshold of infinitesimal disturbances. The results obtained from these two different methods of evaluating the convective flow are compared in the neighborhood of the critical Rayleigh number. We find that for small Prandtl numbers the discrepancy of the two methods is noticeable. © 2009 The Physical Society of Japan.

Human Herpesvirus Infections Of The Central Nervous System: Laboratory Diagnosis Based On Dna Detection By Nested Pcr In Plasma And Cerebrospinal Fluid Samples.

Infections of the central nervous systems (CNS) present a diagnostic problem for which an accurate laboratory diagnosis is essential. Invasive practices, such as cerebral biopsy, have been replaced by obtaining a polymerase chain reaction (PCR) diagnosis using cerebral spinal fluid (CSF) as a reference method. Tests on DNA extracted from plasma are noninvasive, thus avoiding all of the collateral effects and patient risks associated with CSF collection. This study aimed to determine whether plasma can replace CSF in nested PCR analysis for the detection of CNS human herpesvirus (HHV) diseases by analysing the proportion of patients whose CSF nested PCR results were positive for CNS HHV who also had the same organism identified by plasma nested PCR. In this study, CSF DNA was used as the gold standard, and nested PCR was performed on both types of samples. Fifty-two patients with symptoms of nervous system infection were submitted to CSF and blood collection. For the eight HHV, one positive DNA result-in plasma and/or CSF nested PCR-was considered an active HHV infection, whereas the occurrence of two or more HHVs in the same sample was considered a coinfection. HHV infections were positively detected in 27/52 (51.9%) of the CSF and in 32/52 (61.5%) of the plasma, difference not significant, thus nested PCR can be performed on plasma instead of CSF. In conclusion, this findings suggest that plasma as a useful material for the diagnosis of cases where there is any difficulty to perform a CSF puncture.

Central nitrergic system regulation of neuroendocrine secretion, fluid intake and blood pressure induced by angiotensin-II

Background: Nitric oxide (NO) synthesis has been described in several circumventricular and hypothalamic structures in the central nervous system that are implicated in mediating central angiotensin-II (ANG-II) actions during water deprivation and hypovolemia. Neuroendocrine and cardiovascular responses, drinking behavior, and urinary excretions were examined following central angiotensinergic stimulation in awake freely-moving rats pretreated with intracerebroventricular injections of N omega-nitro-L-arginine methyl ester (L-NAME, 40 mu g), an inhibitor of NO synthase, and L-arginine (20 ug), a precursor of NO. Results: Injections of L-NAME or ANG-II produced an increase in plasma vasopressin (VP), oxytocin (OT) and atrial natriuretic peptide (ANP) levels, an increase in water and sodium intake, mean arterial blood pressure and sodium excretion, and a reduction of urinary volume. L-NAME pretreatment enhanced the ANG-II response, while L-arginine attenuated VP and OT release, thirst, appetite for sodium, antidiuresis, and natriuresis, as well as pressor responses induced by ANG-II. Discussion and conclusion: Thus, the central nitrergic system participates in the angiotensinergic responses evoked by water deprivation and hypovolemia to refrain neurohypophysial secretion, hydromineral balance, and blood pressure homeostasis.