Numerical analysis for the mixed Navier–Stokes and Darcy Problem with the Beavers–Joseph interface condition

In this article, we consider the coupled Navier–Stokes and Darcy problem with the Beavers–Joseph interface condition. With suitable restrictions of physical parameters α and ν, we prove the existence and local uniqueness of a weak solution. Then we propose a coupled finite element scheme and a decoupled and linearized scheme based on two‐grid finite element. Under suitable further restrictions, their optimal error estimates are obtained. Finally numerical experiments indicate the validity of the theoretical results as well as the efficiency and effectiveness of the decoupled and linearized two‐grid algorithm. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1009–1030, 2015     

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.     

Second order fully discrete defect‐correction scheme for nonstationary conduction‐convection problem at high Reynolds number

This survey enfolds rigorous analysis of the defect‐correction finite element (FE) method for the time‐dependent conduction‐convection problem which based on the Crank‐Nicolson scheme. The method consists of two steps: solve a nonlinear problem with an added artificial viscosity term on a FE grid and correct the solutions on the same grid using a linearized defect‐correction technique. The stability and optimal error estimate of the fully discrete scheme are derived. As a consequence, the effectiveness of the method to deal with high Reynolds number is illustrated in several numerical experiments. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 681–703, 2017     

Two‐level Fourier analysis of multigrid for higher‐order finite‐element discretizations of the Laplacian

In this paper, we employ local Fourier analysis (LFA) to analyze the convergence properties of multigrid methods for higher‐order finite‐element approximations to the Laplacian problem. We find that the classical LFA smoothing factor, where the coarse‐grid correction is assumed to be an ideal operator that annihilates the low‐frequency error components and leaves the high‐frequency components unchanged, fails to accurately predict the observed multigrid performance and, consequently, cannot be a reliable analysis tool to give good performance estimates of the two‐grid convergence factor. While two‐grid LFA still offers a reliable prediction, it leads to more complex symbols that are cumbersome to use to optimize parameters of the relaxation scheme, as is often needed for complex problems. For the purposes of this analytical optimization as well as to have simple predictive analysis, we propose a modification that is “between” two‐grid LFA and smoothing analysis, which yields reasonable predictions to help choose correct damping parameters for relaxation. This exploration may help us better understand multigrid performance for higher‐order finite element discretizations, including for Q₂‐Q₁ (Taylor‐Hood) elements for the Stokes equations. Finally, we present two‐grid and multigrid experiments, where the corrected parameter choice is shown to yield significant improvements in the resulting two‐grid and multigrid convergence factors.     

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 multigrid compact finite difference method for solving the one‐dimensional nonlinear sine‐Gordon equation

The aim of this paper is to propose a multigrid method to obtain the numerical solution of the one‐dimensional nonlinear sine‐Gordon equation. The finite difference equations at all interior grid points form a large sparse linear system, which needs to be solved efficiently. The solution cost of this sparse linear system usually dominates the total cost of solving the discretized partial differential equation. The proposed method is based on applying a compact finite difference scheme of fourth‐order for discretizing the spatial derivative and the standard second‐order central finite difference method for the time derivative. The proposed method uses the Richardson extrapolation method in time variable. The obtained system has been solved by V‐cycle multigrid (VMG) method, where the VMG method is used for solving the large sparse linear systems. The numerical examples show the efficiency of this algorithm for solving the one‐dimensional sine‐Gordon equation. Copyright © 2014 John Wiley & Sons, Ltd.     

Two‐grid mixed finite‐element methods for nonlinear Schrödinger equations

Two‐grid mixed finite element schemes are developed for solving both steady state and unsteady state nonlinear Schrödinger equations. The schemes use discretizations based on a mixed finite‐element method. The two‐grid approach yields iterative procedures for solving the nonlinear discrete equations. The idea is to relegate all of the Newton‐like iterations to grids much coarser than the final one, with no loss in order of accuracy. Numerical tests are performed. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 63‐73, 2012     

Two‐grid method for two‐dimensional nonlinear Schrödinger equation by finite element method

A conservative two‐grid finite element scheme is presented for the two‐dimensional nonlinear Schrödinger equation. One Newton iteration is applied on the fine grid to linearize the fully discrete problem using the coarse‐grid solution as the initial guess. Moreover, error estimates are conducted for the two‐grid method. It is shown that the coarse space can be extremely coarse, with no loss in the order of accuracy, and still achieve the asymptotically optimal approximation as long as the mesh sizes satisfy in the two‐grid method. The numerical results show that this method is very effective.     

Two‐Grid method for nonlinear parabolic equations by expanded mixed finite element methods

In this article, we develop a two‐grid algorithm for nonlinear reaction diffusion equation (with nonlinear compressibility coefficient) discretized by expanded mixed finite element method. The key point is to use two‐grid scheme to linearize the nonlinear term in the equations. The main procedure of the algorithm is solving a small‐scaled nonlinear equations on the coarse grid and dealing with a linearized system on the fine space using the Newton iteration with the coarse grid solution. Error estimation to the expanded mixed finite element solution is analyzed in detail. We also show that two‐grid solution achieves the same accuracy as long as the mesh sizes satisfy H = O(h^1/2). Two numerical experiments are given to verify the effectiveness of the algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

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     

Superconvergence of a two‐grid method for fourth‐order dispersive‐dissipative wave equations

In this paper, the superconvergence analysis of a two‐grid method (TGM) with low‐order finite elements is presented for the fourth‐order dispersive‐dissipative wave equations for a second order fully discrete scheme. The superclose estimates in the H¹‐norm on the two grids are obtained by the combination technique of the interpolation and Ritz projection. Then, with the help of the interpolated postprocessing technique, the global superconvergence properties are deduced. Finally, numerical results are provided to show the performance of the proposed TGM for conforming bilinear element and nonconforming element, respectively. It shows that the TGM is an effective method to the problem considered of our paper compared with the traditional Galerkin finite element method (FEM).     

WEB‐spline based mesh‐free finite element analysis for the heat equation and the time‐dependent Navier–Stokes equation: A survey

In this article, we give some numerical techniques and error estimates using web‐spline based mesh‐free finite element method for the heat equation and the time‐dependent Navier–Stokes equations on bounded domains. The web‐spline method uses weighted extended B‐splines on a regular grid as basis functions and does not require any grid generation. We demonstrate the method by providing numerical results for the Poisson's and stationary Stokes equation. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

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 Two‐Parameter Stabilized Finite Element Method for Incompressible Flows

Based on two‐grid discretizations, a two‐parameter stabilized finite element method for the steady incompressible Navier–Stokes equations at high Reynolds numbers is presented and studied. In this method, a stabilized Navier–Stokes problem is first solved on a coarse grid, and then a correction is calculated on a fine grid by solving a stabilized linear problem. The stabilization term for the nonlinear Navier–Stokes equations on the coarse grid is based on an elliptic projection, which projects higher‐order finite element interpolants of the velocity into a lower‐order finite element interpolation space. For the linear problem on the fine grid, either the same stabilization approach (with a different stabilization parameter) as that for the coarse grid problem or a completely different stabilization approach could be employed. Error bounds for the discrete solutions are estimated. Algorithmic parameter scalings of the method are also derived. The theoretical results show that, with suitable scalings of the algorithmic parameters, this method can yield an optimal convergence rate. Numerical results are provided to verify the theoretical predictions and demonstrate the effectiveness of the proposed method. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 425–444, 2017     

High‐order approximation of 2D convection‐diffusion equation on hexagonal grids

We derive a fourth‐order finite difference scheme for the two‐dimensional convection‐diffusion equation on an hexagonal grid. The difference scheme is defined on a single regular hexagon of size h over a seven‐point stencil. Numerical experiments are conducted to verify the high accuracy of the derived scheme, and to compare it with the standard second‐order central difference scheme. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

A non‐standard finite difference scheme for an advection‐diffusion‐reaction equation

In this note, a non‐standard finite difference (NSFD) scheme is proposed for an advection‐diffusion‐reaction equation with nonlinear reaction term. We first study the diffusion‐free case of this equation, that is, an advection‐reaction equation. Two exact finite difference schemes are constructed for the advection‐reaction equation by the method of characteristics. As these exact schemes are complicated and are not convenient to use, an NSFD scheme is derived from the exact scheme. Then, the NSFD scheme for the advection‐reaction equation is combined with a finite difference space‐approximation of the diffusion term to provide a NSFD scheme for the advection‐diffusion‐reaction equation. This new scheme could preserve the fixed points, the positivity, and the boundedness of the solution of the original equation. Numerical experiments verify the validity of our analytical results. Copyright © 2014 JohnWiley & Sons, Ltd.     

Superconvergence of the split least‐squares method for second‐order hyperbolic equations

This article studies superconvergence phenomena of the split least‐squares mixed finite element method for second‐order hyperbolic equations. By selecting the least‐squares functional properly, the procedure can be split into two independent symmetric positive definite subprocedures, one of which is for the primitive unknown and the other is for the flux. Based on interpolation operators and an auxiliary projection, superconvergent H¹ error estimates for the primary variable u and L² error estimates for the introduced flux variable σ are obtained under the standard quasiuniform assumptions on finite element partition. A numerical example is given to show the performance of the introduced scheme. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 222‐238, 2014     

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 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.