zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

There's a bug in Schaeffer's expression for the frictional viscosity,

Browse Source 
which was pointed out by Rasmus.

The alpha[celli] term in muf should not be there.
This commit is contained in:
henry
2008-05-20 18:05:07 +01:00
parent b2e05717aa
commit cdb78d6741
1 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModel/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C
View File

@ -71,8 +71,8 @@ frictionalPressure
const dimensionedScalar& eta, const dimensionedScalar& eta,
const dimensionedScalar& p const dimensionedScalar& p
) const ) const
{ {
return return
dimensionedScalar("1e24", dimensionSet(1, -1, -2, 0, 0), 1e24) dimensionedScalar("1e24", dimensionSet(1, -1, -2, 0, 0), 1e24)
*pow(Foam::max(alpha - alphaMinFriction, scalar(0)), 10.0); *pow(Foam::max(alpha - alphaMinFriction, scalar(0)), 10.0);
} }
@ -89,7 +89,7 @@ frictionalPressurePrime
const dimensionedScalar& p const dimensionedScalar& p
) const ) const
{ {
return return
dimensionedScalar("1e25", dimensionSet(1, -1, -2, 0, 0), 1e25) dimensionedScalar("1e25", dimensionSet(1, -1, -2, 0, 0), 1e25)
*pow(Foam::max(alpha - alphaMinFriction, scalar(0)), 9.0); *pow(Foam::max(alpha - alphaMinFriction, scalar(0)), 9.0);
} }
@ -129,8 +129,8 @@ Foam::tmp<Foam::volScalarField> Foam::SchaefferFrictionalStress::muf
{ {
if (alpha[celli] > alphaMax.value()-5e-2) if (alpha[celli] > alphaMax.value()-5e-2)
{ {
muf_[celli] = muf_[celli] =
0.5*alpha[celli]*pf[celli]*sin(phi.value()) 0.5*pf[celli]*sin(phi.value())
/( /(
sqrt(1.0/6.0*(sqr(D[celli].xx() - D[celli].yy()) sqrt(1.0/6.0*(sqr(D[celli].xx() - D[celli].yy())
+ sqr(D[celli].yy() - D[celli].zz()) + sqr(D[celli].yy() - D[celli].zz())