From 51a7b0ad85b9c79b1b68270d5d60b0b8c6a325c0 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Mon, 7 Mar 2011 12:58:15 +0000
Subject: [PATCH] HerschelBulkley viscosity model: Corrected sign

---
 .../viscosityModels/HerschelBulkley/HerschelBulkley.C         | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/transportModels/incompressible/viscosityModels/HerschelBulkley/HerschelBulkley.C b/src/transportModels/incompressible/viscosityModels/HerschelBulkley/HerschelBulkley.C
index 90d7164476..52085f8b69 100644
--- a/src/transportModels/incompressible/viscosityModels/HerschelBulkley/HerschelBulkley.C
+++ b/src/transportModels/incompressible/viscosityModels/HerschelBulkley/HerschelBulkley.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2004-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2004-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -54,7 +54,7 @@ Foam::viscosityModels::HerschelBulkley::calcNu() const
     dimensionedScalar rtone("rtone", dimless/dimTime, 1.0);
     tmp<volScalarField> sr(strainRate());
     return (min(nu0_,(tau0_ + k_* rtone *( pow(tone * sr(), n_)
-        + pow(tone*tau0_/nu0_,n_))) / (max(sr(), dimensionedScalar
+        - pow(tone*tau0_/nu0_,n_))) / (max(sr(), dimensionedScalar
         ("VSMALL", dimless/dimTime, VSMALL)))));
 }