A discontinuous finite volume element method for second‐order elliptic problems

In this article, we propose a new discontinuous finite volume element (DFVE) method for the second‐order elliptic problems. We treat the DFVE method as a perturbation of the interior penalty method and get a superapproximation estimate in a mesh dependent norm between the solution of the DFVE method and that of the interior penalty method. This reveals that the DFVE method is much closer to the interior penalty method than we have known. By using this superapproximation estimate, we can easily get the optimal order error estimates in the L² ‐norm and in the maximum norms of the DFVE method.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 425–440, 2012     

A stabilizer-free C0 weak Galerkin method for the biharmonic equations

In this article, we present and analyze a stabilizer-free C⁰ weak Galerkin(SF-C0WG) method for solving the biharmonic problem. The SF-C0WG method is formulated in terms of cell unknowns which are C⁰ continuous piecewise polynomials of degree k + 2 with k≥0 and in terms of face unknowns which are discontinuous piecewise polynomials of degree k + 1. The formulation of this SF-C0WG method is without the stabilized or penalty term and is as simple as the C¹ conformin...     

The fictitious domain method with L2‐penalty for the Stokes problem with the Dirichlet boundary condition

We consider the fictitious domain method with L²‐penalty for the Stokes problem with the Dirichlet boundary condition. First, we investigate the error estimates for the penalty method at the continuous level. We obtain the convergence of order in H¹‐norm for the velocity and in L²‐norm for the pressure, where is the penalty parameter. The L²‐norm error estimate for the velocity is upgraded to . Moreover, we derive the a priori estimates depending on for the solution of the penalty problem. Next, we apply the finite element approximation to the penalty problem using the P1/P1 element with stabilization. For the discrete penalty problem, we prove the error estimate in H¹‐norm for the velocity and in L²‐norm for the pressure, where h denotes the discretization parameter. For the velocity in L²‐norm, the convergence rate is improved to . The theoretical results are verified by the numerical experiments.     

Strong convergence of the discontinuous Galerkin scheme for the low regularity miscible displacement equations

Strong convergence of the numerical solution to a weak solution is proved for a nonlinear coupled flow and transport problem arising in porous media. The method combines a mixed finite element method for the pressure and velocity with an interior penalty discontinuous Galerkin method in space for the concentration. Using functional tools specific to broken Sobolev spaces, the convergence of the broken gradient of the numerical concentration to the weak solution is obtained in the L² norm. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 489–513, 2017     

A characteristics-mixed covolume method for a convection-dominated transport problem

In this paper, we propose a characteristics-mixed covolume method for approximating the solution to a convection dominated transport problem. The method is a combination of characteristic approximation to handle the convection term in time and mixed covolume method spatial approximation to deal with the diffusion term. The velocity and press are approximated by the lowest order Raviart-Thomas mixed finite element space on rectangles. The projection of a mixed covolume element is introduced. We prove its first order optimal rate of convergence for the approximate velocities in the L² norm as well as for the approximate pressures in the L² norm.     

Multigrid methods for the symmetric interior penalty method on graded meshes

The symmetric interior penalty (SIP) method on graded meshes and its fast solution by multigrid methods are studied in this paper. We obtain quasi‐optimal error estimates in both the energy norm and the L₂ norm for the SIP method, and prove uniform convergence of the W‐cycle multigrid algorithm for the resulting discrete problem. The performance of these methods is illustrated by numerical results. Copyright © 2009 John Wiley & Sons, Ltd.     

Interior penalty method for the indefinite time-harmonic Maxwell equations

Summary In this paper, we introduce and analyze the interior penalty discontinuous Galerkin method for the numerical discretization of the indefinite time-harmonic Maxwell equations in the high-frequency regime. Based on suitable duality arguments, we derive a-priori error bounds in the energy norm and the L²-norm. In particular, the error in the energy norm is shown to converge with the optimal order (h^min{s,}) with respect to the mesh size h, the polynomial degree , and the regularity exponent s of the analytical solution. Under additional regularity assumptions, the L²-error is shown to converge with the optimal order (h⁺¹). The theoretical results are confirmed in a series of numerical experiments.Supported by the EPSRC (Grant GR/R76615).Supported by the Swiss National Science Foundation under project 21-068126.02.Supported in part by the Natural Sciences and Engineering Council of Canada.     

A priori error estimates of an extrapolated space‐time discontinuous galerkin method for nonlinear convection‐diffusion problems

We deal with the numerical solution of a scalar nonstationary nonlinear convection‐diffusion equation. We employ a combination of the discontinuous Galerkin finite element (DGFE) method for the space as well as time discretization. The linear diffusive and penalty terms are treated implicitly whereas the nonlinear convective term is treated by a special higher order explicit extrapolation from the previous time step, which leads to the necessity to solve only a linear algebraic problem at each time step. We analyse this scheme and derive a priori asymptotic error estimates in the L^∞(L²) –norm and the L²(H¹) –seminorm with respect to the mesh size h and time step τ. Finally, we present an efficient solution strategy and numerical examples verifying the theoretical results. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1456–1482, 2010     

Optimal L ∞(L 2)-Error Estimates for the DG Method Applied to Nonlinear Convection–Diffusion Problems with Nonlinear Diffusion

This article is concerned with the analysis of the discontinuous Galerkin finite element method (DGFEM) applied to the space semidiscretization of a nonstationary convection–diffusion problem with nonlinear convection and nonlinear diffusion. Optimal estimates in the L ^∞(L ²)-norm are derived for the symmetric interior penalty (SIPG) scheme in two dimensions. The error analysis is carried out for nonconforming triangular meshes under the assumption that the exact solution of the problem and the solution of a linearized elliptic dual problem are sufficiently regular.     

Continuous interior penalty finite element methods for Sobolev equations with convection‐dominated term

We consider implicit and semi‐implicit time‐stepping methods for continuous interior penalty (CIP) finite element approximations of Sobolev equations with convection‐dominated term. Stability is obtained by adding an interior penalty term giving L² ‐control of the jump of the gradient over element faces. Several $\cal {A}$ ‐stable time‐stepping methods are analyzed and shown to be unconditionally stable and optimally convergent. We show that the contribution from the gradient jumps leading to an extended matrix pattern may be extrapolated from previous time steps, and hence handled explicitly without loss of stability and accuracy. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

A nonconforming finite element method for a two-dimensional curl–curl and grad-div problem

A numerical method for a two-dimensional curl–curl and grad-div problem is studied in this paper. It is based on a discretization using weakly continuous P ₁ vector fields and includes two consistency terms involving the jumps of the vector fields across element boundaries. Optimal convergence rates (up to an arbitrary positive ) in both the energy norm and the L ₂ norm are established on graded meshes. The theoretical results are confirmed by numerical experiments. The work of the first author was supported in part by the National Science Foundation under Grant No. DMS-03-11790 and by the Humboldt Foundation through her Humboldt Research Award. The work of the third author was supported in part by the National Science Foundation under Grant No. DMS-06-52481.     

A CFL‐free explicit characteristic interior penalty scheme for linear advection‐reaction equations

We develop a CFL‐free, explicit characteristic interior penalty scheme (CHIPS) for one‐dimensional first‐order advection‐reaction equations by combining a Eulerian‐Lagrangian approach with a discontinuous Galerkin framework. The CHIPS method retains the numerical advantages of the discontinuous Galerkin methods as well as characteristic methods. An optimal‐order error estimate in the L² norm for the CHIPS method is derived and numerical experiments are presented to confirm the theoretical estimates. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

A splitting positive definite mixed element method for miscible displacement of compressible flow in porous media

A miscible displacement of one compressible fluid by another in a porous medium is governed by a nonlinear parabolic system. A new mixed finite element method, in which the mixed element system is symmetric positive definite and the flux equation is separated from pressure equation, is introduced to solve the pressure equation of parabolic type, and a standard Galerkin method is used to treat the convection‐diffusion equation of concentration of one of the fluids. The convergence of the approximate solution with an optimal accuracy in L²‐norm is proved. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 229–249, 2001     

Fully discrete interior penalty discontinuous Galerkin methods for nonlinear parabolic equations

In this article, we investigate interior penalty discontinuous Galerkin (IPDG) methods for solving a class of two‐dimensional nonlinear parabolic equations. For semi‐discrete IPDG schemes on a quasi‐uniform family of meshes, we obtain a priori bounds on solutions measured in the L² norm and in the broken Sobolev norm. The fully discrete IPDG schemes considered are based on the approximation by forward Euler difference in time and broken Sobolev space. Under a restriction related to the mesh size and time step, an hp ‐version of an a priori l^∞(L²) and l²(H¹) error estimate is derived and numerical experiments are presented.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 288–311, 2012     

A Crank‐Nicolson Petrov‐Galerkin method with quadrature for semi‐linear parabolic problems

We propose and analyze an application of a fully discrete C² spline quadrature Petrov‐Galerkin method for spatial discretization of semi‐linear parabolic initial‐boundary value problems on rectangular domains. We prove second order in time and optimal order H¹ norm convergence in space for the extrapolated Crank‐Nicolson quadrature Petrov‐Galerkin scheme. We demonstrate numerically both L² and H¹ norm optimal order convergence of the scheme even if the nonlinear source term is not smooth. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

Stability and convergence of two‐grid Crank‐Nicolson extrapolation scheme for the time‐dependent natural convection equations

In this paper, we consider the Crank‐Nicolson extrapolation scheme for the 2D/3D unsteady natural convection problem. Our numerical scheme includes the implicit Crank‐Nicolson scheme for linear terms and the recursive linear method for nonlinear terms. Standard Galerkin finite element method is used to approximate the spatial discretization. Stability and optimal error estimates are provided for the numerical solutions. Furthermore, a fully discrete two‐grid Crank‐Nicolson extrapolation scheme is developed, the corresponding stability and convergence results are derived for the approximate solutions. Comparison from aspects of the theoretical results and computational efficiency, the two‐grid Crank‐Nicolson extrapolation scheme has the same order as the one grid method for velocity and temperature in H¹‐norm and for pressure in L²‐norm. However, the two‐grid scheme involves much less work than one grid method. Finally, some numerical examples are provided to verify the established theoretical results and illustrate the performances of the developed numerical schemes.     

Error analysis of the L2 least‐squares finite element method for incompressible inviscid rotational flows

In this article we analyze the L² least‐squares finite element approximations to the incompressible inviscid rotational flow problem, which is recast into the velocity‐vorticity‐pressure formulation. The least‐squares functional is defined in terms of the sum of the squared L² norms of the residual equations over a suitable product function space. We first derive a coercivity type a priori estimate for the first‐order system problem that will play the crucial role in the error analysis. We then show that the method exhibits an optimal rate of convergence in the H¹ norm for velocity and pressure and a suboptimal rate of convergence in the L² norm for vorticity. A numerical example in two dimensions is presented, which confirms the theoretical error estimates. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

A priori error analysis of a discontinuous Galerkin approximation for a kind of compressible miscible displacement problems

A kind of compressible miscible displacement problems which include molecular diffusion and dispersion in porous media are investigated.A symmetric interior penalty discontinuous Galerkin (SIPG) method is applied to the coupled system of flow and transport.Using the induction hypotheses instead of the cut-off operator and the interpolation projection properties,a priori hp error estimates are presented.The error bounds in L2(H1) norm for concentration and in L∞(L2) norm for velocity are optimal in h and suboptimal in p with a loss of power 1/2.     

Adaptive finite element solution of the porous medium equation in pressure formulation

The pressure formulation of the porous medium equation has been commonly used in theoretical studies due to its much better regularities than the original formulation. The goal here is to study its use in the adaptive moving mesh finite element solution. The free boundary is traced explicitly through Darcy's law. The method is shown numerically second‐order in space and first‐order in time in the pressure variable. Moreover, the convergence order of the error in the location of the free boundary is almost second‐order in the maximum norm. However, numerical results also show that the convergence order in the original variable stays between first‐order and second‐order in L¹ norm or between 0.5th‐order and first‐order in L² norm. Nevertheless, the current method can offer some advantages over numerical methods based on the original formulation for situations with large exponents or when a more accurate location of the free boundary is desired.     

A supercloseness result for the discontinuous Galerkin stabilization of convection–diffusion problems on Shishkin meshes

We consider a convection–diffusion problem with Dirichlet boundary conditions posed on a unit square. The problem is discretized using a combination of the standard Galerkin FEM and an h–version of the nonsymmetric discontinuous Galerkin FEM with interior penalties on a layer–adapted mesh with linear/bilinear elements. With specially chosen penalty parameters for edges from the coarse part of the mesh, we prove uniform convergence (in the perturbation parameter) in an associated norm. In the same norm we also establish a supercloseness result. Numerical tests support our theoretical estimates.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007