research/OpenFOAM-4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

MULES: Improve handling of very small geometries

Browse source
This commit is contained in:
Henry Weller 2017-07-24 10:35:00 +01:00
parent 4cfa55f2ea
commit 03c63bf7a5
2 changed files with 10 additions and 10 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

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

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2013-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2013-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -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)
{
@ -361,7 +361,7 @@ void Foam::MULES::limiterCorr
max(min
(
(sumlPhip[celli] + psiMaxn[celli])
/(mSumPhim[celli] - SMALL),
/(mSumPhim[celli] + ROOTVSMALL),
1.0), 0.0
);
@ -369,7 +369,7 @@ void Foam::MULES::limiterCorr
max(min
(
(mSumlPhim[celli] + psiMinn[celli])
/(sumPhip[celli] + SMALL),
/(sumPhip[celli] + ROOTVSMALL),
1.0), 0.0
);
}

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

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -236,8 +236,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)
{
@ -428,7 +428,7 @@ void Foam::MULES::limiter
max(min
(
(sumlPhip[celli] + psiMaxn[celli])
/(mSumPhim[celli] - SMALL),
/(mSumPhim[celli] + ROOTVSMALL),
1.0), 0.0
);
@ -436,7 +436,7 @@ void Foam::MULES::limiter
max(min
(
(mSumlPhim[celli] + psiMinn[celli])
/(sumPhip[celli] + SMALL),
/(sumPhip[celli] + ROOTVSMALL),
1.0), 0.0
);
}