A posteriori error estimate for discontinuous Galerkin finite element method on polytopal mesh

In this work, we derive a posteriori error estimates for discontinuous Galerkin finite element method on polytopal mesh. We construct a reliable and efficient a posteriori error estimator on general polygonal or polyhedral meshes. An adaptive algorithm based on the error estimator and DG method is proposed to solve a variety of test problems. Numerical experiments are performed to illustrate the effectiveness of the algorithm.     

Weak Galerkin finite element methods for electric interface model with nonhomogeneous jump conditions

In this paper, the weak Galerkin finite element method (WG-FEM) is applied to a pulsed electric model arising in biological tissue when a biological cell is exposed to an electric field. A fitted WG-FEM is proposed to approximate the voltage of the pulsed electric model across the physical media involving an electric interface (surface membrane), and heterogeneous permittivity and a heterogeneous conductivity. This method uses totally discontinuous functions in approximation space and allows the usage of finite element partitions consisting of general polygonal meshes. Optimal pointwise-in-time error estimates in L²-norm and H¹-norm are shown to hold for the semidiscrete scheme even if the regularity of the solution is low on the whole domain. Furthermore, a fully discrete approximation based on backward Euler scheme is analyzed and related optimal error estimates are derived.     

Virtual element methods for parabolic problems on polygonal meshes

The virtual element method (VEM) is a recent technology that can make use of very general polygonal/polyhedral meshes without the need to integrate complex nonpolynomial functions on the elements and preserving an optimal order of convergence. In this article, we develop for the first time, the VEM for parabolic problems on polygonal meshes, considering time‐dependent diffusion as our model problem. After presenting the scheme, we develop a theoretical analysis and show the practical behavior of the proposed method through a large array of numerical tests. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2110–2134, 2015     

On the Local Uniqueness of Solutions of Variational Inequalities Under H-Differentiability

In this paper, we give some sufficient conditions for the local uniqueness of solutions to nonsmooth variational inequalities where the underlying functions are H-differentiable and the underlying set is a closed convex set/polyhedral set/box/polyhedral cone. We show how the solution of a linearized variational inequality is related to the solution of the variational inequality. These results extend/unify various similar results proved for C ¹ and locally Lipschitzian variational inequality problems. When specialized to the nonlinear complementarity problem, our results extend/unify those of C ² and C ¹ nonlinear complementarity problems.     

Discontinuous Galerkin and mimetic finite difference methods for coupled Stokes–Darcy flows on polygonal and polyhedral grids

We study locally mass conservative approximations of coupled Darcy and Stokes flows on polygonal and polyhedral meshes. The discontinuous Galerkin (DG) finite element method is used in the Stokes region and the mimetic finite difference method is used in the Darcy region. DG finite element spaces are defined on polygonal and polyhedral grids by introducing lifting operators mapping mimetic degrees of freedom to functional spaces. Optimal convergence estimates for the numerical scheme are derived. Results from computational experiments supporting the theory are presented.     

Domain Decomposition Preconditioners for Discontinuous Galerkin Discretizations of Compressible Fluid Flows

In this article we consider the application of Schwarz-type domain decomposition preconditioners to the discontinuous Galerkin finite element approximation of the compressible Navier-Stokes equations. To discretize this system of conservation laws, we exploit the (adjoint consistent) symmetric version of the interior penalty discontinuous Galerkin finite element method. To define the necessary coarse-level solver required for the definition of the proposed preconditioner, we exploit ideas from composite finite element methods, which allow for the definition of finite element schemes on general meshes consisting of polygonal (agglomerated) elements. The practical performance of the proposed preconditioner is demonstrated for a series of viscous test cases in both two- and three-dimensions.     

Weak Galerkin finite element methods for the biharmonic equation on polytopal meshes

A new weak Galerkin (WG) finite element method is introduced and analyzed in this article for the biharmonic equation in its primary form. This method is highly robust and flexible in the element construction by using discontinuous piecewise polynomials on general finite element partitions consisting of polygons or polyhedra of arbitrary shape. The resulting WG finite element formulation is symmetric, positive definite, and parameter‐free. Optimal order error estimates in a discrete H² norm is established for the corresponding WG finite element solutions. Error estimates in the usual L² norm are also derived, yielding a suboptimal order of convergence for the lowest order element and an optimal order of convergence for all high order of elements. Numerical results are presented to confirm the theory of convergence under suitable regularity assumptions. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1003–1029, 2014     

Semi-Discrete and Fully Discrete Weak Galerkin Finite Element Methods for a Quasistatic Maxwell Viscoelastic Model

This paper considers weak Galerkin finite element approximations on polygonal/polyhedral meshes for a quasistatic Maxwell viscoelastic model. The spatial discretization uses piecewise polynomials of degree $k (k ≥ 1)$ for the stress approximation, degree $k+1$ for the velocity approximation, and degree $k$ for the numerical trace of velocity on the inter-element boundaries. The temporal discretization in the fully discrete method adopts a backward Euler difference scheme. We show the existence and uniqueness of the semi-discrete and fully discrete solutions, and derive optimal a priori error estimates. Numerical examples are provided to support the theoretical analysis.     

Post processing of solution and flux for the nodal mimetic finite difference method

We develop and analyze a post processing technique for the family of low‐order mimetic discretizations based on vertex unknowns for the numerical treatment of diffusion problems on unstructured polygonal and polyhedral meshes. The post processing works in two steps. First, from the nodal degrees of freedom, we reconstruct an elemental‐based vector field that approximates the gradient of the exact solution. Second, we solve a local problem for each mesh vertex associated with a scheme degree of freedom to determine a post processed normal flux that is conservative and divergence preserving. Theoretical results and numerical experiments for two‐dimensional (2D) and 3D benchmark problems show optimal convergence rates. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 336–363, 2015     

Discrete maximum principle for Galerkin approximations of the Laplace operator on arbitrary meshes

We derive a nonlinear stabilized Galerkin approximation of the Laplace operator for which we prove a discrete maximum principle on arbitrary meshes and for arbitrary space dimension without resorting to the well-known acute condition or generalizations thereof. We also prove the existence of a discrete solution and discuss the extension of the scheme to convection–diffusion–reaction equations. Finally, we present examples showing that the new scheme cures local minima produced by the standard Galerkin approach while maintaining first-order accuracy in the H¹-norm. To cite this article: E. Burman, A. Ern, C. R. Acad. Sci. Paris, Ser. I 338 (2004).     

THE INITIAL DIRICHLET BOUNDARY VALUE PROBLEM FOR GENERAL SECOND ORDER PARABOLIC SYSTEMS IN NONSMOOTH MANIFOLDS

Abstract

We consider some unilateral boundary value problems in polygonal and polyhedral domains with unilateral transmission conditions. Regularity results in terms of weighted Sobolev spaces are obtained using a penalization technique, similar regularity results for the penalized problems and by showing uniform estimates with respect to the penalization parameter.     

A model of fracture for elliptic problems with flux and solution jumps

A new model of fracture for elliptic problems combining flux and solution jumps as immersed boundary conditions is proposed and proved to be well-posed. An application of this model to the flow in fractured porous media is also proposed including the cases of “impermeable fracture” and “fully permeable fracture” satisfying the so-called “cubic law”, as well as intermediate cases. A finite volume scheme on general polygonal meshes is built to solve such problems. Since no unknown is required at the fracture interface, the scheme is as cheap as standard schemes for the same problems without fault. The convergence of the scheme can be proved to the weak solution of the problem. With weak regularity assumptions, we also establish for the discrete H¹₀ and L² norms some error estimates in $\text{O}\text{(h)}$, where h is the maximum diameter of the control volumes of the mesh. To cite this article: Ph. Angot, C. R. Acad. Sci. Paris, Ser. I 337 (2003).     

Regularity and a priori error analysis of a Ventcel problem in polyhedral domains

We consider the regularity of a mixed boundary value problem for the Laplace operator on a polyhedral domain, where Ventcel boundary conditions are imposed on one face of the polyhedron and Dirichlet boundary conditions are imposed on the complement of that face in the boundary. We establish improved regularity estimates for the trace of the variational solution on the Ventcel face and use them to derive a decomposition of the solution into a regular and a singular part that belongs to suitable weighted Sobolev spaces. This decomposition, in turn, via interpolation estimates both in the interior as well as on the Ventcel face, allows us to perform an a priori error analysis for the finite element approximation of the solution on anisotropic graded meshes. Numerical tests support the theoretical analysis. Copyright © 2016 John Wiley & Sons, Ltd.     

Boundary element methods for Maxwell's equations on non-smooth domains

Summary. Variational boundary integral equations for Maxwell's equations on Lipschitz surfaces in are derived and their well-posedness in the appropriate trace spaces is established. An equivalent, stable mixed reformulation of the system of integral equations is obtained which admits discretization by Galerkin boundary elements based on standard spaces. On polyhedral surfaces, quasioptimal asymptotic convergence of these Galerkin boundary element methods is proved. A sharp regularity result for the surface multipliers on polyhedral boundaries with plane faces is established. Received January 5, 2001 / Revised version received August 6, 2001 / Published online December 18, 2001 Correspondence to: C. Schwab     

On MAC schemes on triangular delaunay meshes,their convergence and application to coupled flow problems

We study the convergence of two generalized marker‐and‐cell covolume schemes for the incompressible Stokes and Navier–Stokes equations introduced by Cavendish, Hall, Nicolaides, and Porsching. The schemes are defined on unstructured triangular Delaunay meshes and exploit the Delaunay–Voronoi duality. The study is motivated by the fact that the related discrete incompressibility condition allows to obtain a discrete maximum principle for the finite volume solution of an advection–diffusion problem coupled to the flow. The convergence theory uses discrete functional analysis and compactness arguments based on recent results for finite volume discretizations for the biharmonic equation. For both schemes, we prove the strong convergence in L² for the velocities and the discrete rotations of the velocities for the Stokes and the Navier–Stokes problem. Further, for one of the schemes, we also prove the strong convergence of the pressure in L². These predictions are confirmed by numerical examples presented in the article. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1397–1424, 2014     

Space-time discontinuos Galerkin method for solving nonstationary convection-diffusion-reaction problems

The paper presents the theory of the discontinuous Galerkin finite element method for the space-time discretization of a linear nonstationary convection-diffusion-reaction initial-boundary value problem. The discontinuous Galerkin method is applied separately in space and time using, in general, different nonconforming space grids on different time levels and different polynomial degrees p and q in space and time discretization, respectively. In the space discretization the nonsymmetric interior and boundary penalty approximation of diffusion terms is used. The paper is concerned with the proof of error estimates in “L ²(L ²)”-and “ ”-norms, where ɛ ⩾ 0 is the diffusion coefficient. Using special interpolation theorems for the space as well as time discretization, we find that under some assumptions on the shape regularity of the meshes and a certain regularity of the exact solution, the errors are of order O(h ^p + τ ^q ). The estimates hold true even in the hyperbolic case when ɛ = 0.     

A space‐time continuous Galerkin method with mesh modification for viscoelastic wave equations

In this article, we consider the space‐time continuous Galerkin (STCG) method for the viscoelastic wave equations. It allows variable temporal step‐sizes, and the changing of the spatial grids in two adjacent time levels. The existence, uniqueness, and stability of the approximate solutions are demonstrated and the error estimates with global and local spatial mesh sizes in norm are derived without any restrictive assumptions on the space‐time meshes. If the meshes in each time level satisfy some reasonable assumptions, then we can get the optimal order error estimates both in time and space. Finally, we give a numerical example on unstructured meshes to confirm the theoretical findings. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1183–1207, 2017     

Error estimates for the numerical approximation of a quaslinear Neumann problem under minimal regularity of the data

The finite element based approximation of a quasilinear elliptic equation of non monotone type with Neumann boundary conditions is studied. Minimal regularity assumptions on the data are imposed. The consideration is restricted to polygonal domains of dimension two and polyhedral domains of dimension three. Finite elements of degree k ≥ 1 are used to approximate the equation. Error estimates are established in the L ²(Ω) and H ¹(Ω) norms for convex and non-convex domains. The issue of uniqueness of a solution to the approximate discrete equation is also addressed.     

The hp-Version of the Boundary Element Method in ℝ3 The Basic Approximation Results

This paper deals with the basic approximation properties of the h–p version of the boundary element method (BEM) in ℝ³. We extend the results on the exponential convergence of the h–p version of the boundary element method on geometric meshes from problems in polygonal domains to problems in polyhedral domains. In 2D elliptic boundary value problems the solutions have only corner singularities whereas in 3D problems they contain additional edge and corner-edge singularities. The solutions of the corresponding boundary integral equations inherit those singularities. The detailed investigations in our analysis take care of the various types of those singularities. While edge singularities can be analysed using standard one-dimensional approximation results the corner-edge singularities demand a new analysis. © 1997 by B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Graded Mesh Refinement and Error Estimates for Finite Element Solutions of Elliptic Boundary Value Problems in Non-smooth Domains

This paper is concerned with the effective numerical treatment of elliptic boundary value problems when the solutions contain singularities. The paper deals first with the theory of problems of this type in the context of weighted Sobolev spaces and covers problems in domains with conical vertices and non-intersecting edges, as well as polyhedral domains with Lipschitz boundaries. Finite element schemes on graded meshes for second-order problems in polygonal/polyhedral domains are then proposed for problems with the above singularities. These schemes exhibit optimal convergence rates with decreasing mesh size. Finally, we describe numerical experiments which demonstrate the efficiency of our technique in terms of ‘actual’ errors for specific (finite) mesh sizes in addition to the asymptotic rates of convergence.