La technologie de projection plasma sous basse pression a attiré l’attention de nombreux chercheurs comme une nouvelle technique qui permet d’établir un pont entre la projection thermique conventionnelle et le dépôt physique en phase vapeur. Ainsi, cette technologie étend les limites de projection thermique classique et augmente également la vitesse de croissance des dépôts par rapport aux procédés PVD ou CVD classiques. Cette technique peut évaporer totalement ou partiellement les poudres injectées et mener à la réalisation de revêtements à microstructure colonnaire et/ou plus denses, difficiles à réaliser avec des procédés de projection thermique conventionnels. La projection plasma de suspension a été effectuée pour la première fois sous basse pression. L’injection axiale de suspension avec une torche tri-cathodes permet d’augmenter l’échange enthalpique entre le jet de plasma et les poudres après l’évaporation du solvant. La spectroscopie à l’émission optique (OES) a été utilisée pour estimer la température électronique et vérifier l’existence de phase vapeur d’YSZ dans le jet de plasma. Finalement, des revêtements plus denses furent réalisés (comparés à ceux préparés par projection plasma de suspension à pression atmosphérique présentant des particules fondues, agglomérées et de la condensation de vapeur. Des tests de nano-indentation instrumentée ont été effectués sur la surface polie des dépôts réalisés. Les résultats montrent des valeurs de 5,8 GPa pour la dureté et 114,5 GPa pour le module d’élasticité, augmentant de 61% et 31%, respectivement, en comparaison avec les valeurs obtenues par SPS sous atmosphère ambiante. Les essais de projection de poudre YSZ agglomérée ont été réalisés avec une torche F4-VB dans le but de synthétiser une phase vapeur d’YSZ. On observe que les dépôts peuvent se former derrière les échantillons en céramique, sans vis-à-vis du plasma, par condensation de vapeur. En face de cette torche, des revêtements composites ont été obtenus par un mélange de poudres fondues et condensation de vapeur, simultanément. Cependant la quantité de phase vapeur est très faible dans le jet de plasma. Pour comprimer ce jet sous basse pression et afin d’améliorer l’échange d’enthalpie entre le jet de plasma et les poudres injectées, une buse rallongée a été mise en place sur la torche F4-VB. Les revêtements présentent ainsi une microstructure plus dense. Ceci est attribué à la haute vitesse des particules fondues vers le substrat suite à l’utilisation de la buse modifiée. Ce type de revêtement montre une valeur maximale de microdureté Vickers de 1273 Hv100 g. Par ailleurs, la réalisation de dépôts de carbures a été effectuée. Les résultats montrent la possibilité de former des carbures par projection plasma sous basse pression. Les revêtements composites (TiC/Ti) sont déposés par projection plasma réactif sous basse pression en utilisant le méthane comme gaz porteur. La température électronique Te calculée est d’environ 6200 K selon les résultats d’OES, ce qui est supérieur à la température d’ébullition du Ti et de TiC. Le revêtement de Ti pur présente une microstructure dense alors que TiC/Ti présente une microstructure lamellaire. Cependant, la quantité de TiC dans les revêtements est d’environ 20 vol.%. La microdureté Vickers, effectuée sur surface polie, a tendance à diminuer de 846±152 à 773±86 Hv100 g avec l’augmentation de la distance de projection. / As a new technology, the very low pressure plasma spraying has attracted attentions of many researchers, making it possible to establish a bridge between the conventional atmospheric plasma spraying (APS) and the vapor deposition (PVD or CVD). As a result, this new technology enlarges the limitation of APS and increases the deposition rate in comparison with the PVD or CVD. It is possible to evaporate partially the injected material and even evaporate completely and finally realize the columnar or dense coatings from the vapor or the mixture of vapor and liquid. The suspension plasma spraying is performed for the first time at low pressure. Taking consideration of the configuration of the three-cathode torch, the axial injection of the suspension is conducted which can increase the enthalpy change between the plasma jet and the sprayed material. The data of optical emission spectroscopy (OES) could be used to calculate the electron temperature and verify the existence of vapor of YSZ in the plasma jet. Finally, the dense coating was prepared by suspension plasma spraying at low pressure, which is composed of the melted particles, the agglomerated particles and the vapor deposition. The test of nano-indentation is conducted on the polished surface. It shows a value of 5.8 GPa for the microhardness and 114.5 GPa for the elastic modulus, increasing 61% and 31%, respectively, compared with the values obtained by SPS in the ambient atmosphere. In this study, another torch F4-VB is also conducted even of it has a low power in compared with that of O3CP torch. The powder feed rate is reduced to about 1.5 g·min-1 to achieve the vapor of YSZ taking the low power input of the torch into consideration. The columnar structure coating is realized from vapor deposition out of line of sight of projection upon the ceramic tubes. The composite structure coating is deposited by the mixture of melted particles and the vapor deposition simultaneously in front of this torch. But the quantity of vapor of YSZ is low in plasma jet. In order to compress the larger plasma jet and then improve the enthalpy change between the plasma jet and sprayed particles, an extended nozzle is prepared. It shows that the coating has a dense structure, which can be attributed to higher velocity of the melted powders. The coating shows a maximum value of microhardness Vickers up to 1273 Hv100 g. The composite coating of TiC/ Ti is realized by reactive plasma spraying using the methane as the carrier gas. The electron temperature Te is calculated to be 6200 K, which is over the boiling point of TiC and Ti. The coating Ti shows a dense structure and the composite coating TiC/Ti shows a lamellar structure. But the quantity of TiC in the composite coating is very low, about 20 vol.%. The Vickers microhardness is performed on the polished surface. It shows a decreasing tendency from 846 ± 152 to 773 ± 86 Hv100g with the increase of spraying distance. The tribological test is also implemented showing a high value of the coefficient of friction of 0.78 to 0.85, which can lead to a high abrasion. In order to synthesize a larger quantity of TiC in the composite coating, a higher power input torch should be put into action in the future.

Cette thèse répond au besoin de mieux maîtriser la durabilité des structures en béton armé exposées aux embruns d.eau de mer ou aux projections de sels de déverglaçage, en proposant une modélisation de la pénétration des ions chlorure dans les matériaux cimentaires lors des cycles de séchage- imbibition. Ce modèle tient compte de la cristallisation des sels, de la fixation des ions chlorure par les constituants solides de la pâte de ciment et plus généralement des interactions physico- chimiques entre les ions et la matrice cimentaire. Ces dernières sont basées sur l'équilibre chimique de la solution interstitielle avec les aluminates tricalciques (C3A), la portlandite (Ca(OH)2), et les sels de Friedel (3CaO.Al2O3.CaCl2.10H2O) et l'adsorption des ions chlorure sur les feuillets de C-S- H par substitution des l'ion hydroxyle. Cette description physico-chimique des interactions a permis de mettre en évidence que les isothermes de fixation des ions chlorure, mesurés expérimentalement, sont l'addition d'une partie dépendant seulement des teneurs initiales en composants solides issues de l'hydratation du ciment et d'une partie dépendant seulement de la nature intrinsèque des C-S-H produits. La pertinence de cette description est démontrée en comparant les résultats du modèle avec ceux de nombreux d'essais effectués sur plusieurs matériaux en conditions saturées issus de la littérature. L'extension de ce modèle d'interaction aux conditions non saturées n'introduit aucune donnée supplémentaire. Le reste du travail l'a démontré à travers la modélisation et la simulation de nombreux essais de séchage, d'imbibition ou de cycles d'imbibition-séchage sur des matériaux de construction ou cimentaires de caractéristiques différentes. Cette dernière partie a permis de mettre en évidence plusieurs phénomènes importants à prendre en compte tels que l'effet de la force ionique de la solution sur les courbes capillaires du matériau ou le colmatage des pores dû à la cristallisation du chlorure de sodium. Enfin le modèle a servi à expliquer les mécanismes physiques à l'origine de certains phénomènes observés lors du séchage et de l'imbibition de matériaux poreux en présence de sels.

L'objectif de cette thèse est de développer des méthodes numériques dans un code de simulation numérique direct contenant une méthode de suivi Lagrangien de phase disperse. Ces méthodes numériques ont permis d'optimiser le suivi Lagrangien parallèle des gouttes afin de pouvoir simuler de très grands domaines. Un algorithme de gestion des collisions entre les gouttes a aussi été développé et mis en place ainsi qu'une procédure permettant d'intégrer des corps immergés de formes quelconques. Ces méthodes numériques ont permis d'effectuer des simulations de référence et d'obtenir des statistiques détaillées de diverses configurations telles que le croisement de jets de gouttes qui se collisionnent ou bien la propagation d'un front de flamme dans un spray évaporant. La comparaison de ces caractéristiques avec les mêmes simulations où la prise en charge du spray était effectuée par une méthode Eulérienne multi-fluide ont été faites. Les très bonnes concordances qualitatives et quantitatives entre les prises en charge Lagrangiennes et Eulériennes ont permis de montrer les capacités du modèle Eulérien à rendre compte de phénomènes physiques complexes. Ce modèle Eulérien permettra de bien meilleurs scaling lors de son portage sur des architectures massivement parallèles.

Mes différentes activités de recherche qui s'insèrent dans une problématique axée autour du transport multi-échelle en milieu poreux sont détaillées ici au travers de 3 volets: (1) Influence des hétérogénéités sur les processus de transfert (2) Prise en compte des effets de couplage en transport réactif (3) Impact des processus biologiques sur l'hydrodynamique et le transport Un survol des défis qui restent à surmonter et des perspectives à venir de ce travail de recherche concluent le document.

Les modèles de transport réactif sont un outil basique pour la modélisation de l'intéraction entre les réactions chimiques et l'écoulement du fluide dans un milieu poreux. Nous présentons un modèle de transport réactif multi-espèces totalement réduit incluant des réactions cinétiques et en équilibre. Une formulation structurée ainsi que différentes approches numériques sont proposées. Les méthodes de décomposition de domaine offrent la possibilité de diviser des problèmes de grande taille dans des problèmes plus petits dont la solution se fait en parallèle. Partant d'un point de vue géométrique, nous présentons la classe des méthodes de Schwarz ayant prouvé une haute performance dans de nombreuses applications. Des questions quant à la réalisation d'une décomposition de domaine et des conditions de transmission au niveau discret sont traitées dans le contexte des volumes finis. Nous proposons et validons numériquement un schéma de volumes finis hybrides pour l'opérateur d'advection-diffusion étant particulièrement adapté à l'utilisation dans le contexte d'une décomposition de domaine. Nous étudions théoriquement et numériquement des méthodes de Schwarz relaxation d'ondes en détail pour un système de deux espèces couplées de type transport réactif avec des termes de couplage linéaire et non-linéaire. Des résultats qualifiant le problème comme bien posé ainsi que la convergence des méthodes de décomposition de domaine sont développés et la sensibilité du comportement de convergence de l'algorithme de Schwarz par rapport au terme de couplage est étudiée. Finalement, nous appliquons une méthode de Schwarz relaxation d'ondes au modèle de transport réactif multi-espèces présenté.

Mes activités de recherche et d'encadrement, depuis ma soutenance de thèse en octobre 2008, gravitent à la confluence de deux problématiques souvent considérées comme distinctes dont le couplage est plus qu'évident en informatique musicale : l'écoute artificielle et la programmation synchrone des signaux musicaux en temps réel. Le but est d'enrichir chaque domaine respectivement et de fournir des outils nécessaires aux artistes pour étendre leurs capacités de création musicale sur ordinateur (à celui de leurs homologues humains). Le premier aspect de tout système de musique interactive est l'acte de l'écoute ou formellement parlant, l'extraction d'informations musicaux en temps réel. J'ai concentré mes efforts sur les approches axées sur les applications telles que la détection de pitch polyphonique et les techniques d'alignement en temps réel d'une part. Puis, d'autre part, j'ai abordé des problèmes fondamentaux liés à la représentation des signaux. Pour le premier volet, nous avons proposé deux approches reconnues comme l'état de l'art et renommées par leur utilisation par la communauté artistique et leur reconnaissance au sein de la communauté MIR. Le deuxième volet théorique porte sur le problème fondamental de la quantification et la qualification du contenu de l'information arrivant progressivement dans un système en utilisant des méthodes de la géométrie de l'information. Le deuxième et plus récent aspect de mon travail se concentre sur la programmation synchrone réactive de la musique, couplée à des systèmes d'écoute (ou la formalisation de réactions dans les systèmes interactifs aux perceptions artificielles en temps réel). Notre approche est motivée par les pratiques actuelles en art numérique et les exigences d'évaluation en temps réel qu'elles demandent. L'incarnation majeure de cette approche est le système Antescofo doté d'une machine d'écoute et d'un langage dynamique et qui s'est imposé dans le monde entier dans le répertoire de la musique mixte, reconnu du public depuis 2009. L'écriture du temps, l'interaction et la tentative d'assurer leur exécution correcte en temps réel posent des problèmes difficiles pour l'informatique et les concepteurs de systèmes. Ce couplage fort fera l'objet de mon travail dans les années à venir sous la bannière des systèmes de musique Cyber-physique. Le couplage entre la machine d'écoute et des langages temps réel nécessite d'augmenter les approches actuelles avec des modèles formels de temps, dépassant les approches fonctionnelles qui dominent pour développer une formalisation des exigences de réaction et d'exécution. Il faudra en conséquence repenser nos méthodes usuelles de traitement audio multiforme et leurs modèles de calculs sous-jacents en ce qui concerne le temps.

Le stockage géologique du dioxyde de carbone (CO2) est une des options envisagées à moyen terme pour limiter les émissions de gaz à effet de serre. Or, le CO2 n'est pas un gaz inerte et aura tendance à acidifier l'eau en place dans les sites de stockage. Cette acidification de l'eau est alors susceptible de modifier la structure de cette roche ainsi que les propriétés de transport des espèces chimiques. Le but de cette étude est de quantifier l'impact du transport réactif sur la répartition d'une espèce chimique et sur la modification de la structure du milieu poreux de l'échelle du pore à celle du réservoir. Nous nous sommes focalisés dans cette étude sur le transport réactif monophasique d'une espèce dissoute aux temps longs. Pour ce faire, nous avons opté pour une approche multi-échelles considérant successi- vement (i) l'échelle locale, où les phénomènes d'écoulement, de réaction et de transport sont connus ; (ii) l'échelle du pore, où le transport réactif est représenté par des équations issues de l'écriture moyennée des équations locales ; (iii) l'échelle de Darcy (ou échelle de la carotte), où la structure de la roche est retranscrite par un réseau tridimensionnel de pores interconnectés par des canaux ; et (iv) l'échelle du réservoir, où les phénomènes physiques, au sein de chaque maille constituant le modèle réservoir, sont pris en compte par l'introduction de coefficients macroscopiques issus de l'étude de ces même phénomènes à l'échelle de Darcy, comme par exemple la perméabilité, la vitesse de réaction apparente, la vitesse apparente du soluté et sa dispersion.

Les techniques de planification de mouvement actuelles sont maintenant capables de résoudre des problèmes mettant en jeu des mécanismes complexes plongés au sein d'environnements encombrés. Néanmoins, l'adaptation de ces planificateurs à des scènes comprenant à la fois des obstacles statiques et des obstacles mobiles s'est avérée limitée jusqu'ici. Une des raisons en est le coût associé à la mise à jour des structures de données qui sont précalculées pour capturer la connexité de l'espace libre. Notre contribution principale concerne la proposition d'un nouveau planificateur capable de traiter des problèmes comprenant à la fois obstacles statiques et obstacles mobiles. Ce planificateur hybride combine deux grandes familles de techniques. D'une part les techniques dites PRM, initialement conçues pour résoudre des problèmes à requêtes multiples et que nous avons étendu à des problèmes de scènes dynamiques. D'autre part, de nouvelles techniques de diffusion, alors que celles-ci sont généralement dédiées aux problèmes simple requête ne nécessitant aucune opération de prétraitement. Les principaux développements accompagnant la construction de ce planificateur sont les suivants : - La proposition d'une architecture originale pour le planificateur dédié aux environnements changeants. Cette architecture inclut notamment plusieurs mécanismes dits d' "évaluation paresseuse" qui permettent de minimiser les test de collision et ainsi d'assurer de bonnes performances. - Le développement d'une nouvelle méthode de diffusion permettant de reconnecter localement certaines portions du réseau invalidées par la présence des obstacles mobiles. Cette méthode, appelée RRT à Domaine Dynamique correspond en fait une extension des planificateur bien connus à bases de RRTs. Un des intérêt propre à notre approche est d'équilibrer automatiquement deux comportements propres au planificateur : l'exploration vers des régions encore inconnues et l'affinage du modèle des régions de l'espac e déjà explorées. - Deux méthodes originales de création de réseaux cycliques qui servent à initialiser notre planificateur. La première assume que les obstacles mobiles sont confinés dans une région donnée, pour construire un réseau adapté aux différents types de changements de position possibles. La seconde est une méthode qui construit des réseaux appelés "réseaux de rétraction". A l'aide d'une structure de donnée de faible taille, cette structure parvient à capturer les différentes variétés de chemins de l'espace, à travers notamment chacune des classes d'homotopie de l'espace libre. Toutes ces méthodes sont implémentées au sein de la plate-forme de travail Move3D développée au LAAS-CNRS et sont évaluées sur différents types de systèmes mécaniques plongés au sein d'environnements 3D.

Le développement ainsi que l’amélioration des différentes techniques de purification sont des défis importants pour la chimie d’aujourd’hui. Certaines des méthodes actuelles, tel que le greffage d’un réactif sur un support solide permettant d’accéder à un produit pur par simple filtration du milieu, comportent toutefois certains inconvénients. En effet, les propriétés de solubilité de ces polymères rendent la mise en œuvre des réactions plus difficiles. C’est dans ce contexte que le groupe du Pr. Charette a rapporté l’utilisation de réactifs liés à un sel de tétraarylphosphonium (TAP). Ces sels peuvent être solubilisés dans un solvant tel que le dichlorométhane et aisément retirés du milieu réactionnel par précipitation à l’aide d’éther diéthylique (Chapitre 1). L’un des objectifs de cette thèse a donc été lié à la découverte de deux méthodes complémentaires qui, jusqu’à présent, sont considérées comme des méthodes de choix pour la synthèse des sels de TAP fonctionnalisés (Chapitre 2). L’une d’entre elles est utilisée par Soluphase inc., une entreprise qui commercialise ces sels de TAP. L’efficacité des sels en tant que support dans la synthèse de petites molécules a été démontrée lors de la synthèse d’un produit naturel, la (–)-coniine (Chapitre 3). L’isolement des intermédiaires synthétiques instables par simple précipitation à l’aide d’un support de TAP a permis de rendre cette synthèse plus efficace que celle déjà connue. Dans le deuxième volet de cette thèse, plusieurs problèmes reliés à la synthèse de dérivés cyclopropaniques 1,1-disubstitués ont été étudiés. Ces derniers font partie intégrale de plusieurs produits naturels et de médicaments. Cependant, leur formation par une réaction de cyclopropanation d’alcènes utilisant des réactifs diazoïques possédant deux groupements de type accepteur n’est pas tâche facile (Chapitre 4). En effet, cette réaction souffre d’un faible contrôle diastéréosélectif. Par le fait même, très peu de méthodologies de synthèse ont rapporté l’utilisation de ce type de réactifs diazoïques dans des réactions de cyclopropanation stéréosélectives. L’étude du mécanisme de la réaction de cyclopropanation catalysée au Rh(II) a proposé des indices favorisant un modèle ayant des précédents dans la littérature (Chapitre 5). Ces études nous ont mené à la découverte de la «capacité trans-directrice» du groupement amide lors des réactions de cyclopropanation d’oléfines. Nous avons donc utilisé cette propriété afin de résoudre plusieurs problèmes rencontrés dans la littérature. Nous avons montré qu’elle permet l’accès à des dérivés cyclopropaniques possédant deux groupements carboxyliques géminaux avec des sélectivités élevées (Chapitre 6). Ces produits étaient accessibles que par des séquences synthétiques nécessitant plusieurs étapes. De plus, nous avons démontrés que ces nouveaux dérivés cyclopropaniques sont des outils synthétiques fort utiles dans la synthèse de produits naturels d’intérêt biologique. Cette formidable «capacité trans-directrice» du groupement amide nous a permi de résoudre le problème de la synthèse asymétrique de dérivés carboxyliques α-cyano cyclopropaniques (Chapitre 7). De plus, ce projet nous a menées à la découverte de l’effet de divers additif achiraux permettant d’augmenter la sélectivité dans certaines réactions. Cette réaction possède une vaste étendue et l’utilité de ces nouveaux dérivés cyclopropaniques a été démontrée par plusieurs transformations de groupements fonctionnels. / The development of new purification techniques is an important challenge for today’s academic and industrial chemists. Present methods in which a given reagent is linked to a solid support to facilitate recovery still have some issues. Indeed, the solubility properties of these supports, mainly polymers, are not always reliable. In this context, the group of Pr. Charette has recently reported the use of a tetraarylphosphonium salt (TAP) as a solubility control group. TAPs are soluble in solvents such dichloromethane and can be quantitatively precipitated by the addition of diethyl ether (Chapter 1). However, the preparation of these TAPs still remains a synthetic challenge. One of the goal of this thesis lead to the discovery of two complementary methods that are considered among the most versatile ways to access functionalized TAPs (Chapter 2). One of these methods is utilized by Soluphase Inc., which markets these salts. The efficiency of the TAP moiety as a solubility control group in the synthesis of small molecules has been demonstrated by the synthesis of (–)-coniine (Chapter 3). All of the synthetic intermediates were isolated by a simple precipitation/filtration sequence. The TAP-supported synthesis has proven to be more efficient than the unsupported one. The second part of this thesis has focused on the synthesis of different 1,1-disubstituted cyclopropanes. Such cyclopropanes are a common motif in many natural products and synthetic drugs. However, the cyclopropanation reaction between an alkene and a metal carbene bearing two acceptor groups still remains a synthetic challenge (Chapter 4). Indeed, this reaction suffers from a low level of diastereocontrol. Therefore, little data have been reported to date on the stereoselective cyclopropanation with such carbenes. Studying the Rh(II)-catalyzed cyclopropanation of alkenes, we have brought strong evidence in favor of a postulated model of this reaction (Chapter 5). This study has also led to the discovery of the «trans-directing ability» of the amide group in Rh(II)-catalyzed cyclopropanation. We have thus utilized this «trans-directing ability» of the amide group in the Rh(II)-catalyzed cyclopropanation to solve different problems observed in the literature. We first showed that it could enable the synthesis of 1,1-dicarboxy cyclopropanes (Chapter 6). Multi-step syntheses were previously necessary to access such products. We have also demonstrated that these new cyclopropanes are useful synthetic tools for the rapid synthesis of a variety of natural products and biologically active molecules. We have also used this «trans-directing ability» of the amide group in Rh(II)-catalyzed cyclopropanation to enable the synthesis of 1-cyano-1-carboxy cyclopropanes (Chapter 7). Achiral additives were found to increase the selectivity of different metal catalyzed cyclopropanation reactions. The wide scope and synthetic utility of these new cyclopropanes has been further demonstrated by several functional group transformations.

Le projet de cette thèse concerne l'étude de l'impact de la variabilité initiale des réservoirs géologiques sur le stockage souterrain du CO₂. Nous nous sommes inspirés de réservoirs naturels : le Dogger carbonaté du Bassin de Paris et les réservoirs gréseux de Ravenscar (mer du Nord). Chacun des réservoirs modèles reconstruits a été défini par la liste de ses minéraux (calcite, dolomite, albite, clinochlore, k-feldspath, kaolinite, quartz ... par exemple pour les réservoirs carbonatés) et par ses propriétés pétrophysiques (porosité, perméabilité). La variabilité de ces paramètres a été simulée à l'aide de méthodes géostatistiques multivariables opérant sur leur portée, palier, et seuillage. Nous avons défini des paramètres opérationnels utiles pour résumer les propriétés principales du stockage (à savoir la quantité de carbone stockée dans le réservoir, la quantité de carbone précipité sous forme de dolomie, la quantité de carbone évacuée par l'aquifère mobile, la quantité de carbone contenue dans l'eau du réservoir ainsi que la variation de la porosité dans le réservoir). Des simulations déterministes de transfert réactif ont ensuite été effectuées en suivant différents scénarii de stockage : (Cas de l'injection du CO₂ initial dans le réservoir qui nous a permis de définir l'état initiale du réservoir, cas d'un aquifère mobile se déplaçant dans le réservoir contenant du CO₂ en saturation résiduel et le cas du déplacement d'une eau acide dans un réservoir contenant un aquifère neutre). Nous avons cherché par la suite à étudier l'influence de la variabilité initiale sur les paramètres opérationnels de stockage. Nous avons ainsi mis en évidence que, par rapport au cas homogène, le champ hétérogène favorise le stockage du carbone sous forme de dolomie et diminue la quantité de carbone entrainée par l'aquifère. Les effets des grandes portées de la pétrophysique et de l'aquifère mobile ont a également été étudié. Enfin, nous avons mis en place un modèle de plans d'expériences pour étudier l'effet de la variabilité de toutes les entrées sur les différentes réponses sélectionnées. L'effet de la dispersion numérique a été étudié et nous avons montré leur importance sur certains résultats. Les plans d'expériences ont mis en exergue l'importance des taches de la dolomie (associé ici à la portée du champ de dolomie seuillée) sur le carbone stocké.