diff --git a/src/finiteVolume/fvMatrices/solvers/MULES/CMULESTemplates.C b/src/finiteVolume/fvMatrices/solvers/MULES/CMULESTemplates.C index b9ffcc9bc4dc844cda584039514867af76673735..e55889eaa365b10d489af1396bcbb7e856d8b189 100644 --- a/src/finiteVolume/fvMatrices/solvers/MULES/CMULESTemplates.C +++ b/src/finiteVolume/fvMatrices/solvers/MULES/CMULESTemplates.C @@ -259,6 +259,16 @@ void Foam::MULES::limiterCorr psiMinn[pfCelli] = min(psiMinn[pfCelli], psiPf[pFacei]); } } + else + { + forAll(phiCorrPf, pFacei) + { + const label pfCelli = pFaceCells[pFacei]; + + psiMaxn[pfCelli] = max(psiMaxn[pfCelli], psiMax); + psiMinn[pfCelli] = min(psiMinn[pfCelli], psiMin); + } + } forAll(phiCorrPf, pFacei) { diff --git a/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C b/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C index 675e4a3bd664d2b909420b6b060de64946c4f88a..4ad56866257f02236dac5904b1c1f3f5c7fd8b3f 100644 --- a/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C +++ b/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C @@ -301,6 +301,16 @@ void Foam::MULES::limiter psiMinn[pfCelli] = min(psiMinn[pfCelli], psiPf[pFacei]); } } + else + { + forAll(phiCorrPf, pFacei) + { + const label pfCelli = pFaceCells[pFacei]; + + psiMaxn[pfCelli] = max(psiMaxn[pfCelli], psiMax); + psiMinn[pfCelli] = min(psiMinn[pfCelli], psiMin); + } + } forAll(phiCorrPf, pFacei) {