Weak and strong plumes are simulated in 3D with a non-equilibrium gas–particle model.
Multiphase effects have a secondary role in the analysed large-scale plume dynamics.
Integral model assumptions and the entrainment hypothesis hold for weak plumes.
In the strong plume case, a supersonic fountaining structure is formed.
Atmospheric stratification hinders the achievement of self-similarity.