Solution of Evolutionary Partial Differential Equations Using Adaptive Finite Differences with Pseudospectral Post-processing

A coordinate transformation approach is described that enables pseudospectral methods to be applied efficiently to unsteady differential problems with steep solutions. The work is an extension of a method presented by Mulholland, Huang, and Sloan for the adaptive pseudospectral solution of steady problems. A coarse grid is generated by a moving mesh finite difference method that is based on equidistribution, and this grid is used to construct a time-dependent coordinate transformation. A sequence of spatial transformations may be generated at discrete points in time, or a single transformation may be generated as a continuous function of space and time. The differential problem is transformed by the coordinate transformation and then solved using a method that combines pseudospectral discretisation in space with a suitable integrator in time. Numerical results are presented for unsteady problems in one space dimension.

[1]  Keith Miller,et al.  Moving Finite Elements. I , 1981 .

[2]  Bernard J. Matkowsky,et al.  Adaptive approximation of solutions to problems with multiple layers by Chebyshev pseudo-spectral methods , 1995 .

[3]  Weizhang Huang,et al.  Moving mesh partial differential equations (MMPDES) based on the equidistribution principle , 1994 .

[4]  Weizhang Huang,et al.  Moving Mesh Methods Based on Moving Mesh Partial Differential Equations , 1994 .

[5]  Joe F. Thompson,et al.  Numerical grid generation , 1985 .

[6]  Weizhang Huang,et al.  A Simple Adaptive Grid Method in Two Dimensions , 1994, SIAM J. Sci. Comput..

[7]  Joseph E. Flaherty,et al.  On the stability of mesh equidistribution strategies for time-dependent partial differential equations , 1986 .

[8]  J. M. Sanz-Serna,et al.  On simple moving grid methods for one-dimensional evolutionary partial differential equations , 1988 .

[9]  Ivo Babuška,et al.  An adaptive method of lines with error control for parabolic equations of the reaction-diffusion type , 1986 .

[10]  J. Lambert Numerical Methods for Ordinary Differential Systems: The Initial Value Problem , 1991 .

[11]  Weizhang Huang,et al.  Pseudospectral Solution of Near-Singular Problems using Numerical Coordinate Transformations Based on Adaptivity , 1998, SIAM J. Sci. Comput..

[12]  E. Dorfi,et al.  Simple adaptive grids for 1-d initial value problems , 1987 .

[13]  G. Whitham,et al.  Linear and Nonlinear Waves , 1976 .

[14]  I. Christie,et al.  A simple adaptive technique for nonlinear wave problems , 1986 .

[15]  B. Fornberg,et al.  A review of pseudospectral methods for solving partial differential equations , 1994, Acta Numerica.

[16]  T. A. Zang,et al.  Spectral methods for fluid dynamics , 1987 .