zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

MULES: Adjust limiter only at boundaries for which the field value is fixed

Browse Source 
Resolves bug-report https://bugs.openfoam.org/view.php?id=2470
This commit is contained in:
Henry Weller
2017-02-23 12:03:50 +00:00
parent 4a22ab94b5
commit cb2bc60fa5
2 changed files with 20 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
src/finiteVolume/fvMatrices/solvers/MULES/CMULESTemplates.C
View File

@ -259,6 +259,16 @@ void Foam::MULES::limiterCorr
psiMinn[pfCelli] = min(psiMinn[pfCelli], psiPf[pFacei]); 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) forAll(phiCorrPf, pFacei)
{ {

10
src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
View File

@ -301,6 +301,16 @@ void Foam::MULES::limiter
psiMinn[pfCelli] = min(psiMinn[pfCelli], psiPf[pFacei]); 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) forAll(phiCorrPf, pFacei)
{ {