34 std::array
const required_medium_properties = {
40 std::array
const required_liquid_properties = {
44 std::array<MaterialPropertyLib::PropertyType, 0>
const
45 required_solid_properties{};
47 std::array<MaterialPropertyLib::PropertyType, 0>
const
48 required_gas_properties{};
51 mesh, media_map, required_medium_properties, required_solid_properties,
52 required_liquid_properties, required_gas_properties);
56 std::string
const& name,
58 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&& jacobian_assembler,
59 std::vector<ProcessVariable>
const& variables,
60 std::vector<std::unique_ptr<ParameterLib::ParameterBase>>
const& parameters,
61 unsigned const integration_order,
63 std::vector<std::unique_ptr<MeshLib::Mesh>>
const& meshes,
64 std::map<
int, std::shared_ptr<MaterialPropertyLib::Medium>>
const& media)
69 DBUG(
"Create LiquidFlowProcess.");
81 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>
83 process_variables.push_back(std::move(per_process_variables));
90 std::vector<double>
const b =
93 int const mesh_space_dimension =
95 if (
static_cast<int>(b.size()) != mesh_space_dimension)
98 "specific body force (gravity vector) has {:d} components, but the "
99 "space dimension is {:d}.",
100 b.size(), mesh_space_dimension);
103 Eigen::VectorXd specific_body_force(b.size());
107 std::copy_n(b.data(), b.size(), specific_body_force.data());
110 std::unique_ptr<ProcessLib::SurfaceFluxData> surfaceflux;
111 auto calculatesurfaceflux_config =
114 if (calculatesurfaceflux_config)
117 *calculatesurfaceflux_config, meshes);
123 auto const is_linear =
128 INFO(
"LiquidFlow process is set to be linear.");
131 DBUG(
"Check the media properties of LiquidFlow process ...");
133 DBUG(
"Media properties verified.");
135 auto const* aperture_size_parameter = &ParameterLib::findParameter<double>(
138 auto const aperture_config =
143 aperture_size_parameter = &ParameterLib::findParameter<double>(
145 *aperture_config,
"parameter", parameters, 1);
149 std::move(media_map),
152 mesh_space_dimension,
153 std::move(specific_body_force),
155 *aperture_size_parameter,
158 return std::make_unique<LiquidFlowProcess>(
159 std::move(name), mesh, std::move(jacobian_assembler), parameters,
160 integration_order, std::move(process_variables),
161 std::move(process_data), std::move(secondary_variables),
162 std::move(surfaceflux), is_linear);
void INFO(fmt::format_string< Args... > fmt, Args &&... args)
void DBUG(fmt::format_string< Args... > fmt, Args &&... args)
std::optional< ConfigTree > getConfigSubtreeOptional(std::string const &root) const
T getConfigParameter(std::string const ¶m) const
ConfigTree getConfigSubtree(std::string const &root) const
void checkConfigParameter(std::string const ¶m, std::string_view const value) const
std::vector< Node * > const & getNodes() const
Get the nodes-vector for the mesh.
std::vector< Element * > const & getElements() const
Get the element-vector for the mesh.
unsigned getDimension() const
Returns the dimension of the mesh (determined by the maximum dimension over all elements).
static PROCESSLIB_EXPORT const std::string constant_one_parameter_name
Handles configuration of several secondary variables from the project file.
void checkMaterialSpatialDistributionMap(MeshLib::Mesh const &mesh, MaterialPropertyLib::MaterialSpatialDistributionMap const &media_map, ContainerMedium const &required_properties_medium, ContainerSolid const &required_properties_solid_phase, ContainerLiquid const &required_properties_liquid_phase, ContainerGas const &required_properties_gas_phase)
MaterialSpatialDistributionMap createMaterialSpatialDistributionMap(std::map< int, std::shared_ptr< Medium > > const &media, MeshLib::Mesh const &mesh)
Eigen::Map< const Vector > toVector(std::vector< double > const &data, Eigen::VectorXd::Index size)
Creates an Eigen mapped vector from the given data vector.
std::vector< Eigen::MatrixXd > getElementRotationMatrices(int const space_dimension, int const mesh_dimension, std::vector< Element * > const &elements)
Element rotation matrix computation.
int getSpaceDimension(std::vector< Node * > const &nodes)
Computes dimension of the embedding space containing the set of given points.
void checkMPLProperties(MeshLib::Mesh const &mesh, MaterialPropertyLib::MaterialSpatialDistributionMap const &media_map)
std::unique_ptr< Process > createLiquidFlowProcess(std::string const &name, MeshLib::Mesh &mesh, std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&jacobian_assembler, std::vector< ProcessVariable > const &variables, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const ¶meters, unsigned const integration_order, BaseLib::ConfigTree const &config, std::vector< std::unique_ptr< MeshLib::Mesh > > const &meshes, std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)
std::vector< std::reference_wrapper< ProcessVariable > > findProcessVariables(std::vector< ProcessVariable > const &variables, BaseLib::ConfigTree const &pv_config, std::initializer_list< std::string > tags)
void createSecondaryVariables(BaseLib::ConfigTree const &config, SecondaryVariableCollection &secondary_variables)
static std::unique_ptr< ProcessLib::SurfaceFluxData > createSurfaceFluxData(BaseLib::ConfigTree const &calculatesurfaceflux_config, std::vector< std::unique_ptr< MeshLib::Mesh > > const &meshes)