An axis treatment for flow equations in cylindrical coordinates based on parity conditions

Abstract A novel axis treatment using parity conditions is presented for flow equations in cylindrical coordinates that are represented in azimuthal Fourier modes. The correct parity states of scalars and the velocity vector are derived such that symmetry conditions for each Fourier mode of the respective variable can be determined. These symmetries are then used to construct finite-difference and filter stencils at and near the axis, and an interpolation scheme for the computation of terms premultiplied by 1/r. A grid convergence study demonstrates that the axis treatment retains the formal accuracy of the spatial discretization scheme employed. Two further test cases are considered for evaluation of the axis treatment, the propagation of an acoustic pulse and direct numerical simulation of a fully turbulent supersonic axisymmetric wake. The results demonstrate the applicability of the axis treatment for non-axisymmetric flows.

[1]  R. Lewis,et al.  Low-storage, Explicit Runge-Kutta Schemes for the Compressible Navier-Stokes Equations , 2000 .

[2]  Richard D. Sandberg,et al.  DNS OF COMPRESSIBLE PIPE FLOW EXITING INTO A COFLOW , 2012, Proceeding of Seventh International Symposium on Turbulence and Shear Flow Phenomena.

[3]  John D. Anderson,et al.  A computational fluid dynamic technique valid at the centerline for non-axisymmetric problems in cylindrical coordinates , 1979 .

[4]  Richard D. Sandberg,et al.  Numerical investigation of transitional supersonic axisymmetric wakes , 2006, Journal of Fluid Mechanics.

[5]  J. Kim,et al.  Optimized Compact Finite Difference Schemes with Maximum Resolution , 1996 .

[6]  Christophe Bailly,et al.  A shock-capturing methodology based on adaptative spatial filtering for high-order non-linear computations , 2009, J. Comput. Phys..

[7]  Parviz Moin,et al.  Direct computation of Mach wave radiation in an axisymmetric supersonic jet , 1997 .

[8]  T. Colonius,et al.  Numerical Treatment of Polar Coordinate Singularities , 2000 .

[9]  Christopher A. Kennedy,et al.  Reduced aliasing formulations of the convective terms within the Navier-Stokes equations for a compressible fluid , 2008, J. Comput. Phys..

[10]  Navier-Stokes Spectral Solver in a Finite Cylinder , 2010 .

[11]  Sanjiva K. Lele,et al.  A highly accurate technique for the treatment of flow equations at the polar axis in cylindrical coordinates using series expansions , 2002 .

[12]  Jie Shen,et al.  Efficient Spectral-Galerkin Methods III: Polar and Cylindrical Geometries , 1997, SIAM J. Sci. Comput..

[13]  Neil D. Sandham,et al.  Large-eddy simulation of low-frequency unsteadiness in a turbulent shock-induced separation bubble , 2009 .

[14]  M. Visbal.  VERY HIGH-ORDER SPATIALLY IMPLICIT SCHEMES FOR COMPUTATIONAL ACOUSTICS ON CURVILINEAR MESHES , 2001 .

[15]  Paul Bellan,et al.  Physical constraints on the coefficients of Fourier expansions in cylindrical coordinates , 1990 .

[16]  P. Moin,et al.  A direct numerical simulation study on the mean velocity characteristics in turbulent pipe flow , 2008, Journal of Fluid Mechanics.

[17]  Eric Serre,et al.  A three-dimensional pseudospectral method for rotating flows in a cylinder , 2001 .

[18]  C. Tam,et al.  Dispersion-relation-preserving finite difference schemes for computational acoustics , 1993 .

[19]  Jae Wook Kim.  Optimised boundary compact finite difference schemes for computational aeroacoustics , 2007, J. Comput. Phys..

[20]  Richard D. Sandberg,et al.  DNS of a canonical compressible nozzle flow , 2011 .

[21]  Duck-Joo Lee,et al.  Generalized Characteristic Boundary Conditions for Computational Aeroacoustics, Part 2 , 2000 .

[22]  Jae Wook Kim,et al.  High-order compact filters with variable cut-off wavenumber and stable boundary treatment , 2010 .