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 THESE, NOUS NOUS SOMMES INTERESSES A LA SIMULATION NUMERIQUE DE PHENOMENES D'INTERACTION FLUIDE-STRUCTURE. PLUS PRECISEMENT NOUS AVONS MIS AU POINT DES METHODES D'INTEGRATION DECALEE, QUE NOUS AVONS ENSUITE APPLIQUEES A DES CAS AEROELASTIQUES 2D NON-LINEAIRES. LES METHODES COURAMMENT UTILISEES POUR CES SIMULATIONS ONT D'ABORD ETE DETAILLEES. PARMI CELLES-CI, LES METHODES DECALEES PERMETTENT L'EMPLOI MODULAIRE DE PROCEDURES CLASSIQUES DANS CHAQUE SOUS-DOMAINE. SUR UN PROBLEME MODELE 1D LINEAIRE, NOUS AVONS MIS EN EVIDENCE L'IMPORTANCE DE L'ALGORITHME DE COUPLAGE UTILISE. NOUS AVONS DEMONTRE QUE DES CONSIDERATIONS ENERGETIQUES PERMETTENT D'ASSURER NON SEULEMENT LA STABILITE DU SYSTEME LINEAIRE COUPLE MAIS AUSSI L'INCONDITIONNELLE STABILITE DE NOUVEAUX SCHEMAS DECALES IMPLICITES, OU LE FLUIDE EST SOUS-CYCLE, CE QUI PERMET UNE SIMULATION PLUS RAPIDE, SURTOUT SUR MACHINES PARALLELES. LES IDEES GENERALES DEGAGEES DE CES MODELES ONT ENSUITE ETE APPLIQUEES A DES PROBLEMES 1D NON-LINEAIRES. NOUS AVONS ELABORE UNE NOUVELLE METHODE A INTERFACE DOUBLE, PLUS SOUPLE ET PLUS PRECISE, PUIS NOUS L'AVONS APPLIQUEE A DES PROBLEMES D'ECOULEMENTS 2D INSTATIONNAIRES NON-VISQUEUX AUTOUR DE STRUCTURES LINEAIRES MOBILES ET/OU DEFORMABLES EN FLOTTEMENT TRANSSONIQUE OU SUPERSONIQUE: DES RESULTATS PRECIS FURENT OBTENUS EN DE BIEN MOINDRES TEMPS DE CALCUL

L'ECOULEMENT INSTATIONNAIRE TURBULENT AUTOUR D'UN PROFIL D'AILE A FORTE INCIDENCE A ETE ETUDIE POUR UN FLUIDE INCOMPRESSIBLE. LA CONFIGURATION DU PROFIL D'AILE A ETE PLACEE DANS UN ECOULEMENT UNIFORME STATIONNAIRE OU OSCILLANT EN TAMIS. LA MISE EN UVRE D'UN LOGICIEL DE CALCUL EST EFFECTUEE EN MAILLAGE CURVILIGNE DE TYPE H, PRECIS AU SECOND ORDRE TEMPOREL ET SPATIAL, PERMETTANT LA RESOLUTION DES EQUATIONS DE NAVIER-STOKES COMPLETES POUR UN SILLAGE DE GEOMETRIE COMPLEXE, AFIN DE PREDETERMINER LA NAISSANCE D'INSTABILITES, DU DECOLLEMENT INSTATIONNAIRE ET DES STRUCTURES ORGANISEES DU SILLAGE A DES NOMBRES DE REYNOLDS MODERES. DANS CE CONTEXTE, UNE ANALYSE PHYSIQUE DES PROPRIETES DE L'ECOULEMENT A LA SORTIE DU DOMAINE A PERMIS DE DEVELOPPER DES CONDITIONS LIMITES NON-REFLECTIVES, EN ECOULEMENT INCOMPRESSIBLE ET INSTATIONNAIRE, DANS LE SYSTEME DES EQUATIONS DE NAVIER-STOKES. CE TYPE DE CONDITIONS LIMITES A PERMIS D'UTILISER UNE DISTANCE RAISONNABLE A L'AVAL ET D'EFFECTUER DES ECONOMIES NOTABLES SUR LE TEMPS DE CALCUL. AINSI IL MINIMISE LES EFFETS DE RETROACTION ET DE CONFINEMENT ARTIFICIEL. LE DEVELOPPEMENT DU MODELE DE TURBULENCE INSTATIONNAIRE A DEUX EQUATIONS EST EFFECTUE SUR LES EQUATIONS DE NAVIER-STOKES EN MOYENNE DE PHASE POUR LES NOMBRES DE REYNOLDS ELEVES. CETTE METHODOLOGIE PERMET DE CALCULER LES EFFETS FORTEMENT INSTATIONNAIRES (ECOULEMENT TURBULENT EN TAMIS) POUR UNE GEOMETRIE COMPLEXE. UN TRAITEMENT DE MODELE A FAIBLE NOMBRE DE REYNOLDS EST INCLUS. UNE MODELISATION POUR LE TERME DE PRODUCTION D'ENERGIE CINETIQUE TURBULENTE, EN FONCTION DU ROTATIONNEL LOCAL, EST PROPOSEE D'APRES DES CONSIDERATIONS PHYSIQUES, AFIN DE SURMONTER LE PROBLEME DE LA SURPRODUCTION D'ENERGIE TURBULENTE DANS LA REGION AU VOISINAGE DU POINT D'ARRET. L'ENSEMBLE DES RESULTATS OBTENUS MONTRENT L'EFFICACITE DE CETTE APPROCHE POUR LA PREDETERMINATION DES ECOULEMENTS INSTATIONNAIRES, DECOLLES EN SITUATION NATURELLE ET FORCEE

CETTE ETUDE EST CONSACREE A LA SIMULATION NUMERIQUE D'ECOULEMENTS INSTATIONNAIRES D'UN FLUIDE COMPRESSIBLE, VISQUEUX ET NON-VISQUEUX. DANS UN PREMIER TEMPS, ON SE PLACE DANS UN CADRE MONODIMENSIONNEL ET L'ON S'INTERESSE A LA MISE AU POINT DE SCHEMAS NUMERIQUES SUFFISAMMENT PRECIS POUR PERMETTRE LA SIMULATION D'ECOULEMENTS DE FLUIDES EN REGIME INSTATIONNAIRE. ON PRESENTE ENSUITE L'EXTENSION DE CES SCHEMAS AU CAS BIDIMENSIONNEL LORSQUE L'ECOULEMENT DU FLUIDE EST MODELISE PAR LE SYSTEME DES EQUATIONS DE NAVIER-STOKES. LA DISCRETISATION SPATIALE EST DU TYPE MIXTE ELEMENTS/VOLUMES FINIS SUR DES MAILLAGES TRIANGULAIRES NON-STRUCTURES. ON ETUDIE ENFIN LA MISE EN UVRE DE SIMULATIONS NUMERIQUES SUR UNE ARCHITECTURE MASSIVEMENT PARALLELE DE TYPE SYNCHRONE. DANS LE CAS OU L'ON UTILISE DES GRILLES DE CALCUL IRREGULIERES, UNE ATTENTION PLUS PARTICULIERE DOIT ETRE PORTEE A LA GESTION DES COMMUNICATIONS SI L'ON NE VEUT PAS VOIR LES PERFORMANCES DE CE TYPE DE MACHINE SE DEGRADER DE FACON DRAMATIQUE

L'ETUDE NUMERIQUE D'ECOULEMENTS INSTATIONNAIRES, EXTERNES, DECOLLES, DE FLUIDE VISQUEUX INCOMPRESSIBLE AUTOUR D'UN PROFIL D'AILE A GRANDE INCIDENCE A ETE EFFECTUE PAR LA TECHNIQUE DES DIFFERENCES FINIES ET LA SIMULATION DES GRANDES ECHELLES. LES EQUATIONS DE NAVIER-STOKES SONT RESOLUES EN FORMULATION VITESSE-TOURBILLON (3D) ET FONCTION DE COURANT-TOURBILLON (2D) SUR DES MAILLAGES CURVILIGNES ORTHOGONAUX. LES METHODES NUMERIQUES ONT ETE EVALUEES POUR UN FAIBLE NOMBRE DE REYNOLDS PAR COMPARAISON AVEC L'EXPERIENCE. DIFFERENTS MODELES DE SOUS-MAILLE ONT ETE TESTES EN DEUX DIMENSIONS ET A GRAND NOMBRE DE REYNOLDS. UN NOUVEAU MODELE APPELE GRADIENT DU TOURBILLON EST PROPOSE DANS LE CAS D'ECOULEMENTS BI-DIMENSIONNELS. UN CALCUL 3D EST EFFECTUE A GRAND NOMBRE DE REYNOLDS ET ON COMPARE LES RESULTATS DES SIMULATIONS DES GRANDES ECHELLES 2D ET 3D AVEC DES RESULTATS EXPERIMENTAUX.

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.