[Julian K.]

Hi!

I am simulating the impaction of particles on a moving circular cylinder using Lagrangian particle tracking. I found out, that CFX does not consider the particle's radius while checking if a particle has collided with the cylinder's surface. Thus, the particle is captured, when its centre coordiantes of the coincide with the cylinder surface.
However, my case requires, that the particle is captured when the particle surface and the cylinder surface touch. Hence, I need to implement this simple algorithm:

if (x <= r_part + r_cyl)
then impaction
else no impaction

Where x is the distance between the centre coordiantes of the particle and the cylinder, r_part is the particle's radius and r_cyl is the radius of the cylinder.

I assume, that I need to implement the code via a User FORTRAN. Since I have absolutely no experience with this, I would like you to ask for some help. I found a file called 'pt_termination.F' (CFX\Samples\UserFortran) with the following code:

IF (DIST.GT.0.75) THEN
   FATE = DEAD
ENDIF

I wonder if this is the right point to start for me?! Do you have any other suggestions for me?

Thanks!