Le bruit émis par les nouvelles architectures de moteurs aéronautiques a été considérablement réduit dans les dernières années. Les différentes sources de bruit ont été identifiées et pour la plupart réduites. Cependant, la contribution relative du bruit de combustion au bruit global a augmenté progressivement avec la décroissance des autres sources. Deux mécanismes de génération de bruit de combustion ont été identifiés : le bruit direct qui est produit par des fluctuations du dégagement de chaleur dû à la combustion, et le bruit indirect qui est généré par l’accélération des spots d’entropie. Dans ce travail, les mécanismes de génération et propagation du bruit entropique sont étudiés par des simulations numériques aux grandes échelles (en anglais LES) et par des modèles analytiques. Dans un premier temps, une configuration simplifiée du phénomène est étudiée : des spots d’entropie sont créés par des résistances chauffantes et ensuite accélérés par une tuyère pour générer du bruit indirect. Cette configuration a été simulée et ses résultats validés par des campagnes expérimentales. Ensuite, la simulation numérique est utilisée pour mieux comprendre les mécanismes de génération du bruit indirect et ses interactions avec des effets visqueux et non visqueux. Dans une seconde partie, une configuration de turbine haute pression à un seul étage est utilisée pour étudier le bruit indirect d’une façon plus réaliste. Dans les deux parties de cette thèse, les résultats numériques sont comparés à des théories analytiques pour mieux comprendre les avantages et inconvénients d’une méthode par rapport à l’autre. / Combustion noise is increasing its relative contribution to aircraft noise, while other sources are being reduced and new low-NOx emission combustion chambers being built. Two mechanisms are responsible for this noise source: direct noise in which acoustic waves are generated by the flame and propagate to the outlet of the aero-engine, and indirect noise, where entropy waves generate noise as they are accelerated and decelerated in the turbine stages. In this work, the analytical models used for the propagation of waves through non-homogeneous flows, including the generation of indirect noise, are revised and extended. In the first part, the quasi-1D case is studied, extending the analytical method to non-zero frequencies and validating the results with numerical methods and experimental data. In the second part, the 2D method for the case of compact turbine blades is studied and validated using numerical simulations of a rotating blade and of a complete turbine stage. Finally, in the third part of this thesis, these models are combined with reactive and compressible Large Eddy Simulations (LES) of combustion chambers to build a hybrid approach, named CHORUS, able to predict combustion noise.

Résumé : Les nanotubes de carbone mono-parois (SWCNTs) sont très peu dispersibles dans les solvants et ils ont besoin d'être chimiquement modifiés avant leur utilisation dans beaucoup d'applications. Ce travail se concentre sur la synthèse du matériau des SWCNTs chimiquement modifié par une approche in situ. Les objectifs principaux de cette recherche sont : I) explorer le procédé chimique in situ pendant la synthèse des SWCNTs et 2) examiner de manière approfondie l'effet de l'environnement réactif sur les SWCNTs. Les effets du type de catalyseur et son contenu sur le produit fini des SWCNTs, synthétisé par plasma thermique inductif (PTI), ont été étudiés pour remplacer le cobalt (Co) toxique dans la matière première. À cet égard, trois mélanges de catalyseurs différents (c.-à-d. Ni-Y203, Ni-Co-Y203, et Ni-Mo-Y203) ont été utilisés. Les résultats expérimentaux ont montré que le type de catalyseur affecte la qualité des SWCNTs. Une qualité similaire peut être produite lorsque la même quantité de Co est remplacée par le Ni. En outre, des résultats observés dans les travaux expérimentaux ont été explicités par les résultats des calculs thermodynamiques. La therrnogravimétrie (TG) a été utilisée tout au long du travail pour caractériser les échantillons de SWCNTs. La TG a tout d'abord été normalisée par l'étude des effets des trois principaux paramètres instrumentaux (rampe de température, RT, la masse initiale de l'échantillon, MI, et le débit de gaz, D) sur le T, et largeur à mi-hauteur (LMH) obtenu à partir de graphiques TG et TG dérivés de noir de carbone, respectivement. Par conséquent, un plan factoriel à deux niveaux a été prévu. L'analyse statistique a montré que l'effet de RT, MI, et à un degré moindre D est significatif sur la LMH et négligeable sur Tonss. Une méthodologie a ensuite été développée sur la base de la synthèse des SWCNT en utilisant le système PTI, à travers une approche chimique in situ. L'ammoniac (NH3) a été choisi et injecté à contre-courant dans le réacteur PTI à trois débits différents et en utilisant quatre types de buses différentes. La simulation numérique a indiqué un meilleur mélange du NH3 dans le réacteur PTI lorsqu'une buse particulière a été utilisée. Les résultats expérimentaux montrent l'augmentation d'intensité de D-bande dans les spectres Raman d'échantillons SWCNTs lors de l'injection du NH3. Le NH3 pourrait augmenter la teneur en azote du produit fini de SWCNTs jusqu'à 10 fois. L'échantillon des SWCNTs traitée avec 15% vol de NH3 a montré une dispersion accrue dans le diméthylformamide et l'isopropanol. Les nanostructures de carbone en forme d'oignon et plane, ont aussi été observées. Une caractérisation complémentaire sur l'échantillon des SWCNTs traités par NH3 à 15% vol., a indiqué une modification de la surface des nanotubes, où des tubes métalliques ont montré une plus grande réactivité avec NH3 que les semi-conducteurs. Le modèle, y compris le champ d'écoulement thermique du réacteur et la cinétique de décomposition thermique de NH3 a suggéré une modification de surface des SWCNTs en deux étapes dans laquelle les nanotubes réagissent premièrement avec les espèces intermédiaires de H et de NH2. Le NH3 s'adsorbe ensuite chimiquement sur les nanotubes. Le modèle a également suggéré que les espèces intermédiaires comme le NNH et le N2H2 jouent un rôle principalement en conduisant la décomposition du NH3 plutôt que la modification chimique des SWCNTs. // Abstract : Pristine single-walled carbon nanotubes (SWCNTs) are poorly dispersible and insoluble in many solvents and need to be chemically modified prior to their use in many applications. This work is focused on the investigation of the synthesis of chemically modified SWCNTs material through an in situ approach. The main objectives of the presented research are: I) to explore the in situ chemical process during the synthesis of SWCNT and 2) to closely examine the effect of a reactive environment on SWCNTs. Effects of the catalyst type and content on the SWCNTs final product, synthesized by induction thermal plasma (1TP), were studied to replace toxic cobalt (Co) in the feedstock. In this regard, three different catalyst mixtures (i.e. Ni-Y203, Ni-Co-Y203, and Ni-Mo-Y203) were used. Experimental results showed that the catalyst type affects the quality of the SWCNT final product. Similar quality SWCNTs can be produced when the same amount of Co was replaced by Ni. Moreover, the results observed in this experimental work were further explained by thermodynamic calculation results. Thermogravimetry (TG) was used throughout the work to characterize the SWCNTs product. TG was firstly standardized by studying the effects of three main instrumental parameters (temperature ramp, TR, initial mass of the sample, 1M, and gas flow rate, FR) on the Lise, and full-width half maximum (FWHM) obtained from TG and derivative TG graphs of carbon black, respectively. Therefore, a two-level factorial statistical design was performed. The statistical analysis showed that the effect of TR, IM, and to a lower extent, FR, is significant on FWHM and insignificant on T01, 1. A methodology was then developed based upon the SWCNTs synthesis using the 1TP system, through an in situ chemistry approach. Ammonia (NH3) was selected and counter-currently injected into the ITP reactor at three different flow rates and by four different nozzle designs. Numerical simulation indicated a better mixing of NH3 in the ITP reactor when a certain nozzle was used. The experimental results showed the increase of D-band intensity in the Raman spectra of SWCNT samples upon the NH3 injection. NH3 could increase the nitrogen content of the SWCNTs final product up to 10 times. The SWCNTs sample treated with 15 vol% NH3 showed an enhanced dispersibility in Dimethylformamide and Isopropanol. Onion-like and planar carbon nanostructures were also observed. Complementary characterization on the SWCNT samples treated by 15 vol% NH3 indicated the surface modification of nanotubes. Metallic tubes showed a higher reactivity with NH3 than semiconducting ones. The model including the reactor thermo-flow field and NH3 thermal decomposition kinetics suggested a two-step SWCNT surface modification in which nanotubes firstly react with H and NH2 intermediates and later, NH3 chemisorbs on the nanotubes. The model also suggested that the intermediate species, like NNH and N2H2, play a rote primarily in driving the NH3 decomposition rather than the chemical modification of SWCNTs. [symboles non conformes]

Nous présentons dans cette étude le développement de la convection à partir de différentes perturbations de l'état conductif d'une couche fluide confinée dans une cavité cylindrique, chauffée par le bas et avec une surface supérieure libre. La discrétisation spatiale du domaine repose sur la méthode des éléments spectraux et les itérations temporelles sont assurées par une méthode splitting.Au déclenchement de la convection, les structures convectives correspondent à des modes de Fourier, et les seuils critiques dépendent du rapport de forme de la cavité, et des nombres de Biotet de Marangoni qui caractérisent la surface libre. Les transitions d'écoulements au-delà du seuil primaire sont caractérisées quantitativement en fonction du nombre de Rayleigh pour différentes valeurs du nombre de Biot et Ma = 0. Les résultats présentés sont obtenus en résolvant l'ensemble des équations non-linéaires de conservation à travers une méthode de continuation. Lorsque la convection se déclenche sous la forme d'un mode axisymétrique m = 0, l'évolution non-linéaire montre la coexistence de différentes structures convectives, des structures axisymétriques avec écoulement montant ou descendant au centre de la cavité et des structures correspondant à des combinaisons de modes qui apparaissent sur des branches secondaires sous-critiques.L'action d'un champ magnétique constant est ensuite étudiée pour des fluides conducteurs dans une même configuration comprenant une surface supérieure libre. Nous montrons l'effet stabilisateur du champ magnétique sur les seuils primaires ainsi que son action sélective sur les différents modes de convection. Nous analysons l'évolution des structures convectives au delà de ces seuils et montrons comment le champ magnétique modifie les transitions entre ces structures.En soumettant le bain fondu à un champ magnétique tournant, le mouvement de rotation du fluide se superpose aux mouvements de convection thermique et on observe une diminution des fluctuations de température et un retard du déclenchement de l'instabilité de Rayleigh-Bénard(lorsque les deux parois haut/bas du bain sont rigides). La rotation influe sur ce déclenchement qui de stationnaire devient oscillatoire, à l'exception du mode m = 0 de Fourier, pour qui la transition reste stationnaire jusqu'à une certaine valeur critique du nombre de Taylor magnétique.La dynamique de l'écoulement axisymétrique de part et d'autre de cette valeur critique sera étudiée en détail. / The growth of thermal convection out of different perturbations of the conductive base state is investigated using a spectral element time-stepping code. The fluid is subject to a vertical heat transfer in a cylindrical cavity with an upper free surface corresponding to the so-called Rayleigh-Bénard-Marangoni situation and the heat exchange through the free surface is evaluated via the Biot number. The results of the stability diagrams show that the evolution of the primary thresholds are largely influenced by the Biot number, the Marangoni number, and the aspect ratio of the cavity. Flow transitions are elucidated in quantitative detail as a function of the Rayleigh number for different Biot numbers in the tension free limit Ma = 0. The results shown are obtained by solving the full nonlinear field equations numerically among a continuation method. When an axisymmetric m = 0 Fourier mode is obtained at onset, the non-linear evolution shows the coexistence of different convective structures, the axisymmetric structures with up-ow or down-ow at the center and mixed-mode structures which appear on secondary subcritical branches. The action of a constant magnetic field is then considered for melts in the same type of configuration with an upper free surface. We show the global stabilizing effect of the magnetic field on the primary bifurcation thresholds and the selective effect on the different instability modes. We also analyze the evolution of the convective structures above the thresholds and show how the magnetic field modifies the transitions between these structures. When applying a magnetic body forcing in the azimuthal direction (RMF), one can damp the unavoidable thermal fluctuations inside the melt and delay the transition to the Rayleigh-Bénard instability (for rigid-rigid circular plates at top and bottom). The rotation effect also changes the transitions from steady to oscillatory, except for the m = 0 Fourier mode where the transitionis first steady until a critical Taylor number and then becomes oscillatory. The dynamics of the transitions to the axisymmetric flow, below and above this value of critical magnetic Taylor number, is particularly interesting and will be described.

Cette thèse est consacrée à l'étude numérique des interactions non-linéaires entre des vagues et un corps rigide perçant la surface libre. La méthode développée repose sur un modèle d'éléments de frontière qui réduit la dimensionnalité du problème d'une dimension. Dans un premier temps, un modèle2D est appliqué à des géométries simples et permet de démontrer la pertinence de l'approche envisagée pour la prédiction des mouvements d'une structure flottante soumise à des vagues monochromatiques régulières. Dans un second temps, en nous inspirant d'un modèle potentiel non-linéaire 3D développé par Grilli textit{et al.}~cite{grilli2001fully}, nous proposons une généralisation de la méthode pour des maillages triangulaires non-structurés de surfaces 3D. Le modèle développé permet de traiter des configurations arbitraires de plusieurs cylindres verticaux en interaction avec les vagues. Nous présentons des cas de validation de nature académique qui permettent d'apprécier le comportement du modèle numérique. Puis nous nous tournons vers l'application visée par EDF R&D, qui concerne le dimensionnement d'éoliennes off-shore flottantes. Un flotteur de type semi-submersible est évalué à l'aide du modèle non-linéaire / This PhD work is devoted to the study of nonlinear interactions between waves and floating rigid structures. The developed model relies on a boundary element method which reduces the dimensionality of the problem by one. First, a 2D model is applied to basic geometries and allows us to demonstrate the validity of the method for predicting the motion of a floating structrure subject to incoming monochromatic regular waves. Secondly, getting inspired by the 3D fully nonlinear potential flow model of Grilli textit{et al.}~cite{grilli2001fully}, we propose a novel model which generalizes the method for unstructured triangular meshes of 3D surfaces. The proposed model is able to deal with arbitrary configurations of multiple vertical cylinders interacting with the waves. We present academic validation test cases which show how the model works and behaves. Finally, we study situations of interest for EDF R&D related to floating off-shore wind turbines. A semi-submersible platform is evaluated with the nonlinear model

La première partie de ce travail a porté sur la mesure du gradient de vitesse et la direction de l'écoulement local de liquide dans un empilement de sphères. Cette étude a été menée grâce à des micro-électrodes tri-segmentées implantées à la paroi de sphères supports. Ces expériences originales ont en particulier permis de caractériser les recirculations dans le pore suivant les conditions opératoires. La seconde partie a porté sur la simulation numérique de ce même écoulement local. Un écoulement monophasique dans une géométrie complexe représentant exactement le milieu poreux a été simulé. Ces simulations ont été confrontées aux expériences; l'accord est satisfaisant, surtout à faible vitesse. L'inadéquation des modèles de turbulence usuels a été montrée dans notre cas. Enfin, des mesures par Vélocimétrie Laser Doppler ont montré l'intérêt de cette technique, bien que des améliorations soient nécessaires pour obtenir des résultats exploitables sur l'ensemble du milieu / The first part of this work deals with local measurement of flow direction and velocity gradients in a structured packing of spheres. Tri-segmented micro-electrodes flush mounted at the wall of supporting spheres are used for the measurements. Among other results, these original experiments allowed to charatcerize the flow recirculation within the pore.In the second part, the local liquid flow in the packing was numerically computed. More specifically, the local flow of liquid in a complex geometry was modelled and simulated. Numerical results are in agreement with experiments, especially for low flow rates. Usual turbulence models were shown to be inappropriate for the present case. Finally, some Laser Doppler Velocimetry measurements were carried out and showed the interesting potential of this technique, but improvements are yet to be made to get velocities in the bulk of the bed

Cette thèse s’inscrit dans le cadre du développement d’un procédé de décarbonation directe du méthane par voie plasma pour la coproduction de noir de carbone et d’hydrogène. Ce procédé est particulièrement intéressant dans la contexte d'un mix électrique faiblement carboné en offrant une solution pour diminuer drastiquement les émissions des procédés actuels de production de noir de carbone et d’hydrogène qui comptent parmi les procédés actuels les plus polluants en termes d'émissions de CO2. A court terme, la viabilité économique de ce procédé passe par la valorisation simultanée de ces deux produits : noir de carbone et hydrogène. À plus long terme, il pourrait représenter une réelle alternative à la capture et le stockage du CO2.Cette étude a pour but de proposer des méthodes numériques fiables et robustes afin de mieux comprendre, contrôler, voire optimiser les caractéristiques morphologiques des noirs de carbone issus de ce procédé, caractéristiques qui jouent un rôle primordial sur la qualité et les applications des noirs de carbone. Elle traite ainsi de l’évolution de systèmes carbonés en partant du combustible sous sa forme moléculaire jusqu’à la formation de nanoparticules puis de microstructures solides et aborde les phénomènes de nucléation, de croissance chimique, de croissance par coagulation, de maturation et d’agrégation.Des outils et des méthodes numériques ont ainsi pu être développés afin de simuler la formation de particules solides au sein d’un écoulement fluide. Ceux-ci purent être implémentés avec succès à un code CFD. Enfin des simulations numériques du procédé en question ont été réalisées en intégrant les phénomènes de transferts thermiques et de turbulence spécifiques aux plasmas thermiques. / This thesis takes part of the development of a direct decarbonation process of methane by plasma to produce both carbon black and hydrogen. This process is particularly interesting in an electrical mix context with low carbon emission. It proffers a solution to reduce drastically CO2 emissions rejected by the current carbon black and hydrogen ways of production, which are ones of the most polluting industrial processes.This study aims to develop reliable and robust numerical methods for a better understanding and a greater control of the morphologic features of the carbon black generated. These features play a key role in the quality and applications of the carbon black produced. This research retraces the evolution of the carbon structure from the molecules of the fuel to the formation of nanoparticles and solid microstructures. It tackles different phenomenon such as: nucleation, chemical growth, coagulation, maturity and aggregation.Numerical tools and methods were developed thereby and enable to simulate carbon particle formation. They were successfully implemented in a commercial CFD software. Eventually numerical simulation of the plasma process were performed, integrating heat transfers and turbulence.

La méthode LS-STAG est une méthode cartésienne pour le calcul d’écoulements incompressibles en géométries complexes, qui propose une discrétisation précise des équations de Navier-Stokes dans les cut-cells, cellules polyédriques de forme complexe créées par l’intersection du maillage cartésien avec la frontière du solide immergé. Originalement développée pour les géométries 2D, son extension aux géométries 3D se heurte au défi posé par le grand nombre de types de cut-cells (108) à considérer. Récemment, la méthode LS-STAG a été étendue aux géométries complexes 3D dont la frontière est parallèle à l’un des axes du repère cartésien, où sont uniquement présentes les contreparties extrudées des cut-cells 2D. Cette étude a notamment souligné deux points à élucider pour le développement d’une méthode totalement 3D : premièrement, le calcul des flux diffusifs par un simple schéma à deux points s’est révélé insuffisamment précis dans les cut-cells 3D-extrudées du fait de la non orthogonalité. Ensuite, l’implémentation de ces flux à la paroi, qui s’effectue en imposant une discrétisation distincte pour chaque type de cut-cell extrudée, se révèle trop complexe pour être étendue avec succès aux nombreux types supplémentaires de cut-cells 3D, et doit être simplifiée et rationalisée. Dans cette thèse, le premier point est résolu en utilisant l’outil des schémas diamants, d’abord étudié en 2D pour l’équation de la chaleur puis les équations de Navier-Stokes dans l’approximation de Boussinesq, puis étendu en 3D. En outre, les schémas diamants ont permis de revisiter intégralement la discrétisation du tenseur des contraintes des équations de Navier-Stokes, où disparaît le traitement au cas par cas selon la disposition de la frontière solide dans les cut-cells. Cela a permis d’aboutir à une discrétisation systématique, précise et algorithmiquement efficace pour les écoulements en géométries totalement 3D. La validation numérique de la méthode LS-STAG avec schémas diamants est présentée pour une série de cas tests en géométries complexes 2D et 3D. Sa précision est d’abord évaluée par comparaison avec des solutions analytiques en 2D, puis en 3D par la simulation d’un écoulement de Stokes entre deux sphères concentriques. La robustesse de la méthode est notamment mise en évidence par l’étude d’écoulements autour d’une sphère en rotation, dans les régimes laminaires (stationnaire et instationnaire), ainsi que pour un régime faiblement turbulent. / The LS-STAG method is a cartesian method for the computations of incompressible flows in complex geometries, which consists in an accurate discretisation of the Navier-Stokes equations in cut-cells, polyhedral cells with complex shape made by the intersection of cartesian mesh and the immersed boundary. Originally developed for 2D geometries, where only three types of generic cut-cells appear, its extension to 3D geometries has to deal with the large amount of cut-cells types (108). Recently, the LS-STAG method had been extended to 3D complex geometries whose boundary is parallel to an axis of the cartesian coordinate system, where there are only the extruded counterparts of 2D cut-cells. This study highlighted two points to deal with in order to develop a totally 3D method: firstly, the computation of diffusive fluxes by a simple 2-points scheme has shown to be insufficiently accurate in 3D-extruded cut-cells due to the non-orthogonality. In addition to that, implementation of these fluxes on the immersed boundary, which is done with a case by case discretisation according to the type of the cut-cells, appears to be too difficult for its successful extension to the several extra types of 3D cut-cells, and needs to be simplified and rationalized. In this thesis, the first point is solved by using the diamond scheme tool, firstly studied in 2D for the heat equation then for the Navier-Stokes equations in Boussinesq approximation, and finally extended to 3D. Moreover, the diamond schemes have been used to fully revisit the discretisation of shear stresses from Navier-Stokes equations, where the case by case procedure is removed. These modifications have permitted to come up with a systematic discretisation that is accurate and algorithmically efficient for flows in totally 3D geometries. The numerical validation of the LS-STAG method with diamond schemes is presented for a series of test cases in 2D and 3D complex geometries. The precision is firstly assessed by comparison with analytical solutions in 2D, then in 3D by the simulation of Stokes flow between two concentric spheres. The robustess of the method is highlighted by the simulations of flows past a rotating sphere, in laminar modes (steady and unsteady), as well as in a weakly turbulent mode.

Pour stopper un emballement thermique dans un réacteur chimique, un moyen efficace consiste à introduire une faible quantité d'un inhibiteur liquide appelé « killer » dans la cuve agitée. Tout au long de cette thèse, l'approche expérimentale a été fortement couplée à la modélisation numérique par Computational Fluid Dynamics (CFD). La première partie du manuscrit porte sur l'hydrodynamique des réacteurs partiellement chicanés incluant la prise en compte du vortex central qui se forme à leur surface. L'utilisation d'une approche numérique multiphasique, non-homogène a permis de modéliser la déformation de la surface-libre, et la faisabilité de cette méthode innovante a été démontrée par un très bon accord entre prédictions numériques et données expérimentales. Dans une deuxième partie, l'introduction d'un jet de liquide sur la surface libre a été couplée à l'hydrodynamique du réacteur. Les résultats numériques, obtenus avec une approche Eulerienne-Lagrangienne, ont également montré un bon accord avec les données expérimentales. Ces résultats ont permis de modéliser la trajectoire du jet, de quantifier sa pénétration dans la cuve agitée, et de définir de nouveaux critères de mélange. Enfin, les méthodes numériques validées à l'échelle pilote ont été étendues à l'échelle industrielle et ont permis de proposer des améliorations concrètes pour une meilleure sécurité des réacteurs industriels étudiés.

L’élimination des inclusions indésirables par flottation en cuve agitée est un procédé largement utilisé dans l’industrie de l’aluminium. Ce procédé consiste en l’injection de gaz au sein de la cuve par le biais d’un rotor ; les bulles ainsi formées captent, lors de leur ascension dans le métal liquide, les inclusions et les drainent en surface où elles sont éliminées mécaniquement.Dans le but de prédire l’évolution numérale de la population d’inclusions dans la cuve de flottation au cours du temps, un modèle 2D a été développé. Afin de rendre compte du transport convectif des inclusions au sein de l’écoulement diphasique métal liquide – bulles, des phénomènes de décantation, d’agglomération et de flottation, ce modèle couple la Mécanique des Fluides Numérique et la technique des Bilans de Population. Une attention particulière a été portée aux modèles d’agglomération et de flottation pour qu’ils soient applicables sur une grande plage de conditions de turbulence.Les simulations numériques réalisées, nous ont notamment permis de mettre en évidence l’importance de veiller à une bonne répartition des bulles dans la cuve, en particulier dans les zones de forte turbulence. En parallèle de ce travail de modélisation, une analyse de la littérature et des travaux expérimentaux sur métal liquide ont été réalisés, afin de mieux comprendre le rôle du chlore dans le procédé de flottation. Un effet de seuil a ainsi pu être mis en évidence : une certaine quantité de chlore injectée dans le métal est nécessaire pour que la flottation devienne efficace. La compréhension du rôle du chlore reste cependant partielle, mais des pistes de recherche prometteuses ont put être identifiées / Inclusions removal by flotation in mechanically agitated vessels is widely used in liquid aluminium treatments. This process consists of gas injection in a molten aluminium reactor using an impeller. Inclusions are attached to the bubbles while rising into the liquid bulk and released in the froth at the vessel surface.With the aim of both a better understanding of the physical processes acting during flotation and the optimization of the refining process, a 2D mathematical modelling of the flotation tank was set up. Transport phenomena, sedimentation, agglomeration and flotation of inclusions were considered here. The model couples Population Balance with convective transport of the inclusions, in the multiphase flow field. Simulations showed that a good dispersion of the bubbles into the tank, especially in the high turbulence zones, is essential for an efficient process. This situation was achieved with a strong agitation. Gas flow affected the process efficiency weakly.Theoretical and experimental studies on chlorine effects were done. It was noticed that a critical amount of chlorine into the melt is required for the flotation to be efficient. Complete understanding of chlorine role is not achieved yet, but promising ideas for future research have been exhibited

Le passage du régime laminaire au régime turbulent s’accompagne d’importantes modifications des propriétés physiques de l’écoulement. Une prévision précise du point du début de la transition laminaire/turbulent revêt donc une importance considérable dans de nombreux domaines pratiques. Lorsque l’intensité des perturbations extérieures est significative, c'est-à-dire dans le cas de couches limites se développant sur une paroi présentant des rugosités ou soumises à une forte turbulence résiduelle (sillage impactant), les mécanismes de formation et d’amplification des instabilités sont profondément modifiés. Ces perturbations sont les modes de Klebanoff (également appelés stries) qui s’amplifient et déclenchent la transition, qualifiée dans ce cas de Bypass. Ces stries sont très énergétiques, caractérisées par des fluctuations de vitesse très importantes (de l’ordre de 10% de la vitesse extérieure), alors que la couche limite conserve son caractère laminaire. La thèse proposée concerne la modélisation de ces stries via la résolution d’une équation de transport pour l’énergie cinétique dite laminaire. Dans un premier temps, le travail du candidat portera sur la modélisation des termes de production et de dissipation de l’énergie cinétique laminaire. Ceux-ci sont liés au processus de réceptivité de la couche limite vis-à-vis des perturbations extérieures et à la dynamique des modes de Klebanoff dans la zone laminaire. Pour ce faire, la thèse s’appuiera sur des études réalisées depuis plusieurs années au sein de l’unité ITAC sur la théorie des perturbations optimales ainsi que sur les travaux numériques et expérimentaux prévus dans le cadre d’un projet de recherche interne Onera. Classiquement cette équation de transport est couplée avec celles correspondant à l’énergie cinétique turbulente et à la dissipation, le mécanisme d’échange entre les énergies cinétiques laminaire et turbulente devra être soigneusement étudié : ce dernier pilote la transition vers la turbulence. Une attention particulière sera portée aux couches limites décollées et plus précisément à la prise en compte de la transition dans ces bulbes. Cette nouvelle modélisation innovante permettra l’amélioration d’une première approche pour le calcul de la transition bypass dans le solveur elsA, développé à l’Onera, et constituera une étape importante vers la mise en place de techniques de prévision de la transition pratiques et performantes. / This work aims to develop a new bypass-transition prediction model based on the Klebanoff modes dynamics. To represent these mode dynamics the Laminar Kinetic Energy (LKE) concept has been chosen, in order to model these mode energy with a new variable. A new deffinition is given to the LKE and a transport equation consequently derived to describe the Klebanoff modes growth and destabilisation. This equation is incorporated in a k-omega turbulence model as done by Walters & Cokljat, to give a three-equation kL-kT-omega formulation. This new model is written in a Reynolds-averaged Navier-Stokes (RANS) pattern and only uses local variables, it thus can be used in an industrial context.