sim_psy2 function

Simulation of a two-bubble process