Flow around a pair of 2D cylinders using a hybrid Eulerian-Lagrangian solver | Synapse