The WR-HSS iteration method for a system of linear differential equations and its applications to the unsteady discrete elliptic problem

TL;DR

This paper introduces the WR-HSS iteration method, combining waveform relaxation and Hermitian/skew-Hermitian splitting, for efficiently solving non-self-adjoint positive definite linear differential equations and their application to unsteady elliptic problems.

Contribution

It develops a new WR-HSS iterative method with proven unconditional convergence for non-self-adjoint linear differential equations, extending the ADI and HSS techniques.

Findings

The WR-HSS method converges unconditionally.

The contraction factor depends only on the Hermitian part.

Applications demonstrate effectiveness and theoretical correctness.

Abstract

We consider the numerical method for non-self-adjoint positive definite linear differential equations, and its application to the unsteady discrete elliptic problem, which is derived from spatial discretization of the unsteady elliptic problem with Dirichlet boundary condition. Based on the idea of the alternating direction implicit (ADI) iteration technique and the Hermitian/skew-Hermitian splitting (HSS), we establish a waveform relaxation (WR) iteration method for solving the non-self-adjoint positive definite linear differential equations, called the WR-HSS method. We analyze the convergence property of the WR-HSS method, and prove that the WR-HSS method is unconditionally convergent to the solution of the system of linear differential equations. In addition, we derive the upper bound of the contraction factor of the WR-HSS method in each iteration which is only dependent on the Hermitian part of the corresponding non-self-adjoint positive definite linear differential operator. Finally, the applications of the WR-HSS method to the unsteady discrete elliptic problem demonstrate its effectiveness and the correctness of the theoretical results.