exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Fixed FPE for NVDTVD schemes; forced vanAlbada flux limiter to be possitive

Browse source
This commit is contained in:
Oliver Borm (boroli) 2011-01-04 10:36:48 +01:00
parent e5a96484f3
commit ba01e50738
3 changed files with 3 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/LimitedScheme/NVDTVD.H
View file

@ -113,7 +113,7 @@ public:
} }
// Stabilise for division // Stabilise for division
gradf = stabilise(gradf, VSMALL); gradf = stabilise(gradf, SMALL);
return 2*(gradcf/gradf) - 1; return 2*(gradcf/gradf) - 1;
} }

2
src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/LimitedScheme/NVDVTVDV.H
View file

@ -114,7 +114,7 @@ public:
} }
// Stabilise for division // Stabilise for division
gradf = stabilise(gradf, VSMALL); gradf = stabilise(gradf, SMALL);
return 2*(gradcf/gradf) - 1; return 2*(gradcf/gradf) - 1;
} }

2
src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/vanAlbada/vanAlbada.H
View file

@ -78,7 +78,7 @@ public:
faceFlux, phiP, phiN, gradcP, gradcN, d faceFlux, phiP, phiN, gradcP, gradcN, d
); );
return r*(r + 1)/(sqr(r) + 1); return max(r*(r + 1)/(sqr(r) + 1),0);
} }
}; };