[h.bagherzadeh]

Hi
I have problem in implementation of periodic boundary condition (PBC) for both particles and fluid flow in a 2D channel. My code works correctly without particles but when I add particles to flow, it diverges after a while. I have applied SIMPLE algorithm to solve Naiver stokes equations along with methodology suggested by Patankar (Fully Developed Flow and Heat Transfer in Ducts Having Streamwise-Periodic Variations of Cross-Sectional Area_1977) for PBC. I have no idea what's wrong with my code? I would be grateful if anyone could help.