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

Corrected the force calculation on the multisphere particle

Browse Source 
Corrected the loop determination for calculating the forcfe acting on a
multispehre particle.
This commit is contained in:
Achuth1992
2019-05-20 16:10:30 +02:00
parent c6815c9156
commit 3eb0fdbc8a
1 changed files with 15 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files

33
src/lagrangian/cfdemParticle/subModels/forceModel/ShirgaonkarIB/ShirgaonkarIB.C
View File

@ -65,8 +65,8 @@ ShirgaonkarIB::ShirgaonkarIB
propsDict_(dict.subDict(typeName + "Props")),
verbose_(false),
twoDimensional_(false),
multisphere_(false), // drag for a multisphere particle
depth_(1),
multisphere_(false), // drag for a multisphere particle
velFieldName_(propsDict_.lookup("velFieldName")),
U_(sm.mesh().lookupObject<volVectorField> (velFieldName_)),
pressureFieldName_(propsDict_.lookup("pressureFieldName")),
@ -171,38 +171,35 @@ void ShirgaonkarIB::calcForce() const
for(int index=0; index< particleCloud_.numberOfParticles(); index++)
{
int prev;
int prev = 0;
for(int subCell=0;subCell<particleCloud_.voidFractionM().cellsPerParticle()[index][0];subCell++)
{
cellID = particleCloud_.cellIDs()[index][subCell];
prev = 0;
if (cellID > -1 && index == 0) // Force on 1st particle
if (cellID > -1) // Force on 1st particle
{
dragMS += h[cellID]*h.mesh().V()[cellID];
if (index == 0) dragMS += h[cellID]*h.mesh().V()[cellID];
}
else if(cellID > -1 && index > 0)
{
// Check for cellID in previous particles
for(int i=index-1; i>= 0; i--)
if(index > 0)
{
for(int j=0; j< particleCloud_.voidFractionM().cellsPerParticle()[i][0];j++)
// Check for cellID in previous particles
for(int i=index-1; i>= 0; i--)
{
label cell = particleCloud_.cellIDs()[i][j];
if(cellID == cell){prev++;}
for(int j=0; j< particleCloud_.voidFractionM().cellsPerParticle()[i][0];j++)
{
label cell = particleCloud_.cellIDs()[i][j];
if(cell > -1 && cell == cellID){prev++;}
}
}
}
if(prev == 0){ dragMS += h[cellID]*h.mesh().V()[cellID];}
}
if(prev == 0){ dragMS += h[cellID]*h.mesh().V()[cellID];}
}
}
Info << "Drag force on particle clump = " << dragMS[0]<<","<<dragMS[1]<<","<<dragMS[2] << endl;
//Info << "Drag force on particle clump = " << dragMS[0] << ", " << dragMS[1] << ", " << dragMS[2] << endl;
printf("Drag force on particle clump = %f, %f, %f\n",dragMS[0],dragMS[1],dragMS[2]);
}