SolidDensity.cpp

Go to the documentation of this file.

 
#include "SolidDensity.h"
 
namespace ProcessLib::ThermoRichardsMechanics
{
template <int DisplacementDim>
void SolidDensityModel<DisplacementDim>::eval(
    SpaceTimeData const& x_t,
    MediaData const& media_data,
    PorosityData const& poro_data,
    TemperatureData<DisplacementDim> const& T_data,
    SolidDensityData& out) const
{
    namespace MPL = MaterialPropertyLib;
    MPL::VariableArray variables;
    variables.temperature = T_data.T;
 
    out.rho_SR = media_data.solid.property(MPL::PropertyType::density)
                     .template value<double>(variables, x_t.x, x_t.t, x_t.dt);
    out.dry_density_solid =
        (1 - poro_data.phi) * out.rho_SR;  // TODO only for output
}
 
template struct SolidDensityModel<2>;
template struct SolidDensityModel<3>;
}  // namespace ProcessLib::ThermoRichardsMechanics