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

CE TRAVAIL A POUR OBJECTIF SIMULER ET ETUDIER LES ECOULEMENTS COMPRESSIBLES TURBULENTS AXISYMETRIQUES D'ARRIERE-CORPS. POUR PRENDRE EN COMPTE L'EFFET DE LA TURBULENCE, LE MODELE K-E DE CHIEN ET LE MODELE DES EQUATIONS DU TENSEUR DE REYNOLDS DE HANJALIC ET LAUNDER ONT ETE INCORPORES DANS UN CODE NAVIER-STOKES EXISTANT. CE CODE EST BASE SUR UNE GENERALISATION DU SOLVEUR DE RIEMANN APPROCHE DE ROE, A ADAPTER POUR LES MODELES AUX EQUATIONS DE TRANSPORT TURBULENT. UNE VALIDATION DE LA PARTIE PARABOLIQUE SUR DES COUCHES LIMITES COMPRESSIBLES ET DE LA PARTIE HYPERBOLIQUE SUR DES ECOULEMENTS DERRIERE UNE MARCHE DESCENDANTE PLANE ET D'ARRIERE-CORPS AXISYMETRIQUE A ETE EFFECTUEE. LE SCHEMA DE ROE NE PERMETTANT PAS DE CONSERVER LA POSITIVITE DU TENSEUR DE REYNOLDS, LE SCHEMA HLLE A ETE DONC INTRODUIT. AINSI, NOUS AVONS DISCUTE, A TRAVERS UNE SIMULATION D'UN ECOULEMENT AXISYMETRIQUE D'ARRIERE-CORPS EN COMPARANT AVEC LES DONNEES EXPERIMENTALES DE HERRIN ET DUTTON, LA PERFORMANCE DES MODELES DE TURBULENCE ET LEUR COMPORTEMENT, AINSI QUE LA PERFORMANCE DES SCHEMAS AVEC LE MEME MODELE AU SECOND ORDRE

Dans le contexte de la prédiction d'écoulements aérodynamiques, l'élaboration d'une stratégie de modélisation hiérarchique permet d'envisager conjointement la simulation précise des phénomènes physiques par une approche "Haute-Fidélité" (HF) et l'intégration de modèles de faible complexité numérique au sein de procédures d'optimisation ou de contrôle. Une loi constitutive est proposée pour la modélisation HF des contraintes turbulentes dans le cadre de l'approche statistique Organised Eddy Simulation grâce à l'introduction d'une viscosité de turbulence tensorielle permettant de capturer les non-linéarités apparaissant entre tenseurs de déformation moyenne et d'anisotropie des contraintes turbulentes dans les régions hors-équilibre. Parallèlement, une approche de modélisation d'ordre réduit fondée sur la projection de Galerkin des équations de Navier-Stokes compressibles, sur une base extraite par Décomposition Orthogonale aux valeurs Propres (POD), est développée.

L'ECOULEMENT INSTATIONNAIRE TURBULENT AUTOUR DU PROFIL D'AILE A FORTE INCIDENCE EST ETUDIE PAR L'ANALYSE PHYSIQUE, SIMULATION NUMERIQUE, AINSI QUE PAR LA COMPARAISON AVEC LES DONNEES EXPERIMENTALES. L'APPROCHE DE MACROSIMULATION DE STRUCTURES ORGANISEES (OES, ORGANIZED EDDY SIMULATION ) EST ADOPTEE EN COUPLAGE AVEC LE MODELE DE TURBULENCE A DEUX EQUATIONS, DANS LE PRESENT TRAVAIL. UNE ETUDE DU CHOIX DU TRAITEMENT PARIETAL APPROPRIE EST EFFECTUEE AU SEIN DU MODELE DE TURBULENCE. UNE ANALYSE DETAILLEE DE LA DEPENDANCE DE LA SOLUTION VIS-A-VIS DE PARAMETRES NUMERIQUES, ET NOTAMMENT DE MAILLAGES EST PRESENTEE. L'ECOULEMENT INSTATIONNAIRE TURBULENT AUTOUR DU PROFIL D'AILE IMMOBILE A DES INCIDENCES DE 12, 14 ET 20 EST ETUDIE. L'INSTATIONNARITE NATURELLE DE L'ECOULEMENT EST PRISE EN COMPTE PAR L'APPROCHE DE MODELISATION DES STRUCTURES ORGANISEES. LE MODELE A DEUX EQUATIONS INCLUANT LE TRAITEMENT PARIETAL A FAIBLE NOMBRE DE REYNOLDS BASE SUR LA FONCTION D'AMORTISSEMENT CHOISIE EST IMPLEMENTE SUR LE MAILLAGE APPROPRIE. L'AUGMENTATION DISCRETE DE L'INCIDENCE PERMET L'OBSERVATION DE L'EVOLUTION DU DECOLLEMENT DE LA COUCHE LIMITE, AINSI QUE L'APPARITION DE L'ALLEE TOURBILLONNAIRE ALTERNEE DE VON KARMAN ET DES TOURBILLONS DE KELVIN-HELMHOLTZ. L'ETUDE DU PROFIL OSCILLANT AUTOUR DE L'INCIDENCE MOYENNE DE 12 PERMET D'ESTIMER L'IMPACT DE L'INSTATIONNARITE FORCEE SUR L'ECOULEMENT AUTOUR DU PROFIL, AINSI QUE SON INTERACTION AVEC L'INSTATIONNARITE NATURELLE. DEUX TYPES D'OSCILLATIONS SONT CONSIDERES : LE MOUVEMENT DE TYPE TANGAGE ET LE MOUVEMENT COMBINE DE TAMIS/TANGAGE. L'APPARITION DU TOURBILLON DE DECROCHAGE ET LE PHENOMENE DE RELAMINARISATION SONT OBSERVES. CE DERNIER SE MANIFESTE PAR L'APPARITION DE BOUCLES D'HYSTERESIS SUR LES COURBES DES PARAMETRES AERODYNAMIQUES. LE MODELE DU PREMIER ORDRE A DEUX EQUATIONS PERMET UNE PREDICTION CORRECTE DES PHENOMENES PROPRES AUX ECOULEMENTS INSTATIONNAIRES AUTOUR DE PROFILS D'AILE ET UN RELATIVEMENT BON ACCORD AVEC L'EXPERIENCE.

DANS CETTE ETUDE, UNE METHODE NUMERIQUE DE HAUTE PRECISION EST MISE AU POINT POUR LA SIMULATION D'ECOULEMENTS COMPRESSIBLES INSTATIONNAIRES TURBULENTS EN AERODYNAMIQUE. UN SCHEMA CENTRE A CINQ POINTS, PRECIS A L'ORDRE TROIS POUR LES EQUATIONS D'EULER, A ETE CONSTRUIT EN CORRIGEANT L'ERREUR DISPERSIVE DU SCHEMA CLASSIQUE D'ORDRE DEUX DE JAMESON. CE SCHEMA A ETE ENSUITE ETENDU AUX MAILLAGES CURVILIGNES BI- ET TRIDIMENSIONNELS PAR UNE METHODE DE VOLUMES FINIS DE HAUTE PRECISION. LE FLUX NUMERIQUE EST CONSTRUIT A L'AIDE D'OPERATEURS DE DISCRETISATION PONDERES, QUI PERMETTENT DE TENIR COMPTE DES DEFORMATIONS DU MAILLAGE. LE SCHEMA OBTENU EST PRECIS A L'ORDRE TROIS SUR DES MAILLAGES MOYENNEMENT DEFORMES ET A L'ORDRE DEUX SUR DES MAILLAGES TRES DEFORMES. LA METHODE DEVELOPPEE A ETE ETENDUE AUX EQUATIONS DE NAVIER-STOKES EN TRAITANT LES TERMES VISQUEUX DE FACON SIMPLE. LE SCHEMA PROPOSE A ETE PROGRAMME DANS UN CODE INDUSTRIEL ET IL A ETE VALIDE PAR DES COMPARAISONS SYSTEMATIQUES AVEC LE SCHEMA DE JAMESON, PROGRAMME DANS LE MEME CODE : LE SCHEMA DE HAUTE PRECISION PERMET DE REDUIRE DE FACON SIGNIFICATIVE L'ERREUR NUMERIQUE AU BORD D'ATTAQUE DE PROFILS ET D'AILES. ENSUITE, LE NOUVEAU SCHEMA A ETE APPLIQUE A L'ETUDE D'ECOULEMENTS AUTOUR DE PROFILS OSCILLANTS, MODELISANT UNE SECTION DE PALE DE ROTOR D'HELICOPTERE EN VOL D'AVANCEMENT. UN PROBLEME DIFFICILE DE DECROCHAGE DU A L'INTERACTION CHOC/COUCHE LIMITE SUR UN PROFIL OSCILLANT EN REGIME TRANSSONIQUE VIENT DEMONTRER LA NECESSITE DE SCHEMAS DE HAUTE PRECISION POUR LA SIMULATION D'ECOULEMENTS INSTATIONNAIRES COMPLEXES. IN FINE, LE NOUVEAU SCHEMA EST APPLIQUE A L'ETUDE DU DECROCHAGE DYNAMIQUE PROFOND D'UN PROFIL OSCILLANT A TRES GRANDE INCIDENCE. POUR CHAQUE APPLICATION, LES INFLUENCES DU MODELE DE TURBULENCE, DU MAILLAGE DE CALCUL ET DE L'ORDRE DE PRECISION DU SCHEMA SONT DISCUTEES.

L'écoulement turbulent autour d'un cylindre circulaire à grand nombre de Reynolds est analysé expérimentalement en utilisant les techniques PIV, PIV Stéréoscopique et PIV haute cadence. Les nombres de Reynolds étudiés correspondent à l'entrée de l'écoulement dans le régime critique, marqué par la chute brutale de traînée. Le cylindre est placé en milieu confiné afin de permettre des comparaisons avec des simulations numériques effectuées sur un domaine réaliste. Une base de données permettant la validation et l'amélioration des modèles de turbulence pour les approches statistiques avancées est ainsi fournie. Une attention particulière est portée à la décomposition du mouvement en une composante organisée et une composante aléatoire. Pour ceci, la décomposition en moyenne de phase, puis la décomposition en modes propres orthogonaux (POD), sont appliquées et comparées. Les contributions de chacune des composantes du mouvement au tenseur des contraintes turbulentes sont ensuite analysées.