CE TRAVAIL DE THESE PORTE SUR LA SIMULATION NUMERIQUE DES ECOULEMENTS COMPRESSIBLES ET TURBULENTS AUTOUR D'ARRIERE-CORPS. POUR CELA, NOUS CONSIDERONS LES EQUATIONS DE NAVIER-STOKES MOYENNEES COMPLETEES PAR UN MODELE DE TURBULENCE. LES MODELES ALGEBRIQUES DE TYPE LONGUEUR DE MELANGE SONT LES PLUS SIMPLES, MAIS LES RESULTATS OBTENUS AVEC DE TELS MODELES SONT EN GENERAL PEU SATISFAISANTS POUR DES ECOULEMENTS COMPLEXES PRESENTANT DE GRANDES ZONES DE RECIRCULATION. DE PLUS, CES MODELES REQUIERENT UNE EVALUATION DE DISTANCES A LA PAROI, CE QUI PEUT PRESENTER UNE CERTAINE DIFFICULTE POUR DES CONFIGURATIONS D'ARRIERE-CORPS AVEC TUYERE ET CULOT. EN CONSEQUENCE, NOUS AVONS MIS EN UVRE UN MODELE DE TURBULENCE A DEUX EQUATIONS DE TRANSPORT PLUS ADAPTE A LA PHYSIQUE DE L'ECOULEMENT. NOUS AVONS TOUT D'ABORD CONSIDERE UN MODELE AVEC INTEGRATION JUSQU'A LA PAROI A SAVOIR LE MODELE K - DE JONES-LAUNDER. CE MODELE PRESENTE L'AVANTAGE QUE LES FONCTIONS INTRODUITES POUR TENIR COMPTE DES EFFETS DE PROCHE PAROI DEPENDENT UNIQUEMENT DU NOMBRE DE REYNOLDS TURBULENT ET NON DE DISTANCES A LA PAROI. AINSI, LE MODELE DE JONES-LAUNDER EST BEAUCOUP PLUS FACILE A APPLIQUER POUR DES CONFIGURATIONS D'ARRIERE-CORPS. NOUS NOUS SOMMES ENSUITE INTERESSE AUX LOIS DE PAROI. CELLES-CI PERMETTENT DE DECRIRE ANALYTIQUEMENT L'ECOULEMENT MOYEN ET LES GRANDEURS TURBULENTES AU VOISINAGE D'UNE PAROI ET AUTORISENT AINSI UNE TAILLE DE MAILLE A LA PAROI NETTEMENT PLUS IMPORTANTE QUE DANS UN CALCUL AVEC INTEGRATION JUSQU'A LA PAROI

On étudie dans le cadre de la thèse une famille de schémas numériques permettant de résoudre les équations de Saint-Venant. Ces schémas utilisent une décomposition d'opérateur de type Lagrange-projection afin de séparer les ondes de gravité et les ondes de transport. Un traitement implicite du système acoustique (relié aux ondes de gravité) permet aux schémas de rester stable avec de grands pas de temps. La correction des flux de pression rend possible l'obtention d'une solution approchée précise quel que soit le régime d'écoulement vis-à-vis du nombre de Froude. Une attention toute particulière est portée sur le traitement du terme source qui permet la prise en compte de l'influence de la topographie. On obtient notamment la propriété dite équilibre permettant de conserver exactement certains états stationnaires, appelés état du "lac au repos". Des versions 1D et 2D sur maillages non-structurés de ces méthodes ont été étudiées et implémentées dans un cadre volumes finis. Enfin, une extension vers des méthodes ordres élevés Galerkin discontinue a été proposée en 1D avec des limiteurs classiques ainsi que combinée avec une boucle MOOD de limitation a posteriori.

LA TRANSITION A LA TURBULENCE D'ECOULEMENTS EN GEOMETRIE COMPLEXE EST ETUDIEE PAR SIMULATION NUMERIQUE COMPRESSIBLE BIDIMENSIONNELLE ET TRIDIMENSIONNELLE. ON PRESENTE UNE SIMULATION DES GRANDES ECHELLES D'UN JET ROND D'AIR, COMPAREE A UNE EXPERIENCE EN LABORATOIRE ET A UNE AUTRE SIMULATION NUMERIQUE INCOMPRESSIBLE. CETTE COMPARAISON NOUS A PERMIS DE VALIDER NOTRE OUTIL NUMERIQUE. NOUS MONTRONS UN BON ACCORD ENTRE LES RESULTATS EXPERIMENTAUX ET NOTRE SIMULATION NUMERIQUE EN ZONE TRANSITIONNELLE. ON APPLIQUE NOTRE OUTIL NUMERIQUE A UNE GEOMETRIE COMPLEXE ET INDUSTRIELLE: L'ENTREE D'AIR EN INCIDENCE. CETTE ETUDE EST EFFECTUEE AVEC UN CODE A GESTION MULTI-DOMAINE DES TRANSFERTS D'INFORMATIONS PAR METHODE AUX CARACTERISTIQUES. LES RESULTATS SONT COMPARES A CEUX D'UNE EXPERIENCE FAITE AU CEAT DE POITIERS. NOUS AVONS MIS EN EVIDENCE LES TOURBILLONS DE KELVIN-HELMHOLTZ, AINSI QUE LES TOURBILLONS LONGITUDINAUX, RESPONSABLES DES FLUCTUATIONS DE PRESSION AU NIVEAU DU PREMIER ETAGE DU COMPRESSEUR. ENFIN, ON A REPRODUIT UNE EXPERIENCE MAGNETOHYDRODYNAMIQUE FAITE AU LABORATOIRE MADYLAM, EN GEOMETRIE BIDIMENSIONNELLE RONDE. L'INFLUENCE DU CHAMP MAGNETIQUE SUR UN ECOULEMENT DE MERCURE SOUMIS AUX FORCES DE LAPLACE CREE UN CISAILLEMENT AU NIVEAU DE L'ELECTRODE. ON PRESENTE UNE EVOLUTION TEMPORELLE DE LA VORTICITE ET DE LA PRESSION, AINSI QUE LES STATISTIQUES OBTENUES EN REGIME PERMANENT. NOUS MONTRONS DES RESULTATS QUALITATIFS SATISFAISANTS, MAIS PAR CONTRE SUR LES STATISTIQUES, NOUS SURESTIMONS LES VALEURS MOYENNES ET FLUCTUANTES. NOUS MONTRONS QUE LES CARACTERISTIQUES DE L'ECOULEMENT NE SONT PAS BIDIMENSIONNELLES, MAIS DOIVENT ETRE TRIDIMENSIONNELLES

L'ETUDE PORTE SUR L'ELABORATION D'UNE STRATEGIE NUMERIQUE EFFICACE PERMETTANT LA MODELISATION DES ECOULEMENTS NON DECOLLES DE FLUIDE COMPRESSIBLE VISQUEUX AUTOUR DE CONFIGURATIONS COMPLEXES (PROFILS D'AILE EN 2D, AILES OU AVIONS COMPLETS EN 3D). LA STRATEGIE S'APPUIE SUR UNE APPROCHE COMBINANT A LA FOIS UN SOLVEUR FLUIDE PARFAIT (POTENTIEL AVEC CORRECTION ENTROPIQUE OU EULER CONSERVATIVE) ET UN SOLVEUR COUCHE LIMITE CARACTERISANT LES EFFETS VISQUEUX PROCHES DE LA PAROI. LE RESULTAT DU CALCUL DE COUCHE LIMITE ENGENDRE L'EVALUATION D'UNE VITESSE NORMALE DITE DE SOUFFLAGE QUI SERT ENSUITE DE CONDITION A LA LIMITE DE PAROI DANS LE CALCUL FLUIDE PARFAIT. EN ITERANT AINSI LES CALCULS, ON OBTIENT UN SOLVEUR BIVALENT PEU COUTEUX, SE SUBSTITUANT A UN SOLVEUR NAVIER-STOKES, POUR DES ECOULEMENTS BIDIMENSIONNELS NON DECOLLES, ET PERMETTANT DE REALISER DES CALCULS D'ECOULEMENTS AUTOUR D'AVIONS COMPLETS. LA METHODOLOGIE AINSI DEVELOPPEE, MET EN EVIDENCE LA ROBUSTESSE DES DEUX SOLVEURS FLUIDE PARFAIT A LA PRISE EN COMPTE DE CONDITIONS AUX LIMITES DE PAROI NON HOMOGENES. LES NOMBREUSES EXPERIENCES NUMERIQUES EFFECTUEES, TANT ENE BIDIMENSIONNEL QU'EN TRIDIMENSIONNEL, FOURNISSENT DES RESULTATS D'UNE PRECISION SATISFAISANTE ET VALIDENT AINSI LA METHODE

L'OBJET DE CETTE THESE EST DE CONCEVOIR UNE METHODE NUMERIQUE CAPABLE DE SIMULER DES JETS LIQUIDES A HAUTE VITESSE AINSI QUE LEUR IMPACT SUR UN MATERIAU A USINER. NOUS PROPOSONS D'UTILISER UNE APPROCHE QUASI-CONSERVATIVE AINSI QU'UN SCHEMA ROBUSTE (LE SCHEMA HLL) AFIN D'OBTENIR UN CODE S'AFFRANCHISSANT DES OSCILLATIONS NUMERIQUES QUI PEUVENT APPARAITRE AUX INTERFACES. LE MEME TYPE DE PROBLEME PEUT SE PRODUIRE LE LONG DES LIGNES DE GLISSEMENT OU NOUS DEVELOPPONS DES CORRECTIONS APPROPRIEES. LA METHODE NUMERIQUE AINSI CONSTRUITE PERMET DES CALCULS D'IMPACTS DE GOUTTELETTES ET DE JETS. ELLE EST COMPLETEMENT EULERIENNE ET CAPTURE AUTOMATIQUEMENT LES INTERFACES, SI BIEN QU'ELLE S'INTEGRE AISEMENT DANS UN ENVIRONNEMENT DE MAILLAGE ADAPTATIF CARTESIEN. CETTE THESE COMPORTE TROIS PARTIES, LA PREMIERE EXPOSE LA MODELISATION ET LES METHODES NUMERIQUES QUI PERMETTENT NOS SIMULATIONS. LA DEUXIEME PRESENTE DES RESULTATS DE CALCULS NUMERIQUES POUR DES IMPACTS DE GOUTTELETTES ET DE JETS. LA TROISIEME PARTIE DECRIT L'IMPLANTATION DE NOTRE SOLVEUR DANS UN ENVIRONNEMENT DE MAILLAGE ADAPTATIF CARTESIEN.