The smoothing properties of the alternating direction implicit method in multigrid iterations

The performance of the alternating direction implicit (ADI) method within the multigrid context is examined. Local mode analysis is used to predict the smoothing properties of ADI for the Laplace problem and the biharmonic equation. The smoothing rate of ADI is shown to compare favourably with that of the Gauss-Seidel method. As a result, fewer iterations are required when ADI is used as the relaxation process in the multigrid method.     

交错网格紧致差分格式和满足等价性的压力Poisson方程

1．引言随着电子计算机的发展，越来越多的实际问题数值模拟成为现实，但还有很多非线性问题数值计算时间太长，内存要求过大．数值方法的改进可使计算量和存储量大大减少，例如，对二维非定常问题，要使误差达到N－4量级，二阶格式计算点数为（N2）3，（包括时间方向），而四阶格式计算点数仅为N3，差N3倍！而计算量差的倍数更多．当N＝16时N3＝4096，当N＝256时，N31678万．紧致差分格式具有精度高，差分式基点少，＜线性）稳定性好，对高频波分辨率高，边界差分点少等优点【’，‘，’，’。’，’，‘’］，本文中的格式基点数为3，…     

An exponential high‐order compact ADI method for 3D unsteady convection–diffusion problems

In this article, we develop an exponential high order compact alternating direction implicit (EHOC ADI) method for solving three dimensional (3D) unsteady convection–diffusion equations. The method, which requires only a regular seven‐point 3D stencil similar to that in the standard second‐order methods, is second order accurate in time and fourth‐order accurate in space and unconditionally stable. The resulting EHOC ADI scheme in each alternating direction implicit (ADI) solution step corresponding to a strictly diagonally dominant matrix equation can be solved by the application of the one‐dimensional tridiagonal Thomas algorithm with a considerable saving in computing time. Numerical experiments for three test problems are carried out to demonstrate the performance of the present method and to compare it with the classical Douglas–Gunn ADI method and the Karaa's high‐order compact ADI method. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

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.     

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.     

How Competitive is the ADI for Tensor Structured Equations?

In [1] we presented an extension of the alternating direction implicit (ADI) method for the solution of Lyapunov equations (1) to higher dimensional problems. The vectorized form of the Lyapunov equation is We considered the generalization of this equation of the form (2) The tensor train structure is one possible generalization of the low rank factorization we find in the right hand side of (1). Therefor we assume B to be of tensor train structure. We showed that in analogy to the low rank ADI case the solution X can be generated in tensor train structure, too. Further we provided an algorithm that computes X using a generalization of the ADI method. Here we compare our new tensor ADI method with an density matrix renormalization group (DMRG) solver for tensor train matrix equations and with matrix equation solvers to investigate the competitiveness of our new solver. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

ADI preconditioned Krylov methods for large Lyapunov matrix equations

In the present paper, we propose preconditioned Krylov methods for solving large Lyapunov matrix equations AX+XA^T+BB^T=0. Such problems appear in control theory, model reduction, circuit simulation and others. Using the Alternating Direction Implicit (ADI) iteration method, we transform the original Lyapunov equation to an equivalent symmetric Stein equation depending on some ADI parameters. We then define the Smith and the low rank ADI preconditioners. To solve the obtained Stein matrix equation, we apply the global Arnoldi method and get low rank approximate solutions. We give some theoretical results and report numerical tests to show the effectiveness of the proposed approaches.     

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 fast numerical solution of the 3D nonlinear tempered fractional integrodifferential equation

In this paper, we investigate the numerical solution of the three-dimensional (3D) nonlinear tempered fractional integrodifferential equation which is subject to the initial and boundary conditions. The backward Euler (BE) method in association with the first-order convolution quadrature rule is employed to discretize this equation for time, and the Galerkin finite element method is applied for space, which is combined with an alternating direction implicit (ADI) algorithm, in order to reduce the computational cost for solving the three-dimensional nonlocal problem. Then a fully discrete BE ADI Galerkin finite element scheme can be obtained by linearizing the non-linear term. Thereafter we prove a positive-type lemma, from which the stability and convergence of the proposed numerical scheme are derived based on the energy method. Numerical experiments are performed to verify the effectiveness of the proposed approach.     

Parameterized families of finite difference schemes for the wave equation

This article is devoted to an analysis of simple families of finite difference schemes for the wave equation. These families are dependent on several free parameters, and methods for obtaining stability bounds as a function of these parameters are discussed in detail. Access to explicit stability bounds such as those derived here may, it is hoped, lead to optimization techniques for so‐called spectral‐like methods, which are difference schemes dependent on many free parameters (and for which maximizing the order of accuracy may not be the defining criterion). Though the focus is on schemes for the wave equation in one dimension, the analysis techniques are extended to two dimensions; implicit schemes such as ADI methods are examined in detail. Numerical results are presented. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 463–480, 2004.     

High-order ADI orthogonal spline collocation method for a new 2D fractional integro-differential problem

We use the generalized L1 approximation for the Caputo fractional derivative, the second-order fractional quadrature rule approximation for the integral term, and a classical Crank-Nicolson alternating direction implicit (ADI) scheme for the time discretization of a new two-dimensional (2D) fractional integro-differential equation, in combination with a space discretization by an arbitrary-order orthogonal spline collocation (OSC) method. The stability of a Crank-Nicolson ADI OSC scheme is rigourously established, and error estimate is also derived. Finally, some numerical tests are given.     

Alternating Direction Implicit Galerkin Finite Element Method for the Two-Dimensional Time Fractional Evolution Equation

New numerical techniques are presented for the solution of the two-dimensional time fractional evolution equation in the unit square. In these methods, Galerkin finite element is used for the spatial discretization, and, for the time stepping, new alternating direction implicit (ADI) method based on the backward Euler method combined with the first order convolution quadrature approximating the integral term are considered. The ADI Galerkin finite element method is proved to be convergent in time and in the $L^2$ norm in space. The convergence order is$\mathcal{O}$($k$|ln $k$|+$h^r$), where $k$ is the temporal grid size and $h$ is spatial grid size in the $x$ and $y$ directions, respectively. Numerical results are presented to support our theoretical analysis.     

On the ADI method for Sylvester equations

This paper is concerned with the numerical solution of large scale Sylvester equations AX−XB=C, Lyapunov equations as a special case in particular included, with C having very small rank. For stable Lyapunov equations, Penzl (2000) [22] and Li and White (2002) [20] demonstrated that the so-called Cholesky factor ADI method with decent shift parameters can be very effective. In this paper we present a generalization of the Cholesky factor ADI method for Sylvester equations. An easily implementable extension of Penz’s shift strategy for the Lyapunov equation is presented for the current case. It is demonstrated that Galerkin projection via ADI subspaces often produces much more accurate solutions than ADI solutions.     

High order difference schemes for the system of two space second order nonlinear hyperbolic equations with variable coefficients

In this paper, we develop implicit difference schemes of O(k⁴ + k²h² + h⁴), where k > 0, h > 0 are grid sizes in time and space coordinates, respectively, for solving the system of two space dimensional second order nonlinear hyperbolic partial differential equations with variable coefficients having mixed derivatives subject to appropriate initial and boundary conditions. The proposed difference method for the scalar equation is applied for the solution of wave equation in polar coordinates to obtain three level conditionally stable ADI method of O(k⁴ + k²h² + h⁴). Some physical nonlinear problems are provided to demonstrate the accuracy of the implementation.     

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.     

A high‐order ADI finite difference scheme for a 3D reaction‐diffusion equation with neumann boundary condition

In this article, we extend the fourth‐order compact boundary scheme in Liao et al. (Numer Methods Partial Differential Equations 18 (2002), 340–354) to a 3D problem and then combine it with the fourth‐order compact alternating direction implicit (ADI) method in Gu et al. (J Comput Appl Math 155 (2003), 1–17) to solve the 3D reaction‐diffusion equation with Neumann boundary condition. First, the reaction‐diffusion equation is solved with a compact fourth‐order finite difference method based on the Padé approximation, which is then combined with the ADI method and a fourth‐order compact scheme to approximate the Neumann boundary condition, to obtain fourth order accuracy in space. The accuracy in the temporal dimension is improved to fourth order by applying the Richardson extrapolation technique, although the unconditional stability of the numerical method is proved, and several numerical examples are presented to demonstrate the accuracy and efficiency of the proposed new algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

解三维抛物型方程的一个ADI格式

构造了一个解三维抛物型方程的高精度ADI格式,格式绝对稳定,截断误差为O（△t^2＋△x^4）;然后应用Richerdson外推法,外推一次得到了具有O（△t^3＋△x^6）阶精度的近似解.     

A class of one-parameter alternating direction implicit methods for two-dimensional wave equations with discrete and distributed time-variable delays

For solving the initial-boundary value problem of two-dimensional wave equations with discrete and distributed time-variable delays, in the present paper, we first construct a class of basic one-parameter methods. In order to raise the computational efficiency of this class methods, we remold the methods as one-parameter alternating direction implicit (ADI) methods. Under the suitable conditions, the remolded methods are proved to be stable and convergent of second order in both of time and space. With several numerical experiments, the computational effectiveness and theoretical exactness of the methods are confirmed. Moreover, it is illustrated that the proposed one-parameter ADI method has the better advantage in computational efficiency than the basic one-parameter methods.     

Analysis of a fourth-order compact ADI method for a linear hyperbolic equation with three spatial variables

This paper is concerned with a three-level alternating direction implicit (ADI) method for the numerical solution of a 3D hyperbolic equation. Stability criterion of this ADI method is given by using von Neumann method. Meanwhile, it is shown by a discrete energy method that it can achieve fourth-order accuracy in both time and space with respect to H ¹- and L ²-norms only if stable condition is satisfied. It only needs solution of a tri-diagonal system at each time step, which can be solved by multiple applications of one-dimensional tri-diagonal algorithm. Numerical experiments confirming the high accuracy and efficiency of the new algorithm are provided.