4D simulation of bubbles

[tfarago]

Motivation

We want to simulate bubbles in a capillary which move in time. All bubbles should move together upwards but some may be faster than others and may even rupture. The bubbles may have different diameters, trajectories and velocities. The aim is to create a simulation of 4D tomography with artifacts stemming from the process under study, i.e. the bubbles will evolve in time and the whole geometry will be rotated to simulate tomographic DAQ.

Key aspects

• some bubbles move faster than 1 px / tomographic rotation -> motion blur
• some bubbles may rupture -> kind of like limited angle tomography

Challenges

• bubbles must not overlap in the 3D space
• bubbles must not collide during motion
• there should be millions of bubbles + phase contrast = over-sampling -> a lot of computations

Sketch

Images are better than words, so (the code is here)