PengRobinson.h

Go to the documentation of this file.

 
#pragma once
 
#include "MaterialLib/MPL/Property.h"
#include "MaterialLib/MPL/VariableType.h"
#include "MaterialLib/PhysicalConstant.h"
#include "ParameterLib/Parameter.h"
 
namespace MaterialPropertyLib
{
class PengRobinson final : public Property
{
public:
    explicit PengRobinson(const double Tc, const double pc, const double omega)
        : Tc_(Tc), pc_(pc), omega_(omega)
    {
        const double gas_constant =
            MaterialLib::PhysicalConstant::IdealGasConstant;
        a_ = 0.457235 * gas_constant * gas_constant * Tc_ * Tc_ / pc_;
        b_ = 0.077796 * gas_constant * Tc_ / pc_;
    }
 
    PropertyDataType value(
        MaterialPropertyLib::VariableArray const& variable_array,
        ParameterLib::SpatialPosition const& pos, double const t,
        double const dt) const override;
    PropertyDataType dValue(
        MaterialPropertyLib::VariableArray const& variable_array,
        Variable const variable, ParameterLib::SpatialPosition const& pos,
        double const t, double const dt) const override;
 
private:
    const double Tc_;
    const double pc_;
    const double omega_;
    double a_;
    double b_;
};
}  // namespace MaterialPropertyLib