OGS
Classes | Functions
ProcessLib::ThermoHydroMechanics Namespace Reference

Classes

struct  ConstitutiveRelationsValues
 
struct  IntegrationPointData
 
struct  IntegrationPointDataForOutput
 
struct  LocalAssemblerInterface
 
struct  SecondaryData
 
class  ThermoHydroMechanicsLocalAssembler
 
class  ThermoHydroMechanicsProcess
 
struct  ThermoHydroMechanicsProcessData
 

Functions

template<int DisplacementDim>
std::unique_ptr< Process > createThermoHydroMechanicsProcess (std::string const &name, MeshLib::Mesh &mesh, std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&jacobian_assembler, std::vector< ProcessVariable > const &variables, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &parameters, std::optional< ParameterLib::CoordinateSystem > const &local_coordinate_system, unsigned const integration_order, BaseLib::ConfigTree const &config, std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)
 
template std::unique_ptr< Process > createThermoHydroMechanicsProcess< 2 > (std::string const &name, MeshLib::Mesh &mesh, std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&jacobian_assembler, std::vector< ProcessVariable > const &variables, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &parameters, std::optional< ParameterLib::CoordinateSystem > const &local_coordinate_system, unsigned const integration_order, BaseLib::ConfigTree const &config, std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)
 
template std::unique_ptr< Process > createThermoHydroMechanicsProcess< 3 > (std::string const &name, MeshLib::Mesh &mesh, std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&jacobian_assembler, std::vector< ProcessVariable > const &variables, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &parameters, std::optional< ParameterLib::CoordinateSystem > const &local_coordinate_system, unsigned const integration_order, BaseLib::ConfigTree const &config, std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)
 

Function Documentation

◆ createThermoHydroMechanicsProcess()

template<int DisplacementDim>
std::unique_ptr< Process > ProcessLib::ThermoHydroMechanics::createThermoHydroMechanicsProcess ( std::string const & name,
MeshLib::Mesh & mesh,
std::unique_ptr< ProcessLib::AbstractJacobianAssembler > && jacobian_assembler,
std::vector< ProcessVariable > const & variables,
std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const & parameters,
std::optional< ParameterLib::CoordinateSystem > const & local_coordinate_system,
unsigned const integration_order,
BaseLib::ConfigTree const & config,
std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const & media )
Input File Parameter
prj__processes__process__type
Input File Parameter
prj__processes__process__THERMO_HYDRO_MECHANICS__coupling_scheme

Process Variables

Input File Parameter
prj__processes__process__THERMO_HYDRO_MECHANICS__process_variables

Primary process variables as they appear in the global component vector:

Input File Parameter
prj__processes__process__THERMO_HYDRO_MECHANICS__process_variables__temperature
Input File Parameter
prj__processes__process__THERMO_HYDRO_MECHANICS__process_variables__pressure
Input File Parameter
prj__processes__process__THERMO_HYDRO_MECHANICS__process_variables__displacement

Process Parameters

Input File Parameter
prj__processes__process__THERMO_HYDRO_MECHANICS__specific_body_force

Definition at line 33 of file CreateThermoHydroMechanicsProcess.cpp.

42{
44 config.checkConfigParameter("type", "THERMO_HYDRO_MECHANICS");
45 DBUG("Create ThermoHydroMechanicsProcess.");
46
47 auto const coupling_scheme =
49 config.getConfigParameterOptional<std::string>("coupling_scheme");
50 const bool use_monolithic_scheme =
51 !(coupling_scheme && (*coupling_scheme == "staggered"));
52
54
56 auto const pv_config = config.getConfigSubtree("process_variables");
57
58 ProcessVariable* variable_T;
59 ProcessVariable* variable_p;
60 ProcessVariable* variable_u;
61 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>
62 process_variables;
63 if (use_monolithic_scheme) // monolithic scheme.
64 {
67 auto per_process_variables = findProcessVariables(
68 variables, pv_config,
69 {
70 "temperature",
72 "pressure",
74 "displacement"});
75 variable_T = &per_process_variables[0].get();
76 variable_p = &per_process_variables[1].get();
77 variable_u = &per_process_variables[2].get();
78 process_variables.push_back(std::move(per_process_variables));
79 }
80 else // staggered scheme.
81 {
82 using namespace std::string_literals;
83 for (auto const& variable_name :
84 {"temperature"s, "pressure"s, "displacement"s})
85 {
86 auto per_process_variables =
87 findProcessVariables(variables, pv_config, {variable_name});
88 process_variables.push_back(std::move(per_process_variables));
89 }
90 variable_T = &process_variables[0][0].get();
91 variable_p = &process_variables[1][0].get();
92 variable_u = &process_variables[2][0].get();
93 }
94
95 if (variable_T->getShapeFunctionOrder() != 1)
96 {
97 OGS_FATAL(
98 "The shape function order of temperature must be 1 but its input "
99 "value in <process_variable><order> is {:d}. Please correct it.",
100 variable_T->getShapeFunctionOrder());
101 }
102 if (variable_p->getShapeFunctionOrder() != 1)
103 {
104 OGS_FATAL(
105 "The shape function order of pressure must be 1 but its input "
106 "value in <process_variable><order> is {:d}. Please correct it.",
107 variable_p->getShapeFunctionOrder());
108 }
109
110 DBUG("Associate displacement with process variable '{:s}'.",
111 variable_u->getName());
112
113 if (variable_u->getNumberOfGlobalComponents() != DisplacementDim)
114 {
115 OGS_FATAL(
116 "Number of components of the process variable '{:s}' is different "
117 "from the displacement dimension: got {:d}, expected {:d}",
118 variable_u->getName(),
119 variable_u->getNumberOfGlobalComponents(),
120 DisplacementDim);
121 }
122
123 DBUG("Associate pressure with process variable '{:s}'.",
124 variable_p->getName());
125 if (variable_p->getNumberOfGlobalComponents() != 1)
126 {
127 OGS_FATAL(
128 "Pressure process variable '{:s}' is not a scalar variable but has "
129 "{:d} components.",
130 variable_p->getName(),
131 variable_p->getNumberOfGlobalComponents());
132 }
133
134 DBUG("Associate temperature with process variable '{:s}'.",
135 variable_T->getName());
136 if (variable_T->getNumberOfGlobalComponents() != 1)
137 {
138 OGS_FATAL(
139 "temperature process variable '{:s}' is not a scalar variable but "
140 "has {:d} components.",
141 variable_T->getName(),
142 variable_T->getNumberOfGlobalComponents());
143 }
144
146 auto solid_constitutive_relations =
147 MaterialLib::Solids::createConstitutiveRelations<DisplacementDim>(
148 parameters, local_coordinate_system, materialIDs(mesh), config);
149
150 auto ice_constitutive_relation =
151 MaterialLib::Solids::createConstitutiveRelationIce<DisplacementDim>(
152 parameters, local_coordinate_system, config);
153
154 // Specific body force
155 Eigen::Matrix<double, DisplacementDim, 1> specific_body_force;
156 {
157 std::vector<double> const b =
159 config.getConfigParameter<std::vector<double>>(
160 "specific_body_force");
161 if (b.size() != DisplacementDim)
162 {
163 OGS_FATAL(
164 "The size of the specific body force vector does not match the "
165 "displacement dimension. Vector size is {:d}, displacement "
166 "dimension is {:d}",
167 b.size(), DisplacementDim);
168 }
169
170 std::copy_n(b.data(), b.size(), specific_body_force.data());
171 }
172
173 auto media_map =
174 MaterialPropertyLib::createMaterialSpatialDistributionMap(media, mesh);
175
176 // Initial stress conditions
177
178 auto initial_stress = ProcessLib::createInitialStress<DisplacementDim>(
179 config, parameters, mesh);
180
181 auto stabilizer = NumLib::createNumericalStabilization(mesh, config);
182
183 ThermoHydroMechanicsProcessData<DisplacementDim> process_data{
184 materialIDs(mesh),
185 std::move(media_map),
186 std::move(solid_constitutive_relations),
187 std::move(ice_constitutive_relation),
188 std::move(initial_stress),
189 specific_body_force,
190 std::move(stabilizer)};
191
192 SecondaryVariableCollection secondary_variables;
193
194 ProcessLib::createSecondaryVariables(config, secondary_variables);
195
196 return std::make_unique<ThermoHydroMechanicsProcess<DisplacementDim>>(
197 std::move(name), mesh, std::move(jacobian_assembler), parameters,
198 integration_order, std::move(process_variables),
199 std::move(process_data), std::move(secondary_variables),
200 use_monolithic_scheme);
201}
OGS_FATAL
#define OGS_FATAL(...)
Definition Error.h:26
DBUG
void DBUG(fmt::format_string< Args... > fmt, Args &&... args)
Definition Logging.h:30
ProcessLib::ProcessVariable::getShapeFunctionOrder
unsigned getShapeFunctionOrder() const
Definition ProcessVariable.h:119
ProcessLib::ProcessVariable::getName
std::string const & getName() const
Definition ProcessVariable.cpp:211
ProcessLib::ProcessVariable::getNumberOfGlobalComponents
int getNumberOfGlobalComponents() const
Returns the number of components of the process variable.
Definition ProcessVariable.h:93
MaterialLib::Solids::createConstitutiveRelations
std::map< int, std::shared_ptr< MaterialLib::Solids::MechanicsBase< DisplacementDim > > > createConstitutiveRelations(std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &parameters, std::optional< ParameterLib::CoordinateSystem > const &local_coordinate_system, MeshLib::PropertyVector< int > const *const material_ids, BaseLib::ConfigTree const &config)
Definition CreateConstitutiveRelation.cpp:110
MaterialLib::Solids::createConstitutiveRelationIce
std::unique_ptr< MaterialLib::Solids::MechanicsBase< DisplacementDim > > createConstitutiveRelationIce(std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &parameters, std::optional< ParameterLib::CoordinateSystem > const &local_coordinate_system, BaseLib::ConfigTree const &config)
Definition CreateConstitutiveRelationIce.cpp:22
MaterialPropertyLib::createMaterialSpatialDistributionMap
MaterialSpatialDistributionMap createMaterialSpatialDistributionMap(std::map< int, std::shared_ptr< Medium > > const &media, MeshLib::Mesh const &mesh)
Definition CreateMaterialSpatialDistributionMap.cpp:18
MeshLib::materialIDs
PropertyVector< int > const * materialIDs(Mesh const &mesh)
Definition Mesh.cpp:268
NumLib::createNumericalStabilization
NumericalStabilization createNumericalStabilization(MeshLib::Mesh const &mesh, BaseLib::ConfigTree const &config)
Definition CreateNumericalStabilization.cpp:20
ProcessLib::findProcessVariables
std::vector< std::reference_wrapper< ProcessVariable > > findProcessVariables(std::vector< ProcessVariable > const &variables, BaseLib::ConfigTree const &pv_config, std::initializer_list< std::string > tags)
Definition ProcessUtils.cpp:57
ProcessLib::createInitialStress
InitialStress createInitialStress(BaseLib::ConfigTree const &config, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &parameters, MeshLib::Mesh const &mesh, bool const mandatory_stress_type)
Definition CreateInitialStress.cpp:25
ProcessLib::createSecondaryVariables
void createSecondaryVariables(BaseLib::ConfigTree const &config, SecondaryVariableCollection &secondary_variables)
Definition CreateSecondaryVariables.cpp:18

References BaseLib::ConfigTree::checkConfigParameter(), MaterialLib::Solids::createConstitutiveRelationIce(), MaterialLib::Solids::createConstitutiveRelations(), ProcessLib::createInitialStress(), MaterialPropertyLib::createMaterialSpatialDistributionMap(), NumLib::createNumericalStabilization(), ProcessLib::createSecondaryVariables(), DBUG(), ProcessLib::findProcessVariables(), BaseLib::ConfigTree::getConfigParameter(), BaseLib::ConfigTree::getConfigParameterOptional(), BaseLib::ConfigTree::getConfigSubtree(), ProcessLib::ProcessVariable::getName(), ProcessLib::ProcessVariable::getNumberOfGlobalComponents(), ProcessLib::ProcessVariable::getShapeFunctionOrder(), and OGS_FATAL.

◆ createThermoHydroMechanicsProcess< 2 >()

template std::unique_ptr< Process > ProcessLib::ThermoHydroMechanics::createThermoHydroMechanicsProcess< 2 > ( std::string const & name,
MeshLib::Mesh & mesh,
std::unique_ptr< ProcessLib::AbstractJacobianAssembler > && jacobian_assembler,
std::vector< ProcessVariable > const & variables,
std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const & parameters,
std::optional< ParameterLib::CoordinateSystem > const & local_coordinate_system,
unsigned const integration_order,
BaseLib::ConfigTree const & config,
std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const & media )

Referenced by ProjectData::parseProcesses().

◆ createThermoHydroMechanicsProcess< 3 >()

template std::unique_ptr< Process > ProcessLib::ThermoHydroMechanics::createThermoHydroMechanicsProcess< 3 > ( std::string const & name,
MeshLib::Mesh & mesh,
std::unique_ptr< ProcessLib::AbstractJacobianAssembler > && jacobian_assembler,
std::vector< ProcessVariable > const & variables,
std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const & parameters,
std::optional< ParameterLib::CoordinateSystem > const & local_coordinate_system,
unsigned const integration_order,
BaseLib::ConfigTree const & config,
std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const & media )

Referenced by ProjectData::parseProcesses().