Analysis of two‐grid discretization scheme for semilinear hyperbolic equations by mixed finite element methods

In this paper, the full discrete scheme of mixed finite element approximation is introduced for semilinear hyperbolic equations. To solve the nonlinear problem efficiently, two two‐grid algorithms are developed and analyzed. In this approach, the nonlinear system is solved on a coarse mesh with width H, and the linear system is solved on a fine mesh with width h≪H. Error estimates and convergence results of two‐grid method are derived in detail. It is shown that if we choose in the first algorithm and in the second algorithm, the two‐grid algorithms can achieve the same accuracy of the mixed finite element solutions. Finally, the numerical examples also show that the two‐grid method is much more efficient than solving the nonlinear mixed finite element system directly.     

A two‐grid method for mixed finite‐element solution of reaction‐diffusion equations

We present a scheme for solving two‐dimensional, nonlinear reaction‐diffusion equations, using a mixed finite‐element method. To linearize the mixed‐method equations, we use a two grid scheme that relegates all the Newton‐like iterations to a grid Δ_H much coarser than the original one Δ_h, with no loss in order of accuracy so long as the mesh sizes obey . The use of a multigrid‐based solver for the indefinite linear systems that arise at each coarse‐grid iteration, as well as for the similar system that arises on the fine grid, allows for even greater efficiency. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 317–332, 1999     

Two grid finite element discretization method for semi‐linear hyperbolic integro‐differential equations

In this paper, we will investigate a two grid finite element discretization method for the semi‐linear hyperbolic integro‐differential equations by piecewise continuous finite element method. In order to deal with the semi‐linearity of the model, we use the two grid technique and derive that once the coarse and fine mesh sizes H, h satisfy the relation h = H² for the two‐step two grid discretization method, the two grid method achieves the same convergence accuracy as the ordinary finite element method. Both theoretical analysis and numerical experiments are given to verify the results.     

A locally conservative Eulerian‐Lagrangian control‐volume method for transient advection‐diffusion equations

Characteristic methods generally generate accurate numerical solutions and greatly reduce grid orientation effects for transient advection‐diffusion equations. Nevertheless, they raise additional numerical difficulties. For instance, the accuracy of the numerical solutions and the property of local mass balance of these methods depend heavily on the accuracy of characteristics tracking and the evaluation of integrals of piecewise polynomials on some deformed elements generally with curved boundaries, which turns out to be numerically difficult to handle. In this article we adopt an alternative approach to develop an Eulerian‐Lagrangian control‐volume method (ELCVM) for transient advection‐diffusion equations. The ELCVM is locally conservative and maintains the accuracy of characteristic methods even if a very simple tracking is used, while retaining the advantages of characteristic methods in general. Numerical experiments show that the ELCVM is favorably comparable with well‐regarded Eulerian‐Lagrangian methods, which were previously shown to be very competitive with many well‐perceived methods. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A new stabilized finite element method for advection‐diffusion‐reaction equations

In this article, a new stabilized finite element method is proposed and analyzed for advection‐diffusion‐reaction equations. The key feature is that both the mesh‐dependent Péclet number and the mesh‐dependent Damköhler number are reasonably incorporated into the newly designed stabilization parameter. The error estimates are established, where, up to the regularity‐norm of the exact solution, the explicit‐dependence of the diffusivity, advection, reaction, and mesh size (or the dependence of the mesh‐dependent Péclet number and the mesh‐dependent Damköhler number) is revealed. Such dependence in the error bounds provides a mathematical justification on the effectiveness of the proposed method for any values of diffusivity, advection, dissipative reaction, and mesh size. Numerical results are presented to illustrate the performance of the method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 616–645, 2016     

A two‐grid method for characteristic expanded mixed finite element solution of miscible displacement problem

A combined method consisting of mixed finite element method (MFEM) for the pressure equation and expanded mixed finite element method with characteristics(CEMFEM) for the concentration equation is presented to solve the coupled system of incompressible miscible displacement problem. To solve the resulting nonlinear system of equations efficiently, the two‐grid algorithm relegates all of the Newton‐like iterations to grids much coarser than the original one, with no loss in order of accuracy. It is shown that coarse space can be extremely coarse and our algorithm achieve asymptotically optimal approximation when the mesh sizes satisfy $H=O(h^{\frac{1}{4}})$. Numerical experiment is provided to confirm our theoretical results.     

A uniform estimate for the MMOC for two‐dimensional advection‐diffusion equations

We prove an optimal‐order error estimate in a weighted energy norm for the modified method of characteristics (MMOC) and the modified method of characteristics with adjusted advection (MMOCAA) for two‐dimensional time‐dependent advection‐diffusion equations, in the sense that the generic constants in the estimates depend on certain Sobolev norms of the true solution but not on the scaling diffusion parameter ε. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

An upwind finite volume scheme and its maximum‐principle‐preserving ADI splitting for unsteady‐state advection‐diffusion equations

We develop an upwind finite volume (UFV) scheme for unsteady‐state advection‐diffusion partial differential equations (PDEs) in multiple space dimensions. We apply an alternating direction implicit (ADI) splitting technique to accelerate the solution process of the numerical scheme. We investigate and analyze the reason why the conventional ADI splitting does not satisfy maximum principle in the context of advection‐diffusion PDEs. Based on the analysis, we propose a new ADI splitting of the upwind finite volume scheme, the alternating‐direction implicit, upwind finite volume (ADFV) scheme. We prove that both UFV and ADFV schemes satisfy maximum principle and are unconditionally stable. We also derive their error estimates. Numerical results are presented to observe the performance of these schemes. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 211–226, 2003     

A posteriori error estimate for finite volume element method of the parabolic equations

In this article, residual‐type a posteriori error estimates are studied for finite volume element (FVE) method of parabolic equations. Residual‐type a posteriori error estimator is constructed and the reliable and efficient bounds for the error estimator are established. Residual‐type a posteriori error estimator can be used to assess the accuracy of the FVE solutions in practical applications. Some numerical examples are provided to confirm the theoretical results. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 259–275, 2017     

A two‐grid stabilized mixed finite element method for Darcy‐Forchheimer model

A two‐grid stabilized mixed finite element method based on pressure projection stabilization is proposed for the two‐dimensional Darcy‐Forchheimer model. We use the derivative of a smooth function, , to approximate the derivative of in constructing the two‐grid algorithm. The two‐grid method consists of solving a small nonlinear system on the coarse mesh and then solving a linear system on the fine mesh. There are a substantial reduction in computational cost. We prove the existence and uniqueness of solution of the discrete schemes on the coarse grid and the fine grid and obtain error estimates for the two‐grid algorithm. Finally, some numerical experiments are carried out to verify the accuracy and efficiency of the method.     

Taylor‐Galerkin B‐spline finite element method for the one‐dimensional advection‐diffusion equation

The advection‐diffusion equation has a long history as a benchmark for numerical methods. Taylor‐Galerkin methods are used together with the type of splines known as B‐splines to construct the approximation functions over the finite elements for the solution of time‐dependent advection‐diffusion problems. If advection dominates over diffusion, the numerical solution is difficult especially if boundary layers are to be resolved. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show the behavior of the method with emphasis on treatment of boundary conditions. Taylor‐Galerkin methods have been constructed by using both linear and quadratic B‐spline shape functions. Results shown by the method are found to be in good agreement with the exact solution. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

An Eulerian‐Lagrangian control volume scheme for two‐dimensional unsteady advection‐diffusion problems

We develop a mass conservative Eulerian‐Lagrangian control volume scheme (ELCVS) for the solution of the transient advection‐diffusion equations in two space dimensions. This method uses finite volume test functions over the space‐time domain defined by the characteristics within the framework of the class of Eulerian‐Lagrangian localized adjoint characteristic methods (ELLAM). It, therefore, maintains the advantages of characteristic methods in general, and of this class in particular, which include global mass conservation as well as a natural treatment of all types of boundary conditions. However, it differs from other methods in that class in the treatment of the mass storage integrals at the previous time step defined on deformed Lagrangian regions. This treatment is especially attractive for orthogonal rectangular Eulerian grids composed of block elements. In the algorithm, each deformed region is approximated by an eight‐node region with sides drawn parallel to the Eulerian grid, which significantly simplifies the integration compared with the approach used in finite volume ELLAM methods, based on backward tracking, while retaining local mass conservation at no additional expenses in terms of accuracy or CPU consumption. This is verified by numerical tests which show that ELCVS performs as well as standard finite volume ELLAM methods, which have previously been shown to outperform many other well‐received classes of numerical methods for the equations considered. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

A new combined finite element‐upwind finite volume method for convection‐dominated diffusion problems

In this article, we develop a combined finite element‐weighted upwind finite volume method for convection‐dominated diffusion problems in two dimensions, which discretizes the diffusion term with the standard finite element scheme, and the convection and source terms with the weighted upwind finite volume scheme. The developed method leads to a totally new scheme for convection‐dominated problems, which overcomes numerical oscillation, avoids numerical dispersion, and has high‐order accuracy. Stability analyses of the scheme are given for the problems with constant coefficients. Numerical experiments are presented to illustrate the stability and optimal convergence of our proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 799–818, 2016     

A uniformly optimal‐order estimate for finite volume method for advection‐diffusion equations

We prove an optimal‐order error estimate in a weighted energy norm for finite volume method for two‐dimensional time‐dependent advection–diffusion equations on a uniform space‐time partition of the domain. The generic constants in the estimates depend only on certain norms of the true solution but not on the scaling parameter. These estimates, combined with a priori stability estimates of the governing partial differential equations with full regularity, yield a uniform estimate of the finite volume method, in which the generic constants depend only on the Sobolev norms of the initial and right side data but not on the scaling parameter. We use the interpolation of spaces and stability estimates to derive a uniform estimate for problems with minimal or intermediate regularity, where the convergence rates are proportional to certain Besov norms of the initial and right‐hand side data. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 17‐43, 2014     

The characteristic finite element alternating direction method with moving meshes for nonlinear convection‐dominated diffusion problems

In modern numerical simulation of prospecting and exploiting oil‐gas resources and in environmental science, it is necessary to consider numerical method of nonlinear convection‐dominated diffusion problems. This thesis, starting from actual conditions such as the three‐dimensional characteristics of large‐scale science‐engineering computation, puts forward a kind of characteristic finite element alternating direction method with moving meshes. Some techniques, such as calculus of variations, operator‐splitting, generalized L² projection, energy method, negative norm estimate, the theory of prior estimates and techniques, are adopted. Optimal order estimates in L² norm are derived to determine the errors in the approximate solution. Thus the important theoretical problem has been solved. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

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     

A collocation method of lines for two‐sided space‐fractional advection‐diffusion equations with variable coefficients

We present the method of lines (MOL), which is based on the spectral collocation method, to solve space‐fractional advection‐diffusion equations (SFADEs) on a finite domain with variable coefficients. We focus on the cases in which the SFADEs consist of both left‐ and right‐sided fractional derivatives. To do so, we begin by introducing a new set of basis functions with some interesting features. The MOL, together with the spectral collocation method based on the new basis functions, are successfully applied to the SFADEs. Finally, four numerical examples, including benchmark problems and a problem with discontinuous advection and diffusion coefficients, are provided to illustrate the efficiency and exponentially accuracy of the proposed method.     

Implicit–explicit multistep characteristic finite element methods for nonlinear convection–diffusion equations

Implicit–explicit multistep characteristic methods are given for convection‐dominated diffusion equations. Multistep difference along characteristics of the one‐order hyperbolic part of the equation is used for discretization in time, and finite element method is used to discrete the space variables. The resulting schemes are consistent, stable and very efficient. Optimal‐rate of convergence is proved. Also, a note is given for a paper published earlier© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

A two‐grid stabilization method for solving the steady‐state Navier‐Stokes equations

We formulate a subgrid eddy viscosity method for solving the steady‐state incompressible flow problem. The eddy viscosity does not act on the large flow structures. Optimal error estimates are obtained for velocity and pressure. The numerical illustrations agree completely with the theoretical results. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Error estimates for two‐level penalty finite volume method for the stationary Navier–Stokes equations

Two‐level penalty finite volume method for the stationary Navier–Stokes equations based on the P₁ ? P₀ element is considered in this paper. The method involves solving one small penalty Navier–Stokes problem on a coarse mesh with mesh size H = ?^{1 / 4}h^{1 / 2}, a large penalty Stokes problem on a fine mesh with mesh size h, where 0 < ? < 1 is a penalty parameter. The method we study provides an approximate solution with the convergence rate of same order as the penalty finite volume solution (u_?h,p_?h), which involves solving one large penalty Navier–Stokes problem on a fine mesh with the same mesh size h. However, our method can save a large amount of computational time. Copyright © 2013 John Wiley & Sons, Ltd.