The main simulation was made in Realflow.

I used the bin particles to add some air bubbles and plankton.

And finally Fusion to comp all the passes together.

the breakdown