Commit 41e4d68f authored by Henry Weller's avatar Henry Weller Committed by Andrew Heather
Browse files

MULES: Improve handling of very small geometries

parent d27457d4
......@@ -202,8 +202,8 @@ void Foam::MULES::limiterCorr
scalarField psiMaxn(psiIf.size(), psiMin);
scalarField psiMinn(psiIf.size(), psiMax);
scalarField sumPhip(psiIf.size(), VSMALL);
scalarField mSumPhim(psiIf.size(), VSMALL);
scalarField sumPhip(psiIf.size(), 0.0);
scalarField mSumPhim(psiIf.size(), 0.0);
forAll(phiCorrIf, facei)
{
......@@ -371,7 +371,7 @@ void Foam::MULES::limiterCorr
max(min
(
(sumlPhip[celli] + psiMaxn[celli])
/(mSumPhim[celli] - SMALL),
/(mSumPhim[celli] + ROOTVSMALL),
1.0), 0.0
);
......@@ -379,7 +379,7 @@ void Foam::MULES::limiterCorr
max(min
(
(mSumlPhim[celli] + psiMinn[celli])
/(sumPhip[celli] + SMALL),
/(sumPhip[celli] + ROOTVSMALL),
1.0), 0.0
);
}
......
......@@ -240,8 +240,8 @@ void Foam::MULES::limiter
scalarField sumPhiBD(psiIf.size(), 0.0);
scalarField sumPhip(psiIf.size(), VSMALL);
scalarField mSumPhim(psiIf.size(), VSMALL);
scalarField sumPhip(psiIf.size(), 0.0);
scalarField mSumPhim(psiIf.size(), 0.0);
forAll(phiCorrIf, facei)
{
......@@ -442,7 +442,7 @@ void Foam::MULES::limiter
max(min
(
(sumlPhip[celli] + psiMaxn[celli])
/(mSumPhim[celli] - SMALL),
/(mSumPhim[celli] + ROOTVSMALL),
1.0), 0.0
);
......@@ -450,7 +450,7 @@ void Foam::MULES::limiter
max(min
(
(mSumlPhim[celli] + psiMinn[celli])
/(sumPhip[celli] + SMALL),
/(sumPhip[celli] + ROOTVSMALL),
1.0), 0.0
);
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment