From 9b10b18ced60b0bec34f5cc0f46d88b82c501824 Mon Sep 17 00:00:00 2001
From: henry <henry>
Date: Tue, 27 Jan 2009 21:19:02 +0000
Subject: [PATCH] Corrected fv2 for the "Ashford" correction of the SA model.

---
 .../LES/compressible/SpalartAllmaras/SpalartAllmaras.C                 | 1 +
 .../LES/incompressible/SpalartAllmaras/SpalartAllmaras.C               | 1 +
 .../RAS/compressible/SpalartAllmaras/SpalartAllmaras.C                 | 3 ++-
 .../RAS/incompressible/SpalartAllmaras/SpalartAllmaras.C               | 3 ++-
 4 files changed, 6 insertions(+), 2 deletions(-)

diff --git a/src/turbulenceModels/LES/compressible/SpalartAllmaras/SpalartAllmaras.C b/src/turbulenceModels/LES/compressible/SpalartAllmaras/SpalartAllmaras.C
index 386168d..568598a 100644
--- a/src/turbulenceModels/LES/compressible/SpalartAllmaras/SpalartAllmaras.C
+++ b/src/turbulenceModels/LES/compressible/SpalartAllmaras/SpalartAllmaras.C
@@ -55,6 +55,7 @@ tmp<volScalarField> SpalartAllmaras::fv1() const
 tmp<volScalarField> SpalartAllmaras::fv2() const
 {
     volScalarField chi = nuTilda_/(mu()/rho());
+    //return scalar(1) - chi/(scalar(1) + chi*fv1());
     return 1.0/pow3(scalar(1) + chi/Cv2_);
 }
 
diff --git a/src/turbulenceModels/LES/incompressible/SpalartAllmaras/SpalartAllmaras.C b/src/turbulenceModels/LES/incompressible/SpalartAllmaras/SpalartAllmaras.C
index f8de6b6..6dd450c 100644
--- a/src/turbulenceModels/LES/incompressible/SpalartAllmaras/SpalartAllmaras.C
+++ b/src/turbulenceModels/LES/incompressible/SpalartAllmaras/SpalartAllmaras.C
@@ -55,6 +55,7 @@ tmp<volScalarField> SpalartAllmaras::fv1() const
 tmp<volScalarField> SpalartAllmaras::fv2() const
 {
     volScalarField chi = nuTilda_/nu();
+    //return scalar(1) - chi/(scalar(1) + chi*fv1());
     return 1.0/pow3(scalar(1) + chi/Cv2_);
 }
 
diff --git a/src/turbulenceModels/RAS/compressible/SpalartAllmaras/SpalartAllmaras.C b/src/turbulenceModels/RAS/compressible/SpalartAllmaras/SpalartAllmaras.C
index 0ee4318..d0757ab 100644
--- a/src/turbulenceModels/RAS/compressible/SpalartAllmaras/SpalartAllmaras.C
+++ b/src/turbulenceModels/RAS/compressible/SpalartAllmaras/SpalartAllmaras.C
@@ -63,7 +63,8 @@ tmp<volScalarField> SpalartAllmaras::fv2
     const volScalarField& fv1
 ) const
 {
-    return 1.0/pow3(scalar(1) + chi/Cv2_);
+    return 1.0 - chi/(1.0 + chi*fv1);
+    //return 1.0/pow3(scalar(1) + chi/Cv2);
 }
 
 
diff --git a/src/turbulenceModels/RAS/incompressible/SpalartAllmaras/SpalartAllmaras.C b/src/turbulenceModels/RAS/incompressible/SpalartAllmaras/SpalartAllmaras.C
index 7a8fc53..7c0253d 100644
--- a/src/turbulenceModels/RAS/incompressible/SpalartAllmaras/SpalartAllmaras.C
+++ b/src/turbulenceModels/RAS/incompressible/SpalartAllmaras/SpalartAllmaras.C
@@ -63,7 +63,8 @@ tmp<volScalarField> SpalartAllmaras::fv2
     const volScalarField& fv1
 ) const
 {
-    return 1.0/pow3(scalar(1) + chi/Cv2_);
+    return 1.0 - chi/(1.0 + chi*fv1);
+    //return 1.0/pow3(scalar(1) + chi/Cv2);
 }
 
 
-- 
2.11.4.GIT