IdealGasLaw.h

Go to the documentation of this file.

 
#pragma once
 
#include <cassert>
#include <string>
 
#include "MaterialLib/Fluid/FluidProperty.h"
#include "MaterialLib/PhysicalConstant.h"
 
namespace MaterialLib
{
namespace Fluid
{
class IdealGasLaw final : public FluidProperty
{
public:
    explicit IdealGasLaw(const double molar_mass) : _molar_mass(molar_mass) {}
 
    std::string getName() const override { return "Ideal gas law"; }
    double getValue(const ArrayType& var_vals) const override
    {
        return _molar_mass *
               var_vals[static_cast<int>(PropertyVariableType::p)] /
               (PhysicalConstant::IdealGasConstant *
                var_vals[static_cast<int>(PropertyVariableType::T)]);
    }
 
    double getdValue(const ArrayType& var_vals,
                     const PropertyVariableType var) const override
    {
        const double T = var_vals[static_cast<int>(PropertyVariableType::T)];
        const double p = var_vals[static_cast<int>(PropertyVariableType::p)];
        switch (var)
        {
            case PropertyVariableType::T:
                return dIdealGasLaw_dT(T, p);
            case PropertyVariableType::p:
                return dIdealGasLaw_dp(T, p);
            default:
                return 0.;
        }
    }
 
private:
    const double _molar_mass;
 
    double dIdealGasLaw_dT(const double T, const double pg) const
    {
        return -_molar_mass * pg / (PhysicalConstant::IdealGasConstant * T * T);
    }
 
    double dIdealGasLaw_dp(const double T, const double /* pg */) const
    {
        return _molar_mass / (PhysicalConstant::IdealGasConstant * T);
    }
};
}  // namespace Fluid
}  // namespace MaterialLib