SF: An Open Source Object-Oriented Platform for Unstructured Shock-Fitting Methods

Recently, shock-fitting methods have been re-gaining popularity as a valid alternative to more standard shock-capturing schemes for tackling simulation of compressible flows including discontinuities, particularly on unstructured grids. This paper presents SF, the first open source framework for developing new-generation shock-fitting methods and couple them to potentially arbitrary Computational Fluid Dynamics solvers. The basic infrastructure and the completely automatized 2D version of the shock-fitting algorithm, which is entirely written in C++ using advanced object-oriented techniques, are described in this chapter, as well as its coupling to a state-of-the-art Residual Distribution solver. The latter is implemented within COOLFluiD, an open source platform for high-performance computing and multi-physics modeling. Numerical results related to the shock-fitted simulations of 2D inviscid (also including thermochemical nonequilibrium effects) and viscous hypersonic flows are discussed, demonstrating the current capabilities of SF.

[1]  Andrea Lani,et al.  Modelling of high-enthalpy, high-Mach number flows , 2009 .

[2]  Andrea Lani,et al.  An object oriented and high performance platform for aerothermodynamics simulation , 2008 .

[3]  William Gropp,et al.  Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries , 1997, SciTools.

[4]  Peter A. Gnoffo,et al.  Updates to Multi-Dimensional Flux Reconstruction for Hypersonic Simulations on Tetrahedral Grids , 2010 .

[5]  Andrea Lani,et al.  Conservative Residual Distribution Method for Viscous Double Cone Flows in Thermochemical Nonequilibrium , 2013 .

[6]  Andrea Lani,et al.  Modeling of non-equilibrium phenomena in expanding flows by means of a collisional-radiative model , 2013 .

[7]  Andrea Lani,et al.  Collisional radiative coarse-grain model for ionization in air , 2013 .

[8]  Peter A. Gnoffo,et al.  Computational Aerothermodynamic Simulation Issues on Unstructured Grids , 2004 .

[9]  Masashi Kanamori,et al.  Shock Wave Detection based on the Theory of Characteristics for CFD Results , 2011 .

[10]  Andrea Lani,et al.  A Residual Distribution Method for Symmetrized Systems in Thermochemical Nonequilibrium , 2011 .

[11]  Robert Haimes,et al.  Shock detection from computational fluid dynamics results , 1999 .

[12]  Hans G. Hornung,et al.  Non-equilibrium dissociating nitrogen flow over spheres and circular cylinders , 1972, Journal of Fluid Mechanics.

[13]  Renato Paciorri,et al.  Shock interaction computations on unstructured, two-dimensional grids using a shock-fitting technique , 2011, J. Comput. Phys..

[14]  H. Deconinck,et al.  An energy-dissipative remedy against carbuncle: Application to hypersonic flows around blunt bodies , 2016 .

[15]  Chul Park,et al.  Assessment of a two-temperature kinetic model for dissociating and weakly ionizing nitrogen , 1986 .

[16]  Renato Paciorri,et al.  Unsteady shock‐fitting for unstructured grids , 2016 .

[17]  J.M.A. Longo,et al.  Assessment of CFD capability for prediction of hypersonic shock interactions , 2012 .

[18]  Zi-Niu Wu,et al.  Review of shock wave detection method in CFD post-processing , 2013 .

[19]  Andrea Lani,et al.  A finite volume implicit time integration method for solving the equations of ideal magnetohydrodynamics for the hyperbolic divergence cleaning approach , 2011, J. Comput. Phys..

[20]  Renato Paciorri,et al.  An unstructured, three-dimensional, shock-fitting solver for hypersonic flows , 2013 .

[21]  Ralph Johnson,et al.  design patterns elements of reusable object oriented software , 2019 .

[22]  Andrea Lani,et al.  A fully-implicit finite-volume method for multi-fluid reactive and collisional magnetized plasmas on unstructured meshes , 2016, J. Comput. Phys..

[23]  Hang Si,et al.  TetGen: A quality tetrahedral mesh generator and a 3D Delaunay triangulator (Version 1.5 --- User's Manual) , 2013 .

[24]  Andrea Lani,et al.  Reusable object-oriented solutions for numerical simulation of PDEs in a high performance environment , 2006, Sci. Program..

[25]  J. Beveridge SELF-REGISTERING OBJECTS IN C++ , 1998 .

[26]  R. Paciorri,et al.  On the Role of Vibrational Excitation in Hypersonic Flow Computations , 1998 .

[27]  Renato Paciorri,et al.  Recognition of shock-wave patterns from shock-capturing solutions , 2012, CompIMAGE.

[28]  Andrea Lani,et al.  Multi-fluid Modeling of Magnetosonic Wave Propagation in the Solar Chromosphere: Effects of Impact Ionization and Radiative Recombination , 2016, 1611.08439.

[29]  Andrea Lani,et al.  COOLFluiD: An open computational platform for multi-physics simulation and research , 2013 .

[30]  Mikhail S. Ivanov,et al.  Computation of weak steady shock reflections by means of an unstructured shock-fitting solver , 2010 .

[31]  Matthew G. Knepley,et al.  petsc: Portable, Extensible Toolkit for Scientific Computation , 2016 .

[32]  H. Deconinck,et al.  Assessment of heat flux prediction capabilities of residual distribution method: Application to atmospheric entry problems , 2015 .

[33]  Andrea Lani,et al.  A GPU-enabled Finite Volume solver for global magnetospheric simulations on unstructured grids , 2014, Comput. Phys. Commun..

[34]  Andrea Lani,et al.  Analysis of non-equilibrium phenomena in inductively coupled plasma generators , 2016 .

[35]  Gerald J. Agin Fitting Ellipses and General Second-Order Curves. , 1981 .

[36]  Renato Paciorri,et al.  A shock-fitting technique for 2D unstructured grids , 2009 .

[37]  Andrea Lani,et al.  The COOLFluiD Framework: Design Solutions for High Performance Object Oriented Scientific Computing Software , 2005, International Conference on Computational Science.

[38]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[39]  Andrea Lani,et al.  An object-oriented implementation of a parallel Monte Carlo code for radiation transport , 2016, Comput. Phys. Commun..

[40]  Andrea Lani,et al.  Numerical investigation of the non equilibrium shock-layer around the EXPERT vehicle , 2007 .