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

JohnsonJacksonFrictionalStress: alphaDeltaMax -> alphaDeltaMin (it is the lower limit)

Browse Source
This commit is contained in:
Henry
2015-03-18 15:40:39 +00:00
parent 2e35cc76e8
commit fa2721ffd5
4 changed files with 8 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/frictionalStressModel/JohnsonJackson/JohnsonJacksonFrictionalStress.C
View File

@ -62,7 +62,7 @@ JohnsonJackson
eta_("eta", dimless, coeffDict_.lookup("eta")),
p_("p", dimless, coeffDict_.lookup("p")),
phi_("phi", dimless, coeffDict_.lookup("phi")),
alphaDeltaMax_("alphaDeltaMax", dimless, coeffDict_.lookup("alphaDeltaMax"))
alphaDeltaMin_("alphaDeltaMin", dimless, coeffDict_.lookup("alphaDeltaMin"))
{
phi_ *= constant::mathematical::pi/180.0;
}
@ -89,7 +89,7 @@ frictionalPressure
return
Fr_*pow(max(alpha1 - alphaMinFriction, scalar(0)), eta_)
/pow(max(alphaMax - alpha1, alphaDeltaMax_), p_);
/pow(max(alphaMax - alpha1, alphaDeltaMin_), p_);
}
@ -107,7 +107,7 @@ frictionalPressurePrime
eta_*pow(max(alpha1 - alphaMinFriction, scalar(0)), eta_ - 1.0)
*(alphaMax-alpha1)
+ p_*pow(max(alpha1 - alphaMinFriction, scalar(0)), eta_)
)/pow(max(alphaMax - alpha1, alphaDeltaMax_), p_ + 1.0);
)/pow(max(alphaMax - alpha1, alphaDeltaMin_), p_ + 1.0);
}
@ -135,7 +135,7 @@ bool Foam::kineticTheoryModels::frictionalStressModels::JohnsonJackson::read()
phi_.read(coeffDict_);
phi_ *= constant::mathematical::pi/180.0;
alphaDeltaMax_.read(coeffDict_);
alphaDeltaMin_.read(coeffDict_);
return true;
}

4
applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/frictionalStressModel/JohnsonJackson/JohnsonJacksonFrictionalStress.H
View File

@ -69,8 +69,8 @@ class JohnsonJackson
//- Angle of internal friction
dimensionedScalar phi_;
//- Limit of alphaMax - alpha1 (~5e-2 or higher)
dimensionedScalar alphaDeltaMax_;
//- Lower limit for (alphaMax - alpha1)
dimensionedScalar alphaDeltaMin_;
public:

2
tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/constant/turbulenceProperties.particles
View File

@ -45,7 +45,7 @@ RAS
eta 2;
p 5;
phi 28.5;
alphaDeltaMax 0.05;
alphaDeltaMin 0.05;
}
}

2
tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/constant/turbulenceProperties.particles
View File

@ -44,7 +44,7 @@ RAS
eta 2;
p 5;
phi 28.5;
alphaDeltaMax 0.05;
alphaDeltaMin 0.05;
}
}