36 namespace incompressible
46 const label patchI = patch().index();
51 const scalarField& nuw = rasModel.nu().boundaryField()[patchI];
56 sqr(
calcUTau(magGradU))/(magGradU + ROOTVSMALL) - nuw
70 rasModel.U().boundaryField()[patch().index()];
73 const scalarField& nuw = rasModel.nu().boundaryField()[patch().index()];
81 scalar magUpara = magUp[facei];
83 scalar ut =
sqrt((nutw[facei] + nuw[facei])*magGradU[facei]);
92 scalar kUu =
min(kappa_*magUpara/ut, 50);
93 scalar fkUu =
exp(kUu) - 1 - kUu*(1 + 0.5*kUu);
96 - ut*y[facei]/nuw[facei]
98 + 1/E_*(fkUu - 1.0/6.0*kUu*
sqr(kUu));
105 scalar uTauNew = ut + f/df;
106 err =
mag((ut - uTauNew)/ut);
109 }
while (ut > VSMALL && err > 0.01 && ++iter < 10);
110 uTau[facei] =
max(0.0, ut);
182 const label patchI = patch().index();
187 const scalarField& nuw = rasModel.nu().boundaryField()[patchI];
199 writeLocalEntries(os);
200 writeEntry(
"value", os);