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.

La propreté inclusionnaire reste un enjeu majeur en élaboration des métaux par voie liquide. La flottation, principal procédé retenu en métallurgie secondaire pour éliminer les particules d’inclusions, consiste à injecter des bulles de gaz au sein du réacteur. Lors de leur ascension, les bulles vont capter les plus grosses inclusions et favoriser la collision et l’agrégation des particules. Dans le but de quantifier les phénomènes influents à l’échelle des inclusions sur la dynamique d’agrégation entre deux inclusions à proximité des bulles, un modèle numérique 3D a été développé. L’écoulement local est modélisé par un cisaillement plan permanent et résolu par une méthode de Boltzmann sur réseau. Le couplage entre les particules et le fluide a été assuré par une méthode de frontière immergée permettant de calculer la perturbation hydrodynamique engendrée par la présence des particules et de mettre à jour les interactions entre particules pour leur suivi lagrangien. Les simulations numériques réalisées ont mis en évidence que les effets hydrodynamiques ont une influence non négligeable sur le comportement des inclusions. Des sections efficaces de collision ont pu être extraites, à partir desquelles ont été calculés des noyaux d’agrégation, données macroscopiques rendant compte des effets à petite échelle. Une première application de ce travail a été menée avec le calcul des fréquences d’agrégation d’un train de bulle dans un réacteur canal à partir de résultats de simulations DNS. Les noyaux d’agrégation ont également été exploités en vue de déterminer, à partir de résultats RANS de l’hydrodynamique d’une poche d’acier, l’évolution de la concentration d’inclusions par un bilan de population global. / Inclusion cleanliness remains a major challenge faced in process metallurgy in liquid phase. Flotation, the main process used in secondary metallurgy to remove inclusions, consists in injecting gas bubbles into the reactor. Rising gas bubbles entrap the biggest inclusions at their surface or in their wake. Besides, they promote collision and aggregation among particles. A 3D numerical model has been developed in order to quantify the roles of the prevailing phenomena on aggregation dynamics between inclusions in the vicinity of bubbles. At inclusion (mesoscopic) scale, the turbulent flow is locally modeled by a steady plane shear flow which is solved using a lattice-Boltzmann method. The coupling between both liquid and solid phases is ensured using an immersed boundary method. This method resolves the hydrodynamic perturbation induced by particles, and hence their interactions that are, in turn, used to update their Lagrangian tracking. The conducted numerical simulations bring out the influence of hydrodynamic effects on inclusion behavior. Collision cross sections have been determined from which ensuing aggregation kernels have been calculated. Such cross sections could provide macroscopic models to represent local particle dynamics. A first application of these results is presented to calculate aggregation frequencies in bubble swarms in a channel flow reactor that was simulated using DNS. Additionally, evolution of inclusion populations in molten steel has been determined from RANS simulation of a liquid steel ladle by means of a global population balance implementing the aggregation kernels determined in the present work.

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.<br /><br />La première partie introduit et analyse la méthode.<br /><br />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.<br /><br />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.

La réduction des émissions de gaz à effet de serre (GES) passe par la réduction des consommations d'énergie. Le stockage de la chaleur dans les parois des bâtiments permet de réduire la consommation d'énergie. Parmi les techniques de stockage, le stockage latent a la capacité de stocker une quantité d'énergie par unité de volume plus importante qu'un système sensible. Le projet INERTRANS a proposé le développement d'une façade associant une isolation translucide et le stockage latent avec un matériau à changement de phase (MCP). La fusion du MCP s'accompagne de la convection naturelle et l'absorption ou transmission du rayonnement. Le couplage de l'ensemble de ces phénomènes n'a pas été étudié dans la littérature. Dans cette thèse un modèle numérique 2D pour l'étude du changement de phase a été développé. Ce modèle utilise la méthode de Boltzmann sur réseau (LB) à temps de relaxation multiple (MRT), pour la résolution du champ de vitesse et la méthode des différences finies, pour la résolution du champ de températures. Le changement de phase a été traité par la formulation enthalpique. L'originalité est l'application de ce modèle au problème de changement de phase avec convection naturelle, d'une part, et au changement de phase avec convection naturelle et rayonnement, d'autre part. Pour vérifier notre modèle sans rayonnement, un cas de référence de la littérature a été simulé. Il s'agit de la fusion des deux MCP, l'étain et l'octadécane, à faible et fort nombre de Prandtl, respectivement. La simulation de l'étain a confirmé un écoulement multicellulaire. La simulation de l'octadécane a montré une forte influence de la convection avec un front de fusion qui se déforme sur toute la cavité. Le nombre de Nusselt pour l'octadécane avec convection est plus de trois fois le Nusselt sans convection. La simulation de l'acide gras de la brique INERTRANS a montré que la convection ne doit pas être négligée, car le flux prédit avec convection peut être jusqu'à trois fois plus grand que le flux prédit sans convection. La fraction fondue est près du double qu'en conduction seule. La méthode LB appliquée aux transferts radiatifs a été étudiée. Il se trouve, qu'à l'état actuel cette méthode n'est pas compétitive par rapport à une méthode classique des ordonnées discrètes (MOD). Enfin, nous avons couplé la MOD pour le calcul du flux radiatif avec la méthode LB pour le calcul du champ de vitesses et des différences finies pour l'équation de l'énergie. Le rayonnement grande longueur d'onde n'a pas d'influence notable sur les transferts thermiques. Le rayonnement courte longueur d'onde augmente les transferts thermiques, pourtant, cet effet n'est pas aussi important que l'augmentation due à la convection pour le matériau choisi. Puisqu'aucune solution de référence n'existe dans la bibliographie, nos résultats peuvent désormais servir d'éléments de comparaison pour de futurs travaux. Une validation expérimentale constituerait une perspective nécessaire.

La réduction des émissions de gaz à effet de serre (GES) passe par la réduction des consommations d’énergie. Le stockage de la chaleur dans les parois des bâtiments permet de réduire la consommation d'énergie. Parmi les techniques de stockage, le stockage latent a la capacité de stocker une quantité d’énergie par unité de volume plus importante qu’un système sensible. Le projet INERTRANS a proposé le développement d’une façade associant une isolation translucide et le stockage latent avec un matériau à changement de phase (MCP). La fusion du MCP s’accompagne de la convection naturelle et l’absorption ou transmission du rayonnement. Le couplage de l’ensemble de ces phénomènes n’a pas été étudié dans la littérature. Dans cette thèse un modèle numérique 2D pour l’étude du changement de phase a été développé. Ce modèle utilise la méthode de Boltzmann sur réseau (LB) à temps de relaxation multiple (MRT), pour la résolution du champ de vitesse et la méthode des différences finies, pour la résolution du champ de températures. Le changement de phase a été traité par la formulation enthalpique. L’originalité est l’application de ce modèle au problème de changement de phase avec convection naturelle, d’une part, et au changement de phase avec convection naturelle et rayonnement, d’autre part. Pour vérifier notre modèle sans rayonnement, un cas de référence de la littérature a été simulé. Il s’agit de la fusion des deux MCP, l’étain et l’octadécane, à faible et fort nombre de Prandtl, respectivement. La simulation de l’étain a confirmé un écoulement multicellulaire. La simulation de l’octadécane a montré une forte influence de la convection avec un front de fusion qui se déforme sur toute la cavité. Le nombre de Nusselt pour l’octadécane avec convection est plus de trois fois le Nusselt sans convection. La simulation de l’acide gras de la brique INERTRANS a montré que la convection ne doit pas être négligée, car le flux prédit avec convection peut être jusqu’à trois fois plus grand que le flux prédit sans convection. La fraction fondue est près du double qu’en conduction seule. La méthode LB appliquée aux transferts radiatifs a été étudiée. Il se trouve, qu’à l’état actuel cette méthode n’est pas compétitive par rapport à une méthode classique des ordonnées discrètes (MOD). Enfin, nous avons couplé la MOD pour le calcul du flux radiatif avec la méthode LB pour le calcul du champ de vitesses et des différences finies pour l’équation de l’énergie. Le rayonnement grande longueur d’onde n’a pas d’influence notable sur les transferts thermiques. Le rayonnement courte longueur d’onde augmente les transferts thermiques, pourtant, cet effet n’est pas aussi important que l’augmentation due à la convection pour le matériau choisi. Puisqu’aucune solution de référence n’existe dans la bibliographie, nos résultats peuvent désormais servir d’éléments de comparaison pour de futurs travaux. Une validation expérimentale constituerait une perspective nécessaire. / Reduction of greenhouse gas emissions requires reduction of energy consumption. Energy storage on building walls allows reduction in energy consumption. Among storage techniques, latent heat storage offers higher energy storage density than sensible heat storage. INERTRANS project has proposed the development of an innovative facade, coupling transparent insulation and energy storage with a fatty acid phase change materials (PCM). Melting of PCM comprises different phenomena, namely, natural convection in the liquid phase and radiation absorption or transmission. The coupling of all this phenomena is not still studied in scientific literature. In this thesis, a 2D numerical model for studying phase change has been developed. This model uses the lattice Boltzmann method (LBM) with multiple relaxation time (MRT) to resolve velocity field, and finite differences for the temperature field. Phase change is treated with the enthalpy formulation. The original contribution is application of this hybrid approach to the phase change with natural convection, on the one hand, and to the phase change with natural convection and radiation, on the other hand. To verify the model without radiation, a test case taken from literature has been simulated. It concerns the melting of two PCM with a low and high Prandtl number, the tin and octadecane, respectively. Tin melting simulation confirms multiple cells flow, starting with four rolls which merges in three then two rolls. Octadecane simulation shows high convection effect, with a melting front deforming all along the cavity height. Nusselt number plot for octadecane melting with convection is more than three times with conduction only. INERTRANS’ fatty acid simulation shows that convection shall not be neglected, because predicted heat flux with convection may be up to three times that predicted with conduction only. Melted fraction is almost twice than with conduction only. The lattice Boltzmann method applied to radiative heat transfer has also been explored. It turns out that in its current state, this method is not competitive compared to a conventional discrete ordinates method (DOM). Finally, we coupled the DOM for radiation heat flux, with the LBM for velocity field calculation and finite differences for the energy equation to solve the coupling between phase change, convection and radiation. Long wavelength radiation has no noticeable effect on heat transfer. Short wavelength radiation increases heat transfer, however, this increase is not as important as that produced by convection for this kind of material. Since no reference solution exists in the literature, our results can now serve as a basis for future work. An experimental validation would be a necessary perspective.

La sélection de variables est une tâche primordiale en fouille de données et apprentissage automatique. Il s’agit d’une problématique très bien connue par les deux communautés dans les contextes, supervisé et non-supervisé. Le contexte semi-supervisé est relativement récent et les travaux sont embryonnaires. Récemment, l’apprentissage automatique a bien été développé à partir des données partiellement labélisées. La sélection de variables est donc devenue plus importante dans le contexte semi-supervisé et plus adaptée aux applications réelles, où l’étiquetage des données est devenu plus couteux et difficile à obtenir. Dans cette thèse, nous présentons une étude centrée sur l’état de l’art du domaine de la sélection de variable en s’appuyant sur les méthodes qui opèrent en mode semi-supervisé par rapport à celles des deux contextes, supervisé et non-supervisé. Il s’agit de montrer le bon compromis entre la structure géométrique de la partie non labélisée des données et l’information supervisée de leur partie labélisée. Nous nous sommes particulièrement intéressés au «small labeled-sample problem» où l’écart est très important entre les deux parties qui constituent les données. / Feature selection is an important task in data mining and machine learning processes. This task is well known in both supervised and unsupervised contexts. The semi-supervised feature selection is still under development and far from being mature. In general, machine learning has been well developed in order to deal with partially-labeled data. Thus, feature selection has obtained special importance in the semi-supervised context. It became more adapted with the real world applications where labeling process is costly to obtain. In this thesis, we present a literature review on semi-supervised feature selection, with regard to supervised and unsupervised contexts. The goal is to show the importance of compromising between the structure from unlabeled part of data, and the background information from their labeled part. In particular, we are interested in the so-called «small labeled-sample problem» where the difference between both data parts is very important.

Au cours de cette thèse, un modèle LBM MPMC avec échanges thermiques est développé. Des tests d'assimilation de données et des mesures par flot optique sont réalisés en vue d'une validation. Le cadre d'application de cette thèse est celui du mélange d'un fluide cryogénique avec l'eau. Dans une première partie, un travail bibliographique rappelant l'équation de Boltzmann, ses diverses hypothèses et simplifications, ainsi que l'aspect algorithmique de la LBM sont exposés. Une comparaison entre opérateur de collision SRT et MRT est réalisée, et une simulation de phénomènes turbulents à différents nombres de Reynolds est étudiée, notamment avec le benchmark de l'instabilité de Von Karman. Dans une seconde partie, le modèle MPMC de Shan & Cehn est rappelé puis étendu au cas où les échanges thermiques entre composants sont présents. Des validations quantitatives sont faites, notamment avec le benchmark du fluide de Couette à deux phases ou à deux composants, du test de cohérence vis-à-vis de la loi de Laplace, ou encore par rapport à un benchmark faisant intervenir la conduction thermique. Des tests qualitatifs de condensation en milieu multicomposants sont proposés pour valider l'aspect des échanges thermiques entre composants en présence d'une transition de phase. Dans la troisième partie de cette thèse, une méthode de validation par assimilation de données est introduite, avec le filtrage de Kalman d'ensemble. Un test d'estimation d'état d'un fluide di-phasique est réalisé, et la compatibilité du filtrage de Kalman d'ensemble par rapport au modèle LBMMPMC est évaluée. Pour la validation du comportement du modèle d'un point de vue de la présence de deux composants, un fluide de substitution (non-cryogénique) au GNL, le butane, a été choisi pour permettre des observations dans des conditions expérimentales accessibles. Puis, une plateforme expérimentale d'injection de butane liquide dans une colonne d'eau sous pression est présentée. Des images d'ombroscopie issues d'expériences de remontée de butane liquide dans de l'eau sont exposées et un algorithme de calcul de flot optique est appliqué à ces images. Une évaluation qualitative des champs de vitesses obtenus par application de cet algorithme est réalisée. / In this thesis, a LBM MPMC model with heat exchange is developed. Data assimilation tests and optical flow measurements are made in order to validate the model. The application context of this thesis is the mixture of a cryogenic fluid with water. In the first part, a bibliographical work reminding the Boltzmann equation and its various assumptions and simplifications, as well as the algorithmic aspect of the LBM are exposed. A comparison between SRT and MRT collision operator is performed, and a simulation of turbulent phenomena at different Reynolds numbers is studied, especially with the benchmark of the instability from Von Karman. In the second part, the MPMC model from Shan & Chen is reminded and extended to the case of the inter-component heat exchanges. Quantitative validations are made, especially with the benchmark of a two-phase or two-component Couette fluid. Consistency is tested against Laplace's law rule, or against a benchmark involving heat conduction. Qualitative testing of condensations in a multi-component medium are proposed to validate the heat exchange between components in the presence of a phase transition. In the third part of this thesis, a validation method for data assimilation is introduced, with the ensemble Kalman filter. A state estimation test of a bi-phase fluid is realized, and compatibility of the ensemble Kalman filtering to the LBM MPMC model is assessed. For validation of the behavior of the model for a two-component case, a substitution fluid (non-cryogenic) for LNG, butane, was selected to permit observations in experimental conditions which are accessible. Then, an experimental platform of injection of liquid butane in a pressurised water column is presented. Shadowgraph images from liquid butane experiments in water are exposed and an optical flow calculation algorithm is applied to these images. A qualitative assessment of the velocity field obtaines by application of this algorithm is performed.

Ce travail de thèse a pour objet l’étude numérique et expérimentale de la déstabilisation de milieux granulaires immergés par fluidisation. Cette instabilité hydromécanique est un mécanisme précurseur de l’érosion régressive, processus de dégradation au coeur de la problématique de l’érosion interne des ouvrages hydrauliques en terre. La compréhension de ces mécanismes d’érosion nécessite une description rigoureuse du couplage et de l’interaction entre le fluide et les particules de sol. A cette fin, un modèle 2D a été utilisé en couplant deux méthodes particulaires, la méthode des éléments discrets (DEM) pour modéliser le comportement mécanique de la phase solide et la méthode Lattice Boltzmann (LBM) pour la phase fluide. Des expériences servant de validation à cette simulation numérique 2D ont également été réalisées en s’appuyant sur une technique de visualisation interne d’un empilement granulaire combinant l’ajustement d’indice de réfraction des deux phases et la fluorescence induite par plan laser. / The subject of this thesis is the numerical analysis and experimental investigation of the destabilization of submerged granular media caused by fluidization. This hydromechanical instability is one of the mechanisms that may trigger the regressive erosion, which is one of the main degradation phenomena driving the internal erosion of earthen hydraulic constructions. Such erosion mechanisms can only be understood through a rigorous description of the coupling and interaction between the eroding fluid and the soil particles. For this purpose, a 2D model has been used coupling two different numerical techniques, namely the discrete element method (DEM) for modelling the mechanical behaviour of the solid phase and the Lattice Boltzmann method (LBM) for the fluid phase. The experimental validation of this numerical 2D simulation has been carried out using two optical techniques for the internal visualization of a granular sample, namely the adjustment of the refraction index of the two phases and the laser-induced fluorescence.

La méthode de Boltzmann sur réseau est une formulation discrète particulière de l'équation de Boltzmann. Depuis ses débuts, il y a trente ans, cette méthode a gagné une certaine popularité, et elle est maintenant utilisée dans presque tous les problèmes habituellement rencontrés en mécanique des fluides notamment pour les écoulements multi-espèces. Dans le cadre de ce travail, une force de friction intermoléculaire est introduite pour modéliser les interactions entre les molécules de différent types causant principalement la diffusion entre les espèces. Les phénomènes de dissipation visqueuse (collision usuelle) et de diffusion moléculaire (force de friction intermoléculaire) sont séparés et peuvent être ajuster indépendamment. Le principal avantage de cette stratégie est sa compatibilité avec des optimisations de la collision usuelle et les opérateurs de collision avancés. Adapter un code mono-espèce pour aboutir à un code multi-espèces est aisé et demande beaucoup moins d’effort comparé aux précédentes tentatives. De plus, il n’ y a pas d’approximation du mélange, chaque espèce a ses propres coefficients de transport pouvant être calculés à l’aide de la théorie cinétique des gaz. En général, la diffusion et la convection sont vus comme deux mécanismes séparés : l’un agissant sur la masse d’une espèce, l’autre sur la quantité de mouvement du mélange. En utilisant une force de friction intermoléculaire, la diffusion et la convection sont couplés par l’intermédiaire la quantité de mouvement de chaque espèce. Les mécanismes de diffusion et de convection sont intimement liés dans de nombreux phénomènes physique tel que la digitation visqueuse.L’instabilité de digitation visqueuse est simulée en considérant dans un milieu poreux deux espèces dans des proportions différentes soit un mélange moins visqueux déplaçant un mélange plus visqueux. Les principaux moteurs de l’instabilité sont la diffusion et le contraste de viscosité entre les espèces. Deux stratégies sont envisagées pour simuler les effets d’un milieu poreux. Les méthodes de rebond partiel et de force de Brinkman bien que basées sur des approches fondamentalement différentes donnent dans notre cas des résultats identiques. Les taux de croissance de l’instabilité calculés à partir de la simulation coïncident avec ceux obtenus à partir d’analyses de stabilité linéaire. L’évolution de la longueur de mélange peut être divisée en deux étapes dominées d’abord par la diffusion puis par la convection. La physique de la digitation visqueuse est ainsi correctement simulée. Toutefois, les effets de diffusion multi-espèces ne sont généralement pas pris en compte lors de la digitation visqueuse de trois espèces et plus. Ces derniers ne sont pas négligeable puisque nous mettons en avant une configuration initialement stable qui se déstabilise. La diffusion inverse entraîne la digitation dont l’impact dépend de la diffusion entre les espèces. / The lattice Boltzmann method (LBM) is a specific discrete formulation of the Boltzmann equation. Since its first premises, thirty years ago, this method has gained some popularity and is now applied to almost all standard problems encountered in fluid mechanics including multi-component flows. In this work, we introduce the inter-molecular friction forces to take into account the interaction between molecules of different kinds resulting primarily in diffusion between components. Viscous dissipation (standard collision) and molecular diffusion (inter-molecular friction forces) phenomena are split, and both can be tuned distinctively. The main advantage of this strategy is optimizations of the collision and advanced collision operators are readily compatible. Adapting an existing code from single component to multiple miscible components is straightforward and required much less effort than the large modifications needed from previously available lattice Boltzmann models. Besides, there is no mixture approximation: each species has its own transport coefficients, which can be calculated from the kinetic theory of gases. In general, diffusion and convection are dealt with two separate mechanisms: one acting respectively on the species mass and the other acting on the mixture momentum. By employing an inter-molecular friction force, the diffusion and convection are coupled through the species momentum. Diffusion and convection mechanisms are closely related in several physical phenomena such as in the viscous fingering instability.A simulation of the viscous fingering instability is achieved by considering two species in different proportions in a porous medium: a less viscous mixture displacing a more viscous mixture. The core ingredients of the instability are the diffusion and the viscosity contrast between the components. Two strategies are investigated to mimic the effects of the porous medium. The gray lattice Boltzmann and Brinkman force models, although based on fundamentally different approaches, give in our case equivalent results. For early times, comparisons with linear stability analyses agree well with the growth rate calculated from the simulations. For intermediate times, the evolution of the mixing length can be divided into two stages dominated first by diffusion then by convection, as found in the literature. The whole physics of the viscous fingering is thus accurately simulated. Nevertheless, multi-component diffusion effects are usually not taken into account in the case of viscous fingering with three and more species. These effects are non-negligible as we showcase an initial stable configuration that becomes unstable. The reverse diffusion induces fingering whose impact depends on the diffusion between species.

La clairance mucociliaire est un processus physico-chimique qui sert à transporter et éliminer le mucus bronchique. Pour cela, des milliards d'appendices de taille micrométrique, que l'on nomme cils, recouvrent l'épithélium respiratoire. Ces cils propulsent le mucus en suivant un motif périodique comprenant une phase de poussée où leur pointe peut pénétrer dans le mucus, et une phase de récupération où ils sont totalement immergés dans le fluide périciliaire. Un dysfonctionnement de ce processus peut engendrer de nombreux problèmes de santé. Il a été observé expérimentalement que les cils ne battent pas aléatoirement, mais synchronisent leurs battements avec leurs voisins, formant ainsi des ondes métachronales. Toutefois, les observations in vivo sont extrêmement difficiles à réaliser, et les propriétés de ces ondes restent mal connues. Dans cette thèse, nous utilisons un solveur Lattice Boltzmann - Frontière Immergée afin de reproduire un épithélium bronchique et étudier l'émergence, ainsi que les capacités de transports et de mélanges, de ces ondes / The mucociliary clearance process is a physico-chemical process which aims is to transport and eliminate bronchial mucus. To do so, billions of micro-sized appendages, called cilia, cover the respiratory epithelium. These cilia propel the mucus by performing a periodical pattern composed of a stroke phase where their tips can enter the mucus layer, and a recovery phase where the cilia are completely immersed in the periciliary liquid layer. A failure of this process may induce numerous health problems. It has been experimentally observed that cilia do not beat randomly, but instead adapt their beatings accordingly to their neighbours, forming metachronal waves. However, in vivo observations are extremely difficult to perfom, and the properties of these waves remain poorly understood. In this thesis, we use a Lattice Boltzmann - Immersed Boundary solver to reproduce a bronchial epithelium and study the emergence, as well as the transport and mixing capacities, of these waves