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

Some testing on the turbulent dispersion model.

Browse Source
This commit is contained in:
tlichtenegger
2020-11-03 09:51:28 +01:00
parent 25efaffdd2
commit 57e98b2f0e
2 changed files with 49 additions and 14 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -58,6 +58,7 @@ turbulentDispersion::turbulentDispersion
forceModel(dict,sm),
propsDict_(dict.subDict(type + "Props")),
interpolate_(propsDict_.lookupOrDefault<bool>("interpolation", false)),
mesh_(sm.mesh()),
ignoreCellsName_(propsDict_.lookupOrDefault<word>("ignoreCellsName","none")),
ignoreCells_(),
existIgnoreCells_(true),
@ -88,17 +89,39 @@ turbulentDispersion::turbulentDispersion
else existIgnoreCells_ = false;
// define a field to indicate if a cell is next to boundary
// label cellI = -1;
// // set wall indicator field
// const fvMesh& mesh = sm.mesh();
// forAll(mesh.boundary(),patchI)
// {
// forAll(mesh.boundary()[patchI], faceI)
// {
// cellI = patch().faceCells()[faceI];
// wallIndicatorField_[cellI] = 1.0;
// }
// }
label cellI = -1;
// set wall indicator field
Info << "Setting wall indicator field." << endl;
forAll(mesh_.boundary(),patchI)
{
word patchName = mesh_.boundary()[patchI].name();
if (patchName.rfind("procBoundary",0) == 0) continue;
Info << "patch = " << mesh_.boundary()[patchI].name() << endl;
forAll(mesh_.boundary()[patchI], faceI)
{
cellI = mesh_.boundary()[patchI].faceCells()[faceI];
wallIndicatorField_[cellI] = 1.0;
// testing
label patchID = -1;
label faceind = -1;
const cell& faces = mesh_.cells()[cellI];
Pout << "cellI = " << cellI << " has faces = " << faces << endl;
forAll (faces, faceJ) // loop over all faces in cellI
{
faceind = faces[faceJ];
patchID = mesh_.boundaryMesh().whichPatch(faceind);
if (patchID < 0) continue;
vector faceINormal = mesh_.Sf()[faceind] / mesh_.magSf()[faceind] ;
// Pout << " faceind = " << faceind << "at " << mesh_.C()[cellI] << ", faceINormal = " << faceINormal << endl ;
}
// testing done
}
}
wallIndicatorField_.write();
// make sure this is the last force model in list so that fluid velocity does not get overwritten
label numLastForceModel = sm.nrForceModels();
@ -165,9 +188,19 @@ void turbulentDispersion::setForce() const
// if particles are pushed through walls, the velocity fluctuations may be regulated here
// check if cell is adjacent to wall and invert corresponding component
// if (wallIndicatorField_[cellI] > 0.5)
// {
// }
label patchID = -1;
if (wallIndicatorField_[cellI] > 0.5)
{
const cell& faces = mesh_.cells()[cellI];
forAll (faces, faceI) // loop over all faces in cellI
{
patchID = mesh_.boundaryMesh().whichPatch(faceI);
if (patchID < 0) continue;
vector faceINormal = mesh_.Sf()[faceI] / mesh_.magSf()[faceI] ;
Info << " faceI = " << faceI << "at " << mesh_.C()[cellI] << ", faceINormal = " << faceINormal << endl ;
}
}
for(int j=0;j<3;j++)
{

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

@ -49,6 +49,8 @@ protected:
bool interpolate_;
const fvMesh& mesh_;
// ignore fluctuations in region
word ignoreCellsName_;