From b21c617a81f9aa47a6fc2fc6fd8aacc0652f470b Mon Sep 17 00:00:00 2001
From: Kutalmis Bercin <kutalmis.bercin@esi-group.com>
Date: Mon, 5 Jul 2021 12:00:32 +0100
Subject: [PATCH] ENH: turbulenceFields: update turbulence-intensity expression
---
.../field/turbulenceFields/turbulenceFieldsTemplates.C | 9 +++++----
1 file changed, 5 insertions(+), 4 deletions(-)
diff --git a/src/functionObjects/field/turbulenceFields/turbulenceFieldsTemplates.C b/src/functionObjects/field/turbulenceFields/turbulenceFieldsTemplates.C
index 5f17e494343..c94b118b843 100644
--- a/src/functionObjects/field/turbulenceFields/turbulenceFieldsTemplates.C
+++ b/src/functionObjects/field/turbulenceFields/turbulenceFieldsTemplates.C
@@ -118,14 +118,15 @@ Foam::functionObjects::turbulenceFields::I
const Model& model
) const
{
- // Assume k is available
- const volScalarField uPrime(sqrt((2.0/3.0)*model.k()));
- const dimensionedScalar U0("U0", dimVelocity, SMALL);
+ // (P:p. 183)
+ // root-mean-square of velocity fluctuations - isotropic turbulence
+ tmp<volScalarField> uPrime = sqrt((2.0/3.0)*model.k());
+ const dimensionedScalar U0(dimVelocity, SMALL);
return tmp<volScalarField>::New
(
"I.tmp",
- uPrime/max(max(uPrime, mag(model.U())), U0)
+ uPrime/max(mag(model.U()), U0)
);
}
--
GitLab