From 52ce8bae96364f299297bc583a69df71d6511484 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Mon, 16 May 2011 10:08:36 +0100
Subject: [PATCH] HerschelBulkley: Changed the definition of the coefficients
and stress to more commonly used forms.
This does not change the form of the function but the coefficients must be set appropriately.
---
.../HerschelBulkley/HerschelBulkley.C | 24 ++++++++++++++++---
.../viscosityModel/viscosityModel.C | 4 ++--
2 files changed, 23 insertions(+), 5 deletions(-)
diff --git a/src/transportModels/incompressible/viscosityModels/HerschelBulkley/HerschelBulkley.C b/src/transportModels/incompressible/viscosityModels/HerschelBulkley/HerschelBulkley.C
index 52085f8b694..c1b4132323e 100644
--- a/src/transportModels/incompressible/viscosityModels/HerschelBulkley/HerschelBulkley.C
+++ b/src/transportModels/incompressible/viscosityModels/HerschelBulkley/HerschelBulkley.C
@@ -52,10 +52,28 @@ Foam::viscosityModels::HerschelBulkley::calcNu() const
{
dimensionedScalar tone("tone", dimTime, 1.0);
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
- ("VSMALL", dimless/dimTime, VSMALL)))));
+
+ // return
+ // (
+ // min
+ // (
+ // nu0_,
+ // (tau0_ + k_*rtone*(pow(tone*sr(), n_) - pow(tone*tau0_/nu0_, n_)))
+ // /max(sr(), dimensionedScalar("VSMALL", dimless/dimTime, VSMALL))
+ // )
+ // );
+
+ return
+ (
+ min
+ (
+ nu0_,
+ (tau0_ + k_*rtone*pow(tone*sr(), n_))
+ /(max(sr(), dimensionedScalar ("VSMALL", dimless/dimTime, VSMALL)))
+ )
+ );
}
diff --git a/src/transportModels/incompressible/viscosityModels/viscosityModel/viscosityModel.C b/src/transportModels/incompressible/viscosityModels/viscosityModel/viscosityModel.C
index 493cb89df69..57be63e30e7 100644
--- a/src/transportModels/incompressible/viscosityModels/viscosityModel/viscosityModel.C
+++ b/src/transportModels/incompressible/viscosityModels/viscosityModel/viscosityModel.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
@@ -57,7 +57,7 @@ Foam::viscosityModel::viscosityModel
Foam::tmp<Foam::volScalarField> Foam::viscosityModel::strainRate() const
{
- return mag(symm(fvc::grad(U_)));
+ return sqrt(2.0)*mag(symm(fvc::grad(U_)));
}
--
GitLab