Simulation Numérique D'écoulements Multiphasiques, Problèmes À Interfaces Et Changement de Phase PDF Download
Are you looking for read ebook online? Search for your book and save it on your Kindle device, PC, phones or tablets. Download Simulation Numérique D'écoulements Multiphasiques, Problèmes À Interfaces Et Changement de Phase PDF full book. Access full book title Simulation Numérique D'écoulements Multiphasiques, Problèmes À Interfaces Et Changement de Phase by Damien Furfaro. Download full books in PDF and EPUB format.
Author: Damien Furfaro
Publisher:
ISBN:
Category :
Languages : en
Pages : 113
Get Book Here
Book Description
This work deals with the numerical simulation of compressible multiphase flows in velocity disequilibrium. A HLLC-type two-phase Riemann solver is developed and validated against exact solutions and experimental data. This solver is robust, simple, accurate and entropy preserving. The numerical method is then implemented in 3D unstructured meshes. Furthermore, a numerical technique consisting in enforcing the correct energy partition at a discrete level in agreement with the multiphase shock relations is built. The multiphase extension of the HLLC-type Riemann solver is realized and allows the simulation of a wide range of applications. Finally, a droplet heat and mass transfer model with large range of validity is derived. It is valid in any situation: evaporation, flashing and condensation. It accounts for coupled heat and mass diffusion in the gas phase, thermodynamics of the multi-component gas mixture and heat diffusion inside the liquid droplet, enabling in this way consideration of both droplets heating and cooling phenomena.
Author: Damien Furfaro
Publisher:
ISBN:
Category :
Languages : en
Pages : 113
Get Book Here
Book Description
This work deals with the numerical simulation of compressible multiphase flows in velocity disequilibrium. A HLLC-type two-phase Riemann solver is developed and validated against exact solutions and experimental data. This solver is robust, simple, accurate and entropy preserving. The numerical method is then implemented in 3D unstructured meshes. Furthermore, a numerical technique consisting in enforcing the correct energy partition at a discrete level in agreement with the multiphase shock relations is built. The multiphase extension of the HLLC-type Riemann solver is realized and allows the simulation of a wide range of applications. Finally, a droplet heat and mass transfer model with large range of validity is derived. It is valid in any situation: evaporation, flashing and condensation. It accounts for coupled heat and mass diffusion in the gas phase, thermodynamics of the multi-component gas mixture and heat diffusion inside the liquid droplet, enabling in this way consideration of both droplets heating and cooling phenomena.
Author: Sami Ammar
Publisher:
ISBN:
Category :
Languages : fr
Pages : 164
Get Book Here
Book Description
Les écoulements multiphasiques avec changement de phase et transfert de chaleur interviennent dans de nombreuses applications industrielles et processus naturels. Malgré le développement rapide des outils numériques, la simulation numérique d'écoulements fluides à plusieurs phases immiscibles demeure un défi en raison de la difficulté inhérente à suivre les interfaces. Habituellement, pour résoudre ce type d'écoulement, des méthodes numériques basées sur les équations de Navier Stokes couplées à des méthodes lagrangiennes ou eulériennes sont utilisées. Cependant, ces dernières décennies, la méthode Boltzamnn sur réseau est apparue comme une méthode prometteuse pour simuler les écoulements à géométrie complexe et à plusieurs phases. Dans ce contexte, l'objectif principal de cette thèse est de développer un modèle multiphasique 3D basé sur la méthode de Boltzmann sur réseau capable d'étudier l'impact d'une gouttelette sur une surface et de traiter les changements de phase. Pour ce faire, le modèle pseudo-potentiel de Shan & Chen pour la simulation d'écoulement à plusieurs phases et plusieurs composants est utilisé comme point de départ. Cette recherche doctorale est organisée autour de trois thèmes afin de mettre en évidence les contributions. Des cas tests sont réalisés afin de valider les améliorations apportées au modèle et des applications physiques sont proposées. Le premier thème apporte une amélioration au modèle pseudo-potentiel en découplant la tension de surface et la densité. Le tenseur de pression du modèle 3D est modifié pour permettre un ajustement de la tension de surface tout en préservant la consistance thermodynamique. Dans ce thème, l'étude de l'impact d'une gouttelette sur une surface sèche et mouillée est analysée. Le thème 2 généralise la méthode de modification du tenseur de pression aux cas d'écoulements à plusieurs composants, dont la viscosité des fluides est différente. Il permet de simuler des écoulements à ratios de densité et de viscosité élevés. Enfin, le thème 3 ajoute les effets thermiques au modèle multiphasique à plusieurs composants. Le changement de phase et les échanges de chaleur entre composants sont étudiés. Pour chacun des thèmes, un article scientifique a été rédigé et soumis à un journal.
Author: Didier Jamet
Publisher:
ISBN:
Category :
Languages : fr
Pages : 234
Get Book Here
Book Description
L'une des principales difficultés rencontrées en simulation numérique directe des écoulements diphasiques en général et des écoulements liquide-vapeur avec changement de phase en particulier, est le suivi des interfaces. L'idée développée dans ce travail consiste à ne pas modéliser une interface liquide-vapeur comme une surface de discontinuité mais comme une zone volumique a travers laquelle les grandeurs physiques varient continument. La théorie du second gradient permet d'obtenir des équations d'évolution du fluide dans l'ensemble d'un système : phases et interfaces. Cela signifie que la résolution d'un seul système d'équations aux dérivées partielles est nécessaire pour résoudre l'ensemble du problème diphasique, les interfaces et leur évolution faisant partie de la solution de ce seul système. On montre dans ce travail qu'il est possible d'épaissir artificiellement une interface sans changer sa tension interfaciale et la chaleur latente de changement d'état. Cela signifie qu'il est possible de suivre l'ensemble des interfaces d'un système diphasique liquide-vapeur avec changement de phase sur un maillage dont la taille est imposée par la plus petite échelle de Kolmogorov des phases par exemple. L'épaississement artificiel d'une zone interfaciale s'obtient en modifiant le comportement thermodynamique du fluide a l'intérieur de la binodale. On montre que cette modification n'a pas d'influence sur la dynamique d'une interface. En revanche, on montre que bien que l'épaisseur d'une interface et sa tension interfaciale varient en fonction du flux de masse et du flux thermique conductif qui la traversent, le changement de thermodynamique nécessaire à l'épaississement d'une interface amplifie de manière importante ces variations. Cela signifie en particulier qu'il faut être prudent dans la manière de choisir l'épaississement artificiel d'une interface si l'on ne souhaite pas trop faire varier sa tension interfaciale dans des problèmes dynamiques.
Author: SAMUEL.. KOKH
Publisher:
ISBN:
Category :
Languages : en
Pages : 327
Get Book Here
Book Description
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.
Author: Didier Jamet
Publisher:
ISBN:
Category :
Languages : fr
Pages : 234
Get Book Here
Book Description
L'une des principales difficultés rencontrées en simulation numérique directe des écoulements diphasiques en général et des écoulements liquide-vapeur avec changement de phase en particulier, est le suivi des interfaces. L'idée développée dans ce travail consiste à ne pas modéliser une interface liquide-vapeur comme une surface de discontinuité mais comme une zone volumique a travers laquelle les grandeurs physiques varient continument. La théorie du second gradient permet d'obtenir des équations d'évolution du fluide dans l'ensemble d'un système : phases et interfaces. Cela signifie que la résolution d'un seul système d'équations aux dérivées partielles est nécessaire pour résoudre l'ensemble du problème diphasique, les interfaces et leur évolution faisant partie de la solution de ce seul système. On montre dans ce travail qu'il est possible d'épaissir artificiellement une interface sans changer sa tension interfaciale et la chaleur latente de changement d'état. Cela signifie qu'il est possible de suivre l'ensemble des interfaces d'un système diphasique liquide-vapeur avec changement de phase sur un maillage dont la taille est imposée par la plus petite échelle de Kolmogorov des phases par exemple. L'épaississement artificiel d'une zone interfaciale s'obtient en modifiant le comportement thermodynamique du fluide a l'intérieur de la binodale. On montre que cette modification n'a pas d'influence sur la dynamique d'une interface. En revanche, on montre que bien que l'épaisseur d'une interface et sa tension interfaciale varient en fonction du flux de masse et du flux thermique conductif qui la traversent, le changement de thermodynamique nécessaire à l'épaississement d'une interface amplifie de manière importante ces variations. Cela signifie en particulier qu'il faut être prudent dans la manière de choisir l'épaississement artificiel d'une interface si l'on ne souhaite pas trop faire varier sa tension interfaciale dans des problèmes dynamiques.
Author: Olivier Basset
Publisher:
ISBN:
Category :
Languages : fr
Pages : 200
Get Book Here
Book Description
Cette thèse porte sur le développement de méthodes numériques pour le calcul d’écoulements incompressibles multi-fluides sur grille de calcul, s’inscrivant ainsi dans le cadre du projet MecaGrid.Une étude de la grille MecaGrid met en évidence son caractère hétérogène et ses conséquences. Plusieurs techniques d’optimisation sont présentées afin d’améliorer son utilisation : répartir la masse de calcul de façon adaptée, et privilégier le travail des processeurs vis-à-vis des communications réseaux pénalisantes.Des méthodes numériques sur maillage non structuré composé de tétraèdres sont choisies pour réaliser la simulation directe d’écoulements multi-fluides avec capture d’interface. Nous adoptons une approche unique qui rappelle celle du Multiscale, dans laquelle la condensation d’une fonction Bulle pyramidale est utilisée comme technique universelle de stabilisation.Les équations de Navier-Stokes incompressible sont résolues par une méthode éléments finis mixtes à interpolation P1+/P1. Un schéma temporel d’Euler implicite est appliqué, en association avec un algorithme de Newton pour linéariser le problème.Les interfaces sont capturées par une technique Level Set à interpolation continue P1 qui consiste à résoudre une équation de transport stabilisée par la condensation d’une Bulle (Residual-Free Bubbles). Le couplage avec une équation d’Hamilton-Jacobi permet de réinitialiser la fonction Level Set au court de son transport. La comparaison avec une méthode Galerkin discontinue proche du Volume of Fluid montre que le Level Set se distingue par sa simplicité et l’absence de diffusion numérique.Enfin, les simulations numériques sont validées par plusieurs cas test reconnus.
Author: Vincent Le maout
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
Les écoulements multiphasiques viscoélastiques en milieux poreux sont au centre de nombreuses recherches fondamentales et appliquées. Initiés par l'industrie pétrolière, leur développement s'étend aujourd'hui sur un large éventail de domaines tels que le génie civil, les géosciences, les matériaux composites et plus récemment les sciences de la vie. La formulation mathématique de ces problèmes repose généralement sur des équations gouvernantes formulées directement à l'échelle macroscopique, ou dérivées par changement d'échelle à partir des relations microscopiques. Généralement, la seconde approche est privilégiée, car elle permet d'établir une connexion claire entre les différentes échelles de description du milieu poreux et d'identifier les hypothèses res-trictives nécessaires à la formulation du problème. Cependant, lorsque le changement d'échelle est effectué, de nombreux paramètres restent encore à caractériser, et des conditions de fermeture supplémentaires sont à consi-dérer afin d'obtenir un système mathématique fermé. Pour les écoulements d'intérêt, présentant à la fois un caractère multiphasique et viscoélastique, la prise en compte des relations macroscopiques empiriques usuelles peut être trop imprécise pour rendre compte correctement de la physique sous-jacente, et peu de données expé-rimentales permettent de renseigner les paramètres manquants à l'échelle du milieu poreux.Une solution à ce problème consiste à modéliser numériquement les écoulements directement à l'échelle micros-copique afin de retrouver les grandeurs manquantes par des expériences de microfluidique in silico. Dans ces travaux de thèse, une formulation mathématique d'écoulement multiphasique viscoélastique a ainsi été dérivée à l'échelle du pore pour deux applications d'intérêt : la modélisation de la récupération améliorée du pétrole et la croissance tumorale. La dérivation d'un modèle unique pour ces deux applications utilise les lois de conservation de grandeurs extensives à l'échelle microscopiques considérées lors des opérations de changement d'échelle. La formulation résultante permet une description multiphasique à travers une méthode d'interface diffuse et la prise en compte de la viscoélasticité des phases est modélisée par une loi constitutive d'Oldroyd-B. Le modèle mathématique, implémenté dans un code éléments finis, permet d'étudier comment la rhéologie viscoélastique des solutions de polymère peut être utilisée afin d'améliorer la mobilisation du pétrole à l'échelle du pore lors d'opérations de récupération tertiaire. L'influence de la viscoélasticité sur les résultats numériques obtenus, tout comme l'efficacité du procédé, est comparée aux résultats expérimentaux existant de la littérature. Dans un second temps, le modèle mathématique a été spécialisé pour simuler la croissance d'agrégats de cellules cancéreuses de quelques centaines de microns. L'analogie entre les tissus biologiques et les fluides viscoélastiques est courante pour des fins de modélisations mathématiques depuis les travaux de Steinberg sur la dynamique des tissus vivants. Dans cette thèse, cette similarité est utilisée afin d'étudier le comportement d'agrégats tumoraux dans différents environnements. La formulation biomécanique du problème permet de simuler le comportement de populations cellulaires soumises à des chargements mécaniques externes, qui montrent une modulation de leur vitesse de croissance selon l'état de contrainte du milieu. L'utilisation du modèle mathématique dans ce contexte permet d'isoler les effets mécaniques des effets biologiques sur le régime de croissance, et propose des explications originales sur l'évolution d'agrégats tumoraux en milieux confinés. Enfin, les capacités prédictives de la formulation sur plusieurs expériences in vitro permettent d'illustrer la pertinence de l'utilisation de modèles d'écoulements multiphasiques viscoélastiques pour les problèmes de croissances tumorales.
Author: Gloria Faccanoni
Publisher:
ISBN:
Category :
Languages : fr
Pages : 422
Get Book Here
Book Description
Cette thèse étudie la modélisation et la simulation numérique d'écoulements diphasiques à interfaces avec changement de phase. Cette transition est localisée en des interfaces qui sont produites dynamiquement. On prend également en compte la diffusion de la chaleur, la tension de surface et les forces de gravité. L'application envisagée est la simulation d'écoulements dans un réacteur à eau pressurisée dans l'industrie nucléaire civile. On s'intéresse ici plus précisément à un éventuel fonctionnement accidentel et en particulier au phénomène de la crise d'ébullition. On modélise les écoulements diphasiques avec changement de phase par un modèle basé sur le système des équations d'Euler fermé par une seule équation d'état obtenue en postulant un équilibre instantané et local des pressions, températures et potentiels chimiques de chaque phase. On en étudie ensuite l'hyperbolicité et le problème de Riemann qui lui est associé. Du point de vue numérique, puisqu'il n'y a pas d'expression analytique pour la loi à l'équilibre dans le cas général, on propose une méthode simple pour approcher cette loi d'état lorsque les propriétés des deux phases sont décrites par des lois très générales, éventuellement sous forme tabulée. Enfin, pour simuler des écoulements diphasiques avec changement de phase, on présente un schéma numérique de type relaxation/projection pour lequel la phase de projection utilise cette approximation de l'équilibre thermodynamique.
Author: Hong Gen Liu
Publisher:
ISBN:
Category :
Languages : fr
Pages : 189
Get Book Here
Book Description
CE TRAVAIL EST CONSACRE A LA RESOLUTION NUMERIQUE DU PROBLEME DE CHANGEMENT DE PHASE SOLIDE-LIQUIDE, MONO ET BI-DIMENSIONNEL, EN REGIME DIFFUSIF. LA PREMIERE PARTIE CONCERNE LA MODELISATION DIRECTE DE CE PROBLEME: LA POSITION DU FRONT EST DEDUITE DE LA RESOLUTION DES EQUATIONS DE TRANSFERT DE CHALEUR DANS LES DEUX PHASES. ON ETUDIE DEUX MODELES: 1) UN MODELE A SUIVI DE FRONT AVEC FORMULATION EN TEMPERATURE, 2) UN MODELE ENTHALPIQUE. DANS TOUS LES CAS, LA DISCRETISATION EST REALISEE EN UTILISANT LA FORMULATION DU VOLUME DE CONTROLE. ON DETAILLE LES DIFFERENTES METHODES DE SIMULATION NUMERIQUE, ET ON PROPOSE DES ELEMENTS DE COMPARAISON PERMETTANT DE CHOISIR LA METHODE LA MIEUX ADAPTEE SUIVANT LES CAS ENVISAGES. LE MODELE A SUIVI DE FRONT FOURNIT UNE DESCRIPTION FINE DES TRANSFERTS AU VOISINAGE DE L'INTERFACE. AVEC LA FORMULATION ENTHALPIQUE, ON EST AMENE A RESOUDRE UN PROBLEME APPROCHE, ET LES PERFORMANCES DE LA METHODE SONT AMELIOREES PAR DES TECHNIQUES DE RAFFINEMENT DE MAILLAGE. UNE COMPARAISON ENTRE LES DEUX MODELES EST PRESENTEE EN 1D ET EN 2D. DANS LA DEUXIEME PARTIE, UNE APPROCHE INVERSE EST CONSIDEREE: EN UTILISANT LA MODELISATION DES TRANSFERTS DE CHALEUR DANS LA PHASE SOLIDE UNIQUEMENT, IL S'AGIT D'IDENTIFIER LA POSITION ET L'EVOLUTION TEMPORELLE DU FRONT, A PARTIR DE MESURES REALISEES DANS CETTE PHASE SOLIDE. CE PROBLEME INVERSE MAL-POSE EST RESOLU PAR UNE METHODE DE REGULARISATION. LE DOMAINE DE VALIDITE DE LA METHODE A ETE ANALYSE, EN FONCTION DU RAPPORT SIGNAL/BRUIT ET DE LA DISTANCE ENTRE LES CAPTEURS DE MESURE ET L'INTERFACE A IDENTIFIER
Author: Adlène Benkenida
Publisher:
ISBN:
Category :
Languages : en
Pages : 270
Get Book Here
Book Description
Ce travail est consacré au développement et à la mise en œuvre d'un outil de simulation directe d'écoulements diphasiques complexes où la topologie des interfaces évolue dans le temps. La stratégie de résolution emploie un maillage fixe sur lequel les interfaces évoluent librement. Les équations du modèle d'écoulement [...] sont obtenues par filtrage spatial et sommation des équations de Navier-Stokes locales et instantanées de chaque phase. L'utilisation d'une approche eulérienne conduit à des difficultés d'estimation de différentes grandeurs diphasiques [...]. Cet obstacle est levé en établissant et en validant une expression des contraintes valable quelle que soit la nature de l'approche eulérienne utilisée et l'orientation de l'interface par rapport au maillage. Pour simplifier les équations du modèle, on suppose l'absence de glissement local entre les deux phases, de changement de phase et de turbulence de petite échelle. La possibilité de relaxer certaines de ces hypothèses est discutée [...]. Le transport des interfaces est assuré par une méthode de capture de front sans reconstruction d'interfaces en utilisant un schéma FCT. On montre [...] que la version [...] de l'algorithme de Zalesak (1979) donne les meilleurs résultats, bien qu'il subsiste une tendance à l'épaississement des zones interfaciales dans les régions de fort étirement. La validation du code est réalisée [...]. Parmi les résultats de ces validations, ceux se rapportant au problème du rebond d'une bulle sur une paroi solide sont les plus originaux et apportent un éclairage nouveau sur ce phénomène, notamment en ce qui concerne l'évolution temporelle du champ hydrodynamique autour de la bulle. L'étude du mouvement d'une bulle de gaz isolée, dite bulle de Taylor, ou de l'interaction de deux bulles de Taylor montant dans un liquide en conduite verticale constitue le thème physique principal de ce travail. Plusieurs mécanismes physiques influençant la vitesse de montée et la dynamique du culot de ces bulles sont analysés (viscosité, tension interfaciale, vitesse superficielle du liquide et accélération de la bulle). Les comparaisons quantitatives avec des résultats théoriques et expérimentaux antérieurs sont très satisfaisantes et des résultats originaux sont obtenus sur plusieurs points.
