About the formulation, verification and validation of the hypersonic flow solver Eilmer

We describe the formulation of the gas dynamics and high-temperature thermochemical modules of the Eilmer code, an open-source Navier-Stokes solver for transient compressible flow in two and three dimensions. The core gas dynamics formulation is based on finite-volume cells, and the thermochemical effects are handled with specialised updating schemes that are coupled into the overall time-stepping scheme. Verification of the code is explored via a number of case studies that use analytic and semi-analytic solutions as comparison. These include both smooth and shocked flows and are used to demonstrate the order of spatial accuracy of the code. Cases include manufactured solutions for rather abstract inviscid and viscous flow, an idealised detonation wave supported by a curved body, and the transient flow of an idealised but high-performance shock tube. Validation of the inviscid gas dynamics and thermochemical models is then explored using data from a selection of experimental studies. These studies include ballistic range experiments with chemically-inert noble gases and high-temperature chemically-reacting air. These comparisons show that the code performs well and they provide a lesson in considering a range of experimental data rather than relying upon isolated data points for validation. These verification and validation cases are described in full detail and will be useful for other code developers of high-temperature compressible flow solvers. Copyright (C) 2013 John Wiley & Sons, Ltd.

[1]  Richard A. Thompson,et al.  A review of reaction rates and thermodynamic and transport properties for the 11-species air model for chemical and thermal nonequilibrium calculations to 30000 K , 1989 .

[2]  R. Kenneth Lobb,et al.  Chapter 26 - Experimental Measurement of Shock Detachment Distance on Spheres Fired in Air at Hypervelocities , 1964 .

[3]  J. Boris,et al.  A numerical technique for solving stiff ordinary differential equations associated with the chemical kinetics of reactive-flow problems , 1977 .

[4]  K. Sawada,et al.  Validation of hypersonic chemical equilibrium flow calculations using ballistic-range data , 2001 .

[5]  Christopher J. Roy,et al.  Verification of Euler/Navier–Stokes codes using the method of manufactured solutions , 2004 .

[6]  F. S. Billig,et al.  Shock-wave shapes around spherical-and cylindrical-nosed bodies. , 1967 .

[7]  P. Roache.  Verification of Codes and Calculations , 1998 .

[8]  R. Schwartz,et al.  Shock Location in Front of a Sphere as a Measure of Real Gas Effects , 1956 .

[9]  G. D. van Albada,et al.  A comparative study of computational methods in cosmic gas dynamics , 1982 .

[10]  Tariq D. Aslam,et al.  Exact Solution for Multidimensional Compressible Reactive Flow for Verifying Numerical Algorithms , 2006 .

[11]  Patrick J. Roache,et al.  Symbolic manipulation and computational fluid dynamics , 1985 .

[12]  Alan H. Karp,et al.  A variable order Runge-Kutta method for initial value problems with rapidly varying right-hand sides , 1990, TOMS.

[13]  J. M. Powers,et al.  Comparison of numerical oblique detonation solutions with an asymptotic benchmark , 1992 .

[14]  Rowan J. Gollan,et al.  Computational Modelling of High-Temperature Gas Effects with Application to Hypersonic Flows , 2008 .

[15]  Peter A. Jacobs,et al.  Rarefied, superorbital flows in an expansion tube , 2004 .

[16]  Peter A. Jacobs,et al.  Single-Block Navier-Stokes Integrator , 1991 .

[17]  Peter A. Jacobs,et al.  Numerical simulation of transient hypervelocity flow in an expansion tube , 1994 .

[18]  Meng-Sing Liou,et al.  A Flux Splitting Scheme with High-Resolution and Robustness for Discontinuities(Proceedings of the 12th NAL Symposium on Aircraft Computational Aerodynamics) , 1994 .

[19]  C. J. Roy.  Grid Convergence Error Analysis for Mixed-Order Numerical Schemes , 2001 .

[20]  P. Roache.  Code Verification by the Method of Manufactured Solutions , 2002 .

[21]  M. Liou,et al.  A New Flux Splitting Scheme , 1993 .

[22]  D. Reda,et al.  Correlation of nosetip boundary-layer transition data measured in ballistics-range experiments , 1979 .

[23]  Peter A. Jacobs,et al.  A simulation technique for radiating shock tube flows , 2009 .

[24]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .

[25]  Michael N. Macrossan,et al.  The equilibrium flux method for the calculation of flows with non-equilibrium chemical reactions , 1989 .

[26]  Richard J. Goozee,et al.  Simulation of a complete reflected shock tunnel showing a vortex mechanism for flow contamination , 2006 .

[27]  F. Homann.  Der Einfluß großer Zähigkeit bei der Strömung um den Zylinder und um die Kugel , 1936 .

[28]  Daniel Potter,et al.  Modelling of radiating shock layers for atmospheric entry at Earth and Mars , 2011 .

[29]  Joseph M. Powers,et al.  Approximate solutions for oblique detonations in the hypersonic limit , 1992 .