26 std::vector<std::unique_ptr<ParameterLib::ParameterBase>>
const& parameters)
28 DBUG(
"Create PorousMediaProperties.");
30 auto const& porous_medium_configs =
34 std::vector<std::unique_ptr<MaterialLib::PorousMedium::Permeability>>
35 intrinsic_permeability_models;
36 std::vector<std::unique_ptr<MaterialLib::PorousMedium::Porosity>>
38 std::vector<std::unique_ptr<MaterialLib::PorousMedium::Storage>>
41 std::vector<int> mat_ids;
42 for (
auto const& porous_medium_config :
44 porous_medium_configs.getConfigSubtreeList(
"porous_medium"))
48 mat_ids.push_back(
id);
50 auto const& porosity_config =
52 porous_medium_config.getConfigSubtree(
"porosity");
53 porosity_models.emplace_back(
59 auto const& permeability_config =
61 porous_medium_config.getConfigSubtree(
"permeability");
62 intrinsic_permeability_models.emplace_back(
64 permeability_config, parameters));
67 auto const& storage_config =
69 porous_medium_config.getConfigSubtree(
"storage");
70 storage_models.emplace_back(
78 auto const material_ids = materialIDs(mesh);
79 int const max_material_id =
82 : *std::max_element(begin(*material_ids), end(*material_ids));
84 if (max_material_id >
static_cast<int>(mat_ids.size() - 1))
87 "The maximum value of MaterialIDs in mesh is {:d}. As the "
88 "given number of porous media definitions in the project file is "
89 "{:d}, the maximum value of MaterialIDs in mesh must be {:d} "
90 "(index starts with zero).",
91 max_material_id, mat_ids.size(), max_material_id - 1);
94 if (max_material_id <
static_cast<int>(mat_ids.size() - 1))
96 "There are {:d} porous medium definitions in the project file but "
97 "only {:d} different values in the MaterialIDs vector/data_array "
99 mat_ids.size(), max_material_id - 1);
101 if (mat_ids.back() !=
static_cast<int>(mat_ids.size()) - 1)
104 "The ids in the porous media definitions in the project file have "
105 "to be sequential, starting with the id zero.");
109 std::move(intrinsic_permeability_models),
110 std::move(storage_models), material_ids};
void DBUG(fmt::format_string< Args... > fmt, Args &&... args)
void WARN(fmt::format_string< Args... > fmt, Args &&... args)
Definition of the Mesh class.
T getConfigAttribute(std::string const &attr) const
ConfigTree getConfigSubtree(std::string const &root) const
void reorderVector(std::vector< ValueType > &v, std::vector< IndexType > const &order)
std::unique_ptr< Permeability > createPermeabilityModel(BaseLib::ConfigTree const &config, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const ¶meters)
std::unique_ptr< Storage > createStorageModel(BaseLib::ConfigTree const &config)
PorousMediaProperties createPorousMediaProperties(MeshLib::Mesh &mesh, BaseLib::ConfigTree const &configs, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const ¶meters)
std::unique_ptr< Porosity > createPorosityModel(BaseLib::ConfigTree const &config, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const ¶meters)