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