From 00ebf8e153c037faac5469549bb5de4ebfae6985 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Mon, 2 Mar 2015 22:27:39 +0000
Subject: [PATCH] BinghamPlastic: Change definition of shear-rate by 1/sqrt(2)
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1552
---
.../mixtureViscosityModels/BinghamPlastic/BinghamPlastic.C | 4 ++--
.../mixtureViscosityModels/BinghamPlastic/BinghamPlastic.H | 4 +++-
.../mixtureViscosityModel/mixtureViscosityModel.H | 5 -----
.../driftFluxFoam/ras/dahl/constant/transportProperties | 2 +-
.../ras/mixerVessel2D/constant/transportProperties | 2 +-
.../driftFluxFoam/ras/tank3D/constant/transportProperties | 2 +-
6 files changed, 8 insertions(+), 11 deletions(-)
diff --git a/applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/BinghamPlastic/BinghamPlastic.C b/applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/BinghamPlastic/BinghamPlastic.C
index bf9aa605ca8..01826f0af47 100644
--- a/applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/BinghamPlastic/BinghamPlastic.C
+++ b/applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/BinghamPlastic/BinghamPlastic.C
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2014 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2014-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -112,7 +112,7 @@ Foam::mixtureViscosityModels::BinghamPlastic::mu
(
tauy
/(
- mag(fvc::grad(U_))
+ sqrt(2.0)*mag(symm(fvc::grad(U_)))
+ 1.0e-4*(tauy + tauySmall)/mup
)
+ mup,
diff --git a/applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/BinghamPlastic/BinghamPlastic.H b/applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/BinghamPlastic/BinghamPlastic.H
index c63bb42d82a..b39f3f53e02 100644
--- a/applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/BinghamPlastic/BinghamPlastic.H
+++ b/applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/BinghamPlastic/BinghamPlastic.H
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2014 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2014-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -27,6 +27,8 @@ Class
Description
Viscosity correction model for Bingham plastics.
+ The strain-rate used is defined as sqrt(2.0)*mag(symm(grad(U)))
+
SourceFiles
BinghamPlastic.C
diff --git a/applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/mixtureViscosityModel/mixtureViscosityModel.H b/applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/mixtureViscosityModel/mixtureViscosityModel.H
index 9b7db43f93d..866459ebce2 100644
--- a/applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/mixtureViscosityModel/mixtureViscosityModel.H
+++ b/applications/solvers/multiphase/driftFluxFoam/mixtureViscosityModels/mixtureViscosityModel/mixtureViscosityModel.H
@@ -33,11 +33,6 @@ Class
Description
An abstract base class for incompressible mixtureViscosityModels.
- The strain rate is defined by:
-
- mag(symm(grad(U)))
-
-
SourceFiles
mixtureViscosityModel.C
mixtureViscosityModelNew.C
diff --git a/tutorials/multiphase/driftFluxFoam/ras/dahl/constant/transportProperties b/tutorials/multiphase/driftFluxFoam/ras/dahl/constant/transportProperties
index 837a0af6a33..6e6819c2663 100644
--- a/tutorials/multiphase/driftFluxFoam/ras/dahl/constant/transportProperties
+++ b/tutorials/multiphase/driftFluxFoam/ras/dahl/constant/transportProperties
@@ -26,7 +26,7 @@ sludge
coeff 0.00023143;
exponent 179.26;
- BinghamCoeff 0.00042189;
+ BinghamCoeff 0.0005966;
BinghamExponent 1050.8;
BinghamOffset 0;
diff --git a/tutorials/multiphase/driftFluxFoam/ras/mixerVessel2D/constant/transportProperties b/tutorials/multiphase/driftFluxFoam/ras/mixerVessel2D/constant/transportProperties
index 837a0af6a33..6e6819c2663 100644
--- a/tutorials/multiphase/driftFluxFoam/ras/mixerVessel2D/constant/transportProperties
+++ b/tutorials/multiphase/driftFluxFoam/ras/mixerVessel2D/constant/transportProperties
@@ -26,7 +26,7 @@ sludge
coeff 0.00023143;
exponent 179.26;
- BinghamCoeff 0.00042189;
+ BinghamCoeff 0.0005966;
BinghamExponent 1050.8;
BinghamOffset 0;
diff --git a/tutorials/multiphase/driftFluxFoam/ras/tank3D/constant/transportProperties b/tutorials/multiphase/driftFluxFoam/ras/tank3D/constant/transportProperties
index 6af13fec407..32ba3d76e84 100644
--- a/tutorials/multiphase/driftFluxFoam/ras/tank3D/constant/transportProperties
+++ b/tutorials/multiphase/driftFluxFoam/ras/tank3D/constant/transportProperties
@@ -26,7 +26,7 @@ sludge
coeff 0.00023143;
exponent 0.17926;
- BinghamCoeff 5.5469e-07;
+ BinghamCoeff 7.845e-07;
BinghamExponent 95.25;
BinghamOffset 0;
--
GitLab