zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: nutUBlended wall function - added protection for uTau. See #696

Browse Source
This commit is contained in:
Andrew Heather
2018-01-08 10:33:56 +00:00
parent f02783b986
commit 4781193e31
1 changed files with 12 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUBlendedWallFunction/nutUBlendedWallFunctionFvPatchScalarField.C
View File

@ -93,17 +93,20 @@ Foam::nutUBlendedWallFunctionFvPatchScalarField::calcUTau
forAll(uTaup, facei) forAll(uTaup, facei)
{ {
scalar ut = sqrt((nutw[facei] + nuw[facei])*magGradU[facei]); scalar ut = sqrt((nutw[facei] + nuw[facei])*magGradU[facei]);
scalar error = GREAT; if (mag(ut) > ROOTVSMALL)
label iter = 0;
while (iter++ < 10 && error > 0.001)
{ {
scalar yPlus = y[facei]*ut/nuw[facei]; scalar error = GREAT;
scalar uTauVis = magUp[facei]/yPlus; label iter = 0;
scalar uTauLog = kappa_*magUp[facei]/log(E_*yPlus); while (iter++ < 10 && error > 0.001)
{
scalar yPlus = y[facei]*ut/nuw[facei];
scalar uTauVis = magUp[facei]/yPlus;
scalar uTauLog = kappa_*magUp[facei]/log(E_*yPlus);
scalar utNew = pow(pow(uTauVis, n_) + pow(uTauLog, n_), 1.0/n_); scalar utNew = pow(pow(uTauVis, n_) + pow(uTauLog, n_), 1.0/n_);
error = mag(ut - utNew)/(ut + ROOTVSMALL); error = mag(ut - utNew)/(ut + ROOTVSMALL);
ut = 0.5*(ut + utNew); ut = 0.5*(ut + utNew);
}
} }
uTaup[facei] = ut; uTaup[facei] = ut;
} }