OGS
ThermoRichardsMechanicsProcess.h
Go to the documentation of this file.
1
11#pragma once
12
14#ifdef OGS_USE_MFRONT
16#endif
19#include "ProcessLib/Process.h"
21
22namespace ProcessLib
23{
24namespace ThermoRichardsMechanics
25{
113template <int DisplacementDim, typename ConstitutiveTraits>
115 : public Process,
116 private AssemblyMixin<
117 ThermoRichardsMechanicsProcess<DisplacementDim, ConstitutiveTraits>>
118{
119 friend class AssemblyMixin<
120 ThermoRichardsMechanicsProcess<DisplacementDim, ConstitutiveTraits>>;
121
122public:
124 std::string name,
125 MeshLib::Mesh& mesh,
126 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&&
127 jacobian_assembler,
128 std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const&
129 parameters,
130 unsigned const integration_order,
131 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
132 process_variables,
134 ConstitutiveTraits>&& process_data,
135 SecondaryVariableCollection&& secondary_variables,
136 bool const use_monolithic_scheme);
137
140
141 bool isLinear() const override;
143
157 const int process_id) const override;
158
159private:
162
163 void constructDofTable() override;
164
166 NumLib::LocalToGlobalIndexMap const& dof_table,
167 MeshLib::Mesh const& mesh,
168 unsigned const integration_order) override;
169
171 std::map<int, std::shared_ptr<MaterialPropertyLib::Medium>> const&
172 media) override;
173
174 void setInitialConditionsConcreteProcess(std::vector<GlobalVector*>& x,
175 double const t,
176 int const process_id) override;
177
178 void assembleConcreteProcess(const double t, double const dt,
179 std::vector<GlobalVector*> const& x,
180 std::vector<GlobalVector*> const& x_prev,
181 int const process_id, GlobalMatrix& M,
182 GlobalMatrix& K, GlobalVector& b) override;
183
185 const double t, double const dt, std::vector<GlobalVector*> const& x,
186 std::vector<GlobalVector*> const& x_prev, int const process_id,
188 GlobalMatrix& Jac) override;
189
190 void preTimestepConcreteProcess(std::vector<GlobalVector*> const& /*x*/,
191 const double /*t*/,
192 const double /*dt*/,
193 const int /*process_id*/) override;
194
195 void postTimestepConcreteProcess(std::vector<GlobalVector*> const& x,
196 std::vector<GlobalVector*> const& x_prev,
197 double const t, double const dt,
198 const int process_id) override;
199
201 const int process_id) const override;
202
203 std::vector<std::string> initializeAssemblyOnSubmeshes(
204 std::vector<std::reference_wrapper<MeshLib::Mesh>> const& meshes)
205 override;
206
207private:
208 std::vector<MeshLib::Node*> base_nodes_;
209 std::unique_ptr<MeshLib::MeshSubset const> mesh_subset_base_nodes_;
212
213 std::vector<std::unique_ptr<LocalAssemblerIF>> local_assemblers_;
214
215 std::unique_ptr<NumLib::LocalToGlobalIndexMap>
217
220 std::unique_ptr<NumLib::LocalToGlobalIndexMap>
222
226
227 void computeSecondaryVariableConcrete(double const t, double const dt,
228 std::vector<GlobalVector*> const& x,
229 GlobalVector const& x_prev,
230 int const process_id) override;
234 std::tuple<NumLib::LocalToGlobalIndexMap*, bool>
235 getDOFTableForExtrapolatorData() const override;
236
240 bool hasMechanicalProcess(int const /*process_id*/) const { return true; }
241};
242
243extern template class ThermoRichardsMechanicsProcess<
244 2, ConstitutiveStress_StrainTemperature::ConstitutiveTraits<2>>;
245extern template class ThermoRichardsMechanicsProcess<
246 3, ConstitutiveStress_StrainTemperature::ConstitutiveTraits<3>>;
247
248#ifdef OGS_USE_MFRONT
249extern template class ThermoRichardsMechanicsProcess<
250 2, ConstitutiveStressSaturation_StrainPressureTemperature::
251 ConstitutiveTraits<2>>;
252extern template class ThermoRichardsMechanicsProcess<
253 3, ConstitutiveStressSaturation_StrainPressureTemperature::
254 ConstitutiveTraits<3>>;
255#endif
256
257} // namespace ThermoRichardsMechanics
258} // namespace ProcessLib
MathLib::SparsityPattern< GlobalIndexType > GlobalSparsityPattern
Global vector based on Eigen vector.
Definition EigenVector.h:25
std::string const name
Definition Process.h:355
Handles configuration of several secondary variables from the project file.
Global assembler for the monolithic scheme of the non-isothermal Richards flow coupled with mechanics...
void assembleConcreteProcess(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b) override
void preTimestepConcreteProcess(std::vector< GlobalVector * > const &, const double, const double, const int) override
void computeSecondaryVariableConcrete(double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_prev, int const process_id) override
ThermoRichardsMechanicsProcessData< DisplacementDim, ConstitutiveTraits > process_data_
ThermoRichardsMechanicsProcess(std::string name, MeshLib::Mesh &mesh, std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&jacobian_assembler, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &parameters, unsigned const integration_order, std::vector< std::vector< std::reference_wrapper< ProcessVariable > > > &&process_variables, ThermoRichardsMechanicsProcessData< DisplacementDim, ConstitutiveTraits > &&process_data, SecondaryVariableCollection &&secondary_variables, bool const use_monolithic_scheme)
NumLib::LocalToGlobalIndexMap const & getDOFTable(const int process_id) const override
void assembleWithJacobianConcreteProcess(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b, GlobalMatrix &Jac) override
std::vector< std::string > initializeAssemblyOnSubmeshes(std::vector< std::reference_wrapper< MeshLib::Mesh > > const &meshes) override
std::tuple< NumLib::LocalToGlobalIndexMap *, bool > getDOFTableForExtrapolatorData() const override
std::unique_ptr< NumLib::LocalToGlobalIndexMap > local_to_global_index_map_with_base_nodes_
void setInitialConditionsConcreteProcess(std::vector< GlobalVector * > &x, double const t, int const process_id) override
void initializeBoundaryConditions(std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media) override
std::unique_ptr< NumLib::LocalToGlobalIndexMap > local_to_global_index_map_single_component_
void initializeConcreteProcess(NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order) override
Process specific initialization called by initialize().
void postTimestepConcreteProcess(std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, const int process_id) override
MathLib::MatrixSpecifications getMatrixSpecifications(const int process_id) const override