Bugfix: make TPFA version of local linearizer work for two-phase runs.

This commit is contained in:
Atgeirr Flø Rasmussen 2022-09-01 13:52:58 +02:00
parent 63ea7b29c1
commit 92e9bd5c37

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))