Efficient numerical simulation of single-phase flow in three-dimensional fractured porous media

Daniel Zegarra Vasquez: Thursday, 15th May 2025 at 10:30 Abstract: Cette présentation est une répétition en vue de la soutenance de ma thèse qui aura lieu le mardi 27 mai 2025. Cette thèse porte sur les écoulements monophasiques dans des milieux poreux souterrains tridimensionnels, caractérisés par des fractures, des discontinuités étroites omniprésentes dans la matrice rocheuse. Dans ce travail, les fractures sont spécifiquement considérées comme des conduits privilégiés pour l’écoulement et sont modélisées via l’approche “Discrete Fracture Matrix” (DFM). Ce modèle conserve la tridimensionnalité de la matrice rocheuse et traite le réseau de fractures comme un objet de codimension un. L’écoulement est régi par des équations aux dérivées partielles (EDPs) de type Darcy couplées, assurant les échanges entre roche et fractures. Alors que la littérature se limite à des DFMs avec quelques milliers de fractures, cette thèse traite des modèles en comptant jusqu’à plusieurs centaines de milliers. L’objectif de ce travail est de concevoir, implémenter et analyser une méthode de simulation adaptée à ces DFMs. Le défi principal réside dans la résolution efficace des systèmes linéaires issus de la discrétisation des EDPs, pouvant atteindre plusieurs centaines de millions de degrés de liberté (DDL). Le chapitre [1] traite de l’analyse mathématique et numérique du système d’EDPs. Il fournit le cadre fonctionnel de la formulation variationnelle mixte et prouve l’existence et l’unicité de la solution. Les éléments finis mixtes sont utilisés pour poser le problème discret, pour lequel l’existence et l’unicité sont démontrées. Des estimations d’erreur a priori sont également établies. La formulation mixte-hybride équivalente (EFMH) conduit à un système linéaire creux, carré, symétrique et défini positif. La convergence numérique d’ordre 1 du schéma EFMH est validée par un test académique avec solution analytique. Le chapitre [2] propose une méthodologie pour évaluer les performances des solveurs linéaires appliqués au système issu de la…

Lire la suite

Geological storage of CO2: an example for the simulation of subsurface flow

Michel Kern: Thursday, 27th March 2025 at 10:30 Abstract: In this talk, I will present the main models used in the simulation of subsurface flow and why they are relevant for the storage of CO2 in deep geological formations. As an example, I will show results from a recent international benchmarking example. The talk is a rehearsal for my “Demi-heure de scicence” presentation on April 3rd.

Lire la suite

Hybrid high-order methods for the wave equation in first-order form

Rekha Khot: Thursday, 20th March 2025 at 10:30 Abstract: In this talk, we will discuss the approximation of the acoustic wave equation in its first-order Friedrichs formulation using hybrid high-order (HHO) methods, proposed and numerically investigated in [Burman-Duran-Ern, 2022]. We first look at energy-error estimates in the time-continuous setting and give several examples of interpolation operators: the classical one in the HHO literature based on L2 orthogonal projections and others from, or inspired from, the hybridizable discontinuous Galerkin (HDG) literature giving improved convergence rates on simplices. The time-discrete setting is based on explicit Runge-Kutta (ERK) schemes in time combined with HHO methods in space. In the fully discrete analysis, the key observation is that it becomes crucial to bound the consistency error in space by means of the stabilization seminorm only. We formulate three abstract properties (A1)-(A3) to lead the analysis. Our main result proves that, under suitable CFL conditions for second- and third-order ERK schemes, the energy error converges optimally in time and quasi-optimally in space, with optimal rates recovered on simplicial meshes. The abstract foundations of our analysis should facilitate its application to other nonconforming hybrid methods such as HDG and weak Galerkin (WG) methods.

Lire la suite

Building (yet other) bridges between polytopal methods

Simon Lemaire: Thursday, 13th March 2025 at 10:30 Abstract: Within the last 20 years or so, a myriad of novel numerical approaches, capable of accommodating general polytopal meshes, have pop up in the literature. The main purpose of this talk is to tidy up the room, and to build connections, in the context of a model variable diffusion problem, between these different approaches. Our study will focus on skeletal methods. As opposed to plain-vanilla finite volume and discontinuous Galerkin discretizations, skeletal methods essentially attach degrees of freedom to the mesh skeleton. Our study will discriminate between primal and mixed formulations of the problem at hand. Somewhat unsurprisingly, we will see that, at the end of the day, all these approaches fall within only two distinct approximation paradigms.

Lire la suite

Shape optimisation using Lipschitz functions

Philip Herbert: Thursday, 06th March 2025 at 10:30 Abstract: In this talk, we discuss a novel method in PDE constrained shape optimisation.  While it is known that many shape optimisation problems have a solution, finding the solution, or an approximation of the solution, may prove non-trivial.  A typical approach to minimisation is to use a first order method; this raises questions when handling shapes – what is a shape derivative, where does it live?  It happens to be convenient to define the derivatives as linear functionals on $W^{1,\infty}$.  We present an analysis of this in a discrete setting along with the existence of directions of steepest descent.  Several numerical experiments will be considered and extensions discussed.

Lire la suite