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

Merge pull request #109 from ParticulateFlow/feature/drag_schiller_naumann

Browse Source 
added coarse graining to SchillerNaumann drag
This commit is contained in:
asanaz
2020-09-16 11:44:53 +02:00
committed by GitHub
parent a83d0290c9 aa82dad8be
commit 59d5a1f1fd
2 changed files with 40 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -66,7 +66,9 @@ SchillerNaumannDrag::SchillerNaumannDrag
propsDict_(dict.subDict(typeName + "Props")), propsDict_(dict.subDict(typeName + "Props")),
verbose_(propsDict_.found("verbose")), verbose_(propsDict_.found("verbose")),
velFieldName_(propsDict_.lookup("velFieldName")), velFieldName_(propsDict_.lookup("velFieldName")),
U_(sm.mesh().lookupObject<volVectorField> (velFieldName_)) U_(sm.mesh().lookupObject<volVectorField> (velFieldName_)),
scaleDia_(1.),
scaleDrag_(1.)
{ {
//Append the field names to be probed //Append the field names to be probed
particleCloud_.probeM().initialize(typeName, typeName+".logDat"); particleCloud_.probeM().initialize(typeName, typeName+".logDat");
@ -85,7 +87,12 @@ 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"))
scaleDia_ = scalar(readScalar(propsDict_.lookup("scale")));
if (propsDict_.found("scaleDrag"))
scaleDrag_ = scalar(readScalar(propsDict_.lookup("scaleDrag")));
} }
@ -99,6 +106,16 @@ SchillerNaumannDrag::~SchillerNaumannDrag()
void SchillerNaumannDrag::setForce() const void SchillerNaumannDrag::setForce() const
{ {
if (scaleDia_ > 1.0)
{
Info << "SchillerNaumann using scale = " << scaleDia_ << endl;
}
else if (particleCloud_.cg() > 1.0)
{
scaleDia_ = particleCloud_.cg();
Info << "SchillerNaumann using scale from liggghts cg = " << scaleDia_ << endl;
}
#include "setupProbeModel.H" #include "setupProbeModel.H"
const volScalarField& nufField = forceSubM(0).nuField(); const volScalarField& nufField = forceSubM(0).nuField();
@ -117,6 +134,7 @@ void SchillerNaumannDrag::setForce() const
vector Us = particleCloud_.velocity(index); vector Us = particleCloud_.velocity(index);
vector Ur = U_[cellI]-Us; vector Ur = U_[cellI]-Us;
scalar ds = 2*particleCloud_.radius(index); scalar ds = 2*particleCloud_.radius(index);
scalar ds_scaled = ds/scaleDia_;
scalar nuf = nufField[cellI]; scalar nuf = nufField[cellI];
scalar rho = rhoField[cellI]; scalar rho = rhoField[cellI];
scalar voidfraction = particleCloud_.voidfraction(index); scalar voidfraction = particleCloud_.voidfraction(index);
@ -127,30 +145,33 @@ void SchillerNaumannDrag::setForce() const
if (magUr > 0) if (magUr > 0)
{ {
// calc particle Re Nr // calc particle Re Nr
Rep = ds*magUr/nuf; Rep = ds_scaled*magUr/nuf;
// calc fluid drag Coeff // calc fluid drag Coeff
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*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;
} }
if(verbose_ && index >100 && index <102) if(verbose_ && index >=100 && index <102)
{ {
Info << "index = " << index << endl; Pout << "cellI = " << cellI << endl;
Info << "Us = " << Us << endl; Pout << "index = " << index << endl;
Info << "Ur = " << Ur << endl; Pout << "Ufluid = " << U_[cellI] << endl;
Info << "ds = " << ds << endl; Pout << "Us = " << Us << endl;
Info << "rho = " << rho << endl; Pout << "Ur = " << Ur << endl;
Info << "nuf = " << nuf << endl; Pout << "ds = " << ds << endl;
Info << "voidfraction = " << voidfraction << endl; Pout << "ds/scale = " << ds_scaled << endl;
Info << "Rep = " << Rep << endl; Pout << "rho = " << rho << endl;
Info << "Cd = " << Cd << endl; Pout << "nuf = " << nuf << endl;
Info << "drag = " << drag << endl; Pout << "voidfraction = " << voidfraction << endl;
Pout << "Rep = " << Rep << endl;
Pout << "Cd = " << Cd << endl;
Pout << "drag = " << drag << endl;
} }
//Set value fields and write the probe //Set value fields and write the probe

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

@ -66,6 +66,10 @@ private:
const volVectorField& U_; const volVectorField& U_;
mutable scalar scaleDia_;
mutable scalar scaleDrag_;
public: public:
//- Runtime type information //- Runtime type information