BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jerome Droniou

Jerome Droniou appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysismath.AP

Catalog footprint

What is connected

17works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

A cost-effective nonlinear extremum-preserving finite volume scheme for highly anisotropic diffusion on Cartesian grids, with application to radiation belt dynamics

We construct a new nonlinear finite volume (FV) scheme for highly anisotropic diffusion equations, that satisfies the discrete minimum-maximum principle. The construction relies on the linearized scheme satisfying less restrictive monotonicity conditions than those of an M-matrix, based on a weakly regular matrix splitting and using the Cartesian structure of the mesh (extension to quadrilateral meshes is also possible). The resulting scheme, obtained by expressing fluxes as nonlinear combinations of linear fluxes, has a larger stencil than other nonlinear positivity preserving or minimum-maximum principle preserving schemes. Its larger "linearized" stencil, closer to the actual complete stencil (that includes unknowns appearing in the convex combination coefficients), enables a faster convergence of the Picard iterations used to compute the solution of the scheme. Steady state dimensionless numerical tests as well as simulations of the highly anisotropic diffusion in electron radiation belts show a second order of convergence of the new scheme and confirm its computational efficiency compared to usual nonlinear FV schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A discrete de Rham method for the Reissner-Mindlin plate bending problem on polygonal meshes

In this work we propose a discretisation method for the Reissner--Mindlin plate bending problem in primitive variables that supports general polygonal meshes and arbitrary order. The method is inspired by a two-dimensional discrete de Rham complex for which key commutation properties hold that enable the cancellation of the contribution to the error linked to the enforcement of the Kirchhoff constraint. Denoting by $k\ge 0$ the polynomial degree for the discrete spaces and by $h$ the meshsize, we derive for the proposed method an error estimate in $h^{k+1}$ for general $k$, as well as a locking-free error estimate for the lowest-order case $k=0$. The theoretical results are validated on a complete panel of numerical tests.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Conditioning of a Hybrid High-Order scheme on meshes with small faces

We conduct a condition number analysis of a Hybrid High-Order (HHO) scheme for the Poisson problem. We find the condition number of the statically condensed system to be independent of the number of faces in each element, or the relative size between an element and its faces. The dependence of the condition number on the polynomial degree is tracked. Next, we consider HHO schemes on cut background meshes, which are commonly used in unfitted discretisations. It is well known that the linear systems obtained on these meshes can be arbitrarily ill-conditioned due to the presence of sliver-cut and small-cut elements. We show that the condition number arising from HHO schemes on such meshes is not as negatively effected as those arising from conforming methods. We describe how the condition number can be improved by aggregating ill-conditioned elements with their neighbours.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Robust Hybrid High-Order method on polytopal meshes with small faces

We design a Hybrid High-Order (HHO) scheme for the Poisson problem that is fully robust on polytopal meshes in the presence of small edges/faces. We state general assumptions on the stabilisation terms involved in the scheme, under which optimal error estimates (in discrete and continuous energy norms, as well as $L^2$-norm) are established with multiplicative constants that do not depend on the maximum number of faces in each element, or the relative size between an element and its faces. We illustrate the error estimates through numerical simulations in 2D and 3D on meshes designed by agglomeration techniques (such meshes naturally have elements with a very large numbers of faces, and very small faces).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An efficient implementation of mass conserving characteristic-based schemes in 2D and 3D

In this paper, we develop the ball-approximated characteristics (B-char) method, which is an algorithm for efficiently implementing characteristic-based schemes in 2D and 3D. Core to the implementation of numerical schemes is the evaluation of integrals, which in the context of characteristic-based schemes with piecewise constant approximations boils down to computing the intersections between two regions. In the literature, these regions are approximated by polytopes (polygons in 2D and polyhedra in 3D) and, due to this, the implementation in 3D is nontrivial. The main novelty in this paper is the approximation of the regions by balls, whose intersections are much cheaper to compute than those of polytopes. Of course, balls cannot fully tessellate a region, and hence some mass may be lost. We perform some adjustments, and also solve an optimisation problem, in order to yield a scheme that is both locally and globally mass conserving. This algorithm can achieve results that are similar to those obtained from an implementation which uses polytopal intersections, with a much cheaper computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convergence analysis of a numerical scheme for a tumour growth model

We consider a one--spatial dimensional tumour growth model [2, 3, 4] that consists of three dependent variables of space and time: volume fraction of tumour cells, velocity of tumour cells, and nutrient concentration. The model variables satisfy a coupled system of semilinear advection equation (hyperbolic), simplified linear Stokes equation (elliptic), and semilinear diffusion equation (parabolic) with appropriate conditions on the time-dependent boundary, which is governed by an ordinary differential equation. We employ a reformulation of the model defined in a larger, fixed time-space domain to overcome some theoretical difficulties related to the time-dependent boundary. This reformulation reduces the complexity of the model by removing the need to explicitly track the time-dependent boundary, but nonlinearities in the equations, noncoercive operators in the simplified Stokes equation, and interdependence between the unknown variables still challenge the proof of suitable a priori estimates. A numerical scheme that employs a finite volume method for the hyperbolic equation, a finite element method for the elliptic equation, and a backward Euler in time--mass lumped finite element in space method for the parabolic equation is developed. We establish the existence of a time interval $(0,T_{\ast})$ over which, using compactness techniques, we can extract a convergent subsequence of the numerical approximations. The limit of any such convergent subsequence is proved to be a weak solution of the continuous model in an appropriate sense, which we call a threshold solution. Numerical tests and justifications that confirm the theoretical findings conclude the paper.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Design and convergence analysis of numerical methods for stochastic evolution equations with Leray-Lions operator

*The gradient discretisation method (GDM) is a generic framework, covering many classical methods (Finite Elements, Finite Volumes, Discontinuous Galerkin, etc.), for designing and analysing numerical schemes for diffusion models. In this paper, we study the GDM for a general stochastic evolution problem based on a Leray--Lions type operator. The problem contains the stochastic $p$-Laplace equation as a particular case. The convergence of the Gradient Scheme (GS) solutions is proved by using Discrete Functional Analysis techniques, Skorohod theorem and the Kolmogorov test. In particular, we provide an independent proof of the existence of weak martingale solutions for the problem. In this way, we lay foundations and provide techniques for proving convergence of the GS approximating stochastic partial differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-conforming finite elements on polytopal meshes

In this work we present a generic framework for non-conforming finite elements on polytopal meshes, characterised by elements that can be generic polygons/polyhedra. We first present the functional framework on the example of a linear elliptic problem representing a single-phase flow in porous medium. This framework gathers a wide variety of possible non-conforming methods, and an error estimate is provided for this simple model. We then turn to the application of the functional framework to the case of a steady degenerate elliptic equation, for which a mass-lumping technique is required; here, this technique simply consists in using a different --piecewise constant-- function reconstruction from the chosen degrees of freedom. A convergence result is stated for this degenerate model. Then, we introduce a novel specific non-conforming method, dubbed Locally Enriched Polytopal Non-Conforming (LEPNC). These basis functions comprise functions dedicated to each face of the mesh (and associated with average values on these faces), together with functions spanning the local $\mathbb{P}^1$ space in each polytopal element. The analysis of the interpolation properties of these basis functions is provided, and mass-lumping techniques are presented. Numerical tests are presented to assess the efficiency and the accuracy of this method on various examples. Finally, we show that generic polytopal non-conforming methods, including the LEPNC, can be plugged into the gradient discretization method framework, which makes them amenable to all the error estimates and convergence results that were established in this framework for a variety of models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The gradient discretisation method for slow and fast diffusion porous media equations

The gradient discretisation method (GDM) is a generic framework for designing and analysing numerical schemes for diffusion models. In this paper, we study the GDM for the porous medium equation, including fast diffusion and slow diffusion models, and a concentration-dependent diffusion tensor. Using discrete functional analysis techniques, we establish a strong $L^2$-convergence of the approximate gradients and a uniform-in-time convergence for the approximate solution, without assuming non-physical regularity assumptions on the data or continuous solution. Being established in the generic GDM framework, these results apply to a variety of numerical methods, such as finite volume, (mass-lumped) finite elements, etc. The theoretical results are illustrated, in both fast and slow diffusion regimes, by numerical tests based on two methods that fit the GDM framework: mass-lumped conforming $\mathbb{P}_1$ finite elements and the Hybrid Mimetic Mixed method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Uniform-in-time convergence of numerical methods for non-linear degenerate parabolic equations

Gradient schemes is a framework that enables the unified convergence analysis of many numerical methods for elliptic and parabolic partial differential equations: conforming and non-conforming Finite Element, Mixed Finite Element and Finite Volume methods. We show here that this framework can be applied to a family of degenerate non-linear parabolic equations (which contain in particular the Richards', Stefan's and Leray--Lions' models), and we prove a uniform-in-time strong-in-space convergence result for the gradient scheme approximations of these equations. In order to establish this convergence, we develop several discrete compactness tools for numerical approximations of parabolic models, including a discontinuous Ascoli-Arzelà theorem and a uniform-in-time weak-in-space discrete Aubin-Simon theorem. The model's degeneracies, which occur both in the time and space derivatives, also requires us to develop a discrete compensated compactness result.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Introduction to discrete functional analysis techniques for the numerical study of diffusion equations with irregular data

We give an introduction to discrete functional analysis techniques for stationary and transient diffusion equations. We show how these techniques are used to establish the convergence of various numerical schemes without assuming non-physical regularity on the data. For simplicity of exposure, we mostly consider linear elliptic equations, and we briefly explain how these techniques can be adapted and extended to non-linear time-dependent meaningful models (Navier--Stokes equations, flows in porous media, etc.). These convergence techniques rely on discrete Sobolev norms and the translation to the discrete setting of functional analysis results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical analysis of a two-phase flow discrete fracture model

We present a new model for two phase Darcy flows in fractured media, in which fractures are modelled as submanifolds of codimension one with respect to the surrounding domain (matrix). Fractures can act as drains or as barriers, since pressure discontinuities at the matrix-fracture interfaces are permitted. Additionally, a layer of damaged rock at the matrix-fracture interfaces is accounted for. The numerical analysis is carried out in the general framework of the Gradient Discretisation Method. Compactness techniques are used to establish convergence results for a wide range of possible numerical schemes; the existence of a solution for the two phase flow model is obtained as a byproduct of the convergence analysis. A series of numerical experiments conclude the paper, with a study of the influence of the damaged layer on the numerical solution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gradient schemes: generic tools for the numerical analysis of diffusion equations

The gradient scheme framework is based on a small number of properties and encompasses a large number of numerical methods for diffusion models. We recall these properties and develop some new generic tools associated with the gradient scheme framework. These tools enable us to prove that classical schemes are indeed gradient schemes, and allow us to perform a complete and generic study of the well-known (but rarely well-studied) mass lumping process. They also allow an easy check of the mathematical properties of new schemes, by developing a generic process for eliminating unknowns via barycentric condensation, and by designing a concept of discrete functional analysis toolbox for schemes based on polytopal meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Finite volume schemes for diffusion equations: introduction to and review of modern methods

We present Finite Volume methods for diffusion equations on generic meshes, that received important coverage in the last decade or so. After introducing the main ideas and construction principles of the methods, we review some literature results, focusing on two important properties of schemes (discrete versions of well-known properties of the continuous equation): coercivity and minimum-maximum principles. Coercivity ensures the stability of the method as well as its convergence under assumptions compatible with real-world applications, whereas minimum-maximum principles are crucial in case of strong anisotropy to obtain physically meaningful approximate solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gradient Schemes for Linear and Non-linear Elasticity Equations

The Gradient Scheme framework provides a unified analysis setting for many different families of numerical methods for diffusion equations. We show in this paper that the Gradient Scheme framework can be adapted to elasticity equations, and provides error estimates for linear elasticity and convergence results for non-linear elasticity. We also establish that several classical and modern numerical methods for elasticity are embedded in the Gradient Scheme framework, which allows us to obtain convergence results for these methods in cases where the solution does not satisfy the full $H^2$-regularity or for non-linear models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On a miscible displacement model in porous media flow with measure data

We establish the existence of a solution to a non-linearly coupled elliptic-parabolic system of PDEs describing the single-phase, miscible displacement of one incompressible fluid by another in a porous medium. We consider a velocity-dependent diffusion-dispersion tensor and model the action of injection and production wells by measures on the domain. We obtain the solution by passing to the limit on problems with regularised well data.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

A unified approach to Mimetic Finite Difference, Hybrid Finite Volume and Mixed Finite Volume methods

We investigate the connections between several recent methods for the discretization of anisotropic heterogeneous diffusion operators on general grids. We prove that the Mimetic Finite Difference scheme, the Hybrid Finite Volume scheme and the Mixed Finite Volume scheme are in fact identical up to some slight generalizations. As a consequence, some of the mathematical results obtained for each of the method (such as convergence properties or error estimates) may be extended to the unified common framework. We then focus on the relationships between this unified method and nonconforming Finite Element schemes or Mixed Finite Element schemes, obtaining as a by-product an explicit lifting operator close to the ones used in some theoretical studies of the Mimetic Finite Difference scheme. We also show that for isotropic operators, on particular meshes such as triangular meshes with acute angles, the unified method boils down to the well-known efficient two-point flux Finite Volume scheme.

0 citations0 reviewsNo code linkedOpen access