I am looking for a simple method to compute potential flow (then non-viscous) around an obstacle. The method I am looking for must NOT use coordinate transformation, panel methods, finite element methods; it should work just on a standard cartesian grid and it should preferably use the finite differences method. There is no need for it to be efficient or fast, it should only be SIMPLE. It is only for demonstration purposes. Thank you.
I looked at the Shortley-Weller algorithm, but I didn't understand what should I do exactly... can it be used to impose Neumann conditions on the obstacle boundary (even an asymmetric arbitrary boundary)? If yes, how can it be done?