From 0915b391bdbd7ca9641bb6606f50087bd8e10bc6 Mon Sep 17 00:00:00 2001
From: Timo Koch <timo.koch@iws.uni-stuttgart.de>
Date: Fri, 17 Feb 2017 16:54:24 +0100
Subject: [PATCH] [test][mixeddimension][doc] Improve doc on 1p2c_richards2c
test
---
.../embedded/1p2c_richards2c/richardstestproblem.hh | 4 +++-
.../embedded/1p2c_richards2c/rootsystemtestproblem.hh | 4 +++-
2 files changed, 6 insertions(+), 2 deletions(-)
diff --git a/test/mixeddimension/embedded/1p2c_richards2c/richardstestproblem.hh b/test/mixeddimension/embedded/1p2c_richards2c/richardstestproblem.hh
index 4350a2bd43..97737e10a8 100644
--- a/test/mixeddimension/embedded/1p2c_richards2c/richardstestproblem.hh
+++ b/test/mixeddimension/embedded/1p2c_richards2c/richardstestproblem.hh
@@ -81,7 +81,7 @@ SET_BOOL_PROP(RichardsTestProblem, VtkAddVelocity, true);
// Set the grid parameter group
SET_STRING_PROP(RichardsTestProblem, GridParameterGroup, "SoilGrid");
-// Use mass fractions
+// Use mole fractions
SET_BOOL_PROP(RichardsTestProblem, UseMoles, true);
}
@@ -235,6 +235,7 @@ public:
sourceValues[conti0EqIdx] = 2* M_PI *rootRadius * Kr *(lowDimVolVars.pressure(wPhaseIdx) - bulkVolVars.pressure(wPhaseIdx))
*bulkVolVars.molarDensity(wPhaseIdx);
+ //! advective transport over root wall
// compute correct upwind concentration
if (sourceValues[conti0EqIdx] > 0)
sourceValues[transportEqIdx] = sourceValues[conti0EqIdx]
@@ -243,6 +244,7 @@ public:
sourceValues[transportEqIdx] = sourceValues[conti0EqIdx]
* bulkVolVars.moleFraction(wPhaseIdx, transportCompIdx);
+ //! diffusive transport over root wall
sourceValues[transportEqIdx] += 2* M_PI *rootRadius * 1.0e-8
*(lowDimVolVars.moleFraction(wPhaseIdx, transportCompIdx) - bulkVolVars.moleFraction(wPhaseIdx, transportCompIdx))
*0.5*(bulkVolVars.molarDensity(wPhaseIdx) + lowDimVolVars.molarDensity(wPhaseIdx));
diff --git a/test/mixeddimension/embedded/1p2c_richards2c/rootsystemtestproblem.hh b/test/mixeddimension/embedded/1p2c_richards2c/rootsystemtestproblem.hh
index 26a95ff07d..b9343bbca6 100644
--- a/test/mixeddimension/embedded/1p2c_richards2c/rootsystemtestproblem.hh
+++ b/test/mixeddimension/embedded/1p2c_richards2c/rootsystemtestproblem.hh
@@ -77,7 +77,7 @@ SET_BOOL_PROP(RootsystemTestProblem, ProblemEnableGravity, true);
// Enable velocity output
SET_BOOL_PROP(RootsystemTestProblem, VtkAddVelocity, true);
-// Use mass fractions
+// Use mole fractions
SET_BOOL_PROP(RootsystemTestProblem, UseMoles, true);
}
@@ -315,6 +315,7 @@ public:
sourceValues[conti0EqIdx] = 2* M_PI *rootRadius * Kr *(bulkVolVars.pressure(wPhaseIdx) - lowDimVolVars.pressure(wPhaseIdx))
*bulkVolVars.molarDensity(wPhaseIdx);
+ //! advective transport over root wall
// compute correct upwind concentration
if (sourceValues[conti0EqIdx] < 0)
sourceValues[transportEqIdx] = sourceValues[conti0EqIdx]
@@ -323,6 +324,7 @@ public:
sourceValues[transportEqIdx] = sourceValues[conti0EqIdx]
* bulkVolVars.moleFraction(wPhaseIdx, transportCompIdx);
+ //! diffusive transport over root wall
sourceValues[transportEqIdx] += 2* M_PI *rootRadius * 1.0e-8
*(bulkVolVars.moleFraction(wPhaseIdx, transportCompIdx) - lowDimVolVars.moleFraction(wPhaseIdx, transportCompIdx))
*0.5*(bulkVolVars.molarDensity(wPhaseIdx) + lowDimVolVars.molarDensity(wPhaseIdx));
--
GitLab