Shallow Water Analogue of the Standing Accretion Shock Instability: Experimental Demonstration and a Two-Dimensional Model

  • Thierry Foglizzo, Frédéric Masset, Jérôme Guilet, and Gilles Durand, Laboratoire AIM Paris-Saclay, CEA/Irfu Université Paris-Diderot CNRS/INSU, F-91191 Gif sur Yvette, France, ICF, UNAM, Cuernavaca, Morelos, 62210, Mexico, and DAMTP, University of Cambridge, Centre for Mathematical Sciences, Cambridge CB3 0WA, United Kingdom
Phys. Rev. Lett. 108, 051103
Like the classical hydraulic jump in a kitchen sink (upper left), the SWASI experiment (upper right and bottom) involves a hydraulic jump associated to the deceleration of a radial flow of water. Water is injected inward from an annular injection reservoir (R) along a hyperbolic potential well, and evacuated through a vertical cylinder (C), whose walls mimic the surface of the neutron star. A pump (P) distributes collected water. The lower picture illustrates the stationary and axisymmetric character of water injection in a stable configuration.Like the classical hydraulic jump in a kitchen sink (upper left), the SWASI experiment (upper right and bottom) involves a hydraulic jump associated to the deceleration of a radial flow of water. Water is injected inward from an annular injection res... Show more

The physics of stellar core collapse involves complex ingredients such as nuclear physics, neutrino interactions, multidimensional hydrodynamics, and general relativity. Current understanding is based on simplified formulations which make the problem numerically tractable at least, if not physically intuitive. A recent breakthrough occurred on the hydrodynamical side with the recognition of a new instability mechanism [1–7] taking place in the inner 200 km of the collapsing core, during a few hundred milliseconds before the explosion, while a spherical accretion shock stalls above the surface of the proto-neutron star. This shock was found to be unstable with respect to global oscillations, with a period of a few tens of milliseconds. An explosion is launched when sufficient neutrino energy diffusing out of the dense central region has been intercepted by the post-shock accreting matter [8]. When the explosion proceeds, the large scale asymmetry resulting from this standing accretion shock instability (SASI) imposes on the neutron star a velocity kick mediated by gravitational forces [9–12], and can also affect its spin in a direction surprisingly opposite to the global rotation of the shock wave [13,14]. The pulsar spin period that could be accounted for by the spiral mode of SASI is still debated [12–15]. These effects are likely to depend on the progenitor mass and rotation rate, the physical ingredients of the numerical model, and could also be sensitive to numerical limitations such as numerical resolution. We shed new light on the physical nature of these processes by adopting an experimental approach for the first time. The outcome is strikingly similar to the astrophysical results obtained through numerical simulations in an adiabatic approximation [13]. The simplicity of the experiment in the familiar environment of shallow water physics should allow physicists to build up physical intuition about this instability.

A shallow water analogue to a shock instability.—The experimental set up named SWASI is designed to produce a shallow water analogue of a shock instability. It is based on the analogy between acoustic waves in a gas and surface gravity waves in shallow water [16]. Their velocity c when the depth H is smaller than their wavelength is defined by c2=gH, where g is the gravitational acceleration in the laboratory. The hyperbolic shape Hgrav(r) of the surface on which the shallow fluid flows is a formal analogue to the gravitational potential Φ=gHgrav, dominated by the mass of the central proto-neutron star. The system of equations defining this flow in the inviscid 2D approximation can be directly compared to the equations describing an isentropic gas with an adiabatic index γ=2 [17], where the compressible gas density plays the same role as the depth of the fluid layer in the experiment:

The hydraulic jump governed by a spiral mode displays a rotating triple point both in the experiment (left) and in the numerical simulation of 2D shallow water equations (right). The altitude of the free surface (upper right) and the vorticity (lower right) are shown. The same shape and dynamics are observed in astrophysical simulations [13]. The angular momentum in the accreted flow is visualized by a horizontal bar in the experiment (lower left). It spins in the direction opposite to the hydraulic jump (movie M4 [21]). The vorticity trails shown in the numerical simulation (lower right) of the experiment illustrate these counterrotating motions (movie M6 [21]).The hydraulic jump governed by a spiral mode displays a rotating triple point both in the experiment (left) and in the numerical simulation of 2D shallow water equations (right). The altitude of the free surface (upper right) and the vorticity (lower... Show more

The Froude number Fr=|v|/c comparing the fluid velocity v to the wave velocity plays the same role as the Mach number in a gas. Shock waves are analogue to hydraulic jumps [16]. Both can be idealized as sharp discontinuities where the mass flux and momentum flux are conserved. In contrast with the kitchen sink experiment [18] (Fig. 1, upper left), the fluid is injected inward in the SWASI experiment (Fig. 1, bottom). The inner boundary is made of a vertical pipe such that water accumulating in the gravitational potential is continuously extracted by spilling over its upper edge (Fig. 1, upper right). The height of the upper edge of the inner cylinder defines a pressure threshold over which water is efficiently evacuated. It is adjusted vertically in order to choose the radius of the hydraulic jump ( rjp20cm). Conveniently simple, this inner boundary could be interpreted as an analogue of a cooling process such as neutrino emission close to the neutrinosphere, which decreases the energy of the gas when the density and temperature are high enough. The experiment size was chosen large enough in order to minimize the effects of viscosity: the viscous drag is negligible except in the thin layer ahead of the jump. The injection radius is 32 cm, and the radius of the inner cylinder is r*=4cm. The shape of the gravitational potential is described by dHgrav/dr=(5.6cm/r)2. Accreted water is pumped back into the annular reservoir through 16 pipes. A layer of sand dampens any inhomogeneity before injection.

Despite its simple set up, this experiment is expected to capture some hydrodynamical properties of the accreting gas in the equatorial plane of the stellar core [14,19], particularly those observed in the adiabatic approximation [13]. The shallow water model can be scaled using the jump radius rjp as distance unit, the free fall velocity vff(2gHgravjp)1/2 as velocity unit, and the free fall time scale tffjprjp/vff as time unit. The dimensionless solution depends on three parameters only [17]: the relative size of the inner boundary r*/rjp, the preshock velocity v1/vff, and the Froude number Fr1 ahead of the shock. The experimental results can be scaled to astrophysical proportions by using the ratio rsh/rjp106 for distances, and the ratio tffsh/tffjp for time scales:

The oscillation period measured in the experiment, shown here as a function of the jump radius, is well captured by both the perturbative analysis and the numerical simulation. Error bars relative to the measurement of the jump radius (±5mm) and the oscillation period (±3.5%) are indicated. The measured flow rate is 1L/s (±5%) and the injection slit size is 0.74 mm (±0.02mm).The oscillation period measured in the experiment, shown here as a function of the jump radius, is well captured by both the perturbative analysis and the numerical simulation. Error bars relative to the measurement of the jump radius (±5mm) and th... Show more

where MNS1.2Msol is the mass of the proto-neutron star and G is the gravitational constant.

Surface gravity waves and advected vorticity perturbations present in shallow water are directly comparable to acoustic and vorticity waves in a compressible gas. These are a possible source of an unstable cycle similar to the advective-acoustic cycle in SASI [2–7] and also comparable to the vortical-acoustic cycle seen in a shocked isothermal gas [20].

The stability of the flow depends on the strength of the hydraulic jump measured by the prejump Froude number Fr1, and the prejump velocity v1 measured in units of the free fall velocity vff along the potential surface. The growth rate of the dominant m=1 mode, deduced from the perturbative analysis of the 2D shallow water model, is displayed by contour lines for 0, 0.035, 0.07, 0.105, 0.14, 0.175 (tffjp)1, for a jump radius rjp/r*=4.9 (tffjp=0.35s). The symbols “+” indicate the set of parameters used in experiments showing a robust m=1 instability, the filled square corresponds to the parameters of Fig. 3. Dashed lines of constant flow rate (0.5, 1.5 and 2.5L/s) from bottom to top correspond to a Reynolds number equal to 408, 1224, and 2040, respectively, measured at the jump radius.The stability of the flow depends on the strength of the hydraulic jump measured by the prejump Froude number Fr1, and the prejump velocity v1 measured in units of the free fall velocity vff along the potential surface. The growth rate of the dominan... Show more

Symmetry breaking.—Indeed, as the water flux is increased, a large scale instability sets in through growing oscillations of the hydraulic jump (Fig. 2, left). This instability is also observed in numerical simulations [21] of the experimental setup in the 2D shallow water approximation (Fig. 2, right). The visual resemblance with astrophysical simulations is supported theoretically by the formal similarity between the set of shallow water equations and the set of adiabatic gas equations used in [13]. A perturbative analysis of the shallow water equations [17] reveals a SASI-like instability dominated by the global mode m=1 in most of the parameter space (Fig. 4). The oscillation period measured in the experiment is in excellent agreement with both the perturbative analysis and the numerical simulations (Fig. 3). Using the scaling factor in Eq. (3), an oscillation frequency of 3 s and a growth rate of 0.2s-1 in the experiment for v1/vff1, rjp/r*=5 and Fr14 would correspond to an oscillation period of 42 ms and a growth time of 70 ms on an astrophysical scale, which are comparable to the values measured in core-collapse models (e.g., Fig. 17 in [4]).

Like 3D astrophysical simulations [13,14], the experiment displays both the sloshing and the spiral modes (movies M1, M2, M3 [21]). These are easily triggered depending on the random initial perturbations, as expected, since they have the same linear growth rates. However they interact during the nonlinear phase of their evolution, ultimately favoring a single right or left spiral mode even though the injected flow contains no angular momentum. This nonlinear behavior, first noted in astrophysical numerical simulations [13], is clearly observed in the experiment if the instability is vigorous (movies M2, M3 and M5 [21]). The shape and dynamical evolution of the hydraulic jump in the nonlinear regime is remarkably akin to the astrophysical numerical simulations [13,14,19,22], with the formation of a triple point when the spiral mode reaches nonlinear amplitudes (Fig. 2, movie M6 [21]). Another striking similarity is the fact that the angular momentum of the accreted fluid is opposite to the direction of rotation of the hydraulic jump. This is made visible in the experiment by using a light horizontal bar freely rotating on a vertical axis, with both ends bent vertically and immersed in the flowing water. The hydraulic jump and the horizontal bar systematically rotate in opposite directions (Fig. 2 and movies M4, M5 [21]).

Effect of the bottom drag for a viscous fluid.—The Reynolds number in the experiment can be defined using the viscosity ν of water, the fluid velocity v and depth H. It is smallest in the outer part of the experiment and inversely scales with radius: ReHv/ν=Q/2πνr where Q is the injected flow rate. According to Fig. 4, the Reynolds number at the shock radius is in the range 500–2000 within the parameters of the experiment, which supports the laminar approximation for the viscous drag [23], proportional to v/H2. This viscous drag is non-negligible only in the most shallow layer ahead of the hydraulic jump. The strength of the viscous drag has been estimated experimentally in a stationary supercritical flow ( Fr>1) by measuring the radial profile of the potential surface and the free surface of the fluid between R=6cm and R=31cm for an injection slit size of 1.3 mm and a flow rate of 1.15L/s. The corresponding velocity profile was best fitted with an effective viscosity coefficient ν¯=0.03cm2/s, such that the viscous drag is written as an additional term ν¯v/H2 in Eq. (2). The perturbative analysis [17] of the 2D shallow water model assumes that this parameter is a constant.

In order to stress the analogy between the dynamics of the post-jump region in the experiment and the dynamics of the post-shock region in the astrophysical simulations, the flow parameters were chosen according to their values at the shock (velocity and Froude number) rather than at injection. Measuring the prejump velocity in units of the inviscid free fall velocity in Fig. 4, the contour lines showing the growth rate of the instability for ν¯=0.03cm2/s are well approximated by those obtained for an inviscid flow.

Modeling approximations.—The hydraulic jump has been modeled as a simple discontinuity where energy is dissipated through a viscous roller of negligible radial extension. In reality, depending on the strength of the jump, small scale waves can be generated which may interfere with the global instability for some parameters. Intermittent behavior of hydraulic jumps is well known in the range 2.5<Fr<4.5, while steadier jumps are expected in the range 4.5<Fr<9 [24]. The latter is more typical of preshock Mach numbers relevant to core-collapse studies. The experiments showed possible interferences between the intermittent character of the jump and the global instability for Fr<3. Figure 4 shows the set of most robust experiments performed for Fr>3.3. The transition to turbulence, expected around Re=2000, would require a different prescription for the viscous drag (proportional to v2/H). Viscosity in the horizontal plane and surface tension have been neglected. The vertical structure of the flow has been idealized in order to obtain a 2D description of the experiment. A more accurate modeling of this vertical structure could be obtained by introducing the Coriolis and Boussinesq correction coefficients [25] for the kinetic energy and the momentum equations, in the St. Venant system and in the boundary conditions. Such quantitative improvement would be at the expense of the simplicity of the analogy with gas equations.

A new tool for supernova physics.—Drawing a parallel between stellar explosions and shallow water physics is unexpected, but is deeply rooted in the universality of the laws of fluid mechanics. The complementarity between the experimental and numerical approaches is beneficial because their inherent limitations are different. The experiment together with its idealized modeling demonstrated the robustness of SASI in a converging flow, even in a fluid without entropy gradients. This emphasizes implicitly the importance of vorticity perturbations, like in the vortical-acoustic cycle in an isothermal flow [20]. This new tool will be used to characterize nonlinear effects determining both the saturation amplitude and the interaction between right and left spiral modes. The saturation amplitude in the inviscid approximation is expected to be governed by the parasitic growth of the Kelvin-Helmholtz instability [26], apparent in the numerical simulations (Fig. 2 lower right). Furthermore, a global rotation of the experiment will be introduced in order to mimic stellar spin. This is expected to accelerate the growth of the spiral mode rotating in the same direction as the progenitor [13,27] and enhance the influence of SASI on the angular momentum budget of the accretor.

This shallow water analogue of a shock instability belongs to the rare category of physical processes which are astrophysically relevant, accessible on Earth at normal temperature and pressure, and low cost. It is both a new research instrument for supernova science and an original tool for public outreach.


  1. J. Blondin, A. Mezzacappa, and C. DeMarino, Astrophys. J. 584, 971 (2003).
  2. N. Ohnishi, K. Kotake, and S. Yamada, Astrophys. J. 641, 1018 (2006).
  3. T. Foglizzo, P. Galletti, L. Scheck, and H.-T. Janka, Astrophys. J. 654, 1006 (2007).
  4. L. Scheck, H.-T. Janka, T. Foglizzo, and K. Kifonidis, Astron. Astrophys. 477, 931 (2008).
  5. T. Foglizzo, Astrophys. J. 694, 820 (2009).
  6. R. Fernandez and C. Thompson, Astrophys. J. 697, 1827 (2009).
  7. J. Guilet and T. Foglizzo, arXiv:1112.1427 [Mon. Not. R. Astron. Soc. (to be published)].
  8. A. Marek and H.-T. Janka, Astrophys. J. 694, 664 (2009).
  9. L. Scheck, T. Plewa, H. T. Janka, K. Kifonidis, and E. Müller, Phys. Rev. Lett. 92, 011103 (2004).
  10. L. Scheck, K. Kifonidis, H.-T. Janka, and E. Müller, Astron. Astrophys. 457, 963 (2006).
  11. J. Nordhaus, T. Brandt, A. Burrows, E. Livne, and C. Ott, Phys. Rev. D 82, 103016 (2010).
  12. A. Wongwathanarat, H.-T. Janka, and E. Müller, Astrophys. J. 725, L106 (2010).
  13. J. Blondin and A. Mezzacappa, Nature (London) 445, 58 (2007).
  14. R. Fernandez, Astrophys. J. 725, 1563 (2010).
  15. E. Rantsiou, A. Burrows, J. Nordhaus, and A. Almgren, Astrophys. J. 732, 57 (2011).
  16. L. Landau and E. Lifchitz, Theoretical Physics 6. Fluid Mechanics (Pergamon Press, New York, 1987), Chap. 108, 2nd ed.
  17. See Supplemental Material at for a description of the set of equations defining the 2D model and the numerical technique used for the simulations; (also
  18. E. J. Watson, J. Fluid Mech. 20, 481 (1964).
  19. J. Blondin and S. Shaw, Astrophys. J. 656, 366 (2007).
  20. T. Foglizzo, Astron. Astrophys. 392, 353 (2002).
  21. See Supplemental Material at for the movies of both the experiment and its numerical modeling; (also
  22. W. Iwakami, K. Kotake, N. Ohnishi, S. Yamada, and K. Sawada, Astrophys. J. 678, 1207 (2008).
  23. L. F. Moody, Trans. ASME 66, 671 (1944).
  24. V. Chow, Open Channel Hydraulics (McGraw-Hill International, New York, 1973).
  25. H. Chanson, The Hydraulics of Open Channel Flow: An Introduction (Elsevier, New York, 2004) Chap. 3.2.2, pp. 27–29, 2nd ed.
  26. J. Guilet, J. Sato, and T. Foglizzo, Astrophys. J. 713, 1350 (2010).
  27. T. Yamasaki and T. Foglizzo, Astrophys. J. 679, 607 (2008).

About the Authors

Image of Thierry Foglizzo
Image of Frédéric Masset
Image of Jérôme Guilet
Image of Gilles Durand

Related Articles

Notes from the Editors: Astroparticle Physics Italian Style

Notes from the Editors: Astroparticle Physics Italian Style

A new film about dark matter and surprises from the Crab nebula were on display in Turin at the International Conference on Topics in Astroparticle and Underground Physics. Read More »

Synopsis: Predicting the Thickness of Sea Ice

Synopsis: Predicting the Thickness of Sea Ice

A new solution to an old equation will make it easier to model the evolution of sea-ice thickness. Read More »

Synopsis: Bright Twins

Synopsis: Bright Twins

An optical-fiber scheme generates a bright “twin beam”—a pair of quantum-correlated beams that could be used in high-precision metrology. Read More »

More Articles