Commit ffb8ce6d authored by Bernd Flemisch's avatar Bernd Flemisch
Browse files

[stokesnc] fix flux variables

Since the molar density is not associated with a component, it was wrong
to set it always to zero. 

git-svn-id: svn:// 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 71da3343
......@@ -121,13 +121,14 @@ protected:
// loop over all components
for (int compIdx=0; compIdx<numComponents; compIdx++){
molarDensity_ = Scalar(0.0);
massFraction_[compIdx] = Scalar(0.0);
diffusionCoeff_[compIdx] = Scalar(0.0);
moleFractionGrad_[compIdx] = Scalar(0.0);
if (phaseCompIdx!=compIdx) //no transport equation parameters needed for the mass balance
molarDensity_ = Scalar(0.0);
// calculate gradients and secondary variables at IPs
for (int scvIdx = 0;
scvIdx < this->fvGeometry_.numScv;
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