diff --git a/src/functionObjects/field/turbulenceFields/turbulenceFields.C b/src/functionObjects/field/turbulenceFields/turbulenceFields.C
index 214aa9896a82014fe8b79fac09f1efbe09d5030e..6c697271de62a7d4849267387f3e127929a13ea2 100644
--- a/src/functionObjects/field/turbulenceFields/turbulenceFields.C
+++ b/src/functionObjects/field/turbulenceFields/turbulenceFields.C
@@ -3,7 +3,7 @@
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2013-2016 OpenFOAM Foundation
- \\/ M anipulation | Copyright (C) 2015-2017 OpenCFD Ltd.
+ \\/ M anipulation | Copyright (C) 2015-2018 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@@ -54,12 +54,14 @@ Foam::functionObjects::turbulenceFields::compressibleFieldNames_
{ compressibleField::cfK, "k" },
{ compressibleField::cfEpsilon, "epsilon" },
{ compressibleField::cfOmega, "omega" },
+ { compressibleField::cfNuTilda, "nuTilda" },
{ compressibleField::cfMut, "mut" },
{ compressibleField::cfMuEff, "muEff" },
{ compressibleField::cfAlphat, "alphat" },
{ compressibleField::cfAlphaEff, "alphaEff" },
{ compressibleField::cfR, "R" },
{ compressibleField::cfDevRhoReff, "devRhoReff" },
+ { compressibleField::cfL, "L" }
};
@@ -72,10 +74,12 @@ Foam::functionObjects::turbulenceFields::incompressibleFieldNames_
{ incompressibleField::ifK, "k" },
{ incompressibleField::ifEpsilon, "epsilon" },
{ incompressibleField::ifOmega, "omega" },
+ { incompressibleField::ifNuTilda, "nuTilda" },
{ incompressibleField::ifNut, "nut" },
{ incompressibleField::ifNuEff, "nuEff" },
{ incompressibleField::ifR, "R" },
{ incompressibleField::ifDevReff, "devReff" },
+ { incompressibleField::ifL, "L" },
};
@@ -193,6 +197,11 @@ bool Foam::functionObjects::turbulenceFields::execute()
processField<scalar>(f, omega(model));
break;
}
+ case cfNuTilda:
+ {
+ processField<scalar>(f, nuTilda(model));
+ break;
+ }
case cfMut:
{
processField<scalar>(f, model.mut());
@@ -223,6 +232,11 @@ bool Foam::functionObjects::turbulenceFields::execute()
processField<symmTensor>(f, model.devRhoReff());
break;
}
+ case cfL:
+ {
+ processField<scalar>(f, L(model));
+ break;
+ }
default:
{
FatalErrorInFunction
@@ -256,6 +270,11 @@ bool Foam::functionObjects::turbulenceFields::execute()
processField<scalar>(f, omega(model));
break;
}
+ case ifNuTilda:
+ {
+ processField<scalar>(f, nuTilda(model));
+ break;
+ }
case ifNut:
{
processField<scalar>(f, model.nut());
@@ -276,6 +295,11 @@ bool Foam::functionObjects::turbulenceFields::execute()
processField<symmTensor>(f, model.devReff());
break;
}
+ case ifL:
+ {
+ processField<scalar>(f, L(model));
+ break;
+ }
default:
{
FatalErrorInFunction
diff --git a/src/functionObjects/field/turbulenceFields/turbulenceFields.H b/src/functionObjects/field/turbulenceFields/turbulenceFields.H
index dc6e168abd4ac738e06ced35791ce94052c022e2..a4ae31d19970328e82d8f0788effc3671d2ce1b2 100644
--- a/src/functionObjects/field/turbulenceFields/turbulenceFields.H
+++ b/src/functionObjects/field/turbulenceFields/turbulenceFields.H
@@ -3,7 +3,7 @@
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2013-2016 OpenFOAM Foundation
- \\/ M anipulation | Copyright (C) 2015-2017 OpenCFD Ltd.
+ \\/ M anipulation | Copyright (C) 2015-2018 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@@ -65,6 +65,7 @@ Usage
k | turbulence kinetic energy
epsilon | turbulence kinetic energy dissipation rate
omega | turbulence specific dissipation rate
+ nuTilda | turbulence modified viscosity
nut | turbulence viscosity (incompressible)
nuEff | effective turbulence viscosity (incompressible)
mut | turbulence viscosity (compressible)
@@ -74,6 +75,7 @@ Usage
R | Reynolds stress tensor
devReff | Deviatoric part of the effective Reynolds stress
devRhoReff | Divergence of the Reynolds stress
+ L | turbulence length scale
\endplaintable
See also
@@ -116,12 +118,14 @@ public:
cfK,
cfEpsilon,
cfOmega,
+ cfNuTilda,
cfMut,
cfMuEff,
cfAlphat,
cfAlphaEff,
cfR,
- cfDevRhoReff
+ cfDevRhoReff,
+ cfL
};
static const Enum<compressibleField> compressibleFieldNames_;
@@ -130,10 +134,12 @@ public:
ifK,
ifEpsilon,
ifOmega,
+ ifNuTilda,
ifNut,
ifNuEff,
ifR,
- ifDevReff
+ ifDevReff,
+ ifL
};
static const Enum<incompressibleField> incompressibleFieldNames_;
@@ -165,6 +171,14 @@ protected:
template<class Model>
tmp<volScalarField> omega(const Model& model) const;
+ //- Return nuTilda calculated from k and omega
+ template<class Model>
+ tmp<volScalarField> nuTilda(const Model& model) const;
+
+ //- Return L calculated from k and epsilon
+ template<class Model>
+ tmp<volScalarField> L(const Model& model) const;
+
private:
diff --git a/src/functionObjects/field/turbulenceFields/turbulenceFieldsTemplates.C b/src/functionObjects/field/turbulenceFields/turbulenceFieldsTemplates.C
index 248fcbf2a019b1b276cf89dfe03ef8aef6c13c3e..ec6b70d92dfe6320d6205d28adb774d0f0e80589 100644
--- a/src/functionObjects/field/turbulenceFields/turbulenceFieldsTemplates.C
+++ b/src/functionObjects/field/turbulenceFields/turbulenceFieldsTemplates.C
@@ -3,7 +3,7 @@
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2012-2016 OpenFOAM Foundation
- \\/ M anipulation |
+ \\/ M anipulation | Copyright (C) 2018 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@@ -91,7 +91,7 @@ Foam::functionObjects::turbulenceFields::omega
(
IOobject
(
- "omega",
+ "omega.tmp",
k.mesh().time().timeName(),
k.mesh()
),
@@ -102,4 +102,43 @@ Foam::functionObjects::turbulenceFields::omega
}
+template<class Model>
+Foam::tmp<Foam::volScalarField>
+Foam::functionObjects::turbulenceFields::nuTilda
+(
+ const Model& model
+) const
+{
+ return tmp<volScalarField>
+ (
+ new volScalarField("nuTilda.tmp", model.k()/omega(model))
+ );
+}
+
+
+template<class Model>
+Foam::tmp<Foam::volScalarField>
+Foam::functionObjects::turbulenceFields::L
+(
+ const Model& model
+) const
+{
+ const scalar Cmu = 0.09;
+
+ // Assume k and epsilon are available
+ const volScalarField k(model.k());
+ const volScalarField epsilon(model.epsilon());
+ const dimensionedScalar eps0("eps0", epsilon.dimensions(), SMALL);
+
+ return tmp<volScalarField>
+ (
+ new volScalarField
+ (
+ "L.tmp",
+ pow(Cmu, 0.75)*pow(k, 1.5)/(epsilon + eps0)
+ )
+ );
+}
+
+
// ************************************************************************* //