OGS
List of all members | Public Member Functions
ProcessLib::ThermoRichardsMechanics::TRMHeatStorageAndFluxModel< DisplacementDim > Struct Template Reference

Detailed Description

template<int DisplacementDim>
struct ProcessLib::ThermoRichardsMechanics::TRMHeatStorageAndFluxModel< DisplacementDim >

Definition at line 33 of file TRMHeatStorageAndFlux.h.

#include <TRMHeatStorageAndFlux.h>

Public Member Functions

void eval (SpaceTimeData const &x_t, MediaData const &media_data, LiquidDensityData const &rho_L_data, SolidDensityData const &rho_S_data, SaturationData const &S_L_data, SaturationDataDeriv const &dS_L_data, PorosityData const &poro_data, LiquidViscosityData const &mu_L_data, PermeabilityData< DisplacementDim > const &perm, TemperatureData< DisplacementDim > const &T_data, DarcyLawData< DisplacementDim > const &darcy_data, TRMHeatStorageAndFluxData< DisplacementDim > &out) const
 

Member Function Documentation

◆ eval()

template<int DisplacementDim>
void ProcessLib::ThermoRichardsMechanics::TRMHeatStorageAndFluxModel< DisplacementDim >::eval ( SpaceTimeData const & x_t,
MediaData const & media_data,
LiquidDensityData const & rho_L_data,
SolidDensityData const & rho_S_data,
SaturationData const & S_L_data,
SaturationDataDeriv const & dS_L_data,
PorosityData const & poro_data,
LiquidViscosityData const & mu_L_data,
PermeabilityData< DisplacementDim > const & perm,
TemperatureData< DisplacementDim > const & T_data,
DarcyLawData< DisplacementDim > const & darcy_data,
TRMHeatStorageAndFluxData< DisplacementDim > & out ) const

Definition at line 18 of file TRMHeatStorageAndFlux.cpp.

27{
28 namespace MPL = MaterialPropertyLib;
29 MPL::VariableArray variables;
30 variables.temperature = T_data.T;
31 variables.porosity = poro_data.phi;
32 variables.liquid_saturation = S_L_data.S_L;
33
34 auto const& liquid_phase = media_data.liquid;
35 auto const& solid_phase = media_data.solid;
36
37 auto const specific_heat_capacity_fluid =
38 liquid_phase.property(MaterialPropertyLib::specific_heat_capacity)
39 .template value<double>(variables, x_t.x, x_t.t, x_t.dt);
40
41 auto const specific_heat_capacity_solid =
42 solid_phase
43 .property(MaterialPropertyLib::PropertyType::specific_heat_capacity)
44 .template value<double>(variables, x_t.x, x_t.t, x_t.dt);
45
46 double const phi = poro_data.phi;
47
48 // TODO real vs. "non-real" values
49 double const volumetric_heat_capacity_liquid =
50 rho_L_data.rho_LR * specific_heat_capacity_fluid;
51
52 // NOTE: Gas phase is not included.
53 double const volumetric_heat_capacity_liquid_and_solid =
54 rho_S_data.rho_SR * specific_heat_capacity_solid * (1 - phi) +
55 S_L_data.S_L * phi * volumetric_heat_capacity_liquid;
56
57 out.M_TT_X_NTN = volumetric_heat_capacity_liquid_and_solid;
58
59 out.K_TT_Laplace = MaterialPropertyLib::formEigenTensor<DisplacementDim>(
60 media_data.medium
61 .property(MaterialPropertyLib::PropertyType::thermal_conductivity)
62 .value(variables, x_t.x, x_t.t, x_t.dt));
63
64 // Unit is J / m^2 / s / K. It's not a heat flux, but related.
65 out.advective_heat_flux_contribution_to_K_liquid =
66 volumetric_heat_capacity_liquid * darcy_data();
67
68 //
69 // temperature equation, pressure part
70 //
71 out.K_Tp_NT_V_dN = -volumetric_heat_capacity_liquid * perm.k_rel /
72 mu_L_data() * (perm.Ki.transpose() * T_data.grad_T);
73 out.K_Tp_X_NTN = -volumetric_heat_capacity_liquid *
74 darcy_data().dot(T_data.grad_T) / perm.k_rel *
75 perm.dk_rel_dS_L * dS_L_data.dS_L_dp_cap;
76}
MaterialPropertyLib::formEigenTensor
Eigen::Matrix< double, GlobalDim, GlobalDim > formEigenTensor(MaterialPropertyLib::PropertyDataType const &values)
Definition FormEigenTensor.cpp:115

References ProcessLib::ThermoRichardsMechanics::TRMHeatStorageAndFluxData< DisplacementDim >::advective_heat_flux_contribution_to_K_liquid, ProcessLib::ThermoRichardsMechanics::PermeabilityData< DisplacementDim >::dk_rel_dS_L, ProcessLib::ThermoRichardsMechanics::SaturationDataDeriv::dS_L_dp_cap, ProcessLib::ConstitutiveRelations::SpaceTimeData::dt, MaterialPropertyLib::formEigenTensor(), ProcessLib::ThermoRichardsMechanics::TemperatureData< DisplacementDim >::grad_T, ProcessLib::ThermoRichardsMechanics::PermeabilityData< DisplacementDim >::k_rel, ProcessLib::ThermoRichardsMechanics::TRMHeatStorageAndFluxData< DisplacementDim >::K_Tp_NT_V_dN, ProcessLib::ThermoRichardsMechanics::TRMHeatStorageAndFluxData< DisplacementDim >::K_Tp_X_NTN, ProcessLib::ThermoRichardsMechanics::TRMHeatStorageAndFluxData< DisplacementDim >::K_TT_Laplace, ProcessLib::ThermoRichardsMechanics::PermeabilityData< DisplacementDim >::Ki, ProcessLib::ThermoRichardsMechanics::MediaData::liquid, MaterialPropertyLib::VariableArray::liquid_saturation, ProcessLib::ThermoRichardsMechanics::TRMHeatStorageAndFluxData< DisplacementDim >::M_TT_X_NTN, ProcessLib::ThermoRichardsMechanics::MediaData::medium, ProcessLib::ThermoRichardsMechanics::PorosityData::phi, MaterialPropertyLib::VariableArray::porosity, MaterialPropertyLib::Medium::property(), MaterialPropertyLib::Phase::property(), ProcessLib::ThermoRichardsMechanics::LiquidDensityData::rho_LR, ProcessLib::ThermoRichardsMechanics::SolidDensityData::rho_SR, ProcessLib::ThermoRichardsMechanics::SaturationData::S_L, ProcessLib::ThermoRichardsMechanics::MediaData::solid, MaterialPropertyLib::specific_heat_capacity, ProcessLib::ThermoRichardsMechanics::TemperatureData< DisplacementDim >::T, ProcessLib::ConstitutiveRelations::SpaceTimeData::t, MaterialPropertyLib::VariableArray::temperature, MaterialPropertyLib::thermal_conductivity, MaterialPropertyLib::Property::value(), and ProcessLib::ConstitutiveRelations::SpaceTimeData::x.

The documentation for this struct was generated from the following files: