47 a_(dict.
lookup<scalar>(
"a")),
48 b_(dict.
lookup<scalar>(
"b")),
49 c_(dict.
lookup<scalar>(
"c")),
50 d_(dict.
lookup<scalar>(
"d")),
51 e_(dict.
lookup<scalar>(
"e"))
57 inline Foam::scalar Foam::SRIFallOffFunction::operator()
63 const scalar logPr =
log10(
max(Pr, small));
65 const scalar X = 1/(1 +
sqr(logPr));
67 const scalar
psi = a_*
exp(-b_/T) +
exp(-T/c_);
69 return d_*
pow(psi, X)*
pow(T, e_);
80 const scalar logPr =
log10(
max(Pr, small));
82 const scalar X = 1/(1 +
sqr(logPr));
84 const scalar
psi = a_*
exp(-b_/T) +
exp(-T/c_);
85 const scalar dpsidT = a_*b_/
sqr(T)*
exp(-b_/T) - 1/c_*
exp(-T/c_);
87 return F*(X/psi*dpsidT + e_/
T);
98 static const scalar logTen =
log(scalar(10));
100 const scalar logPr =
log10(
max(Pr, small));
101 const scalar dlogPrdPr = Pr >= small ? 1/(logTen*Pr) : 0;
103 const scalar X = 1/(1 +
sqr(logPr));
104 const scalar dXdPr = -
sqr(X)*2*logPr*dlogPrdPr;
106 const scalar
psi = a_*
exp(-b_/T) +
exp(-T/c_);
108 return F*
log(psi)*dXdPr;
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
virtual Ostream & write(const char)=0
Write character.
The SRI fall-off function.
dimensionedScalar log(const dimensionedScalar &ds)
A list of keyword definitions, which are a keyword followed by any number of values (e...
dimensionedSymmTensor sqr(const dimensionedVector &dv)
const dimensionedScalar b
Wien displacement law constant: default SI units: [m K].
void write(Ostream &os) const
Write to stream.
const dimensionedScalar c
Speed of light in a vacuum.
scalar ddPr(const scalar T, const scalar Pr, const scalar F) const
stressControl lookup("compactNormalStress") >> compactNormalStress
SRIFallOffFunction(const scalar a, const scalar b, const scalar c, const scalar d, const scalar e)
Construct from components.
dimensionedScalar exp(const dimensionedScalar &ds)
scalar ddT(const scalar T, const scalar Pr, const scalar F) const
const volScalarField & psi
An Ostream is an abstract base class for all output systems (streams, files, token lists...
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
void writeEntry(Ostream &os, const HashTable< T, Key, Hash > &ht)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar log10(const dimensionedScalar &ds)
const dimensionedScalar e
Elementary charge.