A fully discrete, kinetic energy consistent finite-volume scheme for compressible flows

A robust, implicit, low-dissipation method suitable for LES/DNS of compressible turbulent flows is discussed. The scheme is designed such that the discrete flux of kinetic energy and its rate of change are consistent with those predicted by the momentum and continuity equations. The resulting spatial fluxes are similar to those derived using the so-called skew-symmetric formulation of the convective terms. Enforcing consistency for the time derivative results in a novel density weighted Crank-Nicolson type scheme. The method is stable without the addition of any explicit dissipation terms at very high Reynolds numbers for flows without shocks. Shock capturing is achieved by switching on a dissipative flux term which tends to zero in smooth regions of the flow. Numerical examples include a one-dimensional shock tube problem, the Taylor-Green problem, simulations of isotropic turbulence, hypersonic flow over a double-cone geometry, and compressible turbulent channel flow.

[1]  Robert D. Moser,et al.  A numerical study of turbulent supersonic isothermal-wall channel flow , 1995, Journal of Fluid Mechanics.

[2]  Graham V. Candler,et al.  Effect of numerics on navier-stokes computations of hypersonic double-cone flows , 2005 .

[3]  Parviz Moin,et al.  A semi-implicit method for resolution of acoustic waves in low Mach number flows , 2002 .

[4]  E. Tadmor.  Skew-selfadjoint form for systems of conservation laws , 1984 .

[5]  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..

[6]  Joseph Oliger,et al.  Energy and Maximum Norm Es-timates for Nonlinear Conservation Laws , 1994 .

[7]  Björn Sjögreen,et al.  Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for High Order Methods , 2004, J. Sci. Comput..

[8]  Timothy Wadhams,et al.  CODE VALIDATION STUDY OF LAMINAR SHOCKIBOUNDARY LAYER AND SHOCK/SHOCK INTERACTIONS IN HYPERSONIC FLOW Part B: Comparison \\ith Navier-Stokes and DSMC Solutions , 2001 .

[9]  Nikolaus A. Adams,et al.  Large-Eddy Simulation of Shock-Turbulence Interaction , 2004 .

[10]  Ralf Deiterding,et al.  A low-numerical dissipation, patch-based adaptive-mesh-refinement method for large-eddy simulation of compressible flows , 2006 .

[11]  Graham V. Candler,et al.  Data-parallel lower-upper relaxation method for the navier-stokes equations , 1996 .

[12]  Gregory A. Blaisdell,et al.  The effect of the formulation of nonlinear terms on aliasing errors in spectral methods , 1996 .

[13]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

[14]  R. Rogallo.  Numerical experiments in homogeneous turbulence , 1981 .

[15]  A. Harten.  On the symmetric form of systems of conservation laws with entropy , 1983 .

[16]  Graham V. Candler,et al.  The solution of the Navier-Stokes equations using Gauss-Seidel line relaxation , 1989 .

[17]  Eitan Tadmor,et al.  The numerical viscosity of entropy stable schemes for systems of conservation laws. I , 1987 .

[18]  Antony Jameson,et al.  Formulation of Kinetic Energy Preserving Conservative Schemes for Gas Dynamics and Direct Numerical Simulation of One-Dimensional Viscous Compressible Flow in a Shock Tube Using Entropy and Kinetic Energy Preserving Schemes , 2008, J. Sci. Comput..

[19]  Vincent Guinot,et al.  High-Order Fluxes for Conservative Skew-Symmetric-like Schemes in Structured Meshes , 2000 .

[20]  Steven A. Orszag,et al.  On the Elimination of Aliasing in Finite-Difference Schemes by Filtering High-Wavenumber Components , 1971 .

[21]  Neil D. Sandham,et al.  Low-Dissipative High-Order Shock-Capturing Methods Using Characteristic-Based Filters , 1999 .

[22]  M. Mock,et al.  Systems of conservation laws of mixed type , 1980 .

[23]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[24]  P. Moin,et al.  Simulation of spatially evolving turbulence and the applicability of Taylor's hypothesis in compressible flow , 1992 .

[25]  Krishnan Mahesh,et al.  A robust, colocated, implicit algorithm for direct numerical simulation of compressible, turbulent flows , 2005 .

[26]  V. Guinot.  Approximate Riemann Solvers , 2010 .

[27]  Parviz Moin,et al.  Higher entropy conservation and numerical stability of compressible turbulence simulations , 2004 .

[28]  I. Nompelis.  Computational study of hypersonic double-cone experiments for code validation , 2004 .

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

[30]  F. Nicoud,et al.  Large-Eddy Simulation of the Shock/Turbulence Interaction , 1999 .

[31]  Oh Joon Kwon,et al.  An efficient and robust implicit operator for upwind point Gauss-Seidel method , 2007, J. Comput. Phys..

[32]  C. Pierce,et al.  Progress-variable approach for large-eddy simulation of turbulent combustion , 2001 .

[33]  Margot Gerritsen,et al.  Designing an efficient solution strategy for fluid flows. 1. A stable high order finite difference scheme and sharp shock resolution for the Euler equations , 1996 .

[34]  P. Sagaut,et al.  Subgrid-Scale Models for Large-Eddy Simulations of Compressible Wall Bounded Flows , 2000 .

[35]  P. Lax,et al.  Systems of conservation equations with a convex extension. , 1971, Proceedings of the National Academy of Sciences of the United States of America.

[36]  H. C. Yee,et al.  A class of high resolution explicit and implicit shock-capturing methods , 1989 .

[37]  D. Lilly,et al.  A proposed modification of the Germano subgrid‐scale closure method , 1992 .

[38]  S. Corrsin,et al.  Simple Eulerian time correlation of full-and narrow-band velocity signals in grid-generated, ‘isotropic’ turbulence , 1971, Journal of Fluid Mechanics.

[39]  P. Moin,et al.  A dynamic subgrid‐scale model for compressible turbulence and scalar transport , 1991 .

[40]  P. Moin,et al.  A numerical method for large-eddy simulation in complex geometries , 2004 .

[41]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[42]  S. Jaisankar,et al.  Diffusion regulation for Euler solvers , 2007, J. Comput. Phys..

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