Solution of Time-Dependent PDE Through Component-Wise Approximation of Matrix Functions

Document Type

Article

Publication Date

2-10-2011

Department

Mathematics

School

Mathematics and Natural Sciences

Abstract

Block Krylov subspace spectral (KSS) methods are a "best-of-both-worlds" compromise between explicit and implicit time-stepping methods for variable-coefficient PDE, in that they combine the efficiency of explicit methods and the stability of implicit methods, while also achieving spectral accuracy in space and high-orer accuracy in time. Block KSS methods compute each Fourier coeffience of the solution usin techniques developed by Gene Golub and Gérard Meurant for approximating elements of functions of matrices by block Gaussian quadrature in te spectral, rather than physical, domain. This paper demonstrates the superiority of block KSS methods, in terms of accuracy and efficiency, to tother Krylov subspace methods in the literature. It is also described how the ideas behind block KSS methods can be applied to a variety of equations, including problems for which Fourier spectral methods are not normally feasible. In paricular, the versatility of the approach behind block KSS methods is demondrated through application to nonlinear diffusion equations for signal and image processing, and adaptiation to finite element discretization.

Publication Title

IAENG International Journal of Applied Mathematics

Volume

41

Issue

1

First Page

1

Last Page

10

Link to Full Text

Find in your library

Share

COinS