small correction/addition:
– The code works fine without MPI, also for multiple Cuboids. The problem only arises with MPI.
– The fringe region sets the inflow, I don’t add any other external force to the fluid voxels.
– I am working this out first in 2D. Later I will go to 3D.