28namespace ThermoRichardsMechanics
30template <
int DisplacementDim,
typename ConstitutiveTraits>
35 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&&
37 std::vector<std::unique_ptr<ParameterLib::ParameterBase>>
const&
39 unsigned const integration_order,
40 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
43 ConstitutiveTraits>&& process_data,
45 bool const use_monolithic_scheme)
46 :
Process(std::move(name), mesh, std::move(jacobian_assembler), parameters,
47 integration_order, std::move(process_variables),
48 std::move(secondary_variables), use_monolithic_scheme),
52 process_data_(std::move(process_data))
55 DisplacementDim>(LocalAssemblerIF::getReflectionDataForOutput(),
56 _integration_point_writer, integration_order,
60template <
int DisplacementDim,
typename ConstitutiveTraits>
61bool ThermoRichardsMechanicsProcess<DisplacementDim,
62 ConstitutiveTraits>::isLinear()
const
67template <
int DisplacementDim,
typename ConstitutiveTraits>
70 ConstitutiveTraits>::getMatrixSpecifications(
const int )
const
72 auto const& l = *_local_to_global_index_map;
73 return {l.dofSizeWithoutGhosts(), l.dofSizeWithoutGhosts(),
74 &l.getGhostIndices(), &this->_sparsity_pattern};
77template <
int DisplacementDim,
typename ConstitutiveTraits>
79 ConstitutiveTraits>::constructDofTable()
82 _mesh_subset_all_nodes = std::make_unique<MeshLib::MeshSubset>(
83 _mesh, _mesh.getNodes(), process_data_.use_TaylorHood_elements);
86 mesh_subset_base_nodes_ = std::make_unique<MeshLib::MeshSubset>(
87 _mesh, base_nodes_, process_data_.use_TaylorHood_elements);
91 std::vector<MeshLib::MeshSubset> all__meshsubsets_single_component{
92 *_mesh_subset_all_nodes};
93 local_to_global_index_map_single_component_ =
94 std::make_unique<NumLib::LocalToGlobalIndexMap>(
95 std::move(all__meshsubsets_single_component),
100 std::vector<MeshLib::MeshSubset> all__meshsubsets{*mesh_subset_base_nodes_};
103 all__meshsubsets.push_back(*mesh_subset_base_nodes_);
106 const int monolithic_process_id = 0;
107 std::generate_n(std::back_inserter(all__meshsubsets),
108 getProcessVariables(monolithic_process_id)[2]
110 .getNumberOfGlobalComponents(),
111 [&]() {
return *_mesh_subset_all_nodes; });
113 std::vector<int>
const vec_n_components{1, 1, DisplacementDim};
114 _local_to_global_index_map =
115 std::make_unique<NumLib::LocalToGlobalIndexMap>(
116 std::move(all__meshsubsets), vec_n_components,
118 assert(_local_to_global_index_map);
121template <
int DisplacementDim,
typename ConstitutiveTraits>
125 unsigned const integration_order)
127 createLocalAssemblers<DisplacementDim, ConstitutiveTraits>(
132 auto add_secondary_variable = [&](std::string
const& name,
133 int const num_components,
134 auto get_ip_values_function)
136 _secondary_variables.addSecondaryVariable(
140 std::move(get_ip_values_function)));
143 ProcessLib::Reflection::addReflectedSecondaryVariables<DisplacementDim>(
144 LocalAssemblerIF::getReflectionDataForOutput(), _secondary_variables,
145 getExtrapolator(), local_assemblers_);
152 add_secondary_variable);
156 process_data_.solid_materials, local_assemblers_,
157 _integration_point_writer, integration_order);
159 process_data_.element_saturation = MeshLib::getOrCreateMeshProperty<double>(
163 process_data_.element_porosity = MeshLib::getOrCreateMeshProperty<double>(
167 process_data_.element_liquid_density =
168 MeshLib::getOrCreateMeshProperty<double>(
172 process_data_.element_viscosity = MeshLib::getOrCreateMeshProperty<double>(
176 process_data_.element_stresses = MeshLib::getOrCreateMeshProperty<double>(
180 DisplacementDim>::RowsAtCompileTime);
182 process_data_.pressure_interpolated =
183 MeshLib::getOrCreateMeshProperty<double>(
186 process_data_.temperature_interpolated =
187 MeshLib::getOrCreateMeshProperty<double>(
188 const_cast<MeshLib::Mesh&
>(mesh),
"temperature_interpolated",
197 *_local_to_global_index_map);
200template <
int DisplacementDim,
typename ConstitutiveTraits>
203 std::map<
int, std::shared_ptr<MaterialPropertyLib::Medium>>
const&
206 const int process_id = 0;
207 initializeProcessBoundaryConditionsAndSourceTerms(
208 *_local_to_global_index_map, process_id, media);
211template <
int DisplacementDim,
typename ConstitutiveTraits>
215 int const process_id)
217 DBUG(
"SetInitialConditions ThermoRichardsMechanicsProcess.");
220 &LocalAssemblerIF::setInitialConditions, local_assemblers_,
221 getDOFTables(x.size()), x, t, process_id);
224template <
int DisplacementDim,
typename ConstitutiveTraits>
227 std::vector<GlobalVector*>
const& ,
228 std::vector<GlobalVector*>
const& ,
233 "The Picard method or the Newton-Raphson method with numerical "
234 "Jacobian is not implemented for ThermoRichardsMechanics with the full "
235 "monolithic coupling scheme");
238template <
int DisplacementDim,
typename ConstitutiveTraits>
241 const double t,
double const dt, std::vector<GlobalVector*>
const& x,
242 std::vector<GlobalVector*>
const& x_prev,
int const process_id,
246 DisplacementDim, ConstitutiveTraits>>::assembleWithJacobian(t, dt, x,
253template <
int DisplacementDim,
typename ConstitutiveTraits>
258 const int process_id)
261 DisplacementDim, ConstitutiveTraits>>::updateActiveElements(process_id);
264template <
int DisplacementDim,
typename ConstitutiveTraits>
265std::vector<std::string>
268 std::vector<std::reference_wrapper<MeshLib::Mesh>>
const& meshes)
270 INFO(
"TRM process initializeSubmeshOutput().");
271 const int process_id = 0;
272 std::vector<std::string> residuum_names{
"HeatFlowRate",
"MassFlowRate",
277 initializeAssemblyOnSubmeshes(process_id, meshes, residuum_names);
279 return residuum_names;
282template <
int DisplacementDim,
typename ConstitutiveTraits>
285 std::vector<GlobalVector*>
const& x_prev,
286 double const t,
double const dt,
287 const int process_id)
289 DBUG(
"PostTimestep ThermoRichardsMechanicsProcess.");
293 &LocalAssemblerIF::postTimestep, local_assemblers_,
298template <
int DisplacementDim,
typename ConstitutiveTraits>
301 std::vector<GlobalVector*>
const& x,
303 int const process_id)
305 DBUG(
"Compute the secondary variables for ThermoRichardsMechanicsProcess.");
309 &LocalAssemblerIF::computeSecondaryVariable, local_assemblers_,
314template <
int DisplacementDim,
typename ConstitutiveTraits>
316 DisplacementDim, ConstitutiveTraits>::getDOFTableForExtrapolatorData()
const
318 const bool manage_storage =
false;
319 return std::make_tuple(local_to_global_index_map_single_component_.get(),
323template <
int DisplacementDim,
typename ConstitutiveTraits>
325 DisplacementDim, ConstitutiveTraits>::getDOFTable(
const int )
328 return *_local_to_global_index_map;
338 2, ConstitutiveStressSaturation_StrainPressureTemperature::
339 ConstitutiveTraits<2>>;
341 3, ConstitutiveStressSaturation_StrainPressureTemperature::
342 ConstitutiveTraits<3>>;
void INFO(fmt::format_string< Args... > fmt, Args &&... args)
void DBUG(fmt::format_string< Args... > fmt, Args &&... args)
Global vector based on Eigen vector.
bool isAxiallySymmetric() const
std::vector< Element * > const & getElements() const
Get the element-vector for the mesh.
Properties & getProperties()
std::vector< std::size_t > const & getActiveElementIDs() const
std::unique_ptr< ProcessLib::AbstractJacobianAssembler > _jacobian_assembler
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
ThermoRichardsMechanicsProcess(std::string name, MeshLib::Mesh &mesh, std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&jacobian_assembler, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const ¶meters, 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)
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
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
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
Eigen::Matrix< double, kelvin_vector_dimensions(DisplacementDim), 1, Eigen::ColMajor > KelvinVectorType
std::vector< Node * > getBaseNodes(std::vector< Element * > const &elements)
@ BY_LOCATION
Ordering data by spatial location.
void addReflectedIntegrationPointWriters(ReflData const &reflection_data, std::vector< std::unique_ptr< MeshLib::IntegrationPointWriter > > &integration_point_writers, unsigned const integration_order, std::vector< std::unique_ptr< LocAsmIF > > const &local_assemblers)
SecondaryVariableFunctions makeExtrapolator(const unsigned num_components, NumLib::Extrapolator &extrapolator, LocalAssemblerCollection const &local_assemblers, typename NumLib::ExtrapolatableLocalAssemblerCollection< LocalAssemblerCollection >::IntegrationPointValuesMethod integration_point_values_method)
void setIPDataInitialConditions(std::vector< std::unique_ptr< MeshLib::IntegrationPointWriter > > const &_integration_point_writer, MeshLib::Properties const &mesh_properties, LocalAssemblersVector &local_assemblers)
static void executeSelectedMemberOnDereferenced(Method method, Container const &container, std::vector< std::size_t > const &active_container_ids, Args &&... args)
static void executeMemberOnDereferenced(Method method, Container const &container, Args &&... args)