OGS
HydroMechanicsProcess.h
Go to the documentation of this file.
1
11#pragma once
12
15#include "ProcessLib/Process.h"
16
17namespace ProcessLib
18{
19namespace LIE
20{
21namespace HydroMechanics
22{
23class HydroMechanicsLocalAssemblerInterface;
24
25template <int GlobalDim>
26class HydroMechanicsProcess final : public Process
27{
28 static_assert(GlobalDim == 2 || GlobalDim == 3,
29 "Currently LIE::HydroMechanicsProcess "
30 "supports only 2D or 3D.");
31
32public:
34 std::string name,
35 MeshLib::Mesh& mesh,
36 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&&
37 jacobian_assembler,
38 std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const&
39 parameters,
40 unsigned const integration_order,
41 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
42 process_variables,
44 SecondaryVariableCollection&& secondary_variables,
45 bool const use_monolithic_scheme);
46
49 bool isLinear() const override;
51
52 void postTimestepConcreteProcess(std::vector<GlobalVector*> const& x,
53 std::vector<GlobalVector*> const& x_prev,
54 double const t, double const dt,
55 int const process_id) override;
56
57private:
59
60 void constructDofTable() override;
61
63 NumLib::LocalToGlobalIndexMap const& dof_table,
64 MeshLib::Mesh const& mesh,
65 unsigned const integration_order) override;
66
67 void assembleConcreteProcess(const double t, double const dt,
68 std::vector<GlobalVector*> const& x,
69 std::vector<GlobalVector*> const& x_prev,
70 int const process_id, GlobalMatrix& M,
71 GlobalMatrix& K, GlobalVector& b) override;
72
74 const double t, double const dt, std::vector<GlobalVector*> const& x,
75 std::vector<GlobalVector*> const& x_prev, int const process_id,
76 GlobalVector& b, GlobalMatrix& Jac) override;
77 void preTimestepConcreteProcess(std::vector<GlobalVector*> const& x,
78 double const t, double const dt,
79 int const process_id) override;
80
81private:
83
84 std::vector<std::unique_ptr<LocalAssemblerInterface>> _local_assemblers;
85
86 std::vector<MeshLib::Element*> _vec_matrix_elements;
87 std::vector<MeshLib::Element*> _vec_fracture_elements;
88 std::vector<MeshLib::Element*> _vec_fracture_matrix_elements;
89 std::vector<MeshLib::Node*> _vec_fracture_nodes;
90
91 std::unique_ptr<MeshLib::MeshSubset const> _mesh_subset_fracture_nodes;
92 std::unique_ptr<MeshLib::MeshSubset const> _mesh_subset_matrix_nodes;
93
94 std::vector<MeshLib::Node*> _mesh_nodes_p;
95 std::unique_ptr<MeshLib::MeshSubset const> _mesh_subset_nodes_p;
96};
97
98} // namespace HydroMechanics
99} // namespace LIE
100} // namespace ProcessLib
Global vector based on Eigen vector.
Definition EigenVector.h:25
std::unique_ptr< MeshLib::MeshSubset const > _mesh_subset_matrix_nodes
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 preTimestepConcreteProcess(std::vector< GlobalVector * > const &x, double const t, double const dt, int const process_id) override
std::vector< std::unique_ptr< LocalAssemblerInterface > > _local_assemblers
std::unique_ptr< MeshLib::MeshSubset const > _mesh_subset_fracture_nodes
std::vector< MeshLib::Element * > _vec_fracture_matrix_elements
HydroMechanicsProcess(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, HydroMechanicsProcessData< GlobalDim > &&process_data, SecondaryVariableCollection &&secondary_variables, bool const use_monolithic_scheme)
std::unique_ptr< MeshLib::MeshSubset const > _mesh_subset_nodes_p
void initializeConcreteProcess(NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order) override
Process specific initialization called by initialize().
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 postTimestepConcreteProcess(std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id) override
std::string const name
Definition Process.h:362
Handles configuration of several secondary variables from the project file.