Skip to content
Snippets Groups Projects
Commit 1e38155c authored by renchao.lu's avatar renchao.lu
Browse files

[PL/CT] store temperature field in process data.

parent f8991d44
No related branches found
No related tags found
No related merge requests found
......@@ -15,6 +15,7 @@
#include "ChemistryLib/ChemicalSolverInterface.h"
#include "MaterialLib/MPL/CreateMaterialSpatialDistributionMap.h"
#include "MathLib/LinAlg/Eigen/EigenMapTools.h"
#include "ParameterLib/Parameter.h"
namespace MaterialPropertyLib
{
......@@ -38,6 +39,9 @@ struct ComponentTransportProcessData
Eigen::VectorXd const specific_body_force;
bool const has_gravity;
bool const non_advective_form;
/// This optional tag provides a simple means of considering the temperature
/// effect on the solute transport process.
ParameterLib::Parameter<double> const* const temperature;
/**
* When this optional tag is on, the feedback of chemical reactions on the
* porosity will be counted. The change of porosity equals to the summation
......
......@@ -15,6 +15,7 @@
#include "ComponentTransportProcessData.h"
#include "MaterialLib/MPL/CreateMaterialSpatialDistributionMap.h"
#include "MeshLib/IO/readMeshFromFile.h"
#include "ParameterLib/Utils.h"
#include "ProcessLib/Output/CreateSecondaryVariables.h"
#include "ProcessLib/SurfaceFlux/SurfaceFluxData.h"
#include "ProcessLib/Utils/ProcessUtils.h"
......@@ -221,6 +222,10 @@ std::unique_ptr<Process> createComponentTransportProcess(
config.getConfigParameter<bool>("chemically_induced_porosity_change",
false);
auto const temperature = ParameterLib::findOptionalTagParameter<double>(
//! \ogs_file_param_special{prj__processes__process__ComponentTransport__temperature_field}
config, "temperature_field", parameters, 1, &mesh);
auto media_map =
MaterialPropertyLib::createMaterialSpatialDistributionMap(media, mesh);
......@@ -233,6 +238,7 @@ std::unique_ptr<Process> createComponentTransportProcess(
specific_body_force,
has_gravity,
non_advective_form,
temperature,
chemically_induced_porosity_change,
chemical_solver_interface.get(),
hydraulic_process_id,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment