OGS
ThermoRichardsMechanicsProcess.h
Go to the documentation of this file.
1 
11 #pragma once
12 
14 #include "ProcessLib/Process.h"
16 
17 namespace ProcessLib
18 {
19 namespace ThermoRichardsMechanics
20 {
108 template <int DisplacementDim>
110 {
111 public:
113  std::string name,
114  MeshLib::Mesh& mesh,
115  std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&&
116  jacobian_assembler,
117  std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const&
118  parameters,
119  unsigned const integration_order,
120  std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
121  process_variables,
123  SecondaryVariableCollection&& secondary_variables,
124  bool const use_monolithic_scheme);
125 
128 
129  bool isLinear() const override;
131 
145  const int process_id) const override;
146 
147 private:
149 
150  void constructDofTable() override;
151 
153  NumLib::LocalToGlobalIndexMap const& dof_table,
154  MeshLib::Mesh const& mesh,
155  unsigned const integration_order) override;
156 
157  void initializeBoundaryConditions() override;
158 
159  void setInitialConditionsConcreteProcess(std::vector<GlobalVector*>& x,
160  double const t,
161  int const process_id) override;
162 
163  void assembleConcreteProcess(const double t, double const dt,
164  std::vector<GlobalVector*> const& x,
165  std::vector<GlobalVector*> const& xdot,
166  int const process_id, GlobalMatrix& M,
167  GlobalMatrix& K, GlobalVector& b) override;
168 
170  const double t, double const dt, std::vector<GlobalVector*> const& x,
171  std::vector<GlobalVector*> const& xdot, int const process_id,
173  GlobalMatrix& Jac) override;
174 
175  void postTimestepConcreteProcess(std::vector<GlobalVector*> const& x,
176  double const t, double const dt,
177  const int process_id) override;
178 
180  const int process_id) const override;
181 
182 private:
183  std::vector<MeshLib::Node*> base_nodes_;
184  std::unique_ptr<MeshLib::MeshSubset const> mesh_subset_base_nodes_;
186 
187  std::vector<std::unique_ptr<LocalAssemblerIF>> local_assemblers_;
188 
189  std::unique_ptr<NumLib::LocalToGlobalIndexMap>
191 
194  std::unique_ptr<NumLib::LocalToGlobalIndexMap>
196 
200 
201  void computeSecondaryVariableConcrete(double const t, double const dt,
202  std::vector<GlobalVector*> const& x,
203  GlobalVector const& x_dot,
204  int const process_id) override;
208  std::tuple<NumLib::LocalToGlobalIndexMap*, bool>
209  getDOFTableForExtrapolatorData() const override;
210 
211  std::vector<NumLib::LocalToGlobalIndexMap const*> getDOFTables(
212  const int number_of_processes) const;
213 
217  bool hasMechanicalProcess(int const /*process_id*/) const { return true; }
218 
222 };
223 
224 extern template class ThermoRichardsMechanicsProcess<2>;
225 extern template class ThermoRichardsMechanicsProcess<3>;
226 
227 } // namespace ThermoRichardsMechanics
228 } // namespace ProcessLib
MathLib::SparsityPattern< GlobalIndexType > GlobalSparsityPattern
Global vector based on Eigen vector.
Definition: EigenVector.h:28
std::string const name
Definition: Process.h:327
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...
std::tuple< NumLib::LocalToGlobalIndexMap *, bool > getDOFTableForExtrapolatorData() const override
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
void computeSecondaryVariableConcrete(double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_dot, int const process_id) override
std::unique_ptr< NumLib::LocalToGlobalIndexMap > local_to_global_index_map_with_base_nodes_
ThermoRichardsMechanicsProcessData< DisplacementDim > process_data_
ThermoRichardsMechanicsProcess(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, ThermoRichardsMechanicsProcessData< DisplacementDim > &&process_data, SecondaryVariableCollection &&secondary_variables, bool const use_monolithic_scheme)
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
MathLib::MatrixSpecifications getMatrixSpecifications(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().
std::vector< NumLib::LocalToGlobalIndexMap const * > getDOFTables(const int number_of_processes) const
void postTimestepConcreteProcess(std::vector< GlobalVector * > const &x, double const t, double const dt, const int process_id) override
void setInitialConditionsConcreteProcess(std::vector< GlobalVector * > &x, double const t, int const process_id) override
std::unique_ptr< NumLib::LocalToGlobalIndexMap > local_to_global_index_map_single_component_
NumLib::LocalToGlobalIndexMap const & getDOFTable(const int process_id) const override
static const double t