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 element analysis of a thermally nonlinear coupled problem

In Loula and Zhou [Comput Appl Math 20 (2001), 321–339], a thermally coupled nonlinear elliptic system modeling a large class of engineering problems was considered, and some mathematical and numerical analyses (C⁰ Lagrangian finite elements combined with a fixed point algorithm) were given. To continue our work, we propose in this article a mixed method for the potential equation and present the corresponding analyses and numerical implementations. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Analysis of a combined mixed finite element and discontinuous Galerkin method for incompressible two‐phase flow in porous media

We analyze a combined method consisting of the mixed finite element method for pressure equation and the discontinuous Galerkin method for saturation equation for the coupled system of incompressible two‐phase flow in porous media. The existence and uniqueness of numerical solutions are established under proper conditions by using a constructive approach. Optimal error estimates in L²(H¹) for saturation and in L ^∞ (H(div)) for velocity are derived. Copyright © 2013 John Wiley & Sons, Ltd.     

An approximation of incompressible miscible displacement in porous media by mixed finite elements and symmetric finite volume element method of characteristics

A nonlinear system of two coupled partial differential equations models miscible displacement of one incompressible fluid by another in a porous medium. A sequential implicit time‐stepping procedure is defined, in which the pressure and Darcy velocity of the mixture are approximated by a mixed finite element method and the concentration is approximated by a combination of a modified symmetric finite volume element method and the method of characteristics. Optimal order convergence in H¹ and in L² are proved for full discrete schemes. Finally, some numerical experiments are presented. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A two-grid finite element method for the Allen-Cahn equation with the logarithmic potential

In this paper, we present a two-grid finite element method for the Allen-Cahn equation with the logarithmic potential. This method consists of two steps. In the first step, based on a fully implicit finite element method, the Allen-Cahn equation is solved on a coarse grid with mesh size H. In the second step, a linearized system whose nonlinear term is replaced by the value of the first step is solved on a fine grid with mesh size h. We give the energy stabilities of the traditional finite element method and the two-grid finite element method. The optimal convergence order of the two-grid finite element method in H¹ norm is achieved when the mesh sizes satisfy h = O(H²). Numerical examples are given to demonstrate the validity of the proposed scheme. The results show that the two-grid method can save the CPU time while keeping the same convergence rate.     

Fully discrete T‐ψ finite element method to solve a nonlinear induction hardening problem

We study an induction hardening model described by Maxwell's equations coupled with a heat equation. The magnetic induction field is assumed a nonlinear constitutional relation and the electric conductivity is temperature‐dependent. The T ‐ψ method is to transform Maxwell's equations to the vector–scalar potential formulations and to solve the potentials by means of the finite element method. In this article, we present a fully discrete T ‐ψ finite element scheme for this nonlinear coupled problem and discuss its solvability. We prove that the discrete solution converges to a weak solution of the continuous problem. Finally, we conclude with two numerical experiments for the coupled system.     

A weak Galerkin mixed finite element method for the Helmholtz equation with large wave numbers

In this article, a new weak Galerkin mixed finite element method is introduced and analyzed for the Helmholtz equation with large wave numbers. The stability and well‐posedness of the method are established for any wave number k without mesh size constraint. Allowing the use of discontinuous approximating functions makes weak Galerkin mixed method highly flexible in term of little restrictions on approximations and meshes. In the weak Galerkin mixed finite element formulation, approximation functions can be piecewise polynomials with different degrees on different elements and meshes can consist elements with different shapes. Suboptimal order error estimates in both discrete H¹ and L² norms are established for the weak Galerkin mixed finite element solutions. Numerical examples are tested to support the theory.     

Stability of a combined finite element ‐ finite volume discretization of convection‐diffusion equations

We consider a time‐dependent and a stationary convection‐diffusion equation. These equations are approximated by a combined finite element – finite volume method: the diffusion term is discretized by Crouzeix‐Raviart piecewise linear finite elements on a triangular grid, and the convection term by upwind barycentric finite volumes. In the nonstationary case, we use an implicit Euler approach for time discretization. This scheme is shown to be L²‐stable uniformly with respect to the diffusion coefficient. In addition, it turns out that stability is unconditional in the time‐dependent case. These results hold if the underlying grid satisfies a condition that is fulfilled, for example, by some structured meshes. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 402–424, 2012     

Quasi-least-squares finite element method for steady flow and heat transfer with system rotation

Two quasi-least-squares finite element schemes based on L ² inner product are proposed to solve a steady Navier–Stokes equations, coupled to the energy equation by the Boussinesq approximation and augmented by a Coriolis forcing term to account for system rotation. The resulting nonlinear systems are linearized around a characteristic state, resulting in linearized least-squares models that yield algebraic systems with symmetric positive definite coefficient matrices. Existence of solutions are investigated and a priori error estimates are obtained. The performance of the formulation is illustrated by using a direct iteration procedure to treat the nonlinearities and shown theoretical convergent rate for general initial guess.     

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     

An upwind center difference parallel method and numerical analysis for the displacement problem with moving boundary

A nonlinear coupled mathematical system of two‐phase seepage flow displacement is discussed in this paper including an elliptic equation for the pressure and a convection‐dominated diffusion equation for the saturation. In fact, the boundary of an underground region where the fluid flows through is nonstationary. So a moving boundary should be considered. The saturation equation is convection‐dominated, therefore the method of upwind finite difference is introduced for the accurate computation. The upwind approximation could eliminate numerical oscillation and strong stability is shown. Since the computational work of saturation is larger than the pressure, the authors apply a parallel method, decomposing the whole domain into several nonoverlapping subdomains, to simplify the computation. A domain decomposition method coupled with upwind differences is presented for the saturation. The pressure equation is discretized by a five‐point center finite difference method. By using a transformation and defining new inner products and norms, error estimates in l² norm is discussed. Finally, two experimental tests are given to illustrate the efficiency and accuracy of the parallel algorithm.     

Weak Galerkin finite element methods for a fourth order parabolic equation

This paper is devoted to a newly developed weak Galerkin finite element method with the stabilization term for a linear fourth order parabolic equation, where weakly defined Laplacian operator over discontinuous functions is introduced. Priori estimates are developed and analyzed in L² and an H² type norm for both semi‐discrete and fully discrete schemes. And finally, numerical examples are provided to confirm the theoretical results.     

High accuracy error estimates of a Galerkin finite element method for nonlinear time fractional diffusion equation

In this work, an effective and fast finite element numerical method with high-order accuracy is discussed for solving a nonlinear time fractional diffusion equation. A two-level linearized finite element scheme is constructed and a temporal–spatial error splitting argument is established to split the error into two parts, that is, the temporal error and the spatial error. Based on the regularity of the time discrete system, the temporal error estimate is derived. Using the property of the Ritz projection operator, the spatial error is deduced. Unconditional superclose result in H¹-norm is obtained, with no additional regularity assumption about the exact solution of the problem considered. Then the global superconvergence error estimate is obtained through the interpolated postprocessing technique. In order to reduce storage and computation time, a fast finite element method evaluation scheme for solving the nonlinear time fractional diffusion equation is developed. To confirm the theoretical error analysis, some numerical results are provided.     

Error estimates for mixed finite element approximations of the viscoelasticity wave equation

This paper studies mixed finite element approximations to the solution of the viscoelasticity wave equation. Two new transformations are introduced and a corresponding system of first‐order differential‐integral equations is derived. The semi‐discrete and full‐discrete mixed finite element methods are then proposed for the problem based on the Raviart–Thomas–Nedelec spaces. The optimal error estimates in L²‐norm are obtained for the semi‐discrete and full‐discrete mixed approximations of the general viscoelasticity wave equation. Copyright © 2004 John Wiley & Sons, Ltd.     

Composite finite element approximation for nonlinear parabolic problems in nonconvex polygonal domains

We study a new class of finite elements so‐called composite finite elements (CFEs), introduced earlier by Hackbusch and Sauter, Numer. Math., 1997; 75:447‐472, for the approximation of nonlinear parabolic equation in a nonconvex polygonal domain. A two‐scale CFE discretization is used for the space discretizations, where the coarse‐scale grid discretized the domain at an appropriate distance from the boundary and the fine‐scale grid is used to resolve the boundary. A continuous, piecewise linear CFE space is employed for the spatially semidiscrete finite element approximation and the temporal discretizations is based on modified linearized backward Euler scheme. We derive almost optimal‐order convergence in space and optimal order in time for the CFE method in the L^∞(L²) norm. Numerical experiment is carried out for an L‐shaped domain to illustrate our theoretical findings.     

Error estimates for a finite element-finite volume discretization of convection-diffusion equations

We consider a time-dependent linear convection-diffusion equation. This equation is approximated by a combined finite element-finite volume method: the diffusion term is discretized by Crouzeix-Raviart piecewise linear finite elements, and the convection term by upwind barycentric finite volumes on a triangular grid. An implicit Euler approach is used for time discretization. It is shown that the error associated with this scheme, measured by a discrete L^∞-L²- and L²-H¹-norm, respectively, decays linearly with the mesh size and the time step. This result holds without any link between mesh size and time step. The dependence of the corresponding error bound on the diffusion coefficient is completely explicit.     

Nonconforming finite element methods with subgrid viscosity applied to advection‐diffusion‐reaction equations

A nonconforming (Crouzeix–Raviart) finite element method with subgrid viscosity is analyzed to approximate advection‐diffusion‐reaction equations. The error estimates are quasi‐optimal in the sense that keeping the Péclet number fixed, the estimates are suboptimal of order in the mesh size for the L²‐norm and optimal for the advective derivative on quasi‐uniform meshes. The method is also reformulated as a finite volume box scheme providing a reconstruction formula for the diffusive flux with local conservation properties. Numerical results are presented to illustrate the error analysis. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Superconvergence of finite volume element method for a nonlinear elliptic problem

In this article, we consider the finite volume element method for the second‐order nonlinear elliptic problem and obtain the H¹ and W^1,∞ superconvergence estimates between the solution of the finite volume element method and that of the finite element method, which reveal that the finite volume element method is in close relationship with the finite element method. With these superconvergence estimates, we establish the L^p and W^1,p (2 < p ≤ ∞) error estimates for the finite volume element method for the second‐order nonlinear elliptic problem. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Numerical Solution of Hamilton-Jacobi-Bellman Equations by an Upwind Finite Volume Method

In this paper we present a finite volume method for solving Hamilton-Jacobi-Bellman(HJB) equations governing a class of optimal feedback control problems. This method is based on a finite volume discretization in state space coupled with an upwind finite difference technique, and on an implicit backward Euler finite differencing in time, which is absolutely stable. It is shown that the system matrix of the resulting discrete equation is an M-matrix. To show the effectiveness of this approach, numerical experiments on test problems with up to three states and two control variables were performed. The numerical results show that the method yields accurate approximate solutions to both the control and the state variables.     

Weak Galerkin finite element method for a class of time fractional generalized Burgers' equation

In this article, we use the weak Galerkin (WG) finite element method to study a class of time fractional generalized Burgers' equation. The existence of numerical solutions and the stability of fully discrete scheme are proved. Meanwhile, by applying the energy method, an optimal order error estimate in discrete L² norm is established. Numerical experiments are presented to validate the theoretical analysis.