A polynomial based iterative method for linear parabolic equations
A new polynomial based method (PBM) is developed to integrate multi-dimensional linear parabolic initial-boundary-value problems. It is based on L/sub 2/-approximations to f(z) = (1 - exp(-z))/z, f(0) = 1, over ellipses in the complex plane using expansions of f in Chebychev polynomials. The calculation of the Fourier coefficients requires numerical integration over only a single line segment in the complex plane whose length and orientation depend on the step size and the parabolic operator itself. The simplicity with which these coefficients are obtained rests on special properties of the Chebychev polynomials. Most of the work in PBM consists of matrix-vector multiplications, involving a matrix L which arises from the spatial discretization of the differential operator. To be specific, PBM integrates the semi-discrete problem u/sub t/ = L(t)u + b(t), u,b in R/sup n/ and L in R/sup n x n/, and requires only a modest amount of storage (a few vectors of order n). Due to the analyticity of f it has good convergence properties and compares favorably to other standard methods from the classes of Hopscotch, Alternating Direction Implicit (ADI) and Locally One-Dimensional (LOD) schemes, as measured by the CPU-times required on a single CPU of a CRAY X-MP/24. It is also competitive with Crank-Nicolson which we couple with two proven iterative solvers. I recommend PBM on problems which require fourth order spatial accuracy, problems whose solutions contain significant high-frequency components, and problems whose operators cannot be split conveniently in an ADI or LOD fashion (for example, problems with mixed derivatives). 30 refs., 11 figs., 12 tabs.
- Research Organization:
- Illinois Univ., Urbana (USA). Center for Supercomputing Research and Development
- DOE Contract Number:
- FG02-85ER25001
- OSTI ID:
- 5635473
- Report Number(s):
- DOE/ER/25001-91; CSRD-661; UILU-ENG-86-8003-05/87; ON: DE88003535
- Resource Relation:
- Other Information: Thesis (Ph.D.). Portions of this document are illegible in microfiche products. Thesis
- Country of Publication:
- United States
- Language:
- English
Similar Records
Superconvergence analysis for nonlinear parabolic equation with EQ{sub 1}{sup rot} nonconforming finite element
Density of states, Potts zeros, and Fisher zeros of the Q
Related Subjects
BOUNDARY-VALUE PROBLEMS
NUMERICAL SOLUTION
CRAY COMPUTERS
DIFFERENTIAL EQUATIONS
ITERATIVE METHODS
LEAST SQUARE FIT
LINEAR PROGRAMMING
MANY-DIMENSIONAL CALCULATIONS
POLYNOMIALS
TIME DEPENDENCE
COMPUTERS
EQUATIONS
FUNCTIONS
MAXIMUM-LIKELIHOOD FIT
PROGRAMMING
990220* - Computers
Computerized Models
& Computer Programs- (1987-1989)