40 #ifndef cellLimitedGrad_H 41 #define cellLimitedGrad_H 58 template<
class Type,
class Limiter>
109 if (k_ < 0 || k_ > 1)
114 ) <<
"coefficient = " << k_
115 <<
" should be >= 0 and <= 1" 126 const scalar maxDelta,
127 const scalar minDelta,
128 const scalar extrapolate
134 const Type& maxDelta,
135 const Type& minDelta,
136 const Type& extrapolate
155 template<
class Type,
class Limiter>
159 const scalar maxDelta,
160 const scalar minDelta,
161 const scalar extrapolate
166 if (extrapolate > small)
168 r = maxDelta/extrapolate;
170 else if (extrapolate < -small)
172 r = minDelta/extrapolate;
183 template<
class Type,
class Limiter>
187 const Type& maxDelta,
188 const Type& minDelta,
189 const Type& extrapolate
192 for (
direction cmpt=0; cmpt<pTraits<Type>::nComponents; cmpt++)
TypeName("cellLimited")
RunTime type information.
virtual tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh > > calcGrad(const GeometricField< Type, fvPatchField, volMesh > &vsf, const word &name) const
Return the gradient of the given field to the gradScheme::grad.
errorManipArg< error, int > exit(error &err, const int errNo=1)
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
Abstract base class with a fat-interface to all derived classes covering all possible ways in which t...
typeOfRank< typename pTraits< arg1 >::cmptType, direction(pTraits< arg1 >::rank)+direction(pTraits< arg2 >::rank) >::type type
Generic GeometricField class.
Mesh data needed to do the Finite Volume discretisation.
void limiter(scalarField &allLambda, const RdeltaTType &rDeltaT, const RhoType &rho, const volScalarField &psi, const surfaceScalarField &phiBD, const surfaceScalarField &phiCorr, const SpType &Sp, const SuType &Su, const PsiMaxType &psiMax, const PsiMinType &psiMin)
A class for handling words, derived from string.
cellLimitedGrad gradient scheme applied to a runTime selected base gradient scheme.
bool readScalar(const char *buf, doubleScalar &s)
Read whole of buf as a scalar. Return true if successful.
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
word name(const complex &)
Return a string representation of a complex.
Abstract base class for gradient schemes.
#define FatalIOErrorInFunction(ios)
Report an error message using Foam::FatalIOError.
void limitFaceCmpt(scalar &limiter, const scalar maxDelta, const scalar minDelta, const scalar extrapolate) const
Mesh data needed to do the Finite Volume discretisation.
static tmp< gradScheme< Type > > New(const fvMesh &mesh, Istream &schemeData)
Return a pointer to a new gradScheme created on freestore.
A class for managing temporary objects.
const fvMesh & mesh() const
Return mesh reference.
label & setComponent(label &l, const direction)
void component(FieldField< Field, typename FieldField< Field, Type >::cmptType > &sf, const FieldField< Field, Type > &f, const direction d)
void limitFace(Type &limiter, const Type &maxDelta, const Type &minDelta, const Type &extrapolate) const