Cette thèse est consacrée à la construction et à l’étude mathématique et numérique d’un modèle d’écoulement diphasique à une phase incompressible. La première partie présente l’établissement du modèle. Le point de départ en est le modèle à deux fluides bien connu dans la littérature spécialisée que l’on considère ici sous sa forme isotherme et isobare et qui se traduit (en une dimension d’espace) par un système de quatre équations couplées. En utilisant la technique du développement de Chapman-Enskog dans la limite d’un temps de relaxation des vitesses tendant vers 0, on montre que ce système peut se réduire à un système à trois équations de conservation et on obtient une loi de comportement de type Darcy pour le déséquilibre des vitesses. La deuxième partie de ce travail est consacrée à l’analyse mathématique de ce modèle. On montre qu’il est hyperbolique, et on donne la solution exacte du problème de Riemann. Enfin, dans la dernière partie, on s’intéresse à l’approximation numérique de ce système. On développe des méthodes numériques basées sur des solveurs de Riemann exact et approchés pour l’approximation des termes hyperboliques et sur des méthodes d’éléments finis pour l’approximation des termes de déséquilibre des vitesses. On construit ensuite des méthodes implicites en temps pour ce type de discrétisation et on poursuit par la mise au point de schémas implicites à deux pas. On conclut par quelques applications 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

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

Dans ce mémoire de thèse, on s'intéresse à la simulation des écoulements liquidevapeur en transition de phase. Pour décrire ces écoulements, une approche bifluide moyennée à deux pressions indépendantes est retenue. Cette description du mélange liquide-vapeur s'appuie sur le modèle à sept équations de Baer et Nunziato. On étudie les aptitudes de cette modélisation à simuler les transitions de phase apparaissant en ingénierie nucléaire. Dans un premier temps, on élabore un cadre thermodynamique théorique pour décrire les écoulements liquide-vapeur. Dans ce cadre, on réalise la fermeture du modèle de Baer et Nunziato. De nouvelles modélisations sont proposées pour les termes d'interaction entre les phases. Ces nouvelles modélisations dotent le modèle bifluide à deux pressions d'une inégalité d'entropie. On étudie ensuite les propriétés mathématiques de ce modèle. Sa partie convective hyperbolique se présente sous une forme non-conservative. On étudie tout d'abord la définition de ses solutions faibles. Divers régimes d'écoulement sont alors mis à jour pour le mélange diphasique. Ces différents régimes d'écoulement présentent des analogies avec le comportement fluvial et torrentiel des écoulements en rivière. Les stabilités linéaire et non-linéaire de l'équilibre liquidevapeur sont ensuite établies. Pour affiner notre description des interactions diphasiques, on étudie pour finir l'implémentation d'un modèle de turbulence, ainsi que l'implémentation d'une procédure de reconstruction pour la densité d'aire interfaciale. On s'intéresse ensuite à la simulation de ce modèle. Suivant une approche à pas fractionnaires, une méthode numérique est élaborée dans un formalisme Volumes Finis. Pour réaliser l'approximation de la partie convective, diverses adaptations non-conservatives de solveurs de Riemann standard sont tout d'abord proposées. A l'inverse du cadre non-conservatif classique, l'ensemble de ces schémas converge vers une unique solution. Un nouveau schéma de relaxation est ensuite proposé pour approcher la dynamique des transferts interfaciaux. L'ensemble de la méthode numérique se caractérise alors par la préservation des équilibres liquide-vapeur. Dans un premier temps, cette méthode numérique est employée à la comparaison des différentes modélisations bifluides à une et deux pressions. On l'applique ensuite à la simulation des écoulements liquide-vapeur dans les circuits hydrauliques des réacteurs à eau sous pression en configuration accidentelle.

Ce travail de thèse est consacré à la modélisation et simulation numérique d'écoulements diphasiques liquide-gaz mettant en jeu des transferts de chaleur. La simulation de configurations où la prise en compte des effets de compressibilité de la phase gazeuse est indispensable (micropompes, microactionneurs, etc...) a nécessité l'utilisation d'un modèle original, considérant le liquide incompressible et le gaz compressible sous l'hypothèse faible Mach. Lors de cette thèse, ce modèle a été implémenté dans un code diphasique prenant en compte l'interface à l'aide d'une méthode de front-tracking. Des cas tests ont été développés spécifiquement afin de vérifier la conservation de l'énergie pour des configurations de complexité croissante. Les résultats des cas tests ont permis de mettre en évidence la difficulté à assurer la conservation de l'énergie lorsque l'interface n'est pas discontinue mais lissée, comme c'est le cas dans la méthode de front-tracking standard. Une méthode de traitement d'interface hybride a été proposée, rétablissant le caractère discontinu de l'interface avec la reconstruction d'une fonction indicatrice de phase échelon, tandis que le déplacement de l'interface est assuré d'un pas de temps à l'autre à l'aide du front-tracking. Les résultats obtenus avec cette nouvelle méthode hybride sont très satisfaisants, la méthode hybride permettant d'assurer la conservation de l'énergie et de la masse avec précision dans les simulations.

NOTRE OBJECTIF EST DE CONSTRUIRE UN OUTIL NUMERIQUE ADAPTE A L'ETUDE DU REMPLISSAGE D'UN RESERVOIR A CARBURANT. UNE TELLE APPLICATION FAIT INTERVENIR DEUX PROBLEMES COUPLES: L'EVOLUTION DE L'ETAT DU CARBURANT ET DU GAZ - CONSIDERES COMME DEUX FLUIDES INCOMPRESSIBLES ET IMMISCIBLES - ET LE DEPLACEMENT DES INTERFACES GAZ-LIQUIDE DANS LE RESERVOIR. NOUS CONSTRUISONS TOUT D'ABORD UN SCHEMA IMPLICITE ORIGINAL DEDIE A LA RESOLUTION DES EQUATIONS DE NAVIER-STOKES EN REGIME INCOMPRESSIBLE INSTATIONNAIRE. NOUS UTILISONS A CET EFFET UNE METHODE DE COMPRESSIBILITE ARTIFICIELLE COUPLEE A UNE TECHNIQUE DE PAS DE TEMPS DUAL. L'ANALYSE DES PROPRIETES LINEAIRES DE STABILITE ET DE CONVERGENCE DE CE SCHEMA PERMET D'EVALUER SON COMPORTEMENT NUMERIQUE AINSI QUE LES PLAGES DE VARIATION DE SES PARAMETRES DE CONTROLE. PLUSIEURS CAS DE VALIDATION - PLAQUE PLANE ET CAVITES ENTRAINEES - SONT CONSIDERES POUR DIVERS NOMBRES DE REYNOLDS EN REGIMES STATIONNAIRE ET INSTATIONNAIRE. LA METHODE DE SUIVI D'INTERFACE PAR ISO-SURFACES EST ENSUITE PRESENTEE. NOUS DETAILLONS SES CARACTERISTIQUES PUIS PRECISONS CERTAINS APPORTS SPECIFIQUES TELS QUE L'ECRITURE DES CONDITIONS AUX LIMITES, L'UTILISATION D'UN PROCEDE DE REINITIALISATION ET LE TRAITEMENT DE RAPPORTS DE MASSE VOLUMIQUE ELEVES. LA TENSION SUPERFICIELLE EST MODELISEE A L'AIDE DE LA METHODE DES FORCES DE SURFACE CONTINUES. L'ENSEMBLE EST VALIDE PAR LE CALCUL DES INSTABILITES HYDRODYNAMIQUES DE RAYLEIGH-TAYLOR ET DE KELVIN-HELMHOLTZ, ET PAR LA SIMULATION DE LA PROPAGATION D'ONDES CAPILLAIRES. NOUS MONTRONS EGALEMENT QUE LA MASSE DE CHAQUE FLUIDE EST CONSERVEE A L'ORDRE DE PRECISION DU SCHEMA NUMERIQUE. DES APPLICATIONS DANS LE DOMAINE DE L'AERODYNAMIQUE INTERNE D'UN VEHICULE AUTOMOBILE SONT FINALEMENT DECRITES. LA MODELISATION DE L'IMPACT NORMAL DE GOUTTES SUR UNE SURFACE SOLIDE MET EN EVIDENCE LE ROLE ESSENTIEL JOUE PAR LA TENSION SUPERFICIELLE ET LA COURBURE DE LA PAROI. DES MODIFICATIONS APPORTEES AU LOGICIEL INDUSTRIEL DE MECANIQUE DES FLUIDES FIRE (AVL) AUTORISENT LA SIMULATION NUMERIQUE DU REMPLISSAGE DE RESERVOIRS 2D ET 3D