Saturation.cpp

Go to the documentation of this file.

 
#include "Saturation.h"
 
namespace ProcessLib::TH2M
{
namespace ConstitutiveRelations
{
void SaturationModel::eval(SpaceTimeData const& x_t,
                           MediaData const& media_data,
                           CapillaryPressureData const& p_cap,
                           SaturationData& S_L_data,
                           SaturationDataDeriv& dS_L_data) const
{
    namespace MPL = MaterialPropertyLib;
    MPL::VariableArray variables;
    variables.capillary_pressure = p_cap();
 
    auto const& medium = media_data.medium;
 
    S_L_data.S_L = medium.property(MPL::PropertyType::saturation)
                       .template value<double>(variables, x_t.x, x_t.t, x_t.dt);
 
    dS_L_data() = medium.property(MPL::PropertyType::saturation)
                      .template dValue<double>(
                          variables, MPL::Variable::capillary_pressure, x_t.x,
                          x_t.t, x_t.dt);
}
}  // namespace ConstitutiveRelations
}  // namespace ProcessLib::TH2M