Skip to content
Snippets Groups Projects
Commit 4224c19d authored by Martin Schneider's avatar Martin Schneider
Browse files

[fix][dispersion][boc] Use extrusion for area calculation

parent 9a636b28
No related branches found
No related tags found
1 merge request!3803Fix/tpfa dispersion
......@@ -119,7 +119,7 @@ public:
}
// compute the dispersion flux
componentFlux[compIdx] = -1.0 * rhoMassOrMole * vtmv(scvf.unitOuterNormal(), dispersionTensor, gradX) * scvf.area();
componentFlux[compIdx] = -1.0 * rhoMassOrMole * vtmv(scvf.unitOuterNormal(), dispersionTensor, gradX)*Extrusion::area(fvGeometry, scvf);
if (BalanceEqOpts::mainComponentIsBalanced(phaseIdx) && !FluidSystem::isTracerFluidSystem())
componentFlux[phaseIdx] -= componentFlux[compIdx];
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment