Skip to content
Snippets Groups Projects
Commit 385f01a4 authored by renchao.lu's avatar renchao.lu Committed by Dmitri Naumov
Browse files

[PL/CT] Rename as pore diffusion.

parent b90fb0f7
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
ProcessLib/ComponentTransport/ComponentTransportFEM.h
+ 10
− 10
View file @ 385f01a4
...@@ -479,11 +479,11 @@ public: ...@@ -479,11 +479,11 @@ public:
.property(MaterialPropertyLib::PropertyType::decay_rate) .property(MaterialPropertyLib::PropertyType::decay_rate)
.template value<double>(vars, pos, t, dt); .template value<double>(vars, pos, t, dt);
auto const& molecular_diffusion_coefficient = auto const& pore_diffusion_coefficient =
MaterialPropertyLib::formEigenTensor<GlobalDim>( MaterialPropertyLib::formEigenTensor<GlobalDim>(
component component
.property(MaterialPropertyLib::PropertyType:: .property(
molecular_diffusion) MaterialPropertyLib::PropertyType::pore_diffusion)
.value(vars, pos, t, dt)); .value(vars, pos, t, dt));
auto const& K = MaterialPropertyLib::formEigenTensor<GlobalDim>( auto const& K = MaterialPropertyLib::formEigenTensor<GlobalDim>(
...@@ -518,7 +518,7 @@ public: ...@@ -518,7 +518,7 @@ public:
GlobalDimMatrixType const hydrodynamic_dispersion = GlobalDimMatrixType const hydrodynamic_dispersion =
velocity_magnitude != 0.0 velocity_magnitude != 0.0
? GlobalDimMatrixType(porosity * ? GlobalDimMatrixType(porosity *
molecular_diffusion_coefficient + pore_diffusion_coefficient +
solute_dispersivity_transverse * solute_dispersivity_transverse *
velocity_magnitude * I + velocity_magnitude * I +
(solute_dispersivity_longitudinal - (solute_dispersivity_longitudinal -
...@@ -526,7 +526,7 @@ public: ...@@ -526,7 +526,7 @@ public:
velocity_magnitude * velocity * velocity_magnitude * velocity *
velocity.transpose()) velocity.transpose())
: GlobalDimMatrixType(porosity * : GlobalDimMatrixType(porosity *
molecular_diffusion_coefficient + pore_diffusion_coefficient +
solute_dispersivity_transverse * solute_dispersivity_transverse *
velocity_magnitude * I); velocity_magnitude * I);
const double R_times_phi(retardation_factor * porosity); const double R_times_phi(retardation_factor * porosity);
...@@ -807,11 +807,11 @@ public: ...@@ -807,11 +807,11 @@ public:
.property(MaterialPropertyLib::PropertyType::decay_rate) .property(MaterialPropertyLib::PropertyType::decay_rate)
.template value<double>(vars, pos, t, dt); .template value<double>(vars, pos, t, dt);
auto const& molecular_diffusion_coefficient = auto const& pore_diffusion_coefficient =
MaterialPropertyLib::formEigenTensor<GlobalDim>( MaterialPropertyLib::formEigenTensor<GlobalDim>(
component component
.property(MaterialPropertyLib::PropertyType:: .property(
molecular_diffusion) MaterialPropertyLib::PropertyType::pore_diffusion)
.value(vars, pos, t, dt)); .value(vars, pos, t, dt));
auto const& K = MaterialPropertyLib::formEigenTensor<GlobalDim>( auto const& K = MaterialPropertyLib::formEigenTensor<GlobalDim>(
...@@ -833,7 +833,7 @@ public: ...@@ -833,7 +833,7 @@ public:
GlobalDimMatrixType const hydrodynamic_dispersion = GlobalDimMatrixType const hydrodynamic_dispersion =
velocity_magnitude != 0.0 velocity_magnitude != 0.0
? GlobalDimMatrixType(porosity * ? GlobalDimMatrixType(porosity *
molecular_diffusion_coefficient + pore_diffusion_coefficient +
solute_dispersivity_transverse * solute_dispersivity_transverse *
velocity_magnitude * I + velocity_magnitude * I +
(solute_dispersivity_longitudinal - (solute_dispersivity_longitudinal -
...@@ -841,7 +841,7 @@ public: ...@@ -841,7 +841,7 @@ public:
velocity_magnitude * velocity * velocity_magnitude * velocity *
velocity.transpose()) velocity.transpose())
: GlobalDimMatrixType(porosity * : GlobalDimMatrixType(porosity *
molecular_diffusion_coefficient + pore_diffusion_coefficient +
solute_dispersivity_transverse * solute_dispersivity_transverse *
velocity_magnitude * I); velocity_magnitude * I);
......
ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp
+ 1
− 1
View file @ 385f01a4
...@@ -41,7 +41,7 @@ void checkMPLProperties( ...@@ -41,7 +41,7 @@ void checkMPLProperties(
std::array const required_properties_components = { std::array const required_properties_components = {
MaterialPropertyLib::PropertyType::retardation_factor, MaterialPropertyLib::PropertyType::retardation_factor,
MaterialPropertyLib::PropertyType::decay_rate, MaterialPropertyLib::PropertyType::decay_rate,
MaterialPropertyLib::PropertyType::molecular_diffusion}; MaterialPropertyLib::PropertyType::pore_diffusion};
for (auto const& element : mesh.getElements()) for (auto const& element : mesh.getElements())
{ {
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment