Deux approches pour caractériser les milieux poreux en terme de distribution de taille de pores (DTP) sont développées au sein de l'équipe ECPS. Ce travail a pour but de confirmer expérimentalement leurs validités. A l'instar des autres méthodes utilisant l'intrusion du mercure, l'adsorption isotherme ou la thermoporosimétrie, la première méthode consiste à utiliser un fluide à seuil d'écoulement. En effet, l'utilisation de l'écoulement d'un fluide à seuil de type Herschel-Bulkley, au travers d'un poreux, en fonction du gradient de pression permet (en utilisant les solutions analytique et numérique du problème inverse) de déterminer la fonction de distribution de la taille de pores. La seconde méthode utilise l'admittance complexe d'un milieu poreux, mesurée à partir de la réponse en débit à une sollicitation harmonique du gradient de pression. Comme la fréquence de la sollicitation est reliée aux rayons des pores par le biais de la profondeur de pénétration hydrodynamique, l'admittance permet de retrouver la distribution de taille de pores par la résolution numérique du problème inverse associé. Ces deux techniques sont basées sur le modèle de faisceaux de capillaires parallèles employé dans la plupart des autres études qui traitent du même problème. Nos expériences s'appuient sur des milieux poreux calibrés. L'application de ces techniques aux milieux poreux réels se fait actuellement en collaboration avec le TREFLE de Bordeaux. Les résultats expérimentaux obtenus affirment clairement la validité et l'applicabilité de ces deux méthodes pour la caractérisation de la DTP. Il est désormais envisageable de les transférer pour un usage industriel.

L'objectif de ce travail est de contribuer à une meilleure prédiction de la rupture des Anévrismes de l'Aorte Abdominale. Pour ce faire, des simulations par éléments finis ont été mises en oeuvre sur des anévrismes modèles dans des conditions proches de la réalité physiopathologique, i.e. en tenant compte de l'anisotropie de la paroi anévrismale, du caractère poreux du thrombus et des Interactions Fluide-Structure. Dans la première partie, une étude statique en l'absence du thrombus a permis de mettre en évidence l'influence de la géométrie de l'anévrisme et du comportement mécanique (isotrope ou anisotrope) de la paroi sur la distribution des contraintes, i.e. la rupture potentielle de l'anévrisme, ainsi que sur l'évolution du module de Peterson. Dans la seconde partie, une modélisation poro-hyperélastique du thrombus est proposée, en s'appuyant sur des données expérimentales de la littérature. La prise en compte de ce comportement et des Interactions Fluide-Structure montre que la pression intra-thrombus demeure du même ordre de grandeur que la pression intra-luminale, conformément à des mesures in vivo réalisées par ailleurs. Enfin, nous montrons que ceci n'est pas en contradiction avec une réduction du risque de rupture potentielle de l'anévrisme.

La stabilité gravitaire des ouvrages hydrauliques (digues fluviales, barrage en remblai, ...) est un problème complexifié par la présence d'infiltration d'eau au sein même de l'ouvrage et de ses fondations. Le territoire national compte près de 10000 km de digues, d'où un enjeu important en terme de risque d'inondation. Les déstabilisations gravitaires observées sur les ouvrages hydrauliques se distinguent par leur brièveté et par leur faible emprise spatiale. L'objectif de ce travail de thèse est d'étudier spécifiquement ces instabilités brèves et de faible déplacement cumulé le long de la pente dans le cas de matériaux granulaires saturés. Dans un premier temps, une campagne expérimentale a été effectuée sur la base d'un protocole de chargement progressif dans lequel un échantillon de sol est lentement incliné jusqu'à obtenir une déstabilisation massive. Il en ressort une phénoménologie complexe composée d'un nombre de précurseurs d'avalanche dépendant de la fraction volumique de l'échantillon. Une forte influence de la fraction volumique du sol sur le seuil de stabilité est également mise en avant. Une deuxième campagne expérimentale basée sur un protocole d'effondrement contrôlé, à inclinaison imposée, a été menée dans le but d'amplifier la cinématique de déstabilisation observée en chargement progressif. Mis à part les échantillons les plus denses qui produisent de simples avalanches de surface, on observe dans les autres cas un ou plusieurs événements successifs, d'allure circulaire, mobilisant des profondeurs très importantes contrairement aux avalanches de surface.. Outre l'inclinaison, on voit bien que, là encore, la fraction volumique est un paramètre clé dans le contrôle de la cinématique de déstabilisation. Des essais préliminaires laissent cependant entendre que, au-delà de la fraction volumique, c'est la microstructure qui semble être le vrai paramètre de contrôle à travers notamment la forme des grains et leur enchevêtrement.

Mes recherches ont porté sur l'analyse des équations à double non linéarité, le transport neutronique et la mécanique des textiles, et plus récemment sur la méthode Level Set et ses applications au couplage fluide-structure, notamment dans le domaine biomécanique.

Cette thèse est consacrée à l'analyse mathématique de trois problèmes d'écoulements de fluides viscoélastiques de type Oldroyd. Tout d'abord, nous étudions des écoulements stationnaires faiblement compressibles dans un domaine borné avec des conditions au bord de type "rentrante-sortante". Nous étudions aussi le problème d'écoulements stationnaires faiblement compressibles dans un coin convexe. En utilisant une méthode de point fixe (premier et deuxième problèmes) et une décomposition de Helmoltz (deuxième problème), nous montrons des résultats d'existence et d'unicité des solutions. Nous étudions également le cas d'un écoulement non stationnaire. Nous montrons un résultat d'existence locale et un résultat d'existence globale, avec des conditions initiales suffisamment petites, pour des fluides compressibles. Nous démontrons aussi la convergence du modèle d'écoulement viscoélastique compressible à faible nombre de Mach vers le modèle incompressible lorsque les données initiales sont "bien préparées"

Le risque individuel de rupture des anévrismes cérébraux est un enjeu majeur dans la prise en charge clinique des anévrismes asymptomatiques. La rupture anévrismale se produit lorsque la contrainte intra-pariétale dépasse la contrainte à rupture du matériau composant la paroi. Notre étude a pour objectif d'être un pas vers une nouvelle mesure biomécanique du risque individuel de rupture des anévrismes cérébraux. Dans un premier temps, une étude expérimentale fût menée pour caractériser le comportement biomécanique de la paroi anévrismale sur 16 échantillons d'anévrismes prélevés chirurgicalement. L'expérimentation sur les échantillons de poche anévrismale a permis de dégager trois grandes classes de tissus pour chaque sexe (homme et femme) : souple, rigide et intermédiaire. Il apparaît que tous les anévrismes non rompus appartiennent à la catégorie rigide ou intermédiaire et que tous les anévrismes rompus correspondent à la catégorie souple. Ceci permet de mettre en évidence une corrélation entre le risque de rupture et les propriétés du matériau composant la paroi anévrismale. Dans un deuxième temps, des simulations d'interaction fluide/structure (FSI) ont été réalisées pour comparer les déformations d'un anévrisme " patient spécifique " constitué d'un matériau dégradé et non dégradé. Les résultats montrent que les propriétés du matériaux ont un impact majeur sur l'ampleur de la variation de volume anévrismale diastolo-systolique. Les changements en terme de variations de volume en fonction des caractéristiques du tissu sont potentiellement visualisable à l'aide de l'imagerie médicale. Une analyse des incertitudes des paramètres est aussi présentée et montre la robustesse des résultats aux incertitudes des données d'entrée. Il a ensuite été démontré sur 12 cas " patient-spécifique " d'anévrismes différents (forme, taille, localisation et conditions aux limites différentes) qu'il existe toujours une différence significative en terme de variation de volume au cours du cycle cardiaque entre un anévrisme dont la paroi est composé d'un matériau rigide et d'un matériau souple. Cette étude suggère donc que la variation de volume anévrismale pourrait être utilisée comme une base pour une évaluation individuelle du risque de rupture des anévrismes cérébraux.

Le comportement des matériaux multiphasiques couvre une multitude de phénomènes qui suscitent un grand intérêt dans le domaine scientifique et professionnel. Les propriétés mécaniques de ces types de matériau trouvent leur origine dans les phases dont ils sont composés, leur distribution et interaction. Un nouveau modèle hydrodynamique couplé est présenté dans ce travail de thèse, à appliquer à l'analyse de l'hydrodynamique des milieux granulaires saturés. Le modèle associe la méthode des éléments discrets (DEM) pour la modélisation de la phase solide, avec une formulation en volumes finis, à l'échelle des pores (PFV), du problème de l'écoulement. Une importance particulière est donné à la description de l'interaction entre les phases, avec la détermination des forces fluides à appliquer sur chacune des particule, tout en assurant un coût de calcul abordable, qui permet la modélisation de plusieurs milliers des particules en trois dimensions. Le milieux est considéré saturé par un fluide incompressible. Les pores et leur connectivité est basée sur une triangulation régulière des assemblages. L'analogie de cette formulation avec la théorie classique de Biot est présenté. Le modèle est validé par la comparaison des résultats numériques obtenus pour un problème de consolidation d'un sol granulaire avec la solution analytique de Terzaghi. Une approche pour analyser l'hydrodynamique d'un sédiment granulaire est finalement présenté. La reproduction du phénomène de liquéfaction d'un sol est également présentée.

Le Syndrome d'Apnées Obstructives du Sommeil est caractérisé<br />par la survenue fréquente d' épisodes d' obstruction<br />des voies aériennes supérieures. L'intérêt d'une modélisation physique<br />est qu'elle autorise une compréhension<br />plus fine du phénomène, et laisse espérer une amélioration<br />des traitements. <br />Le but a donc été de concevoir, puis<br />de valider, un algorithme de simulation numérique de<br />l'interaction entre les tissus vivants et le flux<br />d'air à l'origine d'un épisode apnéique. Afin d'alléger<br />les calculs et de réduire le temps de simulation, des hypothèses<br />simplificatrices ont été envisagées. D'une part, en ce qui concerne<br />les tissus vivants, la méthode des éléments finis permet une<br />prédiction réaliste de leur déformation. Le cadre<br />des petites perturbations et de l'élasticité linéaire implique<br />de plus un calcul rapide de la réponse mécanique. <br />D'autre part, la simulation<br />de l'écoulement d'air se fait via une formulation<br />asymptotique des équations de Navier-Stokes (équations<br />de Navier-Stokes Réduites / Prandtl), qui facilite la résolution<br />numérique.<br />Afin de valider hypothèses physiques et méthode de résolution<br />numérique, une maquette in-vitro a été utilisée. Celle-ci<br />permet de reproduire, dans des conditions contrôlées, une<br />interaction entre flux d'air et paroi déformable analogue<br />à celle qui se produit à la base de la langue en début d'obstruction.<br />Une mesure précise de la déformation du conduit d'écoulement<br />est obtenue à l'aide d'une caméra digitale.<br />Une série de comparaisons quantitatives <br />a montré qu'en dépit des simplifications effectuées, l'erreur entre prédiction<br />et mesures est faible.<br />Finalement, pour se rapprocher de la réalité clinique,<br />des modèles de voies aériennes supérieures de quatre patients apnéiques ont<br />été construits à partir de radiographies sagittales.<br />Des comparaisons entre simulations à partir de radiographies<br />pré-opératoires et post-opératoires ont montré<br />que les prédictions étaient globalement cohérentes avec les conséquences<br />du geste chirurgical. Elles ont pû également<br />mettre en évidence certaines limites de notre approche,<br />dûes à la complexité du phénomène.

Ce travail de thèse concerne la prédiction des structures en béton soumises à des chargements extrêmes. Il s'intéresse plus particulièrement au comportement du béton sous fort confinement où la contrainte peut atteindre des niveaux de l'ordre du giga Pascal. La modélisation de ce comportement doit être capable de reproduire la compaction irréversible. Pour ce faire, deux lois de comportement élasto-plastique - endommageable ont été développées et implantées dans un code aux éléments discrets. Les paramètres utilisés dans ces lois sont calibrés par les simulations des essais de traction/compression uniaxial, des essais hydrostatiques et triaxiaux. Une fois les paramètres calibrés, la loi montrant le meilleur agrément avec l'expérience a été choisie pour la prédiction de la réponse du béton sous différents niveaux de confinement. Les résultats du modèle sont analysés non seulement à l'échelle macroscopique mais également à l'échelle de l'élément discret. La nécessité de prendre en compte une loi d'interaction de type élasto-plastique-endommageable est aussi montrée. La deuxième partie du travail de thèse développe une méthode de couplage entre le modèle éléments discrets et un modèle d'écoulement compressible en tenant compte des mécanismes physiques fondamentaux d'interaction entre l'écoulement interne et les particules solides d'un matériau poreux. Le problème d'écoulement est résolu par une méthode en volumes finis, où le volume est discrétisé en tétraèdres issus d'une triangulation régulière de Delaunay. Notre modèle est une adaptation aux fluides compressibles d'un modèle développé initialement pour les écoulements incompressibles. Ce couplage a été utilisé pour simuler le comportement triaxial des bétons humides et saturés sous différents niveaux de confinement. Les résultats nous montrent une bonne reproduction du comportement non-drainé du béton saturé sous faible confinement. Pour fort confinement, les simulations ne se rapprochent des résultats expérimentaux qu'au prix d'une compressibilité du fluide plus faible que celle de l'eau. Par ailleurs, la contrainte effective était une variable pertinente pour décrire le comportement du béton humide par un état limite intrinsèque indépendant du degré de saturation.

Contrairement aux liquides ordinaires, les fluides complexes comme le sang exhibent des comportements étranges qui dépendent essentiellement des structures sous-jacentes qui les composent. La simulation des écoulements sanguins continue de poser un formidable défi pour les modélisations théoriques et numériques dont l'intérêt est de développer des méthodes et des outils de simulation pour la communauté médicale. Nous proposons dans cette thèse une contribution à ce projet qui sera majoritairement centré sur les aspects numériques et informatiques. Nous nous sommes particulièrement intéressés à l'interaction entre le sang et la paroi vasculaire, qui joue un rôle important dans les grandes artères comme l'aorte. Nous nous sommes aussi investis dans la simulation du transport des cellules sanguines dans le sang. Pour la résolution des équations aux dérivées partielles décrivant nos modèles d'hémodynamique, nous avons choisi d'utiliser des méthodes numériques dont la précision pourra être accrue de manière arbitraire. Dans ce but, les principaux ingrédients qui ont été mis en oeuvre sont (i) la méthode des éléments finis basée sur des approximations de Galerkin d'ordre arbitraire en espace et géométrie, (ii) la méthode ALE pour la prise en compte de la mobilité des domaines pour des déplacements d'ordre arbitraire, (iii) les couplages implicites et semi-implicites pour l'interaction fluide-structure. Nous proposons également une nouvelle formulation de la méthode de la frontière élargie visant à modéliser le transport de particules déformables immergées dans un fluide. Nos simulations numériques se sont appuyées sur la librairie de calcul Feel++, spécialisée dans la résolution d'EDP. Outre l'implémentation des modèles physiques, nous y avons développé diverses fonctionnalités nécessaires à la mise en oeuvre de nos méthodes : interpolation, méthode de Galerkin non standard, méthode ALE, environnement pour l'interaction fluide-structure. De plus, de par la taille des géométries et la complexité des modèles mis en jeu, le passage au calcul parallèle a été indispensable pour pouvoir réaliser nos simulations. Ainsi, nous avons décrit le développement qui a été effectué dans cette librairie pour permettre le déploiement de nos programmes sur des architectures parallèles.