TY - JOUR
T1 - A spectral solver for the Navier-Stokes equations in the velocity-vorticity formulation for flows with two non-periodic directions
AU - Clercx, H.J.H.
PY - 1997
Y1 - 1997
N2 - A novel pseudospectral solution procedure for the Navier–Stokes equations in the velocity–vorticity formulation, suitable to simulate flows with two nonperiodic directions, is proposed. An influence matrix method, including a tau correction procedure, has been employed to elicit an algorithmic substitute for thea priorilacking boundary conditions for the vorticity. Following O. Daube (J. Comput. Phys.103,402 (1992)) the influence matrix is built by enforcing either the definition of the vorticity or the continuity equation at the boundary of the domain. The influence matrix is nonsingular in both cases. The order of the influence matrix is twice larger than one would expect from analogous solution methods based on finite differences or finite elements. The spatial discretization is based on a 2D Chebyshev expansion on anonstaggeredgrid of collocation points, and the boundary conditions are imposed via the Lanczos tau procedure. The time marching scheme is Adams–Bashforth for the advection term and Crank–Nicolson for the viscous term. The proposed scheme yields machine accurate divergence–free flow fields, and the definition of the vorticity is satisfied within machine accuracy.
AB - A novel pseudospectral solution procedure for the Navier–Stokes equations in the velocity–vorticity formulation, suitable to simulate flows with two nonperiodic directions, is proposed. An influence matrix method, including a tau correction procedure, has been employed to elicit an algorithmic substitute for thea priorilacking boundary conditions for the vorticity. Following O. Daube (J. Comput. Phys.103,402 (1992)) the influence matrix is built by enforcing either the definition of the vorticity or the continuity equation at the boundary of the domain. The influence matrix is nonsingular in both cases. The order of the influence matrix is twice larger than one would expect from analogous solution methods based on finite differences or finite elements. The spatial discretization is based on a 2D Chebyshev expansion on anonstaggeredgrid of collocation points, and the boundary conditions are imposed via the Lanczos tau procedure. The time marching scheme is Adams–Bashforth for the advection term and Crank–Nicolson for the viscous term. The proposed scheme yields machine accurate divergence–free flow fields, and the definition of the vorticity is satisfied within machine accuracy.
U2 - 10.1006/jcph.1997.5799
DO - 10.1006/jcph.1997.5799
M3 - Article
VL - 137
SP - 186
EP - 211
JO - Journal of Computational Physics
JF - Journal of Computational Physics
SN - 0021-9991
ER -