OGS
List of all members | Public Member Functions
ProcessLib::TH2M::ConstitutiveRelations::AdvectionModel< DisplacementDim > Struct Template Reference

Detailed Description

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

Definition at line 41 of file Advection.h.

#include <Advection.h>

Public Member Functions

void eval (ConstituentDensityData const &constituent_density_data, PermeabilityData< DisplacementDim > const &permeability_data, PureLiquidDensityData const &rho_W_LR, ViscosityData const &viscosity_data, AdvectionData< DisplacementDim > &advection_data) const
 
void dEval (ConstituentDensityData const &constituent_density_data, PermeabilityData< DisplacementDim > const &permeability_data, ViscosityData const &viscosity_data, SaturationDataDeriv const &dS_L_dp_cap, PhaseTransitionData const &phase_transition_data, AdvectionDerivativeData< DisplacementDim > &advection_d_data) const
 

Member Function Documentation

◆ dEval()

template<int DisplacementDim>
void ProcessLib::TH2M::ConstitutiveRelations::AdvectionModel< DisplacementDim >::dEval ( ConstituentDensityData const & constituent_density_data,
PermeabilityData< DisplacementDim > const & permeability_data,
ViscosityData const & viscosity_data,
SaturationDataDeriv const & dS_L_dp_cap,
PhaseTransitionData const & phase_transition_data,
AdvectionDerivativeData< DisplacementDim > & advection_d_data ) const

Definition at line 39 of file Advection.cpp.

46{
47 GlobalDimMatrix<DisplacementDim> const k_over_mu_G =
48 permeability_data.Ki * permeability_data.k_rel_G / viscosity_data.mu_GR;
49 GlobalDimMatrix<DisplacementDim> const k_over_mu_L =
50 permeability_data.Ki * permeability_data.k_rel_L / viscosity_data.mu_LR;
51
52 GlobalDimMatrix<DisplacementDim> const dk_over_mu_G_dp_cap =
53 permeability_data.Ki * permeability_data.dk_rel_G_dS_L * dS_L_dp_cap() /
54 viscosity_data.mu_GR;
55 GlobalDimMatrix<DisplacementDim> const dk_over_mu_L_dp_cap =
56 permeability_data.Ki * permeability_data.dk_rel_L_dS_L * dS_L_dp_cap() /
57 viscosity_data.mu_LR;
58
59 advection_d_data.dadvection_C_dp_GR =
60 phase_transition_data.drho_C_GR_dp_GR * k_over_mu_G
61 // + rhoCGR * (dk_over_mu_G_dp_GR = 0)
62 // + rhoCLR * (dk_over_mu_L_dp_GR = 0)
63 + phase_transition_data.drho_C_LR_dp_GR * k_over_mu_L;
64
65 advection_d_data.dadvection_C_dp_cap =
66 //(drho_C_GR_dp_cap = 0) * k_over_mu_G
67 constituent_density_data.rho_C_GR * dk_over_mu_G_dp_cap +
68 (-phase_transition_data.drho_C_LR_dp_LR) * k_over_mu_L +
69 constituent_density_data.rho_C_LR * dk_over_mu_L_dp_cap;
70}
ProcessLib::TH2M::ConstitutiveRelations::GlobalDimMatrix
Eigen::Matrix< double, DisplacementDim, DisplacementDim, Eigen::RowMajor > GlobalDimMatrix
Definition Base.h:32

References ProcessLib::TH2M::ConstitutiveRelations::AdvectionDerivativeData< DisplacementDim >::dadvection_C_dp_cap, ProcessLib::TH2M::ConstitutiveRelations::AdvectionDerivativeData< DisplacementDim >::dadvection_C_dp_GR, ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::dk_rel_G_dS_L, ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::dk_rel_L_dS_L, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_C_GR_dp_GR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_C_LR_dp_GR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_C_LR_dp_LR, 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::ConstituentDensityData::rho_C_GR, and ProcessLib::TH2M::ConstitutiveRelations::ConstituentDensityData::rho_C_LR.

◆ eval()

template<int DisplacementDim>
void ProcessLib::TH2M::ConstitutiveRelations::AdvectionModel< DisplacementDim >::eval ( ConstituentDensityData const & constituent_density_data,
PermeabilityData< DisplacementDim > const & permeability_data,
PureLiquidDensityData const & rho_W_LR,
ViscosityData const & viscosity_data,
AdvectionData< DisplacementDim > & advection_data ) const

Definition at line 17 of file Advection.cpp.

23{
24 GlobalDimMatrix<DisplacementDim> const k_over_mu_G =
25 permeability_data.Ki * permeability_data.k_rel_G / viscosity_data.mu_GR;
26 GlobalDimMatrix<DisplacementDim> const k_over_mu_L =
27 permeability_data.Ki * permeability_data.k_rel_L / viscosity_data.mu_LR;
28
29 advection_data.advection_C_G =
30 constituent_density_data.rho_C_GR * k_over_mu_G;
31 advection_data.advection_C_L =
32 constituent_density_data.rho_C_LR * k_over_mu_L;
33 advection_data.advection_W_G =
34 constituent_density_data.rho_W_GR * k_over_mu_G;
35 advection_data.advection_W_L = rho_W_LR() * k_over_mu_L;
36}

References ProcessLib::TH2M::ConstitutiveRelations::AdvectionData< DisplacementDim >::advection_C_G, ProcessLib::TH2M::ConstitutiveRelations::AdvectionData< DisplacementDim >::advection_C_L, ProcessLib::TH2M::ConstitutiveRelations::AdvectionData< DisplacementDim >::advection_W_G, ProcessLib::TH2M::ConstitutiveRelations::AdvectionData< DisplacementDim >::advection_W_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::ConstituentDensityData::rho_C_GR, ProcessLib::TH2M::ConstitutiveRelations::ConstituentDensityData::rho_C_LR, and ProcessLib::TH2M::ConstitutiveRelations::ConstituentDensityData::rho_W_GR.

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