UN CODE DESTINE AU CALCUL DES ECOULEMENTS TURBULENTS AUTOUR DES BATIMENTS, BAPTISE CHENSI, A ETE DEVELOPPE ET VALIDE SUR DIFFERENTS CAS REELS. LA DISCRETISATION DES TERMES DE CONVECTION DANS LES EQUATIONS DE LA QUANTITE DE MOUVEMENT ET DE TRANSPORT DE K- EST EFFECTUEE AVEC UN SCHEMA EN AMONT PONDERE DANS LE PLAN DE CALCUL, LE COUPLAGE VITESSE-PRESSION UTILISE LA METHODE DE COMPRESSIBILITE ARTIFICIELLE, L'ALGORITHME EST EXPLICITE EN TEMPS, SA CONVERGENCE EST ACCELEREE PAR L'INTRODUCTION D'UN PAS DE TEMPS EVOLUTIF ET LOCAL. LE MODELE DE TURBULENCE EST DU TYPE DE K-. ON UTILISE UNE LOI DE PAROI AFIN D'EVITER LA RESOLUTION DU SYSTEME DES EQUATIONS POUR LA SOUS-COUCHE VISQUEUSE. LE CODE A ETE APPLIQUE AU DEVELOPPEMENT D'UNE COUCHE LIMITE DU A UN CHANGEMENT DE RUGOSITE, A UN ECOULEMENT SUR UNE MARCHE DESCENDANTE, SUR UNE MARCHE MONTANTE ET AU-DESSUS D'UN BLOC BIDIMENSIONNEL. ENFIN LES CAS D'UN ECOULEMENT TRIDIMENSIONNEL AUTOUR D'UN OBSTACLE CUBIQUE ET PARALLELEPIPEDIQUE ONT ETE TRAITES. LA STRUCTURE DE L'ECOULEMENT, LE CHAMP DE VITESSE ET LE CHAMP DE PRESSION SONT CORRECTEMENT PREDITS PAR LE MODELE. CEPENDANT LE CALCUL NUMERIQUE SOUS-ESTIME L'ENERGIE TURBULENTE. L'EXTENSION DU SILLAGE DANS LES SIMULATIONS TRIDIMENSIONNELLES EST SOUS-ESTIMEE PAR LA SIMULATION. LES DIFFERENTES SOURCES D'INCERTITUDE ONT ETE ANALYSEES

La simulation numérique d'un écoulement visqueux incompressible autour d'un obstacle est réalisée par la résolution des équations de navier stokes bimensionnelles incompressibles instationnaires par une méthode spectrale. Pour tenir compte des conditions aux limites on utilise un opérateur de projection qu'on applique au champ de vitesse. cette méthode s'inspire de celle introduite par Basdevant et Sadourny en (1984). les résultats d'expériences a moyens et grands nombres de Reynolds ont montre la validité du modèle en comparaison avec des résultats numériques et physiques.

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

La simulation numérique directe des écoulements turbulents est un outil fondamental en mécanique des fluides car elle permet de connaître précisément les propriétés cinématiques de l'écoulement en tout point et à chaque instant. Cependant, en raison de la discrétisation spatiale d'ordre élevé nécessaire à la résolution de la cascade énergétique, cette approche est souvent limitée à des géométries simples. Ce travail propose une méthode de décomposition de domaine pour la simulation numérique d'écoulements turbulents. Elle repose sur une résolution directe multidomaine en formulation vitesse/pression des équations de Navier-Stokes. La méthode proposée résout efficacement (calcul parallèle) les équations de Navier-Stokes incompressibles sur des géométries complexes, telles que des quadrangles inter-connectés, tout en conservant les propriétés de haute précision. De nombreux tests de validation sont proposés: précisions spatiale et temporelle, cavité entraînée, marche descendante et écoulement autour d'un barreau de section carrée. Une simulation numérique directe d'un écoulement turbulent dans une conduite de section carrée est également réalisée. Enfin, une part importante de ce travail concerne l'étude d'un jet plan impactant sur une surface plane et sur un obstacle parallélépipédique pour le nombre de Reynolds Re=3000. Une analyse détaillée des propriétés statistique, spectrale et structurelle est menée, afin d'apporter une contribution à la compréhension des phénomènes physiques complexes de ces écoulements. Enfin, une étude du front de propagation de la température lors d'un refroidissement rapide d'une lamelle solide par un jet plan laminaire est présentée.

UNE ETUDE STATISTIQUE ET TOPOLOGIQUE D'UN ECOULEMENT TURBULENT COMPLEXE EN AVAL D'UNE MARCHE DESCENDANTE A ETE REALISEE AU MOYEN DES TECHNIQUES DE SIMULATION DIRECTE ET DE SIMULATION DES GRANDES ECHELLES. NOUS AVONS MIS EN EVIDENCE LA FORMATION DES STRUCTURES COHERENTES QUI COMPOSENT LA COUCHE DE MELANGE BIDIMENSIONNELLE. LA SIMULATION DIRECTE A ETE UTILISEE DANS CE CAS. LES ECOULEMENTS TRIDIMENSIONNELS ONT ETE ETUDIES A L'AIDE DE LA SIMULATION DES GRANDES ECHELLES. DANS CE CAS LE MODELE SOUS-MAILLE EST UNE ADAPTATION LOCALE DANS L'ESPACE PHYSIQUE DU CONCEPT DE VISCOSITE TURBULENTE DANS L'ESPACE DE FOURIER. LA STATISTIQUE DE LA TURBULENCE EST VALIDEE AU MOYEN DES RESULTATS EXPERIMENTAUX CORRESPONDANT A UNE MARCHE PETITE. DES CALCULS SUPPLEMENTAIRES SUR UNE MARCHE HAUTE ONT PERMIS L'ETUDE TOPOLOGIQUE DES INSTABILITES PRIMAIRES ET SECONDAIRES LIEES A LA COUCHE DE MELANGE