Analysis of a decoupled time‐stepping algorithm for reduced MHD system modeling magneto‐convection

In this article, we analyze the stability and error estimate of a decoupled algorithm for a magneto‐convection problem. Magneto‐convection is assumed to be modeled by a coupled system of reduced magneto‐hydrodynamic (RMHD) equations and convection‐diffusion equation. The proposed algorithm applies the second‐order backward difference formula in time and finite element in space. To obtain a noniterative decouple algorithm from the fully discrete nonlinear system, we use a second‐order extrapolation in time to the nonlinear terms such that their skew symmetry properties are preserved. We prove the stability of the algorithm and derive error estimates without assuming any stability conditions. The algorithm is unconditionally stable and requires the solution of one RMHD problem and one convection‐diffusion equation per time step. Numerical test is presented that illustrates the accuracy and efficiency of the algorithm.     

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.     

A new two grid variational multiscale method for steady‐state natural convection problem

A two‐grid variational multiscale method based on two local Gauss integrations for solving the stationary natural convection problem is presented in this article. A significant feature of the method is that we solve the natural convection problem on a coarse mesh using finite element variational multiscale method based on two local Gauss integrations firstly, and then find a fine grid solution by solving a linearized problem on a fine grid. In the computation, we introduce two local Gauss integrations as a stabilizing term to replace the projection operator without adding other variables. The stability estimates and convergence analysis of the new method are derived. Ample numerical experiments are performed to validate the theoretical predictions and demonstrate the efficiency of the new method. Copyright © 2016 John Wiley & Sons, Ltd.     

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     

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     

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 robust algebraic multilevel preconditioner for non‐symmetric M‐matrices

Stable finite difference approximations of convection‐diffusion equations lead to large sparse linear systems of equations whose coefficient matrix is an M‐matrix, which is highly non‐symmetric when the convection dominates. For an efficient iterative solution of such systems, it is proposed to consider in the non‐symmetric case an algebraic multilevel preconditioning method formerly proposed for pure diffusion problems, and for which theoretical results prove grid independent convergence in this context. These results are supplemented here by a Fourier analysis that applies to constant coefficient problems with periodic boundary conditions whenever using an ‘idealized’ version of the two‐level preconditioner. Within this setting, it is proved that any eigenvalue λ of the preconditioned system satisfies for some real constant c such that . This result holds independently of the grid size and uniformly with respect to the ratio between convection and diffusion. Extensive numerical experiments are conducted to assess the convergence of practical two‐ and multi‐level schemes. These experiments, which include problems with highly variable and rotating convective flow, indicate that the convergence is grid independent. It deteriorates moderately as the convection becomes increasingly dominating, but the convergence factor remains uniformly bounded. This conclusion is supported for both uniform and some non‐uniform (stretched) grids. Copyright © 2000 John Wiley & Sons, Ltd.     

Decoupled energy‐law preserving numerical schemes for the Cahn–Hilliard–Darcy system

We study two novel decoupled energy‐law preserving and mass‐conservative numerical schemes for solving the Cahn‐Hilliard‐Darcy system which models two‐phase flow in porous medium or in a Hele–Shaw cell. In the first scheme, the velocity in the Cahn–Hilliard equation is treated explicitly so that the Darcy equation is completely decoupled from the Cahn–Hilliard equation. In the second scheme, an intermediate velocity is used in the Cahn–Hilliard equation which allows for the decoupling. We show that the first scheme preserves a discrete energy law with a time‐step constraint, while the second scheme satisfies an energy law without any constraint and is unconditionally stable. Ample numerical experiments are performed to gauge the efficiency and robustness of our scheme. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 936–954, 2016     

Integrated fast and high‐accuracy computation of convection diffusion equations using multiscale multigrid method

We present an explicit sixth‐order compact finite difference scheme for fast high‐accuracy numerical solutions of the two‐dimensional convection diffusion equation with variable coefficients. The sixth‐order scheme is based on the well‐known fourth‐order compact (FOC) scheme, the Richardson extrapolation technique, and an operator interpolation scheme. For a particular implementation, we use multiscale multigrid method to compute the fourth‐order solutions on both the coarse grid and the fine grid. Then, an operator interpolation scheme combined with the Richardson extrapolation technique is used to compute a sixth‐order accurate fine grid solution. We compare the computed accuracy and the implementation cost of the new scheme with the standard nine‐point FOC scheme and Sun–Zhang's sixth‐order method. Two convection diffusion problems are solved numerically to validate our proposed sixth‐order scheme. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

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

Numerical treatment of two‐parameter singularly perturbed parabolic convection diffusion problems with non‐smooth data

In the present work, we consider a parabolic convection‐diffusion‐reaction problem where the diffusion and convection terms are multiplied by two small parameters, respectively. In addition, we assume that the convection coefficient and the source term of the partial differential equation have a jump discontinuity. The presence of perturbation parameters leads to the boundary and interior layers phenomena whose appropriate numerical approximation is the main goal of this paper. We have developed a uniform numerical method, which converges almost linearly in space and time on a piecewise uniform space adaptive Shishkin‐type mesh and uniform mesh in time. Error tables based on several examples show the convergence of the numerical solutions. In addition, several numerical simulations are presented to show the effectiveness of resolving layer behavior and their locations.     

Convergence analysis of a vertex‐centered finite volume scheme for a copper heap leaching model

In this paper a two‐dimensional solute transport model is considered to simulate the leaching of copper ore tailing using sulfuric acid as the leaching agent. The mathematical model consists in a system of differential equations: two diffusion–convection‐reaction equations with Neumann boundary conditions, and one ordinary differential equation. The numerical scheme consists in a combination of finite volume and finite element methods. A Godunov scheme is used for the convection term and an P1‐FEM for the diffusion term. The convergence analysis is based on standard compactness results in L². Some numerical examples illustrate the effectiveness of the scheme. Copyright © 2009 John Wiley & Sons, Ltd.     

A decoupling two‐grid algorithm for the mixed Stokes‐Darcy model with the Beavers‐Joseph interface condition

In this article, a decoupling scheme based on two‐grid finite element for the mixed Stokes‐Darcy problem with the Beavers‐Joseph interface condition is proposed and investigated. With a restriction of a physical parameter α, we derive the numerical stability and error estimates for the scheme. Numerical experiments indicate that such two‐grid based decoupling finite element schemes are feasible and efficient. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1066–1082, 2014     

Two‐grid methods for time‐harmonic Maxwell equations

In this paper, we develop several two‐grid methods for the Nédélec edge finite element approximation of the time‐harmonic Maxwell equations. We first present a two‐grid method that uses a coarse space to solve the original problem and then use a fine space to solve a corresponding symmetric positive definite problem. Then, we present two types of iterative two‐grid methods, one is to add the kernel of the curl ‐operator in the fine space to a coarse mesh space to solve the original problem and the other is to use an inner iterative method for dealing with the kernel of the curl ‐operator in the fine space and the coarse space, separately. We provide the error estimates for the first two methods and present numerical experiments to show the efficiency of our methods.Copyright © 2012 John Wiley & Sons, Ltd.     

A second‐order partitioned method with different subdomain time steps for the evolutionary Stokes‐Darcy system

A second‐order decoupled algorithm for the nonstationary Stokes‐Darcy system, which allows different time steps in two subregions, is proposed and analyzed in this paper. The algorithm, which is a combination of the second‐order backward differentiation formula and second‐order extrapolation method, uncouples the problem into two decoupled problems per time step. We prove the unconditional stability and long‐time stability of the decoupled scheme with different time steps and derive error estimates of this decoupled algorithm using finite element spatial discretization. Numerical experiments are provided to illustrate the accuracy, effectiveness, and stability of the decoupled algorithm and show its advantages of increasing accuracy and efficiency.     

One‐level and two‐level finite element methods for the Boussinesq problem

Previous works on the convergence of numerical methods for the Boussinesq problem were conducted, while the optimal L²‐norm error estimates for the velocity and temperature are still lacked. In this paper, the backward Euler scheme is used to discrete the time terms, standard Galerkin finite element method is adopted to approximate the variables. The MINI element is used to approximate the velocity and pressure, the temperature field is simulated by the linear polynomial. Under some restriction on the time step, we firstly present the optimal L² error estimates of approximate solutions. Secondly, two‐level method based on Stokes iteration for the Boussinesq problem is developed and the corresponding convergence results are presented. By this method, the original problem is decoupled into two small linear subproblems. Compared with the standard Galerkin method, the two‐level method not only keeps good accuracy but also saves a lot of computational cost. Finally, some numerical examples are provided to support the established theoretical analysis.     

The second order projection method in time for the time-dependent natural convection problem

We consider the second-order projection schemes for the time-dependent natural convection problem. By the projection method, the natural convection problem is decoupled into two linear subproblems, and each subproblem is solved more easily than the original one. The error analysis is accomplished by interpreting the second-order time discretization of a perturbed system which approximates the time-dependent natural convection problem, and the rigorous error analysis of the projection schemes is presented. Our main results of the second order projection schemes for the time-dependent natural convection problem are that the convergence for the velocity and temperature are strongly second order in time while that for the pressure is strongly first order in time.     

Local projection stabilized and characteristic decoupled scheme for the fluid–fluid interaction problems

In this article, we propose and analyse a local projection stabilized and characteristic decoupled scheme for the fluid–fluid interaction problems. We use the method of characteristics type to avert the difficulties caused by the nonlinear term, and use the local projection stabilized method to control spurious oscillations in the velocities due to dominant convection, and use a geometric averaging idea to decouple the monolithic problems. The stability analysis is derived and numerical tests are performed to demonstrate the robustness of this new method. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 704–723, 2017     

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.     

A monotone conservative Eulerian–Lagrangian scheme for reaction‐diffusion‐convection equations modeling chemotaxis

A numerical method for convection dominated diffusion problems, that exploits the use of characteristics, is derived and analyzed. It is shown to preserve positivity of solutions and be locally mass conserving. Stability, consistency and order one convergence are verified. Because of the way in which it determines characteristic pre‐images of grid cells, the method can be easily implemented for 1‐, 2‐, or 3‐dimensional problems on rectangular grids.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007