“…In the following paragraph, the two-dimensional, axisymmetric, numerical setup, and the applied methods to reproduce the experimentally observed stabilization process are outlined. The parallelized, finite difference code CIAO solves the Navier-Stokes equations in the low Mach number limit, utilizing spatial and temporal staggering, together with Crank-Nicolson type time advancement [26,27]. The Poisson equation for pressure is solved by the multi-grid HYPRE solver and the scalar equations are solved via the bounded quadratic upstream interpolation for convective kinematics (BQUICK) scheme, based on the QUICK scheme of Leonard [28].…”