equation_of_state.h

Go to the documentation of this file.

#pragma once
 
#include "cmath"
 
#include "adiabatic_index.h"
 
 
namespace exaparty {
  namespace kernels {
  namespace legacy {
    namespace eos {
 
      // @todo we need to make this modular. This is only valid for ideal gases. It needs to be exchangeable.
 
      __attribute__((always_inline, const)) inline static double gas_internal_energy_from_entropy(
        double density, double entropy
      ) {
 
        return entropy * adiabaticIndex::pow_gamma_minus_one(density) * hydro_one_over_gamma_minus_one;
      }
 
      __attribute__((always_inline, const)) inline static double gas_pressure_from_entropy(
        double density, double entropy
      ) {
 
        return entropy * adiabaticIndex::pow_gamma(density);
      }
 
      __attribute__((always_inline, const)) inline static double gas_entropy_from_pressure(
        double density, double pressure
      ) {
 
        return pressure * adiabaticIndex::pow_minus_gamma(density);
      }
 
      __attribute__((always_inline, const)) inline static double gas_soundspeed_from_entropy(
        double density, double entropy
      ) {
 
        return std::sqrt(hydro_gamma * adiabaticIndex::pow_gamma_minus_one(density) * entropy);
      }
 
      __attribute__((always_inline, const)) inline static double gas_entropy_from_internal_energy(
        double density, double u
      ) {
 
        return hydro_gamma_minus_one * u * adiabaticIndex::pow_minus_gamma_minus_one(density);
      }
 
      __attribute__((always_inline, const)) inline static double gas_pressure_from_internal_energy(
        double density, double u
      ) {
 
        return hydro_gamma_minus_one * u * density;
      }
 
      __attribute__((always_inline, const)) inline static double gas_internal_energy_from_pressure(
        double density, double pressure
      ) {
 
        return hydro_one_over_gamma_minus_one * pressure / density;
      }
 
      __attribute__((always_inline, const)) inline static double gas_soundspeed_from_internal_energy(double u) {
 
        return std::sqrt(u * hydro_gamma * hydro_gamma_minus_one);
      }
 
      __attribute__((always_inline, const)) inline static double gas_soundspeed_from_pressure(
        double density, double P
      ) {
 
        return std::sqrt(hydro_gamma * P / density);
      }
 
    } // namespace eos
  }   // namespace legacy
  }
} // namespace exaparty