Alternating direction implicit OSC scheme for the two-dimensional fractional evolution equation with a weakly singular kernel

In this paper, a new kind of alternating direction implicit (ADI) Crank-Nicolson-type orthogonal spline collocation (OSC) method is formulated for the two-dimensional fractional evolution equation with a weakly singular kernel arising in the theory of linear viscoelasticity. The novel OSC method is used for the spatial discretization, and ADI Crank-Nicolson-type method combined with the second order fractional quadrature rule are considered for the temporal component. The stability of proposed scheme is rigourously established, and nearly optimal order error estimate is also derived. Numerical experiments are conducted to support the predicted convergence rates and also exhibit expected super-convergence phenomena.     

Galerkin-Legendre spectral schemes for nonlinear space fractional Schrödinger equation

In the paper, we first propose a Crank-Nicolson Galerkin-Legendre (CN-GL) spectral scheme for the one-dimensional nonlinear space fractional Schrödinger equation. Convergence with spectral accuracy is proved for the spectral approximation. Further, a Crank-Nicolson ADI Galerkin-Legendre spectral method for the two-dimensional nonlinear space fractional Schrödinger equation is developed. The proposed schemes are shown to be efficient with second-order accuracy in time and spectral accuracy in space which are higher than some recently studied methods. Moreover, some numerical results are demonstrated to justify the theoretical analysis.     

Alternating direction implicit orthogonal spline collocation on some non-rectangular regions with inconsistent partitions

The alternating direction implicit (ADI) method is a highly efficient technique for solving multi-dimensional dependent initial-boundary value problems on rectangles. Earlier we have used the ADI technique in conjunction with orthogonal spline collocation (OSC) for discretization in space to solve parabolic problems on rectangles and rectangular polygons. Recently, we extended applications of ADI OSC schemes to the solution of parabolic problems on some non-rectangular regions that allow for consistent nonuniform partitions. However, for many regions, it is impossible to construct such partitions. Therefore, in this paper, we show how to extend our approach further to solve parabolic problems on some non-rectangular regions using inconsistent uniform partitions. Numerical results are presented using piecewise Hermite cubic polynomials for spatial discretizations and our ADI OSC scheme for parabolic problems to demonstrate its performance on several regions.     

A backward Euler alternating direction implicit difference scheme for the three‐dimensional fractional evolution equation

A backward Euler alternating direction implicit (ADI) difference scheme is formulated and analyzed for the three‐dimensional fractional evolution equation. In our method, the Riemann‐Liouville fractional integral term is treated by means of first order convolution quadrature suggested by Lubich. Meanwhile, an ADI technique is adopted to reduce the multidimensional problem to a series of one‐dimensional problems. A fully discrete difference scheme is constructed with space discretization by finite difference method. Two new inner products and corresponding norms are defined to analyze the scheme. The verification of stability and convergence is based on the nonnegative character of the real quadratic form associated with the convolution quadrature. Numerical experiments are reported to demonstrate the efficiency of our scheme.     

Fourth-Order Splitting Methods for Time-Dependant Differential Equations

This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we study both parabolic and hyperbolic equations.We focus on ADI (alternating direction implicit) methods and LOD (locally one-dimensional) methods,which are standard splitting methods of lower order,e.g.second-order.Our aim is to develop higher-order ADI methods,which are performed by Richardson extrapolation,Crank-Nicolson methods and higher-order LOD methods,based on locally higher-order methods.We discuss the new theoretical results of the stability and consistency of the ADI methods.The main idea is to apply a higher- order time discretization and combine it with the ADI methods.We also discuss the dis- cretization and splitting methods for first-order and second-order evolution equations. The stability analysis is given for the ADI method for first-order time derivatives and for the LOD (locally one-dimensional) methods for second-order time derivatives.The higher-order methods are unconditionally stable.Some numerical experiments verify our results.     

Chebyshev-Legendre pseudo-spectral method for the generalised Burgers-Fisher equation

In this paper, we consider numerical approximation of generalised Burgers-Fisher equation using the pseudo-spectral method. For the time discretization we apply Crank-Nicolson /leapfrog scheme. The space discretization is based on Legendre Galerkin formulation while the Chebyshev-Gauss-Lobatto (CGL) nodes are used in practical computation, which is called “Chebyshev-Legendre” method. The stability and convergence are rigorously set up. Numerical experiments are presented to demonstrate the effectiveness of the method and to confirm the theoretical results.     

An a posteriori error estimate for a dual mixed method applied to Stokes system with non-null source terms

Advances in Computational Mathematics - We propose and analyze a time-stepping Crank-Nicolson(CN) alternating direction implicit(ADI) scheme combined with an arbitrary-order orthogonal spline...     

Gauss‐Lobatto‐Legendre‐Birkhoff pseudospectral approximations for the multi‐term time fractional diffusion‐wave equation with Neumann boundaryconditions

A second‐order finite difference/pseudospectral scheme is proposed for numerical approximation of multi‐term time fractional diffusion‐wave equation with Neumann boundary conditions. The scheme is based upon the weighted and shifted Grünwald difference operators approximation of the time fractional calculus and Gauss‐Lobatto‐Legendre‐Birkhoff (GLLB) pseudospectral method for spatial discretization. The unconditionally stability and convergence of the scheme are rigorously proved. Numerical examples are carried out to verify theoretical results.     

Fast evaluation and high accuracy finite element approximation for the time fractional subdiffusion equation

In this article, an efficient algorithm for the evaluation of the Caputo fractional derivative and the superconvergence property of fully discrete finite element approximation for the time fractional subdiffusion equation are considered. First, the space semidiscrete finite element approximation scheme for the constant coefficient problem is derived and supercloseness result is proved. The time discretization is based on the L1‐type formula, whereas the space discretization is done using, the fully discrete scheme is developed. Under some regularity assumptions, the superconvergence estimate is proposed and analyzed. Then, extension to the case of variable coefficients is also discussed. To reduce the computational cost, the fast evaluation scheme of the Caputo fractional derivative to solve the fractional diffusion equations is designed. Finally, numerical experiments are presented to support the theoretical results.     

Two-time level ADI finite volume method for a class of second-order hyperbolic problems

In this paper, we develop a two-time level alternating direction implicit (ADI) method for a class of second-order hyperbolic problems on a rectangular domain. The method builds on the finite volume method with biquadratic basis functions for the discretization in space, and a Crank-Nicolson approach for the time stepping. We obtain a second-order error estimation in the H¹ norm. Numerical experiments are performed to demonstrate the theoretical findings.     

CONVERGENCE OF THE CRANK-NICOLSON/NEWTON SCHEME FOR NONLINEAR PARABOLIC PROBLEM

In this paper, the Crank-Nicolson/Newton scheme for solving numerically secondorder nonlinear parabolic problem is proposed. The standard Galerkin finite element method based on P2 conforming elements is used to the spatial discretization of the problem and the Crank-Nicolson/Newton scheme is applied to the time discretization of the resulted finite element equations. Moreover, assuming the appropriate regularity of the exact solution and the finite element solution, we obtain optimal error estimates of the fully discrete CrankNicolson/Newton scheme of nonlinear parabolic problem. Finally, numerical experiments are presented to show the efficient performance of the proposed scheme.     

非定常 N-S 方程的非线性 Galerkin 算法及其误差估计

本文给出了二维非定常N-S方程的三种数值格式,其中空间变量用谱非线性Galerkin算法进行离散,时间变量用有限差分离散,并研究了这些格式数值解的逼近精度．最后,给出了部分数值计算结果．     

A fourth-order compact ADI method for solving two-dimensional unsteady convection–diffusion problems

In this article, an exponential high-order compact (EHOC) alternating direction implicit (ADI) method, in which the Crank–Nicolson scheme is used for the time discretization and an exponential fourth-order compact difference formula for the steady-state 1D convection–diffusion problem is used for the spatial discretization, is presented for the solution of the unsteady 2D convection–diffusion problems. The method is temporally second-order accurate and spatially fourth order accurate, which requires only a regular five-point 2D stencil similar to that in the standard second-order methods. The resulting EHOC ADI scheme in each ADI solution step corresponds to a strictly diagonally dominant tridiagonal matrix equation which can be inverted by simple tridiagonal Gaussian decomposition and may also be solved by application of the one-dimensional tridiagonal Thomas algorithm with a considerable saving in computing time. The unconditionally stable character of the method was verified by means of the discrete Fourier (or von Neumann) analysis. Numerical examples are given to demonstrate the performance of the method proposed and to compare mostly it with the high order ADI method of Karaa and Zhang and the spatial third-order compact scheme of Note and Tan.     

A Stabilized Galerkin Scheme for the Convection-Diffusion-Reaction Equations

A fully discrete stabilized scheme is proposed for solving the time-dependent convection-diffusion-reaction equations. A time derivative term results in our stabilized algorithm. The finite element method for spatial discretization and the backward Euler or Crank-Nicolson scheme for time discretization are employed. The long-time stability and convergence are established in this article. Finally, some numerical experiments are provided to confirm the theoretical analysis.     

Weak convergence of finite element approximations of linear stochastic evolution equations with additive noise II. Fully discrete schemes

We present an abstract framework for analyzing the weak error of fully discrete approximation schemes for linear evolution equations driven by additive Gaussian noise. First, an abstract representation formula is derived for sufficiently smooth test functions. The formula is then applied to the wave equation, where the spatial approximation is done via the standard continuous finite element method and the time discretization via an I-stable rational approximation to the exponential function. It is found that the rate of weak convergence is twice that of strong convergence. Furthermore, in contrast to the parabolic case, higher order schemes in time, such as the Crank-Nicolson scheme, are worthwhile to use if the solution is not very regular. Finally we apply the theory to parabolic equations and detail a weak error estimate for the linearized Cahn-Hilliard-Cook equation as well as comment on the stochastic heat equation.     

Orthogonal spline collocation methods for the stream function‐vorticity formulation of the Navier–Stokes equations

An orthogonal spline collocation (OSC) spatial discretization is proposed for the solution of the fully coupled stream function‐vorticity formulation of the Navier–Stokes equations in two dimensions. For the time‐stepping, a three‐level leapfrog scheme is employed. This method is algebraically linear, and, at each time step, gives rise to a system of linear equations of the form arising in the OSC approximation of the biharmonic Dirichlet problem and can be solved by a fast direct method. Error estimates in the H^l–norm in space, l = 1,2, are derived for the semi‐discrete method and the fully‐discrete leapfrog scheme which is also shown to be second order accurate in time. Numerical results are presented which confirm the theoretical analysis and exhibit superconvergence phenomena, which provide superconvergent approximations to the components of the velocity. © John Wiley & Sons, Inc. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

带跳随机波动率模型的高阶ADI分裂格式

针对带跳随机波动率模型满足的偏积分微分方程,提出一种新的高阶交替方向隐式（ADI）有限差分格式,该模型是一个具有混合导数和非常数系数的对流扩散型初边值问题.我们将不同的高阶空间离散与时间步ADI分裂格式相结合,得到了一种空间四阶精度、时间二阶精度的有效方法,并采用Fourier方法分析了高阶ADI格式的稳定性.最后,通过对欧式看跌期权定价模型进行数值实验证实了数值方法的高阶收敛性.     

A finite element method for a two-sex model of population dynamics

A finite element scheme is described to approximate the solution of a nonlinear and non-local system of integro-differential equations that models the dynamics of a two-sex population. Crank-Nicolson time discretization is used and error estimates are derived for the appoximation.     

Numerical simulation of a class of fractional subdiffusion equations via the alternating direction implicit method

In this article, a new numerical technique is proposed for solving the two‐dimensional time fractional subdiffusion equation with nonhomogeneous terms. After a transformation of the original problem, standard central difference approximation is used for the spatial discretization. For the time step, a new fractional alternating direction implicit (FADI) scheme based on the L₁ approximation is considered. This FADI scheme is constructed by adding a small term, so it is different from standard FADI methods. The solvability, unconditional stability and H¹ norm convergence are proved. Finally, numerical examples show the effectiveness and accuracy of our proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 531–547, 2016     

Efficient Simulation of Wave Propagation with Implicit Finite Difference Schemes

Finite difference method is an important methodology in the approximation of waves. In this paper, we will study two implicit finite difference schemes for the simulation of waves. They are the weighted alternating direction implicit (ADI) scheme and the locally one-dimensional (LOD) scheme. The approximation errors, stability conditions, and dispersion relations for both schemes are investigated. Our analysis shows that the LOD implicit scheme has less dispersion error than that of the ADI scheme. Moreover, the unconditional stability for both schemes with arbitrary spatial accuracy is established for the first time. In order to improve computational efficiency, numerical algorithms based on message passing interface (MPI) are implemented. Numerical examples of wave propagation in a three-layer model and a standard complex model are presented. Our analysis and comparisons show that both ADI and LOD schemes are able to efficiently and accurately simulate wave propagation in complex media.