Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.

Media, investors and bubbles

We analysed the specific case of how information in the financial press influences economic bubbles. We found considerable flaws in the information market due to several factors: demand, the predominance of what are termed “irrational investors” (herding), and supply, which has the problem that the sources of information are biasedand feeds. A financial bubble is a deviation between real value of a financial asset and its persistent market price in time, which also has a speculative origin fed back by the illusion of the owners of these financial values, who will take benefits because of the future prices, which must be higher than the previous ones. The economical information in the media is submitting three problems. First of all, it is information generated by companies. In second place, the information circuit is fed back. A problem of informative independence becomes created, particularly serious in the case of the banks, which are very were as creditors. And in a third place, some informative biases are manifested for the companies of regulated sectors which are starring the economical information in the media.

Online teaching of inflammatory skin pathology by a French-speaking international university network

Introduction: Developments in technology, webbased teaching and whole slide imaging have broadened the teaching horizon in anatomic pathology. Creating online learning material including many types of media like radiologic images, videos, clinical and macroscopic photographs and whole slides imaging is now accessible to almost every university. Unfortunately, a major limiting factor to maintain and update the learning material is the amount of work, time and resources needed. In this perspective, a French national university network was initiated in 2011 to build mutualised online teaching pathology modules with clinical cases and tests. This network has been extended to an international level in 2012-2014 (Quebec, Switzerland and Ivory Coast). Method: One of the first steps of the international project was to build a learning module on inflammatory skin pathology intended for interns and residents of pathology and dermatology. A pathology resident from Quebec spent 6 weeks in France and Switzerland to develop the contents and build the module on an e-learning Moodle platform (http: //moodle.sorbonne-paris-cite.fr) under the supervision of two dermatopathologists (BV, MB). The learning module contains text, interactive clinical cases, tests with feedback, whole slides images (WSI), images and clinical photographs. For that module, the virtual slides are decentralized in 2 universities (Bordeaux and Paris 7). Each university is responsible of its own slide scanning, image storage and online display with virtual slide viewers. Results: The module on inflammatory skin pathology includes more than 50 web pages with French original content, tests and clinical cases, links to over 45 WSI and more than 50 micro and clinical photographs. The whole learning module is currently being revised by four dermatopathologists and two senior pathologists. It will be accessible to interns and residents in spring 2014. The experience and knowledge gained from that work will be transferred to the next international fellowship intern whose work will be aimed at creating lung and breast pathology learning modules. Conclusion: The challenges of sustaining a project of this scope are numerous. The technical aspect of whole-slide imaging and storage needs to be developed by each university or group. The content needs to be regularly updated, completed and its use and existence needs to be promoted by the different actors in pathology. Of the great benefits of that kind of project are the international partnerships and connections that have been established between numerous Frenchspeaking universities and pathologists with the common goals of promoting education in pathology and the use of technology including whole slide imaging. * The Moodle website is hosted by PRES Sorbonne Paris Cité, and financial supports for hardware have been obtained from UNF3S (http://www.unf3s.org/) and PRES Sorbonne Paris Cité. Financial support for international fellowships has been obtained from CFQCU (http://www.cfqcu.org/).

Fertility discourse in parental leave policies' media coverage : a frame analysis of French-speaking Swiss press articles from 1999 to 2009

This paper analyses the media coverage of parental leave policies (parental and paternity leaves) in Swiss French-speaking press articles from 1999 to 2009. Switzerland is one of the rare European countries which has no statutory parental or paternity leave. The aim is to describe the mediatisation of these policies and to analyse the arguments in favour and against their implementation. We investigate the status of a fertility frame - the mobilisation of discourse relating to fertility issues - among the various arguments used to justify or reject parental leave policies. We proceed with a content analysis of 579 press articles, as well as a frame analysis on a subset in which parental leave policies are the central theme (N=206). Results show that paternity leave is the predominant public issue addressed in the dataset. A mediatisation peak was reached in 2007, following an initiative of a member of the Federal executive to implement a short paternity leave. Parental leave policies are predominantly represented in a positive light. The main positive frame is economic, in which leaves are represented as serving the interests of companies. Involved fatherhood and gender equality are also frequently mentioned as positive frames. The fertility frame is only moderately used in articles covering Swiss news on paternity leaves. Conversely, the fertility frame is largely mobilised in articles covering parental leave in other countries. We discuss some interpretations of this discrepancy and suggest future avenues of research on parental leave policies in Switzerland.

Hierarchy of interactive functions in father-mother-baby three-way games

In developmental research, the family has mainly been studied through dyadic interaction. Three-way interactions have received less attention, partly because of their complexity. This difficulty may be overcome by distinguishing between four hierarchically embedded functions in three-way interactions: (1) participation (inclusion of all participants), (2) organization (partners keeping to their roles), (3) focalization (sharing a common focus) and (4) affective contact (being in tune). We document this hierarchical model on a sample of 80 families observed in the Lausanne Trilogue Play situation across four different sites. Hierarchy between functions was demonstrated by means of Guttman scalability coefficient. Given the importance of the child's development in a threesome, the pertinence of this model for family assessment is discussed.

Secondhand smoke exposure in children is associated with common carotid artery intima-media thickness

Objective: To investigate the association between common carotid artery intima-media thickness (cIMT) and exposure to secondhand smoke (SHS) in children. Methods: Data were available at baseline in the Quebec Adiposity and Lifestyle investigation in Youth (QUALITY) study, an ongoing longitudinal investigation of Caucasian children aged 8e10 years at cohort inception, who had at least one obese parent. Data on exposure to parents, siblings and friends smoking were collected in interviewer-administered child, and self-report parent questionnaires. Blood cotinine was measured with a high sensitivity ELISA. cIMTwas measured by ultrasound. The association between blood cotinine and cIMT was investigated in multivariable linear regression analyses controlling for age, body mass index, and child smoking status. Results: Mean (SD) cIMT (0.5803 (0.04602)) did not differ across age or sex. Overall 26%, 6% and 3% of children were exposed to parents, siblings and friends smoking, respectively. Cotinine ranged from 0.13 ng/ml to 7.38 ng/ml (median (IQR)¼0.18 ng/ml)). Multivariately, a 1 ng/ml increase in cotinine was associated with a 0.090 mm increase in cIMT (p¼0.034). Conclusion: In children as young as age 8e10 years, exposure to SHS relates to cIMT, a marker of pre-clinical atherosclerosis. Given the wide range of health effects of SHS, increased public health efforts are needed to reduced exposure among children in homes an private vehicles.

Diffusion in spatially and temporarily inhomogeneous media: Effects of turbulent mixing

We consider diffusion of a passive substance C in a phase-separating nonmiscible binary alloy under turbulent mixing. The substance is assumed to have different diffusion coefficients in the pure phases A and B, leading to a spatially and temporarily dependent diffusion ¿coefficient¿ in the diffusion equation plus convective term. In this paper we consider especially the effects of a turbulent flow field coupled to both the Cahn-Hilliard type evolution equation of the medium and the diffusion equation (both, therefore, supplemented by a convective term). It is shown that the formerly observed prolonged anomalous diffusion [H. Lehr, F. Sagués, and J.M. Sancho, Phys. Rev. E 54, 5028 (1996)] is no longer seen if a flow of sufficient intensity is supplied.

Diffusion in spatially and temporarily inhomogeneous media

In this paper we consider diffusion of a passive substance C in a temporarily and spatially inhomogeneous two-dimensional medium. As a realization for the latter we choose a phase-separating medium consisting of two substances A and B, whose dynamics is determined by the Cahn-Hilliard equation. Assuming different diffusion coefficients of C in A and B, we find that the variance of the distribution function of the said substance grows less than linearly in time. We derive a simple identity for the variance using a probabilistic ansatz and are then able to identify the interface between A and B as the main cause for this nonlinear dependence. We argue that, finally, for very large times the here temporarily dependent diffusion "constant" goes like t-1/3 to a constant asymptotic value D¿. The latter is calculated approximately by employing the effective-medium approximation and by fitting the simulation data to the said time dependence.

Front propagation in spatially modulated media

A two-dimensional reaction-diffusion front which propagates in a modulated medium is studied. The modulation consists of a spatial variation of the local front velocity in the transverse direction to that of the front propagation. We study analytically and numerically the final steady-state velocity and shape of the front, resulting from a nontrivial interplay between the local curvature effects and the global competition process between different maxima of the control parameter. The transient dynamics of the process is also studied numerically and analytically by means of singular perturbation techniques.