A hybrid numerical simulation of isotropic compressible turbulence

A novel hybrid numerical scheme with built-in hyperviscosity has been developed to address the accuracy and numerical instability in numerical simulation of isotropic compressible turbulence in a periodic domain at high turbulent Mach number. The hybrid scheme utilizes a 7th-order WENO (Weighted Essentially Non-Oscillatory) scheme for highly compressive regions (i.e., shocklet regions) and an 8th-order compact central finite difference scheme for smooth regions outside shocklets. A flux-based conservative and formally consistent formulation is developed to optimize the connection between the two schemes at the interface and to achieve a higher computational efficiency. In addition, a novel numerical hyperviscosity formulation is proposed within the context of compact finite difference scheme for the smooth regions to improve numerical stability of the hybrid method. A thorough and insightful analysis of the hyperviscosity formulation in both Fourier space and physical space is presented to show the effectiveness of the formulation in improving numerical stability, without compromising the accuracy of the hybrid method. A conservative implementation of the hyperviscosity formulation is also developed. Combining the analysis and test simulations, we have also developed a criterion to guide the specification of a numerical hyperviscosity coefficient (the only adjustable coefficient in the formulation). A series of test simulations are used to demonstrate the accuracy and numerical stability of the scheme for both decaying and forced compressible turbulence. Preliminary results for a high-resolution simulation at turbulent Mach number of 1.08 are shown. The sensitivity of the simulated flow to the detail of thermal forcing method is also briefly discussed.

[1]  T. Mexia,et al.  Author ' s personal copy , 2009 .

[2]  Sergio Pirozzoli,et al.  Direct numerical simulations of isotropic compressible turbulence: Influence of compressibility on dynamics and structures , 2004 .

[3]  M. Pino Martín,et al.  Stencil Adaptation Properties of a WENO Scheme in Direct Numerical Simulations of Compressible Turbulence , 2007, J. Sci. Comput..

[4]  Parviz Moin,et al.  Computational issues and algorithm assessment for shock/turbulence interaction problems , 2007 .

[5]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[6]  P. Woodward,et al.  Measures of intermittency in driven supersonic flows. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[7]  A. Pouquet,et al.  A Turbulent Model for the Interstellar Medium. II. Magnetic Fields and Rotation , 1994 .

[8]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[9]  M. Pino Martín,et al.  Optimization of nonlinear error for weighted essentially non-oscillatory methods in direct numerical simulations of compressible turbulence , 2007, J. Comput. Phys..

[10]  D. Pullin,et al.  Hybrid tuned center-difference-WENO method for large eddy simulations in the presence of strong shocks , 2004 .

[11]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[12]  P. Woodward,et al.  Inertial range structures in decaying compressible turbulent flows , 1998 .

[13]  Krishnan Mahesh,et al.  High order finite difference schemes with good spectral resolution , 1997 .

[14]  Shiyi Chen,et al.  On statistical correlations between velocity increments and locally averaged dissipation in homogeneous turbulence , 1993 .

[15]  Shiyi Chen,et al.  Examination of hypotheses in the Kolmogorov refined turbulence theory through high-resolution simulations. Part 1. Velocity field , 1996, Journal of Fluid Mechanics.

[16]  Chi-Wang Shu,et al.  Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy , 2000 .

[17]  W. Sutherland.  LII. The viscosity of gases and molecular force , 1893 .

[18]  Yuxin Ren,et al.  A characteristic-wise hybrid compact-WENO scheme for solving hyperbolic conservation laws , 2003 .

[19]  Steven A. Orszag,et al.  Energy and spectral dynamics in forced compressible turbulence , 1990 .

[20]  A. G. Bashkirov.  Influence of viscous stresses on shock wave stability in gases , 1991 .

[21]  P. Moin,et al.  DIRECT NUMERICAL SIMULATION: A Tool in Turbulence Research , 1998 .

[22]  Ravi Samtaney,et al.  Direct numerical simulation of decaying compressible turbulence and shocklet statistics , 2001 .

[23]  T. A. Zang,et al.  Direct and large-eddy simulations of three-dimensional compressible Navier-Stokes turbulence , 1992 .

[24]  Sergio Pirozzoli,et al.  Conservative Hybrid Compact-WENO Schemes for Shock-Turbulence Interaction , 2002 .

[25]  Johan Larsson,et al.  Stability criteria for hybrid difference methods , 2008, J. Comput. Phys..

[26]  Feng He,et al.  A new family of high-order compact upwind difference schemes with good spectral resolution , 2007, J. Comput. Phys..

[27]  V. Gregory Weirs,et al.  A bandwidth-optimized WENO scheme for the effective direct numerical simulation of compressible turbulence , 2006, J. Comput. Phys..

[28]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[29]  Tim Colonius,et al.  Acoustic Saturation in Bubbly Cavitating Flow Adjacent to an Oscillating Wall , 2000 .

[30]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[31]  Tapan K. Sengupta,et al.  High Accuracy Compact Schemes and Gibbs' Phenomenon , 2004, J. Sci. Comput..

[32]  Jang-Hyuk Kwon,et al.  A high-order accurate hybrid scheme using a central flux scheme and a WENO scheme for compressible flowfield analysis , 2005 .

[33]  Parviz Moin,et al.  EDDY SHOCKLETS IN DECAYING COMPRESSIBLE TURBULENCE , 1991 .

[34]  Tapan K. Sengupta,et al.  A new flux-vector splitting compact finite volume scheme , 2005 .

[35]  John A. Ekaterinaris,et al.  Regular Article: Implicit, High-Resolution, Compact Schemes for Gas Dynamics and Aeroacoustics , 1999 .

[36]  Nikolaus A. Adams,et al.  A High-Resolution Hybrid Compact-ENO Scheme for Shock-Turbulence Interaction Problems , 1996 .

[37]  Tapan K. Sengupta,et al.  Analysis of central and upwind compact schemes , 2003 .