Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

Impaired chemical coupling of cerebral blood flow is compatible with intact neurological outcome in neonates with perinatal risk factors.

Early detection of pathophysiological factors associated with permanent brain damage is a major issue in neonatal medicine. The aim of our study was to evaluate the significance of the CO2 reactivity of cerebral blood flow (CBF) in neonates with perinatal risk factors. Fourteen ventilated neonates with perinatal risk factors (pathological cardiotocogramm, low cord pH, postpartal encephalopathy) were enrolled into this prospective study. The study was performed 18-123 h after birth. CBF was measured using the noninvasive intravenous 133Xe method. Two measurements were taken with a minimal PaCO2-difference of 5 mm Hg. From the two CBF values the CO2 reactivity was calculated. Outcome was evaluated 1 year after birth. The CBF values at a lower PaCO2 ranged from 6.6 to 115. 2 ml/100 g brain issue/min (median = 18.2) and at a higher PaCO2 level from 7.1 to 125.7 ml/100 g brain tissue/min (median = 18.75). The calculated CO2 reactivity ranged from -9.6 to 6.6% (median 1.1%) change in CBF/mm Hg change in PaCO2. CO2 reactivity correlated with lowest pH (r2 = 0.35, p = 0.02). Two infants died, one of neonatal sepsis, the other of heart failure. Neurological outcome at the age of 1 year was normal in 11 patients, 1 had severe cerebral palsy. From the 12 surviving patients the patient with severe neurological deficit showed the highest CBF values (125.7 ml/100 g/min). Impaired chemical coupling of cerebral blood flow is compatible with intact neurological outcome in neonates with perinatal risk factors. CO2 reactivity in these newborns correlates with the lowest pH and may reflect the severity of perinatal asphyxia.

Estimating gene flow across hybrid zones: how reliable are microsatellites?

Due to their high polymorphism, microsatellites have become one of the most valued genetic markers in population biology. We review the first two published studies on hybrid zones of the common shrew based on microsatellites. Both reveal surprisingly high interracial gene flow. It can be shown that these are overestimates. Indeed, in classical population genetics models as F-statistics, mutation is neglected. This constitutes an acceptable assumption as long as migration is higher than mutation. However, in hybrid zones where genetic exchanges can be rare, neglecting mutation will lead to strong overestimates of migration when working with microsatellites which display mutation rates up to 10(-3). As there seems to be no straightforward way to correct for this bias, interracial gene flow estimates based on microsatellites should be taken with caution. This problem should however not conceal the enormous potential of microsatellites to unravel the genetics of hybrid zones.

A study of the effect of sartorius transposition on lymph flow after ilioinguinal node dissection.

Ilioinguinal dissection is associated with a high rate of lymphatic complications. Prolonged lymph flow causes greatest concern and preventive strategies are needed. A retrospective study of 28 consecutive patients undergoing groin dissection for melanoma metastases was performed to evaluate the influence of sartorius muscle transposition on lymph flow. Modification of the surgical technique with transposition of the sartorius muscle was not associated with reduced drainage time (P = 0.66). A 2-staged approach, with initial sentinel lymph node resection and lymph node dissection in a second operation, however, lead to shortened duration of the lymph flow (P = 0.01). Prolonged lymphorrhea was more frequent in older (P = 0.03), obese (P = 0.02) patients affected by diabetes mellitus (P = 0.03) and hypertension (P = 0.04).

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.

No major impact of skin aging on the response of skin blood flow to a submaximal local thermal stimulus.

OBJECTIVE: This study was undertaken to investigate how aging affects dermal microvascular reactivity in skin areas differentially exposed to sunlight, and therefore to different degrees of photoaging. METHODS: We assessed, in young (18-30 years, n = 13) and aged males (≥60 years, n = 13), the thigh, forearm, and forehead's skin vasodilatory response to local heating (LTH) with a LDI. In each subject and at each location, local Tskin was brought from 34°C (baseline) to 39 or 41°C for 30 minutes, to effect submaximal vasodilation, with maximal vasodilation then elicited by further heating to 44°C. RESULTS: The CVCs evaluated at baseline and after maximal vasodilation (CVCmax ) were higher in the forehead than in the two other anatomical locations. On all locations, CVCmax decreased with age but less markedly in the forehead compared to the two other locations. When expressed in % of CVCmax , the plateau increase of CVCs in response to submaximal temperatures (39 and 41°C) did not vary with age, and minimally so with location. CONCLUSION: Skin aging, whether intrinsic or combined with photoaging, reduces the maximal vasodilatory capacity of the dermal microcirculation, but not its reactivity to local heating.

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.

Superior flow for bridge to life with self-expanding venous cannulas.

BACKGROUND: Recently, a compact cardiopulmonary support (CPS) system designed for quick set-up for example, during emergency cannulation, has been introduced. Traditional rectilinear percutaneous cannulas are standard for remote vascular access with the original design. The present study was designed to assess the potential of performance increase by the introduction of next-generation, self-expanding venous cannulas, which can take advantage of the luminal width of the venous vasculature despite a relatively small access orifice. METHODS: Veno-arterial bypass was established in three bovine experiments (69+/-10 kg). The Lifebridge (Lifebridge GmbH, Munich, Germany) system was connected to the right atrium in a trans-jugular fashion with various venous cannulas; and the oxygenated blood was returned through the carotid artery with a 17 F percutaneous cannula. Two different venous cannulas were studied, and the correlation between the centrifugal pump speed (1500-3900 RPM), flow and the required negative pressure on the venous side was established: (A) Biomedicus 19 F (Medtronic, Tolochenaz, Switzerland); (B) Smart canula 18 F/36 F (Smartcanula LLC, Lausanne, Switzerland). RESULTS: At 1500 RPM, the blood flow was 0.44+/-0.26 l min(-1) for the 19 F rectilinear cannula versus 0.73+/-0.34 l min(-1) for the 18/36 F self-expanding cannula. At 2500 RPM the blood flow was 1.63+/-0.62 l min(-1) for the 19F rectilinear cannula versus 2.13+/-0.34 l min(-1) for the 18/36 F self-expanding cannula. At 3500 RPM, the blood flow was 2.78+/-0.47 l min(-1) for the 19 F rectilinear cannula versus 3.64+/-0.39 l min(-1) for the 18/36 F self-expanding cannula (p<0.01 for 18/36 F vs 19 F). At 1500 RPM, the venous line pressure was 18+/-8 mmHg for the 19F rectilinear cannula versus 19+/-5 mmHg for the 18/36 F self-expanding cannula. At 2500 RPM the venous line pressure accounted for -22+/-32 mmHg for the 19 F rectilinear cannula versus 2+/-5 mmHg for the 18/36 F self-expanding cannula. At 3500 RPM, the venous line pressure was -112+/-42 mmHg for the rectilinear cannula versus 28+/-7 mmHg for the 18/36 F self-expanding cannula (p<0.01 for 18 F/36 F vs 19 F). Conclusions: The negative pressure required to achieve adequate venous drainage with the self-expanding venous cannula accounts for approximately 31% of the pressure necessary with the 19 F rectilinear cannula. In addition, a pump flow of more than 4 l min(-1) can be achieved with the self-expanding design and a well-accepted negative inlet pressure for minimal blood trauma of less than 50 mmHg.

Quantification of Myocardial Blood Flow in Absolute Terms Using (82)Rb PET Imaging: The RUBY-10 Study.

OBJECTIVES: The purpose of this study was to compare myocardial blood flow (MBF) and myocardial flow reserve (MFR) estimates from rubidium-82 positron emission tomography ((82)Rb PET) data using 10 software packages (SPs) based on 8 tracer kinetic models. BACKGROUND: It is unknown how MBF and MFR values from existing SPs agree for (82)Rb PET. METHODS: Rest and stress (82)Rb PET scans of 48 patients with suspected or known coronary artery disease were analyzed in 10 centers. Each center used 1 of 10 SPs to analyze global and regional MBF using the different kinetic models implemented. Values were considered to agree if they simultaneously had an intraclass correlation coefficient >0.75 and a difference <20% of the median across all programs. RESULTS: The most common model evaluated was the Ottawa Heart Institute 1-tissue compartment model (OHI-1-TCM). MBF values from 7 of 8 SPs implementing this model agreed best. Values from 2 other models (alternative 1-TCM and Axially distributed) also agreed well, with occasional differences. The MBF results from other models (e.g., 2-TCM and retention) were less in agreement with values from OHI-1-TCM. CONCLUSIONS: SPs using the most common kinetic model-OHI-1-TCM-provided consistent results in measuring global and regional MBF values, suggesting that they may be used interchangeably to process data acquired with a common imaging protocol.

Toxicity of triclosan, penconazole and metalaxyl on Caulobacter crescentus and a freshwater microbial community as assessed by flow cytometry.

Biocides are widely used for domestic hygiene, agricultural and industrial applications. Their widespread use has resulted in their introduction into the environment and raised concerns about potential deleterious effects on aquatic ecosystems. In this study, the toxicity of the biocides triclosan, penconazole and metalaxyl were evaluated with the freshwater bacterium Caulobacter crescentus and with a freshwater microbial community using a combination of single- and double-stain flow cytometric assays. Growth of C.  crescentus and the freshwater community were repressed by triclosan but not by penconazole or metalaxyl at concentrations up to 250 μM. The repressive effect of triclosan was dependent on culture conditions. Caulobacter crescentus was more sensitive to triclosan when grown with high glucose at high cell density than when grown directly in sterilized lake water at low cell density. This suggests that the use of conventional growth conditions may overestimate biocide toxicity. Additional experiments showed that the freshwater community was more sensitive to triclosan than C.  crescentus, with 10 nM of triclosan being sufficient to repress growth and change the phylogenetic composition of the community. These results demonstrate that isolate-based assays may underestimate biocide toxicity and highlight the importance of assessing toxicity directly on natural microbial communities. Because 10 nM of triclosan is within the range of concentrations observed in freshwater systems, these results also raise concerns about the risk of introducing triclosan into the environment.

Sediment budget monitoring of a debris-flow torrent (French Prealps)

The Manival near Grenoble (French Prealps) is a very active debris-flow torrent equipped with a large sediment trap (25 000 m3) protecting an urbanized alluvial fan from debris-flows. We began monitoring the sediment budget of the catchment controlled by the trap in Spring 2009. Terrestrial laser scanner is used for monitoring topographic changes in a small gully, the main channel, and the sediment trap. In the main channel, 39 cross-sections are surveyed after every event. Three periods of intense geomorphic activity are documented here. The first was induced by a convective storm in August 2009 which triggered a debris-flow that deposited ~1,800 m3 of sediment in the trap. The debris-flow originated in the upper reach of the main channel and our observations showed that sediment outputs were entirely supplied by channel scouring. Hillslope debris-flows were initiated on talus slopes, as revealed by terrestrial LiDAR resurveys; however they were disconnected to the main channel. The second and third periods of geomorphic activity were induced by long duration and low intensity rainfall events in September and October 2009 which generate small flow events with intense bedload transport. These events contribute to recharge the debris-flow channel with sediments by depositing important gravel dunes propagating from headwaters. The total recharge in the torrent subsequent to bedload transport events was estimated at 34% of the sediment erosion induced by the August debris-flow.

Simultaneous measurement of regional cerebral blood flow by perfusion CT and stable xenon CT: a validation study.

BACKGROUND AND PURPOSE: Knowledge of cerebral blood flow (CBF) alterations in cases of acute stroke could be valuable in the early management of these cases. Among imaging techniques affording evaluation of cerebral perfusion, perfusion CT studies involve sequential acquisition of cerebral CT sections obtained in an axial mode during the IV administration of iodinated contrast material. They are thus very easy to perform in emergency settings. Perfusion CT values of CBF have proved to be accurate in animals, and perfusion CT affords plausible values in humans. The purpose of this study was to validate perfusion CT studies of CBF by comparison with the results provided by stable xenon CT, which have been reported to be accurate, and to evaluate acquisition and processing modalities of CT data, notably the possible deconvolution methods and the selection of the reference artery. METHODS: Twelve stable xenon CT and perfusion CT cerebral examinations were performed within an interval of a few minutes in patients with various cerebrovascular diseases. CBF maps were obtained from perfusion CT data by deconvolution using singular value decomposition and least mean square methods. The CBF were compared with the stable xenon CT results in multiple regions of interest through linear regression analysis and bilateral t tests for matched variables. RESULTS: Linear regression analysis showed good correlation between perfusion CT and stable xenon CT CBF values (singular value decomposition method: R(2) = 0.79, slope = 0.87; least mean square method: R(2) = 0.67, slope = 0.83). Bilateral t tests for matched variables did not identify a significant difference between the two imaging methods (P >.1). Both deconvolution methods were equivalent (P >.1). The choice of the reference artery is a major concern and has a strong influence on the final perfusion CT CBF map. CONCLUSION: Perfusion CT studies of CBF achieved with adequate acquisition parameters and processing lead to accurate and reliable results.

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.