simple_waves_ideal_gas.hpp

Go to the documentation of this file.

#ifndef SIMPLE_WAVES_IDEAL_GAS_HPP
#define SIMPLE_WAVES_IDEAL_GAS_HPP 1

#include <string>
#include "../one_dimensional/spatial_distribution1d.hpp"
#include "../common/equation_of_state.hpp"
#include "../common/ideal_gas.hpp"

//using namespace std;
using std::string;

double calc_entropy(double d, double p, double g);

class ConstEntropy: public SpatialDistribution1D
{
public:

  ConstEntropy(SpatialDistribution1D const& density,
               double s, double g);

  double operator()(double x) const;

private:

  SpatialDistribution1D const& density_;
  const double s_;
  const double g_;
};

class SoundSpeedDist: public SpatialDistribution1D
{
public:

  SoundSpeedDist(SpatialDistribution1D const& pressure,
                 SpatialDistribution1D const& density,
                 EquationOfState const& eos);

  double operator()(double x) const;

private:

  SpatialDistribution1D const& pressure_;
  SpatialDistribution1D const& density_;
  EquationOfState const& eos_;
};

double calc_riemann_invariant(double v, double c,
                              double g, bool dir);

class ConstRiemannInv: public SpatialDistribution1D
{
public:

  ConstRiemannInv(double rv, bool dir,
                  SpatialDistribution1D const& sound_speed,
                  double g);

  double operator()(double x) const;

private:

  const double rv_;
  const bool dir_;
  SpatialDistribution1D const& sound_speed_;
  const double g_;
};

class SimpleWaveIdealGasInitCond
{
public:

  SimpleWaveIdealGasInitCond
  (SpatialDistribution1D const& density,
   double entropy,
   double adiabatic_index,
   double edge=3);

  IdealGas const& getEOS(void) const;

  SpatialDistribution1D const& getProfile(string pname) const;

private:

  SpatialDistribution1D const& density_;
  ConstEntropy pressure_;
  IdealGas eos_;
  SoundSpeedDist sound_speed_;
  double c_edge_;
  double rv_edge_;
  ConstRiemannInv xvelocity_;
  Uniform yvelocity_;
};

class EntropyProf: public SpatialDistribution1D
{
public:

  EntropyProf(SpatialDistribution1D const& density,
              SpatialDistribution1D const& pressure,
              double adiabatic_index);

  double operator()(double x) const;

private:

  SpatialDistribution1D const& density_;
  SpatialDistribution1D const& pressure_;
  double g_;
};

#endif // SIMPLE_WAVES_IDEAL_GAS_HPP