A posteriori error estimates for the mixed finite element method with Lagrange multipliers

We consider the mixed finite element method with Lagrange multipliers as applied to second‐order elliptic equations in divergence form with mixed boundary conditions. The corresponding Galerkin scheme is defined by using Raviart‐Thomas spaces. We develop a posteriori error analyses yielding a reliable and efficient estimate based on residuals, and a reliable and quasi‐efficient estimate based on local problems, respectively. Several numerical results illustrate the suitability of these a posteriori estimates for the adaptive computation of the discrete solutions. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Least‐squares mixed finite element methods for the RLW equations

A least‐squares mixed finite element (LSMFE) schemes are formulated to solve the 1D regularized long wave (RLW) equations and the convergence is discussed. The L² error estimates of LSMFE methods for RLW equations under the standard regularity assumption on the finite element partition are given.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

On the weighting method for mixed least squares–total least squares problems

It is well known that the standard algorithm for the mixed least squares–total least squares (MTLS) problem uses the QR factorization to reduce the original problem into a standard total least squares problem with smaller size, which can be solved based on the singular value decomposition (SVD). In this paper, the MTLS problem is proven to be closely related to a weighted total least squares problem with its error‐free columns multiplied by a large weighting factor. A criterion for choosing the weighting factor is given; and for the sake of stability in solving the MTLS problem, the Cholesky factorization‐based inverse (Cho‐INV) iteration and Rayleigh quotient iteration are also considered. For large‐scale MTLS problems, numerical tests show that Cho‐INV is superior to the standard QR‐SVD method, especially for the case with big gap between the desired and undesired singular values and the case when the coefficient matrix has much more error‐contaminated columns. Rayleigh quotient iteration behaves more efficient than QR‐SVD for most cases and fails occasionally, and in some cases, it converges much faster than Cho‐INV but still less efficient due to its higher computation cost.     

Analysis of expanded mixed methods for fourth-order elliptic problems

The recently proposed expanded mixed formulation for numerical solution of second-order elliptic problems is here extended to fourth-order elliptic problems. This expanded formulation for the differential problems under consideration differs from the classical formulation in that three variables are treated, i.e., the displacement, the stress, and the moment tensors. It works for the case where the coefficient of the differential equations is small and does not need to be inverted, or for the case in which the stress tensor of the equations does not need to be symmetric. Based on this new formulation, various mixed finite elements for fourth-order problems are considered; error estimates of quasi-optimal or optimal order depending upon the mixed elements are derived. Implementation techniques for solving the linear system arising from these expanded mixed methods are discussed, and numerical results are presented. © 1997 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 483–503, 1997     

A mixed finite volume method for elliptic problems

We derive a novel finite volume method for the elliptic equation, using the framework of mixed finite element methods to discretize the pressure and velocities on two different grids (covolumes), triangular (tetrahedral) mesh and control volume mesh. The new discretization is defined for tensor diffusion coefficient and well suited for heterogeneous media. When the control volumes are created by connecting the center of gravity of each triangle to the midpoints of its edges, we show that the discretization is stable and first order accurate for both scalar and vector unknowns. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

A posteriori error analysis of nonconforming finite volume elements for general second‐order elliptic PDEs

In this article, we study the a posteriori H¹ and L² error estimates for Crouzeix‐Raviart nonconforming finite volume element discretization of general second‐order elliptic problems in ?². The error estimators yield global upper and local lower bounds. Finally, numerical experiments are performed to illustrate the theoretical findings. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Error estimates for finite volume element methods for general second‐order elliptic problems

We treat the finite volume element method (FVE) for solving general second order elliptic problems as a perturbation of the linear finite element method (FEM), and obtain the optimal H¹ error estimate, H¹ superconvergence and L^p (1 < p ≤ ∞) error estimates between the solution of the FVE and that of the FEM. In particular, the superconvergence result does not require any extra assumptions on the mesh except quasi‐uniform. Thus the error estimates of the FVE can be derived by the standard error estimates of the FEM. Moreover we consider the effects of numerical integration and prove that the use of barycenter quadrature rule does not decrease the convergence orders of the FVE. The results of this article reveal that the FVE is in close relationship with the FEM. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 693–708, 2003.     

Discontinuous Galerkin finite volume element methods for second‐order linear elliptic problems

In this article, a one parameter family of discontinuous Galerkin finite volume element methods for approximating the solution of a class of second‐order linear elliptic problems is discussed. Optimal error estimates in L² and broken H¹‐ norms are derived. Numerical results confirm the theoretical order of convergences. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A weak Galerkin method for second order elliptic problems with polynomial reduction

The second order elliptic equation, which is also know as the diffusion-convection equation, is of great interest in many branches of physics and industry. In this paper, we use the weak Galerkin finite element method to study the general second order elliptic equation. A weak Galerkin finite element method is proposed and analyzed. This scheme features piecewise polynomials of degree $k\geq 1$ on each element and piecewise polynomials of degree $k-1\geq 0$ on each edge or face of the element. Error estimates of optimal order of convergence rate are established in both discrete $H^1$ and standard $L^2$ norm. The paper also presents some numerical experiments to verify the efficiency of the method.     

A weak Galerkin finite element method for the second order elliptic problems with mixed boundary conditions

In this paper, a weak Galerkin finite element method is proposed and analyzed for the second-order elliptic equation with mixed boundary conditions. Optimal order error estimates are established in both discrete $H^1$ norm and the standard $L^2$ norm for the corresponding WG approximations. The numerical experiments are presented to verify the efficiency of the method.     

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     

On the mixed finite element method with Lagrange multipliers

In this note we analyze a modified mixed finite element method for second‐order elliptic equations in divergence form. As a model we consider the Poisson problem with mixed boundary conditions in a polygonal domain of R ². The Neumann (essential) condition is imposed here in a weak sense, which yields the introduction of a Lagrange multiplier given by the trace of the solution on the corresponding boundary. This approach allows to handle nonhomogeneous Neumann boundary conditions, theoretically and computationally, in an alternative and usually easier way. Then we utilize the classical Babu?ka‐Brezzi theory to show that the resulting mixed variational formulation is well posed. In addition, we use Raviart‐Thomas spaces to define the associated finite element method and, applying some elliptic regularity results, we prove the stability, unique solvability, and convergence of this discrete scheme, under appropriate assumptions on the mesh sizes. Finally, we provide numerical results illustrating the performance of the algorithm for smooth and singular problems. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 192–210, 2003     

A mixed and discontinuous Galerkin finite volume element method for incompressible miscible displacement problems in porous media

The incompressible miscible displacement problem in porous media is modeled by a coupled system of two nonlinear partial differential equations, the pressure‐velocity equation and the concentration equation. In this article, we present a mixed finite volume element method for the approximation of pressure‐velocity equation and a discontinuous Galerkin finite volume element method for the concentration equation. A priori error estimates in L^∞(L²) are derived for velocity, pressure, and concentration. Numerical results are presented to substantiate the validity of the theoretical results. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

Mixed finite volume method for nonlinear elliptic problems

In this article we construct and analyze a mixed finite volume method for second‐order nonlinear elliptic problems employing H(div; Ω)‐conforming approximations for the vector variable and completely discontinuous approximations for the scalar variable. The main attractive feature of our method is that, although the vector variable is H(div; Ω)‐conforming, one can eliminate it in a local manner to obtain a discontinuous Galerkin method for the scalar variable. Optimal error estimates will be established for both vector and scalar variables. We also present a fully discrete version of this method that is more convenient for computational purposes. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A general mixed covolume framework for constructing conservative schemes for elliptic problems

We present a general framework for the finite volume or covolume schemes developed for second order elliptic problems in mixed form, i.e., written as first order systems. We connect these schemes to standard mixed finite element methods via a one-to-one transfer operator between trial and test spaces. In the nonsymmetric case (convection-diffusion equation) we show one-half order convergence rate for the flux variable which is approximated either by the lowest order Raviart-Thomas space or by its image in the space of discontinuous piecewise constants. In the symmetric case (diffusion equation) a first order convergence rate is obtained for both the state variable (e.g., concentration) and its flux. Numerical experiments are included.

     

High‐order difference schemes for 2D elliptic and parabolic problems with mixed derivatives

We propose a 9‐point fourth‐order finite difference scheme for 2D elliptic problems with a mixed derivative and variable coefficients. The same approach is extended to derive a class of two‐level high‐order compact schemes with weighted time discretization for solving 2D parabolic problems with a mixed derivative. The schemes are fourth‐order accurate in space and second‐ or lower‐order accurate in time depending on the choice of a weighted average parameter μ. Unconditional stability is proved for 0.5 ≤ μ ≤ 1, and numerical experiments supporting our theoretical analysis and confirming the high‐order accuracy of the schemes are presented. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 366–378, 2007     

A C0‐weak Galerkin finite element method for fourth‐order elliptic problems

This article proposes and analyzes a C⁰‐weak Galerkin (WG) finite element method for solving the biharmonic equation in two‐dimensional and three‐dimensional. The new WG method uses continuous piecewise‐polynomial approximations of degree for the unknown u and discontinuous piecewise‐polynomial approximations of degree k for the trace of on the interelement boundaries. Optimal error estimates are obtained in H², H¹, and L² norms. Numerical experiments illustrate and confirm the theoretical results. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1090–1104, 2016     

Some nonconforming mixed box schemes for elliptic problems

In this article, we introduce three schemes for the Poisson problem in 2D on triangular meshes, generalizing the FVbox scheme introduced by Courbet and Croisille [1]. In this kind of scheme, the approximation is performed on the mixed form of the problem, but contrary to the standard mixed method, with a pair of trial spaces different from the pair of test spaces. The latter is made of Galerkin‐discontinuous spaces on a unique mesh. The first scheme uses as trial spaces the P¹ nonconforming space of Crouzeix‐Raviart both for u and for the flux p = ?u. In the two others, the quadratic nonconforming space of Fortin and Soulie is used. An important feature of all these schemes is that they are equivalent to a first scheme in u only and an explicit representation formula for the flux p = ?u. The numerical analysis of the schemes is performed using this property. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 355–373, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10003     

Optimal error estimates of the local discontinuous Galerkin method for nonlinear second‐order elliptic problems on Cartesian grids

In this paper, we study the local discontinuous Galerkin (LDG) methods for two‐dimensional nonlinear second‐order elliptic problems of the type u_xx + u_yy = f(x, y, u, u_x, u_y) , in a rectangular region Ω with classical boundary conditions on the boundary of Ω . Convergence properties for the solution and for the auxiliary variable that approximates its gradient are established. More specifically, we use the duality argument to prove that the errors between the LDG solutions and the exact solutions in the L² norm achieve optimal (p + 1)th order convergence, when tensor product polynomials of degree at most p are used. Moreover, we prove that the gradient of the LDG solution is superclose with order p + 1 toward the gradient of Gauss–Radau projection of the exact solution. The results are valid in two space dimensions on Cartesian meshes using tensor product polynomials of degree p ≥ 1 , and for both mixed Dirichlet–Neumann and periodic boundary conditions. Preliminary numerical experiments indicate that our theoretical findings are optimal.     

Finite volume element method for monotone nonlinear elliptic problems

In this article, we consider the finite volume element method for the monotone nonlinear second‐order elliptic boundary value problems. With the assumptions which guarantee that the corresponding operator is strongly monotone and Lipschitz‐continuous, and with the minimal regularity assumption on the exact solution, that is, u∈H¹(Ω), we show that the finite volume element method has a unique solution, and the finite volume element approximation is uniformly convergent with respect to the H¹ ‐norm. If u∈H^1+ε(Ω),0 < ε ≤ 1, we develop the optimal convergence rate \begin{align*}\mathcal{O}(h^{\epsilon})\end{align*} in the H¹ ‐norm. Moreover, we propose a natural and computationally easy residual‐based H¹ ‐norm a posteriori error estimator and establish the global upper bound and local lower bounds on the error. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013