Merge pull request #713 from atgeirr/fix-fast-assembly-twophase

Bugfix: make TPFA version of local linearizer work for two-phase runs.
This commit is contained in:
Atgeirr Flø Rasmussen 2022-09-01 16:50:55 +02:00 committed by GitHub
commit d034fc3b58

View File

@ -126,6 +126,9 @@ public:
storage = 0.0;
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
if (!FluidSystem::phaseIsActive(phaseIdx)) {
continue;
}
unsigned activeCompIdx = Indices::canonicalToActiveComponentIndex(FluidSystem::solventComponentIndex(phaseIdx));
LhsEval surfaceVolume =
Toolbox::template decay<LhsEval>(fs.saturation(phaseIdx))