Reply To: Usage of ParticleDynamics2D::checkAndRemoveEscaped() in dkt2d.cpp
Dear Junwei Guo,
as forwarded from our developer Davide Dapelo, who implemented the function:
ParticleDynamics2D::checkAndRemoveEscaped() removes a particle IF the intersection between the particle’s indicator and the domain’s inner layer is non-empty.
Most likely you misunderstood the function of indicatorF in ParticleDyamics2D’s constructor, since you defined it as the indicator of the domain’s bottom layer (colBase), BUT in fact it should be the whole domain. Following your definition, colBase’s inner layer is just a zero-height layer, and therefore its intersection with any particle indicator is always empty.
The solution to this problem should simply be to change colBase with a whole domain indicator:
Vector orig0(0,0);
Vector ext0(lengthX,lengthY);
If that does not work (but it should), you should define a tighter indicator.
Best wishes,
Nicolas