IdealGasLawBinaryMixture.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "MaterialLib/MPL/Properties/IdealGasLawBinaryMixture.h"
 
#include "MaterialLib/PhysicalConstant.h"
 
namespace MaterialPropertyLib
{
IdealGasLawBinaryMixture::IdealGasLawBinaryMixture(std::string name)
{
    name_ = std::move(name);
}
 
PropertyDataType IdealGasLawBinaryMixture::value(
    VariableArray const& variable_array,
    ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/,
    double const /*dt*/) const
{
    const double R = MaterialLib::PhysicalConstant::IdealGasConstant;
    const double pGR = variable_array.gas_phase_pressure;
    const double T = variable_array.temperature;
    const double MG = variable_array.molar_mass;
 
    return pGR * MG / R / T;
}
 
PropertyDataType IdealGasLawBinaryMixture::dValue(
    VariableArray const& variable_array, Variable const variable,
    ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/,
    double const /*dt*/) const
{
    const double R = MaterialLib::PhysicalConstant::IdealGasConstant;
    const double pGR = variable_array.gas_phase_pressure;
    const double T = variable_array.temperature;
    const double MG = variable_array.molar_mass;
    const double dMG_dx = variable_array.molar_mass_derivative;
 
    switch (variable)
    {
        case Variable::gas_phase_pressure:
            return MG / R / T + pGR / R / T * dMG_dx;
        case Variable::capillary_pressure:
            return pGR / R / T * dMG_dx;
        case Variable::temperature:
            return pGR / R / T / T * (T * dMG_dx - MG);
        default:
            OGS_FATAL(
                "IdealGasLawBinaryMixture::dValue is implemented only for "
                "derivatives "
                "w.r.t. gas_pressure, capillary_pressure, and temperature.");
            break;
    }
    return 0.;
}
 
PropertyDataType IdealGasLawBinaryMixture::d2Value(
    VariableArray const& /*variable_array*/, Variable const /*variable1*/,
    Variable const /*variable2*/, ParameterLib::SpatialPosition const& /*pos*/,
    double const /*t*/, double const /*dt*/) const
{
    OGS_FATAL("IdealGasLawBinaryMixture::d2Value is not implemented.");
 
    return 0.;
}
 
}  // namespace MaterialPropertyLib