From 3ba4f024b3dce73fec1b8d66dfdb012c21a17bb9 Mon Sep 17 00:00:00 2001
From: Gabi Seitz <gabriele.seitz@iws.uni-stuttgart.de>
Date: Thu, 21 Dec 2017 14:13:49 +0100
Subject: [PATCH] [thermochemreaction] fix simple reaction kinetics
---
test/porousmediumflow/1pncmin/implicit/thermochemreaction.hh | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/test/porousmediumflow/1pncmin/implicit/thermochemreaction.hh b/test/porousmediumflow/1pncmin/implicit/thermochemreaction.hh
index 69a7a06950..3293859c72 100644
--- a/test/porousmediumflow/1pncmin/implicit/thermochemreaction.hh
+++ b/test/porousmediumflow/1pncmin/implicit/thermochemreaction.hh
@@ -228,13 +228,13 @@ public:
Scalar krd = 0.05;
- Scalar rDehydration = (volVars.density(cPhaseIdx)- realSolidDensityAverage)
+ Scalar rDehydration = -(volVars.density(cPhaseIdx)- realSolidDensityAverage)
* krd * (Teq-T)/ Teq;
qMass = rDehydration;
}
- if(Teq -T <= 0) qMass = 0;
+ if(Teq -T == 0) qMass = 0;
return qMass;
}
--
GitLab