From b236147be29c454231ec329292fd65a23fcdc4a0 Mon Sep 17 00:00:00 2001
From: andy <a.heather@opencfd.co.uk>
Date: Fri, 5 Nov 2010 12:59:51 +0000
Subject: [PATCH] BUG: Corrected missing sqrt(2) factors in compressible
kOmegaSST model
---
src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C b/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C
index 09e07ef57b0..96d14423eca 100644
--- a/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C
+++ b/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C
@@ -268,7 +268,7 @@ kOmegaSST::kOmegaSST
/ max
(
a1_*omega_,
- F2()*sqrt(magSqr(symm(fvc::grad(U_))))
+ F2()*sqrt(2.0*magSqr(symm(fvc::grad(U_))))
)
);
mut_.correctBoundaryConditions();
@@ -365,7 +365,7 @@ void kOmegaSST::correct()
// Re-calculate viscosity
mut_ =
a1_*rho_*k_
- /max(a1_*omega_, F2()*sqrt(magSqr(symm(fvc::grad(U_)))));
+ /max(a1_*omega_, F2()*sqrt(2.0*magSqr(symm(fvc::grad(U_)))));
mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity
@@ -446,7 +446,7 @@ void kOmegaSST::correct()
// Re-calculate viscosity
- mut_ = a1_*rho_*k_/max(a1_*omega_, F2()*sqrt(S2));
+ mut_ = a1_*rho_*k_/max(a1_*omega_, F2()*sqrt(2.0*S2));
mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity
--
GitLab