Techniques for improving monotonicity in a fourth-order finite-volume algorithm solving shocks and detonations

Abstract Techniques are proposed to reduce numerical oscillations in a fourth-order, finite-volume algorithm for solving thermally-perfect, reacting fluid flows with strong discontinuities, such as shock or detonation waves. These additional mechanisms have proven necessary for multispecies flows solved at fourth-order accuracy, and contribute towards bounding the variation of the solution in the vicinity of strong discontinuities. There, oscillations can form due to strong gradients in the flow and may be further intensified by numerical procedures introduced to treat the thermally-perfect thermodynamic system and physical constraints on species mass fractions. The new techniques are designed to respect the conservative property of the base algorithm, retain fourth-order accuracy of the solution in regions of smooth flow, and cooperate with the high-order piecewise parabolic method limiter. Extensive numerical tests, ranging from multispecies mixing flows to reacting flows with detonations, are performed to verify that the new techniques meet the design criteria while effectively suppressing oscillations. The proposed techniques are applied to solve the Shu-Osher and double Mach reflection problems and a set of oblique detonation wave problems. The results demonstrate the effectiveness and robustness of the algorithm.

[1]  Landon D. Owen,et al.  A high-order adaptive algorithm for multispecies gaseous flows on mapped domains , 2018, Computers & Fluids.

[2]  Sanford Gordon,et al.  Computer program for calculation of complex chemical equilibrium compositions , 1972 .

[3]  Kyu Hong Kim,et al.  Accurate, efficient and monotonic numerical methods for multi-dimensional compressible flows Part II: Multi-dimensional limiting process , 2005 .

[4]  Anthony Ruiz,et al.  Large-Eddy Simulation of Supercritical-Pressure Round Jets , 2010 .

[5]  Chongam Kim,et al.  Multi-dimensional limiting process for finite volume methods on unstructured grids , 2012 .

[6]  R. J. R. Williams,et al.  An improved reconstruction method for compressible flows with low Mach number features , 2008, J. Comput. Phys..

[7]  Phillip Colella,et al.  A limiter for PPM that preserves accuracy at smooth extrema , 2008, J. Comput. Phys..

[8]  Landon D. Owen,et al.  A parallel adaptive numerical method with generalized curvilinear coordinate transformation for compressible Navier–Stokes equations , 2016 .

[9]  Hiroshi Terashima,et al.  Approach for simulating gas-liquid-like flows under supercritical pressures using a high-order central differencing scheme , 2012, J. Comput. Phys..

[10]  J. M. Powers,et al.  Mapped weighted essentially non-oscillatory schemes: Achieving optimal order near critical points , 2005 .

[11]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[12]  Phillip Colella,et al.  High-order, finite-volume methods in mapped coordinates , 2010, J. Comput. Phys..

[13]  P. Roache.  Fundamentals of computational fluid dynamics , 1998 .

[14]  Elaine S. Oran,et al.  Detonation structures behind oblique shocks , 1994 .

[15]  Mark F. Adams,et al.  Chombo Software Package for AMR Applications Design Document , 2014 .

[16]  Thomas H. Pulliam,et al.  Fundamental Algorithms in Computational Fluid Dynamics , 2014 .

[17]  Rémi Abgrall,et al.  An adaptive shock-capturing algorithm for solving unsteady reactive flows , 2003 .

[18]  Chongam Kim,et al.  Multi-dimensional limiting process for three-dimensional flow physics analyses , 2008, J. Comput. Phys..

[19]  G. Billet.  Improvement of convective concentration fluxes in a one step reactive flow solver , 2005 .

[20]  Xinfeng Gao,et al.  A Fourth-Order Finite-Volume Algorithm for Compressible Flow with Chemical Reactions on Mapped Grids , 2017 .

[21]  Béatrice Rivière,et al.  Discontinuous Galerkin methods for solving elliptic and parabolic equations - theory and implementation , 2008, Frontiers in applied mathematics.

[22]  P. Colella,et al.  A Single Stage Flux-Corrected Transport Algorithm for High-Order Finite-Volume Methods , 2015, 1506.02999.

[23]  P. Woodward.  Piecewise-parabolic methods for astrophysical fluid dynamics , 1986 .

[24]  Dong-Ho Lee,et al.  A new approach of a limiting process for multi-dimensional flows , 2010, J. Comput. Phys..

[25]  Richard P. Brent,et al.  An Algorithm with Guaranteed Convergence for Finding a Zero of a Function , 1971, Comput. J..

[26]  P. Colella.  Multidimensional upwind methods for hyperbolic conservation laws , 1990 .

[27]  Ch. Hirsch,et al.  Fundamentals Of Computational Fluid Dynamics , 2016 .

[28]  Phillip Colella,et al.  A Freestream-Preserving High-Order Finite-Volume Method for Mapped Grids with Adaptive-Mesh Refinement , 2011 .

[29]  Smadar Karni,et al.  Hybrid Multifluid Algorithms , 1996, SIAM J. Sci. Comput..

[30]  Phillip Colella,et al.  A HIGH-ORDER FINITE-VOLUME METHOD FOR CONSERVATION LAWS ON LOCALLY REFINED GRIDS , 2011 .

[31]  Z J Wang,et al.  High-order computational fluid dynamics tools for aircraft design , 2014, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[32]  Rémi Abgrall,et al.  High‐order CFD methods: current status and perspective , 2013 .

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

[34]  Tariq D. Aslam,et al.  Simulations of pulsating one-dimensional detonations with true fifth order accuracy , 2006, J. Comput. Phys..

[35]  Xinfeng Gao,et al.  A high-performance finite-volume algorithm for solving partial differential equations governing compressible viscous flows on structured grids , 2016, Comput. Math. Appl..

[36]  H. Teng,et al.  Numerical study of oblique detonation wave initiation in a stoichiometric hydrogen-air mixture , 2015 .

[37]  Landon D. Owen,et al.  A freestream-preserving fourth-order finite-volume method in mapped coordinates with adaptive-mesh refinement , 2015 .

[38]  J. Lambert.  Numerical Methods for Ordinary Differential Systems: The Initial Value Problem , 1991 .

[39]  Robert J. Kee,et al.  CHEMKIN-III: A FORTRAN chemical kinetics package for the analysis of gas-phase chemical and plasma kinetics , 1996 .

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

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

[42]  Kenneth K. Kuo,et al.  A low-dissipation and time-accurate method for compressible multi-component flow with variable specific heat ratios , 2011, J. Comput. Phys..

[43]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

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

[46]  Rémi Abgrall,et al.  Computations of compressible multifluids , 2001 .

[47]  H. Chelliah,et al.  Numerical prediction of oblique detonation wave structures using detailed and reduced reaction mechanisms , 1997 .

[48]  Matthias Ihme,et al.  An entropy-stable hybrid scheme for simulations of transcritical real-fluid flows , 2017, J. Comput. Phys..

[49]  M. Z. Dauhoo,et al.  A new fourth-order non-oscillatory central scheme for hyperbolic conservation laws , 2008 .

[50]  Luı́s Fernando Figueira da Silva,et al.  Stabilization of an oblique detonation wave by a wedge: a parametric numerical study , 2000 .