35 namespace incompressible
49 fixedValueFvPatchScalarField(p, iF),
66 fixedValueFvPatchScalarField(ptf, p, iF, mapper),
82 fixedValueFvPatchScalarField(p, iF, dict),
96 fixedValueFvPatchScalarField(nwfpsf),
97 UName_(nwfpsf.UName_),
98 nuName_(nwfpsf.nuName_),
99 kappa_(nwfpsf.kappa_),
111 fixedValueFvPatchScalarField(nwfpsf, iF),
112 UName_(nwfpsf.UName_),
113 nuName_(nwfpsf.nuName_),
114 kappa_(nwfpsf.kappa_),
142 scalar magUpara = magUp[facei];
144 scalar utau =
sqrt((nuSgsw[facei] + nuw[facei])*magFaceGradU[facei]);
153 scalar kUu =
min(kappa_*magUpara/utau, 50);
154 scalar fkUu =
exp(kUu) - 1 - kUu*(1 + 0.5*kUu);
157 - utau/(ry[facei]*nuw[facei])
159 + 1/E_*(fkUu - 1.0/6.0*kUu*
sqr(kUu));
162 - 1.0/(ry[facei]*nuw[facei])
164 - 1/E_*kUu*fkUu/utau;
166 scalar utauNew = utau - f/df;
167 err =
mag((utau - utauNew)/utau);
170 }
while (utau > VSMALL && err > 0.01 && ++iter < 10);
173 max(
sqr(
max(utau, 0))/magFaceGradU[facei] - nuw[facei], 0.0);
181 fixedValueFvPatchScalarField::evaluate();
188 writeEntryIfDifferent<word>(os,
"U",
"U", UName_);
189 writeEntryIfDifferent<word>(os,
"nu",
"nu", nuName_);
192 writeEntry(
"value", os);