zhyang-dev/CFDEMcoupling-PFM
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Lookup turbulent kinetic energy in execution loop.

Browse Source
This commit is contained in:
tlichtenegger
2020-11-02 10:27:05 +01:00
parent 7eacdef7eb
commit 25efaffdd2
2 changed files with 5 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
src/lagrangian/cfdemParticle/subModels/forceModel/turbulentDispersion/turbulentDispersion.C
View File

@ -73,8 +73,6 @@ turbulentDispersion::turbulentDispersion
sm.mesh(),
dimensionedScalar("zero", dimensionSet(0,0,0,0,0,0,0), 0.0)
),
turbKinetcEnergyFieldName_(propsDict_.lookupOrDefault<word>("turbKinetcEnergyFieldName","k")),
turbKinetcEnergy_(sm.mesh().lookupObject<volScalarField> (turbKinetcEnergyFieldName_)),
minTurbKinetcEnergy_(propsDict_.lookupOrDefault<scalar>("minTurbKinetcEnergy", 0.0)),
voidfractionFieldName_(propsDict_.lookupOrDefault<word>("voidfractionFieldName","voidfraction")),
voidfraction_(sm.mesh().lookupObject<volScalarField> (voidfractionFieldName_)),
@ -130,13 +128,15 @@ bool turbulentDispersion::ignoreCell(label cell) const
void turbulentDispersion::setForce() const
{
const volScalarField turbKinetcEnergy(particleCloud_.turbulence().k());
vector position(0,0,0);
scalar k = 0.0;
vector flucU(0,0,0);
label cellI = 0;
interpolationCellPoint<scalar> turbKinetcEnergyInterpolator_(turbKinetcEnergy_);
interpolationCellPoint<scalar> turbKinetcEnergyInterpolator_(turbKinetcEnergy);
for(int index = 0;index < particleCloud_.numberOfParticles(); ++index)
{
@ -156,7 +156,7 @@ void turbulentDispersion::setForce() const
}
else
{
k = turbKinetcEnergy_[cellI];
k = turbKinetcEnergy[cellI];
}
if (k < minTurbKinetcEnergy_) k = minTurbKinetcEnergy_;
@ -175,7 +175,7 @@ void turbulentDispersion::setForce() const
}
}
}
}
}
}
vector turbulentDispersion::unitFlucDir() const

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

@ -58,10 +58,6 @@ protected:
mutable volScalarField wallIndicatorField_;
word turbKinetcEnergyFieldName_;
const volScalarField& turbKinetcEnergy_;
scalar minTurbKinetcEnergy_;
word voidfractionFieldName_;