OGS
List of all members | Public Types | Public Member Functions | Static Public Attributes | Private Types | Private Member Functions | Static Private Member Functions | Private Attributes | Static Private Attributes
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim > Class Template Reference

Detailed Description

template<typename ShapeFunctionDisplacement, typename ShapeFunctionPressure, int DisplacementDim>
class ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >

Definition at line 52 of file ThermoHydroMechanicsFEM.h.

#include <ThermoHydroMechanicsFEM.h>

Inheritance diagram for ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >:
[legend]
Collaboration diagram for ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >:
[legend]

Public Types

using ShapeMatricesTypeDisplacement
 
using ShapeMatricesTypePressure
 
using GlobalDimMatrixType
 
using GlobalDimVectorType
 
using Invariants = MathLib::KelvinVector::Invariants<KelvinVectorSize>
 
using SymmetricTensor = Eigen::Matrix<double, KelvinVectorSize, 1>
 

Public Member Functions

 ThermoHydroMechanicsLocalAssembler (ThermoHydroMechanicsLocalAssembler const &)=delete
 
 ThermoHydroMechanicsLocalAssembler (ThermoHydroMechanicsLocalAssembler &&)=delete
 
 ThermoHydroMechanicsLocalAssembler (MeshLib::Element const &e, std::size_t const, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, ThermoHydroMechanicsProcessData< DisplacementDim > &process_data)
 
std::size_t setIPDataInitialConditions (std::string_view const name, double const *values, int const integration_order) override
 Returns number of read integration points.
 
void assemble (double const, double const, std::vector< double > const &, std::vector< double > const &, std::vector< double > &, std::vector< double > &, std::vector< double > &) override
 
void assembleWithJacobian (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &local_x_prev, std::vector< double > &local_rhs_data, std::vector< double > &local_Jac_data) override
 
void initializeConcrete () override
 
void setInitialConditionsConcrete (Eigen::VectorXd const local_x, double const t, int const process_id) override
 
void preTimestepConcrete (std::vector< double > const &, double const, double const) override
 
void postTimestepConcrete (Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, double const t, double const dt, int const) override
 
void computeSecondaryVariableConcrete (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev) override
 
Eigen::Map< const Eigen::RowVectorXd > getShapeMatrix (const unsigned integration_point) const override
 Provides the shape matrix at the given integration point.
 
std::vector< double > const & getIntPtDarcyVelocity (const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const override
 
std::vector< double > getSigma () const override
 
std::vector< double > const & getIntPtFluidDensity (const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const override
 
std::vector< double > const & getIntPtViscosity (const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const override
 
- Public Member Functions inherited from ProcessLib::LocalAssemblerInterface
virtual ~LocalAssemblerInterface ()=default
 
virtual void setInitialConditions (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, double const t, int const process_id)
 
virtual void initialize (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table)
 
virtual void preAssemble (double const, double const, std::vector< double > const &)
 
virtual void assembleForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)
 
virtual void assembleWithJacobianForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, int const process_id, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data)
 
virtual void computeSecondaryVariable (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_prev, int const process_id)
 
virtual void preTimestep (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, double const t, double const delta_t)
 
virtual void postTimestep (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
 
void postNonLinearSolver (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
 
virtual Eigen::Vector3d getFlux (MathLib::Point3d const &, double const, std::vector< double > const &) const
 
virtual Eigen::Vector3d getFlux (MathLib::Point3d const &, double const, std::vector< std::vector< double > > const &) const
 Fits to staggered scheme.
 
- Public Member Functions inherited from NumLib::ExtrapolatableElement
virtual ~ExtrapolatableElement ()=default
 

Static Public Attributes

static int const KelvinVectorSize
 
static constexpr auto & N_u_op
 

Private Types

using BMatricesType
 
using IpData
 

Private Member Functions

ConstitutiveRelationsValues< DisplacementDim > updateConstitutiveRelations (Eigen::Ref< Eigen::VectorXd const > const local_x, Eigen::Ref< Eigen::VectorXd const > const local_x_prev, ParameterLib::SpatialPosition const &x_position, double const t, double const dt, IpData &ip_data, IntegrationPointDataForOutput< DisplacementDim > &ip_data_output) const
 
std::size_t setSigma (double const *values)
 
std::vector< double > const & getIntPtSigma (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
 
std::vector< double > const & getIntPtSigmaIce (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
 
std::vector< double > getEpsilonM () const override
 
virtual std::vector< double > const & getIntPtEpsilonM (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
 
std::vector< double > getEpsilon () const override
 
virtual std::vector< double > const & getIntPtEpsilon (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
 
std::vector< double > const & getIntPtIceVolume (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
 
unsigned getNumberOfIntegrationPoints () const override
 
int getMaterialID () const override
 
std::vector< double > getMaterialStateVariableInternalState (std::function< std::span< double >(typename MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables &)> const &get_values_span, int const &n_components) const override
 
MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables const & getMaterialStateVariablesAt (unsigned integration_point) const override
 

Static Private Member Functions

template<typename SolutionVector >
static constexpr auto localDOF (SolutionVector const &x)
 

Private Attributes

ThermoHydroMechanicsProcessData< DisplacementDim > & _process_data
 
std::vector< IpData, Eigen::aligned_allocator< IpData > > _ip_data
 
std::vector< IntegrationPointDataForOutput< DisplacementDim >, Eigen::aligned_allocator< IntegrationPointDataForOutput< DisplacementDim > > > _ip_data_output
 
NumLib::GenericIntegrationMethod const & _integration_method
 
MeshLib::Element const & _element
 
bool const _is_axially_symmetric
 
SecondaryData< typename ShapeMatricesTypeDisplacement::ShapeMatrices::ShapeType > _secondary_data
 

Static Private Attributes

static const int temperature_index = 0
 
static const int temperature_size = ShapeFunctionPressure::NPOINTS
 
static const int pressure_index = ShapeFunctionPressure::NPOINTS
 
static const int pressure_size = ShapeFunctionPressure::NPOINTS
 
static const int displacement_index = ShapeFunctionPressure::NPOINTS * 2
 
static const int displacement_size
 

Member Typedef Documentation

◆ BMatricesType

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
using ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::BMatricesType
private
Initial value:
BMatrixPolicyType<ShapeFunctionDisplacement, DisplacementDim>

Definition at line 262 of file ThermoHydroMechanicsFEM.h.

◆ GlobalDimMatrixType

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
using ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::GlobalDimMatrixType
Initial value:
typename ShapeMatricesTypePressure::GlobalDimMatrixType
EigenFixedShapeMatrixPolicy::GlobalDimMatrixType
MatrixType< GlobalDim, GlobalDim > GlobalDimMatrixType
Definition ShapeMatrixPolicy.h:84

Definition at line 63 of file ThermoHydroMechanicsFEM.h.

◆ GlobalDimVectorType

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
using ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::GlobalDimVectorType
Initial value:
typename ShapeMatricesTypePressure::GlobalDimVectorType
EigenFixedShapeMatrixPolicy::GlobalDimVectorType
VectorType< GlobalDim > GlobalDimVectorType
Definition ShapeMatrixPolicy.h:85

Definition at line 66 of file ThermoHydroMechanicsFEM.h.

◆ Invariants

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
using ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::Invariants = MathLib::KelvinVector::Invariants<KelvinVectorSize>

Definition at line 71 of file ThermoHydroMechanicsFEM.h.

◆ IpData

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
using ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::IpData
private
Initial value:
IntegrationPointData<BMatricesType, ShapeMatricesTypeDisplacement,
ShapeMatricesTypePressure, DisplacementDim,
ShapeFunctionDisplacement::NPOINTS>
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler::ShapeMatricesTypeDisplacement
ShapeMatrixPolicyType< ShapeFunctionDisplacement, DisplacementDim > ShapeMatricesTypeDisplacement
Definition ThermoHydroMechanicsFEM.h:56
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler::BMatricesType
BMatrixPolicyType< ShapeFunctionDisplacement, DisplacementDim > BMatricesType
Definition ThermoHydroMechanicsFEM.h:262
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler::ShapeMatricesTypePressure
ShapeMatrixPolicyType< ShapeFunctionPressure, DisplacementDim > ShapeMatricesTypePressure
Definition ThermoHydroMechanicsFEM.h:60

Definition at line 264 of file ThermoHydroMechanicsFEM.h.

◆ ShapeMatricesTypeDisplacement

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
using ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::ShapeMatricesTypeDisplacement
Initial value:
ShapeMatrixPolicyType<ShapeFunctionDisplacement, DisplacementDim>

Definition at line 56 of file ThermoHydroMechanicsFEM.h.

◆ ShapeMatricesTypePressure

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
using ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::ShapeMatricesTypePressure
Initial value:
ShapeMatrixPolicyType<ShapeFunctionPressure, DisplacementDim>

Definition at line 60 of file ThermoHydroMechanicsFEM.h.

◆ SymmetricTensor

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
using ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::SymmetricTensor = Eigen::Matrix<double, KelvinVectorSize, 1>

Definition at line 73 of file ThermoHydroMechanicsFEM.h.

Constructor & Destructor Documentation

◆ ThermoHydroMechanicsLocalAssembler() [1/3]

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::ThermoHydroMechanicsLocalAssembler ( ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim > const & )
delete

◆ ThermoHydroMechanicsLocalAssembler() [2/3]

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::ThermoHydroMechanicsLocalAssembler ( ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim > && )
delete

◆ ThermoHydroMechanicsLocalAssembler() [3/3]

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::ThermoHydroMechanicsLocalAssembler ( MeshLib::Element const & e,
std::size_t const ,
NumLib::GenericIntegrationMethod const & integration_method,
bool const is_axially_symmetric,
ThermoHydroMechanicsProcessData< DisplacementDim > & process_data )

Definition at line 36 of file ThermoHydroMechanicsFEM-impl.h.

43 : _process_data(process_data),
44 _integration_method(integration_method),
45 _element(e),
46 _is_axially_symmetric(is_axially_symmetric)
47{
48 unsigned const n_integration_points =
49 _integration_method.getNumberOfPoints();
50
51 _ip_data.reserve(n_integration_points);
52 _ip_data_output.resize(n_integration_points);
53 _secondary_data.N_u.resize(n_integration_points);
54
55 auto const shape_matrices_u =
56 NumLib::initShapeMatrices<ShapeFunctionDisplacement,
57 ShapeMatricesTypeDisplacement,
58 DisplacementDim>(e, is_axially_symmetric,
59 _integration_method);
60
61 auto const shape_matrices_p =
62 NumLib::initShapeMatrices<ShapeFunctionPressure,
63 ShapeMatricesTypePressure, DisplacementDim>(
64 e, is_axially_symmetric, _integration_method);
65
66 auto const& solid_material =
67 MaterialLib::Solids::selectSolidConstitutiveRelation(
68 _process_data.solid_materials, _process_data.material_ids,
69 e.getID());
70
71 // Consistency check: if frozen liquid phase is given, then the constitutive
72 // relation for ice must also be given, and vice versa.
73 auto const& medium = _process_data.media_map.getMedium(_element.getID());
74 if (medium->hasPhase("FrozenLiquid") !=
75 (_process_data.ice_constitutive_relation != nullptr))
76 {
77 OGS_FATAL(
78 "Frozen liquid phase is {:s} and the solid material constitutive "
79 "relation for ice is {:s}. But both must be given (or both "
80 "omitted).",
81 medium->hasPhase("FrozenLiquid") ? "specified" : "not specified",
82 _process_data.ice_constitutive_relation != nullptr
83 ? "specified"
84 : "not specified");
85 }
86 for (unsigned ip = 0; ip < n_integration_points; ip++)
87 {
88 _ip_data.emplace_back(solid_material);
89 auto& ip_data = _ip_data[ip];
90 auto const& sm_u = shape_matrices_u[ip];
91 ip_data.integration_weight =
92 _integration_method.getWeightedPoint(ip).getWeight() *
93 sm_u.integralMeasure * sm_u.detJ;
94
95 ip_data.N_u = sm_u.N;
96 ip_data.dNdx_u = sm_u.dNdx;
97
98 ip_data.N = shape_matrices_p[ip].N;
99 ip_data.dNdx = shape_matrices_p[ip].dNdx;
100
101 _secondary_data.N_u[ip] = shape_matrices_u[ip].N;
102 }
103}
OGS_FATAL
#define OGS_FATAL(...)
Definition Error.h:26
MathLib::WeightedPoint::getWeight
double getWeight() const
Definition WeightedPoint.h:80
MeshLib::Element::getID
std::size_t getID() const
Returns the ID of the element.
Definition Element.h:89
NumLib::GenericIntegrationMethod::getWeightedPoint
MathLib::WeightedPoint const & getWeightedPoint(unsigned const igp) const
Definition GenericIntegrationMethod.h:36
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler::_process_data
ThermoHydroMechanicsProcessData< DisplacementDim > & _process_data
Definition ThermoHydroMechanicsFEM.h:392
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler::_ip_data_output
std::vector< IntegrationPointDataForOutput< DisplacementDim >, Eigen::aligned_allocator< IntegrationPointDataForOutput< DisplacementDim > > > _ip_data_output
Definition ThermoHydroMechanicsFEM.h:398
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler::_secondary_data
SecondaryData< typename ShapeMatricesTypeDisplacement::ShapeMatrices::ShapeType > _secondary_data
Definition ThermoHydroMechanicsFEM.h:405
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler::_ip_data
std::vector< IpData, Eigen::aligned_allocator< IpData > > _ip_data
Definition ThermoHydroMechanicsFEM.h:394
MaterialLib::Solids::selectSolidConstitutiveRelation
auto & selectSolidConstitutiveRelation(SolidMaterialsMap const &constitutive_relations, MeshLib::PropertyVector< int > const *const material_ids, std::size_t const element_id)
Definition SelectSolidConstitutiveRelation.h:36
NumLib::initShapeMatrices
std::vector< typename ShapeMatricesType::ShapeMatrices, Eigen::aligned_allocator< typename ShapeMatricesType::ShapeMatrices > > initShapeMatrices(MeshLib::Element const &e, bool const is_axially_symmetric, IntegrationMethod const &integration_method)
Definition InitShapeMatrices.h:53
ProcessLib::ThermoHydroMechanics::SecondaryData::N_u
std::vector< ShapeMatrixType, Eigen::aligned_allocator< ShapeMatrixType > > N_u
Definition ThermoHydroMechanicsFEM.h:47

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_element, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_integration_method, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data_output, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_process_data, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_secondary_data, MeshLib::Element::getID(), NumLib::GenericIntegrationMethod::getNumberOfPoints(), MathLib::WeightedPoint::getWeight(), NumLib::GenericIntegrationMethod::getWeightedPoint(), NumLib::initShapeMatrices(), ProcessLib::ThermoHydroMechanics::SecondaryData< ShapeMatrixType >::N_u, OGS_FATAL, and MaterialLib::Solids::selectSolidConstitutiveRelation().

Member Function Documentation

◆ assemble()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
void ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::assemble ( double const ,
double const ,
std::vector< double > const & ,
std::vector< double > const & ,
std::vector< double > & ,
std::vector< double > & ,
std::vector< double > &  )
inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 97 of file ThermoHydroMechanicsFEM.h.

103 {
104 OGS_FATAL(
105 "ThermoHydroMechanicsLocalAssembler: assembly without Jacobian is "
106 "not implemented.");
107 }

References OGS_FATAL.

◆ assembleWithJacobian()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
void ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::assembleWithJacobian ( double const t,
double const dt,
std::vector< double > const & local_x,
std::vector< double > const & local_x_prev,
std::vector< double > & local_rhs_data,
std::vector< double > & local_Jac_data )
overridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 553 of file ThermoHydroMechanicsFEM-impl.h.

559{
560 assert(local_x.size() ==
561 pressure_size + displacement_size + temperature_size);
562
563 auto const x =
564 Eigen::Map<Eigen::VectorXd const>(local_x.data(), local_x.size());
565 auto const x_prev = Eigen::Map<Eigen::VectorXd const>(local_x_prev.data(),
566 local_x_prev.size());
567
568 auto const [T, p, u] = localDOF(local_x);
569 auto const [T_prev, p_prev, u_prev] = localDOF(local_x_prev);
570
571 auto local_Jac = MathLib::createZeroedMatrix<
572 typename ShapeMatricesTypeDisplacement::template MatrixType<
573 temperature_size + displacement_size + pressure_size,
574 temperature_size + displacement_size + pressure_size>>(
575 local_Jac_data, displacement_size + pressure_size + temperature_size,
576 displacement_size + pressure_size + temperature_size);
577
578 auto local_rhs = MathLib::createZeroedVector<
579 typename ShapeMatricesTypeDisplacement::template VectorType<
580 displacement_size + pressure_size + temperature_size>>(
581 local_rhs_data, displacement_size + pressure_size + temperature_size);
582
583 typename ShapeMatricesTypePressure::NodalMatrixType MTT;
584 MTT.setZero(temperature_size, temperature_size);
585
586 typename ShapeMatricesTypePressure::NodalMatrixType KTT;
587 KTT.setZero(temperature_size, temperature_size);
588
589 typename ShapeMatricesTypePressure::NodalMatrixType KTp;
590 KTp.setZero(temperature_size, pressure_size);
591
592 typename ShapeMatricesTypePressure::NodalMatrixType dKTT_dp;
593 dKTT_dp.setZero(temperature_size, pressure_size);
594
595 typename ShapeMatricesTypePressure::NodalMatrixType laplace_p;
596 laplace_p.setZero(pressure_size, pressure_size);
597
598 typename ShapeMatricesTypePressure::NodalMatrixType laplace_T;
599 laplace_T.setZero(pressure_size, temperature_size);
600
601 typename ShapeMatricesTypePressure::NodalMatrixType storage_p;
602 storage_p.setZero(pressure_size, pressure_size);
603
604 typename ShapeMatricesTypePressure::NodalMatrixType storage_T;
605 storage_T.setZero(pressure_size, temperature_size);
606
607 typename ShapeMatricesTypeDisplacement::template MatrixType<
608 displacement_size, pressure_size>
609 Kup;
610 Kup.setZero(displacement_size, pressure_size);
611
612 auto const& medium = _process_data.media_map.getMedium(_element.getID());
613 bool const has_frozen_liquid_phase = medium->hasPhase("FrozenLiquid");
614
615 unsigned const n_integration_points =
616 _integration_method.getNumberOfPoints();
617
618 std::vector<GlobalDimVectorType> ip_flux_vector;
619 double average_velocity_norm = 0.0;
620 ip_flux_vector.reserve(n_integration_points);
621
622 for (unsigned ip = 0; ip < n_integration_points; ip++)
623 {
624 auto const& N_u = _ip_data[ip].N_u;
625 ParameterLib::SpatialPosition const x_position{
626 std::nullopt, _element.getID(), ip,
627 MathLib::Point3d(
628 NumLib::interpolateCoordinates<ShapeFunctionDisplacement,
629 ShapeMatricesTypeDisplacement>(
630 _element, N_u))};
631
632 auto const crv = updateConstitutiveRelations(
633 x, x_prev, x_position, t, dt, _ip_data[ip], _ip_data_output[ip]);
634
635 auto const& w = _ip_data[ip].integration_weight;
636
637 auto const& dNdx_u = _ip_data[ip].dNdx_u;
638
639 auto const& N = _ip_data[ip].N;
640 auto const& dNdx = _ip_data[ip].dNdx;
641
642 auto const T_int_pt = N.dot(T);
643
644 auto const x_coord =
645 NumLib::interpolateXCoordinate<ShapeFunctionDisplacement,
646 ShapeMatricesTypeDisplacement>(
647 _element, N_u);
648 auto const B =
649 LinearBMatrix::computeBMatrix<DisplacementDim,
650 ShapeFunctionDisplacement::NPOINTS,
651 typename BMatricesType::BMatrixType>(
652 dNdx_u, N_u, x_coord, _is_axially_symmetric);
653
654 auto const& b = _process_data.specific_body_force;
655 auto const velocity = _ip_data_output[ip].velocity;
656
657 //
658 // displacement equation, displacement part
659 //
660
661 if (has_frozen_liquid_phase)
662 {
663 local_Jac
664 .template block<displacement_size, displacement_size>(
665 displacement_index, displacement_index)
666 .noalias() += B.transpose() * crv.J_uu_fr * B * w;
667
668 local_rhs.template segment<displacement_size>(displacement_index)
669 .noalias() -= B.transpose() * crv.r_u_fr * w;
670
671 local_Jac
672 .template block<displacement_size, temperature_size>(
673 displacement_index, temperature_index)
674 .noalias() -= B.transpose() * crv.J_uT_fr * N * w;
675 }
676
677 local_Jac
678 .template block<displacement_size, displacement_size>(
679 displacement_index, displacement_index)
680 .noalias() += B.transpose() * crv.C * B * w;
681
682 local_Jac
683 .template block<displacement_size, temperature_size>(
684 displacement_index, temperature_index)
685 .noalias() -= B.transpose() * crv.C *
686 crv.solid_linear_thermal_expansion_coefficient * N *
687 w;
688
689 local_rhs.template segment<displacement_size>(displacement_index)
690 .noalias() -= (B.transpose() * _ip_data[ip].sigma_eff -
691 N_u_op(N_u).transpose() * crv.rho * b) *
692 w;
693
694 //
695 // displacement equation, pressure part (K_up)
696 //
697 Kup.noalias() +=
698 B.transpose() * crv.alpha_biot * Invariants::identity2 * N * w;
699
700 //
701 // pressure equation, pressure part (K_pp and M_pp).
702 //
703 laplace_p.noalias() += dNdx.transpose() * crv.K_over_mu * dNdx * w;
704
705 storage_p.noalias() +=
706 N.transpose() *
707 (_ip_data[ip].porosity * crv.fluid_compressibility +
708 (crv.alpha_biot - _ip_data[ip].porosity) * crv.beta_SR) *
709 N * w;
710
711 laplace_T.noalias() +=
712 dNdx.transpose() * crv.K_pT_thermal_osmosis * dNdx * w;
713 //
714 // RHS, pressure part
715 //
716 double const fluid_density = _ip_data_output[ip].fluid_density;
717 local_rhs.template segment<pressure_size>(pressure_index).noalias() +=
718 dNdx.transpose() * fluid_density * crv.K_over_mu * b * w;
719 //
720 // pressure equation, temperature part (M_pT)
721 //
722 storage_T.noalias() += N.transpose() * crv.beta * N * w;
723
724 //
725 // pressure equation, displacement part.
726 //
727 // Reusing Kup.transpose().
728
729 //
730 // temperature equation, temperature part.
731 //
732 KTT.noalias() +=
733 dNdx.transpose() * crv.effective_thermal_conductivity * dNdx * w;
734
735 ip_flux_vector.emplace_back(velocity * fluid_density * crv.c_f);
736 average_velocity_norm += velocity.norm();
737
738 if (has_frozen_liquid_phase)
739 {
740 local_Jac
741 .template block<temperature_size, temperature_size>(
742 temperature_index, temperature_index)
743 .noalias() -= N.transpose() * crv.J_TT_fr * N * w;
744 }
745
746 MTT.noalias() +=
747 N.transpose() * crv.effective_volumetric_heat_capacity * N * w;
748
749 //
750 // temperature equation, pressure part
751 //
752 KTp.noalias() +=
753 dNdx.transpose() * T_int_pt * crv.K_pT_thermal_osmosis * dNdx * w;
754
755 // linearized darcy
756 dKTT_dp.noalias() -= fluid_density * crv.c_f * N.transpose() *
757 (dNdx * T).transpose() * crv.K_over_mu * dNdx * w;
758
759 /* TODO (Joerg) Temperature changes due to thermal dilatation of the
760 * fluid, which are usually discarded as being very small.
761 * Zhou et al. (10.1016/S0020-7683(98)00089-4) states that:
762 * "Biot (1956) neglected this term and it is included here for
763 * completeness"
764 * Keeping the code here in the case these are needed for the named
765 * effects in the future.
766 if (fluid_compressibility != 0)
767 {
768 auto const C_el = _ip_data[ip].computeElasticTangentStiffness(
769 t, x_position, dt, static_cast<double>(T_int_pt));
770 auto const solid_skeleton_compressibility =
771 1 / solid_material.getBulkModulus(t, x_position, &C_el);
772 double const fluid_volumetric_thermal_expansion_coefficient =
773 MaterialPropertyLib::getLiquidThermalExpansivity(
774 liquid_phase, vars, fluid_density, x_position, t, dt);
775
776 KTT.noalias() +=
777 dNdx.transpose() *
778 (-T_int_pt * fluid_volumetric_thermal_expansion_coefficient *
779 K_pT_thermal_osmosis / fluid_compressibility) *
780 dNdx * w;
781
782 local_rhs.template segment<temperature_size>(temperature_index)
783 .noalias() +=
784 dNdx.transpose() *
785 (-T_int_pt * fluid_volumetric_thermal_expansion_coefficient /
786 fluid_compressibility) *
787 fluid_density * K_over_mu * b * w;
788 MTu part for rhs and Jacobian:
789 (-T_int_pt *
790 Invariants::trace(solid_linear_thermal_expansion_coefficient) /
791 solid_skeleton_compressibility) *
792 N.transpose() * identity2.transpose() * B * w;
793 KTp part for rhs and Jacobian:
794 dNdx.transpose() *
795 (T_int_pt * fluid_volumetric_thermal_expansion_coefficient *
796 K_over_mu / fluid_compressibility) *
797 dNdx * w;
798 }
799 */
800 }
801
802 NumLib::assembleAdvectionMatrix(
803 _process_data.stabilizer, _ip_data, ip_flux_vector,
804 average_velocity_norm / static_cast<double>(n_integration_points), KTT);
805
806 // temperature equation, temperature part
807 local_Jac
808 .template block<temperature_size, temperature_size>(temperature_index,
809 temperature_index)
810 .noalias() += KTT + MTT / dt;
811
812 // temperature equation, pressure part
813 local_Jac
814 .template block<temperature_size, pressure_size>(temperature_index,
815 pressure_index)
816 .noalias() += KTp + dKTT_dp;
817
818 // displacement equation, pressure part
819 local_Jac
820 .template block<displacement_size, pressure_size>(displacement_index,
821 pressure_index)
822 .noalias() -= Kup;
823
824 // pressure equation, temperature part.
825 local_Jac
826 .template block<pressure_size, temperature_size>(pressure_index,
827 temperature_index)
828 .noalias() -= storage_T / dt - laplace_T;
829
830 // pressure equation, pressure part.
831 local_Jac
832 .template block<pressure_size, pressure_size>(pressure_index,
833 pressure_index)
834 .noalias() += laplace_p + storage_p / dt;
835
836 // pressure equation, displacement part.
837 local_Jac
838 .template block<pressure_size, displacement_size>(pressure_index,
839 displacement_index)
840 .noalias() += Kup.transpose() / dt;
841
842 // pressure equation (f_p)
843 local_rhs.template segment<pressure_size>(pressure_index).noalias() -=
844 laplace_p * p + laplace_T * T + storage_p * (p - p_prev) / dt -
845 storage_T * (T - T_prev) / dt + Kup.transpose() * (u - u_prev) / dt;
846
847 // displacement equation (f_u)
848 local_rhs.template segment<displacement_size>(displacement_index)
849 .noalias() += Kup * p;
850
851 // temperature equation (f_T)
852 local_rhs.template segment<temperature_size>(temperature_index).noalias() -=
853 KTT * T + MTT * (T - T_prev) / dt;
854
855 local_rhs.template segment<temperature_size>(temperature_index).noalias() -=
856 KTp * p;
857}
ProcessLib::BMatrixPolicyType::BMatrixType
MatrixType< _kelvin_vector_size, _number_of_dof > BMatrixType
Definition BMatrixPolicy.h:55
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler::updateConstitutiveRelations
ConstitutiveRelationsValues< DisplacementDim > updateConstitutiveRelations(Eigen::Ref< Eigen::VectorXd const > const local_x, Eigen::Ref< Eigen::VectorXd const > const local_x_prev, ParameterLib::SpatialPosition const &x_position, double const t, double const dt, IpData &ip_data, IntegrationPointDataForOutput< DisplacementDim > &ip_data_output) const
Definition ThermoHydroMechanicsFEM-impl.h:224
MathLib::createZeroedVector
Eigen::Map< Vector > createZeroedVector(std::vector< double > &data, Eigen::VectorXd::Index size)
Definition EigenMapTools.h:153
MathLib::createZeroedMatrix
Eigen::Map< Matrix > createZeroedMatrix(std::vector< double > &data, Eigen::MatrixXd::Index rows, Eigen::MatrixXd::Index cols)
Definition EigenMapTools.h:32
NumLib::detail::assembleAdvectionMatrix
void assembleAdvectionMatrix(IPData const &ip_data_vector, NumLib::ShapeMatrixCache const &shape_matrix_cache, std::vector< FluxVectorType > const &ip_flux_vector, Eigen::MatrixBase< Derived > &laplacian_matrix)
Definition AdvectionMatrixAssembler.h:28
NumLib::interpolateXCoordinate
double interpolateXCoordinate(MeshLib::Element const &e, typename ShapeMatricesType::ShapeMatrices::ShapeType const &N)
Definition InitShapeMatrices.h:71
NumLib::interpolateCoordinates
std::array< double, 3 > interpolateCoordinates(MeshLib::Element const &e, typename ShapeMatricesType::ShapeMatrices::ShapeType const &N)
Definition InitShapeMatrices.h:82
ProcessLib::LinearBMatrix::computeBMatrix
BMatrixType computeBMatrix(DNDX_Type const &dNdx, N_Type const &N, const double radius, const bool is_axially_symmetric)
Fills a B-matrix based on given shape function dN/dx values.
Definition LinearBMatrix.h:46
ProcessLib::TES::fluid_density
double fluid_density(const double p, const double T, const double x)
Definition TESOGS5MaterialModels.h:19
EigenFixedShapeMatrixPolicy::NodalMatrixType
MatrixType< ShapeFunction::NPOINTS, ShapeFunction::NPOINTS > NodalMatrixType
Definition ShapeMatrixPolicy.h:74
MathLib::KelvinVector::Invariants::identity2
static Eigen::Matrix< double, KelvinVectorSize, 1 > const identity2
Kelvin mapping of 2nd order identity tensor.
Definition KelvinVector.h:107

References NumLib::detail::assembleAdvectionMatrix(), ProcessLib::LinearBMatrix::computeBMatrix(), MathLib::createZeroedMatrix(), MathLib::createZeroedVector(), NumLib::interpolateCoordinates(), and NumLib::interpolateXCoordinate().

◆ computeSecondaryVariableConcrete()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
void ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::computeSecondaryVariableConcrete ( double const t,
double const dt,
Eigen::VectorXd const & local_x,
Eigen::VectorXd const & local_x_prev )
overridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 918 of file ThermoHydroMechanicsFEM-impl.h.

922{
923 auto const p = local_x.template segment<pressure_size>(pressure_index);
924 auto const T =
925 local_x.template segment<temperature_size>(temperature_index);
926
927 unsigned const n_integration_points =
928 _integration_method.getNumberOfPoints();
929
930 double fluid_density_avg = 0;
931 double viscosity_avg = 0;
932
933 using KV = MathLib::KelvinVector::KelvinVectorType<DisplacementDim>;
934 KV sigma_avg = KV::Zero();
935
936 for (unsigned ip = 0; ip < n_integration_points; ip++)
937 {
938 auto& ip_data = _ip_data[ip];
939
940 fluid_density_avg += _ip_data_output[ip].fluid_density;
941 viscosity_avg += _ip_data_output[ip].viscosity;
942 sigma_avg += ip_data.sigma_eff;
943 }
944
945 fluid_density_avg /= n_integration_points;
946 viscosity_avg /= n_integration_points;
947 sigma_avg /= n_integration_points;
948
949 (*_process_data.element_fluid_density)[_element.getID()] =
950 fluid_density_avg;
951 (*_process_data.element_viscosity)[_element.getID()] = viscosity_avg;
952
953 Eigen::Map<KV>(&(*_process_data.element_stresses)[_element.getID() *
954 KV::RowsAtCompileTime]) =
955 MathLib::KelvinVector::kelvinVectorToSymmetricTensor(sigma_avg);
956
957 NumLib::interpolateToHigherOrderNodes<
958 ShapeFunctionPressure, typename ShapeFunctionDisplacement::MeshElement,
959 DisplacementDim>(_element, _is_axially_symmetric, p,
960 *_process_data.pressure_interpolated);
961
962 NumLib::interpolateToHigherOrderNodes<
963 ShapeFunctionPressure, typename ShapeFunctionDisplacement::MeshElement,
964 DisplacementDim>(_element, _is_axially_symmetric, T,
965 *_process_data.temperature_interpolated);
966}
MathLib::KelvinVector::kelvinVectorToSymmetricTensor
Eigen::Matrix< double, 4, 1 > kelvinVectorToSymmetricTensor(Eigen::Matrix< double, 4, 1, Eigen::ColMajor, 4, 1 > const &v)
Definition KelvinVector.cpp:142
MathLib::KelvinVector::KelvinVectorType
Eigen::Matrix< double, kelvin_vector_dimensions(DisplacementDim), 1, Eigen::ColMajor > KelvinVectorType
Definition KelvinVector.h:47
NumLib::interpolateToHigherOrderNodes
void interpolateToHigherOrderNodes(MeshLib::Element const &element, bool const is_axially_symmetric, Eigen::MatrixBase< EigenMatrixType > const &node_values, MeshLib::PropertyVector< double > &interpolated_values_global_vector)
Definition Interpolation.h:107

References NumLib::interpolateToHigherOrderNodes(), and MathLib::KelvinVector::kelvinVectorToSymmetricTensor().

◆ getEpsilon()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
std::vector< double > ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getEpsilon ( ) const
inlineoverrideprivatevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 323 of file ThermoHydroMechanicsFEM.h.

324 {
325 constexpr int kelvin_vector_size =
326 MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
327
328 return transposeInPlace<kelvin_vector_size>(
329 [this](std::vector<double>& values)
330 { return getIntPtEpsilon(0, {}, {}, values); });
331 }
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler::getIntPtEpsilon
virtual std::vector< double > const & getIntPtEpsilon(const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
Definition ThermoHydroMechanicsFEM.h:333
MathLib::KelvinVector::kelvin_vector_dimensions
constexpr int kelvin_vector_dimensions(int const displacement_dim)
Kelvin vector dimensions for given displacement dimension.
Definition KelvinVector.h:24
ProcessLib::transposeInPlace
std::vector< double > transposeInPlace(StoreValuesFunction const &store_values_function)
Definition TransposeInPlace.h:23

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtEpsilon(), MathLib::KelvinVector::kelvin_vector_dimensions(), and ProcessLib::transposeInPlace().

◆ getEpsilonM()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
std::vector< double > ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getEpsilonM ( ) const
inlineoverrideprivatevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 303 of file ThermoHydroMechanicsFEM.h.

304 {
305 constexpr int kelvin_vector_size =
306 MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
307
308 return transposeInPlace<kelvin_vector_size>(
309 [this](std::vector<double>& values)
310 { return getIntPtEpsilonM(0, {}, {}, values); });
311 }
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler::getIntPtEpsilonM
virtual std::vector< double > const & getIntPtEpsilonM(const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
Definition ThermoHydroMechanicsFEM.h:313

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtEpsilonM(), MathLib::KelvinVector::kelvin_vector_dimensions(), and ProcessLib::transposeInPlace().

◆ getIntPtDarcyVelocity()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
std::vector< double > const & ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtDarcyVelocity ( const double t,
std::vector< GlobalVector * > const & x,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & dof_table,
std::vector< double > & cache ) const
overridevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 862 of file ThermoHydroMechanicsFEM-impl.h.

868{
869 unsigned const n_integration_points =
870 _integration_method.getNumberOfPoints();
871
872 cache.clear();
873 auto cache_matrix = MathLib::createZeroedMatrix<Eigen::Matrix<
874 double, DisplacementDim, Eigen::Dynamic, Eigen::RowMajor>>(
875 cache, DisplacementDim, n_integration_points);
876
877 for (unsigned ip = 0; ip < n_integration_points; ip++)
878 {
879 cache_matrix.col(ip).noalias() = _ip_data_output[ip].velocity;
880 }
881
882 return cache;
883}

References MathLib::createZeroedMatrix().

◆ getIntPtEpsilon()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
virtual std::vector< double > const & ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtEpsilon ( const double ,
std::vector< GlobalVector * > const & ,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & ,
std::vector< double > & cache ) const
inlineoverrideprivatevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 333 of file ThermoHydroMechanicsFEM.h.

338 {
339 return ProcessLib::getIntegrationPointKelvinVectorData<DisplacementDim>(
340 _ip_data, &IpData::eps, cache);
341 }
ProcessLib::getIntegrationPointKelvinVectorData
std::vector< double > const & getIntegrationPointKelvinVectorData(IntegrationPointDataVector const &ip_data_vector, MemberType IpData::*const member, std::vector< double > &cache)
Definition SetOrGetIntegrationPointData.h:68

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::eps, and ProcessLib::getIntegrationPointKelvinVectorData().

Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getEpsilon().

◆ getIntPtEpsilonM()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
virtual std::vector< double > const & ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtEpsilonM ( const double ,
std::vector< GlobalVector * > const & ,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & ,
std::vector< double > & cache ) const
inlineoverrideprivatevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 313 of file ThermoHydroMechanicsFEM.h.

318 {
319 return ProcessLib::getIntegrationPointKelvinVectorData<DisplacementDim>(
320 _ip_data, &IpData::eps_m, cache);
321 }

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::eps_m, and ProcessLib::getIntegrationPointKelvinVectorData().

Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getEpsilonM().

◆ getIntPtFluidDensity()

template<typename ShapeFunctionDisplacement , typename ShapeFunction , int DisplacementDim>
std::vector< double > const & ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim >::getIntPtFluidDensity ( const double t,
std::vector< GlobalVector * > const & x,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & dof_table,
std::vector< double > & cache ) const
overridevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 888 of file ThermoHydroMechanicsFEM-impl.h.

894{
895 return ProcessLib::getIntegrationPointScalarData(
896 _ip_data_output,
897 &IntegrationPointDataForOutput<DisplacementDim>::fluid_density, cache);
898}
ProcessLib::getIntegrationPointScalarData
std::vector< double > const & getIntegrationPointScalarData(IntegrationPointDataVector const &ip_data_vector, MemberType IpData::*const member, std::vector< double > &cache)
Definition SetOrGetIntegrationPointData.h:176

References ProcessLib::getIntegrationPointScalarData().

◆ getIntPtIceVolume()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
std::vector< double > const & ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtIceVolume ( const double ,
std::vector< GlobalVector * > const & ,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & ,
std::vector< double > & cache ) const
inlineoverrideprivatevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 343 of file ThermoHydroMechanicsFEM.h.

348 {
349 return ProcessLib::getIntegrationPointScalarData(
350 _ip_data, &IpData::phi_fr, cache);
351 }

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data, ProcessLib::getIntegrationPointScalarData(), and ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::phi_fr.

◆ getIntPtSigma()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
std::vector< double > const & ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtSigma ( const double ,
std::vector< GlobalVector * > const & ,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & ,
std::vector< double > & cache ) const
inlineoverrideprivatevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 283 of file ThermoHydroMechanicsFEM.h.

288 {
289 return ProcessLib::getIntegrationPointKelvinVectorData<DisplacementDim>(
290 _ip_data, &IpData::sigma_eff, cache);
291 }

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data, ProcessLib::getIntegrationPointKelvinVectorData(), and ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::sigma_eff.

Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getSigma().

◆ getIntPtSigmaIce()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
std::vector< double > const & ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtSigmaIce ( const double ,
std::vector< GlobalVector * > const & ,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & ,
std::vector< double > & cache ) const
inlineoverrideprivatevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 293 of file ThermoHydroMechanicsFEM.h.

298 {
299 return ProcessLib::getIntegrationPointKelvinVectorData<DisplacementDim>(
300 _ip_data, &IpData::sigma_eff_ice, cache);
301 }

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data, ProcessLib::getIntegrationPointKelvinVectorData(), and ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::sigma_eff_ice.

◆ getIntPtViscosity()

template<typename ShapeFunctionDisplacement , typename ShapeFunction , int DisplacementDim>
std::vector< double > const & ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim >::getIntPtViscosity ( const double t,
std::vector< GlobalVector * > const & x,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & dof_table,
std::vector< double > & cache ) const
overridevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 903 of file ThermoHydroMechanicsFEM-impl.h.

909{
910 return ProcessLib::getIntegrationPointScalarData(
911 _ip_data_output,
912 &IntegrationPointDataForOutput<DisplacementDim>::viscosity, cache);
913}

References ProcessLib::getIntegrationPointScalarData().

◆ getMaterialID()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
int ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getMaterialID ( ) const
inlineoverrideprivatevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 358 of file ThermoHydroMechanicsFEM.h.

359 {
360 return _process_data.material_ids == nullptr
361 ? 0
362 : (*_process_data.material_ids)[_element.getID()];
363 }

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_element, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_process_data, and MeshLib::Element::getID().

◆ getMaterialStateVariableInternalState()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
std::vector< double > ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getMaterialStateVariableInternalState ( std::function< std::span< double >(typename MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables &)> const & get_values_span,
int const & n_components ) const
inlineoverrideprivatevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 365 of file ThermoHydroMechanicsFEM.h.

370 {
371 return ProcessLib::getIntegrationPointDataMaterialStateVariables(
372 _ip_data, &IpData::material_state_variables, get_values_span,
373 n_components);
374 }
ProcessLib::getIntegrationPointDataMaterialStateVariables
std::vector< double > getIntegrationPointDataMaterialStateVariables(IntegrationPointDataVector const &ip_data_vector, MemberType member, std::function< std::span< double >(MaterialStateVariables &)> get_values_span, int const n_components)
Definition SetOrGetIntegrationPointData.h:247
ProcessLib::ThermoHydroMechanics::IntegrationPointData::material_state_variables
std::unique_ptr< typename MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables > material_state_variables
Definition IntegrationPointData.h:71

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data, ProcessLib::getIntegrationPointDataMaterialStateVariables(), and ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::material_state_variables.

◆ getMaterialStateVariablesAt()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables const & ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getMaterialStateVariablesAt ( unsigned integration_point) const
inlineoverrideprivatevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 378 of file ThermoHydroMechanicsFEM.h.

379 {
380 return *_ip_data[integration_point].material_state_variables;
381 }

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data.

◆ getNumberOfIntegrationPoints()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
unsigned ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getNumberOfIntegrationPoints ( ) const
inlineoverrideprivatevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 353 of file ThermoHydroMechanicsFEM.h.

354 {
355 return _integration_method.getNumberOfPoints();
356 }

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_integration_method, and NumLib::GenericIntegrationMethod::getNumberOfPoints().

◆ getShapeMatrix()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
Eigen::Map< const Eigen::RowVectorXd > ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getShapeMatrix ( const unsigned integration_point) const
inlineoverridevirtual

Provides the shape matrix at the given integration point.

Implements NumLib::ExtrapolatableElement.

Definition at line 221 of file ThermoHydroMechanicsFEM.h.

223 {
224 auto const& N_u = _secondary_data.N_u[integration_point];
225
226 // assumes N is stored contiguously in memory
227 return Eigen::Map<const Eigen::RowVectorXd>(N_u.data(), N_u.size());
228 }

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_secondary_data, and ProcessLib::ThermoHydroMechanics::SecondaryData< ShapeMatrixType >::N_u.

◆ getSigma()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
std::vector< double > ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getSigma ( ) const
inlineoverridevirtual

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 239 of file ThermoHydroMechanicsFEM.h.

240 {
241 constexpr int kelvin_vector_size =
242 MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
243
244 return transposeInPlace<kelvin_vector_size>(
245 [this](std::vector<double>& values)
246 { return getIntPtSigma(0, {}, {}, values); });
247 }
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler::getIntPtSigma
std::vector< double > const & getIntPtSigma(const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
Definition ThermoHydroMechanicsFEM.h:283

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtSigma(), MathLib::KelvinVector::kelvin_vector_dimensions(), and ProcessLib::transposeInPlace().

◆ initializeConcrete()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
void ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::initializeConcrete ( )
inlineoverridevirtual

Set initial stress from parameter.

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 115 of file ThermoHydroMechanicsFEM.h.

116 {
117 unsigned const n_integration_points =
118 _integration_method.getNumberOfPoints();
119
120 for (unsigned ip = 0; ip < n_integration_points; ip++)
121 {
122 auto& ip_data = _ip_data[ip];
123
124 ParameterLib::SpatialPosition const x_position{
125 std::nullopt, _element.getID(), ip,
126 MathLib::Point3d(
127 NumLib::interpolateCoordinates<
128 ShapeFunctionDisplacement,
129 ShapeMatricesTypeDisplacement>(_element, ip_data.N_u))};
130
132 if (_process_data.initial_stress.value)
133 {
134 ip_data.sigma_eff =
135 MathLib::KelvinVector::symmetricTensorToKelvinVector<
136 DisplacementDim>((*_process_data.initial_stress.value)(
137 std::numeric_limits<double>::quiet_NaN() /* time
138 independent
139 */
140 ,
141 x_position));
142 }
143
144 double const t = 0; // TODO (naumov) pass t from top
145 ip_data.solid_material.initializeInternalStateVariables(
146 t, x_position, *ip_data.material_state_variables);
147
148 ip_data.pushBackState();
149 }
150 }
MathLib::KelvinVector::symmetricTensorToKelvinVector
Eigen::Matrix< double, Eigen::MatrixBase< Derived >::RowsAtCompileTime, 1 > symmetricTensorToKelvinVector(Eigen::MatrixBase< Derived > const &v)
Definition KelvinVector.h:201

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_element, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_integration_method, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_process_data, MeshLib::Element::getID(), NumLib::GenericIntegrationMethod::getNumberOfPoints(), NumLib::interpolateCoordinates(), and MathLib::KelvinVector::symmetricTensorToKelvinVector().

◆ localDOF()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
template<typename SolutionVector >
static constexpr auto ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::localDOF ( SolutionVector const & x)
inlinestaticconstexprprivate

Definition at line 385 of file ThermoHydroMechanicsFEM.h.

386 {
387 return NumLib::localDOF<
388 ShapeFunctionPressure, ShapeFunctionPressure,
389 NumLib::Vectorial<ShapeFunctionDisplacement, DisplacementDim>>(x);
390 }
NumLib::localDOF
auto localDOF(ElementDOFVector const &x)
Definition LocalDOF.h:64

References NumLib::localDOF().

Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::postTimestepConcrete().

◆ postTimestepConcrete()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
void ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::postTimestepConcrete ( Eigen::VectorXd const & local_x,
Eigen::VectorXd const & local_x_prev,
double const t,
double const dt,
int const  )
inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 169 of file ThermoHydroMechanicsFEM.h.

173 {
174 unsigned const n_integration_points =
175 _integration_method.getNumberOfPoints();
176
177 auto const [T_prev, p_prev, u_prev] = localDOF(local_x_prev);
178
179 for (unsigned ip = 0; ip < n_integration_points; ip++)
180 {
181 auto& ip_data = _ip_data[ip];
182 auto const& N_u = ip_data.N_u;
183 auto const& dNdx_u = ip_data.dNdx_u;
184
185 ParameterLib::SpatialPosition const x_position{
186 std::nullopt, _element.getID(), ip,
187 MathLib::Point3d(
188 NumLib::interpolateCoordinates<
189 ShapeFunctionDisplacement,
190 ShapeMatricesTypeDisplacement>(_element, N_u))};
191
192 updateConstitutiveRelations(local_x, local_x_prev, x_position, t,
193 dt, _ip_data[ip], _ip_data_output[ip]);
194
195 auto const x_coord =
196 NumLib::interpolateXCoordinate<ShapeFunctionDisplacement,
197 ShapeMatricesTypeDisplacement>(
198 _element, N_u);
199 auto const B = LinearBMatrix::computeBMatrix<
200 DisplacementDim, ShapeFunctionDisplacement::NPOINTS,
201 typename BMatricesType::BMatrixType>(dNdx_u, N_u, x_coord,
202 _is_axially_symmetric);
203
204 ConstitutiveRelationsValues<DisplacementDim> crv;
205
206 MathLib::KelvinVector::KelvinVectorType<DisplacementDim> const
207 eps_prev = B * u_prev;
208
209 _ip_data[ip].eps0 =
210 _ip_data[ip].eps0_prev +
211 (1 - _ip_data[ip].phi_fr_prev / _ip_data[ip].porosity) *
212 (eps_prev - _ip_data[ip].eps0_prev);
213 _ip_data[ip].pushBackState();
214 }
215 }

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_element, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_integration_method, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data_output, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_is_axially_symmetric, ProcessLib::LinearBMatrix::computeBMatrix(), MeshLib::Element::getID(), NumLib::GenericIntegrationMethod::getNumberOfPoints(), NumLib::interpolateCoordinates(), NumLib::interpolateXCoordinate(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::localDOF(), and ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::updateConstitutiveRelations().

◆ preTimestepConcrete()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
void ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::preTimestepConcrete ( std::vector< double > const & ,
double const ,
double const  )
inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 156 of file ThermoHydroMechanicsFEM.h.

158 {
159 unsigned const n_integration_points =
160 _integration_method.getNumberOfPoints();
161
162 for (unsigned ip = 0; ip < n_integration_points; ip++)
163 {
164 _ip_data_output[ip].velocity.setConstant(
165 DisplacementDim, std::numeric_limits<double>::quiet_NaN());
166 }
167 }

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_integration_method, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data_output, and NumLib::GenericIntegrationMethod::getNumberOfPoints().

◆ setInitialConditionsConcrete()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
void ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::setInitialConditionsConcrete ( Eigen::VectorXd const local_x,
double const t,
int const process_id )
overridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 179 of file ThermoHydroMechanicsFEM-impl.h.

183{
184 if (!_process_data.initial_stress.isTotalStress())
185 {
186 return;
187 }
188
189 // TODO: For staggered scheme, overload
190 // LocalAssemblerInterface::setInitialConditions to enable local_x contains
191 // the primary variables from all coupled processes.
192 auto const p = local_x.template segment<pressure_size>(pressure_index);
193
194 double const dt = 0.0;
195
196 MPL::VariableArray vars;
197 auto const& medium = _process_data.media_map.getMedium(_element.getID());
198
199 int const n_integration_points = _integration_method.getNumberOfPoints();
200 for (int ip = 0; ip < n_integration_points; ip++)
201 {
202 auto const& N = _ip_data[ip].N;
203 auto const& N_u = _ip_data[ip].N_u;
204 ParameterLib::SpatialPosition const x_position{
205 std::nullopt, _element.getID(), ip,
206 MathLib::Point3d(
207 NumLib::interpolateCoordinates<ShapeFunctionDisplacement,
208 ShapeMatricesTypeDisplacement>(
209 _element, N_u))};
210 auto const alpha_b =
211 medium->property(MPL::PropertyType::biot_coefficient)
212 .template value<double>(vars, x_position, t, dt);
213
214 auto& sigma_eff = _ip_data[ip].sigma_eff;
215 sigma_eff.noalias() += alpha_b * N.dot(p) * Invariants::identity2;
216 _ip_data[ip].sigma_eff_prev.noalias() = sigma_eff;
217 }
218}

References NumLib::interpolateCoordinates().

◆ setIPDataInitialConditions()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
std::size_t ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::setIPDataInitialConditions ( std::string_view const name,
double const * values,
int const integration_order )
overridevirtual

Returns number of read integration points.

Implements ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface< DisplacementDim >.

Definition at line 109 of file ThermoHydroMechanicsFEM-impl.h.

112{
113 if (integration_order !=
114 static_cast<int>(_integration_method.getIntegrationOrder()))
115 {
116 OGS_FATAL(
117 "Setting integration point initial conditions; The integration "
118 "order of the local assembler for element {:d} is different from "
119 "the integration order in the initial condition.",
120 _element.getID());
121 }
122
123 if (name == "sigma")
124 {
125 if (_process_data.initial_stress.value)
126 {
127 OGS_FATAL(
128 "Setting initial conditions for stress from integration "
129 "point data and from a parameter '{:s}' is not possible "
130 "simultaneously.",
131 _process_data.initial_stress.value->name);
132 }
133
134 return ProcessLib::setIntegrationPointKelvinVectorData<DisplacementDim>(
135 values, _ip_data, &IpData::sigma_eff);
136 }
137 if (name == "epsilon_m")
138 {
139 return ProcessLib::setIntegrationPointKelvinVectorData<DisplacementDim>(
140 values, _ip_data, &IpData::eps_m);
141 }
142 if (name == "epsilon")
143 {
144 return ProcessLib::setIntegrationPointKelvinVectorData<DisplacementDim>(
145 values, _ip_data, &IpData::eps);
146 }
147 if (name.starts_with("material_state_variable_"))
148 {
149 name.remove_prefix(24);
150
151 // Using first ip data for solid material. TODO (naumov) move solid
152 // material into element, store only material state in IPs.
153 auto const& internal_variables =
154 _ip_data[0].solid_material.getInternalVariables();
155 if (auto const iv = std::find_if(
156 begin(internal_variables), end(internal_variables),
157 [&name](auto const& iv) { return iv.name == name; });
158 iv != end(internal_variables))
159 {
160 DBUG("Setting material state variable '{:s}'", name);
161 return ProcessLib::setIntegrationPointDataMaterialStateVariables(
162 values, _ip_data, &IpData::material_state_variables,
163 iv->reference);
164 }
165
166 int const element_id = _element.getID();
167 DBUG(
168 "The solid material of element {:d} (material ID {:d}) does not "
169 "have an internal state variable called {:s}.",
170 element_id, (*_process_data.material_ids)[element_id], name);
171 }
172
173 return 0;
174}
DBUG
void DBUG(fmt::format_string< Args... > fmt, Args &&... args)
Definition Logging.h:30
ProcessLib::setIntegrationPointDataMaterialStateVariables
std::size_t setIntegrationPointDataMaterialStateVariables(double const *values, IntegrationPointDataVector &ip_data_vector, MemberType member, std::function< std::span< double >(MaterialStateVariables &)> get_values_span)
Definition SetOrGetIntegrationPointData.h:270
ProcessLib::setIntegrationPointKelvinVectorData
std::size_t setIntegrationPointKelvinVectorData(double const *values, IntegrationPointDataVector &ip_data_vector, MemberType IpData::*const member)
Definition SetOrGetIntegrationPointData.h:133

References DBUG(), OGS_FATAL, ProcessLib::setIntegrationPointDataMaterialStateVariables(), and ProcessLib::setIntegrationPointKelvinVectorData().

◆ setSigma()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
std::size_t ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::setSigma ( double const * values)
inlineprivate

Definition at line 277 of file ThermoHydroMechanicsFEM.h.

278 {
279 return ProcessLib::setIntegrationPointKelvinVectorData<DisplacementDim>(
280 values, _ip_data, &IpData::sigma_eff);
281 }

References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data, ProcessLib::setIntegrationPointKelvinVectorData(), and ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::sigma_eff.

◆ updateConstitutiveRelations()

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
ConstitutiveRelationsValues< DisplacementDim > ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::updateConstitutiveRelations ( Eigen::Ref< Eigen::VectorXd const > const local_x,
Eigen::Ref< Eigen::VectorXd const > const local_x_prev,
ParameterLib::SpatialPosition const & x_position,
double const t,
double const dt,
IpData & ip_data,
IntegrationPointDataForOutput< DisplacementDim > & ip_data_output ) const
private

Definition at line 223 of file ThermoHydroMechanicsFEM-impl.h.

230{
231 assert(local_x.size() ==
232 pressure_size + displacement_size + temperature_size);
233
234 auto const [T, p, u] = localDOF(local_x);
235 auto const [T_prev, p_prev, u_prev] = localDOF(local_x_prev);
236
237 auto const& solid_material =
238 MaterialLib::Solids::selectSolidConstitutiveRelation(
239 _process_data.solid_materials, _process_data.material_ids,
240 _element.getID());
241
242 auto const& medium = _process_data.media_map.getMedium(_element.getID());
243 auto const& liquid_phase = medium->phase("AqueousLiquid");
244 auto const& solid_phase = medium->phase("Solid");
245 auto* const frozen_liquid_phase = medium->hasPhase("FrozenLiquid")
246 ? &medium->phase("FrozenLiquid")
247 : nullptr;
248 MaterialPropertyLib::VariableArray vars;
249
250 auto const& N_u = ip_data.N_u;
251 auto const& dNdx_u = ip_data.dNdx_u;
252
253 auto const& N = ip_data.N;
254 auto const& dNdx = ip_data.dNdx;
255
256 auto const T_int_pt = N.dot(T);
257 auto const T_prev_int_pt = N.dot(T_prev);
258 double const dT_int_pt = T_int_pt - T_prev_int_pt;
259
260 auto const x_coord =
261 NumLib::interpolateXCoordinate<ShapeFunctionDisplacement,
262 ShapeMatricesTypeDisplacement>(_element,
263 N_u);
264 auto const B =
265 LinearBMatrix::computeBMatrix<DisplacementDim,
266 ShapeFunctionDisplacement::NPOINTS,
267 typename BMatricesType::BMatrixType>(
268 dNdx_u, N_u, x_coord, _is_axially_symmetric);
269
270 ConstitutiveRelationsValues<DisplacementDim> crv;
271
272 auto& eps = ip_data.eps;
273 eps.noalias() = B * u;
274 MathLib::KelvinVector::KelvinVectorType<DisplacementDim> const eps_prev =
275 B * u_prev;
276
277 vars.temperature = T_int_pt;
278 double const p_int_pt = N.dot(p);
279 vars.liquid_phase_pressure = p_int_pt;
280
281 vars.liquid_saturation = 1.0;
282
283 auto const solid_density =
284 solid_phase.property(MaterialPropertyLib::PropertyType::density)
285 .template value<double>(vars, x_position, t, dt);
286
287 auto const porosity =
288 medium->property(MaterialPropertyLib::PropertyType::porosity)
289 .template value<double>(vars, x_position, t, dt);
290 vars.porosity = porosity;
291 ip_data.porosity = porosity;
292
293 crv.alpha_biot =
294 medium->property(MaterialPropertyLib::PropertyType::biot_coefficient)
295 .template value<double>(vars, x_position, t, dt);
296 auto const& alpha = crv.alpha_biot;
297
298 auto const C_el = ip_data.computeElasticTangentStiffness(
299 t, x_position, dt, static_cast<double>(T_int_pt));
300 auto const solid_skeleton_compressibility =
301 1 / solid_material.getBulkModulus(t, x_position, &C_el);
302
303 crv.beta_SR = (1 - alpha) * solid_skeleton_compressibility;
304
305 // Set mechanical variables for the intrinsic permeability model
306 // For stress dependent permeability.
307 {
308 auto const& identity2 = Invariants::identity2;
309 auto const sigma_total =
310 (ip_data.sigma_eff - alpha * p_int_pt * identity2).eval();
311 vars.total_stress.emplace<SymmetricTensor>(
312 MathLib::KelvinVector::kelvinVectorToSymmetricTensor(sigma_total));
313 }
314 // For strain dependent permeability
315 vars.volumetric_strain = Invariants::trace(ip_data.eps);
316 vars.equivalent_plastic_strain =
317 ip_data.material_state_variables->getEquivalentPlasticStrain();
318
319 auto const intrinsic_permeability =
320 MaterialPropertyLib::formEigenTensor<DisplacementDim>(
321 medium->property(MaterialPropertyLib::PropertyType::permeability)
322 .value(vars, x_position, t, dt));
323
324 auto const fluid_density =
325 liquid_phase.property(MaterialPropertyLib::PropertyType::density)
326 .template value<double>(vars, x_position, t, dt);
327 ip_data_output.fluid_density = fluid_density;
328 vars.density = fluid_density;
329
330 auto const drho_dp =
331 liquid_phase.property(MaterialPropertyLib::PropertyType::density)
332 .template dValue<double>(
333 vars, MaterialPropertyLib::Variable::liquid_phase_pressure,
334 x_position, t, dt);
335
336 crv.fluid_compressibility = 1 / fluid_density * drho_dp;
337
338 double const fluid_volumetric_thermal_expansion_coefficient =
339 MaterialPropertyLib::getLiquidThermalExpansivity(
340 liquid_phase, vars, fluid_density, x_position, t, dt);
341
342 // Use the viscosity model to compute the viscosity
343 ip_data_output.viscosity =
344 liquid_phase.property(MaterialPropertyLib::PropertyType::viscosity)
345 .template value<double>(vars, x_position, t, dt);
346 crv.K_over_mu = intrinsic_permeability / ip_data_output.viscosity;
347
348 auto const& b = _process_data.specific_body_force;
349
350 // Consider also anisotropic thermal expansion.
351 crv.solid_linear_thermal_expansion_coefficient =
352 MPL::formKelvinVector<DisplacementDim>(
353 solid_phase
354 .property(
355 MaterialPropertyLib::PropertyType::thermal_expansivity)
356 .value(vars, x_position, t, dt));
357
358 MathLib::KelvinVector::KelvinVectorType<DisplacementDim> const
359 dthermal_strain =
360 crv.solid_linear_thermal_expansion_coefficient * dT_int_pt;
361
362 crv.K_pT_thermal_osmosis =
363 (solid_phase.hasProperty(
364 MaterialPropertyLib::PropertyType::thermal_osmosis_coefficient)
365 ? MaterialPropertyLib::formEigenTensor<DisplacementDim>(
366 solid_phase
367 .property(MaterialPropertyLib::PropertyType::
368 thermal_osmosis_coefficient)
369 .value(vars, x_position, t, dt))
370 : Eigen::MatrixXd::Zero(DisplacementDim, DisplacementDim));
371
372 GlobalDimVectorType const velocity = -crv.K_over_mu * dNdx * p -
373 crv.K_pT_thermal_osmosis * dNdx * T +
374 crv.K_over_mu * fluid_density * b;
375 ip_data_output.velocity = velocity;
376
377 //
378 // displacement equation, displacement part
379 //
380 auto& eps_m = ip_data.eps_m;
381 auto& eps_m_prev = ip_data.eps_m_prev;
382 eps_m.noalias() = eps_m_prev + eps - eps_prev - dthermal_strain;
383 vars.mechanical_strain
384 .emplace<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
385 eps_m);
386
387 crv.C = ip_data.updateConstitutiveRelation(vars, t, x_position, dt,
388 T_prev_int_pt);
389
390 crv.rho = solid_density * (1 - porosity) + porosity * fluid_density;
391
392 crv.beta =
393 porosity * fluid_volumetric_thermal_expansion_coefficient +
394 (alpha - porosity) *
395 Invariants::trace(crv.solid_linear_thermal_expansion_coefficient);
396
397 //
398 // pressure equation, displacement part.
399 //
400 // Reusing Kup.transpose().
401
402 //
403 // temperature equation, temperature part.
404 //
405 crv.c_f =
406 liquid_phase
407 .property(MaterialPropertyLib::PropertyType::specific_heat_capacity)
408 .template value<double>(vars, x_position, t, dt);
409 crv.effective_thermal_conductivity =
410 MaterialPropertyLib::formEigenTensor<DisplacementDim>(
411 medium
412 ->property(
413 MaterialPropertyLib::PropertyType::thermal_conductivity)
414 .value(vars, x_position, t, dt));
415
416 // Thermal conductivity is moved outside and zero matrix is passed instead
417 // due to multiplication with fluid's density times specific heat capacity.
418 crv.effective_thermal_conductivity.noalias() +=
419 fluid_density * crv.c_f *
420 NumLib::computeHydrodynamicDispersion(
421 _process_data.stabilizer, _element.getID(),
422 GlobalDimMatrixType::Zero(DisplacementDim, DisplacementDim),
423 velocity, 0. /* phi */, 0. /* dispersivity_transversal */,
424 0. /*dispersivity_longitudinal*/);
425
426 double const c_s =
427 solid_phase
428 .property(MaterialPropertyLib::PropertyType::specific_heat_capacity)
429 .template value<double>(vars, x_position, t, dt);
430
431 // Also modified by freezing terms.
432 crv.effective_volumetric_heat_capacity =
433 porosity * fluid_density * crv.c_f +
434 (1.0 - porosity) * solid_density * c_s;
435
436 if (frozen_liquid_phase)
437 {
438 MaterialPropertyLib::VariableArray vars_ice;
439 double const phi_fr =
440 (*medium)[MaterialPropertyLib::PropertyType::volume_fraction]
441 .template value<double>(vars, x_position, t, dt);
442 ip_data.phi_fr = phi_fr;
443
444 auto const frozen_liquid_value =
445 [&](MaterialPropertyLib::PropertyType const p)
446 {
447 return (*frozen_liquid_phase)[p].template value<double>(
448 vars, x_position, t, dt);
449 };
450
451 double const rho_fr =
452 frozen_liquid_value(MaterialPropertyLib::PropertyType::density);
453
454 double const c_fr = frozen_liquid_value(
455 MaterialPropertyLib::PropertyType::specific_heat_capacity);
456
457 double const l_fr = frozen_liquid_value(
458 MaterialPropertyLib::PropertyType::specific_latent_heat);
459
460 double const dphi_fr_dT =
461 (*medium)[MaterialPropertyLib::PropertyType::volume_fraction]
462 .template dValue<double>(
463 vars, MaterialPropertyLib::Variable::temperature,
464 x_position, t, dt);
465
466 double const phi_fr_prev = [&]()
467 {
468 MaterialPropertyLib::VariableArray vars_prev;
469 vars_prev.temperature = T_prev_int_pt;
470 return (*medium)[MaterialPropertyLib::PropertyType::volume_fraction]
471 .template value<double>(vars_prev, x_position, t, dt);
472 }();
473 ip_data.phi_fr_prev = phi_fr_prev;
474
475 // alpha_T^I
476 MathLib::KelvinVector::KelvinVectorType<
477 DisplacementDim> const ice_linear_thermal_expansion_coefficient =
478 MPL::formKelvinVector<DisplacementDim>(
479 frozen_liquid_phase
480 ->property(
481 MaterialPropertyLib::PropertyType::thermal_expansivity)
482 .value(vars, x_position, t, dt));
483
484 MathLib::KelvinVector::KelvinVectorType<DisplacementDim> const
485 dthermal_strain_ice =
486 ice_linear_thermal_expansion_coefficient * dT_int_pt;
487
488 // alpha_{phi_I} -- linear expansion coeff. due to water-to-ice
489 // transition (phase change), and related phase_change_strain term
490 MathLib::KelvinVector::KelvinVectorType<DisplacementDim> const
491 phase_change_expansion_coefficient =
492 MPL::formKelvinVector<DisplacementDim>(
493 frozen_liquid_phase
494 ->property(MaterialPropertyLib::PropertyType::
495 phase_change_expansivity)
496 .value(vars, x_position, t, dt));
497
498 MathLib::KelvinVector::KelvinVectorType<DisplacementDim> const
499 dphase_change_strain = phase_change_expansion_coefficient *
500 (phi_fr - phi_fr_prev) / porosity;
501
502 // eps0 ia a 'history variable' -- a solid matrix strain accrued
503 // prior to the onset of ice forming
504 auto& eps0 = ip_data.eps0;
505 auto const& eps0_prev = ip_data.eps0_prev;
506
507 // definition of eps_m_ice
508 auto& eps_m_ice = ip_data.eps_m_ice;
509 auto const& eps_m_ice_prev = ip_data.eps_m_ice_prev;
510
511 eps_m_ice.noalias() = eps_m_ice_prev + eps - eps_prev -
512 (eps0 - eps0_prev) - dthermal_strain_ice -
513 dphase_change_strain;
514
515 vars_ice.mechanical_strain
516 .emplace<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
517 eps_m_ice);
518 auto const C_IR = ip_data.updateConstitutiveRelationIce(
519 *_process_data.ice_constitutive_relation, vars_ice, t, x_position,
520 dt, T_prev_int_pt);
521 crv.effective_volumetric_heat_capacity +=
522 -phi_fr * fluid_density * crv.c_f + phi_fr * rho_fr * c_fr -
523 l_fr * rho_fr * dphi_fr_dT;
524
525 // part of dMTT_dT derivative for freezing
526 double const d2phi_fr_dT2 =
527 (*medium)[MaterialPropertyLib::PropertyType::volume_fraction]
528 .template d2Value<double>(
529 vars, MaterialPropertyLib::Variable::temperature,
530 MaterialPropertyLib::Variable::temperature, x_position, t,
531 dt);
532
533 crv.J_uu_fr = phi_fr * C_IR;
534
535 auto const& sigma_eff_ice = ip_data.sigma_eff_ice;
536 crv.r_u_fr = phi_fr * sigma_eff_ice;
537
538 crv.J_uT_fr = phi_fr * C_IR * ice_linear_thermal_expansion_coefficient;
539
540 crv.J_TT_fr = ((rho_fr * c_fr - fluid_density * crv.c_f) * dphi_fr_dT +
541 l_fr * rho_fr * d2phi_fr_dT2) *
542 dT_int_pt / dt;
543 }
544 return crv;
545}
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler::GlobalDimVectorType
typename ShapeMatricesTypePressure::GlobalDimVectorType GlobalDimVectorType
Definition ThermoHydroMechanicsFEM.h:66
MaterialPropertyLib::formKelvinVector
MathLib::KelvinVector::KelvinVectorType< GlobalDim > formKelvinVector(MaterialPropertyLib::PropertyDataType const &values)
A function to form a Kelvin vector from strain or stress alike property like thermal expansivity for ...
Definition FormKelvinVector.cpp:82
MaterialPropertyLib::formEigenTensor
Eigen::Matrix< double, GlobalDim, GlobalDim > formEigenTensor(MaterialPropertyLib::PropertyDataType const &values)
Definition FormEigenTensor.cpp:115
MaterialPropertyLib::getLiquidThermalExpansivity
double getLiquidThermalExpansivity(Phase const &phase, VariableArray const &vars, const double density, ParameterLib::SpatialPosition const &pos, double const t, double const dt)
Definition GetLiquidThermalExpansivity.cpp:21
NumLib::computeHydrodynamicDispersion
Eigen::MatrixXd computeHydrodynamicDispersion(NumericalStabilization const &stabilizer, std::size_t const element_id, Eigen::MatrixXd const &pore_diffusion_coefficient, Eigen::VectorXd const &velocity, double const porosity, double const solute_dispersivity_transverse, double const solute_dispersivity_longitudinal)
Definition HydrodynamicDispersion.h:71
ProcessLib::Graph::eval
auto eval(Function &f, Tuples &... ts) -> typename detail::GetFunctionReturnType< decltype(&Function::eval)>::type
Definition Apply.h:274
MathLib::KelvinVector::Invariants::trace
static double trace(Eigen::Matrix< double, KelvinVectorSize, 1 > const &v)
Trace of the corresponding tensor.
Definition KelvinVector-impl.h:62

References MaterialPropertyLib::biot_coefficient, ProcessLib::LinearBMatrix::computeBMatrix(), ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::computeElasticTangentStiffness(), NumLib::computeHydrodynamicDispersion(), MaterialPropertyLib::density, MaterialPropertyLib::VariableArray::density, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::dNdx, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::dNdx_u, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::eps, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::eps0, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::eps0_prev, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::eps_m, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::eps_m_ice, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::eps_m_ice_prev, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::eps_m_prev, MaterialPropertyLib::VariableArray::equivalent_plastic_strain, ProcessLib::ThermoHydroMechanics::IntegrationPointDataForOutput< DisplacementDim >::fluid_density, MaterialPropertyLib::formEigenTensor(), MaterialPropertyLib::formKelvinVector(), MaterialPropertyLib::getLiquidThermalExpansivity(), NumLib::interpolateXCoordinate(), MathLib::KelvinVector::kelvinVectorToSymmetricTensor(), MaterialPropertyLib::liquid_phase_pressure, MaterialPropertyLib::VariableArray::liquid_phase_pressure, MaterialPropertyLib::VariableArray::liquid_saturation, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::material_state_variables, MaterialPropertyLib::VariableArray::mechanical_strain, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::N, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::N_u, MaterialPropertyLib::permeability, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::phi_fr, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::phi_fr_prev, MaterialPropertyLib::porosity, MaterialPropertyLib::VariableArray::porosity, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::porosity, MaterialLib::Solids::selectSolidConstitutiveRelation(), ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::sigma_eff, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::sigma_eff_ice, MaterialPropertyLib::specific_heat_capacity, MaterialPropertyLib::specific_latent_heat, MaterialPropertyLib::temperature, MaterialPropertyLib::VariableArray::temperature, MaterialPropertyLib::thermal_conductivity, MaterialPropertyLib::thermal_expansivity, MaterialPropertyLib::thermal_osmosis_coefficient, MaterialPropertyLib::VariableArray::total_stress, ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::updateConstitutiveRelation(), ProcessLib::ThermoHydroMechanics::IntegrationPointData< BMatricesType, ShapeMatrixTypeDisplacement, ShapeMatricesTypePressure, DisplacementDim, NPoints >::updateConstitutiveRelationIce(), ProcessLib::ThermoHydroMechanics::IntegrationPointDataForOutput< DisplacementDim >::velocity, MaterialPropertyLib::viscosity, ProcessLib::ThermoHydroMechanics::IntegrationPointDataForOutput< DisplacementDim >::viscosity, MaterialPropertyLib::volume_fraction, and MaterialPropertyLib::VariableArray::volumetric_strain.

Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::postTimestepConcrete().

Member Data Documentation

◆ _element

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
MeshLib::Element const& ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_element
private

Definition at line 401 of file ThermoHydroMechanicsFEM.h.

Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::ThermoHydroMechanicsLocalAssembler(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getMaterialID(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::initializeConcrete(), and ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::postTimestepConcrete().

◆ _integration_method

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
NumLib::GenericIntegrationMethod const& ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_integration_method
private

Definition at line 400 of file ThermoHydroMechanicsFEM.h.

Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::ThermoHydroMechanicsLocalAssembler(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getNumberOfIntegrationPoints(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::initializeConcrete(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::postTimestepConcrete(), and ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::preTimestepConcrete().

◆ _ip_data

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
std::vector<IpData, Eigen::aligned_allocator<IpData> > ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data
private

Definition at line 394 of file ThermoHydroMechanicsFEM.h.

Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::ThermoHydroMechanicsLocalAssembler(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtEpsilon(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtEpsilonM(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtIceVolume(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtSigma(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getIntPtSigmaIce(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getMaterialStateVariableInternalState(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getMaterialStateVariablesAt(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::initializeConcrete(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::postTimestepConcrete(), and ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::setSigma().

◆ _ip_data_output

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
std::vector<IntegrationPointDataForOutput<DisplacementDim>, Eigen::aligned_allocator< IntegrationPointDataForOutput<DisplacementDim> > > ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_ip_data_output
private

Definition at line 398 of file ThermoHydroMechanicsFEM.h.

Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::ThermoHydroMechanicsLocalAssembler(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::postTimestepConcrete(), and ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::preTimestepConcrete().

◆ _is_axially_symmetric

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
bool const ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_is_axially_symmetric
private

Definition at line 402 of file ThermoHydroMechanicsFEM.h.

Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::postTimestepConcrete().

◆ _process_data

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
ThermoHydroMechanicsProcessData<DisplacementDim>& ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_process_data
private

Definition at line 392 of file ThermoHydroMechanicsFEM.h.

Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::ThermoHydroMechanicsLocalAssembler(), ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getMaterialID(), and ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::initializeConcrete().

◆ _secondary_data

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
SecondaryData< typename ShapeMatricesTypeDisplacement::ShapeMatrices::ShapeType> ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::_secondary_data
private

Definition at line 405 of file ThermoHydroMechanicsFEM.h.

Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::ThermoHydroMechanicsLocalAssembler(), and ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::getShapeMatrix().

◆ displacement_index

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
const int ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::displacement_index = ShapeFunctionPressure::NPOINTS * 2
staticprivate

Definition at line 413 of file ThermoHydroMechanicsFEM.h.

◆ displacement_size

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
const int ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::displacement_size
staticprivate
Initial value:
=
ShapeFunctionDisplacement::NPOINTS * DisplacementDim

Definition at line 414 of file ThermoHydroMechanicsFEM.h.

◆ KelvinVectorSize

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
int const ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::KelvinVectorSize
static
Initial value:
=
MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim)

Definition at line 69 of file ThermoHydroMechanicsFEM.h.

◆ N_u_op

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
auto& ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::N_u_op
staticconstexpr
Initial value:
= MathLib::eigenBlockMatrixView<
DisplacementDim,
typename ShapeMatricesTypeDisplacement::NodalRowVectorType>
MathLib::eigenBlockMatrixView
constexpr Eigen::CwiseNullaryOp< EigenBlockMatrixViewFunctor< D, M >, typename EigenBlockMatrixViewFunctor< D, M >::Matrix > eigenBlockMatrixView(const Eigen::MatrixBase< M > &matrix)
Definition EigenBlockMatrixView.h:49
EigenFixedShapeMatrixPolicy::NodalRowVectorType
RowVectorType< ShapeFunction::NPOINTS > NodalRowVectorType
Definition ShapeMatrixPolicy.h:78

Definition at line 75 of file ThermoHydroMechanicsFEM.h.

◆ pressure_index

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
const int ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::pressure_index = ShapeFunctionPressure::NPOINTS
staticprivate

Definition at line 411 of file ThermoHydroMechanicsFEM.h.

◆ pressure_size

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
const int ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::pressure_size = ShapeFunctionPressure::NPOINTS
staticprivate

Definition at line 412 of file ThermoHydroMechanicsFEM.h.

◆ temperature_index

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
const int ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::temperature_index = 0
staticprivate

Definition at line 409 of file ThermoHydroMechanicsFEM.h.

◆ temperature_size

template<typename ShapeFunctionDisplacement , typename ShapeFunctionPressure , int DisplacementDim>
const int ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::temperature_size = ShapeFunctionPressure::NPOINTS
staticprivate

Definition at line 410 of file ThermoHydroMechanicsFEM.h.

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