Abstract: A new methodology for the numerical solution of Partial Differential Equations (PDEs) in smoothly bounded domains will be presented. Our algorithms are based on the use of a certain ``Fourier Continuation'' (FC) method for the resolution of the Gibbs phenomenon in conjunction with well-known alternating direction (AD) strategies. This FC-AD methodology will be demonstrated through a variety of examples including the Heat, Laplace, Wave Equations and both Dirichlet and Neumann problems. The high-order algorithms possess the desirable property of unconditional stability for general domains in computational time that grows in an essentially linear manner with the number of unknowns. The significant improvements that these new algorithms can provide over the computing times and memory required by alternative general-domain solvers will be discussed. In particular, the debilitating spatial ``pollution error'', which arises as finite-difference and finite-element solvers are applied to the solution of wave propagation problems, is essentially absent from our calculations due to the Fourier basis used in the FC-AD calculations.
This talk will be accessible to graduate students.