Méthode numérique lagrangienne pour la simulation d'écoulements à surface libre

Méthode numérique lagrangienne pour la simulation d'écoulements à surface libre

Author: Jean-Christophe Marongiu

Publisher:

Published: 2007

Total Pages: 170

ISBN-13:

DOWNLOAD EBOOK

La méthode SPH (Smoothed Particle Hydrodynamics) est une méthode numérique sans maillage utilisée dans cette étude pour la discrétisation spatiale des équations de la mécanique des fluides (essentiellement les équations d 'Euler). La méthode SPH rencontre depuis quelques années un certain succès dans la simulation d’écoulements à surface libre car son formalisme lagrangien facilite le traitement et le suivi des interfaces. Cette étude a pour but d’appliquer cette méthode pour la simulation des écoulements â surface libre se produisant dans les turbines Pelton. Le formalisme SPH standard est tout d’abord testé, il permet de valider la faisabilité de ce choix mais montre également les limites de la méthode SPH standard, en terme de précision et de fiabilité notamment. Le choix s'est alors porté vers une formulation hybride SPH-ALE (Arbitrary Lagrange Euler) qui entretient une certaine filiation avec le formalisme des volumes finis. SPH-ALE utilise en effet une formulation conservative des équations du mouvement et est capable théoriquement de décrire l'écoulement quelque soit le déplacement des points de discrétisation. Par ailleurs, sur un plan purement numérique, ce formalisme permet l'utilisation de schémas numériques décentrés, en particulier les schémas de type Godunov et leurs variantes d’ordre supérieur. Cette méthode hybride se révèle en pratique nettement supérieure A la méthode standard pour les applications visées. La stabilité des simulations est largement renforcée, et la précision des résultats est fortement améliorée. En particulier le champ de pression retrouve une forme satisfaisante sans lissage numérique particulier. La méthode hybride facilite également le traitement des conditions limites. Alors que ce point constitue une difficulté majeure pour la méthode SPH standard, la méthode SPH-ALE permet de traiter les conditions limites à travers des flux aux frontières qui peuvent être eux-aussi décentrés. La mise en place d’un traitement cohérent et rigoureux des conditions limites constitue la principale contribution de ce travail de thèse. La méthode SPH-ALE est finalement testée sur des cas représentatifs des applications visées et fournit des résultats satisfaisants. En particulier le champ de pression en paroi solide est prédit correctement. En conclusion, les développements effectués dans cette étude ont été guidés par l'application en turbine Pelton qui était visée. La nécessité de manipuler des géométries complexes et d'obtenir un niveau de précision correct ont conduit à privilégier et à développer la méthode hybride SPH-ALE. Ce travail ouvre des perspectives prometteuses de développement rapide grâce au lien existant entre SPH-ALE et la méthode des volumes finis.


Simulation numérique des écoulements à surface libre avec turbulence

Simulation numérique des écoulements à surface libre avec turbulence

Author: FRANCOIS.. BOUTTES

Publisher:

Published: 1994

Total Pages: 341

ISBN-13:

DOWNLOAD EBOOK

La méthode des éléments finis (m.e.f.) A été appliquée avec succès pour représenter l'hydrodynamique a surface libre, mais avec une modélisation de la turbulence très simplifiée. Ce travail a pour but de simuler, a l'aide de la m.e.f., des écoulements a surface libre avec des représentations de la turbulence plus complexes. Il commence par une présentation générale des modèles de turbulence. Nous avons choisi d'utiliser plus particulièrement deux types de modèle: le modèle k-, des modèles aux longueurs de mélange développés a partir de la simplification du modèle précédent. Une simplification supplémentaire est effectuée en bidimensionnalisant les modèles. Un travail de ce type a déjà été effectué par rodi et al pour le modèle k-. Nous l'avons complété pour les modèles aux longueurs de mélange. Dans le cas ou les courbures de l'écoulement horizontal sont importantes, nous proposons des variantes des deux types de modèle précédents et par ailleurs, une modification des frottements de fond intervenant dans l'équation de l'écoulement (saint-venant). Cette dernière est notamment importante pour les études sedimentologiques. La méthode de résolution des équations du modèle k- utilisée dans cette thèse est la méthode des éléments finis ; il s'agit alors de résoudre un système non linéaire complexe ; dans ce travail, nous avons fait un choix d'interpolation particulier et teste plusieurs méthodes de résolution. Le couplage avec le modèle hydrodynamique est effectue par un algorithme explicite. Sont présentés plusieurs critères pour contrôler la convergence du modèle k-. Nous terminons par une validation des modèles de turbulence utilises sur de nombreux exemples. Sur ces derniers, nous avons vérifié les méthodes employées pour résoudre les équations du modèle k-


Analyse mathématique et simulation numérique de modèles d'écoulement de fluides incompressibles en surface libre et milieu poreux déformable

Analyse mathématique et simulation numérique de modèles d'écoulement de fluides incompressibles en surface libre et milieu poreux déformable

Author:

Publisher:

Published:

Total Pages:

ISBN-13:

DOWNLOAD EBOOK

Milieu poreux déformable. On dérive un modèle mathématique pour un milieu poreux déformable en utilisant la conservation de la masse, la loi de Darcy généralisée et le tenseur de déformation dépendant de la fonction de retrait de vertisol de E. Braudeau. On démontre l'existence globale d'une solution faible de ce problème. On utilise pour la preuve un théorème de compacité de Dubinskii, généralisé pour des hypothèses plus faibles. On montre que les paramètres de Van Genuchten associés à la conductivité hydraulique et au potentiel matriciel sont liés à la régularité de la solution faible. Un logiciel calculant la teneur en eau du milieu poreux déformable, écrit en fortran 90, est mis en place. La méthode des éléments finis et un schéma d'Euler implicite sont utilisés pour le code de calcul. Modèle de Saint-Venant visqueux. On étudie l'existence de solutions pour un modèle de Saint-Venant non linéaire décrivant mathématiquement l'écoulement d'un cours d'eau pour une formulation hauteur-débit. Utilisant la technique du point fixe de Banach dans le cas du système linéarisé on démontre l'existence globale et l'unicité de solutions. Dans le cas non linéaire, l'utilisation du théorème de point fixe de Schauder et d'un résultat de compacité de J. Simon nous permet de prouver l'existence d'une solution pour des temps petits. Des résultats d'une simulation numérique du système de Saint-Venant sont présentés dans le cas d'une rupture de barrage. Des comparaisons sont données en fonction de la viscosité et de la non-linéarité du problème.


Analyse a posteriori et adaptation de maillage pour des problèmes d'écoulements souterrains et à surface libre

Analyse a posteriori et adaptation de maillage pour des problèmes d'écoulements souterrains et à surface libre

Author: Tarek Ghoudi

Publisher:

Published: 2018

Total Pages: 0

ISBN-13:

DOWNLOAD EBOOK

L'objectif de cette thèse est d'analyser et de développer des outils numériques adaptatifs efficaces pour les problèmes d'écoulements souterrains et à surface libre, en proposant une nouvelle méthode d'adaptation basée sur les estimateurs d'erreur a posteriori.Dans la première partie, nous appréhendons le cadre mathématique en présentant une analyse détaillée sur l'existence et l'unicité ainsi que la convergence d'erreur dans les différentes normes pour les équations d'écoulement et de transport en milieux poreux.La deuxième partie est dédiée à une nouvelle stratégie d'adaptation de maillage qui consiste à coupler deux stratégies d'adaptation, à savoir la méthode Adapt et la stratégie Newest Vertex bisection (NVB). La première consiste à diviser un triangle en quatre-sous triangles, et la deuxième, consiste à le diviser en joignant le somment opposé à l'arête la plus longue par son milieu. La conformité de notre méthode émane de cette dernière méthode.La troisième partie est consacrée à la démarche empruntée pour valider notre nouvelle stratégie.Celle-ci s'appuie sur le schéma volumes finis "vertex centered", en considérant l'équation elliptique de second ordre à coefficients discontinus afin de prendre en compte les hétérogénéités du milieu. Des tests numériques viendront corroborer l'efficacité de notre méthode, ainsi que la convergence de l'erreur exacte et de l'estimateur, dont le rapport définit l'indice d'efficacité qui est proche de 1.Dans la quatrième partie, nous avons essayé d'optimiser cette stratégie d'adaptation, en concevant une stratégie multi-niveaux pour aboutir enfin à une nouvelle méthode numérique volumes finis-semi lagrangienne. Son principe est de résoudre le problème par une phase volume finis(phase correcteur) précédée d'une phase lagrangienne utilisant la méthode des caractéristiques(phase prédicteur). La phase volumes finis utilise comme flux numérique aux interfaces le vraiflux physique évalué en un état approché aux interfaces, obtenu dans la phase prédicteur par la méthode des caractéristiques. Elle fait également appel à des processus d'interpolation devant être judicieusement choisis.


Simulation d'écoulements à surface libre entre un liquide et un gaz compressible

Simulation d'écoulements à surface libre entre un liquide et un gaz compressible

Author: Sébastien Lortie

Publisher:

Published: 2018

Total Pages: 89

ISBN-13:

DOWNLOAD EBOOK

Ce mémoire présente une formulation éléments finis qui permet de résoudre des écoulements compressibles instationnaires subsoniques. La motivation de ce travail provient de la volonté d'améliorer les prédictions d'écoulements diphasiques séparés par une interface. En particulier, on s'intéresse au déferlement de vagues qui emprisonnent une poche d'air. D'abord, on développe une formulation des équations de Navier-Stokes compressible qui permet de gérer plusieurs régimes d'écoulements compressibles dans le même domaine. Pour résoudre ce système d'équations, on utilise une formulation par éléments finis. Pour les problèmes instationnaires, on utilise la méthode des " Backward Differentiation Formula " (BDF) avec un estimateur d'erreur de troncature pour effectuer l'intégration temporelle. En plus, on utilise une formulation ALE qui laisse le maillage de se déplacer dans le temps pour lui permettre de suivre les déformations d'une ou plusieurs frontières du domaine (comme c'est le cas pour un écoulement avec interface). Ces outils mentionnés ci-haut sont déjà tous utilisés pour le calcul d'écoulements à interface, mais pour des fluides incompressibles seulement. Donc, le but de ce mémoire est d'implémenter dans les formulations les termes compressibles qui s'ajoutent aux équations incompressibles. Pour vérifier que l'implémentation est bien faite, on effectue la vérification par la méthode des solutions manufacturées. On présente d'abord un problème stationnaire compressible pour montrer que les termes stationnaires sont bien résolus suivi par une solution manufacturée instationnaire compressible qui fait intervenir tous les termes ajoutés. Ainsi, on peut vérifier que les méthodes BDF se comportent de la bonne façon pour l'estimation des dérivées temporelles de la pression et de la température (essentiellement les termes ajoutés par la formulation compressible). Ensuite, puisque le problème de la vague déferlante se caractérise par la compression d'une poche de gaz, on procède à la validation du code avec le problème de Bagnold. Ce problème en 1D consiste à évaluer la pression de l'air dans un contenant compressé par un piston. Une fois la validation terminée, on se lance dans le calcul final du cas de déferlement de la vague. Les résultats obtenus sont encourageants comparativement à la formulation tout incompressible, mais il reste du travail à faire pour obtenir un modèle qui représente le mieux possible la réalité


Modélisation et contrôle d’écoulements à surface libre par la méthode de Boltzmann sur réseau

Modélisation et contrôle d’écoulements à surface libre par la méthode de Boltzmann sur réseau

Author: Olivier Marcou

Publisher:

Published: 2010

Total Pages: 242

ISBN-13:

DOWNLOAD EBOOK

Ce travail de thèse traite de la modélisation et la simulation des systèmes complexes et s’inscrit dans la thématique du contrôle et la gestion des ressources en eau. Nous proposons une approche originale basée sur les modèles de Boltzmann sur réseau pour représenter les écoulements au sein des canaux d’irrigation, habituellement décrits par les équations non linéaires de Saint-Venant. Nous avons adapté un modèle bi-fluide et étudié les conditions aux bords qui permettent de reproduire la géométrie d’un canal à surface libre. Des méthodes de détermination des grandeurs hydrauliques d’intérêt ont été développées. Nous nous sommes intéressés au comportement des vannes de fond en régime noyé, et nous démontrons que le modèle permet de décrire celui-ci spontanément dans différentes situations de fonctionnement. Des validations ont été effectuées par le biais d’une comparaison avec des expérimentations menées sur un micro-canal expérimental. Nous avons également introduit dans le modèle des phénomènes de sédimentation et étudié l’influence de la présence d’un dépôt de sédiments sur l’écoulement. Là aussi, les résultats de simulation et d’expérimentation convergent.


COUPLAGE DE METHODES NUMERIQUES EN SIMULATION DIRECTE D'ECOULEMENT INCOMPRESSIBLES

COUPLAGE DE METHODES NUMERIQUES EN SIMULATION DIRECTE D'ECOULEMENT INCOMPRESSIBLES

Author: MOHAMED LEMINE.. OULD SALIHI

Publisher:

Published: 1998

Total Pages: 240

ISBN-13:

DOWNLOAD EBOOK

CE TRAVAIL EST CONSACRE AU DEVELOPPEMENT DES METHODES LAGRANGIENNES COMME ALTERNATIVES OU COMPLEMENTS AUX METHODES EULERIENNES CONVENTIONNELLES POUR LA SIMULATION D'ECOULEMENTS INCOMPRESSIBLES EN PRESENCE D'OBSTACLES. ON S'INTERESSE EN PARTICULIER A DES TECHNIQUES OU DES SOLVEURS EULERIENS ET LAGRANGIENS COHABITENT DANS LE MEME DOMAINE DE CALCUL MAIS TRAITENT DIFFERENTS TERMES DIFFERENTS DES EQUATIONS DE NAVIER-STOKES, AINSI QU'A DES TECHNIQUES DE DECOMPOSITION DE DOMAINES OU DIFFERENTS SOLVEURS SONT UTILISES DANS CHAQUES SOUS-DOMAINES. LORSQUE LES METHODES EULERIENNES ET LAGRANGIENNES COHABITENT DANS LE MEME DOMAINE DE CALCUL (METHODE V.I.C.), LES FORMULES DE PASSAGE PARTICULES-GRILLES PERMETTENT DE REPRESENTER LA VORTICITE AVEC LA MEME PRECISION SUR UNE GRILLE FIXE ET SUR LA GRILLE LAGRANGIENNE. LES METHODES V.I.C. AINSI OBTENUES COMBINENT STABILITE ET PRECISION ET FOURNISSENT UNE ALTERNATIVE AVANTAGEUSE AUX METHODES DIFFERENCES-FINIES POUR DES ECOULEMENTS CONFINES. LORSQUE LE DOMAINE DE CALCUL EST DECOMPOSE EN SOUS-DOMAINES DISTINCTS TRAITES PAR METHODES LAGRANGIENNES ET PAR METHODES EULERIENNES, L'INTERPOLATION D'ORDRE ELEVE PERMET DE REALISER DES CONDITIONS D'INTERFACE CONSISTANTES ENTRE LES DIFFERENTS SOUS-DOMAINES. ON DISPOSE ALORS DE METHODES DE CALCUL AVEC DECOMPOSITION EN SOUS-DOMAINES, DE TYPE EULER/LAGRANGE OU LAGRANGE/LAGRANGE, ET RESOLUTION EN FORMULATION (VITESSE-TOURBILLON)/(VITESSE-TOURBILLON) OU (VITESSE-PRESSION)/(VITESSE-TOURBILLON). LES DIFFERENTES METHODES DEVELOPPEES ICI SONT TESTEES SUR PLUSIEURS TYPES D'ECOULEMENTS (CAVITE ENTRAINEE, REBOND DE DIPOLES DE VORTICITE, ECOULEMENT DANS UNE CONDUITE ET SUR UNE MARCHE, ECOULEMENT AUTOUR D'OBSTACLES) ET COMPAREES A DES METHODES DE DIFFERENCES-FINIES D'ORDRE ELEVE.


SIMULATION NUMERIQUE D'ECOULEMENTS PAR DES METHODES DE TRANSFORMATION DE DOMAINES

SIMULATION NUMERIQUE D'ECOULEMENTS PAR DES METHODES DE TRANSFORMATION DE DOMAINES

Author: AHMAD.. MAHMOUD

Publisher:

Published: 1999

Total Pages: 216

ISBN-13:

DOWNLOAD EBOOK

CE TRAVAIL FAIT LE POINT SUR LES METHODES DE TRANSFORMATION DE DOMAINES UTILISEES POUR LA SIMULATION NUMERIQUE D'ECOULEMENTS BI- ET TRIDIMENSIONNELS. DANS CE CADRE, NOUS AVONS DEVELOPPE DES METHODES DE TRANSFORMATION PERMETTANT DE CALCULER DES ECOULEMENTS ETABLIS DANS DES FILIERES A SECTION CONSTANTE, MAIS QUELCONQUE, ET DES ECOULEMENTS COMPLEXES, POUR DES MODELES RHEOLOGIQUES DE TYPE NEWTONIEN, NON-NEWTONIEN PUREMENT VISQUEUX ET VISCOELASTIQUE. A LA DIFFERENCE DES ANALYSES PROPOSEES DANS LA LITTERATURE, NOUS CONSIDERONS UNE TRANSFORMATION INCONNUE ENTRE LE DOMAINE PHYSIQUE ET UN DOMAINE TRANSFORME (DE CALCUL) PLUS SIMPLE, ANALYSE A LAQUELLE SE RATTACHENT LES METHODES DE TUBES DE COURANT. PLUSIEURS PROBLEMES ONT ETE ABORDES ET RESOLUS : 1) LA SIMULATION D'ECOULEMENTS ETABLIS DE FLUIDES OBEISSANT A DES LOIS DE COMPORTEMENT NEWTONIENNE, NON-NEWTONIENNE PUREMENT VISQUEUSE ET VISCOELASTIQUE INTEGRALE POUR DES CONDUITES DE SECTIONS COMPLEXES, POUR LESQUELLES UNE APPROCHE PAR DECOMPOSITION DE DOMAINES A ETE MISE EN OEUVRE. 2) LE CALCUL D'ECOULEMENTS AXISYMETRIQUES A SURFACE LIBRE INCONNUE CORRESPONDANT A L'EXTRUSION EN SORTIE DE FILIERE ANNULAIRE, POUR UN LIQUIDE NEWTONIEN. 3) LA SIMULATION DE PROBLEMES D'ECOULEMENTS AXISYMETRIQUES COMPORTANT DES ZONES DE SINGULARITE : SURFACES AVEC TRANSITION ADHERENCE-GLISSEMENT OU LIBRES. DANS LE CONTEXTE D'UNE ANALYSE DE TUBES DE COURANT, UNE FORMULATION ORIGINALE PAR ELEMENTS FINIS UTILISANT LA METHODE DE GALERKIN A ETE PROPOSEE POUR LES CAS 2D ET 3D. LES RESULTATS, OBTENUS AVEC DES ELEMENTS D'APPROXIMATION D'HERMITE SONT CONFORMES A CEUX DE LA LITTERATURE, POUR LE CAS NEWTONIEN ET EN ACCORD AVEC LES DONNEES EXPERIMENTALES SUR UN POLYETHYLENE LINEAIRE DE BASSE DENSITE REPRESENTE PAR UNE LOI DE COMPORTEMENT INTEGRALE.