25template <
int DisplacementDim>
29 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&& jacobian_assembler,
30 std::vector<std::unique_ptr<ParameterLib::ParameterBase>>
const& parameters,
31 unsigned const integration_order,
32 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
36 bool const use_monolithic_scheme)
37 :
Process(std::move(
name), mesh, std::move(jacobian_assembler), parameters,
38 integration_order, std::move(process_variables),
39 std::move(secondary_variables), use_monolithic_scheme),
42 if (use_monolithic_scheme)
45 "Monolithic scheme is not implemented for the HMPhaseField "
53template <
int DisplacementDim>
59template <
int DisplacementDim>
62 const int process_id)
const
65 if (process_id ==
_process_data._mechanics_related_process_id)
68 return {l.dofSizeWithoutGhosts(), l.dofSizeWithoutGhosts(),
74 return {l.dofSizeWithoutGhosts(), l.dofSizeWithoutGhosts(),
78template <
int DisplacementDim>
83 if (process_id ==
_process_data._mechanics_related_process_id)
92template <
int DisplacementDim>
101 std::vector<MeshLib::MeshSubset> all_mesh_subsets_single_component{
104 std::make_unique<NumLib::LocalToGlobalIndexMap>(
105 std::move(all_mesh_subsets_single_component),
116template <
int DisplacementDim>
120 unsigned const integration_order)
131 DisplacementDim>::RowsAtCompileTime,
138 DisplacementDim>::RowsAtCompileTime,
161template <
int DisplacementDim>
163 std::map<
int, std::shared_ptr<MaterialPropertyLib::Medium>>
const& media)
180template <
int DisplacementDim>
182 const double ,
double const ,
183 std::vector<GlobalVector*>
const& ,
184 std::vector<GlobalVector*>
const& ,
int const ,
188 "HMPhaseFieldLocalAssembler: assembly for the Picard non-linear solver "
189 "is not implemented.");
192template <
int DisplacementDim>
194 const double t,
double const dt, std::vector<GlobalVector*>
const& x,
195 std::vector<GlobalVector*>
const& x_prev,
int const process_id,
198 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_tables;
204 "Assemble the Jacobian equations of phase field in "
205 "HMPhaseFieldProcess for the staggered scheme.");
210 "Assemble the Jacobian equations of pressure in "
211 "HMPhaseFieldProcess for the staggered scheme.");
216 "Assemble the Jacobian equations of deformation in "
217 "HMPhaseFieldProcess for the staggered scheme.");
228 process_id, &b, &Jac);
230 if (process_id ==
_process_data._mechanics_related_process_id)
238template <
int DisplacementDim>
240 std::vector<GlobalVector*>
const& x,
double const t,
double const dt,
241 const int process_id)
243 DBUG(
"PreTimestep HMPhaseFieldProcess {}.", process_id);
247 DBUG(
"Store the value of phase field at previous time step.");
258template <
int DisplacementDim>
260 std::vector<GlobalVector*>
const& x,
261 std::vector<GlobalVector*>
const& x_prev,
const double t,
const double dt,
262 int const process_id)
266 DBUG(
"PostTimestep HMPhaseFieldProcess.");
272 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_tables;
274 dof_tables.emplace_back(
276 dof_tables.emplace_back(
297template <
int DisplacementDim>
299 std::vector<GlobalVector*>
const& x,
300 std::vector<GlobalVector*>
const& x_prev,
const double t,
double const dt,
301 const int process_id)
303 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_tables;
311 INFO(
"Update fracture width and porous properties");
319 INFO(
"PostNonLinearSolver for Hydro Process");
328template <
int DisplacementDim>
MathLib::EigenMatrix GlobalMatrix
MathLib::EigenVector GlobalVector
GlobalMatrix::IndexType GlobalIndexType
void INFO(fmt::format_string< Args... > fmt, Args &&... args)
void DBUG(fmt::format_string< Args... > fmt, Args &&... args)
void copyValues(std::vector< double > &u) const
void set(IndexType rowId, double v)
set entry
bool isAxiallySymmetric() const
std::vector< Element * > const & getElements() const
Get the element-vector for the mesh.
std::unique_ptr< GlobalVector > _x_previous_timestep
bool isLinear() const override
void postTimestepConcreteProcess(std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &, const double t, const double dt, int const process_id) override
void constructDofTable() override
HMPhaseFieldProcess(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, HMPhaseFieldProcessData< DisplacementDim > &&process_data, SecondaryVariableCollection &&secondary_variables, bool const use_monolithic_scheme)
std::unique_ptr< NumLib::LocalToGlobalIndexMap > _local_to_global_index_map_single_component
void initializeBoundaryConditions(std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media) override
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 initializeConcreteProcess(NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order) override
Process specific initialization called by initialize().
MathLib::MatrixSpecifications getMatrixSpecifications(const int process_id) const override
NumLib::LocalToGlobalIndexMap const & getDOFTable(const int process_id) const override
void postNonLinearSolverConcreteProcess(std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, const double t, double const dt, int const process_id) override
HMPhaseFieldProcessData< DisplacementDim > _process_data
std::vector< std::unique_ptr< LocalAssemblerInterface > > _local_assemblers
GlobalSparsityPattern _sparsity_pattern_with_single_component
void updateConstraints(GlobalVector &lower, GlobalVector &upper, int const process_id) override
void preTimestepConcreteProcess(std::vector< GlobalVector * > const &x, double const t, double const dt, 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
MeshLib::PropertyVector< double > * _nodal_forces
void postNonLinearSolver(std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
virtual void postTimestep(std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
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)
std::vector< NumLib::LocalToGlobalIndexMap const * > getDOFTables(int const number_of_processes) const
void constructDofTableOfSpecifiedProcessStaggeredScheme(const int specified_process_id)
void initializeProcessBoundaryConditionsAndSourceTerms(const NumLib::LocalToGlobalIndexMap &dof_table, const int process_id, std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)
Process(std::string name_, MeshLib::Mesh &mesh, std::unique_ptr< 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, SecondaryVariableCollection &&secondary_variables, const bool use_monolithic_scheme=true)
std::unique_ptr< MeshLib::MeshSubset const > _mesh_subset_all_nodes
std::vector< std::size_t > const & getActiveElementIDs() const
SecondaryVariableCollection _secondary_variables
VectorMatrixAssembler _global_assembler
std::unique_ptr< NumLib::LocalToGlobalIndexMap > _local_to_global_index_map
NumLib::Extrapolator & getExtrapolator() const
GlobalSparsityPattern _sparsity_pattern
Handles configuration of several secondary variables from the project file.
void assembleWithJacobian(std::size_t const mesh_item_id, LocalAssemblerInterface &local_assembler, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, 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)
Eigen::Matrix< double, kelvin_vector_dimensions(DisplacementDim), 1, Eigen::ColMajor > KelvinVectorType
void copy(PETScVector const &x, PETScVector &y)
void setLocalAccessibleVector(PETScVector const &x)
PropertyVector< T > * getOrCreateMeshProperty(Mesh &mesh, std::string const &property_name, MeshItemType const item_type, int const number_of_components)
@ 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.
void showEnergyAndWork(const double t, double &_elastic_energy, double &_surface_energy, double &_pressure_work)
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
virtual void approximateFractureWidth(std::size_t mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, double const t, double const dt)=0
virtual void computeEnergy(std::size_t mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, double const t, double &elastic_energy, double &surface_energy, double &pressure_work)=0
virtual std::vector< double > const & getIntPtWidth(const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const =0
virtual std::vector< double > const & getIntPtEpsilon(const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const =0