A finite element modified method of characteristics for convective heat transport

We propose a finite element modified method of characteristics for numerical solution of convective heat transport. The flow equations are the incompressible Navier‐Stokes equations including density variation through the Boussinesq approximation. The solution procedure consists of combining an essentially non‐oscillatory modified method of characteristics for time discretization with finite element method for space discretization. These numerical techniques associate the geometrical flexibility of the finite elements with the ability offered by modified method of characteristics to solve convection‐dominated flows using time steps larger than its Eulerian counterparts. Numerical results are shown for natural convection in a squared cavity and heat transport in the strait of Gibraltar. Performance and accuracy of the method are compared to other published data. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A monotone scheme for Hamilton–Jacobi equations via the nonstandard finite difference method

A usual way of approximating Hamilton–Jacobi equations is to couple space finite element discretization with time finite difference discretization. This classical approach leads to a severe restriction on the time step size for the scheme to be monotone. In this paper, we couple the finite element method with the nonstandard finite difference method, which is based on Mickens' rule of nonlocal approximation. The scheme obtained in this way is unconditionally monotone. The convergence of the new method is discussed and numerical results that support the theory are provided. Copyright © 2009 John Wiley & Sons, Ltd.     

On compactness of the velocity field in the incompressible limit of the full Navier–Stokes–Fourier system on large domains

The incompressible limit for the full Navier–Stokes–Fourier system is studied on a family of domains containing balls of the radius growing with a speed that dominates the inverse of the Mach number. It is shown that the velocity field converges strongly to its limit locally in space, in particular, the effect of the sound waves is eliminated by means of the local decay estimates for the acoustic wave equation. Copyright © 2008 John Wiley & Sons, Ltd.     

A compact finite difference method for the solution of the generalized Burgers–Fisher equation

In this article, numerical solutions of the generalized Burgers–Fisher equation are obtained using a compact finite difference method with minimal computational effort. To verify this, a combination of a sixth‐order compact finite difference scheme in space and a low‐storage third‐order total variation diminishing Runge–Kutta scheme in time have been used. The computed results with the use of this technique have been compared with the exact solution to show the accuracy of it. The approximate solutions to the equation have been computed without transforming the equation and without using linearization. Comparisons indicate that there is a very good agreement between the numerical solutions and the exact solutions in terms of accuracy. The present method is seen to be a very good alternative to some existing techniques for realistic problems. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

A superconvergent nonconforming mixed finite element method for the Navier–Stokes equations

The superconvergence for a nonconforming mixed finite element approximation of the Navier–Stokes equations is analyzed in this article. The velocity field is approximated by the constrained nonconforming rotated Q₁ (CNRQ₁) element, and the pressure is approximated by the piecewise constant functions. Under some regularity assumptions, the superconvergence estimates for both the velocity in broken H¹‐norm and the pressure in L²‐norm are obtained. Some numerical examples are presented to demonstrate our theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 646–660, 2016     

Iterative methods in penalty finite element discretizations for the Steady Navier–Stokes equations

Three penalty finite element methods are designed to solve numerically the steady Navier–Stokes equations, where the Stokes, Newton, and Oseen iteration methods are used, respectively. Moreover, the stability analysis and error estimate for these nine algorithms are provided. Finally, the numerical tests confirm the theoretical results of the presented algorithms. Meanwhile, the numerical investigations are provided to show that the proposed methods are efficient for solving the steady Navier–Stokes equations with the different viscosity. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 74‐94, 2014     

Economical Iterative Algorithms for Solving Stationary Problems of Mathematical Physics

We construct and investigate additive iterative methods of complete approximation for solving stationary problems of mathematical physics. We prove the convergence of the proposed methods and obtain error estimates without the requirement of commutativity of the decomposition operators. We provide the results of a computational experiment for a three-dimensional boundary-value problem. We consider possible generalizations of algorithms for equations with mixed derivatives and Navier–Stokes equation systems.     

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.     

Finite volume method based on stabilized finite elements for the nonstationary Navier–Stokes problem

A finite volume method based on stabilized finite element for the two‐dimensional nonstationary Navier–Stokes equations is investigated in this work. As in stabilized finite element method, macroelement condition is introduced for constructing the local stabilized formulation of the nonstationary Navier–Stokes equations. Moreover, for P₁ ? P₀ element, the H¹ error estimate of optimal order for finite volume solution (u_h,p_h) is analyzed. And, a uniform H¹ error estimate of optimal order for finite volume solution (u_h, p_h) is also obtained if the uniqueness condition is satisfied. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

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.     

A three‐level linearized finite difference scheme for the camassa–holm equation

The Camassa–Holm (CH) system is a strong nonlinear third‐order evolution equation. So far, the numerical methods for solving this problem are only a few. This article deals with the finite difference solution to the CH equation. A three‐level linearized finite difference scheme is derived. The scheme is proved to be conservative, uniquely solvable, and conditionally second‐order convergent in both time and space in the discrete L_∞ norm. Several numerical examples are presented to demonstrate the accuracy and efficiency of the proposed method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 451–471, 2014     

A higher-order moment method of the lattice Boltzmann model for the conservation law equation

In this paper, we proposed a higher-order moment method in the lattice Boltzmann model for the conservation law equation. In contrast to the lattice Bhatnagar–Gross–Krook (BGK) model, the higher-order moment method has a wide flexibility to select equilibrium distribution function. This method is based on so-called a series of partial differential equations obtained by using multi-scale technique and Chapman–Enskog expansion. According to Hirt’s heuristic stability theory, the stability of the scheme can be controlled by modulating some special moments to design the third-order dispersion term and the fourth-order dissipation term. As results, the conservation law equation is recovered with higher-order truncation error. The numerical examples show the higher-order moment method can be used to raise the accuracy of the truncation error of the lattice Boltzmann scheme for the conservation law equation.     

A parallel finite element variational multiscale method based on fully overlapping domain decomposition for incompressible flows

Based on fully overlapping domain decomposition and a recent variational multiscale method, a parallel finite element variational multiscale method for convection dominated incompressible flows is proposed and analyzed. In this method, each processor computes a local finite element solution in its own subdomain using a global mesh that is locally refined around its own subdomain, where a stabilization term based on two local Gauss integrations is adopted to stabilize the numerical form of the Navier–Stokes equations. Using the technical tool of local a priori estimate for the finite element solution, error bounds of the discrete solution are estimated. Algorithmic parameter scalings are derived. Numerical tests are also given to verify the theoretical predictions and demonstrate the effectiveness of the method. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 856–875, 2015     

A finite difference method for fractional partial differential equation

An implicit unconditional stable difference scheme is presented for a kind of linear space–time fractional convection–diffusion equation. The equation is obtained from the classical integer order convection–diffusion equations with fractional order derivatives for both space and time. First-order consistency, unconditional stability, and first-order convergence of the method are proven using a novel shifted version of the classical Grünwald finite difference approximation for the fractional derivatives. A numerical example with known exact solution is also presented, and the behavior of the error is examined to verify the order of convergence.     

Error analysis of a fully discrete finite element variational multiscale method for time‐dependent incompressible Navier–Stokes equations

A finite element variational multiscale method based on two local Gauss integrations is applied to solve numerically the time‐dependent incompressible Navier–Stokes equations. A significant feature of the method is that the definition of the stabilization term is derived via two local Guass integrations at element level, making it more efficient than the usual projection‐based variational multiscale methods. It is computationally cheap and gives an accurate approximation to the quantities sought. Based on backward Euler and Crank–Nicolson schemes for temporal discretization, we derive error bounds of the fully discrete solution which are first and second order in time, respectively. Numerical tests are also given to verify the theoretical predictions and demonstrate the effectiveness of the proposed method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A stabilized finite element method based on local polynomial pressure projection for the stationary Navier–Stokes equations

This article considers a stabilized finite element approximation for the branch of nonsingular solutions of the stationary Navier–Stokes equations based on local polynomial pressure projection by using the lowest equal-order elements. The proposed stabilized method has a number of attractive computational properties. Firstly, it is free from stabilization parameters. Secondly, it only requires the simple and efficient calculation of Gauss integral residual terms. Thirdly, it can be implemented at the element level. The optimal error estimate is obtained by the standard finite element technique. Finally, comparison with other methods, through a series of numerical experiments, shows that this method has better stability and accuracy.     

A block‐centered finite difference method for Darcy–Forchheimer model with variable Forchheimer number

A block‐centered finite difference scheme is introduced to solve the nonlinear Darcy–Forchheimer equation with variable Forchheimer number, in which the velocity and pressure can be approximated simultaneously. For variable Forchheimer number the second‐order error estimates for both pressure and velocity are established on nonuniform rectangular grid. An iteration process is given to solve the nonlinear system. Numerical experiments using the scheme show the consistency of the convergence rates of the presented methods with the theoretical analysis. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1603–1622, 2015     

A two-level finite element method for the Navier–Stokes equations based on a new projection

We consider a fully discrete two-level approximation for the time-dependent Navier–Stokes equations in two dimension based on a time-dependent projection. By defining this new projection, the iteration between large and small eddy components can be reflected by its associated space splitting. Hence, we can get a weakly coupled system of large and small eddy components. This two-level method applies the finite element method in space and Crank–Nicolson scheme in time. Moreover,the analysis and some numerical examples are shown that the proposed two-level scheme can reach the same accuracy as the classical one-level Crank–Nicolson method with a very fine mesh size h by choosing a proper coarse mesh size H. However, the two-level method will involve much less work.     

A nonstandard finite difference method for solving a mathematical model of HIV-TB co-infection

We design and analyze an unconditionally convergent nonstandard finite-difference method to study transmission dynamics of a mathematical model of HIV-TB co-infection. The dynamics of this model are studied using the qualitative theory of dynamical systems. These qualitative features of the continuous model are preserved by the numerical method that we propose in this paper. This method also preserves positivity of the solution which is one of the essential requirements when modelling epidemic diseases. Furthermore, we show that the numerical method is unconditionally stable. Competitive numerical results confirming theoretical investigations are provided. Comparisons are also made with other conventional approaches that are routinely used to solve these types of problems.     

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