Hybrid mimetic finite-difference and virtual element formulation for coupled poromechanics

We present a hybrid mimetic finite-difference and virtual element formulation for coupled single-phase poromechanics on unstructured meshes. The key advantage of the scheme is that it is convergent on complex meshes containing highly distorted cells with arbitrary shapes. We use a local pressure-jump stabilization method based on unstructured macro-elements to prevent the development of spurious pressure modes in incompressible problems approaching undrained conditions. A scalable linear solution strategy is obtained using a block-triangular preconditioner designed specifically for the saddle-point systems arising from the proposed discretization. The accuracy and efficiency of our approach are demonstrated numerically on two-dimensional benchmark problems.Comment: 25 pages, 17 figure

The Hybrid High-Order Method for Polytopal Meshes: Design, Analysis, and Applications

International audienceHybrid High-Order (HHO) methods are new generation numerical methods for models based on Partial Differential Equations with features that set them apart from traditional ones. These include: the support of polytopal meshes including non star-shaped elements and hanging nodes; the possibility to have arbitrary approximation orders in any space dimension; an enhanced compliance with the physics; a reduced computational cost thanks to compact stencil and static condensation. This monograph provides an introduction to the design and analysis of HHO methods for diffusive problems on general meshes, along with a panel of applications to advanced models in computational mechanics. The first part of the monograph lays the foundation of the method considering linear scalar second-order models, including scalar diffusion, possibly heterogeneous and anisotropic, and diffusion-advection-reaction. The second part addresses applications to more complex models from the engineering sciences: non-linear Leray-Lions problems, elasticity and incompressible fluid flows

Operator-adapted wavelets for finite-element differential forms

We introduce in this paper an operator-adapted multiresolution analysis for finite-element differential forms. From a given continuous, linear, bijective, and self-adjoint positive-definite operator L, a hierarchy of basis functions and associated wavelets for discrete differential forms is constructed in a fine-to-coarse fashion and in quasilinear time. The resulting wavelets are L-orthogonal across all scales, and can be used to derive a Galerkin discretization of the operator such that its stiffness matrix becomes block-diagonal, with uniformly well-conditioned and sparse blocks. Because our approach applies to arbitrary differential p-forms, we can derive both scalar-valued and vector-valued wavelets block-diagonalizing a prescribed operator. We also discuss the generality of the construction by pointing out that it applies to various types of computational grids, offers arbitrary smoothness orders of basis functions and wavelets, and can accommodate linear differential constraints such as divergence-freeness. Finally, we demonstrate the benefits of the corresponding operator-adapted multiresolution decomposition for coarse-graining and model reduction of linear and non-linear partial differential equations

Reactive Flows in Deformable, Complex Media

Many processes of highest actuality in the real life are described through systems of equations posed in complex domains. Of particular interest is the situation when the domain is changing in time, undergoing deformations that depend on the unknown quantities of the model. Such kind of problems are encountered as mathematical models in the subsurface, material science, or biological systems.The emerging mathematical models account for various processes at different scales, and the key issue is to integrate the domain deformation in the multi-scale context. The focus in this workshop was on novel techniques and ideas in the mathematical modelling, analysis, the numerical discretization and the upscaling of problems as described above

Error estimates for the gradient discretisation of degenerate parabolic equation of porous medium type

International audienceThe gradient discretisation method (GDM) is a generic framework for the spatial discretisation of partial differential equations. The goal of this contribution is to establish an error estimate for a class of degenerate parabolic problems, obtained under very mild regularity assumptions on the exact solution. Our study covers well-known models like the porous medium equation and the fast diffusion equations, as well as the strongly degenerate Stefan problem. Several schemes are then compared in a last section devoted to numerical results