Commit d21d8d30 authored by Alexander Kissinger's avatar Alexander Kissinger
Browse files

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
parent e324a4b2
......@@ -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());
......
Supports Markdown
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