Dans cette thèse, l'application de l'assimilation de données (AD) à la MFN est étudiée, avec comme objectif global de contribuer à l'amélioration de la prévision numérique d'écoulements complexes. L'AD consiste à fusionner les outils de prévision numérique avec des données expérimentales afin d'améliorer l'estimation des paramètres d'entrée du code MFN. Les aspects méthodologiques de l'AD et son application pour des études physiques sont tous deux examinés dans cette thèse. Dans un premier temps, l'AD est utilisée pour une étude théorique de la turbulence de grille. Un modèle spectral pour les écoulements turbulents homogènes et anisotropes est également proposé. Plusieurs méthodes d'AD sont ensuite implémentées pour un code MFN et appliquées à la reconstruction d'écoulements instationnaires et compressibles en présence d'incertitudes sur des paramètres d'entrée de grandes dimensions afin d'évaluer les forces et faiblesses respectives de ces techniques. Des stratégies pour le placement optimal de réseaux de capteurs sont élaborées afin d'améliorer les performances du processus d'AD. Enfin, l'AD est appliquée à l'identification de sources de polluants et à la reconstruction de conditions météorologiques pour des écoulements en milieu urbain prédits par Simulation des Grandes Echelles. / In this thesis, we investigate the use of various data assimilation (DA) techniques in the context of CFD, with the ultimate goal of enhancing the prediction of real-world flows. DA consists in merging numerical predictions and experimental observations in order to improve the estimation of the CFD solver inputs. Both methodological aspects of DA and its potential application to physics investigations are explored for various flow configurations. First, DA is considered for the theoretical analysis of grid turbulence decay. Fundamental aspects of anisotropic homogeneous turbulence are also investigated through spectral modelling. Various DA methodologies are deployed in conjunction with a Navier-Stokes solver and are assessed for the reconstruction of unsteady compressible flows with large control vectors. Sensor placement strategies are developed to enhance the performances of the DA process. Finally, a first application of DA to Large Eddy Simulations of full-scale urban flows is proposed with the aim of identifying source and wind parameters from concentration measurements.

L'étude des mécanismes de formation d'un spray à partir d'un écoulement liquide-gaz est un sujet de recherche fondamental, et dont les applications industrielles sont nombreuses, notamment dans le domaine de la combustion et de la propulsion. La simulation numérique de tels écoulements apparaît comme un complément indispensable aux expériences et aux théories pour aboutir à une compréhension et à une prédiction fine de tels phénomènes physiques. Durant cette étude nous avons développé une technique numérique de suivi d'interface couplée à un solveur des équations de Navier-Stokes afin d'étudier précisément la dynamique de l'interface qui sépare la phase liquide de la phase gazeuse. Nous décrivons en détail la méthode Level Set qui a été utilisée durant cette étude pour suivre le mouvement de l'interface, ainsi que les techniques numériques utilisées pour la résolution des équations de Navier-Stokes. Notamment différents schémas numériques sont testés pour que la précision du code de calcul développé soit optimal. Le couplage entre les équations du mouvement et le suivi d'interface peut être effectué par différentes techniques. L'utilisation de la Ghost Fluid Method pour réaliser ce couplage permet de préserver numériquement le caractère discontinu de l'interface liquide-gaz tout en proposant une discrétisation des équations précises et robustes. Les codes développés (2D, 3D parallélisation MPI) sont ensuite utilisés pour étudier des collisions de gouttes. Les comparaisons entre les résultats expérimentaux et les simulations numériques sont très probantes et mettent en valeur la prédictibilité des calculs effectués. D'autres configurations plus complexes sont étudiées dans le cadre d'une étude de faisabilité. L'étude de l'atomisation primaire d'un jet liquide continu en un spray de gouttelettes semble notamment être un domaine d'investigation particulièrement prometteur pour ce type de simulations. Enfin l'étude d'interfaces dites « réactives » comme la vaporisation d'un liquide et la propagation de flammes prémalangées a également été abordé en utilisant la Ghost Fluid Method pour imposer les conditions de saut spécifiques à ce type de phénomènes.

L'objet de cette thèse a été de développer un code de calcul pour la simulation d'écoulements diphasiques de fluides non miscibles, incompressibles et isothermes afin de l'appliquer au phénomène de fragmentation d'un jet liquide, et plus particulièrement à la désintégration assistée d'une nappe de liquide par deux écoulements d'air portés à haute vitesse.<br /><br />Les choix des hypothèses physiques et des schémas numériques associés ont été fait afin de respecter au mieux la physique complexe de brisure d'un jet liquide. La résolution des équations de Navier-Stokes incompressibles est faite de façon directe par le biais d'une méthode de projection. La méthode naissante et prometteuse Level-Set assure quant à elle le suivi numérique de la surface de séparation entre deux fluides non miscibles. Enfin, la méthode Ghost Fluid permet un traitement correct des discontinuités à la traversée de l'interface en préservant au niveau discret les conditions de saut entre équilibre des forces de capillarité, de pression et de viscosité. Les bonnes aptitudes de tels schémas numériques sont montrées à travers une batterie de cas tests académiques.<br /><br />Les mécanismes physiques mis en jeu lors de l'atomisation primaire d'un jet liquide ont fait l'objet de nombreuses études théoriques et expérimentales. Les instabilités se développant à la surface du liquide sont multiples et clairement tridimensionnelles. Actuellement et malgré beaucoup d'efforts de recherche, aucune théorie ou modèle ne rend compte rigoureusement de ce phénomène. Or, il est montré que l'outil numérique peut apporter une lumière nouvelle. Par exemple, l'impact de la couche limite gazeuse sur la fréquence d'oscillation est étudié. Nous avons également pu retrouver par la simulation la dynamique ligamentaire intrinsèque à la désintégration assistée d'un jet liquide.

Les caractéristiques d’un nuage de poussières avant l’activation de la source d’inflammation ont une grande influence sur la sensibilité et la sévérité de son explosion. Pour cette raison, le procédé de dispersion d’une poussière combustible qui se développe à l’intérieur des chambres d’explosion standardisés telles que le tube de Hartmann modifié et la sphère 20 L a été décrit par deux approches complémentaires. Dans un premier temps, une caractérisation expérimentale de l’évolution du nuage de poussière a identifié les niveaux de turbulence de l’écoulement de gaz, ainsi que les principales variations de la distribution de la taille des particules combustibles. Ces résultats ont été complétés par l’identification des niveaux de ségrégation de la poudre dispersée à l’intérieur des chambres d’explosion. Par la suite, une série de simulations CFD basée sur une approche Euler-Lagrange a été développée pour prédire le comportement du nuage. Cette étude a été réalisée grâce à la caractérisation des principaux mécanismes d’interaction tels que le transfert de quantité de mouvement (couplage bidirectionnel) et le phénomène de fragmentation des poudres. Deux types de particules microniques ont été spécifiquement étudiés : une poussière métallique, l’aluminium et un composé organique, l’amidon de blé. Ces résultats ont démontré qu’il est possible de définir ab initio, par simulation numérique, les conditions de tests les plus pertinentes (les plus réalistes ou les plus pénalisantes) en vue de la quantification de ces risques majeurs / The pre-ignition stage of an explosibility test determines the ignitability of a dust cloud as well as the main characteristics of the flame propagation. For this reason, the dispersion process of a combustible dust that develops inside the explosion chambers of the modified Hartmann tube and the 20 L sphere has been described by two complementary approaches. Initially, an experimental characterization of the evolution of the dust cloud identified the evolution of the turbulence levels of the gas flow along with the main variations of the particle size distribution of the combustible dust. These results were complemented by the identification of the segregation levels of the dispersed powder inside the explosion chambers. Thereafter, a set of CFD simulations based on the Euler-Lagrange formulation was developed to predict the behavior of the combustible dust cloud. This study was accomplished through the characterization of the main interaction mechanisms such as the momentum exchange (two-way coupling) and the dust fragmentation phenomenon. In this manner, the research work constituted for this thesis allowed determining the most appropriate conditions to ignite a dust cloud formed by a metallic (Aluminum micro-Al 42) or organic powders (wheat starch) in a case study. Thus, this research work presents a methodology that can be extended for the analyses of combustible dusts and the further development of the standard test methods

Cette thèse traite de la modélisation des réacteurs chimiques par la « modélisation compartimentale », qui consiste à diviser le système en un réseau d’une dizaine à quelques centaines de volumes interconnectés, appelés compartiments. La structure du réseau est déduite à partir d’informations provenant d’expériences de traçage, d’informations techniques sur le réacteur chimique, de simulations de mécanique des fluides numérique et des objectifs de la modélisation. Cette méthode procure un bon compromis entre temps de calcul et finesse des résultats. Quand ils sont correctement menés, les modèles à compartiments donnent des prédictions similaires, en termes de réactions chimiques, à ceux issus des simulations de mécanique des fluides numérique réactive avec un temps de calcul plus court et une représentation physique plus concrète du comportement du réacteur. Chaque étude issue de la littérature est consacrée à un réacteur spécifique avec une approche particulière qui ne peut pas être directement transposée sur un autre réacteur. L’objectif de cette thèse est d’apporter une contribution au développement d’une méthodologie la plus générale possible et de développer un outil de génération automatique et de résolution du système d’équations différentielles qui doit être résolu. Dans le premier chapitre, un état de l’art est réalisé, définissant le champ d’application de notre méthode, dans le but d’identifier les méthodes de découpage les plus pertinentes et les différentes méthodes pour calculer les échanges entre les compartiments. Dans un second chapitre, une méthode générale pour de la modélisation compartimentale est développée. Une approche polyvalente est proposée, consistant à découper le réacteur en tranches identiques. Le calcul des échanges entre compartiments, dus à la convection et la turbulence, est présenté en détail, avec la description des trois méthodes de calcul des échanges turbulents. Une interface a été développée permettant de construire n’importe quel réseau de compartiments. À partir de cette interface, les équations sont écrites et automatiquement résolues. La méthode est appliquée dans un troisième chapitre sur un cas défavorable au découpage en tranches. Cela a permis de tester les limites de cette approche. En particulier, deux points ont été étudiés : (1) l’applicabilité du découpage en tranches identiques et (2) la comparaison entre les méthodes de calcul des échanges turbulents. Le premier test a prouvé la robustesse de l’approche par division mais le second test n’a pas permis d’établir si une méthode de calcul est meilleure qu’une autre. Finalement, la méthode a été valorisée et transférée en implémentant les algorithmes développés dans un logiciel commercial. Ce logiciel permet de simuler la dispersion d’espèces réactives et non réactives (traceurs), dans un modèle contenant plusieurs centaines de compartiments organisés en tranches identiques / This PhD deals with modelling of chemical reactors with the “compartmental modelling” approach, which consists in dividing the system into a network from a dozen to several hundreds of interconnected volumes, called compartments. The structure of the network is deduced from tracer experiments, technical information about the chemical reactor and computational fluid dynamics flow simulations. This method provides a good compromise between computation time and results accuracy. When they are properly set-up, compartmental models give similar predictions, in terms of chemical reactions, as those of CFD simulations with a shorter calculation time and a more concrete representation of the reactor behavior. Every study from the literature is devoted to a specific reactor with a particular approach that cannot be straightforwardly transposed to other reactors. The aim of this PhD is to provide a contribution to the development of the most general possible methodology and to develop an automatic tool of generation and resolution of the differential equations system which must be solved. In the first chapter, a state of the art is proposed, defining the field of application of our method, in order to identify the most relevant division methods and the different methods to calculate the exchange between compartments. In the second chapter, a general methodology for compartmental modelling is developed. A versatile approach is proposed, consisting in dividing the reactor in identical slices. The calculation of exchange between compartments, both due to convection and turbulence, is presented in detail, with the description of three calculation methods for turbulent exchange. An interface has been developed, allowing to build any network of compartments. From this interface, the equations are written and solved automatically. The methodology is applied in the third chapter to an unfavorable case for slice cutting. This has allowed to test the limit of this approach. In particular, two points have been studied: (1) the applicability of division into identical slices and (2) the comparison between the turbulent exchange calculation methods. The first test has proved the robustness of the division approach but the second test has not allowed to establish whether one calculation method is better than another. Finally, the methodology has been promoted and transferred by implementing the developed algorithms within a commercial software. This software allows to simulate the dispersion of reactive and non-reactive (tracers) species, in model containing hundreds of compartments organized in identical slices

Les tourbillons sont des structures importantes pour une large gamme d'écoulements de fluides, notamment les sillages, l'interaction fluide-structure, les décollements de couche limite et la turbulence. Cependant, les méthodes numériques classiques n'arrivent généralement pas à donner une représentation précise des tourbillons. Ceci est principalement lié à la dissipation numérique des schémas qui, si elle n'est pas spécifiquement calibrée pour le calcul des écoulements tourbillonnaires, conduit à une diffusion artificielle très rapide des tourbillons dans les calculs. Parmi d'autres approches, la méthode "Vorticity Confinement" (VC) de J. Steinhoff permet de compenser la dissipation des schémas au sein des tourbillons en introduisant une anti-dissipation non-linéaire, mais elle n’est précise qu’au premier ordre. D’autre part, des progrès significatifs ont récemment été accomplis dans le développement de méthodes numériques d’ordre élevé. Celles-ci permettent de réduire ce problème de dissipation excessive, mais la diffusion des tourbillons reste importante pour de nombreuses applications. La présente étude vise à développer des extensions d’ordre élevé de la méthode VC pour réduire cette dissipation excessive des tourbillons, tout en préservant la précision d'ordre élevé des schémas. Tout d'abord, les schémas de confinement sont analysés dans le cas de l'équation de transport linéaire, à partir de discrétisations couplées et découplées en espace et en temps. Une analyse spectrale de ces schémas est effectuée analytiquement et numériquement en raison de leur caractère non linéaire. Elle montre des propriétés dispersives et dissipatives améliorées par rapport aux schémas linéaires de base à tous les ordres de précision. Dans un second temps, des schémas VC précis au troisième et cinquième ordre sont développés pour les équations de Navier-Stokes compressibles. Les termes correctifs restent conservatifs, invariants par rotation et indépendants du schéma de base, comme la formulation originale VC2. Les tests numériques valident l'ordre de précision et la capacité des extensions VC d’ordre élevé à réduire la dissipation dans les tourbillons. Enfin, les schémas avec VC sont appliqués au calcul des écoulements turbulents, dans une approche de simulation de grandes échelles implicite (ILES). Les schémas numériques avec VC présentent une résolvabilité améliorée par rapport à leur version linéaire de base, et montrent leur capacité à décrire de façon cohérente ces écoulements tourbillonnaires complexes. / Vortices are flow structures of primary interest in a wide range of fluid dynamics applications including wakes, fluid-structure interaction, flow separation and turbulence. Albeit their importance, standard Computational Fluid Dynamics (CFD) methods very often fail to provide an accurate representation of vortices. This is primarily related to the schemes’ numerical dissipation which, if inadequately tuned for the calculation of vortical flows, results in the artificial spreading and diffusion of vortices in numerical simulations. Among other approaches, the Vorticity Confinement (VC) method of J. Steinhoff allows balancing the baseline dissipation within vortices by introducing non-linear anti-dissipation in the discretization of the flow equations, but remains at most first-order accurate. At the same time, remarkable progress has recently been made on the development of high-order numerical methods. These allow reducing the problem of excess dissipation, but the diffusion of vortices remains important for many applications. The present study aims at developing high-order extensions of the VC method to reduce the excess dissipation of vortices, while preserving the accuracy of high-order methods. First, the schemes are analyzed in the case of the linear transport equation, based on time-space coupled and uncoupled formulations. A spectral analysis of nonlinear schemes with VC is performed analytically and numerically, due to their nonlinear character. These schemes exhibit improved dispersive and dissipative properties compared to their linear counterparts at all orders of accuracy. In a second step, third- and fifth-order accurate VC schemes are developed for the compressible Navier-Stokes equations. These remain conservative, rotationally invariant and independent of the baseline scheme, as the original VC2 formulation. Numerical tests validate the increased order of accuracy and the capability of high-order VC extensions to balance dissipation within vortices. Finally, schemes with VC are applied to the calculation of turbulent flows, in an implicit Large Eddy Simulation (ILES) approach. In these applications, numerical schemes with VC exhibit improved resolvability compared to their baseline linear version, while they are capable of producing consistent results even in complex vortical flows.

Dans le secteur aéronautique, la performance des alliages d’aluminium connait une amélioration continue, grâce notamment à l’optimisation des procédés d’élaboration. Dans ce cadre, le travail de recherche vise à prédire le comportement des inclusions dans un bain d’aluminium brassé par induction afin d’améliorer la propreté inclusionnaire des alliages coulés. Un modèle numérique a été développé pour simuler le comportement magnétohydrodynamique du bain d’aluminium dans le creuset suivi par la modélisation du comportement d’une population d’inclusions non-métalliques. Le modèle 2D axisymétrique en régime transitoire s’appuie sur le code de CFD commercial ANSYS Fluent, bien que de nombreuses fonctions utilisateurs aient été introduites pour simuler les phénomènes spécifiques comme l’induction électromagnétique et la résolution des bilans de population. Le modèle MHD résout dans un unique maillage les phénomènes d’induction électromagnétiques, l’écoulement turbulent du bain d’Al, la déformation de la surface libre et les effets de la présence d’une couche de métal oxydée en surface du bain. Une méthode dite de vitesse de glissement (entre les particules et le fluide) a été choisie pour simuler à la fois le transport macroscopique des inclusions dans le bain d’Al et les interactions mésoscopique entre les inclusions (c.à.d. les mécanismes d’agrégation et de fragmentation). Des campagnes expérimentales à l’échelle d’un four de laboratoire et d’une installation industrielle accompagnent le travail numérique pour le valider. Les résultats de modélisation MHD exprimés sous la forme du profil de déformation du bain sont en accord raisonnable avec les mesures faites au laboratoire. Les résultats numériques démontrent également l'effet du frottement induit par la couche d'oxyde sur le profil du bain, ainsi que sur l'écoulement à proximité de la surface du dôme. Pour des conditions opératoires du four industriel en mode de maintien, l’évolution temporelle de la population au sein du bain est calculée. Il apparaît que la séparation magnétique est très intense, particulièrement dans la peau électromagnétique, et est ainsi responsable du transport et de la capture d’une grande fraction de la population d’inclusions à la paroi du four. / With an objective of improving processing and development of aerospace aluminum alloys, the current dissertation presents experimental and numerical tools which help comprehend the behavior of a non-metallic inclusion population in an Al bath stirred by induction. The mechanisms occurring in the metallurgical reactor were separated into two interlinked issues – (i) Magnetohydrodynamics (MHD) of the induction furnace, and (ii) Inclusion population dynamics in the Al bath, which were modeled using the ANSYS Fluent software and in-house User Defined Functions. For a 2D axisymmetric geometry, numerical simulations were performed in a single framework and calculated: (i) the electromagnetic forces using the A-V formulation, (ii) the free surface deformation using the Volume Of Fluid method, (iii) the turbulent stirring of the bath using a RANS-based k-omega model and (iv) the friction force due to the oxide layer by imposing a pseudo-wall condition on the bath free surface. The steady state MHD results and the physical properties of the inclusion population were used as input data for the transient inclusion behavior modeling. A combination of the Drift Concentration Method and the Population Balance Method was developed to respectively model the mean transport of inclusions within the bath at the macroscopic scale and the inclusion interactions (turbulent aggregation and fragmentation) at the mesoscopic scale. The performance of the MHD numerical tool was evaluated by comparing the model results with experimental results at laboratory and industrial scales. The simulation results in the form of the average bath surface profile were found to be consistent with the laboratory measurements. The results also illustrated the impact of the friction due to the oxide layer on the bath surface deformation as well as on the flow near the dome interface. The inclusion behavior simulations were performed for the holding mode operation of an industrial IMF. The deduced removal frequency compared the relative importance of each phenomenon. It was found that the electromagnetic migration, especially in the electromagnetic skin, dominates the inclusion dynamics and is responsible for the capture of a large fraction of the inclusion population.

Les simulations dépendant du temps sont toujours un challenge dans l'industrie, notamment à cause des problèmes posés par les géométries mobiles en termes de CPU et de précision. Cette thèse présente des contributions à certains aspects des simulations en géométrie mobile. Un algorithme de bouger de maillage fondé sur une déformation de maillage sur un grand pas de temps et des changements de connectivité (swaps) est étudié. Une méthode d'élasticité et une méthode d'interpolation directe sont comparées en 3D, démontrant l'efficacité de l'algorithme. Cet algorithme est couplé à un solver ALE, dont les schémas et l'implémentation en 3D sont décrits en détail. Une interpolation linéaire est utilisée pou traiter les swaps. Des cas de validation montrent que les swaps n'influent pas notablement sur la précision de la solution. Plusieurs examples complexes en 3D démontrent la puissance de cette approche, pour des mouvement imposés ou pour des problèmes d'Interaction Fluide-Structure. L'adaptation de maillage anisotrope a démontré son efficacité pour améliorer la précision des calculs stationnaires pour un coût raisonnable. On considère l'extension de ces méthodes aux problèmes instationnaires, en mettant à jour l'algorithme de point fixe précédent grâce à une ananlyse de l'erreur espace-temps fondée sur le modèle de maillage continu. Une parallélisation efficace permet de réaliser des simulations adaptatives instationnaires avec une précision inégalée. Cet algorithme est étendu au cas des géométries mobiles en corrigeant la métrique optimale instationnaire. Finalement, plusieurs exemples 3D de simulations adaptatives en géométries mobiles démontrent l'efficacité de l'approche. / Time dependent simulations are still a challenge for industry, notably due to problems raised by moving boundaries, both in terms of CPU cost and accuracy. This thesis presents contributions to several aspects of simulations with moving meshes. A moving-mesh algorithm based on a large deformation time step and connectivity changes (swaps) is studied. An elasticity method and an Inverse Distance Weighted interpolation method are compared on many 3D examples, demonstrating the efficiency of the algorithm in handling large geometry displacement without remeshing. This algorithm is coupled with an Arbitrary-Lagrangian-Eulerian (ALE) solver, whose schemes and implementation in 3D are described in details. A linear interpolation scheme is used to handle swaps. Validation test cases showed that the use of swaps does not impact notably the accuracy of the solution, while several other complex 3D examples demonstrate the capabilities of the approach both with imposed motion and Fluid-Structure Interaction problems. Metric-based mesh adaptation has proved its efficiency in improving the accuracy of steady simulation at a reasonable cost. We consider the extension of these methods to unsteady problems, updating the previous fixed-point algorithm thanks to a new space-time error analysis based on the continuous mesh model. An efficient p-thread parallelization enables running 3D unsteady adaptative simulations with a new level of accuracy. This algorithm is extended to moving mesh problems, notably by correcting the optimal unsteady metric. Finally several 3D examples of adaptative moving mesh simulations are exhibited, that prove our concept by improving notably the accuracy of the solution for a reasonable time cost.

La simulation numérique de phénomènes physiques complexes requiert généralement l’utilisation d’un maillage. En mécanique des fluides numérique, cela consisteà représenter un objet dans un gros volume de contrôle. Cet objet étant celui dont l’on souhaite simuler le comportement. Usuellement, l’objet et la boîte englobante sont représentés par des maillage de surface linéaires et la zone intermédiaire est remplie par un maillage volumique. L’objectif de cette thèse est de s’intéresser à deux manières différentes de représenter cet objet. La première approche dite immergée consiste à mailler intégralement le volume de contrôle et ensuite à simuler le comportement autour de l’objet sans avoir à mailler explicitement dans le volume ladite géometrie. L’objet étant implicitement pris en compte par le schéma numérique. Le couplage de cette méthode avec de l’adaptation de maillage linéaire est notamment étudié. La deuxième approche dite d’ordre élevé consiste quant à elle consiste à augmenter le degré polynomial du maillage de surface de l’objet. La première étape consiste donc à générer le maillage de surface de degré élevé et ensuite àpropager l’information de degré élevé dans les éléments volumiques environnants si nécessaire. Dans ce cadre-là, il s’agit de s’assurer de la validité de telles modifications et à considérer l’extension des méthodes classiques de modification de maillages linéaires. / The numerical simulation of complex physical phenomenons usually requires a mesh. In Computational Fluid Dynamics, it consists in representing an object inside a huge control volume. This object is then the subject of some physical study. In general, this object and its bounding box are represented by linear surface meshes and the intermediary zone is filled by a volume mesh. The aim of this thesis is to have a look on two different approaches for representing the object. The first approach called embedded method consist in integrally meshing the bounding box volume without explicitly meshing the object in it. In this case, the presence of the object is implicitly simulated by the CFD solver. The coupling of this method with linear mesh adaptation is in particular discussed.The second approach called high-order method consist on the contrary by increasing the polynomial order of the surface mesh of the object. The first step is therefore to generate a suitable high-order mesh and then to propagate the high-order information in the neighboring volume if necessary. In this context, it is mandatory to make sure that such modifications are valid and then the extension of classic mesh modification techniques has to be considered.

Le développement des énergies renouvelables passe par l’exploitation de nouvelles sources d’énergie. La filière hydrolienne, dédiée à la récupération de l’énergie des courants de marée, est proche de l’industrialisation. Cependant, les conditions hydrodynamiques turbulentes des sites hydroliens sont encore mal connues. Cette thèse propose d’examiner à l’échelle locale l’effet des rugosités du fond marin sur la génération de tourbillons hautement énergétiques par la simulation numérique en mécanique des fluides de type méthode de Boltzmann sur réseau. Cette méthode est particulièrement adaptée à la simulation d’écoulements instationnaires sur un domaine de simulation complexe. Dans un premier temps, les phénomènes physiques de détachements tourbillonnaires sur des macro-rugosités canoniques sont décrits. L’appariement de structures tourbillonnaires est mis en évidence dans le processus de formation de tourbillons hautement énergétiques. Dans un deuxième temps, la simulation permet d’observer de tels phénomènes dans le cas d’écoulements environnementaux intégrant une bathymétrie réelle. Ces simulations, validées par rapport à des mesures in situ, mènent à une meilleure compréhension des effets du fond marin sur la turbulence en milieu hydrolien. En particulier, l’importance des failles géologiques dans la génération de turbulence dans la zone d’étude est mise en évidence. / Renewable energy development calls for exploitation of new energy resources. Tidal stream power harvesting is now close to the industrialisation step. Still, turbulent hydrodynamic conditions at tidal sites are not well understood. This thesis aims to investigate the local scale effect of sea bottom roughnesses on energetic vortex generation with computational fluid simulations using the lattice Boltzmann method. This method is highly indicated for unsteady flow simulations of complex domains. First, the physical phenomena involved in vortex emission around canonical macroroughnesses are described. Vortex merging is identified in the generation process of energetic vortices. Then, such physical events are reproduced in the case of environmental flow simulations using a real seabed morphology. These simulations are validated on in situ measured data, and lead to a better understanding of the sea bottom effect on tidal stream site turbulence. They demonstrate the role of geological faults on the local turbulence.