Une connaissance plus approfondie du comportement des roches argileuses, et plus particulièrement de leur tenue, s'avère aujourd'hui indispensable dans le cadre de plusieurs thèmes de recherche d'importance croissante : comme couverture des réservoirs pétroliers pour la production d’hydrocarbures et la séquestration de CO2, ou en tant que roche hôte pour le stockage des déchets nucléaires en couches géologiques profondes. Ces roches sont particulièrement difficiles à caractériser expérimentalement du fait de leur très faible perméabilité et de leur sensibilité à l'eau. L'objectif de cette thèse est d'améliorer la compréhension des couplages hydromécaniques intervenant dans ces matériaux. Tout d’abord, une contribution à la caractérisation de l’évolution de la structure des roches argileuses soumises à un chargement mécanique est apportée. A cet effet, un essai de compression uniaxiale sur l’argilite de Tournemire associant des observations au MEB et une méthode de corrélation d’images numériques est proposé. Par la suite, des simulations numériques des essais œdométriques ont été réalisées pour appréhender l'influence des propriétés poromécaniques, en particulier de la faible perméabilité sur le comportement de l'échantillon. L’étude des équations de diffusion et les simulations numériques effectuées nous permettent de proposer une méthode de détermination de la perméabilité intrinsèque de ces matériaux. Dans la dernière partie, une étude du comportement d’une roche couverture d’un réservoir pétrolier à grande échelle est réalisée sur un cas synthétique 3D. L’objectif de ce chapitre est d’étudier les variations du champ de contrainte au sein de la couverture au cours de la production du réservoir / Knowledge of shale poromechanical behavior is essential for various environmental issues such as CO2 deep geological storage, high level radioactive waste storage, oil field abandonment… Furthermore, due to their good confinement properties, there are considered as potential host rocks for high level radioactive waste storage. However, their experimental characterization is very complex because of their very low permeability and their saturating fluid sensitivity. The main objective of this study is to improve knowledge of hydromechanical coupling of such a material. To study the structure of argillites under stress, SEM in-situ uniaxial compression test is performed on parallelepiped sample; compression orientation is perpendicular to bedding planes direction. Correlation techniques are used to estimate strain distributions. Second part underlines the key points of shale experimental characterization using the framework of Biot’s mechanics of fluid saturated porous solids. Shales are well known as a more or less transverse isotropy material. The complete methodology to conduct œdometric tests on such sensitive and weakly permeable material is described. Measurements realized on Tournemire argillite are proposed and a comparison with the poroelastic parameters of Meuse / Haute-Marne is given. Furthermore, to investigate hydromechanical coupling, a transverse isotropic poroelastic model is used to estimate the influence of anisotropy on diffusion characteristic time. The experimental data are used in numerical simulations to estimate the intrinsic permeability of the samples. In the last part, the caprock behavior of a petroleum field is studied at large scale with a 3D synthetic model (finite element analysis). The purpose of this last chapter is to study changes in stress field during the oil production

This paper presents the theoretical background of a phenomenological biphasic material approach at large strains based on the theory of porous media as well as its numerical realization within the context of an adaptive mixed finite element formulation. The study is aimed at the simulation of coupled multiphysics problems with special focus on biomechanics. As the materials of interest can be considered as a mixture of two immiscible components (solid and fluid phases), they can be modeled as saturated porous media. For the numerical treatment of according problems within a finite element approach, weak formulations of the balance equations of momentum and volume of the mixture are developed. Within this context, a generalized Lagrangean approach is preferred assuming the initial configuration of the solid phase as reference configuration of the mixture. The transient problem results in weak formulations with respect to the displacement and pore pressure fields as well as their time derivatives. Therefore special linearization techniques are applied, and after spatial discretization a global system for the incremental solution of the initial boundary value problem within the framework of a stable mixed U/p-c finite element approach is defined. The global system is solved using an iterative solver with hierarchical preconditioning. Adaptive mesh evolution is controlled by a residual a posteriori error estimator. The accuracy and the efficiency of the numerical algorithms are demonstrated on a typical example.

[pt] Este trabalho estuda um modelo acoplado de deformação e fluxo, baseado na teoria de Biot, para solos saturados e não saturados, porém levando-se em consideração a não linearidade física do solo. Na condição não saturada o fluxo do fluido pode ser considerado bifásico ou monofásico. No caso de fluxo bifásico o modelo calcula simultaneamente a pressão de água e do ar, possibilitando assim, a avaliação dos valores da sucção. O grau de saturação é então determinado a partir dos valores obtidos de sucção. Uma vez conhecido o grau de saturação é possível avaliar a permeabilidade relativa do solo para o fluxo de água e ar, tornando o problema não linear em relação à permeabilidade. No problema de fluxo monofásico a pressão de ar é mantida constante, assumindo-se o valor da pressão atmosférica. O comportamento não linear da relação tensão-deformação em solos não saturados é avaliado pelo modelo básico de Barcelona (BBM), baseado na teoria da plasticidade e conceitos de comportamento de solos no estado crítico. Este modelo pode ser compreendido como uma generalização do modelo Cam Clay freqüentemente usado para representação do comportamento mecânico de solos saturados. A influência da sucção na resistência ao cisalhamento de solos não saturados é verificada por meio de algumas aplicações numéricas. Neste trabalho, a não linearidade física de solos secos e saturados também é simulada através do modelo elastoplástico de Lade & Kim. / [en] This dissertation studies a coupled flow and deformation model for saturated and non-saturated soils, based on Biot`s theory but taking into account the nonlinear material behavior of soil. In the non-saturated condition the fluid flow can be considered either as a two-phase or an one-phase problem. In the two-phase case, the model calculates the pressure of air and water simultaneously, thus permitting the evaluation of the pore suction values and the degree of saturation. Once known, the degree of saturation is used for the determination of the relative permeability of the soil considering the flow of water and air, solving a problem that is nonlinear with respect to permeability. In the one-phase flow, the pressure of air is kept constant and equal to the atmospheric pressure. The stress-strain nonlinear behavior of non-saturated soils is modeled through the basic Barcelona model (BBM model) based on plasticity theory and concepts from the critical state behavior of soils. This model may be understood as an extension of the modified Cam Clay model frequently used for representation of the mechanical behavior of saturated soils. The suction influence on the shear resistance of non-saturated soils is examined through some numerical applications. In this work, the material nonlinearity of dry or saturated soils is also simulated through the Lade & Kim elastoplastic model.

<p>Injection to the subcutaneous (SC) tissue is one of the preferred methods for drug delivery of pharmaceuticals, from small molecules to monoclonal antibodies. Delivery to SC has become widely popular in part thanks to the low cost, ease of use, and effectiveness of drug delivery through the use of auto-injector devices. However, injection physiology, from initial plume formation to the eventual uptake of the drug in the lymphatics, is highly dependent on SC mechanics, poroelastic properties in particular. Yet, the poroelastic properties of SC have been understudied. In this thesis, I present a two-pronged approach to understanding the poroelastic properties of SC. Experimentally, mechanical and fluid transport properties of SC were measured with confined compression experiments and compared against gelatin hydrogels used as SC-phantoms. It was found that SC tissue is a highly non-linear material that has viscoelastic and porohyperelastic dissipation mechanisms. Gelatin hydrogels showed a similar, albeit more linear response, suggesting a micromechanical mechanism may underline the nonlinear behavior. The second part of the thesis focuses on the multiscale modeling of SC to gain a fundamental understanding of how geometry and material properties of the microstructure drive the macroscale response. SC is composed of adipocytes (fat cells) embedded in a collagen network. The geometry can be characterized with Voroni-like tessellations. Adipocytes are fluid-packed, highly deformable and capable of volume change through fluid transport. Collagen is highly nonlinear and nearly incompressible. Representative volume element (RVE) simulations with different Voroni tesselations shows that the different materials, coupled with the geometry of the packing, can contribute to different material response under the different kinds of loading. Further investigation of the effect of geometry showed that cell packing density nonlinearly contributes to the macroscale response. The RVE models can be homogenized to obtain macroscale models useful in large scale finite element simulations of injection physiology. Two types of homogenization were explored: fitting to analytical constitutive models, namely the Blatz-Ko material model, or use of Gaussian process surrogates, a data-driven non-parametric approach to interpolate the macroscale response.</p>

The scope of this work is the theoretical and computational modeling of the interaction between a Newtonian fluid and a cellular biological medium attached on the surface of a vessel. First and foremost, an extensive and comprehensive review is presented with regard to the available approaches for modeling momentum transfer within cellular biological media, including single-scale-single-phase approaches, Biot's poroelasticity, mixture theory, upscaling methods and multiscale computational equation free methods. Thereafter, at the cellular biological medium level, a theoretical model is developed for the description of momentum transfer within a poroelastic biomaterial, taking into account the interaction between the extracellular fluid and the solid skeleton that consists of cells and extracellular matrix (ECM). A continuum based formulation of momentum transport in a fluid-solid system at the finer spatial scale is used as starting point, and then the method of local spatial averaging with a weight function is implemented in order to establish the partial differential equations that describe the dynamics of fluid flow and matrix deformation at the coarser (macroscopic) spatial scale. In the special case of a homogeneous medium and under certain other conditions, the derived equations become similar to those which are postulated in the theory of interacting continua (mixture theory) and Biot's theory of poroelasticity. At the vessel level, the contribution of this work is twofold. First, a benchmark problem is developed for the validation of numerical methods used to solve problems that involve interactions between a fluid and a poroelastic material. Specifically, an analytical solution is developed for the problem of plane Couette-Poiseuille flow past a poroelastic layer. Second, a computational study is performed for plane Poiseuille flow past and through a semi-elliptical poroelastic biomaterial, which is attached to the surface of a straight vessel. Fluid flow in the clear fluid region is described by the Navier-Stokes equations, and momentum transfer within the biomaterial is described by the upscaled biphasic equations established in this work. The effect of the Reynolds and Darcy number that characterize the flow past and through the biomaterial, respectively, is investigated for obstacles with different configuration with respect to flow (semicircle, oblate semi-ellipse, prolate semi-ellipse). The distribution of the von Mises stress within the biomaterial is determined and, also, the drag and lift forces exerted by the fluid on the biomaterial are calculated. / Στόχος της παρούσας εργασίας είναι η θεωρητική και υπολογιστική μοντελοποίηση της αλληλεπίδρασης μεταξύ ενός Νευτώνειου ρευστού και ενός κυτταρικού βιολογικού υλικού το οποίο βρίσκεται προσκολημμένο στην επιφάνεια ενός αγγείου. Αρχικά παρουσιάζεται μια εκτεταμένη και περιεκτική ανασκόπηση των διαθέσιμων προσεγγίσεων για τη μοντελοποίηση της μεταφοράς ορμής σε κυτταρικά βιολογικά υλικά, συμπεριλαμβανομένων των προσεγγίσεων μιας κλίμακας και μιας φάσης, της θεωρίας ποροελαστικότητας του Biot, της θεωρίας αλληλεπιδρώντων συνεχών, των τεχνικών αλλαγής κλίμακας προς τα άνω, και τέλος, των υπολογιστικών τεχνικών πολλαπλών κλιμάκων χωρις τον ορισμό καταστατικών εξισώσεων. Στην συνέχεια, στο επίπεδο του κυτταρικού βιολογικού υλικού, αναπτύσεται ένα θεωρητικό μοντέλο για την περιγραφή της μεταφοράς ορμής εντός ενός ποροελαστικού υλικού, λαμβάνοντας υπόψη την αλληλεπίδραση μεταξύ του εξωκυτταρικού ρευστού και της στερεής μήτρας που αποτελείται από τα κύτταρα και το δίκτυο εξωκυτταρικών πολυμερών. Ως σημείο εκκίνησης στην μικρότερη κλίμακα παρατήρησης, χρησιμοποιείται μια περιγραφή της μεταφοράς ορμής που βασίζεται σε ένα συνεχές μοντέλο και έπειτα εφαρμόζεται η μέθοδος χωρικής στάθμισης μέσω συνάρτησης βάρους προκειμένου να εξαχθούν οι μερικές διαφορικές εξισώσεις που περιγράφουν την δυναμική της ροής του εξωκυτταρικού ρευστού και της παραμόρφωσης της στερεής μήτρας στην μακροσκοπική κλίμακα. Για την ειδική περίπτωση ενός ομογενούς μέσου και υπό την ισχύ ορισμένων πρόσθετων συνθηκών, οι εξαχθείσες εξισώσεις λαμβάνουν μορφή παρόμοια με αυτή των αντίστοιχων εξισώσεων οι οποίες ισχύουν στην θεωρία αλληλεπιδρώντων συνεχών καθώς και στην θεωρία ποροελαστικότητας του Biot. Στο επίπεδο του αγγείου, η συνεισφορά της παρούσας εργασίας λαμβάνει χώρα σε δύο άξονες. Κατά πρώτον, αναπτύσσεται ένα πρότυπο πρόβλημα το οποίο μπορεί να χρησιμεύσει για την επαλήθευση αριθμητικών μεθόδων που χρησιμοποιούνται για την επίλυση προβλημάτων στα οποία ενέχεται η αλληλεπίδραση ενός ρευστού με ένα ποροελαστικό υλικό. Συγκεκριμένα, εξάγεται μια αναλυτική λύση σε κλειστή μορφή για το πρόβλημα της επίπεδης ροής Couette-Poiseuille μέσα και γύρω από ένα ποροελαστικό στρώμα. Κατά δεύτερον, διεξάγεται μια υπολογιστική μελέτη της επίπεδης ροής Poiseuille μέσα και γύρω από ένα ημιελλειπτικό ποροελαστικό βιολογικό υλικό, το οποίο βρίσκεται προσκολημμένο στην επιφάνεια ενός ευθύγραμμου αγγείου. Στην περιοχή καθαρού ρευστού, η ροή περιγράφεται από τις εξισώσεις Navier-Stokes , ενώ η μεταφορά ορμής εντός του βιολογικού υλικού περιγράφεται με τις εξισώσεις που εξήχθησαν σε αυτή την εργασία μέσω της μεθόδου χωρικής στάθμισης. Η επίδραση των αριθμών Reynolds και Darcy, οι οποίοι χαρακτηρίζουν τη ροή γύρω και μέσα από το βιολογικό υλικό αντίστοιχα, διερευνάται για εμπόδια με διάφορες γωμετρικές διαμορφώσεις (ημικύκλιο, και ημιέλλειψη). Προσδιορίζεται η χωρική κατανομή της τάσης von Mises εντός του βιολογικού υλικού και, επιπρόσθετα, υπολογίζονται η οπισθέλκουσα και η ανυψωτική δύναμη που ασκούνται από το ρευστό στο υλικό.

Pour les roches saturées, la comparaison entre les mesures ultrasoniques (1 MHz) au laboratoire et les mesures sismiques (100 Hz) ou de diagraphie (10 kHz) sur le terrain n’est pas directe à cause de la dispersion des vitesses des ondes. Les mécanismes impliqués dans la dépendance en fréquence sont les écoulements de fluides à différentes échelles provoqués par le passage de l’onde. La dispersion et l’atténuation des modules élastiques de roches carbonatées ont été étudiées expérimentalement. Les calcaires étudiés sont : un Lavoux, un Indiana intact et craqué thermiquement, un calcaire Urgonien de Provence (Rustrel), et un coquina pré-sel du Congo. Les mesures ont été faites sur une large gamme de fréquence, en combinant les techniques d’oscillations forcées (10-3 to 102 Hz) et ultrasoniques (1 MHz) dans une presse triaxiale, pour différentes pressions effectives. Le forçage peut être hydrostatique pour mesurer un module d’incompressibilité, ou axial pour mesurer le module de Young et le coefficient de Poisson. Pour étudier l’effet de la viscosité, les mesures ont été faites en condition sèche, puis saturée en glycérine et en eau. Le drainage global et le mécanisme d’écoulement crack-pore ont été caractérisés, en termes d’amplitude de dispersion, d’atténuation viscoélastique, et de fréquence de coupure. Pour nos échantillons, la théorie de Biot-Gassmann s’est montrée valide pour les fréquences sismiques (10-100 Hz) sauf pour l’Indiana craqué thermiquement. La dispersion liée à des écoulements crackspores a été observée pour tous les échantillons sauf le Lavoux. Les fréquences de coupures de ceux-ci sont toutes dans la gamme des fréquences des diagraphies (10 kHz) pour des conditions de saturation en eau. Un modèle simple, combinant poroélasticité et milieux effectifs, a été développé pour prédire la dispersion des modules sur toute la gamme de fréquence, et s’est montré généralement en accord avec les résultats expérimentaux. / For fluid-saturated rocks, comparing ultrasonic measurements (1 MHz) in the laboratory and seismic (100 Hz) or logging (10 kHz) measurements in the field is not straightforward due to dispersion of the bodywave velocities. The frequency-dependent mechanisms involved are wave-induced fluid flows that occur at different scales. The dispersion and the attenuation of the elastic moduli of four fluid-saturated carbonate rocks have been studied experimentally. The selected limestones were a Lavoux, an intact and thermally cracked Indiana, a Urgonian limestone from Provence (Rustrel), and a presalt coquina from offshore Congo. Measurements were done over a large frequency range, by the combination of forced oscillations (10-3 to 102 Hz) and ultrasonic measurements (1 MHz) in a triaxial cell, at various effective pressures. The forced oscillations were either hydrostatic to deduce the bulk modulus, or axial to deduce Young’s modulus and Poisson’s ratio. The measurements were done in dry-, glycerinand water-saturated conditions to investigate the effect of viscosity. For all our samples, the global drainage and the squirt-flow mechanisms were characterized experimentally, in terms of amplitude of dispersion, amount of viscoelastic attenuation, and cut-off frequencies. Biot- Gassmann’s theory was found to be valid at seismic frequencies (10-100 Hz) for all the samples except the thermally cracked Indiana. Squirt-flow transitions were observed for all the samples, except the Lavoux. The cut-off frequencies were all in the range of logging frequencies (10 kHz), for watersaturated conditions. A simple model, combining poroelasticity and the noninteraction approximation effective medium, was developed to predict the dispersion of the moduli over the whole frequency range, and was generally in agreement with the experimental results.