Open system behavior is predicated on a fundamental relationship between the timescale over which mass is transported and the timescale over which it is chemically transformed. This relationship describes the basis for the multidisciplinary field of reactive transport (RT). In the 20 years since publication of Review in Mineralogy and Geochemistry volume 34: Reactive Transport in Porous Media, RT principles have expanded beyond early applications largely based in contaminant hydrology to become broadly utilized throughout the Earth Sciences. RT is now employed to address a wide variety of natural and engineered systems across diverse spatial and temporal scales, in tandem with advances in computational capability, quantitative imaging and reactive interface characterization techniques. The present volume reviews the diversity of reactive transport applications developed over the past 20 years, ranging from the understanding of basic processes at the nano- to micrometer scale to the prediction of Earth global cycling processes at the watershed scale. Key areas of RT development are highlighted to continue advancing our capabilities to predict mass and energy transfer in natural and engineered systems.

On s'intéresse dans ce travail à la simulation des écoulements diphasiques. Différents modèles, tous hyperboliques, sont considérés suivant les configurations étudiées. Dans un premier temps, plusieurs schémas Volumes Finis sont comparés pour l'approximation du modèle HEM (Homogeneous Equilibrium Model), notamment en présence de faibles densités. Ensuite on démontre l'existence et l'unicité de la solution faible entropique d'une loi de conservation scalaire gouvernant l'évolution de la saturation d'un écoulement diphasique dans un milieu poreux. On propose alors deux schémas Volumes Finis tenant compte du caractère résonnant de cette équation. La troisième partie concerne les écoulements en eaux peu profondes et l'approximation des termes sources raides. Une méthode pemettant le maintien d'états au repos ainsi que le recouvrement et l'apparition de zones sèches, est présentée et comparée aux méthodes habituellement utilisées dans l'industrie. Enfin, une classe de modèles hyperboliques non conservatifs se basant sur l'approche bifluide à deux vitesses et deux pressions est proposée. Une étude des solutions discontinues du système convectif permet d'exhiber une classe de fermetures sur la vitesse interfaciale et sur la pression interfaciale, tout en permettant de définir de manière unique les produits non conservatifs. L'approximation se fait à l'aide d'une méthode de "splitting" d'opérateur. On utilise deux schémas Volumes Finis, le schéma de Rusanov et le schéma de Godunov approché VFRoe-ncv pour l'étape de convection. Plusieurs cas tests sont présentés et commentés : tubes à choc, conditions limites de paroi, robinet d'eau, sédimentation

Cette thèse a pour but de présenter les travaux liés à l'étude numérique des écoulements critiques diphasiques dans la tuyère primaire d'un éjecteur utilisé pour améliorer le rendement d'un cycle de réfrigération. En premier lieu, les bases d'un modèle d'écoulement multiphasique ainsi que les termes de source associés sont présentés. Les développements entourant les points singuliers ainsi que les conditions critiques sont également présentés. Deux schémas numériques unidimensionnels sont proposés pour l'étude d'un écoulement critique diphasique particulaire. Les comparaisons des flux massiques critiques ainsi que des profils de pression numériques et expérimentaux montrent que le modèle représente bien l'évolution de l'écoulement diphasique. L'existence d'un point critique situé à un endroit différent du col géométrique a été vérifiée. L'étude numérique présente également une validation des conditions critiques utilisées ainsi que les effets des déséquilibres sur la position du point critique et sur l'écoulement. Des travaux supplémentaires sont cependant requis pour mieux comprendre le phénomène entourant les chocs en écoulements diphasiques et pour accroître la plage de taux de vide couverte par les schémas numériques.

On s'interesse dans ce travail à la simulation des écoulements diphasiques. Différents modèles, tous hyperboliques, sont considérés suivant les configurations étudiées. Dans un premier temps, plusieurs schémas Volumes Finis sont comparés pour l'approximation du modèle HEM (Homogeneous Equilibrium Model), notamment en présence de faibles densités. Ensuite on démontre l'existence et l'unicité de la solution faible entropique d'une loi de conservation scalaire gouvernant l'évolution de la saturation d'un écoulement diphasique dans un milieu poreux. On propose alors deux schémas Volumes Finis tenant compte du caractère résonnant de cette équation. La troisième partie concerne les écoulements en eaux peu profondes et l'approximation des termes sources raides. Une méthode pemettant le maintien d'états au repos ainsi que le recouvrement et l'apparition de zones sèches, est présentée et comparée aux méthodes habituellement utilisées dans l'industrie. Enfin, une classe de modèles hyperboliques non conservatifs se basant sur l'approche bifluide à deux vitesses et deux pressions est proposée. Une étude des solutions discontinues du système convectif permet d'exhiber une classe de fermetures sur la vitesse interfaciale et sur la pression interfaciale, tout en permettant de définir de manière unique les produits non conservatifs. L'approximation se fait à l'aide d'une méthode de "splitting" d'opérateur. On utilise deux schémas Volumes Finis, le schéma de Rusanov et le schéma de Godunov approché VFRoe-ncv pour l'étape de convection. Plusieurs cas tests sont présentés et commentés : tubes à choc, conditions limites de paroi, robinet d'eau, sédimentation

CE TRAVAIL EST CONSACRE A LA SIMULATION NUMERIQUE DIRECTE DES ECOULEMENTS DIPHASIQUES COMPRESSIBLES. L'APPROCHE RETENUE EST L'UTILISATION DES METHODES DE CAPTURE D'INTERFACE. CES TECHNIQUES FONT APPEL A UN MAILLAGE EULERIEN FIXE TANT POUR DECRIRE LES VARIABLES DE L'ECOULEMENT QUE POUR L'INTERFACE ENTRE LES FLUIDES. DANS UNE PREMIERE PARTIE NOUS RAPPELONS LES METHODES CLASSIQUES DE CAPTURE D'INTERFACE EN DISTINGUANT CELLES BASEES SUR DES FONCTIONS COULEURS DISCONTINUES DE CELLES BASEES SUR LES FONCTIONS LEVEL SET. CETTE ETUDE MET EN EVIDENCE CERTAINES DIFFICULTES : PAR EXEMPLE CONCILIER LES SCHEMAS NUMERIQUES CONSERVATIFS AVEC UNE DESCRIPTION CONVENABLE DE LA PRESSION. EN PARTICULIER, NOUS DETAILLONS UN MODELE QUASI-CONSERVATIF POUR DEUX STIFFENED GAS ET LE SOLVEUR NUMERIQUE BASE SUR LE SCHEMA DE ROE QUI LUI EST ASSOCIE. NOUS PROPOSONS DANS UNE DEUXIEME PARTIE UNE EXTENSION DE L'APPROCHE PRECEDENTE : LE MODELE A CINQ EQUATIONS, AFIN DE PERMETTRE LE PLUS LARGE CHOIX POSSIBLE D'EQUATIONS D'ETAT POUR LES FLUIDES EN PRESENCE. NOUS PRESENTONS TROIS VARIANTES DE CE SYSTEME ET CONSTRUISONS POUR L'UNE D'ELLES UN SOLVEUR INSPIRE DU SCHEMA DE ROE. LA TROISIEME PARTIE DE NOTRE TRAVAIL EST CONSACREE A UN RAFFINEMENT DES METHODES DE CAPTURE D'INTERFACE ET PLUS SPECIFIQUEMENT AUX PROBLEMES DE DIFFUSION NUMERIQUE DE L'INTERFACE. UNE TECHNIQUE BASEE SUR LE RAIDISSEMENT DE LA FONCTION COULEUR ET UNE AUTRE QUI ALLIE LEVEL SET ET METHODE QUASI-CONSERVATIVE SONT PROPOSEES. ENFIN, UNE DERNIERE PARTIE PORTE SUR L'ETUDE DU CHANGEMENT DE PHASE DYNAMIQUE. NOUS ETUDIONS GRACE AUX OUTILS DE LA THERMODYNAMIQUE NON-CLASSIQUE LES STRUCTURES D'UNE INTERFACE QUI REALISE LA TRANSITION DE PHASE. IL EN RESSORT QU'UNE RELATION DE SAUT A L'INTERFACE ENTRE LES FLUIDES OU RELATION CINETIQUE FOURNIT LE FLUX DE MASSE A TRAVERS L'INTERFACE DU AU CHANGEMENT DE PHASE. NOUS PRESENTONS DE MEME UNE IMPLEMENTATION NUMERIQUE SIMPLE DANS LE CADRE DU MODELE A CINQ EQUATIONS.