A stable high-order finite difference scheme for the compressible Navier-Stokes equations, far-field boundary conditions

We construct a stable high-order finite difference scheme for the compressible Navier-Stokes equations, that satisfy an energy estimate. The equations are discretized with high-order accurate finite difference methods that satisfy a Summation-By-Parts rule. The boundary conditions are imposed with penalty terms known as the Simultaneous Approximation Term technique. The main result is a stability proof for the full three-dimensional Navier-Stokes equations, including the boundary conditions. We show the theoretical third-, fourth-, and fifth-order convergence rate, for a viscous shock, where the analytic solution is known. We demonstrate the stability and discuss the non-reflecting properties of the outflow conditions for a vortex in free space. Furthermore, we compute the three-dimensional vortex shedding behind a circular cylinder in an oblique free stream for Mach number 0.5 and Reynolds number 500.

[1]  Magnus Svärd,et al.  On the order of accuracy for difference approximations of initial-boundary value problems , 2006, J. Comput. Phys..

[2]  Wendy Kress,et al.  High Order Finite Difference Methods in Space and Time , 2003 .

[3]  Haecheon Choi,et al.  Control of laminar vortex shedding behind a circular cylinder using splitter plates , 1996 .

[4]  John C. Strikwerda,et al.  Initial boundary value problems for incompletely parabolic systems , 1976 .

[5]  A. Majda,et al.  Radiation boundary conditions for acoustic and elastic wave calculations , 1979 .

[6]  Jan Nordström,et al.  The use of characteristic boundary conditions for the Navier-Stokes equations , 1995 .

[7]  Magnus Svärd,et al.  High-order accurate computations for unsteady aerodynamics , 2007 .

[8]  R. D. Richtmyer,et al.  Survey of the stability of linear finite difference equations , 1956 .

[9]  J. Nordström,et al.  Summation by Parts Operators for Finite Difference Approximations of Second-Derivatives with Variable Coefficients , 2004, Journal of Scientific Computing.

[10]  Parviz Moin,et al.  B-Spline Method and Zonal Grids for Simulations of Complex Turbulent Flows , 1997 .

[11]  Magnus Svärd,et al.  Accuracy requirements for transient aerodynamics , 2003 .

[12]  B. Strand.  Summation by parts for finite difference approximations for d/dx , 1994 .

[13]  David Gottlieb,et al.  Optimal time splitting for two- and three-dimensional navier-stokes equations with mixed derivatives , 1981 .

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

[15]  H. Kreiss,et al.  Time-Dependent Problems and Difference Methods , 1996 .

[16]  Magnus Svärd,et al.  Stable and Accurate Artificial Dissipation , 2004, J. Sci. Comput..

[17]  D. Gottlieb,et al.  A Stable and Conservative Interface Treatment of Arbitrary Spatial Accuracy , 1999 .

[18]  Charles R. Johnson,et al.  Matrix analysis , 1985, Statistical Inference for Engineers and Data Scientists.

[19]  Magnus Svärd,et al.  Well-Posed Boundary Conditions for the Navier-Stokes Equations , 2005, SIAM J. Numer. Anal..

[20]  Jan Nordström,et al.  High-order finite difference methods, multidimensional linear problems, and curvilinear coordinates , 2001 .

[21]  Magnus Svärd,et al.  On Coordinate Transformations for Summation-by-Parts Operators , 2004, J. Sci. Comput..

[22]  H. Kreiss,et al.  Finite Element and Finite Difference Methods for Hyperbolic Partial Differential Equations , 1974 .

[23]  Gilbert Strang,et al.  Accurate partial difference methods , 1964 .

[24]  Bertil Gustafsson,et al.  Incompletely parabolic problems in fluid dynamics , 1978 .

[25]  C. Williamson.  Defining a Universal and Continuous Strouhal-Reynolds Number Relationship for the Laminar Vortex She , 1988 .

[26]  T. Pulliam,et al.  A diagonal form of an implicit approximate-factorization algorithm , 1981 .

[27]  D. Gottlieb,et al.  Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes , 1994 .

[28]  Magnus Svärd,et al.  Steady-State Computations Using Summation-by-Parts Operators , 2005, J. Sci. Comput..

[29]  Jan Nordström,et al.  Boundary and Interface Conditions for High-Order Finite-Difference Methods Applied to the Euler and Navier-Stokes Equations , 1999 .

[30]  H. Kreiss,et al.  Comparison of accurate methods for the integration of hyperbolic equations , 1972 .

[31]  Jan S. Hesthaven,et al.  A Stable Penalty Method for the Compressible Navier-Stokes Equations: I. Open Boundary Conditions , 1996, SIAM J. Sci. Comput..

[32]  H. Kreiss,et al.  Initial-Boundary Value Problems and the Navier-Stokes Equations , 2004 .