BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Robert Eymard appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We consider the Stefan problem, firstly with regular data and secondly with irregular data. In both cases is given a proof for the convergence of an approximation obtained by regularising the problem. These proofs are based on weak formulations and on compactness results in some Sobolev spaces with negative exponents.
Given a densely defined skew-symmetric operators A 0 on a real or complex Hilbert space V , we parametrize all m-dissipative extensions in terms of contractions $Φ$ : H-$\rightarrow$ H + , where Hand H + are Hilbert spaces associated with a boundary quadruple. Such an extension generates a unitary C 0-group if and only if $Φ$ is a unitary operator. As corollary we obtain the parametrization of all selfadjoint extensions of a symmetric operator by unitary operators from Hto H +. Our results extend the theory of boundary triples initiated by von Neumann and developed by V. I. and M. L. Gorbachuk, J. Behrndt and M. Langer, S. A. Wegner and many others, in the sense that a boundary quadruple always exists (even if the defect indices are different in the symmetric case).
The Representation Theorem of Lions (RTL) is a version of the Lax--Milgram Theorem where completeness of one of the spaces is not complete. In this paper, RTL is deduced from an operator-theoretical version on normed space. The main point of the paper is a theory of derivations, based on RTL, for which well-posedness is proved. One application concerns non-autonomous evolution equations with a new initial-value and a periodic boundary condition for the time variable.
We approximate the solution to some linear and degenerate quasi-linear problem involving a linear elliptic operator (like the semi-discrete in time implicit Euler approximation of Richards and Stefan equations) with measure right-hand side and heterogeneous anisotropic diffusion matrix. This approximation is obtained through the addition of a weighted p--Laplace term. A well chosen diffeomorphism between R and (--1, 1) is used for the estimates of the approximated solution, and is involved in the above weight. We show that this approximation converges to a weak sense of the problem for general right-hand-side, and to the entropy solution in the case where the right-hand-side is in L 1 .
We show that, using the Crouzeix-Raviart scheme, a cheap algebraic transformation, applied to the coupled velocity-pressure linear systems issued from the transient or steady Stokes or Navier-Stokes problems, leads to a linear system only involving as many auxiliary variables as the velocity components. This linear system, which is symmetric positive definite in the case of the transient Stokes problem and symmetric invertible in the case of the steady Stokes problem, with the same stencil as that of the velocity matrix, provides the exact solution of the initial coupled linear system. Numerical results show the increase of performance when applying direct or iterative solvers to the resolution of these linear systems.
In this paper we study the conforming Galerkin approximation of the problem: find u $\in$ U such that a(u, v) = <L, v> for all v $\in$ V, where U and V are Hilbert or Banach spaces, a is a continuous bilinear or sesquilinear form and L $\in$ V' a given data. The approximate solution is sought in a finite dimensional subspace of U, and test functions are taken in a finite dimensional subspace of V. We provide a necessary and sufficient condition on the form a for convergence of the Galerkin approximation, which is also equivalent to convergence of the Galerkin approximation for the adjoint problem. We also characterize the fact that U has a finite dimensional Schauder decomposition in terms of properties related to the Galerkin approximation. In the case of Hilbert spaces, we prove that the only bilinear or sesquilinear forms for which any Galerkin approximation converges (this property is called the universal Galerkin property) are the essentially coercive forms. In this case, a generalization of the Aubin-Nitsche Theorem leads to optimal a priori estimates in terms of regularity properties of the right-hand side L, as shown by several applications. Finally, a section entitled "Supplement" provides some consequences of our results for the approximation of saddle point problems.
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.
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.
We study the behaviour of solutions to a class of nonlinear degenerate parabolic problems when the data are perturbed. The class includes the Richards equation, Stefan problem and the parabolic $p$-Laplace equation. We show that, up to a subsequence, weak solutions of the perturbed problem converge uniformly-in-time to weak solutions of the original problem as the perturbed data approach the original data. We do not assume uniqueness or additional regularity of the solution. However, when uniqueness is known, our result demonstrates that the weak solution is uniformly temporally stable to perturbations of the data. Beginning with a proof of temporally-uniform, spatially-weak convergence, we strengthen the latter by relating the unknown to an underlying convex structure that emerges naturally from energy estimates on the solution. The double degeneracy --- shown to be equivalent to a maximal monotone operator framework --- is handled with techniques inspired by a classical monotonicity argument and a simple variant of the compensated compactness phenomenon.
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.
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.
We present a new finite volume scheme for anisotropic heterogeneous diffusion problems on unstructured irregular grids, which simultaneously gives an approximation of the solution and of its gradient. In the case of simplicial meshes, the approximate solution is shown to converge to the continuous ones as the size of the mesh tends to 0, and an error estimate is given. In the general case, we propose a slightly modified scheme for which we again prove the convergence, and give an error estimate. An easy implementation method is then proposed, and the efficiency of the scheme is shown on various types of grids.
Finite volume methods for problems involving second order operators with full diffusion matrix can be used thanks to the definition of a discrete gradient for piecewise constant functions on unstructured meshes satisfying an orthogonality condition. This discrete gradient is shown to satisfy a strong convergence property on the interpolation of regular functions, and a weak one on functions bounded for a discrete $H^1$ norm. To highlight the importance of both properties, the convergence of the finite volume scheme on a homogeneous Dirichlet problem with full diffusion matrix is proven, and an error estimate is provided. Numerical tests show the actual accuracy of the method.