Activity report
RNSR: 201521772E
Research center
In partnership with:
Ecole des Ponts ParisTech
Team name:
Simulation for the Environment: Reliable and Efficient Numerical Algorithms
In collaboration with:
Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique (CERMICS)
Digital Health, Biology and Earth
Earth, Environmental and Energy Sciences
Creation of the Team: 2015 June 01, updated into Project-Team: 2017 April 01


  • A2.1.3. Object-oriented programming
  • A2.1.4. Functional programming
  • A2.4.3. Proofs
  • A6.1.1. Continuous Modeling (PDE, ODE)
  • A6.1.4. Multiscale modeling
  • A6.1.5. Multiphysics modeling
  • A6.2.1. Numerical analysis of PDE and ODE
  • A6.2.5. Numerical Linear Algebra
  • A6.2.8. Computational geometry and meshes
  • A6.3.1. Inverse problems
  • A6.3.4. Model reduction
  • A6.3.5. Uncertainty Quantification
  • B3.1. Sustainable development
  • B3.3.1. Earth and subsoil
  • B3.4.2. Industrial risks and waste
  • B3.4.3. Pollution
  • B4.1. Fossile energy production (oil, gas)
  • B4.2.1. Fission
  • B5.5. Materials

1 Team members, visitors, external collaborators

Research Scientists

  • Martin Vohralík [Team leader, Inria, Senior Researcher, HDR]
  • François Clément [Inria, Researcher]
  • Zhaonan Dong [Inria, Researcher, from Oct 2020]
  • Jean-Luc Guermond [Inria, Chair, from Mar 2020 until Apr 2020, HDR]
  • Michel Kern [Inria, Researcher]
  • Géraldine Pichot [Inria, Researcher]

Faculty Member

  • Alexandre Ern [Ecole Nationale des Ponts et Chaussées, Professor, HDR]

Post-Doctoral Fellows

  • Omar Duran [Ecole Nationale des Ponts et Chaussées]
  • Kenan Kergrene-Profit [Inria, until Aug 2020]
  • Koondanibha Mitra [Inria, from Apr 2020]

PhD Students

  • Joëlle Ferzly [IFPEN]
  • Frédéric Marazzato [CEA, until May 2020]
  • Riccardo Milani [EDF, until Dec 2020]
  • Ani Miraçi [Inria, until Dec 2020]
  • Idrissa Niakh [EDF]
  • Stefano Piccardo [Ecole Nationale des Ponts et Chaussées]
  • Morgane Steins [CEA, from Oct 2020]

Technical Staff

  • Florent Hedin [Inria, Engineer, until Oct 2020]
  • Simon Legrand [Inria, Engineer, from Sep 2020]

Interns and Apprentices

  • Théo Kaprélian [Ecole Centrale de Lyon, until Feb 2020]
  • Nathanael Tepakbong Tematio [Inria, from Mar 2020 until Jul 2020]

Administrative Assistant

  • Derya Gök [Inria]

Visiting Scientists

  • Roland Becker [Inria, February 3–7, Professor, HDR]
  • Théophile Chaumont-Frelet [Inria, February 24–28, October 12–23, Researcher]

External Collaborators

  • Erik Burman [University College London, HDR]
  • Guy Chavent [Retired]
  • François Delebecque [Inria]
  • Jérôme Jaffré [Retired, HDR]
  • Caroline Japhet [Université Paris-Nord]
  • Kenan Kergrene-Profit [Université de Technologie de Compiègne, from Sep 2020]
  • Vincent Martin [Université de Technologie de Compiègne, from May 2020]
  • Jean-Elizabeth Roberts [Retired, HDR]
  • Pierre Weis [Retired]

2 Overall objectives

The project-team SERENA is concerned with numerical methods for environmental problems. The main topics are the conception and analysis of models based on partial differential equations, the study of their precise and efficient numerical approximation, and implementation issues with special concern for reliability and correctness of programs. We are in particular interested in guaranteeing the quality of the overall simulation process. SERENA has taken over the project-team POMDAPI2 which ended on May 31, 2015. It has been given an authorization to become a joint project-team between INRIA and ENPC at the Committee of Projects, September 1st, 2016, and was created as project-team on April 10, 2017.

3 Research program

3.1 Multiphysics coupling

Within our project, we start from the conception and analysis of models based on partial differential equations (PDEs). Already at the PDE level, we address the question of coupling of different models; examples are that of simultaneous fluid flow in a discrete network of two-dimensional fractures and in the surrounding three-dimensional porous medium, or that of interaction of a compressible flow with the surrounding elastic deformable structure. The key physical characteristics need to be captured, whereas existence, uniqueness, and continuous dependence on the data are minimal analytic requirements that we seek to satisfy. At the modeling stage, we also develop model-order reduction techniques, such as the use of reduced basis techniques or proper generalized decompositions, to tackle evolutive problems, in particular in the nonlinear case, and we are also interested in developing reduced-order methods for variational inequalities such as those encountered in solid mechanics with contact and possibly also friction.

3.2 Discretization by hybrid high-order and discrete element methods

We consequently design numerical methods for the devised model. Traditionally, we have worked in the context of finite element, finite volume, mixed finite element, and discontinuous Galerkin methods. Novel classes of schemes enable the use of general polygonal and polyhedral meshes with nonmatching interfaces, and we develop them in response to a high demand from our industrial partners (namely EDF, CEA, and IFP Energies Nouvelles). In the lowest-order case, our focus is to design discrete element methods for solid mechanics. The novelty is to devise these methods to treat dynamic elastoplasticity as well as quasi-static and dynamic crack propagation. We also develop structure-preserving methods for the Navier–Stokes equations, i.e., methods that mimic algebraically at the discrete level fundamental properties of the underlying PDEs, such as conservation principles and preservation of invariants. In the higher-order case, we actively contribute to the development of hybrid high-order methods. We contribute to the numerical analysis in nonlinear cases (obstacle problem, Signorini conditions), we apply these methods to challenging problems from solid mechanics involving large deformations and plasticity, and we develop a comprehensive software implementing them. We believe that these methods belong to the future generation of numerical methods for industrial simulations; as a concrete example, the implementation of these methods in an industrial software of EDF has been completed in 2019 in the framework of the PhD thesis of Nicolas Pignet.

3.3 Domain decomposition and Newton–Krylov (multigrid) solvers

We next concentrate an intensive effort on the development and analysis of efficient solvers for the systems of nonlinear algebraic equations that result from the above discretizations. We have in the past developed Newton–Krylov solvers like the adaptive inexact Newton method, and we place a particular emphasis on parallelization achieved via the domain decomposition method. Here we traditionally specialize in Robin transmission conditions, where an optimized choice of the parameter has already shown speed-ups in orders of magnitude in terms of the number of domain decomposition iterations in model cases. We concentrate in the SERENA project on adaptation of these algorithms to the above novel discretization schemes, on the optimization of the free Robin parameter for challenging situations, and also on the use of the Ventcell transmission conditions. Another feature is the use of such algorithms in time-dependent problems in space-time domain decomposition that we have recently pioneered. This allows the use of different time steps in different parts of the computational domain and turns out to be particularly useful in porous media applications, where the amount of diffusion (permeability) varies abruptly, so that the evolution speed varies significantly from one part of the computational domain to another. Our new theme here are Newton–multigrid solvers, where the geometric multigrid solver is tailored to the specific problem under consideration and to the specific numerical method, with problem- and discretization-dependent restriction, prolongation, and smoothing. Using patchwise smoothing, we have in particular recently developed a first multigrid method whose behavior is both in theory and in practice insensitive of (robust with respect to) the approximation polynomial degree. With patchwise techniques, we also achieve mass balance at each iteration step, a highly demanded feature in most of the target applications. The solver itself is then adaptively steered at each execution step by an a posteriori error estimate (adaptive stepsize, adaptive smoothing).

3.4 Reliability by a posteriori error control

The fourth part of our theoretical efforts goes towards guaranteeing the results obtained at the end of the numerical simulation. Here a key ingredient is the development of rigorous a posteriori estimates that make it possible to estimate in a fully computable way the error between the unknown exact solution and its numerical approximation. Our estimates also allow to distinguish the different components of the overall error, namely the errors coming from modeling, from the discretization scheme, from the nonlinear (Newton) solver, and from the linear algebraic (Krylov, domain decomposition, multigrid) solver. A new concept here is that of local stopping criteria, where all the error components are balanced locally within each computational mesh element. This naturally connects all parts of the numerical simulation process and gives rise to novel fully adaptive algorithms. We also theoretically address the question of convergence of the new fully adaptive algorithms. We identify theoretical conditions so that the error diminishes at each adaptive loop iteration by a contraction factor and we in particular derive a guaranteed error reduction factor in model cases. We have also proved a numerical optimality of the derived algorithms in model cases in the sense that, up to a generic constant, the smallest possible computational effort to achieve the given accuracy is needed.

3.5 Safe and correct programming

Finally, we concentrate on the issue of computer implementation of scientific computing programs. Increasing complexity of algorithms for modern scientific computing makes it a major challenge to implement them in the traditional imperative languages popular in the community. As an alternative, the computer science community provides theoretically sound tools for safe and correct programming. We explore here the use of these tools to design generic solutions for the implementation of the class of scientific computing software that we deal with. Our focus ranges from high-level programming via functional programming with OCaml through safe and easy parallelism via skeleton parallel programming with Sklml to proofs of correctness of numerical algorithms and programs via mechanical proofs with Coq.

4 Application domains

4.1 Multiphase flows and transport of contaminants in the subsurface

  • subsurface depollution after chemical leakage
  • nuclear waste disposal in deep underground repositories
  • flow in large scale discrete fracture networks
  • production of oil and gas

4.2 Industrial risks in energy production

  • Stokes and Navier–Stokes flows related to nuclear reactor operation
  • seismic wave propagation for detection and protection
  • electromagnetism for interfaces between dielectrics and negative metamaterials

4.3 Nonlinear mechanics

  • quasi-static and dynamic elastoplastic evolutions with small and large deformations
  • quasi-static and dynamic crack propagation
  • nonlinear contact and friction conditions
  • application to engineering components mainly related to nuclear reactor operation and safety analysis

4.4 Computational quantum chemistry

  • guaranteed bounds for ground-state energy (eigenvalues) and ground-state density matrix (eigenvectors) in first-principle molecular simulation
  • application to Laplace, Gross–Pitaevskii, Kohn–Sham, and Schrödinger models

5 Social and environmental responsibility

5.1 Impact of research results

Via applications with our industrial and environmental partners EDF, CEA, IFP Energies Nouvelles, ANDRA, and BRGM.

6 Highlights of the year

Many new results of the ERC GATIPOR project in the ERC GATIPOR Gallery.

7 New software and platforms

7.1 New software

7.1.1 CELIA3D

  • Name: CELIA3D
  • Keywords: Fluid mechanics, Multi-physics simulation
  • Functional Description: The CELIA3D code simulates the coupling between a compressible fluid flow and a deformable structure. The fluid is handled by a Finite Volume method on a structured Cartesian grid. The solid is handled by a Discrete Element method (Mka3d scheme). The solid overlaps the fluid grid and the coupling is carried out with immersed boundaries (cut cells) in a conservative way.
  • URL: http://cermics.enpc.fr/~monassel/CELIA3D/
  • Contacts: Laurent Monasse, Alexandre Ern
  • Partners: Ecole des Ponts ParisTech, CEA

7.1.2 DiSk++

  • Name: Discontinuous Skeletal C++ Library
  • Keywords: High order methods, Polyhedral meshes, C++
  • Scientific Description: Discontinuous Skeletal methods approximate the solution of boundary-value problems by attaching discrete unknowns to mesh faces (hence the term skeletal) while allowing these discrete unknowns to be chosen independently on each mesh face (hence the term discontinuous). Cell-based unknowns, which can be eliminated locally by a Schur complement technique (also known as static condensation), are also used in the formulation. Salient examples of high-order Discontinuous Skeletal methods are Hybridizable Discontinuous Galerkin methods and the recently-devised Hybrid High-Order methods. Some major benefits of Discontinuous Skeletal methods are that their construction is dimension-independent and that they offer the possibility to use general meshes with polytopal cells and non-matching interfaces. The mathematical flexibility of Discontinuous Skeletal methods can be efficiently replicated in a numerical software: by using generic programming, the DiSk++ library offers an environment to allow a programmer to code mathematical problems in a way completely decoupled from the mesh dimension and the cell shape.
  • Functional Description: The software provides a numerical core to discretize partial differential equations arising from the engineering sciences (mechanical, thermal, diffusion). The discretization is based on the "Hybrid high-order" or "Discontinuous Skeletal" methods, which use as principal unknowns polynomials of arbitrary degree on each face of the mesh. An important feature of these methods is that they make it possible to treat general meshes composed of polyhedral cells. The DiSk ++ library, using generic programming techniques, makes it possible to write a code for a mathematical problem independently of the mesh. When a user writes the code for his problem using the basic operations offered by DiSk ++, that code can be executed without modifications on all types of mesh already supported by the library and those that will be added in the future.
  • URL: https://github.com/wareHHOuse/diskpp
  • Publication: hal-01429292
  • Author: Matteo Cicuttin
  • Contacts: Matteo Cicuttin, Alexandre Ern, Martin Vohralík, Nicolas Pignet
  • Partner: CERMICS

7.1.3 NEF++

  • Name: NEF++
  • Keywords: Fracture network, Finite element modelling, High order methods
  • Scientific Description: NEF++ is able to solve in a very efficient way flow in large scale fractures networks using the HHO method.
  • Functional Description: The software NEF++ allows to solve flow problems in fractured rocks. The discrete fracture network (DFN) framework is used to model the fractured rocks. As the mesh steps can be different from one fracture to another, the NEF++ software is able to manage matching and non matching meshes at the intersections between fractures. NEF++ is based on a hybrid high order (HHO) method. It is based on the C++17 standard. It relies on the Eigen library, which is a C++ template library for linear algebra and on the DiSk++ library for HHO. The linear systems can be solved with direct solvers or with iterative solvers like the preconditioned conjugate gradient or multigrid solvers.
  • Release Contributions: NEF++ is now available as a Python module.
  • News of the Year: NEF++ supports non matching meshes at the intersections between fractures.
  • URL: https://gitlab.inria.fr/nef/NEFpp
  • Publications: hal-03115385, hal-02315491
  • Authors: Géraldine Pichot, Florent Hedin
  • Contacts: Florent Hedin, Géraldine Pichot
  • Participants: Florent Hedin, Géraldine Pichot, Alexandre Ern, Nicolas Pignet

7.1.4 ParaCirce

  • Name: Parallel Circulant Embedding
  • Keywords: 2D, 3D, Hydrogeology, Gaussian random fields, MPI
  • Scientific Description: ParaCirce implements the algorithm proposed by [C. R. Dietrich and G. N. Newsam. A fast and exact method for multidimensional gaussian stochastic simulations. Water Resources Research, 29(8):2861-2869, 1993].
  • Functional Description: ParaCirce implements a parallel Circulant Embedding method for the generation in parallel of 2D or 3D Gaussian Random Fields (second order stationary).
  • Release Contributions: Add domain and function abstractions to improve modularity and simplify API, Handle floating point precision with template arguments, Toggle code in header only library,
  • News of the Year: Add [Graham et al. Analysis of circulant embedding methods for sampling stationary random fields. SIAM Vol. 56, No. 3, pp. 1871-1895, 2018] criteria to validate the method,
  • URL: https://gitlab.inria.fr/slegrand/paracirce
  • Contacts: Géraldine Pichot, Simon Legrand
  • Participants: Géraldine Pichot, Simon Legrand

7.1.5 PRune

  • Name: ParserRUNnEr
  • Keywords: Test, Parameters study
  • Scientific Description: Utility to ease parameters study.
  • Functional Description: Python tool to parse single or multi configurations parameters files and to automatically run a program and store the results in a predefined tree.
  • Release Contributions: Can be used as autonomous software, as well as python module
  • News of the Year: The study is defined in a Json file.
  • URL: https://gitlab.inria.fr/slegrand/Parser_Launcher
  • Contacts: Simon Legrand, Géraldine Pichot
  • Participants: Simon Legrand, Géraldine Pichot

7.1.6 NEF-Flow-polygons

  • Name: NEF-Flow-polygons
  • Keywords: Fracture network, Large scale, High order finite elements
  • Scientific Description: The motivation of this development comes from our interest in solving flow in Discrete Fracture Networks (DFN) with the Hybrid High Order (HHO) method. One option to reduce the computational cost in solving flow in large scale DFN is to mesh the fractures independently from one fracture to another, with a mesh step chosen in a efficient way (ideally from a posteriori estimators). Therefore the meshes are possibly non matching at the intersections. As the HHO method handles naturally elements of arbitrary shape, a vertices insertion strategy in a non matching mesh made of triangles allows to build a new matching mesh made of triangles and polygons. This new mesh is then suitable to be used in flow simulations solved with the HHO method.
  • Functional Description: NEF-flow-polygons is a Matlab software to create a matching mesh made of triangles and polygons from a non matching triangular mesh of large scale discrete fracture networks (DFN). NEF-flow-polygons is interfaced with NEF++ to solve flow in these large DFN with the HHO method.
  • Release Contributions: Implementation of the algorithm for creating geometrically conforming meshes. Set up of Conitnuous Integration (CI) on gitlab-inria.fr.
  • News of the Year: NEF-flow-polygons has been successfully applied to various DFN, among with one with 508,339 fractures meshed with 7,283,887 triangles.
  • URL: https://gitlab.inria.fr/gpichot/nef-flow-polygons
  • Publication: hal-03115385
  • Contact: Géraldine Pichot
  • Participant: Géraldine Pichot

7.1.7 PyNKRT

  • Name: A Python Newton-Krylov solver for Reactive Transport
  • Keywords: Reactive transport, Nonlinear system, Python
  • Functional Description:

    Python library for the simulation of reactive transport in porous media.

    The library couples a transport module with a geochemistry module. The transport module is provided by the ComPASS code (developed by BRGM and Inria) and the chemistry module is developed as part of this project.

  • URL: https://gitlab.inria.fr/charms/pynkrt
  • Contact: Michel Kern
  • Partner: FSTG Marrakech

8 New results

8.1 Space–time domain decomposition for nonlinear diffusion problems

Participants: Elyes Ahmed, Caroline Japhet, Michel Kern.

In 11, we completed the study in  1 by giving the details of a global in time domain-decomposition (DD) method for a problem modelling two-phase flow with discontinuous capillary pressure. We proved the existence of a solution for the local problem (a nonlinea degenerate parabolic problem with nonlinear Robin boundary conditions). The proof is by convergence of a finite volume scheme, and uses new compactness results by [Andreianov, Cancès, and Moussa, J. Funct. Anal. 273 (2017)]. We also studied in detail the convergence of the Optiomized Schwarz Waveform Relaxation (OSWR) algorithm.

Figure 1: Saturation in a model with a low capillary pressure region (left), and convergence of OSWR with one or two Robin parameters (right)

8.2 hp–Version discontinuous Galerkin methods on essentially arbitrarily-shaped elements.

Participants: Andrea Cangiani, Zhaonan Dong, Emmanuil Georgoulis.

In 50, we extend the applicability of the popular interior-penalty discontinuous Galerkin (dG) method discretizing advection-diffusion-reaction problems to meshes comprising extremely general, essentially arbitrarily-shaped elemens. In particular, our analysis allows for curved element shapes, without the use of non-linear elemental maps. The feasibility of the method relies on the definition of a suitable choice of the discontinuity penalization, which turns out to be explicitly dependent on the particular element shape, but essentially independent on small shape variations. This is achieved upon proving extensions of classical trace inverse estimate to arbitrary element shapes. A further new H1L2-type inverse estimate on essentially arbitrary element shapes enables the proof of inf-sup stability of the method in a streamline diffusion-like norm. These inverse estimates may be of independent interest. A priori error bounds for the resulting method are given under very mild structural assumptions restricting the magnitude of the local curvature of element boundaries. Moreover, we have also developed a GPU library for the massive parallel implemetation of the proposed DG methods, see 57.

Figure 2: DG methods with polygonal/cuvred elements for strongly convection-dominated problems with (Pe = 105).

8.3 Residual-based a posteriori error estimates for hp-discontinuous Galerkin discretisations of the biharmonic problem

Participants: Zhaonan Dong, Lorenzo Mascotto, Oliver Sutton.

In 26, we introduce a residual-based a posteriori error estimator for a novel hp-version interior penalty discontinuous Galerkin method for the biharmonic problem in two and three dimensions. We prove that the error estimate provides an upper bound and a local lower bound on the error, and that the lower bound is robust to the local mesh size but not the local polynomial degree. The suboptimality in terms of the polynomial degree is fully explicit and grows at most algebraically. Our analysis does not require the existence of a C1 -conforming piecewise polynomial space and is instead based on an elliptic reconstruction of the discrete solution to the H2 space and a generalised Helmholtz decomposition of the error. This is a hp-version error estimator for the biharmonic problem in two and three dimensions. The practical behavior of the estimator is investigated through numerical examples in two and three dimensions.

Figure 3: Errors and estimators (left) and Adapted meshes for singular solution (right)

8.4 A-posteriori-steered p-robust multigrid with optimal step-sizes and adaptive number of smoothing steps

Participants: Ani Miraçi, Jan Papež, Martin Vohralík.

In 36, we continue our developments started in 35. We in particular develop an iterative algebraic solver of the multigrid type, whose distinctive feature is that it is steered by an a posteriori estimator of the algebraic error. We work in the context of a second-order elliptic diffusion problem discretized by conforming finite elements of arbitrary polynomial degree p1. In contrast to an established practice, our solver employs no pre-smoothing and solely one post-smoothing step, carried over via the overlapping Schwarz (block-Jacobi) method. Moreover, it features an optimal choice of the step-sizes in the smoothing correction on each level by line search. This leads to a simple Pythagorean formula of the algebraic error in the next step in terms of the current error and level-wise and patch-wise error reductions. We show the two following results and their equivalence: the solver contracts the algebraic error independently of the polynomial degree p; and the estimator represents a two-sided p-robust bound on the algebraic error. We consider quasi-uniform or graded bisection simplicial meshes and prove mild dependence on the number of mesh levels for minimal H1-regularity and complete independence for H2-regularity. We also present a simple and effective way for the solver to adaptively choose the number of post-smoothing steps necessary at each individual level, yielding a yet improved error reduction. Numerical tests confirm p-robustness and show the benefits of the adaptive number of smoothing steps. An illustration of the performance of the developed solver is illustrated in Figure 4 (red boxes), in comparision with some established algorithms.

Comparison of iteration numbers isi_{\mathrm {s}} and CPU times of the developed solver (red boxes), in comparision with some established algorithms.
Figure 4: Comparison of iteration numbers is and CPU times of the developed solver (red boxes), in comparision with some established algorithms.

9 Bilateral contracts and grants with industry

9.1 Bilateral contracts with industry

  • Two-part contract with EDF accompanying the PhD thesis of Idrissa Niakh.
  • Two-part contract with IFP Energies Nouvelles accompanying the PhD thesis of Joëlle Ferzly.
  • A new two-part contract with ANDRA accompanying the PhD thesis of Ari Rappaport (starting in January 2021).

10 Partnerships and cooperations

10.1 International initiatives

Informal international partners

University College London, Great Britain, Texas A&M University, USA, UPC (Spain), Institute of Mathematics, Czech Academy of Sciences, Charles University, Prague, Ecole Nationale d'Ingénieurs de Tunis

10.2 International research visitors

  • Title: Curved H(div), H(curl) elements, and magnetohydrodynamics & Approximation of hyperbolic systems
  • International Partner (Institution - Laboratory - Researcher):
    • Texas A&M University (United States) - Department of Mathematics - Jean-Luc Guermond
  • Duration: 2019 - 2023
  • Start year: 2019
  • See also: https://www.math.tamu.edu/~guermond/
  • The program is articulated around two themes: (1) Theoretical aspects in finite elements and applications to multi-physics magneto-hydrodynamics; (2) Finite element approximation of hyperbolic systems and applications. The results from this research will have applications in problems related to porous media flows, magnetohydrodynamics, water management, and compressible and incompressible fluid flows.

10.3 European initiatives

10.3.1 FP7 & H2020 Projects


  • Title: Extreme-scale Mathematically-based Computational Chemistry
  • Duration: 2019 - 2025
  • Coordinator: E. Cancès, L. Grigori, Y. Maday, J.-P. Piquemal
  • Partners:
    • Inria
    • Sorbonne Université
  • Inria contact: L. Grigori
  • Summary: The goal of this ERC Synergy is to develop physical and chemical models in chemistry, condensed matter physics, molecular biology, materials science, and nanosciences, altogether with mathematically-certified and numerically-efficient algorithms, and to implement them in a scalable way on various computer architectures.
  • Web page: https://erc-emc2.eu/


  • Title: Guaranteed fully adaptive algorithms with tailored inexact solvers for complex porous media flows
  • Duration: 2015 - 2021
  • Coordinator: M. Vohralík
  • Partners:
    • Inria
  • Inria contact: M. Vohralík
  • Summary: The subject of this ERC consolidator grant are new approaches to porous media multiphase flows: inexact Newton-multigrid solvers, local stopping criteria, adaptivity, and a posteriori error control. The goal is to guarantee the overall simulation error and to speed-up importantly the present-day simulations.
  • Web page: https://project.inria.fr/gatipor/

10.4 Regional initiatives

GiS: scientific collaboration network between ten public institutions from the Paris (Ile-de-France) region, focused on natural resources and environment. The project-team SERENA is a member.

11 Dissemination

11.1 Promoting scientific activities

11.1.1 Scientific events: organisation

General chair, scientific chair

11.1.2 Scientific events: selection

Member of the conference program committees

  • Géraldine Pichot is a member of the Scientific Committee of the SIAM Conference on Mathematical & Computational Issues in the Geosciences (GS21), June 21 -24, 2021, Politecnico di Milano, Italy. Géraldine Pichot was also a member of the Scientific Commitee of FRAME2020 - FRActured MEdia: Numerical methods for fluid flow and mechanics, September 21-24, 2020.
  • Alexandre Ern is a member of the Scientific Committee of the European Finite Element Fair (EFEF) and of the European Numerical Mathematics and Advanced Applications Conference (ENUMATH).
  • Michel Kern was a member of the Program Comittee of JCAD2019, Journées Calcul et Données .

11.1.3 Journal

Member of the editorial boards

Reviewer - reviewing activities

11.1.4 Invited talks

  • Michel Kern, talk at the Workshop on Recent Developments in Modelling, Analysis, and Simulation of Processes in Porous Media, University of Erlangen–Nürmberg, March 5-6 2020.
  • Martin Vohralík was an invited speaker at Finite Volumes for Complex Applications IX, Bergen, Norway, June 2020, PANM 2020, Hejnice, Czech Republic, June 2020, and InterPore 2020, online, September 2020,

11.1.5 Leadership within the scientific community

  • Alexandre Ern is the coordinator of the Master Mathématiques et applications, Ecole nationale des ponts et chaussées.
  • Michel Kern was Chair of the selection committee for the SIAM Activity Group on Geosciences Career Prize.
  • Michel Kern is a member of the Committee on Committees and Appointments of SIAM.
  • Martin Vohralík is a member of the steering committee of Summer schools CEA–EDF–INRIA.
  • Martin Vohralík is a member of the scientific board of the IFP Energies NouvellesInria joint strategic partnership laboratory.
  • Martin Vohralík is in charge of the topic “Numerical schemes, mesh generation algorithms, and error control” in the ANDRA, BRGM, CEA, EDF, IFP Energies Nouvelles, and Total working group on High-Performance Numerical Simulation in the Geosciences (identification of common challenges and collaboration opportunities).

11.1.6 Research administration

  • François Clément is member of the Comité local d'hygiène, de sécurité et des conditions de travail (CLHSCT) of the Inria Research Center of Paris.
  • Géraldine Pichot was a member of the selection committee of the Maître de Conférence (MCF 26e) position at Université de technologie de Troyes (UTT), Mai 2020.
  • Michel Kern is in charge of the Comission de Développement Technologique of the Inria Research Center of Paris.

11.2 Teaching - Supervision - Juries

11.2.1 Teaching

  • Licence : Alexandre Ern, Numerical analysis, 77h, L3/M1, Ecole Polytechnique, France.
  • Master : Alexandre Ern, Discontinuous Galerkin methods, 20h, M2, Sorbonne University, France.
  • Master: Michel Kern, Models and numerical methods for subsurface flow, 30h, M2, Université Paris Saclay, France.
  • Master: Michel Kern, Inverse Problems, 24h, M1, Mines-ParisTech, France.
  • Master: Michel Kern, Advanced numerical analysis, 30h, M1, Institut Galilée, Université Paris-Nord, France.
  • Master: Martin Vohralík, A posteriori error estimates for efficiency and error control in numerical simulations, 36h, M2, Charles University, Prague, Czech Republic.
  • Master: Martin Vohralík, Advanced finite elements, 21h, M1, Ecole nationale supérieure de techniques avancées, Paris, France.

11.2.2 Supervision

  • PhD defended on May 29, 2020: Frédéric Marazzato, Discrete element methods for fracture and fragmentation, Alexandre Ern.
  • PhD defended on December 16, 2020: Riccardo Milani, Compatible Discrete Operator schemes for Navier–Stokes equations, Alexandre Ern.
  • PhD defended on December 14, 2020: Ani Miraçi, A-posteriori-steered and adaptive p-robust multigrid solvers, Martin Vohralík and Alexandre Ern.
  • PhD in progress: Joëlle Ferzly, A posteriori error estimates and adaptivity for complementarity problems, started 15 October 2019, Martin Vohralík.
  • PhD in progress: Idrissa Niakh, Reduced models for variational inequalities in computational mechanics, started 01 November 2019, Alexandre Ern and Virginie Ehrlacher (ENPC).
  • PhD in progress: Stefano Piccardo, High-fidelity simulation of droplets in complex shear flows, started 01 October 2019, Alexandre Ern and Antonio Huerta (UPC).
  • PhD in progress: Morgane Steins, Hybrid High Order methods for explicit dynamics of structures with adaptive mesh refinement, started 01 October 2020, Alexandre Ern and Martin Vohralík.
  • Internship: Nathanael Tepakbong, from March 2020 until July 2020, ParaCirce: logiciel de génération de champs aléatoires Gaussiens, Géraldine Pichot.
  • Internship: Théo Kaprélian, from September 2019 until February 2020, An hp-adaptive strategy in dimension one taking into account the cost of the resolution, Martin Vohralík.

11.2.3 Juries

  • Géraldine Pichot, Examiner, PhD Quentin Van Aken, Université de technologie de Troyes (UTT), September 2020.
  • Michel Kern, Examiner, PhD Bastien Hamlat, September 2020

11.3 Popularization

11.3.1 Internal or external Inria responsibilities

  • Géraldine Pichot is a member of the editorial committee of Interstices.

11.3.2 Interventions

12 Scientific production

12.1 Major publications

  • 1 articleE. Ahmed, S. Ali Hassan, C. Japhet, M. Kern and M. Vohralík. 'A posteriori error estimates and stopping criteria for space-time domain decomposition for two-phase flow between different rock types'.SMAI Journal of Computational Mathematics5December 2019, 195-227
  • 2 articleS. Boldo, F. Clément, J.-C. Filliâtre, M. Mayero, G. Melquiond and P. Weis. 'Trusting computations: A mechanized proof from partial differential equations to actual program'.Computers and Mathematics with Applications683August 2014, 325--352
  • 3 articleS. Boldo, F. Clément, J.-C. Filliâtre, M. Mayero, G. Melquiond and P. Weis. 'Wave equation numerical resolution: a comprehensive mechanized proof of a C program'.Journal of Automated Reasoning504April 2013, 423--456
  • 4 articleE. Cancès, G. Dusson, Y. Maday, B. Stamm and M. Vohralík. 'Guaranteed and robust a posteriori bounds for Laplace eigenvalues and eigenvectors: conforming approximations'.SIAM Journal on Numerical Analysis555September 2017, 2228-2254
  • 5 unpublishedA. Cangiani, Z. Dong and E. Georgoulis. 'hp -Version discontinuous Galerkin methods on essentially arbitrarily-shaped elements'.January 2021, working paper or preprint
  • 6 articleD. Di Pietro and A. Ern. 'A hybrid high-order locking-free method for linear elasticity on general meshes'.Comput. Methods Appl. Mech. Engrg.2832015, 1--21URL: http://dx.doi.org/10.1016/j.cma.2014.09.009
  • 7 articleA. Ern and J.-L. Guermond. 'Finite element quasi-interpolation and best approximation'.ESAIM Math. Model. Numer. Anal.5142017, 1367--1385URL: https://doi.org/10.1051/m2an/2016066
  • 8 articleA. Ern and M. Vohralík. 'Polynomial-degree-robust a posteriori estimates in a unified setting for conforming, nonconforming, discontinuous Galerkin, and mixed discretizations'.SIAM Journal on Numerical Analysis532April 2015, 1058-1081
  • 9 articleT.-T. Hoang, J. Jaffré, C. Japhet, M. Kern and J. Roberts. 'Space-time domain decomposition methods for diffusion problems in mixed formulations'.SIAM J. Numer. Anal.5162013, 3532--3559URL: http://dx.doi.org/10.1137/130914401
  • 10 articleG. Pichot, J. Erhel and J.-R. De Dreuzy. 'A generalized mixed hybrid mortar method for solving flow in stochastic discrete fracture networks'.SIAM J. Sci. Comput.3412012, B86--B105URL: http://dx.doi.org/10.1137/100804383

12.2 Publications of the year

International journals

  • 11 articleE. Ahmed, C. Japhet and M. Kern. 'Space-time domain decomposition for two-phase flow between different rock types'.Computer Methods in Applied Mechanics and Engineering371November 2020, 113294
  • 12 articleL. Amir and M. Kern. 'Jacobian Free Methods for Coupling Transport with Chemistry in Heterogenous Porous Media'.Water133February 2021, 370
  • 13 articleI. Ben Gharbia, J. Dabaghi, V. Martin and M. Vohralík. 'A posteriori error estimates for a compositional two-phase flow with nonlinear complementarity constraints'.Computational Geosciences24January 2020, 1031-1055
  • 14 articleJ. Blechta, J. Málek and M. Vohralík. 'Localization of the W -1,q norm for local a posteriori efficiency'.IMA Journal of Numerical Analysis402April 2020, 914-950
  • 15 article E. Burman, M. Cicuttin, G. Delay and A. Ern. 'An unfitted Hybrid High-Order method with cell agglomeration for elliptic interface problems'. SIAM Journal on Scientific Computing 2021
  • 16 article E. Burman, G. Delay and A. Ern. 'An unfitted Hybrid High-Order method for the Stokes interface problems'. IMA Journal of Numerical Analysis 2021
  • 17 articleE. Cancès, G. Dusson, Y. Maday, B. Stamm and M. Vohralík. 'Guaranteed a posteriori bounds for eigenvalues and eigenvectors: multiplicities and clusters'.Mathematics of Computation89326July 2020, 2563-2611
  • 18 article E. Cancès, G. Dusson, Y. Maday, B. Stamm and M. Vohralík. 'Post-processing of the planewave approximation of Schrödinger equations. Part I: linear operators'. IMA Journal of Numerical Analysis September 2020
  • 19 articleC. Cancès, F. Nabet and M. Vohralík. 'Convergence and a posteriori error analysis for energy-stable finite element approximations of degenerate parabolic equations'.Mathematics of Computation902021, 517-563
  • 20 article T. Chaumont-Frelet, A. Ern and M. Vohralík. 'On the derivation of guaranteed and p-robust a posteriori error estimates for the Helmholtz equation'. Numerische Mathematik 2021
  • 21 articleT. Chaumont-Frelet, A. Ern and M. Vohralík. 'Polynomial-degree-robust H(curl)-stability of discrete minimization in a tetrahedron'.Comptes Rendus Mathématique3589-10December 2020, 1101-1110
  • 22 articleF. Chouly, A. Ern and N. Pignet. 'A Hybrid High-Order discretization combined with Nitsche's method for contact and Tresca friction in small strain elasticity'.SIAM Journal on Scientific Computing422020, 2300-2324
  • 23 articleJ. Dabaghi, V. Martin and M. Vohralík. 'A posteriori estimates distinguishing the error components and adaptive stopping criteria for numerical approximations of parabolic variational inequalities'.Computer Methods in Applied Mechanics and Engineering367August 2020, 113105
  • 24 articleJ. Dabaghi, V. Martin and M. Vohralík. 'Adaptive inexact semismooth Newton methods for the contact problem between two membranes'.Journal of Scientific Computing84July 2020, 28
  • 25 articleP. Daniel, A. Ern and M. Vohralík. 'An adaptive hp-refinement strategy with inexact solvers and computable guaranteed bound on the error reduction factor'.Computer Methods in Applied Mechanics and Engineering359February 2020, 112607
  • 26 article Z. Dong, L. Mascotto and O. Sutton. 'Residual-based a posteriori error estimates for hp-discontinuous Galerkin discretisations of the biharmonic problem'. SIAM Journal on Numerical Analysis 2021
  • 27 article A. Ern, T. Gudi, I. Smears and M. Vohralík. 'Equivalence of local-and global-best approximations, a simple stable local commuting projector, and optimal hp approximation estimates in H(div)'. IMA Journal of Numerical Analysis 2021
  • 28 articleA. Ern and M. Vohralík. 'Stable broken H1 and H(div) polynomial extensions for polynomial-degree-robust potential and flux reconstruction in three space dimensions'.Mathematics of Computation89322March 2020, 551-594
  • 29 article A. Ern, M. Vohralík and M. Zakerzadeh. 'Guaranteed and robust L 2 -norm a posteriori error estimates for 1D linear advection problems'. ESAIM: Mathematical Modelling and Numerical Analysis 2021
  • 30 article A. Haberl, D. Praetorius, S. Schimanko and M. Vohralík. 'Convergence and quasi-optimal cost of adaptive algorithms for nonlinear operators including iterative linearization and algebraic solver'. Numerische Mathematik 2021
  • 31 articleM. Kern, A. Taakili and M. Zarrouk. 'Preconditioned iterative method for reactive transport with sorption in porous media'.Mathematical Modelling and Analysis2542020, 546 - 568
  • 32 article C. Lusso, F. Bouchut, A. Ern and A. Mangeney. 'Explicit solutions to a free interface model for the static/flowing transition in thin granular flows'. ESAIM: Mathematical Modelling and Numerical Analysis 2020
  • 33 articleG. Mallik, M. Vohralík and S. Yousef. 'Goal-oriented a posteriori error estimation for conforming and nonconforming approximations with inexact solvers'.Journal of Computational and Applied Mathematics366March 2020, 112367
  • 34 article F. Marazzato, A. Ern and L. Monasse. 'A variational discrete element method for quasi-static and dynamic elasto-plasticity'. International Journal for Numerical Methods in Engineering 2020
  • 35 articleA. Miraçi, J. Papež and M. Vohralík. 'A multilevel algebraic error estimator and the corresponding iterative solver with p-robust behavior'.SIAM Journal on Numerical Analysis585October 2020, 2856-2884
  • 36 article A. Miraçi, J. Papež and M. Vohralík. 'A-posteriori-steered p-robust multigrid with optimal step-sizes and adaptive number of smoothing steps'. SIAM Journal on Scientific Computing 2021
  • 37 article A. Miraçi, J. Papež and M. Vohralík. 'Contractive local adaptive smoothing based on Dörfler's marking in a-posteriori-steered p-robust multigrid solvers'. Computational Methods in Applied Mathematics 2021
  • 38 articleJ. Papež, U. Rüde, M. Vohralík and B. Wohlmuth. 'Sharp algebraic and total a posteriori error bounds for h and p finite elements via a multilevel approach. Recovering mass balance in any situation'.Computer Methods in Applied Mechanics and Engineering371August 2020, 113243
  • 39 articleI. Smears and M. Vohralík. 'Simple and robust equilibrated flux a posteriori estimates for singularly perturbed reaction-diffusion problems'.ESAIM: Mathematical Modelling and Numerical Analysis546October 2020, 1951-1973

Conferences without proceedings

  • 40 inproceedings F. Hédin and G. Pichot. 'A hybrid high order (HHO) method with non matching meshes in discrete fracture networks'. ALGORITMY 2020 - Conference on Scientific Computing Podbanske, Slovakia September 2020

Doctoral dissertations and habilitation theses

  • 41 thesis F. Marazzato. 'Discrete element and time-integration methods for elasto-plasticity and dynamic cracking'. Université Paris-Est Marne-la-Vallée (UPEM) May 2020
  • 42 thesis R. Milani. 'Compatible Discrete Operator schemes for the unsteady incompressible Navier–Stokes equations'. Université Paris-Est December 2020
  • 43 thesis A. Miraçi. 'A-posteriori-steered and adaptive p-robust multigrid solvers'. INRIA Paris; Sorbonne Universite December 2020

Reports & preprints

  • 44 misc L. Amir, M. Kern, Z. Mghazli and J. Roberts. 'Intersecting fractures in porous media: theoretical and numerical analysis'. February 2021
  • 45 misc E. Burman, G. Delay and A. Ern. 'A hybridized high-order method for unique continuation subject to the Helmholtz equation'. October 2020
  • 46 misc E. Burman, G. Delay and A. Ern. 'The unfitted HHO method for the Stokes problem on curved domains'. January 2020
  • 47 misc E. Burman, O. Duran and A. Ern. 'Hybrid high-order methods for the acoustic wave equation in the time domain'. August 2020
  • 48 misc E. Burman, O. Duran, A. Ern and M. Steins. 'Convergence analysis of hybrid high-order methods for the wave equation'. December 2020
  • 49 misc E. Burman, O. Duran and A. Ern. 'Unfitted hybrid high-order methods for the wave equation'. December 2020
  • 50 misc A. Cangiani, Z. Dong and E. Georgoulis. 'hp -Version discontinuous Galerkin methods on essentially arbitrarily-shaped elements'. January 2021
  • 51 misc C. Carstensen, A. Ern and S. Puttkammer. 'Guaranteed lower bounds on eigenvalues of elliptic operators with a hybrid high-order method'. October 2020
  • 52 misc T. Chaumont-Frelet, A. Ern and M. Vohralík. 'Stable broken H(curl) polynomial extensions and p-robust quasi-equilibrated a posteriori estimators for Maxwell's equations'. May 2020
  • 53 misc T. Chaumont-Frelet and M. Vohralík. 'Equivalence of local-best and global-best approximations in H(curl)'. June 2020
  • 54 reportF. Clément and V. Martin. 'Lebesgue integration. Detailed proofs to be formalized in Coq'.Inria ParisJanuary 2021, 284
  • 55 misc P. Daniel and M. Vohralík. 'Guaranteed contraction of adaptive inexact hp-refinement strategies with realistic stopping criteria'. February 2020
  • 56 misc Q. Deng and A. Ern. 'SoftFEM: revisiting the spectral finite element approximation of elliptic operators'. November 2020
  • 57 misc Z. Dong, E. Georgoulis and K. Thomas. 'GPU-accelerated discontinuous Galerkin methods on polytopic meshes'. January 2021
  • 58 misc Z. Dong and L. Mascotto. 'On the suboptimality of the p-version discontinuous Galerkin methods for first order hyperbolic problems'. January 2021
  • 59 misc F. Marazzato, A. Ern and L. Monasse. 'Quasi-static crack propagation with a Griffith criterion using a discrete element method'. January 2021
  • 60 misc R. Milani, J. Bonelle and A. Ern. 'Compatible Discrete Operator schemes for the steady incompressible Stokes and Navier-Stokes equations'. February 2020
  • 61 misc J. Papež and M. Vohralík. 'Inexpensive guaranteed and efficient upper bounds on the algebraic error in finite element discretizations'. February 2021

12.3 Cited publications