Numerical Algorithms for Compressible Turbulence.
Abstract
This research involves higher order methods for solving partial differential equations. A common theme of this research is the use of spectral methods. One of the challenges of the spectral Fourier method is approximating smooth functions that are non-periodic. The rate of convergence deteriorates to O(1/N) away from the boundaries and there are spurious oscillations at the boundary points, known as the Gibbs phenomenon. We discuss an analytic but non-periodic function f(x) defined on -1,1, for which the first 2N+1 Fourier coefficients are known, and focuses on two methods that successfully eliminate the Gibbs phenomenon based on the information obtained from the Fourier coefficients. A new method is suggested to improve upon the results that have been obtained thus far. Also discussed is the Gibbs phenomenon for spherical harmonic spectral methods applied to functions that are piecewise analytic on spheres. We prove that knowledge of the first N spherical harmonic coefficients yield an exponentially convergent approximation to a spherical piecewise smooth function, hence completely overcoming the Gibbs phenomenon. Higher order numerical methods are applied to a specific one-dimensional hyperbolic system that describes the shallow water dynamics of a linearly sloping gulf. This problem is challenging due to the constantly changing location of where the water touches the shore. At this point the equation loses its strong hyperbolicity, and conventional numerical methods may not produce the correct solution. Also, any discontinuities that may arise as a result of the nonlinearity cannot be predicted, and this must be considered in selecting a numerical method.
Document Details
- Document Type
- Technical Report
- Publication Date
- Feb 26, 1997
- Accession Number
- ADA322928
Entities
People
- David Gottlieb
Organizations
- Brown University