evaluateNearWall.H
Go to the documentation of this file.
1 {
2  // Evaluate near-wall behaviour
3 
4  const scalar nu = turbulence->nu()().boundaryField()[patchId][faceId];
5  const scalar nut = turbulence->nut()().boundaryField()[patchId][faceId];
6  const scalar nuEff = turbulence->nuEff()().boundaryField()[patchId][faceId];
7  const symmTensor R =
8  turbulence->devSigma()().boundaryField()[patchId][faceId];
9  const scalar epsilon = turbulence->epsilon()()[cellId];
10  const scalar omega = turbulence->omega()()[cellId];
11  const scalar k = turbulence->k()()[cellId];
12  const scalar magUp = mag(U[cellId] - U.boundaryField()[patchId][faceId]);
13 
14  const scalar tauw = flowDirection & R & wallNormal;
15  const scalar uTau = ::sqrt(mag(tauw));
16 
17  const scalar yPlus = uTau*y[cellId]/(nu + rootVSmall);
18  const scalar uPlus = magUp/(uTau + rootVSmall);
19  const scalar nutPlus = nut/nu;
20  const scalar nuEffPlus = nuEff/nu;
21  const scalar kPlus = k/(sqr(uTau) + rootVSmall);
22  const scalar epsilonPlus = epsilon*nu/(pow4(uTau) + rootVSmall);
23  const scalar omegaPlus = omega*nu/(sqr(uTau) + rootVSmall);
24 
25  const scalar Rey = magUp*y[cellId]/nu;
26 
27  Info<< "Rey = " << Rey
28  << ", uTau = " << uTau
29  << ", nut+ = " << nutPlus
30  << ", nuEff+ = " << nuEffPlus
31  << ", y+ = " << yPlus << ", u+ = " << uPlus
32  << ", k+ = " << kPlus << ", epsilon+ = " << epsilonPlus
33  << ", omega+ = " << omegaPlus
34  << endl;
35 }
scalar y
const scalar omega
const symmTensor R
const scalar yPlus
const scalar uTau
const scalar nut
const scalar kPlus
const scalar nutPlus
const scalar omegaPlus
const scalar uPlus
const scalar nuEff
const scalar epsilon
const scalar Rey
const scalar k
const scalar nuEffPlus
const scalar magUp
const scalar epsilonPlus
const scalar tauw
label patchId(-1)
vector wallNormal(Zero)
const label cellId
label faceId(-1)
U
Definition: pEqn.H:72
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:257
SymmTensor< scalar > symmTensor
SymmTensor of scalars.
Definition: symmTensor.H:48
messageStream Info
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensioned< scalar > mag(const dimensioned< Type > &)
dimensionedScalar pow4(const dimensionedScalar &ds)
vector flowDirection
Definition: createFields.H:41
autoPtr< incompressible::momentumTransportModel > turbulence(incompressible::momentumTransportModel::New(U, phi, viscosity))