On the stability of alternating‐direction explicit methods for advection‐diffusion equations

Alternating‐Direction Explicit (A.D.E.) finite‐difference methods make use of two approximations that are implemented for computations proceeding in alternating directions, e.g., from left to right and from right to left, with each approximation being explicit in its respective direction of computation. Stable A.D.E. schemes for solving the linear parabolic partial differential equations that model heat diffusion are well‐known, as are stable A.D.E. schemes for solving the first‐order equations of fluid advection. Several of these are combined here to derive A.D.E. schemes for solving time‐dependent advection‐diffusion equations, and their stability characteristics are discussed. In each case, it is found that it is the advection term that limits the stability of the scheme. The most stable of the combinations presented comprises an unconditionally stable approximation for computations carried out in the direction of advection of the system, from left to right in this case, and a conditionally stable approximation for computations proceeding in the opposite direction. To illustrate the application of the methods and verify the stability conditions, they are applied to some quasi‐linear one‐dimensional advection‐diffusion problems. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

A modified alternating‐direction finite volume method for modeling secondary hydrocarbon migration and accumulation processes

We present a finite volume method for the coupled systems of nonlinear partial differential equations arising in secondary hydrocarbon migration and accumulation processes in basin modeling. We propose a modified alternating direction implicit splitting acceleration technique to solve the scheme in an effective and robust manner. We apply the scheme to simulate model problems with clear physical features and to secondary migration and accumulation problems in real petroleum reservoirs. The numerical experiments show that the scheme is very robust and generates stable and physically reasonable numerical results, even in the extremely long time period that ranges from millions of years to tens of millions of years. This clearly shows the strong potential of the method developed. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 254–270, 2003     

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 quadrature finite element method for semilinear second‐order hyperbolic problems

In this work we propose and analyze a fully discrete modified Crank–Nicolson finite element (CNFE) method with quadrature for solving semilinear second‐order hyperbolic initial‐boundary value problems. We prove optimal‐order convergence in both time and space for the quadrature‐modified CNFE scheme that does not require nonlinear algebraic solvers. Finally, we demonstrate numerically the order of convergence of our scheme for some test problems. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

一阶线性双曲方程的耗散谱元法

讨论了二维一阶线性变系数双曲方程的耗散谱元法,得到拟最优估计.数值结果表明,耗散谱元法对于具有较复杂边界条件的问题同样有效,对于有限光滑问题,耗散谱元法能够得到比传统的谱元法更好的结果.     

A space‐time mixed‐hybrid finite element method for the damped wave equation

A space‐time finite element method is introduced to solve the linear damped wave equation. The scheme is constructed in the framework of the mixed‐hybrid finite element methods, and where an original conforming approximation of H(div;Ω) is used, the latter permits us to obtain an upwind scheme in time. We establish the link between the nonstandard finite difference scheme recently introduced by Mickens and Jordan and the scheme proposed. In this regard, two approaches are considered and in particular we employ a formulation allowing the solution to be marched in time, i.e., one only needs to consider one time increment at a time. Numerical results are presented and compared with the analytical solution illustrating good performance of the present method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

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     

Legendre multiwavelet Galerkin method for solving the hyperbolic telegraph equation

Recently, it is found that telegraph equation is more suitable than ordinary diffusion equation in modeling reaction diffusion for such branches of sciences. In this article a numerical method for solving the one‐dimensional hyperbolic telegraph equation is presented. The method is based upon Legendre multiwavelet approximations. The properties of Legendre multiwavelet are first presented. These properties together with Galerkin method are then utilized to reduce the telegraph equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Developing weak Galerkin finite element methods for the wave equation

In this article, we extend the recently developed weak Galerkin method to solve the second‐order hyperbolic wave equation. Many nice features of the weak Galerkin method have been demonstrated for elliptic, parabolic, and a few other model problems. This is the initial exploration of the weak Galerkin method for solving the wave equation. Here we successfully developed and established the stability and convergence analysis for the weak Galerkin method for solving the wave equation. Numerical experiments further support the theoretical analysis. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 868–884, 2017     

A method based on meshless approach for the numerical solution of the two‐space dimensional hyperbolic telegraph equation

In the current article, we investigate the RBF solution of second‐order two‐space dimensional linear hyperbolic telegraph equation. For this purpose, we use a combination of boundary knot method (BKM) and analog equation method (AEM). The BKM is a meshfree, boundary‐only and integration‐free technique. The BKM is an alternative to the method of fundamental solution to avoid the fictitious boundary and to deal with low accuracy, singular integration and mesh generation. Also, on the basis of the AEM, the governing operator is substituted by an equivalent nonhomogeneous linear one with known fundamental solution under the same boundary conditions. Finally, several numerical results and discussions are demonstrated to show the accuracy and efficiency of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Alternating direction finite volume element methods for 2D parabolic partial differential equations

On the basis of rectangular partition and bilinear interpolation, this article presents alternating direction finite volume element methods for two dimensional parabolic partial differential equations and gives three computational schemes, one is analogous to Douglas finite difference scheme with second order splitting error, the second has third order splitting error, and the third is an extended locally one dimensional scheme. Optimal L² norm or H¹ semi‐norm error estimates are obtained for these schemes. Finally, two numerical examples illustrate the effectiveness of the schemes. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Singly diagonally implicit runge‐kutta method for time‐dependent reaction‐diffusion equation

In this article, an efficient fourth‐order accurate numerical method based on Padé approximation in space and singly diagonally implicit Runge‐Kutta method in time is proposed to solve the time‐dependent one‐dimensional reaction‐diffusion equation. In this scheme, we first approximate the spatial derivative using the second‐order central finite difference then improve it to fourth‐order by applying Padé approximation. A three stage fourth‐order singly diagonally implicit Runge‐Kutta method is then used to solve the resulting system of ordinary differential equations. It is also shown that the scheme is unconditionally stable, and is suitable for stiff problems. Several numerical examples are solved by the scheme and the efficiency and accuracy of the new scheme are compared with two widely used high‐order compact finite difference methods. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1423–1441, 2011     

A method for finding coefficients of a quasilinear hyperbolic equation

The inverse problem of finding the coefficients q(s) and p(s) in the equation u _tt = a ² u _xx + q(u)u _t ? p(u)u _x is investigated. As overdetermination required in the inverse setting, two additional conditions are set: a boundary condition and a condition with a fixed value of the timelike variable. An iteration method for solving the inverse problem is proposed based on an equivalent system of integral equations of the second kind. A uniqueness theorem and an existence theorem in a small domain are proved for the inverse problem to substantiate the convergence of the algorithm.     

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.     

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.     

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     

A splitting positive definite mixed element method for second‐order hyperbolic equations

In this article, we establish a new mixed finite element procedure, in which the mixed element system is symmetric positive definite, to solve the second‐order hyperbolic equations. The convergence of the mixed element methods with continuous‐ and discrete‐time scheme is proved. And the corresponding error estimates are given. Finally some numerical results are presented. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009