ClausiusClapeyron.cpp

Go to the documentation of this file.

 
#include "MaterialLib/MPL/Properties/ClausiusClapeyron.h"
 
#include "MaterialLib/MPL/Medium.h"
#include "MaterialLib/PhysicalConstant.h"
 
namespace MaterialPropertyLib
{
ClausiusClapeyron::ClausiusClapeyron(std::string name,
                                     const double triple_temperature,
                                     const double triple_pressure,
                                     const double critical_temperature,
                                     const double critical_pressure,
                                     const double ref_temperature,
                                     const double ref_pressure)
    : T_triple_(triple_temperature),
      p_triple_(triple_pressure),
      T_critical_(critical_temperature),
      p_critical_(critical_pressure),
      T_ref_(ref_temperature),
      p_ref_(ref_pressure)
{
    name_ = std::move(name);
}
 
void ClausiusClapeyron::checkScale() const
{
    if (!(std::holds_alternative<Phase*>(scale_) ||
          std::holds_alternative<Component*>(scale_)))
    {
        OGS_FATAL(
            "The property 'ClausiusClapeyron' is implemented on 'phase' and "
            "'component' scales only.");
    }
}
 
double ClausiusClapeyron::molarMass(
    std::variant<Medium*, Phase*, Component*> const scale,
    VariableArray const& variable_array,
    ParameterLib::SpatialPosition const& pos, double const t,
    double const dt) const
{
    return std::visit(
        [&variable_array, &pos, t, dt](auto&& s) -> double
        {
            return s->property(PropertyType::molar_mass)
                .template value<double>(variable_array, pos, t, dt);
        },
        scale);
}
 
PropertyDataType ClausiusClapeyron::value(
    VariableArray const& variable_array,
    ParameterLib::SpatialPosition const& pos, double const t,
    double const dt) const
{
    const double T = variable_array.temperature;
 
    const double M = molarMass(scale_, variable_array, pos, t, dt);
 
    if (T >= T_critical_)
    {
        return p_critical_;
    }
    if (T <= T_triple_)
    {
        return p_triple_;
    }
 
    const double dh = variable_array.enthalpy_of_evaporation;
    const double R = MaterialLib::PhysicalConstant::IdealGasConstant;
 
    return p_ref_ * std::exp((1. / T_ref_ - 1. / T) * M * dh / R);
}
 
PropertyDataType ClausiusClapeyron::dValue(
    VariableArray const& variable_array, Variable const variable,
    ParameterLib::SpatialPosition const& pos, double const t,
    double const dt) const
{
    const double T = variable_array.temperature;
    const double M = molarMass(scale_, variable_array, pos, t, dt);
 
    if (T > T_critical_)
    {
        return 0.;
    }
    if (T < T_triple_)
    {
        return 0.;
    }
    if (variable == Variable::gas_phase_pressure)
    {
        return 0.;
    }
 
    const double R = MaterialLib::PhysicalConstant::IdealGasConstant;
    const double dh = variable_array.enthalpy_of_evaporation;
    const double p_vap = std::get<double>(value(variable_array, pos, t, dt));
 
    if (variable == Variable::temperature)
    {
        return p_vap * M * dh / (R * T * T);
    }
    OGS_FATAL(
        "ClausiusClapeyron::dValue is implemented for derivatives with respect "
        "to phase pressure and temperature only.");
}
 
PropertyDataType ClausiusClapeyron::d2Value(
    VariableArray const& /*variable_array*/, Variable const /*variable1*/,
    Variable const /*variable2*/, ParameterLib::SpatialPosition const& /*pos*/,
    double const /*t*/, double const /*dt*/) const
{
    OGS_FATAL("ClausiusClapeyron::d2Value is not implemented.");
}
 
}  // namespace MaterialPropertyLib