OGS
ThermoHydroMechanicsProcess.h
Go to the documentation of this file.
1
11#pragma once
12
14#include "ProcessLib/Process.h"
16
17namespace ProcessLib
18{
19namespace ThermoHydroMechanics
20{
26template <int DisplacementDim>
28{
29public:
31 std::string name,
32 MeshLib::Mesh& mesh,
33 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&&
34 jacobian_assembler,
35 std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const&
36 parameters,
37 unsigned const integration_order,
38 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
39 process_variables,
41 SecondaryVariableCollection&& secondary_variables,
42 bool const use_monolithic_scheme);
43
46
47 bool isLinear() const override;
49
63 const int process_id) const override;
64
65private:
66 void constructDofTable() override;
67
69 NumLib::LocalToGlobalIndexMap const& dof_table,
70 MeshLib::Mesh const& mesh,
71 unsigned const integration_order) override;
72
74 std::map<int, std::shared_ptr<MaterialPropertyLib::Medium>> const&
75 media) override;
76
77 void setInitialConditionsConcreteProcess(std::vector<GlobalVector*>& x,
78 double const t,
79 int const process_id) override;
80
81 void assembleConcreteProcess(const double t, double const dt,
82 std::vector<GlobalVector*> const& x,
83 std::vector<GlobalVector*> const& x_prev,
84 int const process_id, GlobalMatrix& M,
85 GlobalMatrix& K, GlobalVector& b) override;
86
88 const double t, double const dt, std::vector<GlobalVector*> const& x,
89 std::vector<GlobalVector*> const& x_prev, int const process_id,
90 GlobalVector& b, GlobalMatrix& Jac) override;
91
92 void preTimestepConcreteProcess(std::vector<GlobalVector*> const& x,
93 double const t, double const dt,
94 const int process_id) override;
95
96 void postTimestepConcreteProcess(std::vector<GlobalVector*> const& x,
97 std::vector<GlobalVector*> const& x_prev,
98 const double t, const double dt,
99 int const process_id) override;
100
102 const int process_id) const override;
103
104private:
105 std::vector<MeshLib::Node*> _base_nodes;
106 std::unique_ptr<MeshLib::MeshSubset const> _mesh_subset_base_nodes;
108
109 std::vector<std::unique_ptr<LocalAssemblerInterface<DisplacementDim>>>
111
112 std::unique_ptr<NumLib::LocalToGlobalIndexMap>
114
117 std::unique_ptr<NumLib::LocalToGlobalIndexMap>
119
123
124 void computeSecondaryVariableConcrete(double const t, double const dt,
125 std::vector<GlobalVector*> const& x,
126 GlobalVector const& x_prev,
127 const int process_id) override;
131 std::tuple<NumLib::LocalToGlobalIndexMap*, bool>
132 getDOFTableForExtrapolatorData() const override;
133
137 bool hasMechanicalProcess(int const process_id) const
138 {
139 return _use_monolithic_scheme || process_id == 2;
140 }
141
145};
146
147extern template class ThermoHydroMechanicsProcess<2>;
148extern template class ThermoHydroMechanicsProcess<3>;
149
150} // namespace ThermoHydroMechanics
151} // namespace ProcessLib
MathLib::SparsityPattern< GlobalIndexType > GlobalSparsityPattern
Global vector based on Eigen vector.
Definition EigenVector.h:25
std::string const name
Definition Process.h:353
const bool _use_monolithic_scheme
Definition Process.h:370
Handles configuration of several secondary variables from the project file.
std::unique_ptr< NumLib::LocalToGlobalIndexMap > _local_to_global_index_map_with_base_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
NumLib::LocalToGlobalIndexMap const & getDOFTable(const int process_id) const override
void initializeConcreteProcess(NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order) override
Process specific initialization called by initialize().
void preTimestepConcreteProcess(std::vector< GlobalVector * > const &x, double const t, double const dt, const int process_id) override
std::unique_ptr< NumLib::LocalToGlobalIndexMap > _local_to_global_index_map_single_component
ThermoHydroMechanicsProcess(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, ThermoHydroMechanicsProcessData< DisplacementDim > &&process_data, SecondaryVariableCollection &&secondary_variables, bool const use_monolithic_scheme)
void initializeBoundaryConditions(std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media) override
std::vector< std::unique_ptr< LocalAssemblerInterface< DisplacementDim > > > _local_assemblers
void computeSecondaryVariableConcrete(double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_prev, const int process_id) override
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
std::tuple< NumLib::LocalToGlobalIndexMap *, bool > getDOFTableForExtrapolatorData() const override
void postTimestepConcreteProcess(std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, const double t, const double dt, int const process_id) override
MathLib::MatrixSpecifications getMatrixSpecifications(const int process_id) const override
void setInitialConditionsConcreteProcess(std::vector< GlobalVector * > &x, double const t, int const process_id) override
ThermoHydroMechanicsProcessData< DisplacementDim > _process_data