A Hermite pseudospectral solver for two-dimensional incompressible flows on infinite domains

The Hermite pseudospectral method is applied to solve the Navier-Stokes equations on a two-dimensional infinite domain. The feature of Hermite function allows us to adopt larger time steps than other spectral methods, but also leads to some extra computation when the stream function is calculated from the vorticity field. The scaling factor is employed to increase the resolution within the region of our main interest, and the aliasing error is fully removed by the 2/3-rule. Several traditional numerical experiments are performed with high accuracy, and some related future work on physical applications of this program is also discussed.