Le thème de l’étude proposée concerne la modélisation du séchage d’un milieu poreux déformable saturé (solide-liquide). Cette modélisation tient compte de la nature du produit ainsi que de ses conditions initiales telles que la teneur en eau et la forme. L’objectif est de prévoir les contraintes mécaniques apparaissant au sein du matériau au cours du séchage, de contrôler la déformation du produit et sa teneur en eau. L’originalité de cette étude réside dans l'écriture du modèle de séchage avec la prise en compte du gradient de pression comme terme moteur réel du transport au travers de la loi de Darcy auquel s'ajoute l'hypothèse de compressibilité du liquide (loi d’état du liquide). Cette écriture présente l’avantage de ne pas recourir à un coefficient de transport effectif souvent identifié à partir d’expériences spécifiques et à terme de pouvoir relier le comportement d’un milieu supposé diphasique à un milieu triphasique (liquide, solide et gaz). Le système d’équations obtenu est issu d'une prise de moyenne volumique où le V.E.R est considérée déformable. La description physique du fort couplage hydromécanique qui existe au sein du matériau, tout au long du processus de séchage, est prise en compte au travers de la théorie de consolidation de Biot en adoptant la notion contraintes effectives de Terzaghi, les variables de couplage sont la vitesse de déformation du solide et la pression du liquide. La résolution numérique du modèle obtenu est effectuée par la méthode des éléments finis en grandes déformations et le changement d'espace Euler-Lagrange. Le modèle est validé pour un gel d’alumine à différentes conditions de séchage et l’étude de sensibilité montre la forte influence de la perméabilité du matériau comme de ses propriétés rhéologiques au cours de son séchage.

Le thème de cette étude est la modélisation du séchage d'un milieu poreux, déformable et partiellement saturé. La modélisation tient compte de la nature du produit ainsi que des conditions initiales et de sa forme. L'objectif est de prévoir les contraintes liées au retrait volumique, afin decontrôler la déformation du produit. La problématique se pose au niveau de la transition entre milieu saturé et milieu non saturé où il n'existe pas actuellement de modélisation physique continue. Dans ce travail, on propose un modèle, décrivant les transferts de chaleur, de masse et de quantité de mouvement, appliqués au séchage d'un milieu non saturé et déformable. La particularité du modèle est qu'il tient compte du fort couplage entre le transport de l'humidité et le comportement rhéologique du matériau en utilisant la notion de contrainte effective. Le modèle est validé pour une argile à différentes conditions de séchage convectif. Les simulations montrent la faisabilité du modèle proposé et l'étude de sensibilité met en évidence la forte influence de la perméabilité du matériau et de la pression capillaire d'une part, et les propriétés rhéologiques d'autre part.

Le thème de cette étude est la modélisation du séchage d’un milieu poreux, déformable, partiellement saturé (solide, liquide et gaz). Cette modélisation tient compte de la nature du produit ainsi que des conditions initiales et de sa forme. En fait, au cours de leur séchage, les milieux poreux déformables subissent des contraintes liées au retrait volumique. L’objectif est de prévoir ces contraintes afin de contrôler la déformation du produit. La modélisation du séchage a été assurée tant pour les milieux saturés que pour les milieux non saturés. La problématique se pose au niveau de la transition entre milieu saturé et milieu non saturé où il n'existe pas actuellement de modélisation physique continue. Dans ce travail, on propose un modèle, décrivant les transferts de chaleur, de masse et de quantité de mouvement, appliqué au séchage d’un milieu non saturé et déformable. Le gradient de pression est le terme moteur du transport de l’eau dans le milieu au travers de la loi de Darcy. La particularité du modèle est qu’il tient compte du fort couplage entre transport d’eau et comportement rhéologique du matériau en utilisant la notion de contrainte effective. Les variables de couplage sont la vitesse de déformation du solide et la pression intrinsèque de la phase liquide. Le modèle est validé pour une argile à différentes conditions de séchage convectif. Les simulations montrent la faisabilité du modèle décrivant le séchage d’un milieu partiellement saturé et l’étude de sensibilité met en évidence la forte influence de la perméabilité du matériau et de la pression capillaire d’une part, et des propriétés rhéologiques d’autre part.

Ces travaux de thèse portent sur la caractérisation des milieux poreux fortement déformables afin de modéliser le transfert réalisé au cours du séchage. La théorie de consolidation de Biot traduit le couplage hydromécanique se produisant durant la déshydratation d'un milieu poreux déformable. Le modèle de Cáceres s'appuie sur cette théorie et sur la loi classique de Darcy. Il est adapté aux grandes déformations du squelette solide et utilise un tenseur de contraintes de décomposition Terzaghi. Une simulation de la déshydratation par l'intermédiaire de ce modèle nécessite au préalable le renseignement de celui-ci en termes de caractéristiques thermophysiques. Deux dispositifs expérimentaux ont été mis en place : l'un mesurant la relaxation de charge, l'autre réalisant la pressurisation dynamique (DP). La 1ère expérience nous a permis de caractériser la perméabilité, le coefficient de Poisson et le module d'Young en fonction de la teneur en eau du gel. La 2ème expérience a révélé l'importance de la prise en considération de la compressibilité de la phase solide dans la modélisation du séchage du gel d'agar ainsi que, la nécessité de corriger la valeur du module de compressibilité de l'eau contenant du gaz dissout. Suite à cette étape expérimentale, nous avons concentré nos efforts sur la modélisation liée à notre étude. Tout d'abord, la DP étant conforme à l'hypothèse de faibles déformations, un modèle numérique basé sur les équations de Biot a été réalisé afin de valider les paramètres mesurés et de corriger la valeur estimée de la perméabilité. Un second modèle traduit le couplage thermo-hydro-mécanique lors du séchage convectif d'un milieu poreux fortement déformable comme les gels d'agar et d'alumine. Ce dernier s'appuie sur la théorie de Biot et d'une part s'adapte aux fortes déformations et d'autre part utilise l'approche eulérienne. Ce modèle constitue donc un compromis entre les modèles de Biot et de Cáceres et s'appuie sur l'étude thermodynamique de Coussy ce qui est une avancée dans la modélisation des gels déformables.

L'ETUDE PORTE SUR LES TRANSPORTS DE CHALEUR ET DE MASSE AYANT LIEU AU SEIN D'UN MILIEU POREUX LORS DE SON SECHAGE. A PARTIR D'UNE MODELISATION TRES COMPLETE, DES CODES DE CALCUL 1D ET 2D SONT MIS AU POINT POUR SIMULER EFFICACEMENT LES PHENOMENES INTERNES. UN NOUVEAU CHOIX D'INCONNUES CONDUIT A UNE ECRITURE DES EQUATIONS SOUS UNE FORME CONSERVATIVE. CETTE PROPRIETE ET UN SCHEMA D'INTEGRATION NUMERIQUE ADAPTE AUX NON-LINEARITES ET AU COUPLAGE DES EQUATIONS ASSURE LA PRECISION. UNE METHODE DE GESTION DU PAS DE TEMPS PROCURE UNE BONNE VITESSE DE RESOLUTION TOUT EN PRESERVANT LA QUALITE DES SOLUTIONS. CES OUTILS PERMETTENT UNE ANALYSE FINE DE LA TRANSITION ENTRE LE REGIME D'ECOULEMENT CAPILLAIRE ET LE REGIME DE DIFFUSION-SORPTION. ON MET ALORS EN EVIDENCE LES IMPERFECTIONS DU MODELE COMMUNEMENT ADMIS POUR PRENDRE EN COMPTE LE CARACTERE HYGROSCOPIQUE D'UN MATERIAU. LE CODE 1D EST ENSUITE UTILISE EN PARALLELE AVEC DES EXPERIENCES POUR PROPOSER UNE METHODOLOGIE D'ESTIMATION DES PERMEABILITES RELATIVES, PARAMETRES PHENOMENOLOGIQUES CLES DU MODELE, NON ACCESSIBLES PAR MESURE DIRECTE

Cette thèse est consacrée au développement d'une méthode numérique visant à modéliser des écoulements dans des milieux poreux fissurés, ainsi qu'à déterminer leur perméabilité effective à partir des avancements théoriques récents. En parallèle, elle a été aussi l'occasion de continuer sur la voie théorique et d'obtenir de nouveaux résultats sur ce plan.Les équations générales gouvernant l'écoulement dans de tels matériaux sont rappelées, et plus particulièrement, la conservation de la masse à l'intersection entre fissures est établie explicitement [132, 139]. Des solutions générales du potentiel sont proposées sous la forme d'une équation intégrale singulière décrivant l'écoulement dans et autour des fissures entourées par une matrice infinie soumise à un champ lointain [136, 139]. Ces solutions représentent le champ de pression dans le milieu infini en fonction de l'infiltration dans les fissures qui prennent en compte complètement l'interaction et l'intersection entre fissures. En considérant le problème d'une fissure super-conductrice, des solutions analytiques de l'écoulement ont été développées et ces solutions sont utilisées comme cas de référence pour valider la solution numérique. De plus, la solution théorique de ce problème dans le cas 3D permet de comparer le modèle d'écoulement de Poiseuille dans une fissure elliptique d'épaisseur nulle et le modèle d'inclusion ellipsoïdale aplatie soumise à l'écoulement de Darcy [140]. Des outils numériques ont été développés en se basant sur la méthode des équations intégrales singulières afin de résoudre les équations générales du potentiel [132, 180]. Cela permet, d'une part, de modéliser l'écoulement stationnaire dans un domaine poreux contenant un grand nombre de fissures et, d'autre part, de proposer une solution semi-analytique de l'infiltration dans une fissure isolée dépendant de la perméabilité de la matrice, de la conductivité de la fissure et de la variable géométrique de la fissure. Cette dernière est l'élément essentiel pour déterminer de la perméabilité effective d'un milieu poreux fissuré en utilisant des schémas d'homogénéisation. Ce modèle auto-cohérent révèle un seuil de percolation qui est alors applicable pour l'estimation de la perméabilité effective d'un matériau contenant un grand nombre de fissures. L'approche par sous-structuration permet l'extension de la solution générale du potentiel, écrite pour un domaine infini, à celle d'un domaine fini [181]. Une solution analytique de l'écoulement dans et autour d'une fissure partiellement saturée est établie, fondée sur la solution analytique pour la fissure super-conductrice. Celle-ci est alors utilisée pour estimer la perméabilité effective d'un milieu poreux fissuré non-saturé [141]. Le modèle de la perméabilité effective est appliqué dans le contexte du stockage géologique du CO2 en vue d'étudier le comportement d'une zone de faille constituée par un noyau argileux et des zones fissurées. La pression d'injection provoque l'augmentation de la pression interstitielle dans le réservoir. Cette surpression engendrée dans le réservoir peut affecter la perméabilité de zones fissurées ce qui conduit à des phénomènes hydromécaniques couplés. Les résultats de simulations numériques permettent d'évaluer, d'une part, le risque de la remontée de la saumure à l'aquifère supérieur, et d'autre part, le risque de l'initiation d'une rupture sur le plan de la faille.