zhyang-dev/OpenFOAM-6
1
0
Fork 0
mirror of https://github.com/OpenFOAM/OpenFOAM-6.git synced 2025-12-08 06:57:46 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

cubicGradientLimiter: Documented private data

Browse Source
This commit is contained in:
Henry Weller
2018-02-27 17:37:28 +00:00
parent c9c9c62a24
commit 9ea6e22b0a
3 changed files with 62 additions and 38 deletions
Download Patch File Download Diff File Expand all files Collapse all files

17
src/finiteVolume/finiteVolume/gradSchemes/limitedGradSchemes/cellLimitedGrad/gradientLimiters/VenkatakrishnanGradientLimiter.H
View File

@ -83,13 +83,18 @@ class Venkatakrishnan
public:
Venkatakrishnan(Istream&)
{}
// Constructors
inline scalar limiter(const scalar r) const
{
return (sqr(r) + 2*r)/(sqr(r) + r + 2);
}
Venkatakrishnan(Istream&)
{}
// Member Functions
inline scalar limiter(const scalar r) const
{
return (sqr(r) + 2*r)/(sqr(r) + r + 2);
}
};

66
src/finiteVolume/finiteVolume/gradSchemes/limitedGradSchemes/cellLimitedGrad/gradientLimiters/cubicGradientLimiter.H
View File

@ -80,40 +80,54 @@ namespace gradientLimiters
class cubic
{
const scalar rt_;
const scalar a_;
const scalar b_;
// Private Data
//- Limiter transition point at which the limiter function -> 1
// Must be > 1
const scalar rt_;
//- Coefficient of the r^3 term (evaluated from rt)
const scalar a_;
// - Coefficient of the r^2 term (evaluated from rt)
const scalar b_;
public:
cubic(Istream& schemeData)
:
rt_(readScalar(schemeData)),
a_(2.0/sqr(rt_) - 2.0/pow3(rt_)),
b_(-(3.0/2.0)*a_*rt_)
{
if (rt_ < 1)
{
FatalIOErrorInFunction
(
schemeData
) << "coefficient = " << rt_
<< " should be > 1"
<< exit(FatalIOError);
}
}
// Constructors
inline scalar limiter(const scalar r) const
{
if (r < rt_)
cubic(Istream& schemeData)
:
rt_(readScalar(schemeData)),
a_(2.0/sqr(rt_) - 2.0/pow3(rt_)),
b_(-(3.0/2.0)*a_*rt_)
{
return ((a_*r + b_)*r + 1)*r;
if (rt_ < 1)
{
FatalIOErrorInFunction
(
schemeData
) << "coefficient = " << rt_
<< " should be > 1"
<< exit(FatalIOError);
}
}
else
// Member Functions
inline scalar limiter(const scalar r) const
{
return 1;
if (r < rt_)
{
return ((a_*r + b_)*r + 1)*r;
}
else
{
return 1;
}
}
}
};

17
src/finiteVolume/finiteVolume/gradSchemes/limitedGradSchemes/cellLimitedGrad/gradientLimiters/minmodGradientLimiter.H
View File

@ -69,13 +69,18 @@ class minmod
public:
minmod(Istream&)
{}
// Constructors
inline scalar limiter(const scalar r) const
{
return min(r, 1);
}
minmod(Istream&)
{}
// Member Functions
inline scalar limiter(const scalar r) const
{
return min(r, 1);
}
};