diff --git a/ProcessLib/ComponentTransport/ComponentTransportProcessData.h b/ProcessLib/ComponentTransport/ComponentTransportProcessData.h index 6cd51d77049805fd6de8de04b033d557767ec116..e0527fbf63f7d66544bda61994d3f2ab1291ea39 100644 --- a/ProcessLib/ComponentTransport/ComponentTransportProcessData.h +++ b/ProcessLib/ComponentTransport/ComponentTransportProcessData.h @@ -38,6 +38,7 @@ struct ComponentTransportProcessData Eigen::VectorXd const specific_body_force; bool const has_gravity; bool const non_advective_form; + bool const chemically_induced_porosity_change; ChemistryLib::ChemicalSolverInterface* const chemical_solver_interface; const int hydraulic_process_id; diff --git a/ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp b/ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp index 4fc9bf55b24d6ec39f9b1745888d6089fe0beff1..427d5f892d8cac48d4a13f8c0fdff403dc91e025 100644 --- a/ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp +++ b/ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp @@ -216,6 +216,11 @@ std::unique_ptr<Process> createComponentTransportProcess( //! \ogs_file_param{prj__processes__process__ComponentTransport__non_advective_form} config.getConfigParameter<bool>("non_advective_form", false); + bool chemically_induced_porosity_change = + //! \ogs_file_param{prj__processes__process__ComponentTransport__chemically_induced_porosity_change} + config.getConfigParameter<bool>("chemically_induced_porosity_change", + false); + auto media_map = MaterialPropertyLib::createMaterialSpatialDistributionMap(media, mesh); @@ -223,13 +228,15 @@ std::unique_ptr<Process> createComponentTransportProcess( checkMPLProperties(mesh, *media_map); DBUG("Media properties verified."); - ComponentTransportProcessData process_data{std::move(media_map), - specific_body_force, - has_gravity, - non_advective_form, - chemical_solver_interface.get(), - hydraulic_process_id, - first_transport_process_id}; + ComponentTransportProcessData process_data{ + std::move(media_map), + specific_body_force, + has_gravity, + non_advective_form, + chemically_induced_porosity_change, + chemical_solver_interface.get(), + hydraulic_process_id, + first_transport_process_id}; SecondaryVariableCollection secondary_variables;