Add wrapper function to simulation loop.

[andytorrestb]

Roughly drawing inspiration from this example case.

https://github.com/FlowKit/palabos-examples/blob/master/showCases/cylinder2d/cylinder2d.cpp

[andytorrestb]

do we still need to do this?