Pseudospectral algorithms for Navier-Stokes simulation of turbulent flows in cylindrical geometry with coordinate singularities (Q1891116)

A family of algorithms for incompressible three-dimensional Navier-Stokes equations in cylindrical geometry is considered. A Galerkin trigonometric approximation for the uniform variables and a pseudospectral polynomial approximation for the radial variable \(r\) are used. The time discretization is performed by implicit and predictor-corrector time advancement schemes. In order to achieve high spatial accuracy at the coordinate singularity \(r = 0\), the behaviour of analytic functions in the vicinity \(r = 0\) is exploited. Further different time-splitting techniques are proposed, and a new set of collocation nodes is presented which possesses the discrete energy conservation law analogous to the continuous energy conservation law for the original Navier-Stokes equations. This is of special significance for the nonlinear stability of algorithms.