32 std::map<
int, std::shared_ptr<MaterialPropertyLib::Medium>>
const& media)
34 std::array
const required_medium_properties = {
41 std::array
const required_liquid_properties = {
45 for (
auto const& m : media)
47 checkRequiredProperties(*m.second, required_medium_properties);
48 checkRequiredProperties(m.second->phase(
"AqueousLiquid"),
49 required_liquid_properties);
50 checkRequiredProperties(m.second->phase(
"Solid"),
51 required_solid_properties);
57 std::string
const& name,
59 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&& jacobian_assembler,
60 std::vector<ProcessVariable>
const& variables,
61 std::vector<std::unique_ptr<ParameterLib::ParameterBase>>
const& parameters,
62 std::optional<ParameterLib::CoordinateSystem>
const&
63 local_coordinate_system,
64 unsigned const integration_order,
66 std::map<
int, std::shared_ptr<MaterialPropertyLib::Medium>>
const& media)
70 DBUG(
"Create RichardsMechanicsProcess.");
72 auto const coupling_scheme =
75 const bool use_monolithic_scheme =
76 !(coupling_scheme && (*coupling_scheme ==
"staggered"));
85 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>
87 if (use_monolithic_scheme)
97 variable_p = &per_process_variables[0].get();
98 variable_u = &per_process_variables[1].get();
99 process_variables.push_back(std::move(per_process_variables));
103 using namespace std::string_literals;
104 for (
auto const& variable_name : {
"pressure"s,
"displacement"s})
106 auto per_process_variables =
108 process_variables.push_back(std::move(per_process_variables));
110 variable_p = &process_variables[0][0].get();
111 variable_u = &process_variables[1][0].get();
114 DBUG(
"Associate displacement with process variable '{:s}'.",
120 "Number of components of the process variable '{:s}' is different "
121 "from the displacement dimension: got {:d}, expected {:d}",
127 DBUG(
"Associate pressure with process variable '{:s}'.",
132 "Pressure process variable '{:s}' is not a scalar variable but has "
138 auto solid_constitutive_relations =
139 MaterialLib::Solids::createConstitutiveRelations<DisplacementDim>(
140 parameters, local_coordinate_system, config);
144 Eigen::Matrix<double, DisplacementDim, 1> specific_body_force;
146 std::vector<double>
const b =
149 "specific_body_force");
150 if (b.size() != DisplacementDim)
153 "The size of the specific body force vector does not match the "
154 "displacement dimension. Vector size is {:d}, displacement "
156 b.size(), DisplacementDim);
159 std::copy_n(b.data(), b.size(), specific_body_force.data());
164 DBUG(
"Check the media properties of RichardsMechanics process ...");
166 DBUG(
"Media properties verified.");
169 auto const initial_stress = ParameterLib::findOptionalTagParameter<double>(
171 config,
"initial_stress", parameters,
176 std::optional<MicroPorosityParameters> micro_porosity_parameters;
177 if (
auto const micro_porosity_config =
184 micro_porosity_config->getConfigSubtree(
"nonlinear_solver")),
186 micro_porosity_config->getConfigParameter<
double>(
187 "mass_exchange_coefficient")};
190 auto const mass_lumping =
194 auto const explicit_hm_coupling_in_unsaturated_zone =
197 "explicit_hm_coupling_in_unsaturated_zone",
false);
201 std::move(media_map),
202 std::move(solid_constitutive_relations),
205 micro_porosity_parameters,
207 explicit_hm_coupling_in_unsaturated_zone};
213 return std::make_unique<RichardsMechanicsProcess<DisplacementDim>>(
214 std::move(name), mesh, std::move(jacobian_assembler), parameters,
215 integration_order, std::move(process_variables),
216 std::move(process_data), std::move(secondary_variables),
217 use_monolithic_scheme);
template std::unique_ptr< Process > createRichardsMechanicsProcess< 3 >(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, std::optional< ParameterLib::CoordinateSystem > const &local_coordinate_system, unsigned const integration_order, BaseLib::ConfigTree const &config, std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)
std::unique_ptr< Process > createRichardsMechanicsProcess(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, std::optional< ParameterLib::CoordinateSystem > const &local_coordinate_system, unsigned const integration_order, BaseLib::ConfigTree const &config, std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)
template std::unique_ptr< Process > createRichardsMechanicsProcess< 2 >(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, std::optional< ParameterLib::CoordinateSystem > const &local_coordinate_system, unsigned const integration_order, BaseLib::ConfigTree const &config, std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)