Plots:

Each run was carried out until a sink particle formed (except the Light ambient case, for which the run was carried out for 5 crossing times to check stability during which time no sink particle was formed). The plots then sampled 5 evenly spaced times during the simulation between t=0 and t=tfinal, where tfinal was the time when a sink particle formed.

1. rho -

Stability Case:

BP Case:

Matched Case:

3 Intermediate cases, log spaced between the Matched and Light cases:

rho_ambient = rho(Rbe)/3.162:

rho_ambient = rho(Rbe)/10:

rho_ambient = rho(Rbe)/31.62: