Academic literature on the topic 'Méthodes de Boltzmann sur réseaux'

Les méthodes de Boltzmann sur réseaux (LBM) ont été appliquées avec beaucoup de succès aux écoulements hydrodynamiques en milieux poreux. Cependant, la limitation de ces méthodes aux écoulements hydrodynamiques et isothermes, les rendent insuffisantes pour simuler des écoulements de gaz dans des milieux micro-poreux. Dans ce cas, il est en effet fréquent que le libre parcours moyen des molécules du gaz, soit du même ordre de grandeur que la taille des pores dans lesquels il s’écoule. De tels écoulements ne seront alors plus en régime hydrodynamique, mais dans des régimes qualifiés de glissement et de transitionnel ; régimes pour lesquels les LBM standards ne sont plus valides. D’autre part, le caractère isotherme des LBM les rendent inutilisables, par exemple dans le cas où le gaz subit une détente à travers le milieu. Il est nécessaire, pour décrire de tels écoulements et phénomènes, de se placer au niveau cinétique. La démarche proposée repose sur la décomposition de la fonction de distribution sur la base des polynômes d’Hermite et l’emploi de la quadrature de Gauss-Hermite associée à cette projection. L’aspect systématique de ce développement amène naturellement à considérer divers ordres d’approximation de l’équation de Boltzmann-BGK sous diverses quadratures. Il résulte alors de ces différentes approximations toute une famille de discrétisations de l’équation de Boltzmann-BGK, dont les LBM classiques ne sont qu’un membre. La détermination de l’approximation la plus adaptée est réalisée par analyse systématique des résultats obtenus aux différents ordres d’approximation. Ces méthodes sont testées avec succès dans des cas modèles<br>The lattice Boltzmann method (LBM) have been applied very successfully to hydrodynamic flows in porous media. However, the limitation of these methods to isothermal and hydrodynamic flows, make them inadequate to simulate gas flows in micro-porous media. Indeed, in these conditions, the mean free path of the molecules could be of the same magnitude order as the pore size in which gas flows. Such flows will not be in hydrodynamic regime, but in regimes qualified of, slip or transitional ; for which the LBM are no longer valid. On the other hand, the isothermal character of LBM make them unusable, for example, in the case where the gas undergoes expansion through the media. It is then necessary, to take the kinetic point of view to describe such flows and phenomena. The proposed approach is based on the decomposition of the distribution function on the Hermite polynomials basis and the use of Gauss-Hermite quadrature associated with this projection. The systematic nature of this development naturally leads to consider different order of approximation of the Boltzmann-BGK equation in various quadratures. It then follows from these various approximations, a family of discretizations of the Boltzmann-BGK equation, whose classical LBM are a member. Determining the most suitable approximation is achieved by systematic analysis of the results obtained with different approximation orders. These methods are successfully tested in model cases

L’objectif de ce travail est d’étudier les capacités de la méthode de Boltzmann sur réseau (LBM) dans un cadre numériquement contraignant : celui de la simulation aéroacoustique en maillage non-uniforme, à très haut nombre de Reynolds et à nombre de Mach non négligeable (Ma &gt; 0.1), appliquée à l’automobile. La problématique industrielle est celle du calcul du bruit intérieur d’origine aérodynamique, dont le calcul du champ de pression pariétal instationnaire sur le vitrage conducteur est la première étape décisive. Il a été constaté qu’un manque de précision sur la faible part acoustique du champ de pression total sur le vitrage, provenant très probablement d’erreurs au niveau des transitions de résolution du maillage, était la cause d’une surestimation du bruit intérieur. Nous présentons d’abord une construction cohérente et unifiée de la méthode de Boltzmann sur réseau à partir de l’équation de Boltzmann, dans un cadre athermal faiblement compressible. Nous étudions ensuite en détail les propriétés aéroacoustiques de la LBM, en parcourant toutes les grandes familles d’opérateurs de collision de la littérature. Une variante de modèle à temps de relaxation multiples, utilisable pour l’aéroacoustique, est présentée et testée. Un modèle alternatif simplifié de filtrage sélectif, rapide et compact, est développé et validé. La problématique des maillages non-uniformes est abordée. Un recensement exhaustif des études LBM menées dans ce cadre dans la littérature montre qu’aucune ne correspond à nos contraintes. Des algorithmes alternatifs aux transitions sont développés. Enfin, des applications industrielles sont réalisées à l’aide des modèles développés dans le mémoire<br>The main goal of this work is to study the capacities of the Lattice Boltzmann Method in a constrained numerical framework : that of numerical simulation in automotive aeroacoustics with non-uniform meshes, at high Reynolds number and non egligible Mach number (Ma &gt; 0.1). The industrial problem is the computation of the interior aerodynamic noise, which includes as its first decisive step the computation of the unsteady wall pressure field on the car windows. It was observed that a lack of precision on the weak acoustic part of the total pressure field on the driver-side window, which is most probably due to errors at mesh refinement interfaces, caused an overestimation of the interior noise. We first present a coherent and unified construction of the Lattice BoltzmannMethod from the Boltzmann equation, in an athermal weakly compressible framework. Then, we study in details the aeroacoustic properties of the LBM by reviewingall the main families of collisional operators that exist in the literature. A variant of multiple relaxation time operator that can be used for aeroacoustics is presented and tested. A simplified alternative selective filter, fast and compact, is developped and numerically validated. The problem of non-uniform meshes is discussed. An exhaustive review of the LBM studies that have been carried out within that framework shows that none of them corresponds to our constraints. Alternative transition nodes algorithms are developed. Finally, all the developed models of this work are applied to industrial cases

Malgré l'efficacité informatique et la faible dissipation numérique de la méthode de Boltzmann sur réseau (LBM) classique reposant sur un algorithme de propagation-collision, cette méthode est limitée aux maillages cartésiens uniformes. L'adaptation de l'étape de discrétisation à différentes échelles de la mécanique des fluides est généralement réalisée par des schémas LBM à échelles multiples, dans lesquels le domaine de calcul est décomposé en plusieurs sous-domaines uniformes avec différentes résolutions spatiales et temporelles. Pour des raisons de connectivité, le facteur de résolution des sous-domaines adjacents doit être un multiple de deux, introduisant un changement abrupt des échelles spatio-temporelles aux interfaces. Cette spécificité peut déclencher des instabilités numériques et produire des sources de bruit parasite rendant l'exploitation de simulations à finalités aéroacoustiques impossible. Dans la présente thèse, nous avons d'abord élucidé le sujet du raffinement de maillage dans la LBM classique en soulignant les défis et les sources potentielles d'erreur. Par la suite, une méthode de Boltzmann sur réseau hybride (HLBM) est proposée, combinant l'algorithme de propagation-collision avec un algorithme de flux au sens eulérien obtenu à partir d'une discrétisation en volumes finis des équations de Boltzmann à vitesse discrète. La HLBM combine à la fois les avantages de la LBM classique et une flexibilité géométrique accrue. La HLBM permet d'utiliser des maillages cartésiens non-uniformes. La validation de la méthode hybride sur des cas tests 2D à finalité aéroacoustique montre qu'une telle approche constitue une alternative viable aux schémas Boltzmann sur réseau à échelles multiples, permettant de réaliser des raffinements locaux en H. Enfin, un couplage original, basé sur l'algorithme de propagation-collision et une formulation isotherme des équations de Navier-Stokes en volumes finis, est proposé. Une telle tentative présente l'avantage de réduire le nombre d'équations du solveur volumes finis tout en augmentant la stabilité numérique de celui-ci, en raison d'une condition CFL plus favorable. Les deux solveurs sont couplés dans l'espace des moments, où la solution macroscopique du solveur Navier-Stokes est injectée dans l'algorithme de propagation-collision à l'aide de la collision des moments centrés. La faisabilité d'un tel couplage est démontrée sur des cas tests 2D, et les résultas obtenus sont comparés avec la HLBM<br>Despite the inherent efficiency and low dissipative behaviour of the standard lattice Boltzmann method (LBM) relying on a two step stream and collide algorithm, a major drawback of this approach is the restriction to uniform Cartesian grids. The adaptation of the discretization step to varying fluid dynamic scales is usually achieved by multi-scale lattice Boltzmann schemes, in which the computational domain is decomposed into multiple uniform subdomains with different spatial resolutions. For the sake of connectivity, the resolution factor of adjacent subdomains has to be a multiple of two, introducing an abrupt change of the space-time discretization step at the interface that is prone to trigger instabilites and generate spurious noise sources that contaminate the expected physical pressure signal. In the present PhD thesis, we first elucidate the subject of mesh refinement in the standard lattice Boltzmann method and point out challenges and potential sources of error. Subsequently, we propose a novel hybrid lattice Boltzmann method (HLBM) that combines the stream and collide algorithm with an Eulerian flux-balance algorithm that is obtained from a finite-volume discretization of the discrete velocity Boltzmann equations. The interest of a hybrid lattice Boltzmann method is the pairing of efficiency and low numerical dissipation with an increase in geometrical flexibility. The HLBM allows for non-uniform grids. In the scope of 2D periodic test cases, it is shown that such an approach constitutes a valuable alternative to multi-scale lattice Boltzmann schemes by allowing local mesh refinement of type H. The HLBM properly resolves aerodynamics and aeroacoustics in the interface regions. A further part of the presented work examines the coupling of the stream and collide algorithm with a finite-volume formulation of the isothermal Navier-Stokes equations. Such an attempt bears the advantages that the number of equations of the finite-volume solver is reduced. In addition, the stability is increased due to a more favorable CFL condition. A major difference to the pairing of two kinetic schemes is the coupling in moment space. Here, a novel technique is presented to inject the macroscopic solution of the Navier-Stokes solver into the stream and collide algorithm using a central moment collision. First results on 2D tests cases show that such an algorithm is stable and feasible. Numerical results are compared with those of the previous HLBM

L'objectif de cette thèse est de développer et d'analyser des techniques numériques basées sur la méthode e Boltzmann sur réseau (LBM) pour résoudre des systèmes non linéaires de type advection-réaction-diffusion provenant de la physique et de la biologie. Avec la LBM, des problèmes portant sur des quantités moyennées densité, potentiel, vitesse, etc) sont exprimés à l'échelle particulaire. Nous approchons la solution de l'équation e Boltzmann relative au comportement d'un champs de particules puis nous recomposons les quantités moyennées solutions des équations traitées. Dans un premier temps, nous développons un cadre général approprié permettant de traiter plusieurs types de systèmes non linéaires (paraboliques, elliptiques, ou couplées ' variables réelles ou complexes), avec des applications à des modèles tels que Burger-Fisher, écoulement de fluides en milieu poreux, Helmoltz, Patlar-Keller-Segel, ou encore Schrodinger. Pour chaque problème, nous analysons le comportement asymptotique de la méthode, quand le nombre de Knudsen tend vers zéro (par le développement de Chapman-Enskog) et nous effectuons l'analyse numérique de la convergence et de la stabilité de la méthode. Dans un deuxième temps, nous nous intéressons à un problème réaliste d'électrophysiologie cardio-vasculaire. Nous adaptons la méthode LBM développée pour approcher les solutions d'un système de type bidomaine permettant de simuler le comportement de potentiels électriques et les interactions ioniques ans la région du myocarde. L'étude et la modélisation d'un tel type de problème est un enjeu sanitaire majeur ans le traitement des pathologies liées par exemple à l'arythmie cardiaque. Notre but étant d'obtenir des comportements réalistes, nous introduisons au sein de ce système bidomaine des opérateurs de retard afin de tenir compte des temps de retard dans les transmissions de signaux. Une fois l'existence et l'unicité de la solution démontrées, nous proposons une série de simulations avec des paramètres physiques et biologiques réalistes afin de valider la méthode proposée<br>In this thesis, we develop and analyze numerical techniques based on the lattice Boltzmann method LBM) for solving systems of nonlinear advection-diffusion-reaction equations from physics and biology. Wi BM, problems relating to averaged quantities (density, potential, velocities, etc.) are expressed at the particle scale. We approach the solution of Boltzmann equation relating to the behavior of a particle field and then we recompose the averaged quantities solutions of treated systems. In the first part, we develop an appropriate general framework to deal with several types of non-linear systems (parabolic, elliptic, or coupled, with real or complex variables), with applications to models such as Burger-Fisher, fluid flow in a porous medium, Helmoltz, Patlar-Keller-Segel, or Schrodinger. For each problem, we analyze the asymptotic behavior of the method, when the number of Knudsen tends to zero (by the development of Chapman-Enskog) and we perform the numerical analysis of convergence and stability of the method. In the second part, we have taken an interest in a realistic problem of cardio-vascular electrophysiology. We adapt the developed LBM method to approach e solutions of a bidomain type system for simulating the behavior of electrical potentials and ionic interactions in myocardial region. The study and modeling of this type of problem is a major health issue in the treatment of pathologies related, for example, to cardiac arrhythmia. Since our goal is to obtain realistic behaviors, we introduce time-delay operators into this coupled system in order to take into account delay in signal transmissions. Once the existence and uniqueness of solution have been demonstrated, we propose a series of simulations with realistic physical and biological parameters to validate the proposed method

Les méthodes de Boltzmann sur réseaux (LBM) ont été appliquées avec beaucoup de succès aux écoulements hydrodynamiques en milieux poreux. Cependant, la limitation de ces méthodes aux écoulements hydrodynamiques et isothermes, les rendent insuffisantes pour simuler des écoulements de gaz dans des milieux micro-poreux. Dans ce cas, il est en effet fréquent que le libre parcours moyen des molécules du gaz, soit du même ordre de grandeur que la taille des pores dans lesquels il s'écoule. De tels écoulements ne seront alors plus en régime hydrodynamique, mais dans des régimes qualifiés de glissement et de transitionnel ; régimes pour lesquels les LBM standards ne sont plus valides. D'autre part, le caractère isotherme des LBM les rendent inutilisables, par exemple dans le cas où le gaz subit une détente à travers le milieu. Il est nécessaire, pour décrire de tels écoulements et phénomènes, de se placer au niveau cinétique. La démarche proposée repose sur la décomposition de la fonction de distribution sur la base des polynômes d'Hermite et l'emploi de la quadrature de Gauss-Hermite associée à cette projection. L'aspect systématique de ce développement amène naturellement à considérer divers ordres d'approximation de l'équation de Boltzmann-BGK sous diverses quadratures. Il résulte alors de ces différentes approximations toute une famille de discrétisations de l'équation de Boltzmann-BGK, dont les LBM classiques ne sont qu'un membre. La détermination de l'approximation la plus adaptée est réalisée par analyse systématique des résultats obtenus aux différents ordres d'approximation. Ces méthodes sont testées avec succès dans des cas modèles.

Cette thèse introduit et étudie une nouvelle classe de schémas de Boltzmann sur réseau appelés schémas à vitesse relative. Les schémas de Boltzmann sur réseau visent à approcher des problèmes de nature macroscopique en mimant la dynamique microscopique d’équations cinétiques du type Boltzmann. L’algorithme calcule des distributions de particules évoluant au travers de deux phases de transport et de relaxation, les particules se déplaçant en les noeuds d’un réseau cartésien en espace. Les schémas de Boltzmann à plusieurs temps de relaxation (ou schéma MRT de d’Humières), dont la relaxation im- plique un ensemble de moments combinaison linéaire polynomiale des distributions, constituent le cadre initial de la thèse. Les schémas à vitesse relative sont une extension de ces schémas de d’Humières. Ils sont inspirés du schéma cascade de Geier apportant davantage de stabilité que les schémas de d’Hu- mières pour des régimes peu visqueux des équations de Navier-Stokes. La différence avec ces schémas se situe au niveau de la relaxation : elle utilise un ensemble de moments relatifs à un paramètre champ de vitesse fonction du temps et de l’espace. Cette différence se matérialise par une matrice de tran- sition des moments fixes (les schémas de d’Humières correspondent à un paramètre champ de vitesse nul) aux moments mobiles. La structure algébrique de cette matrice est étudiée. Le schéma cascade est ensuite traduit comme un schéma à vitesse relative pour un nouvel ensemble de polynômes définissant les moments. L’étude de la consistance des schémas à vitesse relative par la méthode des équations équivalentes est un point central de la thèse. Les équations limites pour un nombre arbitraire de dimen- sions et de vitesses sont dérivées et illustrées sur des exemples tels que le D2Q9 pour les équations de Navier-Stokes. Ces équations équivalentes sont également un outil pour prédire la stabilité des schémas grâce à l’analyse des termes de diffusion et dispersion. La dernière partie traite de la stabilité suivant le choix du paramètre champ de vitesse. Nous sommes particulièrement intéressés en les deux choix de paramètre nul (d’Humières) et la vitesse du fluide (cascade). Le schéma D2Q9 pour les équations de Navier-Stokes est étudié numériquement par une méthode de Von Neumann puis appuyé sur des cas tests non linéaires. La stabilité des schémas relatifs à la vitesse du fluide est dépendante du choix des polynômes définissant les moments. L’amélioration la plus notable se produit si les polynômes du schéma cascade sont choisis. Nous étudions enfin les stabilités théorique et numérique d’un schéma bidimensionnel minimal. Le contexte physique est la simulation d’une équation d’advection diffusion linéaire. Le choix de la vitesse d’advection comme paramètre champ de vitesse annule certains termes de dispersion des équations équivalentes contrairement aux schémas de d’Humières. Ceci se traduit par un meilleur comportement en termes de stabilité pour de grandes vitesses, appuyé théoriquement à l’aide d’une notion de stabilité à poids<br>In this PhD thesis, a new class of lattice Boltzmann schemes called relative velocity schemes is introduced and studied. The purpose of lattice Boltzmann schemes is to approximate problems of macroscopic nature using the microscopic dynamic of Boltzmann type kinetic equations. They compute particle distributions through two phases of transport and relaxation, the particles moving on the nodes of a cartesian lattice. The multiple relaxation times schemes---MRT of d'Humières---, whose relaxation uses a set of moments, linear combinations of the particle distributions, constitutes the initial framework of the thesis. The relative velocity schemes extend the MRT d'Humières schemes. They originate from the cascaded automaton of Geier which provides more stability for the low viscosity regime of the Navier-Stokes equations. Their difference with the d'Humières schemes is carried by the relaxation : a set of moments relative to a velocity field parameter function of space and time is used. This difference is represented by a shifting matrix sending the fixed moments---The d'Humières schemes are associated with a zero velocity field parameter---On the relative moments. The algebraic structure of this matrix is studied. The cascaded automaton is then interpreted as a relative velocity scheme for a new set of polynomials defining the moments. The consistency study of the relative velocity schemes with the equivalent equations method is a keypoint of the thesis. These equations are derived for an arbitrary number of dimensions and velocities. They are then illustrated on examples like the D2Q9 scheme for the Navier-Stokes equations. These equivalent equations are also a tool to predict the stability behaviour of the schemes by analysing their diffusion and dispersion terms. In a last part, the stability according to the velocity field parameter is studied. Two cases especially interest us : a parameter equal to zero---D'Humières schemes---And equal to the fluid velocity---Cascaded automaton. The D2Q9 scheme for the Navier-Stokes equations is numerically studied with a linear Von Neumann analysis and some non linear test cases. The stability of the relative velocity schemes depends on the choice of the polynomials defining the moments. The most important improvement occurs if the polynomials of the cascaded automaton are chosen. We finally study the theoretical and numerical stability of a minimal bidimensional scheme for a linear advection diffusion equation. If the velocity field parameter is chosen equal to the advection velocity, some dispersion terms of the equivalent equations vanish unlike the d'Humières scheme. This implies a better stability behaviour for high velocities, characterized thanks to theoretical weighted stability notion

La connaissance de la mobilité est un enjeu majeur pour les autorités organisatrices de mobilité et l'aménagement urbain. En raison du manque de définition formelle de la mobilité humaine, l'expression "mobilité des personnes" sera utilisée dans cette ouvrage. Ce sujet sera introduit par une description de l'écosystème en considérant ces acteurs et ces applications.La création d'un modèle d'apprentissage a des prérequis: la compréhension des typologies des ensembles de données disponibles, leurs forces et leurs faiblesses. Cet état de l'art de la connaissance de la mobilité passe par le modèle à quatre étapes qui existe et est utilisé depuis 1970 pour finir sur le renouvellement des méthodologies de ces dernières années.Nos modélisations de la mobilité des personnes sont ensuite présentées. Leur point commun est la mise en avant de l'individu contrairement aux approches classiques qui prennent comme référence la localité. Les modèles que nous proposons s'appuient sur le fait que la prise de décision des individus se fait selon leur perception de l'environnement.Cet ouvrage fini sur l'étude des méthodes d'apprentissage profond des machines de Boltzmann restreintes. Après un état de l'art de cette famille de modèles, nous recherchons des stratégies pour rendre ces modèles viables dans le monde applicatif. Ce dernier chapitre est notre contribution théorique principale, par l'amélioration de la robustesse et la performance de ces modèles<br>Knowledge of mobility is a major challenge for authorities mobility organisers and urban planning. Due to the lack of formal definition of human mobility, the term "people's mobility" will be used in this book. This topic will be introduced by a description of the ecosystem by considering these actors and applications.The creation of a learning model has prerequisites: an understanding of the typologies of the available data sets, their strengths and weaknesses. This state of the art in mobility knowledge is based on the four-step model that has existed and been used since 1970, ending with the renewal of the methodologies of recent years.Our models of people's mobility are then presented. Their common point is the emphasis on the individual, unlike traditional approaches that take the locality as a reference. The models we propose are based on the fact that the intake of individuals' decisions is based on their perception of the environment.This finished book on the study of the deep learning methods of Boltzmann machines restricted. After a state of the art of this family of models, we are looking for strategies to make these models viable in the application world. This last chapter is our contribution main theoretical, by improving robustness and performance of these models

Ces dernières années ont été marquées par l'émergence d'un grand nombre de base données dans de nombreux domaines comme la médecine par exemple. La création de ces bases données a ouvert la voie à de nouvelles applications. Les propriétés des données sont parfois complexes (non linéarité, dynamique, grande dimension ou encore absence d'étiquette) et nécessite des modèles d'apprentissage performants. Parmi les modèles d'apprentissage existant, les réseaux de neurones artificiels ont connu un large succès ces dernières décennies. Le succès de ces modèles repose sur la non linéarité des neurones, l'utilisation de variables latentes et leur grande flexibilité leur permettant de s'adapter à de nombreux problèmes. Les machines de Boltzmann présentées dans cette thèse sont une famille de réseaux de neurones non supervisés. Introduite par Hinton dans les années 80, cette famille de modèle a connu un grand intérêt dans le début du 21e siècle et de nouvelles extensions sont proposées régulièrement.Cette thèse est découpée en deux parties. Une partie exploratoire sur la famille des machines de Boltzmann et une partie applicative. L'application étudiée est l'apprentissage non supervisé des signaux électroencéphalogramme intracrânien chez les rats Parkinsonien pour le contrôle des symptômes de la maladie de Parkinson.Les machines de Boltzmann ont donné naissance aux réseaux de diffusion. Il s'agit de modèles non supervisés qui reposent sur l'apprentissage d'une équation différentielle stochastique pour des données dynamiques et stochastiques. Ce réseau fait l'objet d'un développement particulier dans cette thèse et un nouvel algorithme d'apprentissage est proposé. Son utilisation est ensuite testée sur des données jouet ainsi que sur des données réelles<br>The past recent years have been marked by the emergence of a large amount of database in many fields like health. The creation of many databases paves the way to new applications. Properties of data are sometimes complex (non linearity, dynamic, high dimensions) and require to perform machine learning models. Belong existing machine learning models, artificial neural network got a large success since the last decades. The success of these models lies on the non linearity behavior of neurons, the use of latent units and the flexibility of these models to adapt to many different problems. Boltzmann machines presented in this thesis are a family of generative neural networks. Introduced by Hinton in the 80's, this family have got a large interest at the beginning of the 21st century and new extensions are regularly proposed.This thesis is divided into two parts. A first part exploring Boltzmann machines and their applications. In this thesis the unsupervised learning of intracranial electroencephalogram signals on rats with Parkinson's disease for the control of the symptoms is studied.Boltzmann machines gave birth to Diffusion networks which are also generative model based on the learning of a stochastic differential equation for dynamic and stochastic data. This model is studied again in this thesis and a new training algorithm is proposed. Its use is tested on toy data as well as on real database

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.<br />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é.<br />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.

L'objectif de ce travail de recherche est d'étudier numériquement le transport et le dépôt de particules dans des milieux poreux à l'échelle des pores.Premièrement, un couplage entre la méthode de Boltzmann sur réseau (LBM) et la méthode des éléments discrets (DEM) est réalisé et utilisé pour simuler l'écoulement d'un fluide chargé en particules. La LBM est utilisée pour décrire l'écoulement du fluide autour des fibres tandis que la DEM est utilisée pour traiter la dynamique des particules. Ce couplage est bidirectionnel dans le sens où le mouvement des particules affecte le flux de fluide et réciproquement. Ce modèle nous a permis de prédire l'efficacité de capture et la chute de pression à l'étape initiale du processus de filtration. Le facteur de qualité est également calculé pour déterminer la qualité de filtration.Ensuite, on se focalise sur l'étude de l'efficacité de la capture de fibres de formes de section transversale différentes (circulaire, losange et carrée). Les résultats issus de nos simulations du processus de filtration de la fibre circulaire concordent bien avec les corrélations empiriques. L'impaction des particules sur la face avant de la fibre de forme carrée est plus importante que dans les cas de fibre de formes circulaire et losange. Cependant, en raison d'une chute de pression plus faible, la fibre de section losange présente une meilleure qualité de filtration. Ensuite, les variations du facteur de qualité dues à l'angle d'orientation et au rapport d'aspect des fibres ont été étudiées numériquement pour la forme rectangulaire. Pour chaque cas, on a déterminé la valeur optimale de la zone au vent pour laquelle le facteur de qualité est maximal. La comparaison des valeurs du facteur de qualité obtenues pour les différentes formes de fibre monte une meilleure performance pour la fibre de section carrée orientée avec un angle de π/4.Enfin, l'influence de l'arrangement des fibres sur la qualité de la filtration est analysée en considérant la configuration en quinconce pour les différentes formes. Les simulations conduites pour différentes tailles de particules et différentes valeurs de la densité (particule/air) montent que la fibre de section losange est plus performante en termes de facteur de qualité pour les particules de grande taille et pour les valeurs de densité élevée. La présente étude fournit des pistes pour optimiser le processus de filtration et prédire la qualité de filtration<br>The objective of the present research was to numerically investigate the transport and deposition of particles in porous media at the pore scale. Firstly, a developed coupled lattice Boltzmann method (LBM) and discrete element method (DEM) is used to simulate the fluid-particle flow. LBM is employed to describe the fluid flow around fibers whereas DEM is used to deal with the particle dynamics. The corresponding method is two-way coupling in the sense that particle motion affects the fluid flow and reciprocally. It allowed us to predict the capture efficiency and pressure drop at the initial stage of filtration process. The quality factor is also calculated for determining the filtration performance. Secondly, we focus on the study the capture efficiency of single fiber with circular, diamond and square cross-section, respectively. The results of LBM-DEM for filtration process of single circular fiber agree well with the empirical correlation. The impaction of particles on the front side of square-shaped fiber is more favorable than those on circular and diamond cases. However, diamond fiber exhibits a good filtration performance. Then the variations of quality factor due to the different orientation angle and aspect ratio of rectangular fiber were studied using LBM-DEM. For each case, we have found the optimal value of the windward area to which corresponds a maximum value of the quality factor. The comparison of the performance of the different forms of fibers shows that the largest quality factor is obtained for square fiber oriented with angle π/4.Finally, the influence of the arrangement of fiber on filtration performance is analyzed by considering the staggered configuration. Simulations conducted for several particle size and density show that the diamond with staggered array performs better for large particles and high particle-to-fluid density ratio in terms of quality factor. The present study provide an insight to optimize the filtration process and predict filtration performance