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.

Using spatially distributed parameters and multi-response objective functions to solve parameterization of complex applications of semi-distributed hydrological models

Application of semi-distributed hydrological models to large, heterogeneous watersheds deals with several problems. On one hand, the spatial and temporal variability in catchment features should be adequately represented in the model parameterization, while maintaining the model complexity in an acceptable level to take advantage of state-of-the-art calibration techniques. On the other hand, model complexity enhances uncertainty in adjusted model parameter values, therefore increasing uncertainty in the water routing across the watershed. This is critical for water quality applications, where not only streamflow, but also a reliable estimation of the surface versus subsurface contributions to the runoff is needed. In this study, we show how a regularized inversion procedure combined with a multiobjective function calibration strategy successfully solves the parameterization of a complex application of a water quality-oriented hydrological model. The final value of several optimized parameters showed significant and consistentdifferences across geological and landscape features. Although the number of optimized parameters was significantly increased by the spatial and temporal discretization of adjustable parameters, the uncertainty in water routing results remained at reasonable values. In addition, a stepwise numerical analysis showed that the effects on calibration performance due to inclusion of different data types in the objective function could be inextricably linked. Thus caution should be taken when adding or removing data from an aggregated objective function.

Naevus du repli semi-lunaire de la conjonctive de découverte tardive [Conjunctival semi-lunar fold nevus of presumed late onset].

Conjunctival nevi are lesions observed commonly occasionnally and can be found later in a life. Such a lesion was observed on the conjunctival semilunar fold in a 76-year-old female patient. It was oval in shape, pigmented, and was mainly visible when the affected eye was in abduction. As no accurate anterior clinical history was available, surgical resection was performed. Hitopathological study disclosed a subepithelial nevus, which was completely resected. Regular follow-up was recommended, as is usual after surgical resection of a conjunctival pigmented tumor.

Deep learning via semi-supervised embedding

We show how nonlinear embedding algorithms popular for use with shallow semi-supervised learning techniques such as kernel methods can be applied to deep multilayer architectures, either as a regularizer at the output layer, or on each layer of the architecture. This provides a simple alternative to existing approaches to deep learning whilst yielding competitive error rates compared to those methods, and existing shallow semi-supervised techniques.

Avaliação de um modelo físico de produção de bovinos no semi-árido integrando caatinga, capim-buffel e leucena.: I. fase de cria

O estudo teve como objetivo avaliar um modelo físico do sistema de cria de bovinos para as áreas mais secas do Semi-Árido brasileiro, associando a pastagem natural com pastos cultivados tolerantes à seca. A caatinga foi pastejada no período verde, e o capim buffel no restante do ano. A suplementação com leucena ocorreu no período seco, sob pastejo direto e feno no cocho. Parâmetros de desempenho, monitorados de novembro de 1991 a outubro de 1995, mostraram uma taxa média de parição da ordem de 72,8% ao ano e taxas de mortalidade praticamente nulas. O peso vivo médio dos bezerros aos 205 dias de idade foi de 153,4 kg e a produ- ção de bezerros desmamados foi de 109,5 kg/matriz exposta/ano e de 25,5 kg/hectare/ano. Os resultados são considerados bastante expressivos, considerando-se que mais da metade da área era ocupada com pastagens nativas de caatinga e que a oferta de leucena foi limitada pela ocorrência de uma forte estiagem.

Uses of also in oral semi-informal German

In recent grammars and dictionaries also (`therefore, so, well¿) continues to be preferably presented as an adverb with a conclusive-consecutive connective function that essentially corresponds to its use in formal written German. Its function as a modal particle is documented, however, since the beginnings of what is known as Partikelforschung, though not all its uses have been systematically investigated contrasting oral and written German, either in mode or concept. In this article we analyse the uses of also in semi-informal oral interactions on the basis of empirical data (from a subsample of the VARCOM corpus). Specifically, we will analyse the presence and frequency of also at the beginning of a sentence or sequence, the functions it serves as a logical-semantic connector or discourse and interaction marker and the interrelations between these functions, in order to contrast these results with the description of also provided by current reference works.

SUPLEMENTAÇÃO ALIMENTAR DE CABRAS ANGLO-NUBIANAS NA ÉPOCA CHUVOSA, NA REGIÃO SEMI-ÁRIDA DO NORDESTE BRASILEIRO

Foi realizado um experimento para avaliar a influência da suplementação alimentar na produção de leite de cabras da raça Anglo-nubiana, durante a época chuvosa, na região semi-árida do Nordeste brasileiro. Vinte e quatro cabras foram mantidas em pastagem nativa rebaixada, distribuídas ao acaso, por tipo de parto e produção de leite, em três níveis de suplementação: N0 - sem suplementação, N1 -150 g de concentrado/animal/dia e N2 - 300 g de concentrado/animal/dia. Não houve diferença significativa (p>0,05) entre os níveis de suplementação. A produção de leite foi descrita por uma regressão múltipla em função do período de lactação em semana. A interação semana x níveis de suplementação foi significativa (p<0,01). Os resultados revelam que a suplementação com concentrado, na dose de 300 g/animal/dia, foi a melhor opção para se produzir leite com cabras Anglo-nubianas, durante a época chuvosa.

Efeito da calagem e inoculação de sabiá em solo da mata úmida e do semi-árido de Pernambuco

Realizou-se um experimento, usando solo álico da Zona da Mata e do Semi-Árido do Estado de Pernambuco, com o objetivo de avaliar o efeito da calagem (0, 3 e 6 t ha-1) e da inoculação das estirpes de rizóbio: NFB 539, NFB 577 e NFB 578, previamente selecionadas para sabiá (Mimosa caesalpiniaefolia) em meio de cultura ácido. Adicionaram-se tratamentos sem inoculação, sem e com adição de N mineral (100 kg ha-1), para fins comparativos. O experimento seguiu o esquema fatorial 2x3x5, no delineamento em blocos ao acaso, com quatro repetições. As plantas foram colhidas 110 dias após a emergência (DAE). A inoculação de rizóbio, em todos os níveis de calagem, mostrou efeito significativo nos parâmetros avaliados. No nível de 3 t ha-1 não houve efeito da calagem, porém a adição de 6 t ha-1 de calcário reduziu o peso de matéria seca, o N total na parte aérea, a nodulação e a atividade da nitrogenase. O baixo pH e Al3+ do solo não prejudicou a fixação do N2, e o crescimento das plantas que receberam o inóculo. Portanto, torna-se desnecessária a calagem no cultivo de sabiá em solos ácidos, quando usadas estirpes selecionadas visando à resistência a acidez.

Efeito de inoculante ectomicorrízico produzido por fermentação semi-sólida sobre o crescimento de Eucalyptus dunnii Maiden

O crescimento de Eucalyptus é favorecido por fungos ectomicorrízicos (FECM). De ocorrência comum nas plantações florestais, os FECM diferem, entretanto, em infectividade e eficiência. O objetivo deste estudo foi testar a compatibilidade e a eficiência do inoculante ectomicorrízico produzido por fermentação semi-sólida, em relação a Eucalyptus dunnii Maiden. Neste sentido, o inoculante de Pisolithus sp., produzido por fermentação semi-sólida em vermiculita-solução nutritiva, foi aplicado ao substrato turfa-vermiculita (30:70, v/v) nas concentrações de 0, 1, 3, 5 e 10% (v/v) que foi semeado com E. dunnii, e mantido em casa de vegetação. A inoculação proporcionou colonização radicular e aumento da matéria seca da parte aérea e do conteúdo de P em razão da concentração do inoculante. O peso da matéria seca das plantas com 1% de inoculante não diferiu do peso das testemunhas. A 3% e 5%, a matéria seca aumentou 62%. O maior valor foi obtido com 10% de inoculante, onde o peso das plantas foi 73% superior ao das testemunhas, e o teor de P, 130%. Resultados semelhantes foram obtidos em relação à raiz. Altura e diâmetro foram positivamente afetados pela inoculação, mas não a relação raiz/parte aérea. As doses mais eficientes de inoculante foram aquelas de 3% ou superiores. Este método de produção deverá ser testado com outros FECM, e sua eficiência, avaliada em viveiro e no campo.

A semi-grand canonical Monte Carlo simulation model for ion binding to ionizable surfaces: proton binding of carboxylated latex particles as a case study

In this paper, we present a computer simulation study of the ion binding process at an ionizable surface using a semi-grand canonical Monte Carlo method that models the surface as a discrete distribution of charged and neutral functional groups in equilibrium with explicit ions modelled in the context of the primitive model. The parameters of the simulation model were tuned and checked by comparison with experimental titrations of carboxylated latex particles in the presence of different ionic strengths of monovalent ions. The titration of these particles was analysed by calculating the degree of dissociation of the latex functional groups vs. pH curves at different background salt concentrations. As the charge of the titrated surface changes during the simulation, a procedure to keep the electroneutrality of the system is required. Here, two approaches are used with the choice depending on the ion selected to maintain electroneutrality: counterion or coion procedures. We compare and discuss the difference between the procedures. The simulations also provided a microscopic description of the electrostatic double layer (EDL) structure as a function of p H and ionic strength. The results allow us to quantify the effect of the size of the background salt ions and of the surface functional groups on the degree of dissociation. The non-homogeneous structure of the EDL was revealed by plotting the counterion density profiles around charged and neutral surface functional groups.

Virosome-Formulated Plasmodium falciparum AMA-1 & CSP Derived Peptides as Malaria Vaccine: Randomized Phase 1b Trial in Semi-Immune Adults & Children.

Background: This trial was conducted to evaluate the safety and immunogenicity of two virosome formulated malaria peptidomimetics derived from Plasmodium falciparum AMA-1 and CSP in malaria semi-immune adults and children.Methods: The design was a prospective randomized, double-blind, controlled, age-deescalating study with two immunizations. 10 adults and 40 children (aged 5-9 years) living in a malaria endemic area were immunized with PEV3B or virosomal influenza vaccine Inflexal (R) V on day 0 and 90.Results: No serious or severe adverse events (AEs) related to the vaccines were observed. The only local solicited AE reported was pain at injection site, which affected more children in the Inflexal (R) V group compared to the PEV3B group (p = 0.014). In the PEV3B group, IgG ELISA endpoint titers specific for the AMA-1 and CSP peptide antigens were significantly higher for most time points compared to the Inflexal (R) V control group. Across all time points after first immunization the average ratio of endpoint titers to baseline values in PEV3B subjects ranged from 4 to 15 in adults and from 4 to 66 in children. As an exploratory outcome, we found that the incidence rate of clinical malaria episodes in children vaccinees was half the rate of the control children between study days 30 and 365 (0.0035 episodes per day at risk for PEV3B vs. 0.0069 for Inflexal (R) V; RR = 0.50 [95%-CI: 0.29-0.88], p = 0.02).Conclusion: These findings provide a strong basis for the further development of multivalent virosomal malaria peptide vaccines.