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 std::unique_ptr<ProcessLib::SurfaceFluxData>&& surfaceflux)
38 :
Process(std::move(
name), mesh, std::move(jacobian_assembler), parameters,
39 integration_order, std::move(process_variables),
40 std::move(secondary_variables), use_monolithic_scheme),
41 _process_data(std::move(process_data)),
42 _surfaceflux(std::move(surfaceflux))
49 unsigned const integration_order)
55 const int process_id = 0;
60 ProcessLib::createLocalAssemblers<MonolithicHTFEM>(
67 ProcessLib::createLocalAssemblers<StaggeredHTFEM>(
81 std::vector<GlobalVector*>
const& x,
82 std::vector<GlobalVector*>
const& xdot,
86 std::vector<std::reference_wrapper<NumLib::LocalToGlobalIndexMap>>
90 DBUG(
"Assemble HTProcess.");
98 "Assemble the equations of heat transport process within "
104 "Assemble the equations of single phase fully saturated "
105 "fluid flow process within HTProcess.");
120 const double t,
double const dt, std::vector<GlobalVector*>
const& x,
121 std::vector<GlobalVector*>
const& xdot,
const double dxdot_dx,
125 DBUG(
"AssembleWithJacobian HTProcess.");
127 std::vector<std::reference_wrapper<NumLib::LocalToGlobalIndexMap>>
144 dxdot_dx, dx_dx, process_id, M, K, b, Jac);
148 int const process_id)
150 DBUG(
"Set the coupled term for the staggered scheme to local assembers.");
158 std::tuple<NumLib::LocalToGlobalIndexMap*, bool>
165 const bool manage_storage =
false;
171 std::vector<MeshLib::MeshSubset> all_mesh_subsets_single_component{
174 const bool manage_storage =
true;
176 std::move(all_mesh_subsets_single_component),
185 std::vector<GlobalVector*>
const& x)
const
188 std::vector<GlobalIndexType> indices_cache;
191 std::vector<std::vector<GlobalIndexType>> indices_of_all_coupled_processes{
192 x.size(), r_c_indices.rows};
203 int const process_id)
209 "The condition of process_id = 0 must be satisfied for monolithic "
210 "HTProcess, which is a single process.");
215 DBUG(
"This is the thermal part of the staggered HTProcess.");
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 std::vector< double > const & getIntPtDarcyVelocity(const double, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &) const =0
void setStaggeredCoupledSolutions(std::size_t const, CoupledSolutionsForStaggeredScheme *const coupling_term)
void postTimestepConcreteProcess(std::vector< GlobalVector * > const &x, const double t, const double delta_t, int const process_id) override
void setCoupledTermForTheStaggeredSchemeToLocalAssemblers(int const process_id) override
void initializeConcreteProcess(NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order) override
Process specific initialization called by initialize().
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
std::tuple< NumLib::LocalToGlobalIndexMap *, bool > getDOFTableForExtrapolatorData() const override
std::unique_ptr< ProcessLib::SurfaceFluxData > _surfaceflux
Eigen::Vector3d getFlux(std::size_t element_id, MathLib::Point3d const &p, double const t, std::vector< GlobalVector * > const &x) const override
HTProcessData _process_data
HTProcess(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, HTProcessData &&process_data, SecondaryVariableCollection &&secondary_variables, bool const use_monolithic_scheme, std::unique_ptr< ProcessLib::SurfaceFluxData > &&surfaceflux)
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
std::vector< std::unique_ptr< HTLocalAssemblerInterface > > _local_assemblers
unsigned getShapeFunctionOrder() const
std::vector< std::size_t > const & getActiveElementIDs() const
NumLib::Extrapolator & getExtrapolator() const
std::unique_ptr< MeshLib::MeshSubset const > _mesh_subset_all_nodes
SecondaryVariableCollection _secondary_variables
CoupledSolutionsForStaggeredScheme * _coupled_solutions
std::vector< std::reference_wrapper< ProcessVariable > > const & getProcessVariables(const int process_id) const
VectorMatrixAssembler _global_assembler
unsigned const _integration_order
std::unique_ptr< NumLib::LocalToGlobalIndexMap > _local_to_global_index_map
const bool _use_monolithic_scheme
Handles configuration of several secondary variables from the project file.
void addSecondaryVariable(std::string const &internal_name, SecondaryVariableFunctions &&fcts)
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)
NumLib::LocalToGlobalIndexMap::RowColumnIndices getRowColumnIndices(std::size_t const id, NumLib::LocalToGlobalIndexMap const &dof_table, std::vector< GlobalIndexType > &indices)
@ BY_LOCATION
Ordering data by spatial location.
std::vector< double > getCoupledLocalSolutions(std::vector< GlobalVector * > const &global_solutions, std::vector< std::vector< GlobalIndexType >> const &indices)
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)
int const hydraulic_process_id
int const heat_transport_process_id