BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jerome Droniou

Jerome Droniou contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.AP

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 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
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