sedov_computer.hpp File Reference

This class produces 1d solutions for a sedov blast wave propagating through a density gradient: rho = rho**(-omega) , in planar(1D), cylindrical(2D) or spherical geometry(3D) for the 'standard', 'singular' and 'vaccum' cases. More...

#include <cmath>
#include <fstream>
#include <functional>
#include <iomanip>
#include <iostream>
#include <string>
#include <variant>
#include <vector>
#include <cstdint>

Go to the source code of this file.

Classes
class	SedovComputer

Typedefs
using	Real = double
using	KeyType = uint64_t

Functions
template<class... T, class... Separators>
void	writeAscii (size_t firstIndex, size_t lastIndex, const std::string &path, bool append, const std::vector< std::variant< T *... >> &fields, Separators &&... separators)
void	printHelp (char *binName)
void	writeColumns1D (const std::string &path)

Detailed Description

This class produces 1d solutions for a sedov blast wave propagating through a density gradient: rho = rho**(-omega) , in planar(1D), cylindrical(2D) or spherical geometry(3D) for the 'standard', 'singular' and 'vaccum' cases.

Author

Jose A. Escartin ja.es.nosp@m.cart.nosp@m.in@gm.nosp@m.ail..nosp@m.com

• standard case: a nonzero solution extends from the shock to the origin, where the pressure is finite.
• singular case: a nonzero solution extends from the shock to the origin, where the pressure vanishes.
• vacuum case : a nonzero solution extends from the shock to a boundary point, where the density vanishes making the pressure meaningless.

     This routine is a C++ conversion of one Fortran code based in these two papers:
     - "Evaluation of the sedov-von neumann-taylor blast wave solution", Jim Kamm, la-ur-00-6055
     - "The sedov self-similiar point blast solutions in nonuniform media", David Book, shock waves, 4, 1, 1994
  
     Although the ordinary differential equations are analytic, the sedov expressions appear to become singular for
  

  various combinations of parameters and at the lower limits of the integration range. All these singularies are removable and done so by this routine.

This file is based on the analytical solution presented in SPH-EXA https://github.com/unibas-dmi-hpc/SPH-EXA

Date

2023-11-30