In the framework of research into radioactive waste disposal, it was decided to investigate the thermally induce pore pressure occurring in the Callovo-Oxfordian claystone, a possible host rock in which the ANDRA underground laboratory of Bure (East of France) has been excavated. Thermal pore pressures appear in low permeability soils and rocks because the thermal expansion coefficient of water is significantly higher than that of the solid grains (Campanella and Mitchell; 1968 [1], Ghabezloo and Sulem; 2009 [2]). This phenomenon has clearly been observed in various in-situ heating tests conducted in Opalinus claystone in the Mont-Terri Underground Research Laboratory (URL) in Switzerland (HE-D test) and in Callovo-Oxfordian (COx) claystone in the Bure URL in France (TER test, Wileveau and Su; 2007 [3]) The processes of coring, transportation, storage and specimen trimming induce some desaturation in the sample. Due to the very low permeability (10-20 m2) of the COx claystone, a long period of time is necessary to properly resaturate the sample, a mandatory condition to satisfactorily investigate thermal pressurisation. Particular emphasis was hence put on the previous saturation procedure that was carried out under in-situ effective stress condition. Thermal pressurization has been investigated by performing undrained heating tests while measuring pore pressures changes in a specially adapted thermal isotropic compression cell. Special care was devoted to calibration procedures to account for the effects of the system on the pore pressure measurements. The thermal pressurization coefficient measured appeared to change with temperature, mainly because of the changes with temperature of both the water thermal expansion coefficient of water and the drained compression coefficient of the claystone.

We address the discretization of two-phase Darcy flows in a fractured and deformable porous medium, including frictional contact between the matrix-fracture interfaces. Fractures are described as a network of planar surfaces leading to the so-called mixed- or hybrid-dimensional models. Small displacements and a linear elastic behavior are considered for the matrix. Phase pressures are supposed to be discontinuous at matrix-fracture interfaces, as they provide a better accuracy than continuous pressure models even for high fracture permeabilities. The general gradient discretization framework is employed for the numerical analysis, allowing for a generic stability analysis and including several conforming and nonconforming discretizations. We establish energy estimates for the discretization, and prove existence of a solution. To simulate the coupled model, we employ a Two-Point Flux Approximation (TPFA) finite volume scheme for the flow and second-order ($\mathbb P_2$) finite elements for the mechanical displacement coupled with face-wise constant ($\mathbb P_0$) Lagrange multipliers on fractures, representing normal and tangential stresses, to discretize the frictional contact conditions. This choice allows to circumvent possible singularities at tips, corners, and intersections between fractures, and provides a local expression of the contact conditions. We present numerical simulations of two benchmark examples and one realistic test case based on a drying model in a radioactive waste geological storage structure.

This paper is devoted to the micro-mechanical origins of the high compressibility of brittle tubular particle assemblies. The material is extremely porous due to the presence of a large hole within the tube-shaped particle. The release of the inner void, protected by a fragile shell, gives the material a very strong ability to compress. The compressive response is investigated by means of the Discrete Element Method, DEM, using crushable-elements. To address the complexity of the model, a step-by-step break-down is developed. The paper comprises the comparison of the numerical results with both results obtained by the authors and existing experiments. With the insights provided by the DEM, we have sought to better understand the phenomena that originate at the grain scale, and that govern macroscopic behaviour. Grain breakage was proven to control the compressive behaviour, and thus, the importance of internal pores dominates the inter-particle voids. Then, a novel concept of compressibility analysis has been proposed using the separation of the double porosity and the quantification of the pore collapse through primary grain breakage. Finally, a general, geometrical development of a semi-analytical model has been proposed aiming the prediction of the evolution of double porosity vs axial strain.

In porous media physics, calibrating model parameters through experiments is a challenge. This process is plagued with errors that come from modelling, measurement and computation of the macroscopic observables through random homogenization -- the forward problem -- as well as errors coming from the parameters fitting procedure -- the inverse problem. In this work, we address these issues by considering a least-square formulation to identify parameters of the microscopic model on the basis on macroscopic observables. In particular, we discuss the selection of the macroscopic observables which we need to know in order to uniquely determine these parameters. To gain a better intuition and explore the problem without a too high computational load, we mostly focus on the one-dimensional case. We show that the Newton algorithm can be efficiently used to robustly determine optimal parameters, even if some small statistical noise is present in the system.