zhyang-dev/CFDEMcoupling-PFM
1
0
Fork 0
mirror of https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git synced 2025-12-08 06:37:44 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Update SchillerNaumannDrag.C

Browse Source 
The changes are added.
This commit is contained in:
asanaz
2020-09-16 10:45:51 +02:00
committed by GitHub
parent 2f5f10c6a5
commit aa82dad8be
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/lagrangian/cfdemParticle/subModels/forceModel/SchillerNaumannDrag/SchillerNaumannDrag.C
View File

@ -87,7 +87,7 @@ SchillerNaumannDrag::SchillerNaumannDrag
// read those switches defined above, if provided in dict // read those switches defined above, if provided in dict
forceSubM(0).readSwitches(); forceSubM(0).readSwitches();
particleCloud_.checkCG(false); particleCloud_.checkCG(true);
if (propsDict_.found("scale")) if (propsDict_.found("scale"))
scaleDia_ = scalar(readScalar(propsDict_.lookup("scale"))); scaleDia_ = scalar(readScalar(propsDict_.lookup("scale")));
@ -151,7 +151,7 @@ void SchillerNaumannDrag::setForce() const
Cd = max(0.44,24.0/Rep*(1.0+0.15*pow(Rep,0.687))); Cd = max(0.44,24.0/Rep*(1.0+0.15*pow(Rep,0.687)));
// calc particle's drag // calc particle's drag
drag = 0.125*Cd*rho*M_PI*ds*ds_scaled*magUr*Ur; drag = 0.125*Cd*rho*M_PI*ds*ds*scaleDia_*magUr*Ur*scaleDrag_;
if (modelType_=="B") if (modelType_=="B")
drag /= voidfraction; drag /= voidfraction;