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.

Measuring Party Positions and Issue Salience from Media Coverage: Introducing and Cross-Validating New Indicators

Recent studies have started to use media data to measure party positions and issue salience. The aim of this article is to compare and cross-validate this alternative approach with the more commonly used party manifestos, expert judgments and mass surveys. To this purpose, we present two methods to generate indicators of party positions and issue salience from media coverage: the core sentence approach and political claims analysis. Our cross-validation shows that with regard to party positions, indicators derived from the media converge with traditionally used measurements from party manifestos, mass surveys and expert judgments, but that salience indicators measure different underlying constructs. We conclude with a discussion of specific research questions for which media data offer potential advantages over more established methods.

Magnetism of isolated Mn12 single-molecule magnets detected by magnetic circular dichroism: Observation of spin tunneling with a magneto-optical technique

We report magnetic and magneto-optical measurements of two Mn12 single-molecule magnet derivatives isolated in organic glasses. Field-dependent magnetic circular dichroism (MCD) intensity curves (hysteresis cycles) are found to be essentially identical to superconducting quantum interference device magnetization results and provide experimental evidence for the potential of the optical technique for magnetic characterization. Optical observation of magnetic tunneling has been achieved by studying the decay of the MCD signal at weak applied magnetic field

Biased diffusion in confined media: Test of the Fick-Jacobs approximation and validity criteria

We study biased, diffusive transport of Brownian particles through narrow, spatially periodic structures in which the motion is constrained in lateral directions. The problem is analyzed under the perspective of the Fick-Jacobs equation, which accounts for the effect of the lateral confinement by introducing an entropic barrier in a one-dimensional diffusion. The validity of this approximation, based on the assumption of an instantaneous equilibration of the particle distribution in the cross section of the structure, is analyzed by comparing the different time scales that characterize the problem. A validity criterion is established in terms of the shape of the structure and of the applied force. It is analytically corroborated and verified by numerical simulations that the critical value of the force up to which this description holds true scales as the square of the periodicity of the structure. The criterion can be visualized by means of a diagram representing the regions where the Fick-Jacobs description becomes inaccurate in terms of the scaled force versus the periodicity of the structure.

Detection of metabolite induction in fungal co-cultures on solid media by high-throughput differential ultra-high pressure liquid chromatography-time-of-flight mass spectrometry fingerprinting.

Access to new biological sources is a key element of natural product research. A particularly large number of biologically active molecules have been found to originate from microorganisms. Very recently, the use of fungal co-culture to activate the silent genes involved in metabolite biosynthesis was found to be a successful method for the induction of new compounds. However, the detection and identification of the induced metabolites in the confrontation zone where fungi interact remain very challenging. To tackle this issue, a high-throughput UHPLC-TOF-MS-based metabolomic approach has been developed for the screening of fungal co-cultures in solid media at the petri dish level. The metabolites that were overexpressed because of fungal interactions were highlighted by comparing the LC-MS data obtained from the co-cultures and their corresponding mono-cultures. This comparison was achieved by subjecting automatically generated peak lists to statistical treatments. This strategy has been applied to more than 600 co-culture experiments that mainly involved fungal strains from the Fusarium genera, although experiments were also completed with a selection of several other filamentous fungi. This strategy was found to provide satisfactory repeatability and was used to detect the biomarkers of fungal induction in a large panel of filamentous fungi. This study demonstrates that co-culture results in consistent induction of potentially new metabolites.

Entropic transport in confined media: a challenge for computational studies in biological and soft-matter systems

Transport in small-scale biological and soft-matter systems typically occurs under confinement conditions in which particles proceed through obstacles and irregularities of the boundaries that may significantly alter their trajectories. A transport model that assimilates the confinement to the presence of entropic barriers provides an efficient approach to quantify its effect on the particle current and the diffusion coefficient. We review the main peculiarities of entropic transport and treat two cases in which confinement effects play a crucial role, with the appearance of emergent properties. The presence of entropic barriers modifies the mean first-passage time distribution and therefore plays a very important role in ion transport through micro- and nano-channels. The functionality of molecular motors, modeled as Brownian ratchets, is strongly affected when the motor proceeds in a confined medium that may constitute another source of rectification. The interplay between ratchet and entropic rectification gives rise to a wide variety of dynamical behaviors, not observed when the Brownian motor proceeds in an unbounded medium. Entropic transport offers new venues of transport control and particle manipulation and new ways to engineer more efficient devices for transport at the nanoscale.

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.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Positional stability of holographic optical traps

The potential of digital holography for complex manipulation of micron-sized particles with optical tweezers has been clearly demonstrated. By contrast, its use in quantitative experiments has been rather limited, partly due to fluctuations introduced by the spatial light modulator (SLM) that displays the kinoforms. This is an important issue when high temporal or spatial stability is a concern. We have investigated the performance of both an analog-addressed and a digitally-addressed SLM, measuring the phase fluctuations of the modulated beam and evaluating the resulting positional stability of a holographic trap. We show that, despite imparting a more unstable modulation to the wavefront, our digitally-addressed SLM generates optical traps in the sample plane stable enough for most applications. We further show that traps produced by the analog-addressed SLM exhibit a superior pointing stability, better than 1 nm, which is comparable to that of non-holographic tweezers. These results suggest a means to implement precision force measurement experiments with holographic optical tweezers (HOTs).

Growth regulators, culture media and antibiotics in the in vitro shoot regeneration from mature tissue of citrus cultivars

The objective of this study was to evaluate the effects of 6-benzylaminopurine (BAP) and α-naphthaleneacetic acid (NAA) combinations, basal media and beta-lactam antibiotics on in vitro organogenesis from mature stem segments of 'Pêra', 'Valência' and 'Bahia' sweet oranges and 'Cravo' rangpur lime. For induction of shoot regeneration, the segments of the four cultivars were placed on Murashige and Skoog (MS) medium containing the following BAP/NAA concentrations: 0.0/0.0; 0.25/0.0; 0.25/0.25; 0.5/0.0; 0.5/0.5; 1.0/0.0; 2.0/0.0; 2.0/0.25; 2.0/0.5; and 2.0/1.0 mg L-1. In order to test the influence of the culture media on shoot-bud induction, (MS), Murashige and Tucker (MT), and woody plant medium (WPM) formulations were evaluated, associated with the best combination of plant growth regulators obtained in the previous experiment. The influence of four beta-lactam antibiotics (timentin, cefotaxime sodium salt, meropenem trihydrate and augmentin) on shoot regeneration was determined. Better regeneration responses were achieved when internodal segments were cultured onto MS-based medium with 500 mg L-1 cefotaxime with the following BAP/NAA concentrations: 0.5 + 0.25 mg L-1 for 'Cravo', 1.0 + 0.25 mg L-1 for 'Valência' and 'Bahia', and 1.0 + 0.5 mg L-1 for 'Pêra'. Genotype, growth regulators, basal media and beta-lactam antibiotics affect the morphogenetic response in mature tissues of citrus.

Off-axis low-coherence interferometry for surface topology measurement.

In this communication we introduce a low or reduced coherence interferometry technique that can be used to retrieve surface topology on samples with high roughness. Moreover, we will show that the approach enables surface topology measurement also at the interface of so-called turbid media, where multiple scattering inside tissues can be a major issue, preventing accurate measurements.

Holographic optical tweezers combined with back-focal-plane displacement detection

A major problem with holographic optical tweezers (HOTs) is their incompatibility with laser-based position detection methods, such as back-focal-plane interferometry (BFPI). The alternatives generally used with HOTs, like high-speed video tracking, do not offer the same spatial and temporal bandwidths. This has limited the use of this technique in precise quantitative experiments. In this paper, we present an optical trap design that combines digital holography and back-focal-plane displacement detection. We show that, with a particularly simple setup, it is possible to generate a set of multiple holographic traps and an additional static non-holographic trap with orthogonal polarizations and that they can be, therefore, easily separated for measuring positions and forces with the high positional and temporal resolutions of laser-based detection. We prove that measurements from both polarizations contain less than 1% crosstalk and that traps in our setup are harmonic within the typical range. We further tested the instrument in a DNA stretching experiment and we discuss an interesting property of this configuration: the small drift of the differential signal between traps.

A novel optical intracellular imaging approach for potassium dynamics in astrocytes.

Astrocytes fulfill a central role in regulating K+ and glutamate, both released by neurons into the extracellular space during activity. Glial glutamate uptake is a secondary active process that involves the influx of three Na+ ions and one proton and the efflux of one K+ ion. Thus, intracellular K+ concentration ([K+]i) is potentially influenced both by extracellular K+ concentration ([K+]o) fluctuations and glutamate transport in astrocytes. We evaluated the impact of these K+ ion movements on [K+]i in primary mouse astrocytes by microspectrofluorimetry. We established a new noninvasive and reliable approach to monitor and quantify [K+]i using the recently developed K+ sensitive fluorescent indicator Asante Potassium Green-1 (APG-1). An in situ calibration procedure enabled us to estimate the resting [K+]i at 133±1 mM. We first investigated the dependency of [K+]i levels on [K+]o. We found that [K+]i followed [K+]o changes nearly proportionally in the range 3-10 mM, which is consistent with previously reported microelectrode measurements of intracellular K+ concentration changes in astrocytes. We then found that glutamate superfusion caused a reversible drop of [K+]i that depended on the glutamate concentration with an apparent EC50 of 11.1±1.4 µM, corresponding to the affinity of astrocyte glutamate transporters. The amplitude of the [K+]i drop was found to be 2.3±0.1 mM for 200 µM glutamate applications. Overall, this study shows that the fluorescent K+ indicator APG-1 is a powerful new tool for addressing important questions regarding fine [K+]i regulation with excellent spatial resolution.

The media accountability of independent regulatory agencies

Independent regulatory agencies (IRAs) are increasingly attracting academic and societal attention, as they represent the institutional cornerstone of the regulatory state and play a key role in policy-making. Besides the expected benefits in terms of credibility and efficiency, these regulators are said to bring about a 'democratic deficit', following their statutory separation from democratic institutions. Consequently, a 'multi-pronged system of control' is required. This article focuses on a specific component of this system, that is, the media. The goal is to determine whether media coverage of IRAs meets the necessary prerequisites to be considered a potential 'accountability forum' for regulators. The results of a comparison of two contrasted cases - the British and Swiss competition commissions - mostly support the expectations, because they show that media coverage of IRAs corresponds to that of the most relevant policy issues and follows the regulatory cycle. Furthermore, a systematic bias in media coverage can be excluded.

Textural controls on induced polarization and permeability in granular media

The spatial resolution visualized with hydrological models and the conceptualized images of subsurface hydrological processes often exceed resolution of the data collected with classical instrumentation at the field scale. In recent years it was possible to increasingly diminish the inherent gap to information from point like field data through the application of hydrogeophysical methods at field-scale. With regards to all common geophysical exploration techniques, electric and electromagnetic methods have arguably to greatest sensitivity to hydrologically relevant parameters. Of particular interest in this context are induced polarisation (IP) measurements, which essentially constrain the capacity of a probed subsurface region to store an electrical charge. In the absence of metallic conductors the IP- response is largely driven by current conduction along the grain surfaces. This offers the perspective to link such measurements to the characteristics of the solid-fluid-interface and thus, at least in unconsolidated sediments, should allow for first-order estimates of the permeability structure.¦While the IP-effect is well explored through laboratory experiments and in part verified through field data for clay-rich environments, the applicability of IP-based characterizations to clay-poor aquifers is not clear. For example, polarization mechanisms like membrane polarization are not applicable in the rather wide pore-systems of clay free sands, and the direct transposition of Schwarz' theory relating polarization of spheres to the relaxation mechanism of polarized cells to complex natural sediments yields ambiguous results.¦In order to improve our understanding of the structural origins of IP-signals in such environments as well as their correlation with pertinent hydrological parameters, various laboratory measurements have been conducted. We consider saturated quartz samples with a grain size spectrum varying from fine sand to fine gravel, that is grain diameters between 0,09 and 5,6 mm, as well as corresponding pertinent mixtures which can be regarded as proxies for widespread alluvial deposits. The pore space characteristics are altered by changing (i) the grain size spectra, (ii) the degree of compaction, and (iii) the level of sorting. We then examined how these changes affect the SIP response, the hydraulic conductivity, and the specific surface area of the considered samples, while keeping any electrochemical variability during the measurements as small as possible. The results do not follow simple assumptions on relationships to single parameters such as grain size. It was found that the complexity of natural occurring media is not yet sufficiently represented when modelling IP. At the same time simple correlation to permeability was found to be strong and consistent. Hence, adaptations with the aim of better representing the geo-structure of natural porous media were applied to the simplified model space used in Schwarz' IP-effect-theory. The resulting semi- empiric relationship was found to more accurately predict the IP-effect and its relation to the parameters grain size and permeability. If combined with recent findings about the effect of pore fluid electrochemistry together with advanced complex resistivity tomography, these results will allow us to picture diverse aspects of the subsurface with relative certainty. Within the framework of single measurement campaigns, hydrologiste can than collect data with information about the geo-structure and geo-chemistry of the subsurface. However, additional research efforts will be necessary to further improve the understanding of the physical origins of IP-effect and minimize the potential for false interpretations.¦-¦Dans l'étude des processus et caractéristiques hydrologiques des subsurfaces, la résolution spatiale donnée par les modèles hydrologiques dépasse souvent la résolution des données du terrain récoltées avec des méthodes classiques d'hydrologie. Récemment il est possible de réduire de plus en plus cet divergence spatiale entre modèles numériques et données du terrain par l'utilisation de méthodes géophysiques, notamment celles géoélectriques. Parmi les méthodes électriques, la polarisation provoquée (PP) permet de représenter la capacité des roches poreuses et des sols à stocker une charge électrique. En l'absence des métaux dans le sous-sol, cet effet est largement influencé par des caractéristiques de surface des matériaux. En conséquence les mesures PP offrent une information des interfaces entre solides et fluides dans les matériaux poreux que nous pouvons lier à la perméabilité également dirigée par ces mêmes paramètres. L'effet de la polarisation provoquée à été étudié dans différentes études de laboratoire, ainsi que sur le terrain. A cause d'une faible capacité de polarisation des matériaux sableux, comparé aux argiles, leur caractérisation par l'effet-PP reste difficile a interpréter d'une manière cohérente pour les environnements hétérogènes.¦Pour améliorer les connaissances sur l'importance de la structure du sous-sol sableux envers l'effet PP et des paramètres hydrologiques, nous avons fait des mesures de laboratoire variées. En détail, nous avons considéré des échantillons sableux de quartz avec des distributions de taille de grain entre sables fins et graviers fins, en diamètre cela fait entre 0,09 et 5,6 mm. Les caractéristiques de l'espace poreux sont changées en modifiant (i) la distribution de taille des grains, (ii) le degré de compaction, et (iii) le niveau d'hétérogénéité dans la distribution de taille de grains. En suite nous étudions comment ces changements influencent l'effet-PP, la perméabilité et la surface spécifique des échantillons. Les paramètres électrochimiques sont gardés à un minimum pendant les mesures. Les résultats ne montrent pas de relation simple entre les paramètres pétro-physiques comme par exemples la taille des grains. La complexité des media naturels n'est pas encore suffisamment représenté par les modèles des processus PP. Néanmoins, la simple corrélation entre effet PP et perméabilité est fort et consistant. En conséquence la théorie de Schwarz sur l'effet-PP a été adapté de manière semi-empirique pour mieux pouvoir estimer la relation entre les résultats de l'effet-PP et les paramètres taille de graines et perméabilité. Nos résultats concernant l'influence de la texture des matériaux et celles de l'effet de l'électrochimie des fluides dans les pores, permettront de visualiser des divers aspects du sous-sol. Avec des telles mesures géo-électriques, les hydrologues peuvent collectionner des données contenant des informations sur la structure et la chimie des fluides des sous-sols. Néanmoins, plus de recherches sur les origines physiques de l'effet-PP sont nécessaires afin de minimiser le risque potentiel d'une mauvaise interprétation des données.