Stoneley wave modeling in heterogeneous porous media with viscous pore fluids

We implemented Biot-type porous wave equations in a pseudo-spectral numerical modeling algorithm for the simulation of Stoneley waves in porous media. Fourier and Chebyshev methods are used to compute the spatial derivatives along the horizontal and vertical directions, respectively. To prevent from overly short time steps due to the small grid spacing at the top and bottom of the model as a consequence of the Chebyshev operator, the mesh is stretched in the vertical direction. As a large benefit, the Chebyshev operator allows for an explicit treatment of interfaces. Boundary conditions can be implemented with a characteristics approach. The characteristic variables are evaluated at zero viscosity. We use this approach to model seismic wave propagation at the interface between a fluid and a porous medium. Each medium is represented by a different mesh and the two meshes are connected through the above described characteristics domain-decomposition method. We show an experiment for sealed pore boundary conditions, where we first compare the numerical solution to an analytical solution. We then show the influence of heterogeneity and viscosity of the pore fluid on the propagation of the Stoneley wave and surface waves in general.

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.

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.

A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work we perform direct numerical simulations of drainage by solving Navier-Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Synthesis of medium-chain-length polyhydroxyalkanoates in arabidopsis thaliana using intermediates of peroxisomal fatty acid beta-oxidation.

Polyhydroxyalkanoate (PHA) is a family of polymers composed primarily of R-3-hydroxyalkanoic acids. These polymers have properties of biodegradable thermoplastics and elastomers. Medium-chain-length PHAs (MCL-PHAs) are synthesized in bacteria by using intermediates of the beta-oxidation of alkanoic acids. To assess the feasibility of producing MCL-PHAs in plants, Arabidopsis thaliana was transformed with the PhaC1 synthase from Pseudomonas aeruginosa modified for peroxisome targeting by addition of the carboxyl 34 amino acids from the Brassica napus isocitrate lyase. Immunocytochemistry demonstrated that the modified PHA synthase was appropriately targeted to leaf-type peroxisomes in light-grown plants and glyoxysomes in dark-grown plants. Plants expressing the PHA synthase accumulated electron-lucent inclusions in the glyoxysomes and leaf-type peroxisomes, as well as in the vacuole. These inclusions were similar to bacterial PHA inclusions. Analysis of plant extracts by GC and mass spectrometry demonstrated the presence of MCL-PHA in transgenic plants to approximately 4 mg per g of dry weight. The plant PHA contained saturated and unsaturated 3-hydroxyalkanoic acids ranging from six to 16 carbons with 41% of the monomers being 3-hydroxyoctanoic acid and 3-hydroxyoctenoic acid. These results indicate that the beta-oxidation of plant fatty acids can generate a broad range of R-3-hydroxyacyl-CoA intermediates that can be used to synthesize MCL-PHAs.

Postmortem computed tomography angiography, contrast medium administration and toxicological analyses in urine.

Postmortem angiography methods that use water soluble or lipid soluble liquid contrast compounds may potentially modify the composition of fluid-based biological samples and thus influence toxicological findings. In this study, we investigated whether toxicological investigations performed in urine collected prior to and post angiography using Angiofil? mixed with paraffin oil are characterized by different qualitative or quantitative results. In addition, we studied whether diluting samples with 1% and 3% contrast medium solution may modify molecule concentration. A postmortem angiography group consisting of 50 cases and a postmortem group without angiography consisting of 50 cases were formed. In the first group, toxicological investigations were performed in urine samples collected prior to and post angiography as well as in undiluted and diluted samples. In the second group, analyses were performed in undiluted and diluted urine, bile, gastric content, cerebrospinal and pericardial fluids collected during autopsy. The preliminary results indicate that differences may be observed between urine samples collected prior to and post angiography in the number of identified molecules in relation to specific cases. Analyses performed in diluted samples failed to reveal differences that might potentially alter the interpretation of toxicological results in all analyzed specimens for nearly all molecules, except for tetrahydrocannabinol and its metabolites. Though these findings suggest that toxicology might be effectively performed, in very special cases and for a large number of molecules, in biological samples collected after angiography, it remains recommendable to collect biological fluids for toxicology prior to contrast medium injection.

Simultaneous cell morphometry and refractive index measurement with dual-wavelength digital holographic microscopy and dye-enhanced dispersion of perfusion medium.

Digital holographic microscopy (DHM) allows optical-path-difference (OPD) measurements with nanometric accuracy. OPD induced by transparent cells depends on both the refractive index (RI) of cells and their morphology. This Letter presents a dual-wavelength DHM that allows us to separately measure both the RI and the cellular thickness by exploiting an enhanced dispersion of the perfusion medium achieved by the utilization of an extracellular dye. The two wavelengths are chosen in the vicinity of the absorption peak of the dye, where the absorption is accompanied by a significant variation of the RI as a function of the wavelength.

Television before TV : a transnational history of an experimental medium on display, 1928 - 1939

Dès le milieu des années 1920 et le développement des premiers prototypes fonctionnels, la télévision se fait objet d'exposition. La nouvelle technologie est montrée à des foires industrielles et dans des grands magasins, à des expositions universelles et nationales : exposée, la télévision devient un mass media avant même qu'elle ne diffuse des émissions régulières. En étudiant les présentations publiques de la télévision en Allemagne, en Grande-Bretagne et aux Etats-Unis entre 1928 - date de l'ouverture des premières expositions annuelles dans les trois pays - et 1939 - moment de l'interruption de celles-ci suite à l'éclatement de la Deuxième Guerre mondiale - Television before TV montre que les lieux où le médium rencontre son premier public constituent les mêmes espaces qui, rétrospectivement, permettent de comprendre son avènement. Cette recherche propose ainsi de nouvelles pistes pour l'histoire des médias en interrogeant la définition de la télévision en tant que médium du privé. Elle rappelle l'importance historique de l'espace public pour la « vision à distance » et, dépassant la perspective nationale, reconstitue les différences et similarités techniques, médiatiques et institutionnelles de la télévision allemande, américaine et britannique dans l'entre-deux-guerres. Développée dans un contexte qui est simultanément nationaliste et nourri d'échanges et de compétitions internationales, la télévision sert de comparatif dans les rivalités entre pays, mais stimule également la circulation de savoir-faire et de personnes. Sa forte valeur symbolique comme emblème du progrès scientifique et, dans le cas allemand, comme preuve de la modernité du régime national-socialiste se traduit dans les salles d'expositions par des scénographies souvent innovatrices, dont les photographies et descriptions se propagent au-delà de la presse locale. Au centre d'un réseau discursif et représentationnel, les expositions jouent alors un rôle essentiel pour la construction sociale, politique et culturelle de la télévision. Remplaçant l'analyse de textes audiovisuels par l'étude de la machine exposée, ma recherche reformule la question bazinienne de qu'est-ce que la télévision ? pour demander où le médium se situe-t-il ? Cette prémisse méthodologique permet de contourner une évaluation qualitative de la technologie qui déplorerait ses multiples imperfections, et éclaire d'un nouveau jour des dispositifs télévisuels qui, même sans programme, intègrent un univers voué à la culture de consommation et du loisir.

Efficacité a court et moyen termes des traitements d'urgence-crise pour couples et families: résultats et perspectives [Short and medium-term effectiveness of the emergency-crisis treatments for couples and families: results and perspectives].

From 2007 to 2010, the emergency-crisis unit of the Couple and Family Consultation Unit -UCCF (West Psychiatric Service, Prangins Psychiatric Hospital, Psychiatric Department of CHUV) has carried out a research about the relevance and usefulness of emergency-crisis, systemic-oriented treatments, for deeply distressed couples and families. Besides epidemiologic data, we present results demonstrating the efficiency of those treatments, both at short-term and at a one year's range. The global impact of such treatments in terms of public health, but also economical issues, make us believe that they should be fully included in the new trend of psychiatric ambulatory care, into the social net.

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.

Increased flow of fatty acids toward beta-oxidation in developing seeds of Arabidopsis deficient in diacylglycerol acyltransferase activity or synthesizing medium-chain-length fatty acids.

Synthesis of polyhydroxyalkanoates (PHAs) from intermediates of fatty acid beta-oxidation was used as a tool to study fatty acid degradation in developing seeds of Arabidopsis. Transgenic plants expressing a peroxisomal PHA synthase under the control of a napin promoter accumulated PHA in developing seeds to a final level of 0. 06 mg g(-1) dry weight. In plants co-expressing a plastidial acyl-acyl carrier protein thioesterase from Cuphea lanceolata and a peroxisomal PHA synthase, approximately 18-fold more PHA accumulated in developing seeds. The proportion of 3-hydroxydecanoic acid monomer in the PHA was strongly increased, indicating a large flow of capric acid toward beta-oxidation. Furthermore, expression of the peroxisomal PHA synthase in an Arabidopsis mutant deficient in the enzyme diacylglycerol acyltransferase resulted in a 10-fold increase in PHA accumulation in developing seeds. These data indicate that plants can respond to the inadequate incorporation of fatty acids into triacylglycerides by recycling the fatty acids via beta-oxidation and that a considerable flow toward beta-oxidation can occur even in a plant tissue primarily devoted to the accumulation of storage lipids.

Use of quantum dots in aqueous solution to detect blood fingermarks on non-porous surfaces.

A new and original reagent based on the use of highly fluorescent cadmium telluride (CdTe) quantum dots (QDs) in aqueous solution is proposed to detect weak fingermarks in blood on non-porous surfaces. To assess the efficiency of this approach, comparisons were performed with one of the most efficient blood reagents on non-porous surfaces, Acid Yellow 7 (AY7). To this end, four non-porous surfaces were studied, i.e. glass, transparent polypropylene, black polyethylene, and aluminium foil. To evaluate the sensitivity of both reagents, sets of depleted fingermarks were prepared, using the same finger, initially soaked with blood, which was then successively applied on the same surface without recharging it with blood or latent secretions. The successive marks were then cut in halves and the halves treated separately with each reagent. The results showed that QDs were equally efficient to AY7 on glass, polyethylene and polypropylene surfaces, and were superior to AY7 on aluminium. The use of QDs in new, sensitive and highly efficient latent and blood mark detection techniques appears highly promising. Health and safety issues related to the use of cadmium are also discussed. It is suggested that applying QDs in aqueous solution (and not as a dry dusting powder) considerably lowers the toxicity risks.

The Peroxisomal Enzyme L-PBE Is Required to Prevent the Dietary Toxicity of Medium-Chain Fatty Acids.

Specific metabolic pathways are activated by different nutrients to adapt the organism to available resources. Although essential, these mechanisms are incompletely defined. Here, we report that medium-chain fatty acids contained in coconut oil, a major source of dietary fat, induce the liver ω-oxidation genes Cyp4a10 and Cyp4a14 to increase the production of dicarboxylic fatty acids. Furthermore, these activate all ω- and β-oxidation pathways through peroxisome proliferator activated receptor (PPAR) α and PPARγ, an activation loop normally kept under control by dicarboxylic fatty acid degradation by the peroxisomal enzyme L-PBE. Indeed, L-pbe(-/-) mice fed coconut oil overaccumulate dicarboxylic fatty acids, which activate all fatty acid oxidation pathways and lead to liver inflammation, fibrosis, and death. Thus, the correct homeostasis of dicarboxylic fatty acids is a means to regulate the efficient utilization of ingested medium-chain fatty acids, and its deregulation exemplifies the intricate relationship between impaired metabolism and inflammation.