23 namespace ThermoMechanicalPhaseField
25 template <
int DisplacementDim>
30 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&&
32 std::vector<std::unique_ptr<ParameterLib::ParameterBase>>
const&
34 unsigned const integration_order,
35 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
39 int const mechanics_related_process_id,
40 int const phase_field_process_id,
41 int const heat_conduction_process_id)
42 :
Process(std::move(
name), mesh, std::move(jacobian_assembler), parameters,
43 integration_order, std::move(process_variables),
44 std::move(secondary_variables), false),
45 _process_data(std::move(process_data)),
46 _mechanics_related_process_id(mechanics_related_process_id),
47 _phase_field_process_id(phase_field_process_id),
48 _heat_conduction_process_id(heat_conduction_process_id)
52 template <
int DisplacementDim>
58 template <
int DisplacementDim>
61 const int process_id)
const
63 if (process_id == _mechanics_related_process_id)
65 auto const& l = *_local_to_global_index_map;
66 return {l.dofSizeWithoutGhosts(), l.dofSizeWithoutGhosts(),
67 &l.getGhostIndices(), &this->_sparsity_pattern};
71 auto const& l = *_local_to_global_index_map_single_component;
72 return {l.dofSizeWithoutGhosts(), l.dofSizeWithoutGhosts(),
73 &l.getGhostIndices(), &_sparsity_pattern_with_single_component};
76 template <
int DisplacementDim>
79 const int process_id)
const
81 if (process_id == _mechanics_related_process_id)
83 return *_local_to_global_index_map;
87 return *_local_to_global_index_map_single_component;
90 template <
int DisplacementDim>
93 const int process_id)
const
95 if (process_id == _mechanics_related_process_id)
97 return *_local_to_global_index_map;
101 return *_local_to_global_index_map_single_component;
104 template <
int DisplacementDim>
108 constructDofTableOfSpecifiedProcessStaggeredScheme(
109 _mechanics_related_process_id);
113 std::vector<MeshLib::MeshSubset> all_mesh_subsets_single_component{
114 *_mesh_subset_all_nodes};
115 _local_to_global_index_map_single_component =
116 std::make_unique<NumLib::LocalToGlobalIndexMap>(
117 std::move(all_mesh_subsets_single_component),
121 assert(_local_to_global_index_map_single_component);
125 *_local_to_global_index_map_single_component, _mesh);
128 template <
int DisplacementDim>
132 unsigned const integration_order)
138 _mechanics_related_process_id, _phase_field_process_id,
139 _heat_conduction_process_id);
141 _secondary_variables.addSecondaryVariable(
145 DisplacementDim>::RowsAtCompileTime,
146 getExtrapolator(), _local_assemblers,
149 _secondary_variables.addSecondaryVariable(
152 DisplacementDim>::RowsAtCompileTime,
153 getExtrapolator(), _local_assemblers,
157 _secondary_variables.addSecondaryVariable(
167 *_local_to_global_index_map);
170 template <
int DisplacementDim>
172 DisplacementDim>::initializeBoundaryConditions()
176 initializeProcessBoundaryConditionsAndSourceTerms(
177 getDOFTableByProcessID(_mechanics_related_process_id),
178 _mechanics_related_process_id);
180 initializeProcessBoundaryConditionsAndSourceTerms(
181 getDOFTableByProcessID(_phase_field_process_id),
182 _phase_field_process_id);
184 initializeProcessBoundaryConditionsAndSourceTerms(
185 getDOFTableByProcessID(_heat_conduction_process_id),
186 _heat_conduction_process_id);
189 template <
int DisplacementDim>
192 std::vector<GlobalVector*>
const& x,
193 std::vector<GlobalVector*>
const& xdot,
197 DBUG(
"Assemble the equations for ThermoMechanicalPhaseFieldProcess.");
199 std::vector<std::reference_wrapper<NumLib::LocalToGlobalIndexMap>>
200 dof_table = {std::ref(*_local_to_global_index_map)};
210 template <
int DisplacementDim>
213 const double t,
double const dt, std::vector<GlobalVector*>
const& x,
214 std::vector<GlobalVector*>
const& xdot,
const double dxdot_dx,
215 const double dx_dx,
int const process_id,
GlobalMatrix& M,
218 std::vector<std::reference_wrapper<NumLib::LocalToGlobalIndexMap>>
221 if (process_id == _mechanics_related_process_id)
224 "Assemble the Jacobian equations of "
225 "temperature-deformation in ThermoMechanicalPhaseFieldProcess for "
226 "the staggered scheme.");
229 if (process_id == _phase_field_process_id)
232 "Assemble the Jacobian equations ofphase field in "
233 "ThermoMechanicalPhaseFieldProcess for the staggered scheme.");
238 "Assemble the Jacobian equations of heat conduction in "
239 "ThermoMechanicalPhaseFieldProcess for the staggered scheme.");
241 dof_tables.emplace_back(
242 getDOFTableByProcessID(_heat_conduction_process_id));
243 dof_tables.emplace_back(
244 getDOFTableByProcessID(_mechanics_related_process_id));
245 dof_tables.emplace_back(getDOFTableByProcessID(_phase_field_process_id));
252 dxdot_dx, dx_dx, process_id, M, K, b, Jac);
255 template <
int DisplacementDim>
260 const int process_id)
262 DBUG(
"PreTimestep ThermoMechanicalPhaseFieldProcess.");
264 if (process_id != _mechanics_related_process_id)
274 *x[process_id], t, dt);
277 template <
int DisplacementDim>
282 int const process_id)
289 DBUG(
"PostTimestep ThermoMechanicalPhaseFieldProcess.");
290 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_tables;
291 auto const n_processes = x.size();
292 dof_tables.reserve(n_processes);
293 for (std::size_t process_id = 0; process_id < n_processes; ++process_id)
295 dof_tables.push_back(&getDOFTable(process_id));
304 template <
int DisplacementDim>
308 double const dt,
const int process_id)
310 if (process_id != _mechanics_related_process_id)
315 DBUG(
"PostNonLinearSolver ThermoMechanicalPhaseFieldProcess.");
317 const bool use_monolithic_scheme =
false;
323 use_monolithic_scheme, process_id);
void DBUG(char const *fmt, Args const &... args)
Global vector based on Eigen vector.
bool isAxiallySymmetric() const
unsigned getDimension() const
Returns the dimension of the mesh (determined by the maximum dimension over all elements).
std::vector< Element * > const & getElements() const
Get the element-vector for the mesh.
virtual void postTimestep(std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, double const t, double const dt)
virtual void preTimestep(std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, double const t, double const delta_t)
virtual void initialize(std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table)
void postNonLinearSolver(std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, GlobalVector const &xdot, double const t, double const dt, bool const use_monolithic_scheme, int const process_id)
std::vector< std::size_t > const & getActiveElementIDs() const
Handles configuration of several secondary variables from the project file.
A class to simulate thermo-mechanical fracturing process using phase-field approach in solids describ...
ThermoMechanicalPhaseFieldProcess(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, ThermoMechanicalPhaseFieldProcessData< DisplacementDim > &&process_data, SecondaryVariableCollection &&secondary_variables, int const mechanics_related_process_id, int const phase_field_process_id, int const heat_conduction_process_id)
NumLib::LocalToGlobalIndexMap const & getDOFTable(const int process_id) const override
void assembleWithJacobianConcreteProcess(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, const double dxdot_dx, const double dx_dx, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b, GlobalMatrix &Jac) override
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
void postNonLinearSolverConcreteProcess(GlobalVector const &x, GlobalVector const &xdot, const double t, double const dt, int const process_id) 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 postTimestepConcreteProcess(std::vector< GlobalVector * > const &x, double const t, double const dt, int const process_id) override
NumLib::LocalToGlobalIndexMap & getDOFTableByProcessID(const int process_id) const
void constructDofTable() override
bool isLinear() 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 assembleWithJacobian(std::size_t const mesh_item_id, LocalAssemblerInterface &local_assembler, std::vector< std::reference_wrapper< NumLib::LocalToGlobalIndexMap >> const &dof_tables, const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, const double dxdot_dx, const double dx_dx, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b, GlobalMatrix &Jac)
void assemble(std::size_t const mesh_item_id, LocalAssemblerInterface &local_assembler, std::vector< std::reference_wrapper< NumLib::LocalToGlobalIndexMap >> const &dof_tables, double const t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b)
Eigen::Matrix< double, kelvin_vector_dimensions(DisplacementDim), 1, Eigen::ColMajor > KelvinVectorType
@ 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.
SecondaryVariableFunctions makeExtrapolator(const unsigned num_components, NumLib::Extrapolator &extrapolator, LocalAssemblerCollection const &local_assemblers, typename NumLib::ExtrapolatableLocalAssemblerCollection< LocalAssemblerCollection >::IntegrationPointValuesMethod integration_point_values_method)
static void executeSelectedMemberOnDereferenced(Method method, Container const &container, std::vector< std::size_t > const &active_container_ids, Args &&... args)
static void executeSelectedMemberDereferenced(Object &object, 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)
virtual std::vector< double > const & getIntPtSigma(const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const =0