Two alternating direction implicit spectral methods for two-dimensional distributed-order differential equation

Numerical Algorithms - In this paper, two alternating direction implicit Galerkin-Legendre spectral methods for distributed-order differential equation in two-dimensional space are developed. It is...     

Extrapolated alternating direction implicit iterative methods

This paper describes a new generalised (extrapolated) A.D.I. method for the solution of Laplace's equation. This method uses (i) a fixed acceleration parameter and (ii) the set of acceleration parameters of Douglas. The theory is applied to the 2-dimensional case and optimum numerical results are obtained.     

An extension ofA-stability to alternating direction implicit methods

Completely implicit, noniterative, finite-difference schemes have recently been developed by several authors for nonlinear, multidimensional systems of hyperbolic and mixed hyperbolic-parabolic partial differential equations. The method of Douglas and Gunn or the method of approximate factorization can be used to reduce the computational problem to a sequence of one-dimensional or alternating direction implicit (ADI) steps. Since the eigenvalues of partial differential equations (for example, the equations of compressible fluid dynamics) are often widely distributed with large imaginary parts,A-stable integration formulas provide ideal time-differencing approximations. In this paper it is shown that if anA-stable linear multistep method is used to integrate a model two-dimensional hyperbolic-parabolic partial differential equation, then one can always construct an ADI scheme by the method of approximate factorization which is alsoA-stable, i.e., unconditionally stable. A more restrictive result is given for three spatial dimensions. Since necessary and sufficient conditions forA-stability can easily be determined by using the theory of positive real functions, the stability analysis of the factored partial difference equations is reduced to a simple algebraic test.The main results of this paper were presented at the SIAM National Meeting, Madison, Wis., May 24 to 26, 1978, and section 9 was part of a presentation at the 751st Meeting of the American Mathematical Society, San Luis Obispo, California, Nov. 11 to 12, 1977.     

Solving semidefinite programming problems via alternating direction methods

In this paper, we consider an alternating direction algorithm for the solution of semidefinite programming problems (SDP). The main idea of our algorithm is that we reformulate the complementary conditions in the primal–dual optimality conditions as a projection equation. By using this reformulation, we only need to make one projection and solve a linear system of equation with reduced dimension in each iterate. We prove that the generated sequence converges to the solution of the SDP under weak conditions.     

The solution by fast Fourier transforms of Laplace's equation in a toroidal region with a rectangular cross-section

The fast Fourier transform method is described for Laplace's equation in a toroidal region using the 9-point difference approximation to the Laplacian operator. Numerical results are given which indicate the efficiency and accuracy of the method. Accurate difference approximations are also derived for the determination of the electrostatic field in a toroidal region.     

Numerical solutions of the Laplace’s equation

This paper uses the sinc methods to construct a solution of the Laplace’s equation using two solutions of the heat equation. A numerical approximation is obtained with an exponential accuracy. We also present a reliable algorithm of Adomian decomposition method to construct a numerical solution of the Laplace’s equation in the form a rapidly convergence series and not at grid points. Numerical examples are given and comparisons are made to the sinc solution with the Adomian decomposition method. The comparison shows that the Adomian decomposition method is efficient and easy to use.     

The new alternating direction implicit difference methods for the wave equations

A new second-order alternating direction implicit (ADI) scheme, based on the idea of the operator splitting, is presented for solving two-dimensional wave equations. The scheme is also extended to a high-order compact difference scheme. Both of them have the advantages of unconditional stability, less impact of the perturbing terms on the accuracy, and being convenient to compute the boundary values of the intermediates. Besides this, the compact scheme has high-order accuracy and costs less in computational time. Numerical examples are presented and the results are very satisfactory.     

Layer potentials and regularity for the Dirichlet problem for Laplace's equation in Lipschitz domains

For D, a bounded Lipschitz domain in Rⁿ, n ? 2, the classical layer potentials for Laplace's equation are shown to be invertible operators on L²(?D) and various subspaces of L²(?D). For 1 < p ? 2 and data in L^p(?D) with first derivatives in L^p(?D) it is shown that there exists a unique harmonic function, u, that solves the Dirichlet problem for the given data and such that the nontangential maximal function of ▽u is in L^p(?D). When n = 2 the question of the invertibility of the layer potentials on every L^p(?D), 1 < p < ∞, is answered.     

The new alternating direction implicit difference methods for solving three-dimensional parabolic equations

A new alternating direction implicit (ADI) scheme for solving three-dimensional parabolic equations with nonhomogeneous boundary conditions is presented. The scheme is also extended to high-order compact difference scheme. Both of them have the advantages of unconditional stability and being convenient to compute the boundary values of the intermediates. Besides this, the compact scheme has high-order accuracy and uses less computational time. Numerical examples are presented and the results are very satisfactory.     

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.     

An unconditionally stable alternating direction implicit scheme for the two space dimensional linear hyperbolic equation

We report a new unconditionally stable implicit alternating direction implicit (ADI) scheme of O(k² + h²) for the difference solution of linear hyperbolic equation u_tt + 2αu_t + β²u = u_xx + u_yy + f(x, y, t), αβ ≥ 0, 0 < x, y < 1, t > 0 subject to appropriate initial and Dirichlet boundary conditions, where α > 0 and β ≥ 0 are real numbers. The resulting system of algebraic equations is solved by split method. Numerical results are provided to demonstrate the efficiency and accuracy of the method. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 684–688, 2001     

Convergence analysis of high-order compact alternating direction implicit schemes for the two-dimensional time fractional diffusion equation

High-order compact finite difference method with operator-splitting technique for solving the two dimensional time fractional diffusion equation is considered in this paper. The Caputo derivative is evaluated by the L1 approximation, and the second order derivatives with respect to the space variables are approximated by the compact finite differences to obtain fully discrete implicit schemes. Alternating Direction Implicit (ADI) method is used to split the original problem into two separate one dimensional problems. One scheme is given by replacing the unknowns by the values on the previous level directly and a correction term is added for another scheme. Theoretical analysis for the first scheme is discussed. The local truncation error is analyzed and the stability is proved by the Fourier method. Using the energy method, the convergence of the compact finite difference scheme is proved. Numerical results are provided to verify the accuracy and efficiency of the two proposed algorithms. For the order of the temporal derivative lies in different intervals $\left(0,\frac{1}{2}\right)$ or $\left[\frac{1}{2},1\right)$ , corresponding appropriate scheme is suggested.     

Fourth-order alternating direction implicit compact finite difference schemes for two-dimensional Schrödinger equations

In this paper, alternating direction implicit compact finite difference schemes are devised for the numerical solution of two-dimensional Schrödinger equations. The convergence rates of the present schemes are of order O(h⁴+τ²). Numerical experiments show that these schemes preserve the conservation laws of charge and energy and achieve the expected convergence rates. Representative simulations show that the proposed schemes are applicable to problems of engineering interest and competitive when compared to other existing procedures.     

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.