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
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}
scalar y
U
Definition: pEqn.H:72
scalar nutPlus
scalar uPlus
scalar yPlus
scalar epsilonPlus
scalar kPlus
scalar magUp
scalar uTau
scalar tauw
symmTensor R
scalar epsilon
scalar k
scalar Rey
scalar nut
compressible::turbulenceModel & turbulence
vector wallNormal(Zero)
label cellId
label patchId(-1)
label faceId(-1)
volScalarField & nu
vector flowDirection
Definition: createFields.H:41