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

Ce manuscrit est consacré au développement et à la validation d'un nouveau modèle de collision destiné à améliorer la stabilité des modèles lattice Boltzmann (LB) d'ordre élevés lors de la simulation d'écoulements : (1) isothermes et faiblement compressibles à nombre de Reynolds élevés, ou (2) compressibles et comprenant des discontinuités telles que des ondes de choc. Ce modèle de collision s'appuie sur une étape de régularisation améliorée. Cette dernière inclut désormais un calcul par récurrence des coefficients hors-équilibre du développement en polynômes d'Hermite. Ces formules de récurrence sont directement issues du développement de Chapman-Enskog, et permettent de correctement filtrer les contributions non-hydrodynamiques émergeant lors de l'utilisation de maillages sous-résolus. Cette approche est d'autant plus intéressante quelle est compatible avec un grand nombre de réseaux de vitesses discrètes. Ce modèle LB d'ordre élevé est validé tout d'abord pour des écoulements isothermes à nombre de Reynolds élevé. Un couplage avec une technique de capture de choc permet ensuite d'étendre son domaine de validité aux écoulements compressibles incluant des ondes de choc. Ce travail se conclut avec une étude de stabilité linéaire des modèles considérés, le tout dans le cas d'écoulements isothermes. Ceci permet de quantifier de manière distincte l'impact des discrétisations en vitesse et numérique, sur le comportement spéctrale du jeu d'équations associé. Cette étude permet au final de confirmer le gain en stabilité induit par le nouveau modèle de collision. / This thesis is dedicated to the derivation and the validation of a new collision model as a stabilization technique for high-order lattice Boltzmann methods (LBM). More specifically, it intends to stabilize simulations of: (1) isothermal and weakly compressible flows at high Reynolds numbers, and (2) fully compressible flows including discontinuities such as shock waves. The new collision model relies on an enhanced regularization step. The latter includes a recursive computation of nonequilibrium Hermite polynomial coefficients. These recursive formulas directly derive from the Chapman-Enskog expansion, and allow to properly filter out second- (and higher-) order nonhydrodynamic contributions in underresolved conditions. This approach is even more interesting since it is compatible with a very large number of velocity sets. This high-order LBM is first validated in the isothermal case, and for high-Reynolds number flows. The coupling with a shock-capturing technique allows to further extend its validity domain to the simulation of fully compressible flows including shockwaves. The present work ends with the linear stability analysis(LSA) of the new approach, in the isothermal case. This leads to a proper quantification of the impact induced by each discretization (velocity and numerical) on the spectral properties of the related set of equations. The LSA of the recursive regularized LBM finally confirms the drastic stability gain obtained with this new approach.

La modélisation des écoulements et du transport dans les milieux poreux est un domaine actif pour, notamment, progresser dans la compréhension du transfert des polluants dans les sols. Les sols présentent fréquemment des hétérogénéités comme des macropores (provoqués par la faune, la flore ou des fissures) et un certain nombre de modèles numériques utilisent les concepts de double ou de multi-perméabilité pour tenir compte de tous les types d’écoulements susceptibles de coexister dans de tels systèmes. Cependant, les modèles classiques semblent sous-estimer l’effet de la macroporosité sur l’écoulement et le transfert préférentiels et restreindre la zone d’écoulement préférentiel au seul volume occupé par la macroporosité. Diverses études expérimentales antérieures à cette thèse ont questionné cette hypothèse. Cette étude se propose de comprendre l’établissement de l’écoulement et du transport préférentiel et en particulier les mécanismes d’échange d’eau et de masse entre un macropore et une matrice poreuse environnante en condition saturée. Pour cela, des traçages de l’eau sont réalisés pour un milieu poreux modèle constitué de billes de verre, traversé par un macropore synthétique et mis en place en colonnes de laboratoire. Elution et transfert dans les colonnes sont caractérisés par suivi de la concentration en sortie et par imagerie par résonance magnétique. Un modèle numérique développé sur la base de la méthode de Boltzmann sur réseau est utilisé pour simuler numériquement des écoulements dans un système macroporé et identifier les mécanismes d’écoulements préférentiels à l’échelle de pores. Les données expérimentales montrent que le transfert du traceur est fortement dépendant du débit d’injection ainsi que du coefficient de diffusion dans l’eau. À fort débit, le transfert semble s’effectuer exclusivement dans le macropore, avec très peu d’échange avec la matrice. Pour des débits plus faibles, la percée présente une inflexion suivie d’un pic. Les images IRM montrent alors un échange significatif de traceur entre le macropore et la matrice poreuse environnante. Les simulations numériques sont utilisées pour calculer le champ de vitesse de l’écoulement dans le système en fonction du débit. Les modélisations numériques montrent que l’écoulement préférentiel est étendu dans la matrice poreuse sur une zone de même dimension que le diamètre moyen des grains indépendamment de la taille du macropore et du débit, dans la gamme de débits simulés. Ces résultats expérimentaux et numériques montrent que l’influence du macropore sur les transferts doit être étendue dans la matrice poreuse sur une zone de la taille des grains pour l’écoulement et sur une zone dépendant du coefficient de diffusion du traceur ainsi que du temps de séjour moyen de celui-ci pour le transfert des solutés / Flow and transport modeling through porous media is of primary concern nowadays, especially in order to progress in the understanding of pollutant transfers through soils. Soils present frequently heterogeneities such as macropores (caused by fauna, flora or cracks) and several numerical models use double or multi permeability concepts in order to take into account all flow types that may exist in such porous systems. Nevertheless, classical models seem underestimate the macropore effect on preferential flow and transport by restricting the preferential flow zone only to the volume occupied by the macroporosity. Various experimental studies prior to this thesis have questioned this hypothesis. This study proposes to understand the establishment of preferential flow and transport and in particular the mechanism of flow and solute exchanges between a synthetic macropore and a surrounding porous matrix in saturated condition. For this purpose, water tracing are realized for a model porous media constituted by glass beads, crossed by a synthetic macropore and implemented in laboratory columns. Breakthrough and transport in columns are characterized by monitoring the concentration at the end of the column by magnetic nuclear resonance. A numerical model developed on the basis of lattice-Boltzmann method is used to simul ate flow in macroporous system and identify preferential flow mechanisms at pore scale. Experimental data show that tracer transport is strongly dependent on injection flow rate and the diffusion coefficient in water. At high flow rate, the transport seems to occur exclusively in the macropore, with very little masse exchange with the porous matrix. At lower flow rates, the breakthrough exhibits an inflexion followed by a peak. The MRI images show a significant mass exchange of tracer between the macropore and the surrounding porous matrix. The numerical simulations are used to calculate the flow field in a porous system as a function of flow rate. They show that preferential flow is extended in porous matrix into a zone of same dimension the mean diameter of beads regardless of macropore size or injected flow rate, in the range of simulated flow rates. These experimental and numerical results show that macropore influence on transport should be extended through the surrounding porous matrix into a zone of the same size of grains diameter for flow and into a zone depending on diffusion coefficient as well as mean residence time of the studied tracer for solute transport

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

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.

La circulation sanguine joue un rôle vital en microcirculation, et ce pour le transport de l'oxygène, le dioxide de carbone et d'autres nutriments. Les globules rouges (GR) constituent la majorité des cellules du sang, c'est pourquoi par "écoulement sanguin", nous entendrons "écoulement d'une suspension de GR". Pendant longtemps l'écoulement sanguin était vu comme un phénomène passif où les GR sont considérés comme des cargos d'oxygène. La vision moderne est tout autre: l'écoulement sanguin est bel et bien un phénomène actif. Les GR ainsi que les cellules endothéliales (qui tapissent les faces internes des vaisseaux sanguins) sont impliquées dans un grand nombre de signalisations biochimiques induites par les contraintes hydrodynamiques, la route vers des régulations vasomotrices sans l'intervention du système nerveux. Par exemple, les GR ne transportent pas que l'oxygène, mais également de l'ATP (adenosine triphosphate), qui est libérée suite à des changements de conformation de protéines membranaires induite par les contraintes hydrodynamiques. Cette thèse est dédiée à la circulation sanguine et son couplage avec la signalisation biochimique ayant lieu en microcirculation. Plus précisément, les questions traités dans cette thèse sont i) la dynamique des GR, ii) le problème de la diffiusion-advection d'espèces chimiques au sein des écoulements sanguins, et iii) le rôle de la géométrie des réseaux vasculaires dans le processus de la signalisation biochimique mentionnés plus haut. Dans un premier temps nous analysons la dynamique de GR dans un écoulement de Poiseuille en présence de valeurs réalistes de contraste de viscosité. Dans un deuxième temps nous développons un modèle de diffusion-advection et le couplons aux écoulements sanguins en adoptant la méthode de Boltzmann sur réseaux; nous exploitons ensuite formulation en l'appliquant au problème de la libération de l'ATP par les GR sous écoulement. Enfin nous présentons des résultats préliminaires pour la problématique générale de l'écoulement sanguin mettant en jeu l'ATP libéré par les GR et la signalisation de calcium par les cellules endothéliales. Cette étude constitue un premier pas vers le problème général et ambitieux de la régulation locale mechano-biochimique impliquée dans la microcirculation. / Blood flow in microcirculation is vital for oxygen, carbon dioxide and nutrients transport. Most of blood cells are red blood cells (RBCs), so that by blood flow we mean flow of a suspension of RBCs. For long time blood flow has been mainly considered as a passive phenomenon, in which RBCs are viewed as passive carriers of oxygen. The modern view is completely different: blood flow is more active than we thought. The RBCs as well as vascular endothelial cells covering the internal walls of blood vessels are involved in a number of biochemical signaling processes that are triggered by shear stress eliciting a number of biochemical events, and ultimately resulting into vasomotor regulation without participation of the nerve system. For example, RBCs do not only carry oxygen but also ATP (adenosine triphosphate) , the release of which occurs thanks to changes of RBC membrane protein conformations caused by shear stress. Released ATP reacts with some endothelial membrane receptors leading to vasodilation. This thesis is devoted to blood flow and its coupling to biochemical signaling. More precisely, we investigate i) the dynamics of RBCs, ii) the advection diffusion of chemicals in blood flow and the role of iii) the geometry of vessel networks, in the mentioned signaling processes in microcirculations. Firstly, we study the RBC dynamics in a pipe flow with realistic viscosity contrast values, where a link between shape dynamics and rheology is established. Secondly, we develop an advection-diffusion solver that can handle general moving curved boundaries based on lattice-Boltzmann method (LBM); we then implement it for the study of the problem of ATP release from RBCs under shear flow. Membrane tension and deformation induced by shear stress together with vessel network geometry contribute to ATP release. Finally we demonstrate the capability of applying our model and our numerical tool to the complete problem of blood under flow involving ATP release from RBCs and endothelial calcium signaling as a preliminary step to the ambitious task of mechano-involved local regulation events in microcirculation.

Au cours des deux dernières décennies, le contexte économique a changé de manière significative en raison de la hausse des prix de l'énergie. Le bâtiment étant devenu le principal secteur consommateur d'énergie, la réduction de celle-ci est devenue un objectif économique, sociétal et environnemental. Ce sujet mobilise de nombreux travaux de recherche. Les Matériaux à Changement de Phase (MCP) représentent une solution innovante qui pourrait contribuer à améliorer la performance énergétique des bâtiments. Ils sont principalement utilisés pour la régulation de température, et leur forte capacité de stockage est un moyen de réduire la consommation d'énergie. Notre étude vise à caractériser, via une approche expérimentale et numérique, le comportement d'un PCM (l’Octadécane). Pour cela, nous avons développé et mis en œuvre un modèle numérique qui corrobore les résultats expérimentaux, et ainsi améliore la prédiction de la performance du MCP considéré.Dans ce travail, notre principale préoccupation est de mettre en évidence les erreurs ou simplifications présentes dans le modèle numérique traditionnel pouvant entraîner un écart global par rapport au comportement réel du MCP. Ces différences conduisent à une estimation erronée des temps de fusion et de la quantité d'énergie stockée. L'amélioration significative de notre modèle est la prise en compte de la convection naturelle, de la surfusion, et l'utilisation des courbes réelles d'enthalpie du MCP considéré. La relation température-enthalpie réelle tient compte de la présence d'une fraction d'impuretés solubles dans le matériau. L’originalité de ce travail est de traiter ces phénomènes physiques via la méthode de Boltzmann réseau (connue sous l'acronyme LBM) avec des fonctions de distribution doubles couplée à une formulation enthalpique. Une telle approche permet de passer outre la non-linéarité des équations régissant l'écoulement et le transfert de chaleur. Sa simplicité de mise en œuvre et son caractère local permettent d'affiner le modèle. Ainsi, on peut couvrir les problèmes de changement de phase, y compris ceux pouvant avoir lieu dans des matrices poreuses ou fibreuses. Ce dernier point a été couvert dans cette thèse.Enfin, il s'est avéré que l'approche numérique adoptée ici pour traiter les problèmes de changement de phase corrobore à la fois nos résultats expérimentaux et ceux disponibles dans la littérature. / Over the past two decades, the economic context has changed significantly due to the rise in energy prices. The building sector has become the main consumer of energy. Thereby, reducing the latter is now an economic, societal and environmental necessity. Accordingly, this topic mobilizes many researches. Phase Change Materials (PCMs) represent an innovative solution, which could improve buildings' energy performance. They are primarily used for temperature regulation, and their high storage capacity can reduce energy consumption.Our study aims at characterizing, via a complementary approach of experimental and numerical simulation, the behavior of a PCM (n-Octadecane). For this, we have developed and implemented a numerical model that corroborates the experimental results, and hence improves the prediction of the PCM performance.In this work, our main concern is to highlight the common errors or simplifications taken in the traditional numerical model, which can result in an overall discrepancy compared to the actual behavior of PCMs. Those discrepancies lead to wrong estimation of the fusion times and amount of energy stored. The major improvement of our model is the consideration of the natural convection, the supercooling, and the use of real enthalpy curves of the considered PCM. The actual temperature-enthalpy relationship takes into account the presence of a fraction of soluble impurities in the material. The originality of this work is to handle these physical phenomena via a lattice Boltzmann method (known by the acronym LBM), which leans on double distribution functions and coupled with the enthalpy formulation. Such an approach overcomes the non-linearity in the governing equations of fluid flow and heat transfer. Its simplicity and local character allow adding complexity to the model. Thereby, one can cover up the phase change problems, including those, which may occur in heterogeneous matrices. This last point has been also covered in this thesis.Finally, it turned out that the approach implemented here for phase change problems supports both, our experimental results and those available in the literature.

Les nids des insectes sociaux fascinent car, aussi complexes et structurés qu'ils puissent être, ils sont construits de façon auto-organisée par des insectes qui peuvent être de plusieurs ordres de grandeur plus petits que la structure. Nous abordons ici l'étude des mécanismes de comportement individuel impliqués dans ces processus de morphogenèse à travers l'exemple de l'agrégation de cadavres par une colonie de fourmis Messor Sancta. Nous proposons un modèle de marche aléatoire non réciproque faisant émerger un suivi des structures, appuyé sur des observations individuelles et dont nous étudions les propriétés statistiques. Nous modélisons également l'interaction entre les actes de construction et un écoulement d'air de faible vitesse comme un exemple des couplages écoulement-constructions supposés responsables de l'adaptation des nids aux conditions environnementales.

Dans de nombreux réservoirs géothermiques profonds, l'échangeur thermique est un massif de roches chaudes fracturées. Comment la morphologie des fractures influence-t-elle le flux hydraulique et le champ de température lorsque de l'eau froide est injectée dans ce milieu ? Cette question est abordée à l'échelle de la fracture, de manière numérique, à l'aide de deux méthodes : par différences finies, en 2D, et par méthode de Boltzmann sur réseau, en 3D. Les calculs en différences finies montrent qu'une ouverture auto-affine induit une variabilité spatiale du flux hydraulique et du champ de température. Dans certains cas, il se crée une chenalisation du fluide, qui réduit l'efficacité du transfert thermique entre la roche et le fluide. Ces caractéristiques hydro-thermiques sont essentiellement engendrées par les plus grandes longueurs d'ondes de ces ouvertures auto-affines. Une première application à grande échelle, dimensionnée suivant les caractéristiques du site géothermique de Soultz-sous-Forêts (Alsace, France), est proposée. Une deuxième application consiste à évaluer la perméabilité du sous-sol de Draix (Alpes, France). Pour cela, la géométrie de fractures présentes dans le sous-sol est caractérisée précisément à partir de carottes géologiques. La topographie des fractures, mesurée avec un laser, est utilisée pour reconstituer la morphologie de l'ouverture des fractures. Ces données sont utilisées pour quantifier la perméabilité, ici élevée -- de l'ordre de 10^-9-10^-8 m^2. À l'aide de méthodes de Boltzmann sur réseau, la modélisation du comportement hydro-thermique dans des fractures présentant une topographie avec de fortes pentes est abordée. Au voisinage d'une aspérité de forme simple, le comportement hydro-thermique modélisé est très différent (recirculation) de celui dans le reste de la fracture. L'ensemble de ces modélisations suggère qu'il est nécessaire de prendre en compte la morphologie des fractures pour estimer le comportement hydro-thermique en géothermie.

La reconstruction tridimensionnelle d'un échantillon représentatif d'un matériau céramique biphasique à partir d'une image de structure bidimensionnelle est proposée. L'algorithme préconisé est de nature stochastique fondé sur une minimisation par recuit simulé. La méthode statistique de Boltzmann sur réseau est présentée pour simuler les phénomènes de transport de matière afin d'estimer la perméabilité de l'échantillon reconstruit par un modèle BR D3Q19. La conductivité thermique « effective » du même domaine 3D est évaluée par la méthode des différences finies dans un schéma Dufort-Frankel. Les exemples étudiés montrent l'aptitude de l'outil de la reconstruction associé aux outils de simulations des phénomènes de transport de chaleur et de masse à synthétiser numériquement un matériau hétérogène ou poreux et à déterminer ses paramètres physiques.