Add a IBM + cylinder example

[mathrack]

The code now contains the capability for immersed boundary without reconstruction. One example is needed to illustrate this capability.