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.

Origin of the relative afferent pupillary defect in optic tract lesions.

OBJECTIVE: To determine the percent decussation of pupil input fibers in humans and to explain the size and range of the log unit relative afferent pupillary defect (RAPD) in patients with optic tract lesions. DESIGN: Experimental study. PARTICIPANTS AND CONTROLS: Five patients with a unilateral optic tract lesion. METHODS: The pupil response from light stimulation of the nasal hemifield, temporal hemifield, and full field of each eye of 5 patients with a unilateral optic tract lesion was recorded using computerized binocular infrared pupillography. Six stimulus light intensities, separated by 0.5-log unit steps, were used; 12 stimulus repetitions were given for each stimulus condition. MAIN OUTCOME MEASURES: For each stimulus condition, the pupil response of each eye was characterized by plotting the mean pupil contraction amplitude as a function of stimulus light intensity. The percentage of decussating afferent pupillomotor input fibers was calculated from the ratio of the maximal pupil contractions elicited from each eye. The RAPD was determined pupillographically from full-field stimulation to each eye. RESULTS: In all patients, the pupil response from the functioning temporal hemifield ipsilateral to the tract lesion was greater than that from the functioning contralateral nasal hemifield. This temporal-nasal asymmetry increased with increasing stimulus intensity and was similar in hemifield and full-field stimuli, eventually saturating at maximal light intensity. The log unit RAPD did not correlate with the estimated percentage of decussating pupil fibers, which ranged from 54% to 67%. CONCLUSIONS: In patients with a unilateral optic tract lesion, the pupillary responses from full-field stimulation to each eye are the same as comparing the functioning temporal field with the functioning nasal field. The percentage of decussating fibers is reflected in the ratio of the maximal pupil contraction amplitudes resulting from stimulus input between the two eyes. The RAPD that occurs in this setting reflects the difference in light sensitivity between the intact temporal and nasal hemifields. Its magnitude does not correlate with the difference in the number of crossed and uncrossed axons, but its sidedness contralateral to the side of the optic tract lesion is consistent with the greater percentage of decussating pupillomotor input.

Diffusion in stationary flow from mesoscopic nonequilibrium thermodynamics

We analyze the diffusion of a Brownian particle in a fluid under stationary flow. By using the scheme of nonequilibrium thermodynamics in phase space, we obtain the Fokker-Planck equation that is compared with others derived from the kinetic theory and projector operator techniques. This equation exhibits violation of the fluctuation-dissipation theorem. By implementing the hydrodynamic regime described by the first moments of the nonequilibrium distribution, we find relaxation equations for the diffusion current and pressure tensor, allowing us to arrive at a complete description of the system in the inertial and diffusion regimes. The simplicity and generality of the method we propose makes it applicable to more complex situations, often encountered in problems of soft-condensed matter, in which not only one but more degrees of freedom are coupled to a nonequilibrium bath.

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.

A hybrid algorithm combining metaheuristic with Monte Carlo simulation for solving the Stochastic Flow Shop problem

In this paper, a hybrid simulation-based algorithm is proposed for the StochasticFlow Shop Problem. The main idea of the methodology is to transform the stochastic problem into a deterministic problem and then apply simulation to the latter. In order to achieve this goal, we rely on Monte Carlo Simulation and an adapted version of a deterministic heuristic. This approach aims to provide flexibility and simplicity due to the fact that it is not constrained by any previous assumption and relies in well-tested heuristics.

A hybrid algorithm combining metaheuristic with Monte Carlo simulation for solving the Stochastic Flow Shop problem

In this paper, a hybrid simulation-based algorithm is proposed for the StochasticFlow Shop Problem. The main idea of the methodology is to transform the stochastic problem into a deterministic problem and then apply simulation to the latter. In order to achieve this goal, we rely on Monte Carlo Simulation and an adapted version of a deterministic heuristic. This approach aims to provide flexibility and simplicity due to the fact that it is not constrained by any previous assumption and relies in well-tested heuristics.

Outcomes of optic nerve sheath decompression for nonarteritic ischemic optic neuropathy.

Efficacy of optic nerve sheath decompression (ONSD) in treating non-arteritic ischemic optic neuropathy (NAION) is not clear. We retrospectively analyzed the records of 91 patients with NAION, who were examined during a two-year period, and compared the final Snellen visual acuities of eyes treated with ONSD with those of eyes that did not have surgery. Seven of 18 eyes with ONSD (39%) demonstrated increased visual acuity of two or more lines; whereas 23 of 71 eyes without surgery (32%) had increased acuity. The ONSD group and no surgery group were further subdivided into eyes with progressive visual loss and nonprogressive visual loss. No statistically significant differences in visual outcome between groups were found. We did not find the high frequency of visual improvement that has been reported in some studies of ONSD for NAION.

Local re-inversion coronary MR angiography: arterial spin-labeling without the need for subtraction.

PURPOSE: To implement a double-inversion bright-blood coronary MR angiography sequence using a cylindrical re-inversion prepulse for selective visualization of the coronary arteries. MATERIALS AND METHODS: Local re-inversion bright-blood magnetization preparation was implemented using a nonselective inversion followed by a cylindrical aortic re-inversion prepulse. After an inversion delay that allows for in-flow of the labeled blood-pool into the coronary arteries, three-dimensional radial steady-state free-precession (SSFP) imaging (repetition/echo time, 7.2/3.6 ms; flip angle, 120 degrees, 16 profiles per RR interval; field of view, 360 mm; matrix, 512, twelve 3-mm slices) is performed. Coronary MR angiography was performed in three healthy volunteers and in one patient on a commercial 1.5 Tesla whole-body MR System. RESULTS: In all subjects, coronary arteries were selectively visualized with positive contrast. In addition, a middle-grade stenosis of the proximal right coronary artery was seen in one patient. CONCLUSION: A novel T1 contrast-enhancement strategy is presented for selective visualization of the coronary arteries without extrinsic contrast medium application. In comparison to former arterial spin-labeling schemes, the proposed magnetization preparation obviates the need for a second data set and subtraction.

Trust region versus line search for computing the optical flow

We consider the numerical treatment of the optical flow problem by evaluating the performance of the trust region method versus the line search method. To the best of our knowledge, the trust region method is studied here for the first time for variational optical flow computation. Four different optical flow models are used to test the performance of the proposed algorithm combining linear and nonlinear data terms with quadratic and TV regularization. We show that trust region often performs better than line search; especially in the presence of non-linearity and non-convexity in the model.

Free-breathing whole-heart coronary MRA with 3D radial SSFP and self-navigated image reconstruction.

Respiratory motion is a major source of artifacts in cardiac magnetic resonance imaging (MRI). Free-breathing techniques with pencil-beam navigators efficiently suppress respiratory motion and minimize the need for patient cooperation. However, the correlation between the measured navigator position and the actual position of the heart may be adversely affected by hysteretic effects, navigator position, and temporal delays between the navigators and the image acquisition. In addition, irregular breathing patterns during navigator-gated scanning may result in low scan efficiency and prolonged scan time. The purpose of this study was to develop and implement a self-navigated, free-breathing, whole-heart 3D coronary MRI technique that would overcome these shortcomings and improve the ease-of-use of coronary MRI. A signal synchronous with respiration was extracted directly from the echoes acquired for imaging, and the motion information was used for retrospective, rigid-body, through-plane motion correction. The images obtained from the self-navigated reconstruction were compared with the results from conventional, prospective, pencil-beam navigator tracking. Image quality was improved in phantom studies using self-navigation, while equivalent results were obtained with both techniques in preliminary in vivo studies.

Reduced microwave losses of YBa2Cu3O7_thin films on electro-optic LiNbO3 crystals

We report on the growth of epitaxial YBa2Cu3O7 thin films on X-cut LiNbO3 single crystals. The use of double CeO2/YSZ buffer layers allows a single in-plane orientation of YBa2Cu3O7, and results in superior superconducting properties. In particular, surface resistance Rs values of 1.4 m¿ have been measured at 8 GHz and 65 K. The attainment of such low values of Rs constitutes a key step toward the incorporation of high Tc materials as electrodes in photonic and acoustic devices.

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.

Milk flow, teat morphology and subclinical mastitis prevalence in Gir cows

The aim of this work was to evaluate the association between milk flow, teat morphological measurements and subclinical mastitis prevalence in Gir cows. Eighty cows in the 2nd and 3rd lactations, with 90 to 200 days of lactation, were divided according to milk flow during milking into fast or slow groups. Teat morphometry was assessed by ultrasound scanning of the right anterior teat and external measurements. Milk samples were collected for somatic cells count (SCC) and microbiological culture. The effect of milk flow during milking was evaluated by analysis of variance of milk yield, SCC, morphometry and external measurements. The association of morphometry and external measurements of the teats with the SCC and microorganisms found in milk were analysed. Milk flow was significantly correlated to milk production. Gir cows with slower milk flow had longer teat canal and greater milk yield, in comparison to cows with fast milk flow. Teat-end to floor distance influenced SCC of Gir cows. Prevalence of subclinical mastitis and the type of mastitis-causing pathogens were not affected by milk flow during milking