Option in cfdemCloud to push particle velocities to LIGGGHTS.

This commit is contained in:
Thomas Lichtenegger
2021-04-26 16:10:26 +02:00
parent 3e5be0f42e
commit dc36ae2e66
2 changed files with 24 additions and 7 deletions

View File

@ -88,6 +88,8 @@ cfdemCloud::cfdemCloud
getParticleDensities_(couplingProperties_.lookupOrDefault<bool>("getParticleDensities",false)),
getParticleEffVolFactors_(couplingProperties_.lookupOrDefault<bool>("getParticleEffVolFactors",false)),
getParticleTypes_(couplingProperties_.lookupOrDefault<bool>("getParticleTypes",false)),
streamingMode_(couplingProperties_.lookupOrDefault<bool>("streamingMode",false)),
streamingFluc_(couplingProperties_.lookupOrDefault<bool>("streamingFluc",false)),
maxDEMForce_(0.),
phaseInForcesTime_(couplingProperties_.lookupOrDefault<scalar>("phaseInForcesTime",0.0)),
modelType_(couplingProperties_.lookup("modelType")),
@ -438,15 +440,26 @@ void cfdemCloud::getDEMdata()
void cfdemCloud::giveDEMdata()
{
if(forceM(0).coupleForce())
if (forceM(0).coupleForce())
{
dataExchangeM().giveData("dragforce","vector-atom",DEMForces_);
if(impDEMdrag_)
if (!streamingMode_)
{
if(verbose_) Info << "sending Ksl and uf" << endl;
dataExchangeM().giveData("Ksl","scalar-atom",Cds_);
dataExchangeM().giveData("uf","vector-atom",fluidVel_);
dataExchangeM().giveData("dragforce","vector-atom",DEMForces_);
if (impDEMdrag_)
{
if (verbose_) Info << "sending Ksl and uf" << endl;
dataExchangeM().giveData("Ksl","scalar-atom",Cds_);
dataExchangeM().giveData("uf","vector-atom",fluidVel_);
}
}
else
{
dataExchangeM().giveData("vrec","vector-atom",particleConvVel_);
if (streamingFluc_)
{
dataExchangeM().giveData("vfluc","vector-atom",particleFlucVel_);
}
}
}
if(verbose_) Info << "giveDEMdata done." << endl;

View File

@ -105,6 +105,10 @@ protected:
const bool getParticleTypes_;
const bool streamingMode_;
const bool streamingFluc_;
scalar maxDEMForce_;
scalar phaseInForcesTime_;