ProcessLib::ThermoRichardsMechanics::GravityModel< DisplacementDim > Struct Template Reference

Detailed Description

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

Definition at line 27 of file Gravity.h.

#include <Gravity.h>

Public Member Functions
	GravityModel (Eigen::Vector< double, DisplacementDim > const &specific_body_force)
void	eval (PorosityData const &poro_data, SolidDensityData const &rho_S_data, LiquidDensityData const &rho_L_data, SaturationData const &S_L_data, SaturationDataDeriv const &dS_L_data, GravityData< DisplacementDim > &out) const

Private Attributes
Eigen::Vector< double, DisplacementDim > const	specific_body_force_

Constructor & Destructor Documentation

GravityModel()

template<int DisplacementDim>

ProcessLib::ThermoRichardsMechanics::GravityModel< DisplacementDim >::GravityModel ( Eigen::Vector< double, DisplacementDim > const & specific_body_force )

inlineexplicit

Definition at line 29 of file Gravity.h.

        : specific_body_force_(specific_body_force)
    {
    }

Member Function Documentation

eval()

template<int DisplacementDim>

void ProcessLib::ThermoRichardsMechanics::GravityModel< DisplacementDim >::eval	(	PorosityData const &	poro_data,
		SolidDensityData const &	rho_S_data,
		LiquidDensityData const &	rho_L_data,
		SaturationData const &	S_L_data,
		SaturationDataDeriv const &	dS_L_data,
		GravityData< DisplacementDim > &	out ) const

Definition at line 16 of file Gravity.cpp.

{
    auto const rho_SR = rho_S_data.rho_SR;
    auto const phi = poro_data.phi;
    auto const S_L = S_L_data.S_L;
    auto const rho_LR = rho_L_data.rho_LR;
    auto const b = specific_body_force_;
 
    double const rho = rho_SR * (1 - phi) + S_L * phi * rho_LR;
    out.volumetric_body_force = rho * b;
 
    // There is no minus in the Jacobian block, because volumetric_body_force is
    // subtracted from the residual vector.
    out.J_up_HT_V_N =
        poro_data.phi * rho_L_data.rho_LR * dS_L_data.dS_L_dp_cap * b;
}

References ProcessLib::ThermoRichardsMechanics::SaturationDataDeriv::dS_L_dp_cap, ProcessLib::ThermoRichardsMechanics::GravityData< DisplacementDim >::J_up_HT_V_N, ProcessLib::ThermoRichardsMechanics::PorosityData::phi, ProcessLib::ThermoRichardsMechanics::LiquidDensityData::rho_LR, ProcessLib::ThermoRichardsMechanics::SolidDensityData::rho_SR, ProcessLib::ThermoRichardsMechanics::SaturationData::S_L, and ProcessLib::ThermoRichardsMechanics::GravityData< DisplacementDim >::volumetric_body_force.

Member Data Documentation

specific_body_force_

template<int DisplacementDim>

Eigen::Vector<double, DisplacementDim> const ProcessLib::ThermoRichardsMechanics::GravityModel< DisplacementDim >::specific_body_force_

private

Definition at line 43 of file Gravity.h.

---

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

• ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/Gravity.h
• ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/Gravity.cpp