------------------------------------------------------------
revno: 3910
committer: Raphael Maurin <[email protected]>
timestamp: Sun 2016-06-26 13:42:49 +0200
message:
HydroForceEngine: add a condition to ignore particles fixed on the sideWalls.
modified:
pkg/common/ForceEngine.cpp
--
lp:yade
https://code.launchpad.net/~yade-pkg/yade/git-trunk
Your team Yade developers is subscribed to branch lp:yade.
To unsubscribe from this branch go to
https://code.launchpad.net/~yade-pkg/yade/git-trunk/+edit-subscription
=== modified file 'pkg/common/ForceEngine.cpp'
--- pkg/common/ForceEngine.cpp 2016-05-19 14:19:12 +0000
+++ pkg/common/ForceEngine.cpp 2016-06-26 11:42:49 +0000
@@ -153,6 +153,8 @@
const Real zPos = b->state->pos[2]-zRef;
int Np = floor(zPos/deltaZ); //Define the layer number with 0 corresponding to zRef. Let the z position wrt to zero, that way all z altitude are positive. (otherwise problem with volPart evaluation)
+ if ((b->state->blockedDOFs==State::DOF_ALL)&&(zPos > s->radius)) continue;// to remove contribution from the fixed particles on the sidewalls.
+
// Relative fluid/particle velocity using also the associated fluid vel. fluct.
if ((Np>=0)&&(Np<nCell)){
uRel = Vector3r(vxFluid[Np]+vFluctX[b->id], vFluctY[b->id],vFluctZ[b->id]) - b->state->vel;
_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-dev
More help : https://help.launchpad.net/ListHelp
