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.

Bilateral anterior pubic osteotomy in bladder exstrophy closure.

PURPOSE: We report our clinical experience with anterior pelvic osteotomy in 16 patients who underwent surgery for bladder exstrophy. The technique and its difficulties are discussed. MATERIALS AND METHODS: Anterior pelvic osteotomy of the superior ramus of the public bone is a simple and efficient method to facilitate symphyseal approximation and abdominal wall closure without or with low tension on the suture lines in neonates who undergo surgery for bladder exstrophy. Older children in whom surgery has been delayed can also benefit from this method. Compared to other methods of osteotomy in exstrophy surgery it does not require additional incisions, nor does the patient need to be repositioned on the operating table. A successful operation does not depend on any particular orthopedic skills and it can easily be done by the pediatric urologist. RESULTS: Immediate postoperative results regarding abdominal wall closure were excellent in all 16 patients. However, major postoperative complications developed in 2 patients. Despite antibiotic prophylaxis a severe soft tissue infection developed in 1 child, resulting in complete bladder dehiscence. In another patient an obturator nerve injury resulted in transient palsy, which resolved completely. While the first complication was not related to osteotomy, the second was osteotomy related. CONCLUSIONS: Bilateral superior ramotomy of the pubic bones is a new alternative, easily performed technique to optimize bladder exstrophy surgery in children.

Antagonistes de l'angiotensine II et cancer: un bilan rassurant [Angiotensin II receptor blockers (ARBs) and cancer: a reassuring balance].

The effects of drugs on new cancer and cancer-related death are a major concern. Recently, a meta-analysis raised the possibility that ARBs might have an adverse impact in this respect. This point of view was highly debated until the publication of two other meta-analyses which did not demonstrate any increased risk of new cancer occurrence as well as of cancer related-death with the use of ARBs in patients with hypertension, heart failure and/or nephropathy. This illustrates that the results of meta-analyses should be interpreted cautiously and critically in order to avoid biased conclusions. Overall the bulk of evidence today indicates that ARBs are not associated with an increased cancer risk.

Premature epiphyseal closure in an adolescent treated by retinoids for acne: An unusual cause of anterior knee pain

Introduction.- Retinoids are effective and widely used for the treatment of severe acne. Their use can be, however, associated with numerous side effects. For example, some rare cases of premature epiphyseal closure were reported.Observation.- A sixteen-year-old soccer player consulted for bilateral progressive anterior knee pain, since two months, evoking a femoro-patellar origin. After physiotherapy, the pain decreases on the right but remained on the left. The history taking brought out the use of isotretinoin for more than 6 months (0.5 mg/kg). Magnetic resonance imaging (MRI) findings showed an irregularity of the growth plate and an important metaphyso-epiphyseal oedema, more marked on the left. The diagnosis of retinoid induced premature ephysieal closure was retained. The treatment was stopped, with a resolution of symptoms within two months. The control MRI of the left knee present persisting small sequelar thumbprint-like growth plate lesion. Eighteen months later, neither limb-length discrepancy nor static disorder was noticed.Discussion.- Premature epiphyseal closure is a rare complication of retinoid treatment of acne. Retinoids induce an invasion of the growth plate by osteoclasts and a decrease in proteoglycans synthesis. The knee seems the most involved joint. The clinical presentation is aspecific, sometimes lightly symptomatic. A careful pharmacological history and an appropriate imaging are necessary. MRI is now the gold standard. It shows an irregularity of the growth plate with an oedema on both sides. In chronic phase, a thumbprint-like image may persist. The symptoms resolution arises in few weeks after the treatment interruption. A single case of static disorder was reported until now. The small size of the growth plate interruptions, insufficient to lead to a growth disorder if the medicament is stopped early enough, explains probably it. This complication being rare, a radiological follow-up of the young patients treated by retinoids is not proposed.

Dispersal routes reconstruction and the minimum cost arborescence problem.

We show that the dispersal routes reconstruction problem can be stated as an instance of a graph theoretical problem known as the minimum cost arborescence problem, for which there exist efficient algorithms. Furthermore, we derive some theoretical results, in a simplified setting, on the possible optimal values that can be obtained for this problem. With this, we place the dispersal routes reconstruction problem on solid theoretical grounds, establishing it as a tractable problem that also lends itself to formal mathematical and computational analysis. Finally, we present an insightful example of how this framework can be applied to real data. We propose that our computational method can be used to define the most parsimonious dispersal (or invasion) scenarios, which can then be tested using complementary methods such as genetic analysis.

Energy balance, physical activity, and thermogenic effect of feeding in premature infants.

In order to assess the contribution of the thermogenic effect of feeding and muscular activity to total energy expenditure, nine premature infants were studied for 2 consecutive days during which time repeated measurements of energy expenditure by indirect calorimetry were performed throughout the day, combined with a visual activity score based on body movement. The infants were growing at 16.6 +/- 4.0 g/kg/day (mean +/- SD) and received 110 +/- 8 kcal/kg/day metabolizable energy (milk formula) and 522 +/- 40 mgN/kg/day. Their total energy expenditure was 68 +/- 4 kcal/kg/day indicating that 41 +/- 7 kcal/kg/day was retained for growth. Based on the combination of energy + N balances it was estimated that 80% of the weight gain was fat-free tissue and 20% was fat tissue. The rate of energy expenditure measured minute-by-minute was significantly and linearly correlated with the activity score in both the premeal (r = 0.75;p less than 0.001) and the postmeal periods (r = 0.74; p less than 0.001) with no difference in the regression slope, but with a significant difference in intercept. In preset feeding schedules the latter allowed an estimation of the thermogenic effect without the confounding effect of activity. This was found to be 3.1 +/- 1.8% when expressed as a percentage of metabolizable energy intake. However when the "classical" approach was used as a comparison (integration of extra energy expenditure induced by the meal), the thermogenic effect was found to be greater, i.e. 9.5 +/- 3.8% of the meal's metabolizable energy, due to the superimposed effect of physical activity in the postprandial state.(ABSTRACT TRUNCATED AT 250 WORDS)

Recent misconceptions about the "database search problem": a probabilistic analysis using Bayesian networks

This paper analyses and discusses arguments that emerge from a recent discussion about the proper assessment of the evidential value of correspondences observed between the characteristics of a crime stain and those of a sample from a suspect when (i) this latter individual is found as a result of a database search and (ii) remaining database members are excluded as potential sources (because of different analytical characteristics). Using a graphical probability approach (i.e., Bayesian networks), the paper here intends to clarify that there is no need to (i) introduce a correction factor equal to the size of the searched database (i.e., to reduce a likelihood ratio), nor to (ii) adopt a propositional level not directly related to the suspect matching the crime stain (i.e., a proposition of the kind 'some person in (outside) the database is the source of the crime stain' rather than 'the suspect (some other person) is the source of the crime stain'). The present research thus confirms existing literature on the topic that has repeatedly demonstrated that the latter two requirements (i) and (ii) should not be a cause of concern.

ENaC-mediated alveolar fluid clearance and lung fluid balance depend on the channel-activating protease 1.

Sodium transport via epithelial sodium channels (ENaC) expressed in alveolar epithelial cells (AEC) provides the driving force for removal of fluid from the alveolar space. The membrane-bound channel-activating protease 1 (CAP1/Prss8) activates ENaC in vitro in various expression systems. To study the role of CAP1/Prss8 in alveolar sodium transport and lung fluid balance in vivo, we generated mice lacking CAP1/Prss8 in the alveolar epithelium using conditional Cre-loxP-mediated recombination. Deficiency of CAP1/Prss8 in AEC induced in vitro a 40% decrease in ENaC-mediated sodium currents. Sodium-driven alveolar fluid clearance (AFC) was reduced in CAP1/Prss8-deficient mice, due to a 48% decrease in amiloride-sensitive clearance, and was less sensitive to beta(2)-agonist treatment. Intra-alveolar treatment with neutrophil elastase, a soluble serine protease activating ENaC at the cell surface, fully restored basal AFC and the stimulation by beta(2)-agonists. Finally, acute volume-overload increased alveolar lining fluid volume in CAP1/Prss8-deficient mice. This study reveals that CAP1 plays a crucial role in the regulation of ENaC-mediated alveolar sodium and water transport and in mouse lung fluid balance.

The anti-apoptotic role of PPARbeta contributes to efficient skin wound healing.

PPARalpha and PPARbeta are expressed in the mouse epidermis during fetal development, but their expression progressively disappears after birth. However, the expression of PPARbeta is reactivated in adult mice upon proliferative stimuli, such as cutaneous injury. We show here that PPARbeta protects keratinocytes from growth factor deprivation, anoikis and TNF-alpha-induced apoptosis, by modulating both early and late apoptotic events via the Akt1 signaling pathway and DNA fragmentation, respectively. The control mechanisms involve direct transcriptional upregulation of ILK, PDK1, and ICAD-L. In accordance with the anti-apoptotic role of PPARbeta observed in vitro, the balance between proliferation and apoptosis is altered in the epidermis of wounded PPARbeta mutant mice, with increased keratinocyte proliferation and apoptosis. In addition, primary keratinocytes deleted for PPARbeta show defects in both cell-matrix and cell-cell contacts, and impaired cell migration. Together, these results suggest that the delayed wound closure observed in PPARbeta mutant mice involves the alteration of several key processes. Finally, comparison of PPARbeta and Akt1 knock-out mice reveals many similarities, and suggests that the ability of PPARbeta to modulate the Akt1 pathway has significant impact during skin wound healing.

New microfluidic-based sampling procedure for overcoming the hematocrit problem associated with dried blood spot analysis.

Hematocrit (Hct) is one of the most critical issues associated with the bioanalytical methods used for dried blood spot (DBS) sample analysis. Because Hct determines the viscosity of blood, it may affect the spreading of blood onto the filter paper. Hence, accurate quantitative data can only be obtained if the size of the paper filter extracted contains a fixed blood volume. We describe for the first time a microfluidic-based sampling procedure to enable accurate blood volume collection on commercially available DBS cards. The system allows the collection of a controlled volume of blood (e.g., 5 or 10 μL) within several seconds. Reproducibility of the sampling volume was examined in vivo on capillary blood by quantifying caffeine and paraxanthine on 5 different extracted DBS spots at two different time points and in vitro with a test compound, Mavoglurant, on 10 different spots at two Hct levels. Entire spots were extracted. In addition, the accuracy and precision (n = 3) data for the Mavoglurant quantitation in blood with Hct levels between 26% and 62% were evaluated. The interspot precision data were below 9.0%, which was equivalent to that of a manually spotted volume with a pipet. No Hct effect was observed in the quantitative results obtained for Hct levels from 26% to 62%. These data indicate that our microfluidic-based sampling procedure is accurate and precise and that the analysis of Mavoglurant is not affected by the Hct values. This provides a simple procedure for DBS sampling with a fixed volume of capillary blood, which could eliminate the recurrent Hct issue linked to DBS sample analysis.