diff --git a/src/finiteVolume/fvMatrices/solvers/MULES/CMULESTemplates.C b/src/finiteVolume/fvMatrices/solvers/MULES/CMULESTemplates.C index b9ffcc9bc4..e55889eaa3 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 675e4a3bd6..4ad5686625 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) {