35 namespace tabulatedWallFunctions
40 tabulatedWallFunction,
88 *(
exp(kUPlus) -
pow3(kUPlus)/6 - 0.5*
sqr(kUPlus) - kUPlus - 1);
95 +
exp(kUPlus)*(kUPlus + 1)
96 - 2.0/3.0*
pow3(kUPlus)
102 scalar uPlusNew =
uPlus -
f/(df + ROOTVSMALL);
103 error =
mag((
uPlus - uPlusNew)/uPlusNew);
110 <<
"Newton iterations not converged:" <<
nl 111 <<
" iters = " << iter <<
", error = " << error <<
endl;
124 const dictionary&
dict,
128 tabulatedWallFunction(
dict,
mesh, typeName),
129 kappa_(coeffDict_.
get<scalar>(
"kappa")),
130 E_(coeffDict_.
get<scalar>(
"E"))
152 + 1/
E_*(
exp(kUPlus) -
pow3(kUPlus)/6 - 0.5*
sqr(kUPlus) - kUPlus - 1);
167 if (invertedTable_.log10())
169 os <<
"log10(Re), y+, u+:" <<
endl;
172 scalar
uPlus = invertedTable_[i];
180 os <<
"Re, y+, u+:" <<
endl;
183 scalar
uPlus = invertedTable_[i];
184 scalar
Re = this->
Re(uPlus);
List< ReturnType > get(const UPtrList< T > &list, const AccessOp &aop)
List of values generated by applying the access operation to each list item.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
static const label maxIters_
Maximum number of iterations.
virtual scalar Re(const scalar uPlus) const
Return Reynolds number as a function of u+.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
constexpr char nl
The newline '\n' character (0x0a)
virtual void invertFunction()
Invert the function.
Ostream & endl(Ostream &os)
Add newline and flush stream.
const bool writeData(pdfDictionary.get< bool >("writeData"))
Macros for easy insertion into run-time selection tables.
#define forAll(list, i)
Loop across all elements in list.
const dimensionedScalar e
Elementary charge.
dimensionedScalar exp(const dimensionedScalar &ds)
virtual void writeData(Ostream &os) const
Write to Ostream.
scalar E_
Law-of-the-wall E coefficient.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
int debug
Static debugging option.
OBJstream os(runTime.globalPath()/outputName)
defineTypeNameAndDebug(combustionModel, 0)
virtual scalar yPlus(const scalar uPlus) const
Return y+ as a function of u+.
scalarField Re(const UList< complex > &cmplx)
Extract real component.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
#define WarningInFunction
Report a warning using Foam::Warning.
messageStream Info
Information stream (stdout output on master, null elsewhere)
static const scalar tolerance_
Tolerance.
uniformInterpolationTable< scalar > invertedTable_
Inverted table.
static const Foam::dimensionedScalar A("", Foam::dimPressure, 611.21)
SpaldingsLaw(const dictionary &dict, const polyMesh &mesh)
dimensionedScalar log10(const dimensionedScalar &ds)
scalar kappa_
Von Karman constant.
addToRunTimeSelectionTable(functionObject, pointHistory, dictionary)