Skip to content
Snippets Groups Projects
Commit 5a7e382b authored by Timo Koch's avatar Timo Koch
Browse files

[doc][mineralization] Fix doc and solid terminology

solid phases->solidComponents
parent fe47be42
No related branches found
No related tags found
2 merge requests!2170[doc][mineralization] Fix doc and solid terminology,!2134WIP Feature/timestepper test
...@@ -76,15 +76,15 @@ public: ...@@ -76,15 +76,15 @@ public:
{ {
auto storage = ParentType::computeStorage(problem, scv, volVars); auto storage = ParentType::computeStorage(problem, scv, volVars);
const auto massOrMoleDensity = [](const auto& volVars, const int phaseIdx) const auto massOrMoleSolidDensity = [](const auto& volVars, const int sCompIdx)
{ return useMoles ? volVars.solidComponentMolarDensity(phaseIdx) : volVars.solidComponentDensity(phaseIdx); }; { return useMoles ? volVars.solidComponentMolarDensity(sCompIdx) : volVars.solidComponentDensity(sCompIdx); };
// compute storage term of all components within all fluid phases // compute storage term of all solid components
for (int phaseIdx = 0; phaseIdx < numSolidComps-numInertSolidComps; ++phaseIdx) for (int sCompIdx = 0; sCompIdx < numSolidComps-numInertSolidComps; ++sCompIdx)
{ {
auto eqIdx = Indices::conti0EqIdx + numComponents + phaseIdx; auto eqIdx = Indices::conti0EqIdx + numComponents + sCompIdx;
storage[eqIdx] += volVars.solidVolumeFraction(phaseIdx) storage[eqIdx] += volVars.solidVolumeFraction(sCompIdx)
* massOrMoleDensity(volVars, phaseIdx); * massOrMoleSolidDensity(volVars, sCompIdx);
} }
return storage; return storage;
......
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