Commun. Comput. Phys., 8 (2010), pp. 663-689.


Navier-Stokes Spectral Solver in a Finite Cylinder

F. Auteri 1*, M. Biava 1, L. Quartapelle 1

1 Politecnico di Milano, Dipartimento di Ingegneria Aerospaziale, Via La Masa 34, 20156 Milano, Italy.

Received 19 August 2009; Accepted (in revised version) 14 December 2009
Available online 20 April 2010
doi:10.4208/cicp.190809.141209a

Abstract

A primitive variable spectral method for simulating incompressible viscous flows inside a finite cylinder is presented. One element of originality of the proposed method is that the radial discretization of the Fourier coefficients depends on the Fourier mode, its dimension decreasing with the increase of the azimuthal modal number. This principle was introduced independently by Matsushima and Marcus and by Verkley in polar coordinates and is adopted here for the first time to formulate a 3D cylindrical Galerkin projection method. A second element of originality is the use of a special basis of Jacobi polynomials introduced recently for the radial dependence in the solution of Dirichlet problems. In this basis the radial operators are represented by matrices of minimal sparsity - diagonal stiffness and tridiagonal mass - provided here in closed form for the first time, and lead to a Helmholtz operator characterized by a favorable condition number. Finally, a new method is presented for eliminating the singular behaviour of the solution originated by the rotation of the lid with respect to the cylindrical wall. Thanks to these elements, the resulting Navier-Stokes spectral solver guarantees the differentiability to any order of the solution in the entire computational domain and does not suffer from the time-step stability restriction occurring in spectral methods with a point clustering close to the axis. Several test examples are offered that demonstrate the spectral accuracy of the solution method under different representative conditions.

AMS subject classifications: 65N30, 65N35

Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: Navier-Stokes equations, finite cylindrical domain, spectral methods, Jacobi and Legendre polynomials, primitive variables, projection method.

*Corresponding author.
Email: auteri@aero.polimi.it (F. Auteri)
 

The Global Science Journal