From 0848c825215ef2213ee5b1338f074fe53f0b5451 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 19 Mar 2012 10:06:39 +0000
Subject: [PATCH] ENH: Spalart-Allmaras - update nut_ from nuTilda_*fv1 when
 turbulence inactive

---
 .../incompressible/RAS/SpalartAllmaras/SpalartAllmaras.C   | 7 ++++++-
 1 file changed, 6 insertions(+), 1 deletion(-)

diff --git a/src/turbulenceModels/incompressible/RAS/SpalartAllmaras/SpalartAllmaras.C b/src/turbulenceModels/incompressible/RAS/SpalartAllmaras/SpalartAllmaras.C
index 213896ff7d..1804ba5b01 100644
--- a/src/turbulenceModels/incompressible/RAS/SpalartAllmaras/SpalartAllmaras.C
+++ b/src/turbulenceModels/incompressible/RAS/SpalartAllmaras/SpalartAllmaras.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -367,6 +367,10 @@ void SpalartAllmaras::correct()
 
     if (!turbulence_)
     {
+        // Re-calculate viscosity
+        nut_ = nuTilda_*fv1(this->chi());
+        nut_.correctBoundaryConditions();
+
         return;
     }
 
@@ -401,6 +405,7 @@ void SpalartAllmaras::correct()
     bound(nuTilda_, dimensionedScalar("0", nuTilda_.dimensions(), 0.0));
     nuTilda_.correctBoundaryConditions();
 
+    // Re-calculate viscosity
     nut_.internalField() = fv1*nuTilda_.internalField();
     nut_.correctBoundaryConditions();
 }