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

Thierry Foglizzo,1 Frdric Masset,2,1 Jrme Guilet,3,1 and Gilles Durand1 e e eo

AIM, IRFU/SAp, CEA Saclay, 91191, France 2 ICF, UNAM, 62210, Mexico 3 DAMTP Cambridge, CB3 0WA, UK (Dated: accepted December 13, 2011)

arXiv:1112.3448v1 [astro-ph.HE] 15 Dec 2011

Despite the sphericity of the collapsing stellar core, the birth conditions of neutron stars can be highly non spherical due to a hydrodynamical instability of the shocked accretion ow. Here we report the rst laboratory experiment of a shallow water analogue, based on the physics of hydraulic jumps. Both the experiment and its shallow water modeling demonstrate a robust linear instability and nonlinear properties of symmetry breaking, in a system which is one million times smaller and about hundred times slower than its astrophysical analogue. PACS numbers: 97.60.Bw, 95.30.Lz, 47.20.-k, 01.50.Pa

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 simplied 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 [17] taking place in the inner 200km of the collapsing core, during a few hundred milliseconds before the explosion, while a spherical accretion shock stalls above the surface of the protoneutron 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 sucient neutrino energy diusing 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 [912], and can also aect 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 [1215]. These eects 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 rst 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

FIG. 1. 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 ow 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 conguration.

surface gravity waves in shallow water [16]. Their velocity c when the depth H is smaller than their wavelength is dened by c2 = gH, where g is the gravitational acceleration in the laboratory. The hyperbolic shape Hgrav (r) of the surface on which the shallow uid ows is a formal analogue to the gravitational potential = gHgrav , dominated by the mass of the central proto-neutron star. The system of equations dening this ow in the inviscid 2D approximation can be directly compared to the equations describing an isentropic gas with an adiabatic index

2 = 2 [17], where the compressible gas density plays the same role as the depth of the uid layer in the experiment: H + (Hv) = 0, (1) t v2 v + ( v) v + + c2 + = 0. (2) t 2 The Froude number Fr = |v|/c comparing the uid 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 ux and momentum ux are conserved. In contrast with the kitchen sink experiment [18] (Fig. 1, upper left), the uid 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 denes a pressure threshold over which water is eciently evacuated. It is adjusted vertically in order to choose the radius of the hydraulic jump (rjp 20cm). 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 eects of viscosity: the viscous drag is negligible except in the thin layer ahead of the jump. The injection radius is 32cm, 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 velocjp ity v (2gHgrav )1/2 as velocity unit, and the free fall jp timescale t rjp /v as time unit. The dimensionless solution depends on three parameters only [17]: the relative size of the inner boundary r /rjp , the pre-shock velocity v1 /v and the Froude number Fr1 ahead of the shock. The experimental results can be scaled to astrophysical proportions by using the ratio rsh /rjp 106 for distances, and the ratio tsh /tjp for timescales: tsh tjp rsh rjp
3 2

jp rjp gHgrav GMNS

1 2

1.4 102 ,

(3)

where MNS 1.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 [27] and also comparable to the vortical-acoustic cycle seen in a shocked isothermal gas [20]. Symmetry breaking. Indeed, as the water ux 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 3s and a growth rate of 0.2s1 in the experiment for v1 /v 1, rjp /r = 5 and Fr1 4 would correspond to an oscillation period of 42ms and a growth time of 70ms 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 ow contains no angular momentum. This nonlinear behavior, rst 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 uid 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 owing water. The hydraulic jump and the horizontal bar systematically rotate in opposite directions (Fig. 2 and movies M4, M5 [21]). Eect of the bottom drag for a viscous uid. The Reynolds number in the experiment can be dened using the viscosity of water, the uid velocity v and depth

FIG. 2. 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 ow 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]).

FIG. 4. The stability of the ow depends on the strength of the hydraulic jump measured by the pre-jump Froude number Fr1 , and the pre-jump velocity v1 measured in units of the free-fall velocity v 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 (tjp )1 , for a jump radius rjp /r = 4.9 (tjp = 0.35s). The symbols + indicate the set of parameters used in experiments showing a robust m = 1 instability, the lled square corresponds to the parameters of Fig. 3. Dashed lines of constant ow rate (0.5, 1.5 and 2.5 L/s) from bottom to top correspond to a Reynolds number 408, 1224 and 2040 respectively, measured at the jump radius.

H. It is smallest in the outer part of the experiment and inversely scales with radius: Re Hv/ = Q/2r where Q is the injected ow rate. According to Fig. 4, the Reynolds number at the shock radius is in the range 5002000 within the parameters of the experiment, which supports the laminar approximation for the viscous drag [23], proportional to v/H 2 . 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 ow (Fr> 1) by measuring the radial prole of the potential surface and the free surface of the uid between R = 6cm and R = 31cm for an injection slit size of 1.3mm and a ow rate of 1.15L/s. The corresponding velocity prole was best tted with an eective viscosity coecient = 0.03 cm2 /s, such that the viscous drag is written as an additional term v/H 2 in Eq. (2). The perturbative analysis [17] of the 2D shallow water model assumes that this parameter is a constant.
FIG. 3. 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 ow rate is 1L/s (5%) and the injection slit size is 0.74mm (0.02mm).

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 ow parameters were chosen according to their values at the shock (velocity and Froude number) rather than at injection. Measuring the pre-jump velocity in units of the inviscid free fall velocity in Fig. 4, the contour lines showing the growth rate of the instability for = 0.03

4 cm2 /s are well approximated by those obtained for an inviscid ow. 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 are well known in the range 2.5 < Fr < 4.5, while steadier jumps are expected in the range 4.5 < Fr < 9 [24]. This latter is more typical of pre-shock Mach numbers relevant to corecollapse studies. The experiments showed possible interferences between the intermittent character of the jump and the global instability for Fr < 3. Fig. 4 shows the set of most robust experiments performed for Fr > 3.3. The transition to turbulence, expected around Re= 2000, would require a dierent prescription for the viscous drag (proportional to v 2 /H). Viscosity in the horizontal plane and surface tension have been neglected. The vertical structure of the ow 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 coefcients [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 uid mechanics. The complementarity between the experimental and numerical approaches is benecial because their inherent limitations are dierent. The experiment together with its idealized modeling demonstrated the robustness of the SASI instability in a converging ow, even in a uid without entropy gradients. This emphasizes implicitly the importance of vorticity perturbations, like in the vortical-acoustic cycle in an isothermal ow [20]. This new tool will be used to characterize nonlinear eects 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 inuence 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. The design and building of the experiment beneted from the help of P. Mulet, E. Grgoire, L. Rodriguez, F. e Boissonnet and V. Foglizzo. This work is part of ANR funded projects Vortexplosion ANR-06-JCJC-0119 and SN2NS ANR-10-BLAN-0503. J.G. acknowledges support from STFC.

[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 e-prints (2011), arXiv:1112.1427 [MNRAS, in press]. [8] A. Marek and H.-T. Janka, Astrophys. J. 694, 664 (2009). [9] L. Scheck, H.-T. J. T. Plewa, K. Kifonidis, and E. Mller, u Phys. Rev. Lett. 92, 011103 (2004). [10] L. Scheck, K. Kifonidis, H.-T. Janka, and E. Mller, u 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. Mller, Asu trophys. J. 725, L106 (2010). [13] J. Blondin and A. Mezzacappa, Nature 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. uid mechanics, (Pergamon Press, 1987) Chap. 108, 2nd ed. [17] (Supplementary Material), a description of the set of equations dening the 2D model and the numerical technique used for the simulations are available at http://link.aps.org/supplemental/ (also http://irfu.cea.fr/Projets/SN2NS/outreach.html). [18] E. 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] (Supplementary Material), the movies of both the experiment and its numerical modeling are available at http://link.aps.org/supplemental/ (also http://irfu.cea.fr/Projets/SN2NS/outreach.html). [22] W. Iwakami, K. Kotake, N. Ohnishi, S. Yamada, and K. Sawada, Astrophys. J. 678, 1207 (2008). [23] L. Moody, Trans. Am. Soc. Mech. Eng. 66, 671 (1944). [24] V. Chow, Open Channel Hydraulics (McGraw-Hill International, New York, 1973). [25] H. Chanson, The hydraulics of open channel ow: an introduction, (Elsevier, 2004) Chap. 3.2.2, pp. 2729. [26] J. Guilet, J. Sato, and T. Foglizzo, Astrophys. J. 713,

5
1350 (2010). [27] T. Yamasaki and T. Foglizzo, Astrophys. J. 679, 607 (2008).