simulating flow past a cylinder

dear all

I’m trying to simulate flow past a cylinder using LAMMPS. I’ve been successful in getting it about say 80% right. the only trouble that I’m facing right now is due to the periodic boundary condition along the x axis due to which the cylinder is influencing the part of flow before encountering it. is there any way in which i can delete the atoms exiting from the cell and create fresh atoms to enter the cell, during each step of the run. Is it possible? I’m attaching the code that I’ve written.


in.flowpastcylinder (1.01 KB)