|
[Sponsors] |
March 19, 2014, 18:07 |
limiter function in openfoam
|
#1 |
Senior Member
Join Date: Jan 2013
Posts: 372
Rep Power: 14 |
Dear All,
In the openfoam limiter function, there is a scalar called k_, does anybody know what does it stand for? For example, we can have a look at the limiter for Gamma differencing: Code:
GammaLimiter(Istream& is) : k_(readScalar(is)) { if (k_ < 0 || k_ > 1) { FatalIOErrorIn("GammaLimiter(Istream& is)", is) << "coefficient = " << k_ << " should be >= 0 and <= 1" << exit(FatalIOError); } // Rescale k_ to be >= 0 and <= 0.5 (TVD conformant) // and avoid the /0 when k_ = 0 k_ = max(k_/2.0, SMALL); } Code:
OpenFOAM/OpenFOAM-<version>/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/Gamma/Gamma.H Code:
OpenFOAM/OpenFOAM-<version>/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/LimitedScheme/NVDTVD.H Code:
scalar r ( const scalar faceFlux, const scalar phiP, const scalar phiN, const vector& gradcP, const vector& gradcN, const vector& d ) const { scalar gradf = phiN - phiP; scalar gradcf; if (faceFlux > 0) { gradcf = d & gradcP; } else { gradcf = d & gradcN; } if (mag(gradcf) >= 1000*mag(gradf)) { return 2*1000*sign(gradcf)*sign(gradf) - 1; } else { return 2*(gradcf/gradf) - 1; } } Code:
scalar r ( const scalar faceFlux, const scalar phiP, const scalar phiN, const vector& gradcP, const vector& gradcN, const vector& d ) const { scalar gradcf; if (faceFlux > 0) { gradcf = d & gradcP; scalar gradf = phiN - phiP; } else { gradcf = d & gradcN; scalar gradf = phiP - phiN; } if (mag(gradcf) >= 1000*mag(gradf)) { return 2*1000*sign(gradcf)*sign(gradf) - 1; } else { return 2*(gradcf/gradf) - 1; } } OFFO |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
foamToTecplot360 | thomasduerr | OpenFOAM Post-Processing | 121 | June 11, 2021 10:05 |
Production limiter only used for k - kOmegaSST - OpenFOAM 2.2.x | Bokse | OpenFOAM Programming & Development | 8 | September 21, 2015 08:55 |
probes function not working in Openfoam 2.1.1 | nandiganavishal | OpenFOAM | 7 | July 26, 2012 11:06 |
LiencubiclowRemodel | nzy102 | OpenFOAM Bugs | 14 | January 10, 2012 08:53 |
Version 15 on Mac OS X | gschaider | OpenFOAM Installation | 113 | December 2, 2009 10:23 |