Les solutions exactes autosemblables de l'hydrodynamique, avec conduction thermique non-linéaire, pour des gaz parfaits en milieu semi-infini, constituent une approche possible pour l'étude de la stabilité d'écoulements en fusion par confinement inertiel. Ces solutions et leurs perturbations linéaires, calculées 'a l'aide d'une méthode spectrale de Chebyshev multidomaine, permettent de rendre compte, sans autre approximation, du caractère compressible et instationnaire de ces écoulements. Suite aux résultats (Clarisse et al., 2008; Lombard, 2008) représentatifs de l'ablation précoce d'une cible par un flux laser non-uniforme (conduction électronique, front thermique subsonique précédé d'un choc quasi-parfait), nous explorons ici d'autres configurations. Pour cette même phase d'ablation, mais avec un rayonnement X incident non-uniforme (conduction radiative), nous traitons un écoulement compressible et un autre faiblement compressible. Dans les deux cas, nous retrouvons les comportements des écoulements compressibles obtenus en conduction électronique avec une instabilité maximale pour un nombre d'onde nul. D'autre part, la méthode spectrale est étendue au calcul de solutions autosemblables en tenant compte de l'onde thermique supersonique en amont du choc. Basée sur une analyse des singularités des équations réduites (front infiniment raide), cette méthode permet d'accéder au régime d'onde thermique supersonique propre à l'irradiation initiale d'une cible et de retrouver les solutions ablatives calculées antérieurement dans l'approximation de précurseur thermique négligeable.

L'objectif de cette étude est d'analyser l'écoulement traversant le récepteur solaire haute température de la centrale PEGASE. Cet écoulement est turbulent et soumis à un très fort gradient de température. Nous utilisons la simulation des grandes échelles thermiques appliquée aux équations bas-Mach pour étudier l'écoulement turbulent dans une géométrie de canal plan bi-périodique avec températures imposées aux parois. Pour cette étude, nous avons réalisé des simulations pour deux nombres de Reynolds turbulents différents (Re_tau m = 180 et Re _taum = 395) et pour quatre rapports de température (T2=T1 = 1 ; 1; 04 ; 2 et 5). Après avoir validé notre modèle, nous avons réalisé une étude sur la modélisation sous-maille thermique qui conclut sur la nécessité d'utiliser un modèle sous-maille thermique dynamique pour un écoulement ayant une forte intensité turbulente soumis à un très fort gradient de température. Nous analysons l'influence du gradient de température sur les différents profils de vitesse et de température ainsi que sur les spectres d'énergie. Augmenter le gradient de température crée une dissymétrie de tous les profils. Une relaminarisation du côté chaud est visible. On remarque une nouvelle répartition des fluctuations à travers le canal ainsi que de l'énergie en fonction de la taille des échelles. Une création de fluctuations de vitesse et de température a également été mise évidence du côté froid et du côté chaud du canal. Ces effets sont dus, à la fois à une variation de la viscosité cinématique (effet direct de la température) et à une interaction entre les champs turbulents dynamique et thermique.

Le but de cette thèse est le développement d'un algorithme adaptatif pour la résolution des équations aux dérivées partielles paraboliques en géométrie cartésienne pour des problèmes en dimension un, deux et trois et l'application aux instabilités de flamme dans l'approximation thermodiffusive. Partant d'un schéma de discrétisation de type volumes finis explicite, nous appliquons une décomposition adaptative multi-résolution pour représenter la solution sur un maillage localement raffiné. Les flux numériques sont calculés directement sur la grille adaptative. Afin de suivre l'évolution de la solution au cours du temps, nous utilisons une stratégie d'adaptation dynamique basée sur la représentation des données en multi-résolution. Cette dernière consiste à représenter la solution à l'aide des valeurs sur une grille grossière, plus l'ensemble des différences de prédiction entre les valeurs d'une grille donnée et celles d'une grille plus fine, l'ensemble constistuant une hiérarchie de grilles emboîtées. La compression des données s'obtient en supprimant les différences inférieures à une certaine tolérance fixée. Nous validons cette méthode par la résolution numérique d'équations de référence, comme l'équation de convection-diffusion ou l'équation de Burgers diffusive, afin de montrer la précision de la méthode et son efficacité par rapport au même schéma volumes finis sur grille fine. En particulier, pour l'équation linéaire de convection-diffusion, nous donnons une relation entre la tolérance et le nombre d'échelles qui permet de réduire le temps de calcul et la place mémoire nécessaires tout en maintenant l'ordre de précision du schéma volumes finis. Ce résultat est confirmé par le calcul numérique. Nous utilisons ensuite la méthode adaptative pour étudier les instabilités de flammes pré-mélangées dans l'approximation thermodiffusive. En particulier, pour les flammes planes, nous déterminons la limite d'apparition des flammes pulsantes, limite que nous comparons aux données de la littérature, tant numériques que théoriques. Nos calculs confirment la théorie pour les grandes valeurs de l'énergie d'activation. Nous montrons également numériquement l'existence de flammes planes non-pulsantes lorsque la conduction de la chaleur est très supérieure à la diffusion des réactants. Nous étudions également les ballons de flamme et montrons que, lorsqu'ils interagissent avec une paroi adiabatique, leur comportement présente une analogie avec la capillarité en mécanique des fluides. Le dernier résultat concerne l'interaction d'un ballon de flamme avec un tourbillon. Il ouvre des perspectives sur la simulation adaptative des ballons de flamme dans un fluide en mouvement.

Le transfert d'eau à l'interface sol-atmosphère est abordée par une approche de thermodynamique qui constitue un cadre général pour décrire à la fois l'état thermodynamique de l'eau dans le sol et les mécanismes de transfert mis en jeu : filtration de la phase liquide, diffusion de la vapeur d'eau, changement de phase liquide vapeur de l'eau. Une caractérisation expérimentale d'un sol a été réalisée puis complétée par des données bibliographiques. Des essais expérimentaux de transfert d'eau dans des colonnes du même sol, placées dans une atmosphère asséchante a permis d'établir la cinétique globale de perte d'eau ainsi que l'évolution des profils de teneur en eau au cours du temps. Une simulation numérique du transfert d'eau au voisinage de l'interface est proposée. Une étude de sensibilité aux paramètres physiques du sol est présentée. Les résultats du modèle numérique sont ensuite confrontés aux essais de transfert sur colonnes. La correspondance est bonne moyennant un ajustement des paramètres hydrologiques issus de la littérature. L'évolution des flux d'eau liquide et de vapeur permettent de décrire les mécanismes au voisinage de l'interface dans le cas d'un sol à faible teneur en eau. On montre que le changement de phase qui assure la continuité des transferts entre le sol et l'atmosphère concerne quelques centimètres du sol sous cette interface.

Les nouveaux bâtiments sont sujets à un renforcement de l'isolation thermique, à une augmentation des surfaces vitrées et à un accroissement des apports internes dus à des appareils électriques de plus en plus nombreux (ordinateurs, électroménager...). Cette tendance conduit à une diminution des demandes de chauffage mais peut engendrer des demandes de rafraîchissement en mi-saison et en période estivale. Les besoins en chauffage et en rafraîchissement deviennent plus équilibrés sur une journée et sur l'année en général. De plus, l'eau chaude sanitaire (ECS) occupe une place de moins en moins négligeable dans les besoins énergétiques totaux, ce qui provoque une simultanéité partielle des besoins en chaud et en froid dès que des demandes de rafraîchissement apparaissent. Une thermofrigopompe (TFP), machine frigorifique produisant de l'énergie calorifique et frigorifique de manière simultanée, associée à un stockage tampon d'énergie devient alors une solution intéressante. La TFP de petite à moyenne puissance est un système de production d'eau chaude et d'eau froide utilisant l'air comme source gratuite, destiné au chauffage, au rafraîchissement et à la production d'ECS pour des bâtiments du secteur résidentiel et du petit tertiaire. Elle présente plusieurs spécificités techniques dans la conception du circuit frigorifique, dans la gestion des modes de fonctionnement et des séquences de dégivrage. L'ajustement des productions de chaud et de froid aux charges du bâtiment s'effectue grâce à l'utilisation d'un échangeur d'équilibrage sur air extérieur fonctionnant soit en condenseur, soit en évaporateur. En période de chauffage, cette TFP offre également la possibilité de stocker une certaine quantité d'énergie sur la boucle d'eau froide à l'aide d'un sousrefroidisseur. Cette énergie stockée est utilisée en temps différé à l'évaporateur à eau afin d'améliorer les performances par un relèvement de la température d'évaporation et éventuellement de dégivrer l'évaporateur à air sans arrêter la production de chaleur. L'énergie de dégivrage est apportée par un thermosiphon diphasique formé entre les deux évaporateurs à des températures différentes. La présence d'un sous-refroidisseur impose le contrôle de la haute pression pour assurer une condensation totale du réfrigérant dans le condenseur utile. Le système de contrôle de pression est constitué d'une bouteille contenant du réfrigérant à l'équilibre liquide / vapeur connectée en partie basse à la ligne liquide et en partie haute au refoulement du compresseur et à un point en basse pression du circuit. Une injection de gaz provenant du refoulement du compresseur entraîne une augmentation de la haute pression. Une chasse vers la basse pression provoque une diminution de la haute pression du système. Une étude expérimentale a été menée sur une machine prototype de TFP au R407C pour des conditions variables de fonctionnement en températures d'air et d'eau et pour les différents modes possibles (chauffage, rafraîchissement et production simultanée). Les performances obtenues expérimentalement sont en accord avec celles annoncées par le logiciel de sélection développé par le constructeur du compresseur. L'étude expérimentale a également permis d'observer le dégivrage par thermosiphon, de vérifier le bon fonctionnement du système de gestion de la pression de condensation et des séquences alternées en période hivernale (1 : mode chauffage avec stockage de chaleur sur la boucle d'eau froide par sous-refroidissement. 2 : mode équilibré utilisant une évaporation par l'énergie stockée sur la boucle d'eau froide, plus performant, avec dégivrage éventuel). Le fonctionnement de la TFP a été modélisé avec deux réfrigérants : le HFC R407C et le CO2. Dans le contexte actuel animé par d'éventuelles décisions politiques concernant l'impact des gaz à effet de serre sur le réchauffement climatique, le dioxyde de carbone est un fluide intéressant d'abord pour son faible impact environnemental (ODP nul et GWP100ans = 1) et pour les caractéristiques particulières du cycle thermodynamique utilisé (cycle transcritique). En effet, une grande quantité d'énergie est récupérable par sous-refroidissement du fluide frigorigène et de l'ECS à haute température peut être produite de manière performante. Dans des simulations annuelles, les TFP au R407C et au CO2 sont comparées à des pompes à chaleur réversibles sur les critères de performance, d'économie d'énergie et d'impact environnemental. Le concept de cette TFP offre clairement une amélioration des performances par rapport à une pompe à chaleur réversible classique et l'utilisation du dioxyde de carbone présente de bonnes perspectives de développement.

Le transport de particules solides dans un matériau poreux saturé, sujet à l'écoulement de fluide se produit dans les milieux naturels et dans diverses applications liées au génie civil. Ces particules solides peuvent colmater les interstices du matériau poreux et par conséquent réduire sa perméabilité ou au contraire se détacher du matériau poreux et provoquer son érosion. Ainsi, cette thèse s'attache à étudier à l'échelle des grains du matériau poreux les mécanismes du transport de particules en présence d'écoulement de fluide incompressible et de calculer les interactions fluide-solide pour contribuer à la connaissance des phénomènes de colmatage ou d'érosion. Le travail est réalisé principalement par la voie numérique utilisant un couplage original entre deux codes de calcul : PFC3D pour la matrice granulaire - CFD FINAS pour l'écoulement du fluide, validé par des expériences de visualisation menées sur des modèles de laboratoire.<br/> Les études numérique et expérimentale de la filtration dans un milieu poreux bidimensionnel montrent l'influence des particules sur le champ de vitesse du fluide qui se met à son tour en mouvement et accélère la filtration. L'étude des interactions particules-fluide par l'approche couplée DEM-CFD (Distinct Element Method-Computational Fluid Dynamics) a abouti à proposer une formulation de la traînée en fonction de la porosité du milieu poreux. Le couplage avec la loi de Darcy conduit à une relation pratique qui prédit la perméabilité en fonction de la porosité. L'étude est étendue à la filtration de particules solides dans un milieu poreux constitué d'un empilement désordonné de sphères dans lequel un fluide s'écoule en permanence. L'analyse à l'échelle des grains montre des zones de colmatage qui se développent dans le milieu poreux. L'étendue du colmatage est analysée en fonction de combinaisons des forces prédominantes que sont les forces de gravité et de viscosité. Des courbes de distribution des particules selon la profondeur sont trouvées.

Cette thèse vise à une meilleure compréhension de la dynamique du voilier et à la validation des outils numériques de prédiction de performances et d'optimisation par l'étude expérimentale in situ du problème aéro-élastique d'un gréement. Une instrumentation est développée sur un voilier de 8m de type J80 pour la mesure dynamique des efforts dans le gréement, de la forme des voiles en navigation, du vent et des attitudes du bateau. Un effort particulier est apporté à la mesure des caractéristiques géométriques et mécaniques des éléments du gréement, la calibration des capteurs et au système d'acquisition des données. Les principaux résultats montrent que le voilier instrumenté est un outil adapté pour les mesures instationnaires et soulignent l'amplitude de variation d'effort rencontrée en mer (20 à 50% de l'effort moyen dans une houle modérée). En outre, les variations du signal d'effort sont déphasées avec l'angle d'assiette, créant un phénomène d'hystérésis. Le comportement dynamique d'un voilier en mouvement diffère ainsi de l'approche quasi-statique. Les simulations numériques proviennent du code ARAVANTI, couplage implicite d'un code structure éléments finis ARA et d'un code fluide parfait, limitant son domaine de validité aux allures de près Les résultats de simulation sont très proches des cas stationnaires et concordent bien avec les mesures en instationnaire dans une houle de face. L'expérimentation numérique d'un gréement soumis à des oscillations harmoniques en tangage souligne l'importance de l'approche Interaction Fluide Structure (IFS) et montre que l'énergie échangée par le système avec la houle est reliée à la fréquence réduite et l'amplitude du mouvement. Certaines informations n'étant pas disponibles sur le voilier instrumenté, une expérience contrôlée en laboratoire est développée. Elle consiste en un carré de tissu tenu par deux lattes en oscillation forcée. Les mesures sur cette " voile oscillante " permettent d'étudier les phénomènes IFS avec décollement et sont utilisées pour la validation du couplage ARA-ISIS entre un code fluide Navier-Stokes (RANS) et le même code structure.

Les transferts d'humidité dans les matériaux de construction ont une influence importante sur leur durabilité et sur les performances hygriques et thermiques des bâtiments. De nombreux mécanismes d'endommagement chimiques et physiques de ces matériaux sont en effet dus à l'infiltration d'eau. En conséquence, leur structure poreuse peut évoluer au cours du temps, et des fissures microscopiques comme macroscopiques peuvent s'y développer. La description des matériaux à l'échelle microscopique est cependant une source d'erreur importante dans les codes de simulation actuels des transferts d'humidité et de chaleur, notamment en raison du fait que les milieux sont considérés comme homogènes, et que les effets du vieillissement des matériaux sont négligés. Il importe donc de trouver un moyen d'inclure les effets de l'endommagement dans les simulations de transferts d'humidité et de chaleur à l'échelle du bâtiment. Des méthodes existent pour la prédiction du comportement de milieux soumis à des sollicitations hygriques et mécaniques, mais supposent que l'ensemble des facteurs extérieurs influant sur l'endommagement soient connus tout au long des simulations. Une nouvelle méthodologie a été proposée pour compléter ces approches prédictives, en combinant des mesures expérimentales d'endommagement avec la simulation de transferts couplés d'humidité et de chaleur. Une étude préliminaire a d'abord été menée, consistant à mesurer la perméabilité vapeur équivalente d'éprouvettes de mortier multi-fissurées. Cette démarche a permis d'identifier les besoins expérimentaux et numériques de la suite du travail, visant à modéliser les écoulements dans un réseau discret de fissures sur la base de leur caractérisation. Une méthodologie expérimentale combinant corrélation d'images numériques et émissions acoustiques a ensuite été développée, permettant de disposer de cartographies d'endommagement et de proposer une démarche pour la mesure de réseaux de fissures dans les matériaux de construction en place. La méthode optique, associée à une procédure de traitement d'images, a permis de disposer de données précises de la géométrie de réseaux de fissures. De plus, une méthode a été proposée pour permettre l'interprétation des mesures d'émissions acoustiques en termes de quantification, localisation et identification des phénomènes d'endommagement. Un code de simulation a ensuite été écrit, permettant d'intégrer ces mesures de fissuration dans la modélisation des écoulements couplés d'humidité et de chaleur en milieu poreux. Ce modèle a été validé sur la base de mesures expérimentales : la corrélation d'images numériques a été appliquée durant la fracturation d'éprouvettes de béton, dans lesquelles l'infiltration d'eau a ensuite été suivie par radiographie aux rayons X. Les résultats numériques obtenus sont en bonne conformité avec les mesures expérimentales en termes de prédiction de la concentration d'eau en deux dimensions. Enfin, la méthodologie a été appliquée à une série de cas test, dans le but de modéliser les performances hygrothermiques de parois multi-couches, incluant des matériaux endommagés, soumises à des conditions climatiques réelles. On a ainsi pu estimer les conséquences potentielles de l'endommagement sur l'accumulation d'eau dans des parois, sur l'amplitude des cycles de sorption et de séchage, ainsi que sur les transferts thermiques.

Cette thèse est scindée en deux parties. Dans un premier temps, nous présentons deux schémas volumes finis pour la discrétisation des problèmes de diffusion-convection-réaction sur les surfaces mouvantes. Le premier schéma présente une extension du schéma volumes finis avec flux à deux points sur les surfaces mouvantes. Le deuxième développe une méthode de type O-méthode. Cette dernière consiste à construire à partir des inconnus au centre des mailles, des solutions linéaires autour des nœuds de maillage qui intègrent la continuité des flux aux interfaces de mailles. La méthode permet aussi la construction des décentrages amont d'ordre 2 et ainsi, offre au procédé l'ordre 2 de convergence sur tout maillage non dégénéré. Ensuite, nous modélisons l'écoulement du couplage filme mince-surfactant (surface active agent) sur les surfaces mouvantes et simulons à l'aide des schémas volumes finis précédemment définis. Ici, l'utilisation du calcul tensoriel et de la méthode d'approximation par la lubrification permettent de réduire les équations de Navier-Stokes caractérisant le mouvement du filme mince en dimension 3 en un système d'équations définies sur la surface courbe mouvante dont l'inconnu est la hauteur du fluide. Le surfactant supposé insoluble est modélisé par une équation de diffusion convection à la l'interface fluid-air. Nous simulons l'ensemble en utilisant une méthode dite de capture d'interface (Interface tracking method) dérivée des volumes finis définis plus haut. Plusieurs exemples illustrent à suffisance l'efficacité et la précision des différentes méthodes.

Durant mon parcours en tant que Maître de Conférences depuis 2005, mes diverses activités de recherche se sont situées dans le domaine de la mécanique des fluides et des transferts. Si l'approche dominante dans mes recherches est la modélisation et la simulation numérique, dans les cas où l'expérimentation s'avère l'approche la plus pertinente, je m'associe à mes collègues expérimentateurs pour mener à bien des études. En effet, dans ma démarche de chercheur, j'ai toujours souhaité contribuer à résoudre des problèmes scientifiques et technologiques issus pour la plupart de préoccupations industrielles ou se situant directement en amont de celles-ci. La principale thématique que j'ai abordée est celle du mélange laminaire. Que ce soit pour l'agitation chaotique ou pour l'émulsification, j'essaye, à travers mes travaux, de promouvoir ce mode de mélange sobre, afin d'encourager son usage industriel. L'étude du mélange d'un ou de plusieurs fluides avec ou sans transfert de chaleur est toujours d'un grand intérêt pour les industriels, surtout lorsqu'il s'agit de fluides très visqueux et/ou très délicats pour lesquels les méthodes classiques de mélange turbulent peuvent s'avérer coûteuses en énergie et/ou destructrices des propriétés physiques de ces fluides, ou tout simplement inefficaces. Le mélange par advection chaotique a montré son efficacité pour le mélange d'un champ scalaire (concentration, température) pour des fluides très visqueux ou en micro-fluidique. En effet, il est possible à partir d'un écoulement simple d'obtenir des trajectoires lagrangiennes complexes, en apportant une modulation temporelle adéquate, par exemple. Nous nous sommes intéressé au mélange thermique (à l'échelle macroscopique) au sein d'un fluide très visqueux à grand nombre de Prandtl. Afin de contribuer à l'optimisation de ce procédé, nous avons modélisé l'écoulement instationnaire dans d'un mélangeur type et nous avons révélé les différents mécanismes de mélange thermique qui interviennent pour des fluides newtoniens ou non, avec ou sans thermodépendance. En outre, j'entreprends le développement d'un solveur numérique de mécanique des fluides qui résout les équations tridimensionnelles (3D) de conservation de la masse, de la quantité de mouvement et de l'énergie. Ce code, que j'ai baptisé Tamaris, utilise la méthode des volumes finis non structurés à maillage hybride pouvant contenir différents types de cellules de calcul (hexaèdres, tétraèdres, prismes et pyramides), ce qui permet une grande flexibilité pour les géométries étudiées. Ce code permet l'étude des écoulements instationnaires de fluides incompressibles newtoniens ou non newtoniens, en prenant en compte les transferts thermiques par convection forcée ou naturelle, les effets thermocapillaires, les transferts conjugués (conduction dans les solides) et les changements de phase liquide-solide. Ce code est parallèle (MPI) par décomposition de domaine sans chevauchement. Parmi les autres thématiques que j'ai abordées, il y a l'émulsification laminaire. Les études expérimentales menées ont montré qu'il est possible de réaliser des émulsions très concentrées en phase dispersée (> 90%), grâce à un procédé sobre énergétiquement, qui utilise un écoulement laminaire. Ce procédé a été intensifié, en le transposant vers un fonctionnement en continu, très souhaitable pour un usage industriel. Une autre de mes thématiques concerne les changements de phase liquide-solide, principalement en présence d'un écoulement. L'application principale concerne le contrôle de la température de parois et le transport d'énergie par chaleur latente. Aujourd'hui, ces thématiques d'émulsification laminaire et de changement de phase se rejoignent dans un nouveau projet pour la conception de nouveaux matériaux pour l'isolation thermique active. C'est en ce sens que mes projets de recherche futurs s'orientent naturellement vers la modélisation numérique des écoulements avec interfaces : écoulements diphasiques et écoulements avec changement de phase.