From cdb78d6741c75f51972fcd5d91090f13c2140bb4 Mon Sep 17 00:00:00 2001 From: henry Date: Tue, 20 May 2008 18:05:07 +0100 Subject: [PATCH] There's a bug in Schaeffer's expression for the frictional viscosity, which was pointed out by Rasmus. The alpha[celli] term in muf should not be there. --- .../Schaeffer/SchaefferFrictionalStress.C | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModel/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C b/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModel/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C index 5f6dcd8735..b3ddd44ff2 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModel/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModel/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C @@ -71,8 +71,8 @@ frictionalPressure const dimensionedScalar& eta, const dimensionedScalar& p ) const -{ - return +{ + return dimensionedScalar("1e24", dimensionSet(1, -1, -2, 0, 0), 1e24) *pow(Foam::max(alpha - alphaMinFriction, scalar(0)), 10.0); } @@ -89,7 +89,7 @@ frictionalPressurePrime const dimensionedScalar& p ) const { - return + return dimensionedScalar("1e25", dimensionSet(1, -1, -2, 0, 0), 1e25) *pow(Foam::max(alpha - alphaMinFriction, scalar(0)), 9.0); } @@ -129,8 +129,8 @@ Foam::tmp Foam::SchaefferFrictionalStress::muf { if (alpha[celli] > alphaMax.value()-5e-2) { - muf_[celli] = - 0.5*alpha[celli]*pf[celli]*sin(phi.value()) + muf_[celli] = + 0.5*pf[celli]*sin(phi.value()) /( sqrt(1.0/6.0*(sqr(D[celli].xx() - D[celli].yy()) + sqr(D[celli].yy() - D[celli].zz())