Pseudospectral solution of the two-dimensional Navier-Stokes equations in a disk
An efficient and accurate algorithm for solving the two-dimensional (2D) incompressible Navier-Stokes equations on a disk with no-slip boundary conditions is described. The vorticity-stream function formulation of these equations is used, and spatially the vorticity and stream functions are expressed as Fourier-Chebyshev expansions. The Poisson and Helmholtz equations which arise from the implicit-explicit time marching scheme are solved as banded systems using a postconditioned spectral {tau}-method. The polar coordinate singularity is handled by expanding fields radically over the entire diameter using a parity modified Chebyshev series and building partial regularity into the vorticity. The no-slip boundary condition is enforced by transferring one of the two boundary conditions imposed on the stream function onto the vorticity via a solvability constraint. Significant grains in run times were realized by parallelizing the code in message passage interface (MPI).
- Research Organization:
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- OSTI ID:
- 20015657
- Journal Information:
- SIAM Journal on Scientific Computing, Vol. 21, Issue 1; Other Information: PBD: Sep 1999; ISSN 1064-8275
- Country of Publication:
- United States
- Language:
- English
Similar Records
CCM Continuity Constraint Method: A finite-element computational fluid dynamics algorithm for incompressible Navier-Stokes fluid flows
Incompressible computational fluid dynamics and the continuity constraint method for the three-dimensional Navier-Stokes equations