Learning Forum

[kgardner]

Hello,

I am attempting to implement my own time stepping particle tracker in MATLAB and was able to do this but without any walls in the domain.  I used the Force Balance equation given in the Ansys Documentation:

This is the general force balance for a particle within a moving fluid but I actually only used the Drag force term without gravitation and other external forces present.  This Documentation then describes how we can discretize this equation through a number of different techniques.  The technique that I chose was the "analytic discretization scheme":

I was able to implement this into MATLAB just fine but when I try to add a wall in the domain I am getting confused about what the documentation says.  I am attempting to use the adhesive impact model and am confused on the explanation:

I am confused on the boundary condition term here in particular when it says "For a given boundary condition, the post-collision velocities,  and , can be computed by solving Equation 16–217." I am wondering if anyone can explain this better to me and how to implement this boundary condition to the previous discretization equations for the particle velocity and particle position.  I believe I can write an if statement in MATLAB and say if the particle is in the region of the wall then I insert a new velocity, but I need the velocity expression from hitting the wall.

Thanks,

Karl Gardner