En raison de la tendance positives pour le développement durable, les systèmes solaires(capteurs solaires, concentrateur solaire, etc.) Intègrent (et demandent d'intégrer encore plus)intensivement les résidences et les industries. Dans ce contexte, les systèmes diphasiques comme le caloduc semblent être très efficaces en raison de leurs capacités élevées de transport de chaleur et de leur fonctionnement passif appliqués aux capteurs. Compte-tenu de la complexité des caloducs à structure poreuse dans ce type d'application, la plupart des systèmes existants sur le marché utilisent des thermosiphons. Ainsi, le besoin croissant de solutions de contrôle thermique fiables et plus efficaces croit rapidement pour de tels systèmes.Ce travail de thèse porte sur la caractérisation des performances des caloducs à structure poreuse utilisés dans les applications solaires. Une étude numérique a été réalisée pour modéliser et simuler le comportement d'un caloduc typique à l'aide de la méthode Lattice Boltzmann. Une étude expérimentale a également été réalisée pour caractériser les performances de trois prototypes testes dans différentes conditions (température du condenseur, puissance introduite et angle d'inclinaison). Les effets induits par plusieurs paramètres incluant le taux de remplissage, le fluide de travail et la symétrie de la puissance appliquée sur les performances de ces dispositifs ont également été étudiés. En particulier, l'asymétrie du chauffage induit un assèchement plus précoce, toutes choses étant égales par ailleurs. L'inclinaison optimale est également déterminée là où est équilibrée la chaleur solaire maximale disponible et reçue par le caloduc et l'écoulement de liquide assisté par gravité à l'intérieur de ce dispositif. / Owing to the trend to development sustainability, solar systems (solar collector, solar concentrator, etc.) Are integrating (and asked to integrate even more) intensively residences and industries. In this context, two-phase systems like heat pipe seem highly effective because of their high heat transport capabilities and their passive operation in collectors’ technology. In view of the complexity of the heat pipes with a porous structure in this kind of application,most of the existing systems on the market use thermosyphons. Thus, the growing need of reliable and more efficient thermal control solutions is increasing for such systems. This thesis work focuses on the performance characterization of heat pipes with porous structure used in solar applications. A numerical study has been performed to model and simulate the behavior of a typical heat pipe using the Lattice Boltzmann method. An experimental study has also been done to characterize the performance of three prototypes tested under different conditions (condenser temperature, heat input and inclination angle). The effects induced by several parameters including the filling rate, working fluid and symmetry of the applied heat on the performance of these devices has also been investigated. In particular, heating asymmetry is found to induce dry-out earlier, all other things being equal. Optimal inclination is also determined where is balanced the maximum solar heat available and received by the heat pipe and the gravity-assisted liquid flow inside that device.

Les extensions bifluides des méthodes Lattice Boltzmann à frontière libre utilisent généralement des pseudopotentiels microscopiques pour modéliser l'interface. Nous avons choisi d'orienter nos recherches vers une méthode Lattice Boltzmann à capture d'interface où la fraction massique d'un des deux fluides, inconnue, est transportée. De nombreux travaux ont montré les difficultés des méthodes Lattice Boltzmann à traiter des systèmes bifluides, et ce d'autant plus que le ratio de densité est important. Nous expliquerons l'origine de ces problèmes en mettant en évidence le manque de diffusion numérique pour capturer précisément les discontinuités de contact. Pour régler cet obstacle, nous proposerons une formulation Arbitrary Lagrangian Eulerian (ALE) des méthodes Lattice Boltzmann. Cela permet de séparer le traitement des ondes matérielles de celui des ondes de pression. Une fois l'étape ALE terminée, une phase de projection ramène les variables sur la grille eulérienne de calcul initiale. Nous expliquons comment obtenir une procédure de projection ayant une précision d'ordre 2 et une interface fine et dépourvue d'oscillations. Il sera montré que la fraction massique satisfait un principe du maximum discret et qu'elle reste donc entre 0 et 1. Les simulations numériques sont en accord avec la théorie. Même si notre méthode n'est pour le moment utilisée que pour simuler des écoulements de fluides non visqueux (Equations d'Euler), nous sommes convaincus qu'elle pourra être étendue à des simulations d'écoulements bifluides visqueux. / Two-fluid extensions of Lattice Boltzmann methods with free boundaries usually consider ``microscopic'' pseudopotential interface models. In this paper, we rather propose an interface-capturing Lattice Boltzmann approach where the mass fraction variable is considered as an unknown and is advected. Several works have reported the difficulties of LBM methods to deal with such two-fluid systems especially for high-density ratio configurations. This is due to the mixing nature of LBM, as with Flux vector splitting approaches for Finite Volume methods. We here give another explanation of the lack of numerical diffusion of Lattice Boltzmann approaches to accurately capture contact discontinuities. To fix the problem, we propose an arbitrary Lagrangian-Eulerian (ALE) formulation of Lattice-Boltzmann methods. In the Lagrangian limit, it allows for a proper separated treatment of pressure waves and advection phenomenon. After the ALE solution, a remapping (advection) procedure is necessary to project the variables onto the Eulerian Lattice-Boltzmann grid.We explain how to derive this remapping procedure in order to get second-order accuracy and achieve sharp stable oscillation-free interfaces. It has been shown that mass fractions variables satisfy a local discrete maximum principle and thus stay in the range $[0,1]$. The theory is supported by numerical computations of rising bubbles (without taking into account surface tension at this current state of development).Even if our methods are currently used for inviscid flows (Euler equations) by projecting the discrete distributions onto equilibrium ones at each time step, we believe that it is possible to extend the framework formulation for multifluid viscous problems. This will be at the aim of a next work.

Les études présentées dans ce mémoire portent sur la simulation numérique et l'analyse physique des milieux granulaires immergés dans un fluide. Des développements numériques ont été réalisés pour coupler la méthode Lattice Boltzmann pour la dynamique du fluide avec la méthode Contact Dynamics en 2D et avec la méthode Molecular Dynamics en 3D pour la dynamique des grains. Ces outils numériques ont été utilisés pour étudier l'initiation des avalanches sur un plan incliné en fonction de la compacité initiale et de l'angle d'inclinaison en 3D. Les résultats sont en bon accord quantitatif avec les expériences et ont permis de mettre en évidence la stabilisation de la pente granulaire par une pression négative du fluide interstitielle induite par la dilatance, et l'évolution spatiotemporelle des grandeurs telles que la compacité et la déformation de cisaillement. Ces évolutions dans la phase de fluage qui précède la rupture de pente ont pu être mises à l'échelle par un modèle théorique incorporant la loi de Darcy et l'effet de la dilatance sur l'angle de frottement interne. L'analyse de la texture granulaire a révélé la distortion du réseau des contacts pendant le fluage et la saturation de l'anisotropie comme un critère de rupture. La propagation des avalanches granulaires a été étudiée dans une configuration 2D pour deux géométries différentes : 1) l'effondrement et l'étalement d'une colonne sous son propre poids, 2) l'étalement d'une pente sous l'effet d'une énergie cinétique injectée. Nous avons en particulier montré que la distance et la durée d'étalement obéissent à des lois de puissance en fonction du rapport d'aspect initial ou de l'énergie injectée. Le fluide exerce deux effets contradictoires : réduire les temps de relaxation et lubrifier les contacts. Ces effets ont été analysés dans le régime visqueux en fonction des conditions initiales et la viscosité du fluide.

L'objet de cette thèse est la modélisation computationnelle de l'Imagerie par Résonance Magnétique (IRM), appliquée à l'imagerie des réseaux vasculaires. Les images sont influencées par la géométrie des vaisseaux mais aussi par le flux sanguin. Par ailleurs, outre la qualité des modèles développés, il est important que les calculs soient performants. C'est pourquoi, le calcul parallèle est utilisé pour gérer ce type de problèmes complexes. Dans cette thèse, trois solutions sont proposées. La première concerne les algorithmes parallèles pour la modélisation des réseaux vasculaires. Des algorithmes dédiés à différentes architectures sont proposés. Le premier est basé sur le modèle de « passage de messages » pour les machines à mémoires distribuées. La parallélisation concerne l'irrigation de nouvelles zones de tissu par les vaisseaux existants. Le deuxième algorithme est dédié aux machines à mémoire partagée. Il parallélise également le processus de perfusion mais des processeurs différents se chargent de gérer les différents arbres vasculaires. Le troisième algorithme est une combinaison des approches précédentes offrant une solution pour les architectures parallèles hybrides. Les algorithmes proposés permettent d'accélérer considérablement la croissance des réseaux vasculaires complexes, ce qui rend possible la simulation de structures vasculaires plus précises, en un temps raisonnable et aide à améliorer le modèle vasculaire et à tester plus facilement différents jeux de paramètres. Une nouvelle approche de modélisation computationnelle des flux en IRM est également proposée. Elle combine le calcul de flux par la méthode de Lattice Boltzmann, la simulation IRM par le suivi temporel de magnétisations locales, ainsi qu'un nouvel algorithme de transport des magnétisations. Les résultats montrent qu'une telle approche intègre naturellement l'influence du flux dans la modélisation IRM. Contrairement aux travaux de la littérature, aucun mécanisme additionnel n'est nécessaire pour considérer les artéfacts de flux, ce qui offre une grande facilité d'extension du modèle. Les principaux avantages de cette méthode est sa faible complexité computationnelle, son implémentation efficace, qui facilitent le lancement des simulations en utilisant différents paramètres physiologiques ou paramètres d'acquisition des images. La troisième partie du travail de thèse a consisté à appliquer le modèle d'imagerie de flux à des réseaux vasculaires complexes en combinant les modèles de vaisseaux, de flux et d'acquisition IRM. Les algorithmes sont optimisés à tous les niveaux afin d'être performants sur des architectures parallèles. Les possibilités du modèle sont illustrées sur différents cas. Cette démarche de modélisation peut aider à mieux interpréter les images IRM grâce à l'intégration, dans les modèles, de connaissances variées allant de la vascularisation des organes jusqu'à la formation de l'image en passant par les propriétés des flux sanguins. / This PhD thesis concerns computer modeling of magnetic resonance imaging (MRI). The main attention is centered on imaging of vascular structures. Such imaging is influenced not only by vascular geometries but also by blood flow which has to been taken into account in modeling. Next to the question about the quality of developed models, the challenge lies also in the demand for high performance computing. Thus, in order to manage computationally complex problems, parallel computing is in use. In the thesis three solutions are proposed. The first one concerns parallel algorithms of vascular network modeling. Algorithms for different architectures are proposed. The first algorithm is based on the message passing model and thus, it is suited for distributed memory architectures. It parallelizes the process of connecting new parts of tissue to existing vascular structures. The second algorithm is designed for shared memory machines. It also parallelizes the perfusion process, but individual processors perform calculations concerning different vascular trees. The third algorithm combines message passing and shared memory approaches providing solutions for hybrid parallel architectures. Developed algorithms are able to substantially speed up the time-demanded simulations of growth of complex vascular networks. As a result, more elaborate and precise vascular structures can be simulated in a reasonable period of time. It can also help to extend the vascular model and to test multiple sets of parameters. Secondly, a new approach in computational modeling of magnetic resonance (MR) flow imaging is proposed. The approach combines the flow computation by lattice Boltzmann method, MRI simulation by following discrete local magnetizations in time and a new magnetization transport algorithm together. Results demonstrate that such an approach is able to naturally incorporate the flow influence in MRI modeling. As a result, in the proposed model, no additional mechanism (unlike in prior works) is needed to consider flow artifacts, what implies its easy extensibility. In combination with its low computational complexity and efficient implementation, the solution is a user-friendly and manageable at different levels tool which facilitates running series of simulations with different physiological and imaging parameters. The goal of the third solution is to apply the proposed MR flow imaging model on complex vascular networks. To this aim, models of vascular networks, flow behavior and MRI are combined together. In all the model components, computations are adapted to be performed at various parallel architectures. The model potential and possibilities of simulations of flow and MRI in complex vascular structures are shown. The model aims at explaining and exploring MR image formation and appearance by the combined knowledge from many processes and systems, starting from vascular geometry, through flow patterns and ending on imaging technology.

Dans les structures en béton précontraint, les phénomènes de fluage et de retrait tendent à réduire les efforts de précontrainte initialement prévus pour maintenir le béton dans un état minimisant les forces de traction et donc la fissuration. La compréhension et la prédiction de ces phénomènes par le biais de modèles sont donc primordiales pour la conception et la maintenance à long terme des ouvrages du génie civil tels que les enceintes de confinement des centrales nucléaires.L’objectif de cette thèse est d’élaborer un cadre de modélisation micromécanique pour décrire de manière unifiée le retrait et le fluage dans les matériaux cimentaires. Pour cela, l’étude se base sur l’échelle de la microstructure poreuse du gel de C-S-H où les mécanismes intrinsèques de ces déformations différées du béton opèrent. Une approche d’homogénéisation numérique modélisant ces phénomènes dans des microstructures poreuses à morphologies quelconques est développée. Une description explicite du réseau poreux ainsi que de la phase liquide de l’eau pendant les processus de séchage/humidification est prise en compte. Les mécanismes concernant lesdéformations différées dans la phase solide sont modélisés par la théorie de la microprécontrainte-solidification (MPS). Les simulations à l’échelle microscopique sont réalisées par une approche originale couplant la méthode de Lattice Boltzmann (LBM) et la méthode des éléments finis (FEM). La LBM est utilisée pour décrire la distribution du liquide capillaire à l’échelle du pore,tandis que la FEM est employée pour simuler la déformation du squelette solide sous l’action combinée de l’eau dans l’espace poreux et d’un chargement macroscopique.La démarche proposée permet, au travers des simulations, de mieux comprendre les mécanismes liés à la non saturation et aux effets capillaires dans les milieux poreux. En particulier, la prise en compte de morphologies réalistes de microstructures et des ménisques formés conduit à différents régimes de retrait/gonflement. Ainsi les effets de l’intensité de la pression capillaire,de la tension de surface et des surfaces de chargement sur la réponse élastique du squelette solide sont évalués. Enfin, nous proposons une extension des approches précédentes au cas d’un squelette viscoélastique se déformant sous les effets de la pression capillaire et des tensions de surface. A partir des observations numériques réalisées, nous proposons un modèle pour décrire le fluage et le retrait du gel de C-S-H de manière unifiée / In pre-stressed concrete structures, creep and shrinkage tend to reduce the pre-stress forces which are initially produced so as to maintain concrete in a state minimizing traction forces and then cracks. Understanding and predicting these phenomena through models are thus highly important for the design and durability of civil engineering structures, such as containment buildings in nuclear power plants.The objective of this thesis is to develop a micromechanical modeling framework to describe shrinkage and creep in cementitious materials in a unified manner. For this purpose, the study focuses on the scale of the porous structure of the C-S-H gel where the intrinsic mechanisms of delayed strains are active. A computational homogenization approach is developed to model these phenomena in porous structures with arbitrary morphologies. An explicit description of the porous network and of the liquid phase of water during the drying/humidification process is taken into account. The mechanisms related to delayed strains in the solid phase are modeled by the microprestress-solidification theory (MPS). The simulations at the microscale are conductedbased on an original approach coupling the Lattice Boltzmann method (LBM) and the finite element method (FEM). The LBM is used to describe the distribution of capillary water in the porous structure, whereas the FEM serves as modeling the strain of the solid skeleton under the capillary water effets and a macroscopic load.The proposed method allows, by means of the simulations, to better understand the mechanisms related to the capillary effects in the porous structure. More specifically, taking into account realistic morphologies of microstructures and of the formed menisci lead to different regimes of shrinkage/swelling. Then, the effects of capillary pressure intensity, of surface tension and of morphologies of capillary surfaces on the elastic response of the solid skeleton are evaluated. Finally, the above approaches are extended to the case of a viscoelastic solid deformed under the action of the capillary water. From numerical observations, we propose a model is proposed to describe the creep and shrinkage of C-S-H gel in a unified way

L'ablation par radiofréquence (ARF) de tumeurs abdominales est rendue difficile par l’influence des vaisseaux sanguins et les variations de la conductivité thermique, compliquant la planification spécifique à un patient donné. En fournissant des outils prédictifs, les modèles biophysiques pourraient aider les cliniciens à planifier et guider efficacement la procédure. Nous introduisons un modèle mathématique détaillé des mécanismes impliqués dans l’ARF des tumeurs du foie comme la diffusion de la chaleur et la nécrose cellulaire. Il simule l’étendue de l’ablation à partir d’images médicales, d’après lesquelles des modèles personnalisés du foie, des vaisseaux visibles et des tumeurs sont segmentés. Dans cette thèse, une nouvelle approche pour résoudre ces équations basée sur la méthode de Lattice Boltzmann est introduite. Le modèle est d’abord évalué sur des données de patients qui ont subi une ARF de tumeurs du foie. Ensuite, un protocole expérimental combinant des images multi-modales, anatomiques et fonctionnelles pré- et post-opératoires, ainsi que le suivi de la température et de la puissance délivrée pendant l'intervention est présenté. Il permet une validation totale du modèle qui considère des données les plus complètes possibles. Enfin, nous estimons automatiquement des paramètres personnalisés pour mieux prédire l'étendu de l’ablation. Cette stratégie a été validée sur 7 ablations dans 3 cas cliniques. A partir de l'étude préclinique, la personnalisation est améliorée en comparant les simulations avec les mesures faites durant la procédure. Ces contributions ont abouti à des résultats prometteurs, et ouvrent de nouvelles perspectives pour planifier et guider l’ARF. / The outcome of radiofrequency ablation (RFA) of abdominal tumors is challenged by the presence of blood vessels and time-varying thermal conductivity, which make patient-specific planning extremely difficult. By providing predictive tools, biophysical models may help clinicians to plan and guide the procedure for an effective treatment. We introduce a detailed computational model of the biophysical mechanisms involved in RFA of hepatic tumors such as heat diffusion and cellular necrosis. It simulates the extent of ablated tissue based on medical images, from which patient-specific models of the liver, visible vessels and tumors are segmented. In this thesis, a new approach for solving these partial differential equations based on the Lattice Boltzmann Method is introduced. The model is first evaluated against clinical data of patients who underwent RFA of liver tumors. Then, a comprehensive pre-clinical experiment that combines multi-modal, pre- and post-operative anatomical and functional images, as well as the interventional monitoring of the temperature and delivered power is presented. This enables an end-to-end validation framework that considers the most comprehensive data set for model validation. Then, we automatically estimate patient-specific parameters to better predict the ablated tissue. This personalization strategy has been validated on 7 ablations from 3 clinical cases. From the pre-clinical study, we can go further in the personalization by comparing the simulated temperature and delivered power with the actual measurements during the procedure. These contributions have led to promising results, and open new perspectives in RFA guidance and planning.