Cette thèse contribue à la simulation numérique des écoulements d'un fluide compressible modélisé par les équations de Euler et Navier-Stokes: étude d'un schéma d'ordre élévé basé sur une matrice de masse, modélisation des écoulement turbulents compressibles à très haut Reynolds, développement des estimateurs et correcteurs d'erreurs a posteriori et a priori, et adaptation de maillage anisotrope pour les fonctionnelles d'observation. Concernant la prédiction des écoulements turbulents, on s'est intéressé aux modèles hybrides de type RANS/LES comportant les nouveautés suivantes: traitement des tourbillons de grande échelle utilisant la formulation VMS (Variational Multi-Scale) et du RANS employé sur la paroi sur une distance imposée via une zone de protection conçue afin d'éviter le phénomène assez commun apelé "grid induced model depletion". Le niveau de viscosité du modèle VMS-LES est de plus controllé par un procédé de double filtre dynamique. La seconde partie concerne l'adaptation de maillage anisotrope pour mieux observer une fontionnelle d'observation. Les estimations a priori sont réalises pour le modèle des équations d'Euler et Navier-Stokes en instationnaire en 2D et 3D. A partir de ces estimations on sait définir les maillages optimaux au cours du calcul instationnaire, en fonction de l'état et de l'état adjoint. Le système d'optimalité est discrétisé et résolu à l'aide d'une méthode de point fixe instationnaire global, comportant une stratégie de stockage/recalcul pour le couplage état/ état adjoint. Des applications à la propagation d'ondes de choc et d'ondes acoustiques sont présentées.

Les torches à arcs plasma sont actuellement utilisées dans de nombreuses applications industrielles pour atteindre des températures nettement supérieures (> 5 000 K) à celle obtenues par la combustion (3 000 K max). Une technologie plasma triphasée avec des électrodes en graphite est en cours de développement au Centre PERSEE MINES ParisTech. Cette technologie diffère sensiblement des technologies à courant continu traditionnelles et vise à dépasser certaines limites des systèmes actuels en termes de robustesse, de coûts d'équipement et d'exploitation pour des applications liées à la conversion et la valorisation de biomasse et déchets. Dans le but d'améliorer la compréhension des phénomènes physiques instationnaires intervenant dans les décharges triphasées, une étude menée en parallèle sur le plan théorique et sur le plan expérimental a été conduite en conditions non réactives (azote et gaz de synthèse). Sur un plan expérimental, cette étude s'est appuyée sur des analyses réalisées avec une caméra ultra rapide (100 000 images par seconde) et l'analyse des signaux électriques. Sur un plan théorique, cette étude a consisté à développer un modèle Magnéto-Hydro-Dynamique (MHD) 3D instationnaire de la zone d'arc dans l'environnement du logiciel Code Saturne® et à effectuer une étude paramétrique basée sur le courant, la fréquence et le débit de gaz plasma. Deux configurations : électrodes coplanaires et électrodes parallèles ont été étudiées. Cette étude a permis de mettre en évidence l'influence des phénomènes électromagnétiques et hydrodynamiques sur le déplacement de l'arc. Dans la configuration coplanaire, les jets aux électrodes semblent jouer un rôle prépondérant sur le mouvement des arcs, les transferts de chaleur dans l'espace inter-électrodes et l'amorçage des arcs. Dans la configuration parallèle, le mouvement des canaux chauds semble être le paramètre dominant. La confrontation des résultats théoriques et expérimentaux a montré un très bon accord à la fois au niveau du mouvement des arcs et des signaux électriques.

Cette thèse vise à une meilleure compréhension de la dynamique du voilier et à la validation des outils numériques de prédiction de performances et d'optimisation par l'étude expérimentale in situ du problème aéro-élastique d'un gréement. Une instrumentation est développée sur un voilier de 8m de type J80 pour la mesure dynamique des efforts dans le gréement, de la forme des voiles en navigation, du vent et des attitudes du bateau. Un effort particulier est apporté à la mesure des caractéristiques géométriques et mécaniques des éléments du gréement, la calibration des capteurs et au système d'acquisition des données. Les principaux résultats montrent que le voilier instrumenté est un outil adapté pour les mesures instationnaires et soulignent l'amplitude de variation d'effort rencontrée en mer (20 à 50% de l'effort moyen dans une houle modérée). En outre, les variations du signal d'effort sont déphasées avec l'angle d'assiette, créant un phénomène d'hystérésis. Le comportement dynamique d'un voilier en mouvement diffère ainsi de l'approche quasi-statique. Les simulations numériques proviennent du code ARAVANTI, couplage implicite d'un code structure éléments finis ARA et d'un code fluide parfait, limitant son domaine de validité aux allures de près Les résultats de simulation sont très proches des cas stationnaires et concordent bien avec les mesures en instationnaire dans une houle de face. L'expérimentation numérique d'un gréement soumis à des oscillations harmoniques en tangage souligne l'importance de l'approche Interaction Fluide Structure (IFS) et montre que l'énergie échangée par le système avec la houle est reliée à la fréquence réduite et l'amplitude du mouvement. Certaines informations n'étant pas disponibles sur le voilier instrumenté, une expérience contrôlée en laboratoire est développée. Elle consiste en un carré de tissu tenu par deux lattes en oscillation forcée. Les mesures sur cette " voile oscillante " permettent d'étudier les phénomènes IFS avec décollement et sont utilisées pour la validation du couplage ARA-ISIS entre un code fluide Navier-Stokes (RANS) et le même code structure.

On se propose d'analyser le régime transitoire de décrochage à l'aide de la simulation numérique instationnaire de type DNS. Cette approche permet de reproduire avec fidélité l'écoulement dans la région critique de Bulbe de Décollement Laminaire au bord d'attaque, encore impossible à modéliser ou mesurer avec précision. Après une étape de validation, la sensibilité du BDL au niveau de turbulence extérieure est étudiée et comparée favorablement à celle établie récemment dans la littérature. La phase d'établissement du décollement massif depuis le BDL est reproduite en réalisant de petites variations d'incidence à travers l'angle critique d'apparition du décrochage. Ce raisonnement "aux petites perturbations" permet de reproduire l'éclatement du BDL communément rattaché au décrochage statique. La déstabilisation de la région de BDL est alors étudiée à l'aide d'une base de données instationnaires et moyennes qui permet pour la première fois de rendre compte des déformations 3D du BDL. Conservant ce protocole, et faisant varier la vitesse du profil, on est en mesure d'évaluer l'influence de cette dernière sur le régime transitoire. Des mouvements de rotation de plus forte amplitude angulaire ont permis de mettre en évidence un mécanisme de décrochage sensiblement différent du précédent. La couche de mélange surplombant le BDL s'enroule alors pour donner naissance à un tourbillon énergétique (Leading Edge Vortex), communément associé au décrochage dynamique. Enfin, l'analyse du champ de vitesse moyen a permis de valider un critère empirique d'apparition de l'éclatement du BDL, qui s'est révélé pertinent aussi bien pour les mouvements de faible amplitude que ceux plus amples.

Ces travaux de recherche ont pour but d’évaluer la méthode dite de la "Simulation auxEchelles Adaptées" (SAS pour Scale-Adaptive Simulation). Cette approche coïncide avec uneapproche RANS classique dans les zones pariétales attachées et adapte le niveau de viscositéturbulente dans les zones décollées pour y permettre une résolution partielle des structures turbulentes.Dans une première partie, une analyse théorique du modèle SAS original a été menéeet a permis de développer une correction visant à favoriser l’adaptation du niveau de viscositéturbulente dans les zones sièges d’instabilités de type Kelvin-Helmholtz. Le modèle ainsi corrigéest nommé SAS-αL. Les modèles SAS et SAS-αL ont été implantés dans le code de calculNavier-Stokes elsA de l’ONERA. À l’issue de cette étape, trois cas académiques d’écoulementsturbulents instationnaires, cylindre à grand nombre de Reynolds, marche descendante et cavitétranssonique, ont été simulés grâce aux trois modèles de turbulence SST, SAS et SAS-αL. Outreune comparaison aux bases de données expérimentales disponibles, une attention particulièrea été portée à l’influence de paramètres numériques tels que des schémas numériques d’ordreélevé. Enfin, afin d’étudier la viabilité de l’approche SAS dans un contexte industriel, les troismodèles de turbulence ont été testés sur une configuration issue de l’industrie aéronautique etcorrespondant à la sortie d’air chaud d’un système de dégivrage des nacelles d’avion. La comparaisondes prévisions obtenues avec les modèles SST, SAS et SAS-αL aux données expérimentalesobtenues à l’ONERA a permis de montrer un gain de précision grâce à l’emploi de l’approcheSAS et ce pour un coût de calcul compatible avec un cycle de conception industrielle. / This research work is meant to assess an upgraded URANS approach, namely the Scale-Adaptive Simulation (SAS). This method is similar to a conventional RANS approach (namelythe SSTmodel) in attached areas and is able to adapt the eddy-viscosity level in detached areas toensure the resolution, at least partially, of the turbulent structures. In a first part of this researchwork, an improvement of the SAS approach is suggestedto allowa better sensitivity of themodelto instabilities such as Kelvin-Helmholtz ones. This "improved" model is referred to as SAS-αLmodel. Both SAS and SAS-αL models were implemented in the ONERA Navier-Stokes solverelsA and both of themaswell as the SSTmodelwere tested on academic test cases : a cylinder in acrossflowat a high Reynolds number, a backward-facing step flowcorresponding to theDriver&Seegmiller experiment and the transonic flow over the M219 cavity experimentally investigatedby de Henshaw. The influence of the numerical parameters was deeply investigated and particularattention was paid to the high-order space-discretization schemes effects. The reliabilityof the SAS approach in an industrial framework was assessed on an aeronautic configurationnamely a nacelle de-icing device. Comparisons between the threemodels (SST, SAS and SAS-αL)and an experimental database available at ONERA - The French Aerospace Lab have shown thebetter accuracy of the SAS approach as well as the high potential of the SAS-αL model.

Ce travail tente d'analyser la réponse d'une jet axisymétrique turbulent à une excitation fluidique stationnaire et instationnaire lorsque le contenu fréquentiel et aziumutal (!,m) de la perturbation est maîtrisé. Le dispositif de contrôle utilisé est composé de 16 microjets ronds répartis sur le bord de fuite de la tuyère. L'utilisation des microjets provoque une réduction du champ acoustique rayonné (particulièrement pour le cas de contrôle stationnaire). Le champ aérodynamique est ensuite sondé grâce à des mesures fil chaud et PIV stéréoscopique résolue en temps. L'excitation instationnaire permet d'utiliser les moyennes de phase afin d'effectuer une décomposition triple du champ de vitesse. L'étude de la composante cyclique de la “réponse du jet” montre une synchronisation spatio-temporelle importante sur une grande étendue spatiale. En d'autres mots, le forçage a une grande autorité déterministe sur l'écoulement. De plus, la comparaison de la composante cyclique de la réponse du jet avec la théorie de la stabilité linéaire indique qu'il existe des ondes d'instabilité hydrodynamique au sein du jet. L'analyse du jet contrôlé par injection fluidique stationnaire montre ensuite comment l'effet du contrôle peut être expliqué par la déformation du champ moyen conduisant à la réduction du taux de croissance des ondes d'instabilité dans le jet. Cette déformation est dûe à l'introduction d'un couple de paramètre (nombre d'onde/fréquences) pour lequel le champ moyen de l'écoulement est stable. La réponse du jet étant turbulente, cela implique que les tensions de Reynolds déforment le champ moyen de manière à ce que les modes les plus instables aient des taux de croissance plus faibles. / This work investigates the response of an axisymetric turbulent jet to steady and unsteady fluidic florcing where the azimuthal wavenumber-frequency (!,m) content of the perturbation is well known. The control setup is composed of 16 round microjets azimutally distributed around the nozzle lip. Such actuation can lead to a decrease in the acoustic energy radiated by the jet (especially for the steady case). The aerodynamic fied is investigated using hotwire measurements and time-resolved stereoscopic PIV. Using the unsteady forcing, phase-averaging is possible, and this allows the implementation of a triple decomposition of the measurements. Examination of the cyclic component of the flow response shows that a non-negligible phase-locked fluctuation is obtained over a large spatial extent, in other words, the actuation has good deterministic control authority over the flow. Furthermore, comparison of the cyclic component of the flow response with Linear Stability Theory supports the idea that the jet response comprises linear hydrodynamic instability waves. Subsequent analysis of jets controlled by steady fluidic actuation shows how the control effect can be explained by a mean-flow modification that leads to the reduction of instability-wave growth rates ; the mean flow modification is argued to be due to the introduction of azimuthal wavenumber-frequency pairs to which the mean flow is stable. The response is therefore turbulent, and involves Reynolds stresses which deform the mean-field such that the most unstable modes have lower growth rates.

L’aérodynamique latérale des véhicules automobiles suscite de nos jours de plus en plus d’intérêt de la part des constructeurs. L’automobiliste est en effet soumis quotidiennement à de forts courants d’air latéraux, que ce soit lors du dépassement d’un autre véhicule, ou alorsen passant dans un couloir de vent du à la topographie du terrain (passage devant un espace entre deux immeubles par exemple). Les efforts aérodynamiques mis en jeu dans ces situations peuvent provoquer des mouvements non désirés du véhicule, pouvant avoir des conséquences dramatiques si le conducteur se laisse surprendre. Des études expérimentales reproduisant les effets d’un dérapage dynamique ont mis en évidence des phénomènes transitoires importants mettant à défaut les modèles stationnaires généralement pratiqués par les constructeurs pour qualifier le comportement de leurs véhicules en présence de dérapage. Les mécanismes responsables de ces phénomènes transitoires sont encore mal connus de la communauté scientifique. Ce travail propose d’approfondir ce sujet au travers de l’étude de l’aérodynamique d’un véhicule terrestre fixe soumis à un vent longitudinal et à une rafale de vent latéral. Le but principal est d’identifier les structures tourbillonnaires au moyen de mesures PIV et de calculs numériques des champs de vitesse autour d’une maquette automobile et de les corréler aux efforts aérodynamiques. Un accord entre l’ISAT, composante de l’Université de Bourgogne, et l’Institut Supérieurde l’Aéronautique et de l’Espace (ISAE) de Toulouse a permis de mener l’étude avec les ressources de cet établissement. Le moyen d’essai principal, créé par l’ISAE, est le banc« rafale latérale », constitué d’une soufflerie principale et d’une soufflerie secondaire, dont la sortie à volet déferlants (« Mexican Wave ») est inspirée de l’approche proposée par Ryan et Dominy (2000). L’analyse expérimentale a été effectuée à l’aide de la PIV résolue en temps et stéréoscopique, et d’une balance dard instationnaire à cinq composantes. Un banc« numérique » identique a été constitué à l’aide du logiciel FLUENT©, pour des calculs 3D. De plus, un modèle 2D annexe, basé sur la méthode « meshless », a été développé pour de futures investigations, en raison de sa robustesse pour des problèmes à fortes discontinuités et sa bonne adaptabilité aux problèmes avec frontières mobiles.Une première phase de ce travail a consisté en la mise au point des bancs expérimental et numérique, avec génération d’un champ de dérapage homogène, de 21° dans la zone de mesure. L’évolution du dérapage en chaque point respecte bien la forme d’un créneau imposé par la rafale. Pour l’analyse des efforts, deux géométries de maquette ont été étudiées, à savoir un corps de Windsor à culot droit générant, pour un écoulement longitudinal, des structures de sillage bidimensionnelles, et son homologue à culot incliné de 25°, générant des tourbillons « cigare ». Des pics d’efforts ont été observés à l’arrivée de la rafale, tout comme la littérature le prédit. Pour ce qui est du coefficient du moment de lacet, les sursauts sont de 29 % et 19 % respectivement pour la maquette à culot droit et celle à culot incliné, par rapport aux valeurs stationnaires. Concernant le coefficient de force de dérive, ils sont de 10 % et 14 %, respectivement. Lors de nos essais, ces efforts se sont établis après 5.5 longueurs de maquette. Afin d’expliquer la différence de comportement entre les deux maquettes en termes d’efforts, l’évolution temporelle des tourbillons nommés, dans ce mémoire, ΓA, ΓB, ΓC et Γ1 à été discutée. Il en est ressorti une forte corrélation entre la circulation du tourbillon ΓA, le plusénergétique, naissant à l’avant du flanc sous le vent de la maquette, et les efforts latéraux, de sorte que ce tourbillon serait le meilleur témoin des phénomènes instationnaires mis en jeu dans l’étude de l’effet du vent latéral. [...] / The automotive manufacturers are nowadays more and more interested in crosswind aerodynamics. Indeed, the driver is subjected every day to strong side air flows, for example when overtaking another vehicle or when passing through a lateral wind wall, generated by terrain topography (as in the case of the passage near the empty space between two buildings).The aerodynamic efforts generated in these situations can lead to undesired lateral deviations,which can be dramatic if the driver is surprised. Different experimental studies, reproducing the effects of a dynamic yaw angle, pointed out the issues of the steady methods, commonly used to qualify the vehicle crosswind behaviour. Little is still known about the physics behind these unsteady phenomena. This is the main topic of this work, by studying the aerodynamics of a fixed vehicle subjected to both a longitudinal flow and a side wind gust. The goal is the identification of the near-vehicle vortex structures, by means of PIV measurements and CFD calculations, and their correlation with the evolution of the efforts. An agreement between the ISAT, a department of the University of Burgundy, and the ISAE of Toulouse, permitted to carry out this research with the resources of the latter laboratory. The work focuses on the use of the “rafale latérale” (side gust) test bench, made up with a main wind tunnel connected with an auxiliary one by means of a shutter system,whose opening is held by a “Mexican Wave” law. This approach is inspired by the work of Dominy and Ryan (2000). The experimental analysis was carried out by means of Time-Resolved and stereoscopic PIV, and by a five components unsteady balance as well. Anidentical test bench was numerically reproduced with the 3D CFD software FLUENT©.Moreover, an additional 2D CFD model, based on the meshless method, has been developed for future studies. This kind of approximation method has been chosen for its robustness innon-continuous problems and because of its adaptability when moving boundaries are needed.The first phase of this work consisted on wind tunnels set-up, both for the real test bench and for the CFD model. The yaw angle field is homogeneous, 21° in the measurement region. The yaw angle evolution, at given point, respects the step wise behaviour, imposed by the gust passage. As far as the efforts are concerned, two versions of the Windsor body car model were studied, that is a squareback geometry, generating, for longitudinal flows, 2D wakestructures, and a fastback geometry (rear window inclined by 25°), producing cone-liketrailing vortices. Force overshoots were seen after the gust arrival, as seen in literature. In particular, the yaw moment coefficient overshoots are 29% and 19% higher than the steady yaw angle tests, for the squareback and the fastback geometries, respectively. If the side forceis concerned, the entities of these overshoots are 10% and 14%, respectively. Our testspointed out that efforts establish after the vehicle has driven 5.5 times its length in thecrosswind. In order to explain the different behaviour of the two geometries, it is discussed about the unsteady evolution of the vortices called, in this report, ΓA, ΓB, ΓC et Γ1. A strong correlation between the side efforts and the circulation of the most energetic vortex, ΓA,having its origin in the front leeward side of the vehicle. The ΓA vortex is so the best index for the study of the crosswind unsteady phenomena. The coupled analysis between vortex structures and efforts confirmed the presence of a higher side force for the squareback geometry. The inverted effect has been observed for the yaw moment

Cette thèse s'intéresse aux simulations dépendantes du temps impliquant des géometries fixes ou mobiles. Ce type de simulations est l'objet d'attentes grandissantes de la part des industriels, qui souhaiteraient voir réaliser ce type de calculs de façon systématique au sein de leurs centres de recherche, ce qui n'est clairement pas le cas à l'heure actuelle. Ce travail tente de satisfaire en partie cette demande et vise notamment à améliorer la précision ainsi que l'efficacité en termes de temps de calcul des algorithmes actuellement utilisés dans ce contexte. Les méthodes d'adaptation de maillage anisotrope par prescription d'un champ de métriques, qui ont aujourd'hui atteint une certaine maturité, notamment dans leur application aux simulations stationnaires, constituent une piste très prometteuse pour l'amélioration des calculs évoluant en temps, mais leur extension dans ce contexte est loin d'être triviale. Quant à leur utilisation sur les simulations en géométries mobiles, seules quelques tentatives peuvent être répertoriées, et très peu portent sur des problèmes réalistes en trois dimensions. Cette étude présente plusieurs nouveautés sur ces questions, notamment l'extension de l'adaptation de maillage multi-échelles par champ de métriques aux problèmes instationnaires en géométries fixes et mobiles. Par ailleurs, essentiellement dans une optique de réduction des temps de calculs, une stratégie originale à été adoptée pour réaliser des calculs impliquant des maillages mobiles. Notamment, il est démontré par la pratique dans cette thèse qu'il possible de déplacer des objets en trois dimensions sur de grandes distances en maintenant le nombre de sommets du maillage constant, c'est-à-dire en limitant les types d'opérations de modification de maillage autorisés. Il en résulte un gain conséquent en terme de temps de calcul aussi bien au niveau du déplacement de maillage qu'au niveau de la résolution numérique. Par ailleurs, un nouveau schéma est proposé qui permet de gérer les changements de connectivité du maillage de manière cohérente avec la description Arbitrary-Lagrangian-Eulerian des équations physiques. La plupart de ces nouvelles méthodes ont été appliquées à la simulation d'écoulements fluides compressibles autour de géometries complexes en deux et trois dimensions d'espace.