From 827e3659b97b0bbe118879e595dda222eba9c0e0 Mon Sep 17 00:00:00 2001
From: Mark Olesen <Mark.Olesen@Germany>
Date: Fri, 31 Jul 2009 18:15:54 +0200
Subject: [PATCH] consistency update: kappa=0.41, E=9.8
- this would be an argument for providing default values at the top-level
compressible/incompressible turbulenceModel
---
.../multiphase/bubbleFoam/createRASTurbulence.H | 4 ++--
.../solvers/multiphase/settlingFoam/createFields.H | 4 ++--
.../applyBoundaryLayer/applyBoundaryLayer.C | 14 +++++++-------
3 files changed, 11 insertions(+), 11 deletions(-)
diff --git a/applications/solvers/multiphase/bubbleFoam/createRASTurbulence.H b/applications/solvers/multiphase/bubbleFoam/createRASTurbulence.H
index f9794a914f6..1ae6219fd6f 100644
--- a/applications/solvers/multiphase/bubbleFoam/createRASTurbulence.H
+++ b/applications/solvers/multiphase/bubbleFoam/createRASTurbulence.H
@@ -82,7 +82,7 @@
(
"kappa",
wallFunctionDict,
- 0.4187
+ 0.41
)
);
@@ -92,7 +92,7 @@
(
"E",
wallFunctionDict,
- 9.0
+ 9.8
)
);
diff --git a/applications/solvers/multiphase/settlingFoam/createFields.H b/applications/solvers/multiphase/settlingFoam/createFields.H
index 9f0ca00596b..6acb19994ae 100644
--- a/applications/solvers/multiphase/settlingFoam/createFields.H
+++ b/applications/solvers/multiphase/settlingFoam/createFields.H
@@ -304,7 +304,7 @@
(
"kappa",
wallFunctionDict,
- 0.4187
+ 0.41
)
);
@@ -314,7 +314,7 @@
(
"E",
wallFunctionDict,
- 9.0
+ 9.8
)
);
diff --git a/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C b/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
index 8aae7e124c9..5714cccbffc 100644
--- a/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
+++ b/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
@@ -41,6 +41,11 @@ Description
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+// turbulence constants - file-scope
+static const scalar Cmu(0.09);
+static const scalar kappa(0.41);
+
+
int main(int argc, char *argv[])
{
argList::validOptions.insert("ybl", "scalar");
@@ -48,7 +53,6 @@ int main(int argc, char *argv[])
argList::validOptions.insert("writenut", "");
# include "setRootCase.H"
-
# include "createTime.H"
# include "createMesh.H"
@@ -117,10 +121,6 @@ int main(int argc, char *argv[])
phi = fvc::interpolate(U) & mesh.Sf();
phi.write();
- // Set turbulence constants
- dimensionedScalar kappa("kappa", dimless, 0.4187);
- dimensionedScalar Cmu("Cmu", dimless, 0.09);
-
// Read and modify turbulence fields if present
IOobject epsilonHeader
@@ -182,11 +182,11 @@ int main(int argc, char *argv[])
Info<< "Reading field epsilon\n" << endl;
volScalarField epsilon(epsilonHeader, mesh);
- scalar ck0 = ::pow(Cmu.value(), 0.25)*kappa.value();
+ scalar ck0 = ::pow(Cmu, 0.25)*kappa;
k = sqr(nut/(ck0*min(y, ybl)));
k.correctBoundaryConditions();
- scalar ce0 = ::pow(Cmu.value(), 0.75)/kappa.value();
+ scalar ce0 = ::pow(Cmu, 0.75)/kappa;
epsilon = ce0*k*sqrt(k)/min(y, ybl);
epsilon.correctBoundaryConditions();
--
GitLab