A Caltech Library Service

High-Order Unconditionally-Stable FC-AD PDE Solvers for General Domains


Lyon, Mark Edward (2009) High-Order Unconditionally-Stable FC-AD PDE Solvers for General Domains. Dissertation (Ph.D.), California Institute of Technology. doi:10.7907/3FFW-GK56.


A new methodology is introduced for the numerical solution of Partial Differential Equations in general spatial domains. The methodology is based on the use of the well-known Alternating Direction Implicit (ADI) approach of Peaceman and Rachford in conjunction with one-dimensional and high-order accurate Fourier representations of non-periodic data, obtained by way of a certain "continuation method" introduced recently for the resolution of the Gibbs phenomenon. We construct a number of high-order convergent PDE solvers on the basis of this strategy. Unlike previous alternating direction methods for general domains of order higher than one, the new algorithms possess the desirable property of unconditional stability for general spatial domains; the computational time required for these methods to advance one time-step, in turn, grows in an essentially linear manner with the number of spatial discretization points. In particular, the new methodology yields significant advantages over traditional low-order methods for computations involving wave propagation and "large domains," as well as PDEs including diffusive terms. In all, we treat Dirichlet problems for the Heat Equation, the Poisson Equation, and the Wave Equation in two- and three-dimensional spatial domains with smooth boundaries. A stability analysis we present hinges upon the numerical evaluation of certain singular value decompositions. Numerical results arising from the implementations of two- and three-dimensional versions of the method for linear parabolic, hyperbolic and elliptic PDEs, exhibit unconditional stability and high-order convergence, in agreement with our theoretical results.

Item Type:Thesis (Dissertation (Ph.D.))
Subject Keywords:ADI; Fourier continuation methed; Fourier series; implicit solvers; PDE solvers
Degree Grantor:California Institute of Technology
Division:Engineering and Applied Science
Major Option:Applied And Computational Mathematics
Thesis Availability:Public (worldwide access)
Research Advisor(s):
  • Bruno, Oscar P.
Thesis Committee:
  • Bruno, Oscar P. (chair)
  • Pullin, Dale Ian
  • Pierce, Niles A.
  • Schroeder, Peter
Defense Date:29 July 2008
Record Number:CaltechETD:etd-07312008-102751
Persistent URL:
Default Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:2992
Deposited By: Imported from ETD-db
Deposited On:01 Aug 2008
Last Modified:26 Nov 2019 19:14

Thesis Files

PDF (lyon_thesis_A100Final.pdf) - Final Version
See Usage Policy.


Repository Staff Only: item control page