From 9fa37cd3eea18fb820faf7ab05143be29723d110 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Fri, 15 Apr 2011 16:07:14 +0100
Subject: [PATCH] kOmegaSST: Rationalisation of S2
---
.../compressible/RAS/kOmegaSST/kOmegaSST.C | 4 ++--
.../incompressible/RAS/kOmegaSST/kOmegaSST.C | 8 ++++----
2 files changed, 6 insertions(+), 6 deletions(-)
diff --git a/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C b/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C
index 54ebb5aa9a1..1c9d53d3918 100644
--- a/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C
+++ b/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C
@@ -390,7 +390,7 @@ void kOmegaSST::correct()
}
tmp<volTensorField> tgradU = fvc::grad(U_);
- volScalarField S2(magSqr(symm(tgradU())));
+ volScalarField S2(2*magSqr(symm(tgradU())));
volScalarField GbyMu((tgradU() && dev(twoSymm(tgradU()))));
volScalarField G("RASModel::G", mut_*GbyMu);
tgradU.clear();
@@ -448,7 +448,7 @@ void kOmegaSST::correct()
// Re-calculate viscosity
- mut_ = a1_*rho_*k_/max(a1_*omega_, F2()*sqrt(2.0*S2));
+ mut_ = a1_*rho_*k_/max(a1_*omega_, F2()*sqrt(S2));
mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity
diff --git a/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C b/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C
index 9e34f054b7c..6baeaf5a0ed 100644
--- a/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C
+++ b/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C
@@ -347,8 +347,8 @@ void kOmegaSST::correct()
y_.correct();
}
- const volScalarField S2(magSqr(symm(fvc::grad(U_))));
- volScalarField G("RASModel::G", nut_*2*S2);
+ const volScalarField S2(2*magSqr(symm(fvc::grad(U_))));
+ volScalarField G("RASModel::G", nut_*S2);
// Update omega and G at the wall
omega_.boundaryField().updateCoeffs();
@@ -368,7 +368,7 @@ void kOmegaSST::correct()
- fvm::Sp(fvc::div(phi_), omega_)
- fvm::laplacian(DomegaEff(F1), omega_)
==
- gamma(F1)*2*S2
+ gamma(F1)*S2
- fvm::Sp(beta(F1)*omega_, omega_)
- fvm::SuSp
(
@@ -402,7 +402,7 @@ void kOmegaSST::correct()
// Re-calculate viscosity
- nut_ = a1_*k_/max(a1_*omega_, F2()*sqrt(2*S2));
+ nut_ = a1_*k_/max(a1_*omega_, F2()*sqrt(S2));
nut_.correctBoundaryConditions();
}
--
GitLab