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

Extended fines deposition model, fixed tutorial case.

Browse Source
This commit is contained in:
tlichtenegger
2019-08-22 11:37:41 +02:00
parent a6be6be3cc
commit e94964e9bf
6 changed files with 90 additions and 44 deletions
Download Patch File Download Diff File Expand all files Collapse all files

100
src/lagrangian/cfdemParticle/subModels/forceModel/Fines/FinesFields.C
View File

@ -45,6 +45,8 @@ FinesFields::FinesFields
propsDict_(dict.subDict("FinesFieldsProps")),
smoothing_(propsDict_.lookupOrDefault<bool>("smoothing",false)),
verbose_(propsDict_.lookupOrDefault<bool>("verbose",false)),
clogKin_(propsDict_.lookupOrDefault<bool>("kineticClogging",false)),
clogStick_(propsDict_.lookupOrDefault<bool>("stickyClogging",false)),
velFieldName_(propsDict_.lookupOrDefault<word>("velFieldName","U")),
U_(sm.mesh().lookupObject<volVectorField> (velFieldName_)),
voidfractionFieldName_(propsDict_.lookupOrDefault<word>("voidfractionFieldName","voidfraction")),
@ -221,33 +223,65 @@ FinesFields::FinesFields
alphaDynMax_(0.1),
alphaMax_(readScalar(propsDict_.lookup ("alphaMax"))),
critVoidfraction_(readScalar(propsDict_.lookup ("critVoidfraction"))),
depRate_(readScalar(propsDict_.lookup ("depRate"))),
depositionLength_(readScalar(propsDict_.lookup ("depositionLength"))),
exponent_(-1.33),
nCrit_(readScalar(propsDict_.lookup ("nCrit"))),
poresizeWidth_(readScalar(propsDict_.lookup ("poresizeWidth"))),
nCrit_(0.0),
poresizeWidth_(0.0),
prefactor_(10.5),
ratioHydraulicPore_(1.5)
ratioHydraulicPore_(1.5),
tauRelease_(readScalar(propsDict_.lookup ("tauRelease"))),
uBind_(0.0),
uMin_(0.001)
{
Sds_.write();
if (propsDict_.found("prefactor"))
{
prefactor_=readScalar(propsDict_.lookup ("prefactor"));
}
if (propsDict_.found("exponent"))
{
exponent_=readScalar(propsDict_.lookup ("exponent"));
}
if (propsDict_.found("dFine"))
{
dFine_.value()=readScalar(propsDict_.lookup ("dFine"));
}
else
{
FatalError <<"Please specify dFine.\n" << abort(FatalError);
}
if (propsDict_.found("diffCoeff"))
{
diffCoeff_.value()=readScalar(propsDict_.lookup ("diffCoeff"));
}
if (propsDict_.found("rhoFine"))
{
rhoFine_.value()=readScalar(propsDict_.lookup ("rhoFine"));
}
else
{
FatalError <<"Please specify rhoFine.\n" << abort(FatalError);
}
if (propsDict_.found("nuAve"))
{
nuAve_.value()=readScalar(propsDict_.lookup ("nuAve"));
}
if (propsDict_.found("alphaDynMax"))
{
alphaDynMax_=readScalar(propsDict_.lookup ("alphaDynMax"));
}
if (clogKin_)
{
nCrit_ = readScalar(propsDict_.lookup ("nCrit"));
poresizeWidth_ = readScalar(propsDict_.lookup ("poresizeWidth"));
}
if (clogStick_)
{
uBind_ = readScalar(propsDict_.lookup ("uBind"));
}
if(verbose_)
{
@ -307,50 +341,52 @@ void FinesFields::update()
void FinesFields::calcSource()
{
if(!clogKin_ & !clogStick_) return;
Sds_.primitiveFieldRef() = 0;
deltaAlpha_.primitiveFieldRef() = 0.0;
scalar f(0.0);
scalar critpore(0.0);
scalar dmean(0.0);
scalar d1(0.0);
scalar d2(0.0);
scalar fKin = 0.0;
scalar fStick = 0.0;
scalar critpore = 0.0;
scalar dmean = 0.0;
scalar d1 = 0.0;
scalar d2 = 0.0;
scalar tauDeposition = 0.0;
forAll(Sds_,cellI)
{
// calculate everything in units auf dSauter
critpore = nCrit_*dFine_.value()/dSauter_[cellI];
// pore size from hydraulic radius
dmean = 2 * (1 - alphaP_[cellI]) / ( (1 + poresizeWidth_*poresizeWidth_/3) * 3 * alphaP_[cellI] );
// Sweeney and Martin, Acta Materialia 51 (2003): ratio of hydraulic to pore throat radius
dmean /= ratioHydraulicPore_;
d1 = dmean * (1 - poresizeWidth_);
d2 = dmean * (1 + poresizeWidth_);
if(clogKin_)
{
// calculate everything in units auf dSauter
critpore = nCrit_*dFine_.value()/dSauter_[cellI];
// pore size from hydraulic radius
dmean = 2 * (1 - alphaP_[cellI]) / ( (1 + poresizeWidth_*poresizeWidth_/3) * 3 * alphaP_[cellI] );
// Sweeney and Martin, Acta Materialia 51 (2003): ratio of hydraulic to pore throat radius
dmean /= ratioHydraulicPore_;
d1 = dmean * (1 - poresizeWidth_);
d2 = dmean * (1 + poresizeWidth_);
f = (critpore*critpore*critpore - d1 * d1 * d1) / (d2 * d2 * d2 - d1 * d1 * d1);
if (f < 0)
{
f = 0.0;
fKin = (critpore*critpore*critpore - d1 * d1 * d1) / (d2 * d2 * d2 - d1 * d1 * d1);
if (fKin < 0) fKin = 0.0;
else if (fKin > 1.0) fKin = 1.0;
}
else if (f > 1.0)
if(clogStick_)
{
f = 1.0;
fStick = 1.0 / ( 1.0 + mag(U_[cellI])/uBind_);
}
// at this point, voidfraction is still calculated from the true particle sizes
deltaAlpha_[cellI] = f * (alphaMax_ - alphaP_[cellI]) - alphaSt_[cellI];
// too much volume occupied: release it (50% per time step)
deltaAlpha_[cellI] = max(fKin,fStick) * (alphaMax_ - alphaP_[cellI]) - alphaSt_[cellI];
// too much volume occupied: release it
if (deltaAlpha_[cellI] < 0.0)
{
Sds_[cellI] = 0.5*deltaAlpha_[cellI];
}
// volume too occupy available: deposit at most 80% of dyn hold up
else if (depRate_ * deltaAlpha_[cellI] > 0.8 * alphaDyn_[cellI])
{
Sds_[cellI] = 0.8 * alphaDyn_[cellI];
Sds_[cellI] = deltaAlpha_[cellI] / tauRelease_;
}
// volume to occupy available
else
{
Sds_[cellI] = depRate_ * deltaAlpha_[cellI];
tauDeposition = depositionLength_ / max(mag(U_[cellI]),uMin_);
Sds_[cellI] = deltaAlpha_[cellI] / tauDeposition;
}
}
}

13
src/lagrangian/cfdemParticle/subModels/forceModel/Fines/FinesFields.H
View File

@ -46,6 +46,10 @@ private:
bool verbose_;
bool clogKin_;
bool clogStick_;
word velFieldName_;
const volVectorField& U_;
@ -90,7 +94,6 @@ private:
volScalarField Sds_;
//volVectorField massFluxDyn_;
surfaceScalarField massFluxDyn_;
volVectorField uDyn_;
@ -111,7 +114,7 @@ private:
scalar critVoidfraction_;
scalar depRate_;
scalar depositionLength_;
scalar exponent_;
@ -123,6 +126,12 @@ private:
scalar ratioHydraulicPore_;
scalar tauRelease_;
scalar uBind_;
scalar uMin_;
void calcSource();
void integrateFields();

11
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/constant/couplingProperties
View File

@ -38,7 +38,7 @@ meshMotionModel noMeshMotion;
regionModel allRegion;
IOModel "basicIO";
IOModel "off";
dataExchangeModel twoWayMPI;//twoWayM2M;//twoWayFiles;//oneWayVTK;//
@ -55,7 +55,7 @@ forceModels
dSauter
Fines
FanningDynFines
ErgunStatFines
BeetstraDragPoly
gradPForce
viscForce
);
@ -140,13 +140,14 @@ viscForceProps
interpolation;
}
ErgunStatFinesProps
BeetstraDragProps
{
fines true;
dFine 0.000388;
velFieldName "U";
granVelFieldName "Us";
densityFieldName "rho";
voidfractionFieldName "voidfraction";
phi 1;
}
FanningDynFinesProps
@ -207,7 +208,7 @@ dividedProps
{
alphaMin 0.25;
scaleUpVol 1.0;
weight 1.0; //1.33;
weight 1.0;
verbose;
}

2
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/system/blockMeshDict
View File

@ -13,7 +13,7 @@ FoamFile
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#include "../../geometry"
#include "../geometry"
convertToMeters 1;
vertices #codeStream

4
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/system/controlDict
View File

@ -59,7 +59,7 @@ functions
volInt
{
type volRegion;
type volFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
writeControl timeStep;
log true;
@ -76,7 +76,7 @@ functions
inflow
{
type surfaceRegion;
type surfaceFieldValue;
libs ("libfieldFunctionObjects.so");
writeControl timeStep;
log true;

4
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/system/fvOptions
View File

@ -37,7 +37,7 @@ source1
{
active yes;
selectionMode all;
Tmin 288;
Tmax 298;
min 288;
max 298;
}
}