LiquidDensity.cpp

Go to the documentation of this file.

 
#include "LiquidDensity.h"
 
namespace ProcessLib::ThermoRichardsMechanics
{
template <int DisplacementDim>
void LiquidDensityModel<DisplacementDim>::eval(
    SpaceTimeData const& x_t, MediaData const& media_data,
    CapillaryPressureData<DisplacementDim> const& p_cap_data,
    TemperatureData<DisplacementDim> const& T_data,
    LiquidDensityData& out) const
{
    namespace MPL = MaterialPropertyLib;
    MPL::VariableArray variables;
    variables.liquid_phase_pressure = -p_cap_data.p_cap;
    variables.temperature = T_data.T;
 
    auto const& liquid_phase = media_data.liquid;
 
    out.rho_LR = liquid_phase.property(MPL::PropertyType::density)
                     .template value<double>(variables, x_t.x, x_t.t, x_t.dt);
 
    out.drho_LR_dp = liquid_phase.property(MPL::PropertyType::density)
                         .template dValue<double>(
                             variables, MPL::Variable::liquid_phase_pressure,
                             x_t.x, x_t.t, x_t.dt);
 
    out.drho_LR_dT =
        liquid_phase.property(MPL::PropertyType::density)
            .template dValue<double>(variables, MPL::Variable::temperature,
                                     x_t.x, x_t.t, x_t.dt);
}
 
template struct LiquidDensityModel<2>;
template struct LiquidDensityModel<3>;
}  // namespace ProcessLib::ThermoRichardsMechanics