This repository contains information and code to reproduce the results presented in the article

@online{oblapenko2024entropyconservative,
  title={Entropy-conservative high-order methods for high-enthalpy gas flows},
  author={Oblapenko, Georgii and Torrilhon, Manuel},
  year={2024},
  eprint={2403.16882},
  archivePrefix={arXiv},
  primaryClass={physics.flu-dyn},
  eprintclass={physics.flu-dyn},
  journal={arXiv preprint 2403.16882},
  doi={10.48550/arXiv.2403.16882}
}

If you find these results useful, please cite the article mentioned above. If you use the implementations provided here, please also cite this repository as

@misc{oblapenko2024entropyconservativeRepro,
  title={Reproducibility repository for "Entropy-conservative high-order methods for high-enthalpy gas flows"},
  author={Oblapenko, Georgii and Torrilhon, Manuel},
  year={2024},
  howpublished={\url{https://github.com/knstmrd/paper-ec_trixi_inte}},
  doi={10.5281/zenodo.10869716}
}

Abstract

A framework for numerical evaluation of entropy-conservative volume fluxes in gas flows with internal energies is developed, for use with high-order discretization methods. The novelty of the approach lies in the ability to use arbitrary expressions for the internal degrees of freedom of the constituent gas species. The developed approach is implemented in an open-source discontinuous Galerkin code for solving hyperbolic equations. Numerical simulations are carried out for several model 2-D flows and the results are compared to those obtained with the finite volume-based solver DLR TAU.

Reproducing the results

Installation

To download the code using git, use

git clone git@github.com:knstmrd/paper-ec_trixi_inte.git

To instantiate the environment execute the following two commands:

cd paper-ec_trixi_inte
julia --project=. -e 'using Pkg; Pkg.instantiate()'

Note that the results are obtained using Julia 1.9.4. Thus, you might need to install the old Julia 1.9.4 release first and replace the julia calls from this README with /YOUR/PATH/TO/julia-1.9.4/bin/julia

Running the code

The scripts for the different test cases are located in the examples directory.

To execute them provide the path (if inside the paper-ec_trixi_inte directory):

julia --project=. examples/elixir_euler_vibr_aho_cylinder.jl

Authors

• Georgii Oblapenko (Corresponding Author)
• Manuel Torrilhon

Disclaimer

Everything is provided as is and without warranty. Use at your own risk!

Acknowledgments

This work has been supported by the German Research Foundation within the research unit DFG-FOR5409.