Spatiotemporal adaptive multiscale multiphysics simulations of two-phase flow

We present a spatiotemporal adaptive multiscale algorithm, which is based on the Multiscale Finite Volume method. The algorithm offers a very efficient framework to deal with multiphysics problems and to couple regions with different spatial resolution. We employ the method to simulate two-phase flow through porous media. At the fine scale, we consider a pore-scale description of the flow based on the Volume Of Fluid method. In order to construct a global problem that describes the coarse-scale behavior, the equations are averaged numerically with respect to auxiliary control volumes, and a Darcy-like coarse-scale model is obtained. The space adaptivity is based on the idea that a fine-scale description is only required in the front region, whereas the resolution can be coarsened elsewhere. Temporal adaptivity relies on the fact that the fine-scale and the coarse-scale problems can be solved with different temporal resolution (longer time steps can be used at the coarse scale). By simulating drainage under unstable flow conditions, we show that the method is able to capture the coarse-scale behavior outside the front region and to reproduce complex fluid patterns in the front region.

Detection and characterization of hydraulically active fractures in a carbonate aquifer: Results from self-potential, temperature and fluid conductivity logging along a 260-m-deep borehole in the Combioula hydrothermal system in the southwestern Swiss Alps

A geophysical and geochemical study has been conducted in a fractured carbonate aquifer located at Combioula in the southwestern Swiss Alps with the objective to detect and characterize hydraulically active fractures along a 260-m-deep borehole. Hydrochemical analyses, borehole diameter, temperature and fluid electrical conductivity logging data were integrated in order to relate electrokinetic self-potential signals to groundwater flow inside the fracture network. The results show a generally good, albeit locally variable correlation of variations of the self-potential signals with variations in temperature, fluid electrical conductivity and borehole diameter. Together with the hydrochemical evidence, which was found to be critical for the interpretation of the self-potential data, these measurements not only made it possible to detect the hydraulically active fractures but also to characterize them as zones of fluid gain or fluid loss. The results complement the available information from the corresponding litholog and illustrate the potential of electrokinetic self-potential signals in conjunction with temperature, fluid electrical conductivity and hydrochemical analyses for the characterization of fractured aquifers, and thus may offer a perspective for an effective quantitative characterization of this increasingly important class of aquifers and geothermal reservoirs.

Fluorescence flow cytometer to determine urine particle reference intervals in doping control samples.

Background: Urine is still the matrix of choice to fight against doping, because it can be collected non-invasively during anti-doping tests. Most of the World Anti-Doping Agency's accredited laboratories have more than 20 years experience in analyzing this biological fluid and the majority of the compounds listed in the 2010 Prohibited List - International Standard are eliminated through the urinary apparatus. Storing and transporting urine samples for doping analyses does not include a specific protocol to prevent microbial and thermal degradation. The use of a rapid and reliable screening method could enable determine reference intervals for urine specimens in doping control samples and evaluate notably the prevalence of microbial contamination known to be responsible for the degradation of chemical substances in urine.Methods: The Sysmex(R) UF-500i is a recent urine flow cytometer analyzer capable of quantifying BACT and other urinary particles such as RBC, WBC, EC, DEBRIS, CAST, PATH. CAST, YLC, SRC as well as measuring urine conductivity. To determine urine anti-doping reference intervals, 501 samples received in our laboratory over a period of two months were submitted to an immediate examination. All samples were collected and then transported at room temperature. Analysis of variance was performed to test the effects of factors such as gender, test type [in-competition, out-of-competition] and delivery time.Results: The data obtained showed that most of the urine samples were highly contaminated with bacteria. The other urine particles were also very different according to the factors.Conclusions: The Sysmex(R) UF-500i was capable of providing a snapshot of urine particles present in the samples at the time of the delivery to the laboratory. These particles, BACT in particular, gave a good idea of the possible microbial degradation which had and/or could have occurred in the sample. This information could be used as the first quality control set up in WADA (World Anti-Doping Agency) accredited laboratories to determine if steroid profiles, endogenous and prohibited substances have possibly been altered. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Extension of Murray's law using a non-Newtonian model of blood flow.

BACKGROUND: So far, none of the existing methods on Murray's law deal with the non-Newtonian behavior of blood flow although the non-Newtonian approach for blood flow modelling looks more accurate. MODELING: In the present paper, Murray's law which is applicable to an arterial bifurcation, is generalized to a non-Newtonian blood flow model (power-law model). When the vessel size reaches the capillary limitation, blood can be modeled using a non-Newtonian constitutive equation. It is assumed two different constraints in addition to the pumping power: the volume constraint or the surface constraint (related to the internal surface of the vessel). For a seek of generality, the relationships are given for an arbitrary number of daughter vessels. It is shown that for a cost function including the volume constraint, classical Murray's law remains valid (i.e. SigmaR(c) = cste with c = 3 is verified and is independent of n, the dimensionless index in the viscosity equation; R being the radius of the vessel). On the contrary, for a cost function including the surface constraint, different values of c may be calculated depending on the value of n. RESULTS: We find that c varies for blood from 2.42 to 3 depending on the constraint and the fluid properties. For the Newtonian model, the surface constraint leads to c = 2.5. The cost function (based on the surface constraint) can be related to entropy generation, by dividing it by the temperature. CONCLUSION: It is demonstrated that the entropy generated in all the daughter vessels is greater than the entropy generated in the parent vessel. Furthermore, it is shown that the difference of entropy generation between the parent and daughter vessels is smaller for a non-Newtonian fluid than for a Newtonian fluid.

Analysis of flow dynamics in right ventricular outflow tract.

BACKGROUND: The mechanism behind early graft failure after right ventricular outflow tract (RVOT) reconstruction is not fully understood. Our aim was to establish a three-dimensional computational fluid dynamics (CFD) model of RVOT to investigate the hemodynamic conditions that may trigger the development of intimal hyperplasia and arteriosclerosis. METHODS: Pressure, flow, and diameter at the RVOT, pulmonary artery (PA), bifurcation of the PA, and left and right PAs were measured in 10 normal pigs with a mean weight of 24.8 ± 0.78 kg. Data obtained from the experimental scenario were used for CFD simulation of pressure, flow, and shear stress profile from the RVOT to the left and right PAs. RESULTS: Using experimental data, a CFD model was obtained for 2.0 and 2.5-L/min pulsatile inflow profiles. In both velocity profiles, time and space averaged in the low-shear stress profile range from 0-6.0 Pa at the pulmonary trunk, its bifurcation, and at the openings of both PAs. These low-shear stress areas were accompanied to high-pressure regions 14.0-20.0 mm Hg (1866.2-2666 Pa). Flow analysis revealed a turbulent flow at the PA bifurcation and ostia of both PAs. CONCLUSIONS: Identified local low-shear stress, high pressure, and turbulent flow correspond to a well-defined trigger pattern for the development of intimal hyperplasia and arteriosclerosis. As such, this real-time three-dimensional CFD model may in the future serve as a tool for the planning of RVOT reconstruction, its analysis, and prediction of outcome.

Fluid dynamics simulation of right ventricular outflow tract oversizing.

OBJECTIVES: Repair of the right ventricular outflow tract (RVOT) in paediatric cardiac surgery remains challenging due to the high reoperation rate. Intimal hyperplasia and consequent arteriosclerosis is one of the most important limitation factors for graft durability. Since local shear stress and pressure are predictive elements for intimal hyperplasia and wall degeneration, we sought to determine in an oversized 12-mm RVOT model, with computed fluid dynamics simulation, the local haemodynamical factors that may explain intimal hyperplasia. This was done with the aim of identifying the optimal degree of oversizing for a 12-mm native RVOT. METHODS: Twenty domestic pigs, with a weight of 24.6 ± 0.89 kg and a native RVOT diameter of 12 ± 1.7 mm, had valve conduits of 12, 16, 18 and 20 mm implanted. Pressure and flow were measured at 75, 100 and 125% of normal flow at RVOT at the pulmonary artery, pulmonary artery bifurcation and at the left and right pulmonary arteries. Three-dimensional computed fluid dynamics (CFD) simulation in all four geometries in all flow modalities was performed. Local shear stress and pressure conditions were investigated. RESULTS: Corresponding to 75, 100 and 125% of steady-state flow, three inlet velocity profiles were obtained, 0.2, 0.29 and 0.36 m/s, respectively. At inflow velocity profiles, low shear stress areas, ranged from 0 to 2 Pa, combined with high-pressure areas ranging from 11.5 to 12.1 mmHg that were found at distal anastomosis, at bifurcation and at the ostia of the left and right pulmonary arteries in all geometries. CONCLUSIONS: In all three oversized geometries, the local reparation of shear stress and pressure in the 16-mm model showed a similar local profile as in the native 12 mm RVOT. According to these findings, we suggest oversizing the natural 12-mm RVOT by not more than 4 mm. The elements responsible for wall degeneration and intimal hyperplasia remain very similar to the conditions present in native RVOT.

FOXC2 and fluid shear stress stabilize postnatal lymphatic vasculature.

Biomechanical forces, such as fluid shear stress, govern multiple aspects of endothelial cell biology. In blood vessels, disturbed flow is associated with vascular diseases, such as atherosclerosis, and promotes endothelial cell proliferation and apoptosis. Here, we identified an important role for disturbed flow in lymphatic vessels, in which it cooperates with the transcription factor FOXC2 to ensure lifelong stability of the lymphatic vasculature. In cultured lymphatic endothelial cells, FOXC2 inactivation conferred abnormal shear stress sensing, promoting junction disassembly and entry into the cell cycle. Loss of FOXC2-dependent quiescence was mediated by the Hippo pathway transcriptional coactivator TAZ and, ultimately, led to cell death. In murine models, inducible deletion of Foxc2 within the lymphatic vasculature led to cell-cell junction defects, regression of valves, and focal vascular lumen collapse, which triggered generalized lymphatic vascular dysfunction and lethality. Together, our work describes a fundamental mechanism by which FOXC2 and oscillatory shear stress maintain lymphatic endothelial cell quiescence through intercellular junction and cytoskeleton stabilization and provides an essential link between biomechanical forces and endothelial cell identity that is necessary for postnatal vessel homeostasis. As FOXC2 is mutated in lymphedema-distichiasis syndrome, our data also underscore the role of impaired mechanotransduction in the pathology of this hereditary human disease.