I am trying to model a flow past a cylinder using the projection method. I am a bit confused on how to define the circular cylinder in my computational domain.
I know that the cylinder has the coordinates $(x-a)^2 + (y-b)^2 = r^2$ but what is the best way to define that in terms of x(i,j), y(i,j) coordinates such that I will not have any points inside the cylinder and so that I can easily apply boundary conditions to the cylinder in my computational domain.
I want to avoid doing a coordinates transformation. It seems that the NS eqtns seem to become very complicated.