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.

Influence de deux milieux séquentiels pour la culture d'embryons sur la production d'hormone gonadotrophine chorionique et de progestérone par des cellules JEG-3 et BeWo [Effects of different embryo development media on hCG and progesterone release by JEG-3 and BeWo cells]

Current in vitro fertilisation (IVF) practice requires synchronisation between the¦environment of cultured oocytes and embryos and the surroundings to what they would have¦been exposed to in vivo. Commercial, sequential media follow this requirement but their exact¦composition is not available. We have compared two widely used IVF culture media systems using¦the two choriocarcinoma cell lines JEG-3 and BeWo. The two hormones hCG and progesterone¦were determined in the culture supernatants as endpoints. In both cell lines, but in a more¦pronounced way in JEG-3, progesterone rather than hCG production was stimulated, and a¦higher hormone release was observed in the fertilisation than in the cleavage media. Differences¦between manufacturers were small and did not favour one system over the other. We conclude¦that both sequential media systems can be equally well used in current IVF laboratory practice.¦© 2012 Elsevier Masson SAS. All rights reserved.

Politicians in the Media: Determinants of Legislators' Presence and Prominence in Swiss Newspapers

In times of increasing "mediatization" of politics, when voters and their elected representatives primarily communicate through the media, the question of who gets into the news and why becomes of the utmost importance. This article examines the determinants of Swiss legislators' presence and prominence in the print media by focusing on three competing approaches drawn from communication studies. The first approach regards the media as a "mirror" of political reality and argues that the media focus on the most active deputies in parliament. Second, news values theory predicts that "authoritative" politicians in leadership positions get the most media coverage. Third, theories of "news bias" hold that the media privilege legislators who are in line with their own editorial interests. Overall, the statistical analyses show an important leadership effect and provide strong support for the second explanation. While deputies in official functions get the most extensive news coverage, media access can also be won by parliamentary activity. The least support is shown for the news bias theory, although some newspapers try to localize parliamentary news coverage by focusing on deputies from their own media market.

The Media and Socio-Cultural Regulation of Religion

(Résumé de l'ouvrage) Originale, insolite, renaissante, l'action religieuse émergente bouscule les habitudes, ébranle les certitudes, construit ici, maintenant, l'autre monde. Peut-on courir le risque? Voilà que la question se pose et se résout en rumeurs publiques, poursuites judiciaires et tensions scolaires, lesquelles mettent à nu des mécanismes inédits d'institutionnalisation de l'expérience religieuse en modernité. As new religious movements seek to carve out their own niche in society, public controversy and opposing beliefs can spark bitter debates, and can even lead to calls for state intervention. How then do new or borderline religious groups negotiate or mediate the building of public space?

Delayed otogenic hydrocephalus after acute otitis media in pediatric patients: the changing presentation of a serious otologic complication.

CONCLUSIONS: The clinical presentation of otogenic dural sinus thrombosis (DST) as a complication of acute otitis media (AOM) can be masked by antibiotic treatment. Morning episodes of vomiting and/or headache, visual impairment and a history of AOM seem to be indicative of otogenic hydrocephalus. We therefore advocate that the MRI scans of patients with similar symptoms should be carefully studied to facilitate the early diagnosis of a potentially life-threatening complication. OBJECTIVE: To describe the frequency, pathognomonic signs, clinical course and outcome of otogenic hydrocephalus and DST as complications of AOM in pediatric patients. MATERIAL AND METHODS: We undertook a retrospective chart review of all pediatric patients (age 1-14 years) treated for otitis media and its complications at an academic medical center between 1999 and 2003. The main outcome measures were otologic and ophthalmologic findings and CT and MRI scans at the beginning of treatment and 3 months later. RESULTS: We report on five cases with otogenic DST following AOM. All but one of them presented initially with diplopia caused by otogenic hydrocephalus. In four cases the otologic complaints had already disappeared by the time of MRI confirmation of the diagnosis. Only one child was referred with severe otologic symptoms. Management included systemic antibiotics, short-term heparin anticoagulation and surgical decompression. In our cases, even after intensive i.v. antibiotic treatment, only surgery led to a significant improvement in the clinical condition.

Hybrid Multiscale Finite Volume method for two-phase flow in porous media

We present a novel hybrid (or multiphysics) algorithm, which couples pore-scale and Darcy descriptions of two-phase flow in porous media. The flow at the pore-scale is described by the Navier?Stokes equations, and the Volume of Fluid (VOF) method is used to model the evolution of the fluid?fluid interface. An extension of the Multiscale Finite Volume (MsFV) method is employed to construct the Darcy-scale problem. First, a set of local interpolators for pressure and velocity is constructed by solving the Navier?Stokes equations; then, a coarse mass-conservation problem is constructed by averaging the pore-scale velocity over the cells of a coarse grid, which act as control volumes; finally, a conservative pore-scale velocity field is reconstructed and used to advect the fluid?fluid interface. The method relies on the localization assumptions used to compute the interpolators (which are quite straightforward extensions of the standard MsFV) and on the postulate that the coarse-scale fluxes are proportional to the coarse-pressure differences. By numerical simulations of two-phase problems, we demonstrate that these assumptions provide hybrid solutions that are in good agreement with reference pore-scale solutions and are able to model the transition from stable to unstable flow regimes. Our hybrid method can naturally take advantage of several adaptive strategies and allows considering pore-scale fluxes only in some regions, while Darcy fluxes are used in the rest of the domain. Moreover, since the method relies on the assumption that the relationship between coarse-scale fluxes and pressure differences is local, it can be used as a numerical tool to investigate the limits of validity of Darcy's law and to understand the link between pore-scale quantities and their corresponding Darcy-scale variables.

Error estimate and control in the MSFV method for multiphase flow in porous media

n this paper the iterative MSFV method is extended to include the sequential implicit simulation of time dependent problems involving the solution of a system of pressure-saturation equations. To control numerical errors in simulation results, an error estimate, based on the residual of the MSFV approximate pressure field, is introduced. In the initial time steps in simulation iterations are employed until a specified accuracy in pressure is achieved. This initial solution is then used to improve the localization assumption at later time steps. Additional iterations in pressure solution are employed only when the pressure residual becomes larger than a specified threshold value. Efficiency of the strategy and the error control criteria are numerically investigated. This paper also shows that it is possible to derive an a-priori estimate and control based on the allowed pressure-equation residual to guarantee the desired accuracy in saturation calculation.

Calculation of correction factors for ionization chamber measurements with small fields in low-density media.

The quantity of interest for high-energy photon beam therapy recommended by most dosimetric protocols is the absorbed dose to water. Thus, ionization chambers are calibrated in absorbed dose to water, which is the same quantity as what is calculated by most treatment planning systems (TPS). However, when measurements are performed in a low-density medium, the presence of the ionization chamber generates a perturbation at the level of the secondary particle range. Therefore, the measured quantity is close to the absorbed dose to a volume of water equivalent to the chamber volume. This quantity is not equivalent to the dose calculated by a TPS, which is the absorbed dose to an infinitesimally small volume of water. This phenomenon can lead to an overestimation of the absorbed dose measured with an ionization chamber of up to 40% in extreme cases. In this paper, we propose a method to calculate correction factors based on the Monte Carlo simulations. These correction factors are obtained by the ratio of the absorbed dose to water in a low-density medium □D(w,Q,V1)(low) averaged over a scoring volume V₁ for a geometry where V₁ is filled with the low-density medium and the absorbed dose to water □D(w,QV2)(low) averaged over a volume V₂ for a geometry where V₂ is filled with water. In the Monte Carlo simulations, □D(w,QV2)(low) is obtained by replacing the volume of the ionization chamber by an equivalent volume of water, according to the definition of the absorbed dose to water. The method is validated in two different configurations which allowed us to study the behavior of this correction factor as a function of depth in phantom, photon beam energy, phantom density and field size.

The prediction of cardiovascular events is improved by combination of Framingham risk score and subclinical atherosclerosis imaging

Purpose: In primary prevention of cardiovascular disease (CVD), it is accepted that the intensity of risk factor treatment should be guided by the magnitude of absolute risk. Risk factors tools like Framingham risk score (FHS) or noninvasive atherosclerosis imaging tests are available to detect high risk subjects. However, these methods are imperfect and may misclassify a large number of individuals. The purpose of this prospective study was to evaluate whether the prediction of future cardiovascular events (CVE) can be improved when subclinical imaging atherosclerosis (SCATS) is combined with the FRS in asymptomatic subjects. Methods: Overall, 1038 asymptomatic subjects (413 women, 625 men, mean age 49.1±12.8 years) were assessed for their cardiovascular risk using the FRS. B-mode ultrasonography on carotid and femoral arteries was performed by two investigators to detect atherosclerotic plaques (focal thickening of intima-media > 1.2 mm) and to measure carotid intima-media thickness (C-IMT). The severity of SCATS was expressed by an ATS-burden Score (ABS) reflecting the number of the arterial sites with >1 plaques (range 0-4). CVE were defined as fatal or non fatal acute coronary syndrome, stroke, or angioplasty for peripheral artery disease. Results: during a mean follow-up of 4.9±3.1 years, 61 CVE were recorded. Event rates the rate of CVE increased significantly from 2.7% to 39.1% according to the ABS (p<0.001) and from 4% to 24.6% according to the quartiles of C-IMT. Similarly, FRS predicted CVE (p<0.001). When computing the angiographic markers of SCATS in addition of FRS, we observed an improvement of net reclassification rate of 16.6% (p< 0.04) for ABS as compared to 5.5% (p = 0.26) for C-IMT. Conclusion: these results indicate that the detection of subjects requiring more attention to prevent CVE can be significantly improved when using both FRS and SCATS imaging.