28namespace RichardsMechanics
30template <
int DisplacementDim>
34 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&& jacobian_assembler,
35 std::vector<std::unique_ptr<ParameterLib::ParameterBase>>
const& parameters,
36 unsigned const integration_order,
37 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
41 bool const use_monolithic_scheme)
42 :
Process(std::move(name), mesh, std::move(jacobian_assembler), parameters,
43 integration_order, std::move(process_variables),
44 std::move(secondary_variables), use_monolithic_scheme),
47 process_data_(std::move(process_data))
56 DisplacementDim>(LocalAssemblerIF::getReflectionDataForOutput(),
57 _integration_point_writer, integration_order,
61template <
int DisplacementDim>
67template <
int DisplacementDim>
70 const int process_id)
const
74 if (_use_monolithic_scheme || process_id == 1)
76 auto const& l = *_local_to_global_index_map;
77 return {l.dofSizeWithoutGhosts(), l.dofSizeWithoutGhosts(),
78 &l.getGhostIndices(), &this->_sparsity_pattern};
82 auto const& l = *local_to_global_index_map_with_base_nodes_;
83 return {l.dofSizeWithoutGhosts(), l.dofSizeWithoutGhosts(),
84 &l.getGhostIndices(), &sparsity_pattern_with_linear_element_};
87template <
int DisplacementDim>
91 _mesh_subset_all_nodes =
92 std::make_unique<MeshLib::MeshSubset>(_mesh, _mesh.getNodes());
95 mesh_subset_base_nodes_ =
96 std::make_unique<MeshLib::MeshSubset>(_mesh, base_nodes_);
100 std::vector<MeshLib::MeshSubset> all_mesh_subsets_single_component{
101 *_mesh_subset_all_nodes};
102 local_to_global_index_map_single_component_ =
103 std::make_unique<NumLib::LocalToGlobalIndexMap>(
104 std::move(all_mesh_subsets_single_component),
108 if (_use_monolithic_scheme)
111 std::vector<MeshLib::MeshSubset> all_mesh_subsets{
112 *mesh_subset_base_nodes_};
115 const int monolithic_process_id = 0;
116 std::generate_n(std::back_inserter(all_mesh_subsets),
117 getProcessVariables(monolithic_process_id)[1]
119 .getNumberOfGlobalComponents(),
120 [&]() {
return *_mesh_subset_all_nodes; });
122 std::vector<int>
const vec_n_components{1, DisplacementDim};
123 _local_to_global_index_map =
124 std::make_unique<NumLib::LocalToGlobalIndexMap>(
125 std::move(all_mesh_subsets), vec_n_components,
127 assert(_local_to_global_index_map);
132 const int process_id = 1;
133 std::vector<MeshLib::MeshSubset> all_mesh_subsets;
134 std::generate_n(std::back_inserter(all_mesh_subsets),
135 getProcessVariables(process_id)[0]
137 .getNumberOfGlobalComponents(),
138 [&]() {
return *_mesh_subset_all_nodes; });
140 std::vector<int>
const vec_n_components{DisplacementDim};
141 _local_to_global_index_map =
142 std::make_unique<NumLib::LocalToGlobalIndexMap>(
143 std::move(all_mesh_subsets), vec_n_components,
148 std::vector<MeshLib::MeshSubset> all_mesh_subsets_base_nodes{
149 *mesh_subset_base_nodes_};
150 local_to_global_index_map_with_base_nodes_ =
151 std::make_unique<NumLib::LocalToGlobalIndexMap>(
152 std::move(all_mesh_subsets_base_nodes),
157 *local_to_global_index_map_with_base_nodes_, _mesh);
159 assert(_local_to_global_index_map);
160 assert(local_to_global_index_map_with_base_nodes_);
164template <
int DisplacementDim>
168 unsigned const integration_order)
177 LocalAssemblerIF::getReflectionDataForOutput(), _secondary_variables,
178 getExtrapolator(), local_assemblers_);
180 auto add_secondary_variable = [&](std::string
const& name,
181 int const num_components,
182 auto get_ip_values_function)
184 _secondary_variables.addSecondaryVariable(
188 std::move(get_ip_values_function)));
196 add_secondary_variable);
200 process_data_.solid_materials, local_assemblers_,
201 _integration_point_writer, integration_order);
215 DisplacementDim>::RowsAtCompileTime);
217 process_data_.pressure_interpolated =
228 *_local_to_global_index_map);
231template <
int DisplacementDim>
233 std::map<
int, std::shared_ptr<MaterialPropertyLib::Medium>>
const& media)
235 if (_use_monolithic_scheme)
237 const int monolithic_process_id = 0;
238 initializeProcessBoundaryConditionsAndSourceTerms(
239 *_local_to_global_index_map, monolithic_process_id, media);
245 const int hydraulic_process_id = 0;
246 initializeProcessBoundaryConditionsAndSourceTerms(
247 *local_to_global_index_map_with_base_nodes_, hydraulic_process_id,
251 const int mechanical_process_id = 1;
252 initializeProcessBoundaryConditionsAndSourceTerms(
253 *_local_to_global_index_map, mechanical_process_id, media);
256template <
int DisplacementDim>
260 int const process_id)
267 DBUG(
"SetInitialConditions RichardsMechanicsProcess.");
270 &LocalAssemblerIF::setInitialConditions, local_assemblers_,
271 getActiveElementIDs(), getDOFTables(x.size()), x, t, process_id);
274template <
int DisplacementDim>
276 const double t,
double const dt, std::vector<GlobalVector*>
const& x,
277 std::vector<GlobalVector*>
const& x_prev,
int const process_id,
280 DBUG(
"Assemble the equations for RichardsMechanics");
283 t, dt, x, x_prev, process_id, M, K, b);
286template <
int DisplacementDim>
289 const double t,
double const dt, std::vector<GlobalVector*>
const& x,
290 std::vector<GlobalVector*>
const& x_prev,
int const process_id,
294 if (_use_monolithic_scheme)
297 "Assemble the Jacobian of RichardsMechanics for the monolithic"
306 "Assemble the Jacobian equations of liquid fluid process in "
307 "RichardsMechanics for the staggered scheme.");
312 "Assemble the Jacobian equations of mechanical process in "
313 "RichardsMechanics for the staggered scheme.");
317 auto const dof_tables = getDOFTables(x.size());
319 assembleWithJacobian(t, dt, x, x_prev, process_id, b, Jac);
321 auto copyRhs = [&](
int const variable_id,
auto& output_vector)
323 if (_use_monolithic_scheme)
325 transformVariableFromGlobalVector(b, variable_id, *dof_tables[0],
327 std::negate<double>());
331 transformVariableFromGlobalVector(b, 0, *dof_tables[process_id],
333 std::negate<double>());
336 if (_use_monolithic_scheme || process_id == 0)
338 copyRhs(0, *hydraulic_flow_);
340 if (_use_monolithic_scheme || process_id == 1)
342 copyRhs(1, *nodal_forces_);
346template <
int DisplacementDim>
348 std::vector<GlobalVector*>
const& x,
double const t,
double const dt,
349 const int process_id)
351 DBUG(
"PreTimestep RichardsMechanicsProcess.");
354 &LocalAssemblerIF::preTimestep, local_assemblers_,
355 getActiveElementIDs(), *_local_to_global_index_map, *x[process_id], t,
362template <
int DisplacementDim>
364 std::vector<GlobalVector*>
const& x,
365 std::vector<GlobalVector*>
const& x_prev,
double const t,
double const dt,
366 const int process_id)
368 if (hasMechanicalProcess(process_id))
370 DBUG(
"PostTimestep RichardsMechanicsProcess.");
373 &LocalAssemblerIF::postTimestep, local_assemblers_,
374 getActiveElementIDs(), getDOFTables(x.size()), x, x_prev, t, dt,
379template <
int DisplacementDim>
380std::vector<std::vector<std::string>>
382 std::vector<std::reference_wrapper<MeshLib::Mesh>>
const& meshes)
384 INFO(
"RichardsMechanics process initializeSubmeshOutput().");
385 std::vector<std::vector<std::string>> residuum_names{
386 {
"MassFlowRate",
"NodalForces"}};
389 initializeAssemblyOnSubmeshes(meshes, residuum_names);
391 return residuum_names;
394template <
int DisplacementDim>
397 std::vector<GlobalVector*>
const& x,
399 int const process_id)
406 DBUG(
"Compute the secondary variables for RichardsMechanicsProcess.");
409 &LocalAssemblerIF::computeSecondaryVariable, local_assemblers_,
410 getActiveElementIDs(), getDOFTables(x.size()), t, dt, x, x_prev,
414template <
int DisplacementDim>
416 DisplacementDim>::getDOFTableForExtrapolatorData()
const
418 const bool manage_storage =
false;
419 return std::make_tuple(local_to_global_index_map_single_component_.get(),
416 DisplacementDim>::getDOFTableForExtrapolatorData()
const {
…}
423template <
int DisplacementDim>
426 const int process_id)
const
428 if (hasMechanicalProcess(process_id))
430 return *_local_to_global_index_map;
434 return *local_to_global_index_map_with_base_nodes_;
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::unique_ptr< ProcessLib::AbstractJacobianAssembler > _jacobian_assembler
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 &x, double const t, double const dt, const int process_id) override
MathLib::MatrixSpecifications getMatrixSpecifications(const int process_id) const override
void constructDofTable() override
void initializeConcreteProcess(NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order) override
Process specific initialization called by initialize().
void computeSecondaryVariableConcrete(double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_prev, int const process_id) override
RichardsMechanicsProcess(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, RichardsMechanicsProcessData< DisplacementDim > &&process_data, SecondaryVariableCollection &&secondary_variables, bool const use_monolithic_scheme)
bool isLinear() const override
void setInitialConditionsConcreteProcess(std::vector< GlobalVector * > &x, double const t, int const process_id) override
void assembleWithJacobianConcreteProcess(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, GlobalVector &b, GlobalMatrix &Jac) override
void initializeBoundaryConditions(std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media) override
void postTimestepConcreteProcess(std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, const int process_id) override
std::vector< std::vector< std::string > > initializeAssemblyOnSubmeshes(std::vector< std::reference_wrapper< MeshLib::Mesh > > const &meshes) override
NumLib::LocalToGlobalIndexMap const & getDOFTable(const int process_id) const override
Handles configuration of several secondary variables from the project file.
Eigen::Matrix< double, kelvin_vector_dimensions(DisplacementDim), 1, Eigen::ColMajor > KelvinVectorType
PropertyVector< T > * getOrCreateMeshProperty(Mesh &mesh, std::string const &property_name, MeshItemType const item_type, int const number_of_components)
std::vector< Node * > getBaseNodes(std::vector< Element * > const &elements)
@ BY_LOCATION
Ordering data by spatial location.
GlobalSparsityPattern computeSparsityPattern(LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh)
Computes a sparsity pattern for the given inputs.
void addReflectedSecondaryVariables(ReflData const &reflection_data, SecondaryVariableCollection &secondary_variables, NumLib::Extrapolator &extrapolator, std::vector< std::unique_ptr< LocAsmIF > > const &local_assemblers)
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)
void createLocalAssemblersHM(std::vector< MeshLib::Element * > const &mesh_elements, NumLib::LocalToGlobalIndexMap const &dof_table, std::vector< std::unique_ptr< LocalAssemblerInterface > > &local_assemblers, ProviderOrOrder const &provider_or_order, ExtraCtorArgs &&... extra_ctor_args)
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)