A Caltech Library Service

Computing Nearly Singular Solutions Using Pseudo-Spectral Methods

Hou, Thomas Y. and Li, Ruo (2007) Computing Nearly Singular Solutions Using Pseudo-Spectral Methods. . (Submitted)

[img] PDF - Submitted Version
See Usage Policy.


Use this Persistent URL to link to this item:


In this paper, we investigate the performance of pseudo-spectral methods in computing nearly singular solutions of fluid dynamics equations. We consider two different ways of removing the aliasing errors in a pseudo-spectral method. The first one is the traditional 2/3 dealiasing rule. The second one is a high (36th) order Fourier smoothing which keeps a significant portion of the Fourier modes beyond the 2/3 cut-off point in the Fourier spectrum for the 2/3 dealiasing method. Both the 1D Burgers equation and the 3D incompressible Euler equations are considered. We demonstrate that the pseudo-spectral method with the high order Fourier smoothing gives a much better performance than the pseudo-spectral method with the 2/3 dealiasing rule. Moreover, we show that the high order Fourier smoothing method captures about 12∼15 more effective Fourier modes in each dimension than the 2/3 dealiasing method. For the 3D Euler equations, the gain in the effective Fourier codes for the high order Fourier smoothing method can be as large as 20% over the 2/3 dealiasing method. Another interesting observation is that the error produced by the high order Fourier smoothing method is highly localized near the region where the solution is most singular, while the 2/3 dealiasing method tends to produce oscillations in the entire domain. The high order Fourier smoothing method is also found be very stable dynamically. No high frequency instability has been observed.

Item Type:Report or Paper (Discussion Paper)
Related URLs:
URLURL TypeDescription
Additional Information:Submitted on 12 Jan 2007. February 2, 2008. We would like to thank Prof. Lin-Bo Zhang from the Institute of Computational Mathematics in Chinese Academy of Sciences (CAS) for providing us with the computing resource to perform this large scale computational project. Additional computing resource was provided by the Center of High Performance Computing in CAS. We also thank Prof. Robert Kerr for providing us with his Fortran subroutine that generates his initial data. This work was in part supported by NSF under the NSF FRG grant DMS-0353838 and ITR Grant ACI-0204932. Part of this work was done while Hou visited the Academy of Systems and Mathematical Sciences of CAS in the summer of 2005 as a member of the Oversea Outstanding Research Team for Complex Systems. Li was supported by the National Basic Research Program of China under the grant 2005CB321701. Finally, we would like to thank Professors Alfio Quarteroni, Jie Shen, and Eitan Tadmor for their valuable comments on our draft manuscript.
Funding AgencyGrant Number
National Basic Research Program of China2005CB321701
Record Number:CaltechAUTHORS:20160322-083351617
Persistent URL:
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:65572
Deposited By: Ruth Sustaita
Deposited On:30 Mar 2016 23:56
Last Modified:03 Oct 2019 09:48

Repository Staff Only: item control page