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

Correct limit for uDyn, boundary conditions.

Browse Source
This commit is contained in:
Thomas Lichtenegger
2016-07-14 14:58:09 +02:00
parent 57e19765e0
commit 45e11369df
1 changed files with 5 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -437,6 +437,8 @@ void FinesFields::updateDragCoeff()
Cd = 0.44;
DragCoeff_.boundaryField()[patchI][faceI] = Cd * beta.boundaryField()[patchI][faceI];
}
DragCoeff_ = max( DragCoeff_, dimensionedScalar("SMALL", dimensionSet(1,-3,-1,0,0), SMALL) );
}
@ -497,7 +499,7 @@ void FinesFields::updateUDyn()
}
}
forAll(uDyn_.boundaryField(), patchI)
/* forAll(uDyn_.boundaryField(), patchI)
forAll(uDyn_.boundaryField()[patchI], faceI)
{
scalar mU(mag(U_.boundaryField()[patchI][faceI]));
@ -507,8 +509,8 @@ void FinesFields::updateUDyn()
uDyn_.boundaryField()[patchI][faceI] *= mU / muDyn;
}
}
// uDyn_ = 0.2*U_;
*/
uDyn_.correctBoundaryConditions();
alphaDyn_.correctBoundaryConditions();
massFluxDyn_ = rhoFine_ * alphaDyn_ * uDyn_;
}