From 6ed5e2353676524599bdb1c86ed8d2dbd646aedd Mon Sep 17 00:00:00 2001
From: Kutalmis Bercin <kutalmis.bercin@esi-group.com>
Date: Mon, 9 Aug 2021 16:34:43 +0100
Subject: [PATCH] BUG: nut{k,U}WallFunction: use viscous predictions in yPlus
(fixes #1773,#1838)
---
.../nutUWallFunctionFvPatchScalarField.C | 30 ++++++++++++--
.../nutkWallFunctionFvPatchScalarField.C | 40 +++++++++++++++----
2 files changed, 60 insertions(+), 10 deletions(-)
diff --git a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUWallFunction/nutUWallFunctionFvPatchScalarField.C b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUWallFunction/nutUWallFunctionFvPatchScalarField.C
index 003bbc48764..b513f54afa6 100644
--- a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUWallFunction/nutUWallFunctionFvPatchScalarField.C
+++ b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUWallFunction/nutUWallFunctionFvPatchScalarField.C
@@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2016 OpenFOAM Foundation
- Copyright (C) 2019 OpenCFD Ltd.
+ Copyright (C) 2019-2021 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@@ -183,7 +183,8 @@ Foam::tmp<Foam::scalarField>
Foam::nutUWallFunctionFvPatchScalarField::yPlus() const
{
const label patchi = patch().index();
- const turbulenceModel& turbModel = db().lookupObject<turbulenceModel>
+
+ const auto& turbModel = db().lookupObject<turbulenceModel>
(
IOobject::groupName
(
@@ -191,10 +192,33 @@ Foam::nutUWallFunctionFvPatchScalarField::yPlus() const
internalField().group()
)
);
+
+ const scalarField& y = turbModel.y()[patchi];
+
+ tmp<scalarField> tnuw = turbModel.nu(patchi);
+ const scalarField& nuw = tnuw();
+
+ tmp<scalarField> tnuEff = turbModel.nuEff(patchi);
+ const scalarField& nuEff = tnuEff();
+
const fvPatchVectorField& Uw = U(turbModel).boundaryField()[patchi];
const scalarField magUp(mag(Uw.patchInternalField() - Uw));
+ const scalarField magGradUw(mag(Uw.snGrad()));
- return calcYPlus(magUp);
+ tmp<scalarField> tyPlus = calcYPlus(magUp);
+ scalarField& yPlus = tyPlus.ref();
+
+ forAll(yPlus, facei)
+ {
+ if (yPlusLam_ > yPlus[facei])
+ {
+ // viscous sublayer
+ yPlus[facei] =
+ y[facei]*sqrt(nuEff[facei]*magGradUw[facei])/nuw[facei];
+ }
+ }
+
+ return tyPlus;
}
diff --git a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkWallFunction/nutkWallFunctionFvPatchScalarField.C b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkWallFunction/nutkWallFunctionFvPatchScalarField.C
index 826a8f354fc..cf8159e4282 100644
--- a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkWallFunction/nutkWallFunctionFvPatchScalarField.C
+++ b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkWallFunction/nutkWallFunctionFvPatchScalarField.C
@@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2016, 2019 OpenFOAM Foundation
- Copyright (C) 2019 OpenCFD Ltd.
+ Copyright (C) 2019-2021 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@@ -33,7 +33,6 @@ License
#include "wallFvPatch.H"
#include "addToRunTimeSelectionTable.H"
-
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
Foam::tmp<Foam::scalarField> Foam::nutkWallFunctionFvPatchScalarField::
@@ -142,7 +141,7 @@ yPlus() const
{
const label patchi = patch().index();
- const turbulenceModel& turbModel = db().lookupObject<turbulenceModel>
+ const auto& turbModel = db().lookupObject<turbulenceModel>
(
IOobject::groupName
(
@@ -153,12 +152,39 @@ yPlus() const
const scalarField& y = turbModel.y()[patchi];
- const tmp<volScalarField> tk = turbModel.k();
+ tmp<volScalarField> tk = turbModel.k();
const volScalarField& k = tk();
- tmp<scalarField> kwc = k.boundaryField()[patchi].patchInternalField();
- tmp<scalarField> nuw = turbModel.nu(patchi);
+ tmp<scalarField> tkwc = k.boundaryField()[patchi].patchInternalField();
+ const scalarField& kwc = tkwc();
+
+ tmp<scalarField> tnuw = turbModel.nu(patchi);
+ const scalarField& nuw = tnuw();
+
+ tmp<scalarField> tnuEff = turbModel.nuEff(patchi);
+ const scalarField& nuEff = tnuEff();
+
+ const fvPatchVectorField& Uw = U(turbModel).boundaryField()[patchi];
+ const scalarField magGradUw(mag(Uw.snGrad()));
+
+ const scalar Cmu25 = pow025(Cmu_);
+
+ auto tyPlus = tmp<scalarField>::New(patch().size(), Zero);
+ auto& yPlus = tyPlus.ref();
+
+ forAll(yPlus, facei)
+ {
+ // inertial sublayer
+ yPlus[facei] = Cmu25*y[facei]*sqrt(kwc[facei])/nuw[facei];
+
+ if (yPlusLam_ > yPlus[facei])
+ {
+ // viscous sublayer
+ yPlus[facei] =
+ y[facei]*sqrt(nuEff[facei]*magGradUw[facei])/nuw[facei];
+ }
+ }
- return pow025(Cmu_)*y*sqrt(kwc)/nuw;
+ return tyPlus;
}
--
GitLab