Letteratura scientifica selezionata sul tema "Méthode de Boltzmann sur réseau – Aérodynamique"

La simulation numérique appliquée à la mécanique des fluides est devenue un outil de conception indispensable pour l’industrie aéronautique. Alors que la plupart des simulations industrielles sont réalisées à l’aide d’une approche RANS (Reynolds Averaged Navier-Stokes), celle-ci montre ses limites dès lors qu’il s’agit de caractériser finement des écoulements turbulents instationnaires ou d’étudier des phénomènes aéroacoustiques large-bande. Ainsi, les industriels expriment un besoin grandissant d'outils de simulation haute-fidélité performants. Deux méthodes numériques se montrent particulièrement prometteuses pour la réalisation de telles simulations dans un futur proche : les méthodes Navier-Stokes et la méthode de Boltzmann sur réseau (LBM). Les travaux menés dans le cadre de cette thèse ont contribué à fournir une meilleure compréhension des avantages et des inconvénients respectifs de ces deux méthodes, démontrant que les méthodes de Boltzmann sur réseau et Navier-Stokes se complètent plutôt qu'elles ne se concurrencent. Pour cela, l’étude s’est divisée en deux grandes parties. En premier lieu, une comparaison exhaustive et rigoureuse des méthodes de Boltzmann sur réseau et Navier-Stokes a été réalisée. Différents aspects des méthodes numériques ont été discutés comme leur dissipation et dispersion intrinsèque, leur performance dans un environnement de calcul parallèle (HPC) ainsi que leur capacité à simuler efficacement différents problèmes canoniques de la LES à un niveau de précision donné. Cette étude a permis d’apporter un nouveau regard sur les propriétés des méthodes de Boltzmann sur réseau et Navier-Stokes et de fournir des éléments d’aide à la décision afin d’orienter le choix des ingénieurs vers l’utilisation d’une méthode par rapport à l’autre selon le type d’application visée et le niveau de fidélité requis. Dans un second temps, la possibilité de la mise en place d’un couplage entre les méthodes de Boltzmann sur réseau et Navier-Stokes a été explorée. En effet, de nombreux aspects de la LBM posent encore problèmes ou restent peu efficaces. En particulier le traitement numérique de la zone de proche paroi reste mal défini due à la forme cartésienne des maillages imposée par la méthode. A l'inverse, les approches Navier-Stokes classiques sont particulièrement performantes dans le voisinage de la paroi de par l'utilisation de maillages curvilignes à très grand rapport d'aspect et de méthodes d'intégration temporelle implicites. Ainsi, une méthode numérique hybride innovante a été développée reposant sur un couplage zonal des méthodes de Boltzmann sur réseau et Navier-Stokes, puis étendu au cas de maillages recouvrants (approche Chimère). De nombreuses validations permettent de montrer l’intérêt de cette stratégie. Notamment, cette nouvelle méthode hybride permet de réduire le coût de simulations aéroacoustiques directes en offrant une précision optimale
Computational Fluid Dynamics has become an important design tool for the aeronautical industry. While most industrial simulations are carried out using a RANS (Reynolds Averaged Navier-Stokes) approach, this approach is showing its limitations when it comes to finely characterising unsteady turbulent flows or studying broadband aeroacoustic phenomena. In this context, manufacturers are increasingly looking for high-performance, high-fidelity simulation tools. Two numerical methods are showing particular promise for performing industrial-scale high-fidelity flow simulations in the near future: the Navier-Stokes method and the lattice Boltzmann method (LBM). These two approaches are often presented as competing, but each has its own specific features and requirements. The research carried out as part of this thesis has helped to provide a better understanding of the respective advantages and disadvantages of these two methods, revealing that the lattice Boltzmann and Navier-Stokes methods complement each other rather than compete. The study is divided into two main parts. Firstly, a comprehensive and rigorous comparison of the lattice Boltzmann and Navier-Stokes methods was conducted. The numerical methods were examined in various aspects, such as their intrinsic dissipation and dispersion, their performance in a parallel computing environment (HPC) and their ability to efficiently simulate various canonical LES problems at a given level of accuracy. This study has offered a new perspective on the properties of the lattice Boltzmann and Navier-Stokes methods, providing several decision aids to help the CFD community choose one method over the other based on the type of application and the fidelity level required. Secondly, this PhD explored the possibility of coupling the lattice Boltzmann and Navier-Stokes methods. Indeed, while the LBM offers many benefits, there are still some issues and inefficiencies, especially regarding the numerical treatment of the near-wall zone. The Cartesian shape of the meshes imposed by the method is one of the main reasons for this problem. In contrast, classical Navier-Stokes approaches are particularly effective in the vicinity of the wall thanks to the use of curvilinear meshes with very high aspect ratios and implicit time integration methods. Therefore, an innovative hybrid numerical method was developed based on a zonal coupling of the lattice Boltzmann and Navier-Stokes methods. This approach was then extended to the case of overset meshes (Chimera approach). Numerous validations have demonstrated the value of this strategy. In particular, this new hybrid method makes it possible to reduce the cost of direct aeroacoustic simulations while maintaining optimum accuracy

Cette étude porte sur le développement d'outils d'analyse pour l'identification des évènements aérodynamiques à l'origine des émissions sonores dans les écoulements turbulents. Une application aux écoulements autour d'un véhicule automobile et à l'intérieur d'un habitacle automobile est proposée dans ces travaux. Ces écoulements permettent d'aborder aussi bien le rayonnement direct que le rayonnement des panneaux soumis à une excitation aéroacoustique. Pour ces analyses, des simulations numériques directes basées sur la méthode Boltzmann sur réseau ont été mises en œuvre pour analyser les phénomènes aéroacoustiques présents dans ces configurations d'écoulements tels que l'aérateur d'une voiture, l'écoulement autour d'un rétroviseur, l'écoulement autour d'une voiture (vitrages). Des premières investigations à partir de ces bases de données ont été menées en utilisant des méthodes de causalité permettant de relier les évènements aérodynamiques aux émissions sonores. Ainsi, pour l'analyse du rayonnement direct, la Décomposition Orthogonale aux valeurs Propres Etendue (EPOD) est retenue. L'utilisation de cette méthode a nécessité dans un premier temps de s'assurer de la fiabilité des grandeurs numériques disponibles en terme de convergence statistique. On a ainsi pu mettre en évidence l'importance d'une telle étude préliminaire de faisabilité pour de telles analyses aéroacoustiques. Ensuite, la procédure initiale de l'EPOD est améliorée par la prise en compte d'un temps de retard entre le bruit rayonné et le champ aérodynamique. Différentes variables (champ de vitesse, tenseur de Lighthill, pression) dans la région aérodynamique pour représenter les évènements ''sources'', ont également été testées. Il a ainsi été montré l'intérêt d'utiliser le champ de pression pour de telles analyses. Des évènements aérodynamiques à l'origine des émissions sonores émises à certaines fréquences ont ainsi pu être caractérisés. L'étude du rayonnement acoustique des panneaux est tout d'abord effectuée par un modèle éléments finis d'une plaque sous une excitation déterministe. Cette description déterministe est obtenue à l'aide d'une démarche de synthèse partant d'un modèle statistique décrivant les composantes turbulente et acoustique du champ d'excitation aéroacoustique. La réponse vibroacoustique de la plaque vis-à-vis des champs homogènes d'excitation est étudiée et comparée favorablement aux prédictions théoriques. L'application à un champ d'excitation inhomogène a également été effectuée avec succès. Ensuite, une analyse de contribution des différentes composantes (aérodynamique et acoustique) du champ d'excitation au champ acoustique rayonné est effectuée à l'aide de la fonction de cohérence. Les limitations de cette méthode de causalité qui sont dues à la nature étendue et incohérente des sources ont été démontrées. Enfin, une démarche basée sur la décomposition du champ de pression pariétale par la POD est développée afin de séparer les composantes acoustique et turbulente du champ total d'excitation. Une première application de cette méthode à des données synthétiques a montré l'efficacité de cette nouvelle approche. Cette méthode est finalement appliquée avec succès à un champ de pression pariétale issu de la simulation numérique de l'écoulement autour d'un véhicule réel. Cette dernière étape de ce travail offre de nombreuses perspectives relatives aussi bien à la modélisation des champs de pression pariétale qu'à des analyses des couplages aéroacoustiques de ces champs de pression.

Cette thèse porte sur l'étude numérique des fluctuations extrêmes de la force de traînée exercée par un écoulement turbulent sur un corps immergé.Ce type d'évènement, très rare, est difficile à caractériser par le biais d'un échantillonnage direct, puisqu'il est alors nécessaire de simuler l'écoulement sur des durées extrêmement longues. Cette thèse propose une approche différente, basée sur l'application d'algorithmes d'échantillonnage d'événements rares. L'objectif de ces algorithmes, issus de la physique statistique, est de modifier la statistique d'échantillonnage des trajectoires d'un système dynamique, de manière à favoriser l'occurrence d'événements rares. Si ces techniques ont été appliquées avec succès dans le cas de dynamiques relativement simples, l'intérêt de ces algorithmes n'est à ce jour pas clair pour des dynamiques déterministes extrêmement complexes, comme c'est le cas pour les écoulement turbulents.Cette thèse présente tout d'abord une étude de la dynamique et de la statistique associée aux fluctuations extrêmes de la force de traînée sur un obstacle carré fixe immergé dans un écoulement turbulent à deux dimensions. Ce cadre simplifié permet de simuler la dynamique sur des durées très longues, permettant d'échantillonner un grand nombre de fluctuations dont l'amplitude est assez élevée pour être qualifiée d'extrême.Dans un second temps, l'application de deux algorithmes d’échantillonnage est présentée et discutée.Dans un premier cas, il est illustré qu'une réduction significative du temps de calcul d'extrêmes peut être obtenue. En outre, des difficultés liées à la dynamique de l'écoulement sont mises en lumière, ouvrant la voie au développement de nouveaux algorithmes spécifiques aux écoulements turbulents
This thesis discusses the numerical simulation of extreme fluctuations of the drag force acting on an object immersed in a turbulent medium.Because such fluctuations are rare events, they are particularly difficult to investigate by means of direct sampling. Indeed, such approach requires to simulate the dynamics over extremely long durations.In this work an alternative route is introduced, based on rare events algorithms.The underlying idea of such algorithms is to modify the sampling statistics so as to favour rare trajectories of the dynamical system of interest.These techniques recently led to impressive results for relatively simple dynamics. However, it is not clear yet if such algorithms are useful for complex deterministic dynamics, such as turbulent flows.This thesis focuses on the study of both the dynamics and statistics of extreme fluctuations of the drag experienced by a square cylinder mounted in a two-dimensional channel flow.This simple framework allows for very long simulations of the dynamics, thus leading to the sampling of a large number of events with an amplitude large enough so as they can be considered extreme.Subsequently, the application of two different rare events algorithms is presented and discussed.In the first case, a drastic reduction of the computational cost required to sample configurations resulting in extreme fluctuations is achieved.Furthermore, several difficulties related to the flow dynamics are highlighted, paving the way to novel approaches specifically designed to turbulent flows

Ce travail de thèse s'inscrit dans une problématique visant à étudier numériquement le bruit d'origine aérodynamique généré par les écoulements turbulents autour des véhicules en utilisant la méthode Boltzmann sur Réseau (LBM). Les objectifs de cette thèse sont l'étude des capacités aéroacoustiques de la LBM ainsi que l'élaboration d'un code de calcul tridimensionnel et parallèle.
Dans un premier temps, les élements historiques et théoriques de la LBM sont présentés ainsi que le développement permettant de passer de l'équation de Boltzmann aux équations de Navier-Stokes. La construction des modèles à vitesses discrètes est également décrite. Deux modèles basés sur des opérateurs de collision différents sont présentés : le modèle LBM-BGK et le modèle LBM-MRT. Pour l'étude des capacités aéroacoustiques de la LBM, une analyse de von Neumann est réalisée pour les modèles LBM-BGK et LBM-MRT ainsi que pour l'équation de Boltzmann à vitesse discrète (DVBE). Une comparaison avec les schémas Navier-Stokes d'ordre élevé est alors menée. Pour remédier aux instabilités numériques de la méthode Boltzmann sur Réseau intervenant lors de la propagation dans des directions particulières à M>0.1, des filtres sélectifs sont utilisés et leur effet sur la dissipation est étudié.
Dans un second temps, le code de calcul L-BEAM est présenté. La structure générale et les différentes techniques de calculs sont décrites. Un algorithme de transition de résolution est développé. La modélisation de la turbulence est abordée et le modèle de Meyers-Sagaut est implémenté dans le code. Enfin, des cas tests numériques sont utilisés pour valider le code et la simulation d'un écoulement turbulent complexe est réalisée.

Un nouveau logiciel CFD, LaBS, basé sur la méthode de lattice Boltzmann sur Réseau a été développé dans le cadre d'un projet entre universités et industries. LaBS est utilisé pour la simulation numérique des écoulements thermiques avec un nouveau modèle de frontière immergée pour les conditions limites thermiques. Ce modèle est basé sur la méthode de reconstruction de la fonction de distribution et est évalué pour des conditions limites coincidentes et non-coincidentes avec le maillage, sur le phénomène de diffusion thermique et de convection naturelle.Renault s'intéresse aux situations d'arrêt péage ou de contact coupé, pour lesquelles sont considérés un véhicule roulant à une vitesse soutenue, sur une autoroute par exemple, et qui subit un arrêt ou un ralentissement brutal (avec ou sans contact coupé).Dans ce genre de situation le refroidissement du compartiment moteur qui était assuré par le phénomène de convection forcé durant le roulage laisse place au phénomène de convection naturelle durant les phases de base vitesse ou de vitesse nulle.Le phénomène de convection naturelle est un phénomène lent, qui peut prendre plusieurs minutes à évacuer la chaleur accumulée dans le compartiment moteur. La présence de température élevée pendant une durée trop importante dans le compartiment moteur peut endommager certains composants qui possèdent des seuils de température critique.Pour anticiper ce problème de surchauffe du compartiment moteur, dans lequel un grand nombre de pièces à géométries complexes sont présentes, le phénomène de convection naturelle est étudié avec le nouveau modèle de frontière immergée thermique.%Ce modèle est d'abord testé sur des cas test académique pour validation et est ensuite appliqué au cas d'une voiture réelle.La modélisation des écoulements thermiques avec la méthodes de lattice Boltzmann sur Réseau (LBM) peut-être classée en trois catégories: l'approche multi-vitesse, l'approche hybride et l'approche à deux fonctions de distribution (DDF: Double-Distribution-Function).L'approche multi-vitesse, utilise une équation pour résoudre le champ de vitesse, de densité et de température qui sont résolus avec la LBM. Tandis que l'approche hybride et l'approche DDF utilise un jeux de deux équations, un pour résoudre le champ de vitesse et de densité et l'autre pour résoudre le champ de température.L'approche hybride résout le champ de vitesse et de densité avec la LBM et utilise une méthode de différence finie ou de volume fini pour résoudre le champ de température. L'approche DDF résout quand à elle les deux équations avec la LBM.Le modèle thermique utilisé dans LaBS est basé sur l'approche DDF où les deux équations sont couplées par l'hypothèse de Boussinesq. Le champ de vitesse et de densité est résolu avec un réseau de dix-neuf vitesses discrètes (D3Q19) et champs de température est résolut soit par un réseau à dix-neuf vitesses discrètes (D3Q19) soit par un réseau à sept vitesses discrètes (D3Q7).Le nouveau modèle de frontière immergée décompose la fonction de distribution aux noeuds frontière en sa partie à l'équilibre et hors équilibre. La partie hors équilibre est calculée à partir d'une formulation théorique issus du développement de Chapman-Enskog.La validation du modèle DDF implémenté dans LaBS est faite sur un ensemble de cas test de complexité croissante. Les résultats obtenus avec LaBS sont comparés aux solutions analytiques ou encore à des articles de référence et sont en accord avec les résultats attendus. Ils montrent que qualitativement les résultats sont aussi bons pour le modèle D3Q19/D3Q19 que pour le modèle D3Q19/D3Q7 mais que quantitativement le modèle D3Q19/D3Q19 reste meilleur
A new three-dimensional CFD solver, LaBS, based on the lattice Boltzmann alogorithms has been developed in a framework of university and industry consortium. In this thesis, this solver is used to simulate thermal flows, with a new thermal boundary condition for immersed solid boundary. The new proposed thermal boundary condition is based on the reconstruction method of the distribution function and is evaluated for immersed solid with coincident and non-coincident wall on the case of diffusion and natural convection phenomena.Renault case study, deals with a vehicle moving at constant speed (highway) that suddently slows down and stops (with or without a cut off contact). In such situation the cooling of the engine compartment first driven by forced convection during taxiing stage, abruptly switches to natural convection in low velocity stages. As natural convection is a slow process, it can take several minutes to remove the accumulated heat in the engine compartment. Such duration could be damaging for some components of the engine compartement which do not tolerate high temperature.In order to anticipate overheating of the engine compartment, where a lot of automotive parts with complex geometry are present and to avoid the above mentioned damages, the phenomenon of natural convection is here studied with the new thermal boundary condition.%The new proposed thermal boundary condition is first tested on academic case studies for validation, and then applied to the case of a real car.The modelling of thermal flows with the lattice Boltzmann method (LBM) can be classified into three categories: the multispeed approach, the hybrid approach and the double-distribution-function (DDF) approach. The multispeed approach, uses only one equation to resolve velocity, density and temperature field, which is solved by the LBM. Whereas the hybrid approach and the DDF approach utilize two sets of equations, one to resolve velocity field and density field and another to resolve temperature field. The hybrid approach solves velocity field and density field by the LBM method and the temperature field by finite-different or finite-volume methods. On the other hand the DDF approach solves the two equations with LBM.The thermal model used in the solver LaBS is based on the coupled DDF approach. In this model, the flow field is solved by a D3Q19 velocity model while the temperature field is solved by a D3Q19 or a D3Q7 velocity model. The coupling between the momentum and the energy transport is made by the boussinesq approximation. The new proposed thermal boundary condition decomposes the distribution function at the boundary node into its equilibrium and non-equilibrium part. The non-equilibrium part is calculated from the theoretical solution based on Chapman-Enskog developement.LaBS thermal model based on the coupled DDF approach is evaluated on a set of cases with increasing complexity. The results obtained with LaBS are compared with analytical solutions or with reference articles and are in a good agreement with the results expected. Results show that the model D3Q19/D3Q7 is qualitatively as good as the model D3Q19/D3Q19 but quantitatively the model D3Q19/D3Q19 remains the best

Cette thèse comporte trois parties: étude du schéma de Boltzmann sur réseau, schéma adjoint de Boltzmann sur réseau pour l'identification de paramètres et construction d'une couche parfaitement absorbante pour ce schéma.

La première partie introduit et analyse la méthode.

La deuxième partie décrit une approche variationnelle pour l'assimilation de paramètres relatifs à la méthode du gaz de Boltzmann sur réseau. Une méthode adjointe discrète en temps est développée. L'algorithme est d'abord testé sur un écoulement de type Poiseuille linéaire (problème de Stokes), puis il est appliqué à un problème non linéaire. Des résultats encourageants sont obtenus pour un et deux paramètres inconnus.

Finalement la troisième partie décrit une adaptation des couches absorbantes de Bérenger. Il en résulte un modèle d'automate de Boltzmann à neuf vitesses discrètes. Une analyse des ondes réfléchies est ensuite réalisée entre deux milieux de Boltzmann à une dimension, ce qui permet d'obtenir un équivalent des formules de Fresnel pour les schémas de Boltzmann et de proposer des modifications du schéma à l'interface pour annuler les ondes réfléchies. En deux dimensions, la même analyse d'ondes réfléchies met en évidence l'apparition de modes de Knudsen et des ondes transverses qui rendent l'analyse complexe.

L'anévrisme cérébral est une région fragile de la paroi d'un vaisseau sanguin dans le cerveau, qui peut se rompre et provoquer des saignements importants et des accidents vasculaires cérébraux. La segmentation de l'anévrisme cérébral est une étape primordiale pour l'aide au diagnostic, le traitement et la planification chirurgicale. Malheureusement, la segmentation manuelle prend encore une part importante dans l'angiographie clinique et elle est devenue couteuse en temps de traitement étant donné la gigantesque quantité de données générées par les systèmes d'imagerie médicale. Les méthodes de segmentation automatique d'image constituent un moyen essentiel pour faciliter et accélérer l'examen clinique et pour réduire l'interaction manuelle et la variabilité inter-opérateurs. L'objectif principal de ce travail de thèse est de développer des méthodes automatiques pour la segmentation et la mesure des anévrismes. Le présent travail de thèse est constitué de trois parties principales. La première partie concerne la segmentation des anévrismes géants qui contiennent à la fois la lumière et le thrombus. La méthode consiste d'abord à extraire la lumière et le thrombus en utilisant une procédure en deux étapes, puis à affiner la forme du thrombus à l'aide de la méthode des courbes de niveaux. Dans cette partie, la méthode proposée est également comparée à la segmentation manuelle, démontrant sa bonne précision. La deuxième partie concerne une approche LBM pour la segmentation des vaisseaux dans des images 2D+t et de l'anévrisme cérébral dans les images en 3D. La dernière partie étudie un modèle de segmentation 4D en considérant les images 3D+t comme un hypervolume 4D et en utilisant un réseau LBM D4Q81, dans lequel le temps est considéré de la même manière que les trois autres dimensions pour la définition des directions de mouvement des particules dans la LBM, considérant les données 3D+t comme un hypervolume 4D et en utilisant un réseau LBM D4Q81. Des expériences sont réalisées sur des images synthétiques d'hypercube 4D et d'hypersphere 4D. La valeur de Dice sur l'image de l'hypercube avec et sans bruit montre que la méthode proposée est prometteuse pour la segmentation 4D et le débruitage
Cerebral aneurysm is a fragile area on the wall of a blood vessel in the brain, which can rupture and cause major bleeding and cerebrovascular accident. The segmentation of cerebral aneurysm is a primordial step for diagnosis assistance, treatment and surgery planning. Unfortunately, manual segmentation is still an important part in clinical angiography but has become a burden given the huge amount of data generated by medical imaging systems. Automatic image segmentation techniques provides an essential way to easy and speed up clinical examinations, reduce the amount of manual interaction and lower inter operator variability. The main purpose of this PhD work is to develop automatic methods for cerebral aneurysm segmentation and measurement. The present work consists of three main parts. The first part deals with giant aneurysm segmentation containing lumen and thrombus. The methodology consists of first extracting the lumen and thrombus using a two-step procedure based on the LBM, and then refining the shape of the thrombus using level set technique. In this part the proposed method is also compared with manual segmentation, demonstrating its good segmentation accuracy. The second part concerns a LBM approach to vessel segmentation in 2D+t images and to cerebral aneurysm segmentation in 3D medical images through introducing a LBM D3Q27 model, which allows achieving a good segmentation and high robustness to noise. The last part investigates a true 4D segmentation model by considering the 3D+t data as a 4D hypervolume and using a D4Q81 lattice in LBM where time is considered in the same manner as for other three dimensions for the definition of particle moving directions in the LBM model

Cette thèse comporte trois parties: étude du schéma de Boltzmann sur réseau, schéma adjoint de Boltzmann sur réseau pour l'identification de paramètres et construction d'une couche parfaitement absorbante pour ce schéma. La première partie introduit et analyse la méthode. La deuxième partie décrit une approche variationnelle pour l'assimilation de paramètres relatifs à la méthode du gaz de Boltzmann sur réseau. Une méthode adjointe discrète en temps est développée. L'algorithme est d'abord testé sur un écoulement de type Poiseuille linéaire (problème de Stokes), puis il est appliqué à un problème non linéaire. Des résultats encourageants sont obtenus pour un et deux paramètres inconnus. Finalement la troisième partie décrit une adaptation des couches absorbantes de Bérenger. Il en résulte un modèle d'automate de Boltzmann à neuf vitesses discrètes. Une analyse des ondes réfléchies est ensuite réalisée entre deux milieux de Boltzmann à une dimension, ce qui permet d'obtenir un équivalent des formules de Fresnel pour les schémas de Boltzmann et de proposer des modifications du schéma à l'interface pour annuler les ondes réfléchies. En deux dimensions, la même analyse d'ondes réfléchies met en évidence l'apparition de modes de Knudsen et des ondes transverses qui rendent l'analyse complexe
This thesis is composed of three parts. Firstly a study of Lattice Boltzmann scheme (LBE) is performed. Then Adjoint Lattice Boltzmann scheme (ALBE) is introduced for parameters identification. Finally a new Lattice Boltzmann scheme (BRB) is proposed to modelise B\'erenger's Perfectly Matched Layer (PML) method. The first part introduces and analyzes the LBE method. The second part describes a variational approach for parameters identification adapted to LBE. A time discrete adjoint method is developed. At first the ALBE method is applied to Stokes' problem and then to a nonlinear problem. Good results have been obtained in the cases of one or tow unknown parameters. Finally the third part describes an adaptation of PML for LBE. The LB scheme is obtained with 9 discretes velocities. An analysis of reflected waves between two one dimensional LB media is performed. It provides us an equivalent for the Fresnel formula for LBE interface. That gives us same ideas to modify the LB scheme at the interface to vanish reflected waves. In the two dimensional case, the same analysis of reflected waves shows the existence of Knudsen modes and transverse waves, which make the analysis more difficult

Dans un contexte mondial où l’accès à l’énergie est un problème de premier plan, l’exploitation des courants de marée avec des hydroliennes revête un intérêt certain. Les écoulements dans les zones à fort potentiel énergétique propices à l’installation d’hydroliennes sont souvent fortement turbulents. Or la turbulence ambiante impacte fortement l’hydrodynamique avoisinante et le fonctionnement de la turbine. Une prédiction fine de la turbulence et du sillage est fondamentale pour l'optimisation d'une ferme d'hydroliennes. Un modèle de simulation de l'écoulement autour de la turbine doit donc être précis et tenir compte de la turbulence ambiante. Un outil basé sur la méthode de Boltzmann sur réseau (LBM) est utilisé à ces fins, en association avec une approche de simulation des grandes échelles (LES). La LBM est une méthode instationnaire de modélisation d’écoulement fluide. Une méthode de génération de turbulence synthétique est implémentée afin de prendre en compte la turbulence ambiante des sites hydroliens. Les géométries complexes, potentiellement en mouvement, sont modélisées avec la méthode des frontières immergées (IBM). La mise en place d’un modèle de paroi est réalisée afin de réduire le cout en calcul du modèle. Ces outils sont ensuite utilisés pour modéliser en LBM-LES une hydrolienne dans un environnement turbulent. Les calculs, réalisés à deux taux de turbulence différents, sont comparés avec des résultats expérimentaux et des résultats NS-LES. Les modélisations LBM-LES sont ensuite utilisées pour analyser le sillage de l'hydrolienne. Il est notamment observé qu'un faible taux de turbulence impacte de manière significative la propagation des tourbillons de bout de pale
In a global context where access to energy is a major problem, the exploitation of tidal currents with tidal turbines is of particular interest. Flows in areas with high energy potential suitable for the installation of tidal turbines are often highly turbulent. However, the ambient turbulence has a strong impact on the surrounding hydrodynamics and the turbine operation. A precise prediction of turbulence and wake is fundamental to the optimization of a tidal farm. A numerical model of the flow around the turbine must therefore be accurate and take into account the ambient turbulence. A tool based on the Lattice Boltzmann Method (LBM) is used for this purpose, in combination with a Large Eddy Simulation (LES) approach. The LBM is an unsteady method for modelling fluid flows. A synthetic turbulence method is implemented to take into account the ambient turbulence of tidal sites. Complex geometries, potentially in motion, are modelled using the Immersed Boundary Method (IBM). The implementation of a wall model is carried out in order to reduce the cost of the simulations. These tools are then used to model a turbine in a turbulent environment. The calculations, performed at two different turbulence rates, are compared with experimental and NS-LES results. The LBM-LES models are then used to analyze the wake of the turbine. In particular, it is observed that a low turbulence rate has a significant impact on the propagation of tip-vortices

Ce travail de thèse traite de la modélisation et la simulation des systèmes complexes et s’inscrit dans la thématique du contrôle et la gestion des ressources en eau. Nous proposons une approche originale basée sur les modèles de Boltzmann sur réseau pour représenter les écoulements au sein des canaux d’irrigation, habituellement décrits par les équations non linéaires de Saint-Venant. Nous avons adapté un modèle bi-fluide et étudié les conditions aux bords qui permettent de reproduire la géométrie d’un canal à surface libre. Des méthodes de détermination des grandeurs hydrauliques d’intérêt ont été développées. Nous nous sommes intéressés au comportement des vannes de fond en régime noyé, et nous démontrons que le modèle permet de décrire celui-ci spontanément dans différentes situations de fonctionnement. Des validations ont été effectuées par le biais d’une comparaison avec des expérimentations menées sur un micro-canal expérimental. Nous avons également introduit dans le modèle des phénomènes de sédimentation et étudié l’influence de la présence d’un dépôt de sédiments sur l’écoulement. Là aussi, les résultats de simulation et d’expérimentation convergent
This PhD work considers the general problem of modelling and simulation of complex systems and deals with the domain of control and management of water resources. We propose here an original approach based on Lattice Boltzmann models (LB) for modelling free surface flows in irrigation canals, usually described with the non-linear shallow water equations. We adapted a bi-fluid model and studied the boundary conditions which allow to reproduce the geometry of a free-surface irrigation canal. Methods for estimating the desired hydraulic quantities were developed. We studied the behavior of submerged underflow gates, and we show that the model is able to spontaneously and correctly describe how the gates function in quite different situations. Validations were realized by comparing results from simulations and experimentations performed on a laboratory micro-canal facility. We also introduced sedimentation phenomena in the model and studied the influence of a sedimentation deposit on the flow. Comparisons between experimental and simulation results were also performed and converged