From d21d8d30dc2dac914c0f33db252c79887c80b72a Mon Sep 17 00:00:00 2001
From: Alexander Kissinger <alexander.kissinger@iws.uni-stuttgart.de>
Date: Tue, 28 Oct 2014 10:14:46 +0000
Subject: [PATCH] Wrong assignment of the Neumann boundary flux for the total
flux equation. Pointed out by Bernd.
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@13567 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
test/implicit/1p2c/1p2cconvectionproblem.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/test/implicit/1p2c/1p2cconvectionproblem.hh b/test/implicit/1p2c/1p2cconvectionproblem.hh
index cd9556f69f..886bccc587 100644
--- a/test/implicit/1p2c/1p2cconvectionproblem.hh
+++ b/test/implicit/1p2c/1p2cconvectionproblem.hh
@@ -393,7 +393,7 @@ public:
if(globalPos[0] < eps_)
{
- values[pressureIdx] = -darcyVelocity_*elemVolVars[scvIdx].molarDensity()*elemVolVars[scvIdx].moleFraction(pressureIdx);
+ values[pressureIdx] = -darcyVelocity_*elemVolVars[scvIdx].molarDensity();
values[massOrMoleFracIdx] = -darcyVelocity_*elemVolVars[scvIdx].molarDensity()*elemVolVars[scvIdx].moleFraction(massOrMoleFracIdx);
values[temperatureIdx] = -darcyVelocity_*elemVolVars[scvIdx].density()
*IapwsH2O::liquidEnthalpy(temperatureHigh_, elemVolVars[scvIdx].pressure());
--
GitLab