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