LiquidViscosity.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "LiquidViscosity.h"
 
namespace ProcessLib::ThermoRichardsMechanics
{
template <int DisplacementDim>
void LiquidViscosityModel<DisplacementDim>::eval(
    SpaceTimeData const& x_t, MediaData const& media_data,
    LiquidDensityData const& rho_L_data,
    TemperatureData<DisplacementDim> const& T_data,
    LiquidViscosityData& out) const
{
    namespace MPL = MaterialPropertyLib;
    MPL::VariableArray variables;
    variables.temperature = T_data.T;
    variables.density = rho_L_data.rho_LR;
 
    *out = media_data.liquid.property(MPL::PropertyType::viscosity)
               .template value<double>(variables, x_t.x, x_t.t, x_t.dt);
}
 
template struct LiquidViscosityModel<2>;
template struct LiquidViscosityModel<3>;
}  // namespace ProcessLib::ThermoRichardsMechanics