OGS
ProcessLib::TH2M::ConstitutiveRelations::FT2Model< DisplacementDim > Struct Template Reference

Detailed Description

template<int DisplacementDim>
struct ProcessLib::TH2M::ConstitutiveRelations::FT2Model< DisplacementDim >

Definition at line 68 of file TEquation.h.

#include <TEquation.h>

Public Member Functions

void eval (DarcyVelocityData< DisplacementDim > const &darcy_velocity_data, FluidDensityData const &fluid_density_data, FluidEnthalpyData const &fluid_enthalpy_data, FT2Data< DisplacementDim > &fT_2) const
 
void dEval (DarcyVelocityData< DisplacementDim > const &darcy_velocity_data, FluidDensityData const &fluid_density_data, FluidEnthalpyData const &fluid_enthalpy_data, PermeabilityData< DisplacementDim > const &permeability_data, PhaseTransitionData const &phase_transition_data, SpecificBodyForceData< DisplacementDim > const &specific_body_force, ViscosityData const &viscosity_data, FT2DerivativeData< DisplacementDim > &dfT_2) const
 

Member Function Documentation

◆ dEval()

template<int DisplacementDim>
void ProcessLib::TH2M::ConstitutiveRelations::FT2Model< DisplacementDim >::dEval ( DarcyVelocityData< DisplacementDim > const & darcy_velocity_data,
FluidDensityData const & fluid_density_data,
FluidEnthalpyData const & fluid_enthalpy_data,
PermeabilityData< DisplacementDim > const & permeability_data,
PhaseTransitionData const & phase_transition_data,
SpecificBodyForceData< DisplacementDim > const & specific_body_force,
ViscosityData const & viscosity_data,
FT2DerivativeData< DisplacementDim > & dfT_2 ) const

Definition at line 69 of file TEquation.cpp.

78{
79 auto const k_over_mu_G =
80 permeability_data.Ki * permeability_data.k_rel_G / viscosity_data.mu_GR;
81 auto const k_over_mu_L =
82 permeability_data.Ki * permeability_data.k_rel_L / viscosity_data.mu_LR;
83
84 dfT_2.dp_GR_Npart = phase_transition_data.drho_GR_dp_GR *
85 fluid_enthalpy_data.h_G * darcy_velocity_data.w_GS +
86 fluid_density_data.rho_GR * fluid_enthalpy_data.h_G *
87 k_over_mu_G * phase_transition_data.drho_GR_dp_GR *
88 specific_body_force();
89 dfT_2.dp_GR_gradNpart =
90 fluid_density_data.rho_GR * fluid_enthalpy_data.h_G * k_over_mu_G -
91 fluid_density_data.rho_LR * fluid_enthalpy_data.h_L * k_over_mu_L;
92
93 // From p_LR = p_GR - p_cap it follows for
94 // drho_LR/dp_GR = drho_LR/dp_LR * dp_LR/dp_GR
95 // = drho_LR/dp_LR * (dp_GR/dp_GR - dp_cap/dp_GR)
96 // = drho_LR/dp_LR * (1 - 0)
97 double const drho_LR_dp_cap = -phase_transition_data.drho_LR_dp_LR;
98
99 dfT_2.dp_cap_Npart =
100 -drho_LR_dp_cap * fluid_enthalpy_data.h_L * darcy_velocity_data.w_LS -
101 fluid_density_data.rho_LR * fluid_enthalpy_data.h_L * k_over_mu_L *
102 drho_LR_dp_cap * specific_body_force();
103 dfT_2.dp_cap_gradNpart =
104 fluid_density_data.rho_LR * fluid_enthalpy_data.h_L * k_over_mu_L;
105
106 dfT_2.dT = phase_transition_data.drho_GR_dT * fluid_enthalpy_data.h_G *
107 darcy_velocity_data.w_GS +
108 fluid_density_data.rho_GR * phase_transition_data.dh_G_dT *
109 darcy_velocity_data.w_GS +
110 phase_transition_data.drho_LR_dT * fluid_enthalpy_data.h_L *
111 darcy_velocity_data.w_LS +
112 fluid_density_data.rho_LR * phase_transition_data.dh_L_dT *
113 darcy_velocity_data.w_LS;
114 // TODO (naumov) + k_over_mu_G * drho_GR_dT * b + k_over_mu_L *
115 // drho_LR_dT * b
116}

References ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::dh_G_dT, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::dh_L_dT, ProcessLib::TH2M::ConstitutiveRelations::FT2DerivativeData< DisplacementDim >::dp_cap_gradNpart, ProcessLib::TH2M::ConstitutiveRelations::FT2DerivativeData< DisplacementDim >::dp_cap_Npart, ProcessLib::TH2M::ConstitutiveRelations::FT2DerivativeData< DisplacementDim >::dp_GR_gradNpart, ProcessLib::TH2M::ConstitutiveRelations::FT2DerivativeData< DisplacementDim >::dp_GR_Npart, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_GR_dp_GR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_GR_dT, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_LR_dp_LR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_LR_dT, ProcessLib::TH2M::ConstitutiveRelations::FT2DerivativeData< DisplacementDim >::dT, ProcessLib::TH2M::ConstitutiveRelations::FluidEnthalpyData::h_G, ProcessLib::TH2M::ConstitutiveRelations::FluidEnthalpyData::h_L, ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::k_rel_G, ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::k_rel_L, ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::Ki, ProcessLib::TH2M::ConstitutiveRelations::ViscosityData::mu_GR, ProcessLib::TH2M::ConstitutiveRelations::ViscosityData::mu_LR, ProcessLib::TH2M::ConstitutiveRelations::FluidDensityData::rho_GR, ProcessLib::TH2M::ConstitutiveRelations::FluidDensityData::rho_LR, ProcessLib::TH2M::ConstitutiveRelations::DarcyVelocityData< DisplacementDim >::w_GS, and ProcessLib::TH2M::ConstitutiveRelations::DarcyVelocityData< DisplacementDim >::w_LS.

◆ eval()

template<int DisplacementDim>
void ProcessLib::TH2M::ConstitutiveRelations::FT2Model< DisplacementDim >::eval ( DarcyVelocityData< DisplacementDim > const & darcy_velocity_data,
FluidDensityData const & fluid_density_data,
FluidEnthalpyData const & fluid_enthalpy_data,
FT2Data< DisplacementDim > & fT_2 ) const

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