Response of nucleardissociating shocks to vorticity perturbations
Abstract
In the context of corecollapse supernova explosions (CCSNe), the interaction of standing accretion shocks with upstream vorticity perturbations is investigated by linear theory analysis. The endothermic effect associated to the nuclear dissociation, which takes place right behind the shock wave, affects the amplitude of the perturbations amplified/generated across the front. For upstream disturbances whose characteristic size is much larger than the postshock dissociationlayer thickness, the effect of nuclear dissociation can be reduced to that of considering the global endothermic effect that scales with the inflow energy flux. The present study focuses on perturbation fields that are not isotropic, which mimic the perturbations in collapsing convective shells of massive stars. The linear interaction of the shock with bidimensional monofrequency vorticity perturbations is theoretically addressed, with the limit of highlystretched vortices being analyzed in detail. The exact spatial distribution of the rotational and acoustic perturbations generated in the postshock flow are provided along with the transient evolution of the shock front. It is found that nuclear dissociation contributes to stabilize the shock oscillations, but increases the amplitude of the density perturbations downstream. An extension of this work that addresses the interaction with tridimensional isotropic turbulent flows can be found in reference Huete, C., et al. 2018, MNRAS, 475, 3305â3323, which analyzes the effect of the postshock flow on the critical conditions that ultimately trigger explosion.

August 2018
Keywords: Shock Wave, Supernovae, Turbulence
1 Introduction
A distinguishing feature of massive stars is that the nuclear fusion occurring in their cores continues even after the hydrogen fuel is exhausted. The high temperatures induced by the strong gravity allow heavier elements such as helium and carbon to fuse sequentially. The energy released due to the mass defect between products and reactants keep the star stable, extending the lifetime of the star. Nonetheless, the characteristic time of these advanced fusion stages decreases very rapidly with the nuclear mass, with heavier nuclei burning on a timescale orders of magnitude shorter than the hydrogen sequence, which counts in million of years [1]. In stars with initial masses , the sequential nuclear fusion lasts until the formation of iron nuclei, a point beyond which nuclear fusion is no longer exothermic. As a result, the iron nuclei accumulates in the center, forming a core supported by the pressure of degenerate electrons. When the core reaches its maximum mass of , pressure begins to decay and the hydrodynamical stability breaks down, triggering a collapse of the iron core (e.g., [2] for a recent review).
The core collapse accelerates until the central density becomes as high as nuclear density (), a point where nuclear matter stiffens. This abruptly halts the collapse of the inner iron core, leading to the formation of a shock wave at the boundary of the inner core. The shock has to expel the stellar envelope and thus power corecollapse supernova (CCSN) explosion, leaving behind a stable neutron star (NS). The propagation of the shock, however, does not progress smoothly. The inherent pressure and temperature rise across the shock produces heavynuclei breaking as it propagates outwards, with associated energy consumption. In addition, the hot material behind the shock cools rapidly due to copious neutrino emission. As a result, the shock quickly loses its energy and turns in a stalled accretion shock within milliseconds after formation (see the sketch in Figure 1). Despite the decades of effort, the details of how to revive the shock and power CCSN explosion remain unclear (e.g., [3, 4] for recent reviews).
Along with the nontrivial, yet conventional, gasdynamics effects, CCSNe is a very rich problem that comprises many different phenomena. The newlyborn NS cools and contracts, releasing potential binding energy as neutrino radiation. The shockcompressed matter is sufficiently opaque to absorb a small fraction of these neutrinos (e.g., [5, 6]). The deposited neutrino energy plays a key role in powering the explosion^{1}^{1}1Corecollapse supernovae with explosion energies as high as have also been observed. Also known as hypernovae, such explosions are relatively rare and are believed to be powered by the rotational kinetic energy of rapidly rotating protoneutron stars (e.g., [7, 8]).. Neutrino heating leads to negative gradient of entropy, driving vigorous neutrinodriven turbulent convection in the postshock region ([9, 10] for recent reviews), exerting additional pressure behind the shock [11]. In addition, due to the standing accretion shock instability (SASI) [12, 13, 14], the shock undergoes largescale nonradial oscillations [15, 16], which expands the size of the region subject to net neutrino heating, resulting in higher heating efficiencies.
Recently, it was shown that the convective instabilities that develop in the innermost nuclearburning shells of massive stars can have an important impact on the explosion [17, 18, 19, 20, 21, 23]. The convective motion in oxygen and silicon shells may persist even during iron core collapse. As the core collapses, these shells contract and, due to the conservation of angular momentum, the velocities amplify as , resulting in an increase by a factor of several (e.g., [25, 28]). Upon reaching the shock, these perturbations interact with the supernova shock and generate additional turbulence in the postshock region. This augments the turbulent pressure behind the shock (e.g., [20]), resulting in more favorable conditions for producing explosion.
The details of how these perturbations interact with the shock affects the properties of the resulting supernova explosion. Using the linear interaction analysis (LIA), Abdikamalov et al. [25, 26] studied the effect of acoustic, entropy, and vorticity perturbations, which are the three components of a generic weak hydrodynamic turbulent flow [27]. Huete et al. [28] improved their models by taking into account the perturbation of nuclear dissociation energy at the shock. They employed longtime asymptotic expressions to compute the turbulent amplification ratios across the shock front for incoming vorticity waves. The impact of these modes on the explosion condition can be assessed using the notion of the critical (i.e., the minimum) neutrino luminosity necessary for driving the explosion [20]. The effect of the entropicbuoyant turbulent perturbations generated by incident vorticity waves was found to reduce the critical luminosity by 17–24 per cent, which approximately agrees with the results of threedimensional simulations of CCSNe [29]. The present study is an extension of [28] and it focuses on the linear interaction of the shock with bidimensional singlemode vorticity perturbations. Employing the mathematical formalism used in describing RichtmyerMeshkovtype flows [30, 31, 32] and perturbed nonreactive and reactive shocks [33, 34, 35, 36, 37, 38], the exact spatial distribution of the rotational and acoustic perturbations generated in the postshock flow are provided along with the transient evolution of the shock front towards the permanent oscillatory mode, which, akin to nonideal gases [39, 40], may change the character of the decay when nuclear dissociation is sufficiently high. The effects of the nuclear dissociation energy, the shock strength, and the characteristic frequency are analyzed in the perturbed flow. Finally, due to accelerated pace of stellar collapse, the inner parts of the star collapses faster than the outer parts (e,g, [41]). As a result, the convective vortices undergo substantial stretching during this phase. Such highly elongated vortices are also considered in this work.
The paper is structured as follows: the problem formulation is shown in Section 2, where the baseflow equations are presented and the linearperturbation version are provided. The resulting Euler equations are integrated in Section 3, where the transient and the longtime response of the shock front is computed. The complete distribution of the perturbations in the postshock flow is also shown. The transient evolution of the shock front in the limit of highly elongated vortices is addressed in Section 4. A summary of the results is provided in Section 5.
2 Problem Formulation
2.1 Baseflow equations
Consider a shock wave located at radius and assume that the shock thickness is much smaller than the shock radius (). In this thinshock limit, one can relate hydrodynamic quantities in the postshock region with those in the preshock region using the conservation equations for the mass, momentum and energy:
(1a)  
(1b)  
(1c) 
Here, the flow ahead of the shock () is denoted with subscript 1, while the flow behind () is marked with subscript 2. The variable is the bulk velocity measured in a reference frame at rest with respect to the center of the star, while variables , and represent the density, pressure, and internal energy. We model stellar matter as a perfect gas with the polytropic index on both sides of the shock. Since we consider scenario of infalling flow and expanding shock, it is natural to define velocities positive, , , for the gas moving inwards and for the shock expanding outwards.
In the presence of nuclear dissociation in a thin layer behind the shock front, the change of the internal energy across the shock is
(1b) 
where is the specific energy employed in the nuclear dissociation process. For stalled shock in CCSNe, can be parametrized as , where freefall speed and is a dimensionless parameter [42, 43]. In this scenario, scales as , which results in typically ranging between and [28]. For flows with vanishing Bernoulli parameter above the shock, one can express in terms of the preshock Mach number (see [28] for the details of the derivation),
(1c) 
where is the sound speed in the preshock region and is the preshock speed in the reference frame at rest with respect to the shock.
The fluid properties behind the shock can be conveniently expressed as functions of and preshock Mach number :
(1d) 
(1e) 
(1f) 
where is the mean Mach number in the postshock region and is the speed of the postshock flow in the reference frame at rest with respect to the shock. The function
(1g) 
contains the effect of nuclear dissociation. In the limit of vanishing nuclear dissociation, , expressions (1d)(1f) reduce to the classical RankineHugoniot relations.
The progress of the nuclear dissociation taking place behind the shock can be quantified in terms of variable , such that refers to values right behind the adiabatic shock and refers to the corresponding properties once nuclear dissociation has been accomplished. Normalized with the flow properties right behind the shock, which are directly obtained from (1d)(1g) with , the inner properties are computed as a function of in Figure 2. The effect of nuclear dissociation is to increase the density ratio with the corresponding velocity decrease, and also to slightly increase the downstream pressure. The rate of the heavynuclei breaking would provide the characteristic reaction time, whose combination with the characteristic velocity can be used to scale dissociation layer thickness .
2.2 Linear perturbation analysis
The upstream flow is seldom uniform. The infalling matter is immersed in a strong gravity field that triggers convective instabilities. Vortex cells are then formed upstream and they perturb the shock front. The postshock flow is correspondingly perturbed. In order to study how the shock wave reacts to perturbations, a canonical case is selected to be studied: the interaction of harmonic vorticity perturbations with the shock wave in the planar fastreaction limit . Such an interaction is sketched in Figure 3, where the disturbed shock induces pressure, density and velocity changes in the flow downstream.
The amplitude of the velocity perturbations is assumed to be much smaller than the background flow properties. The small dimensionless amplitude factor , is used to scale the preshock and postshock perturbation variables. The incident shear wave in the frame comoving with the infalling fluid particles is expressed as a divergencefree velocity field
(1ha)  
(1hb) 
where and are the orderunity streamwise and transverse velocity perturbation components. The angle brackets represent the timeaveraged value of the flow variables. The dimensionless vorticity function, associated to the rotational velocity perturbation (1ha)(1hb), is
(1hi) 
where is the perturbation wavenumber in the preshock region, which can be expressed in terms of the wavelengths and sketched in Figure 3, or in terms of the incident shear angle . For the analysis of the postshock flow, it is most natural to use a reference frame comoving with the postshock flow. For this reason, hereafter, the dimensionless coordinates and and the dimensionless time are used to describe the solution in the postshock region.
As a result of the interaction, the shock ripples and the fluid downstream is correspondingly altered with acoustic and entropicvorticity waves. The former travels at the speed of sound downstream and the latter moves with the fluid particles. In the postshock region, the dimensionless pressure, density and velocity perturbations are defined as orderunity functions
(1hj) 
The adiabatic Euler equations governing the postshock flow are written in terms of these variables. Since and are proportional to and , the conservation equations for mass, momentum, momentum, and energy
(1hk) 
can be combined to yield an equation for pressure perturbation :
(1hl) 
This is a periodicallysymmetric twodimensional wave equation, the solution of which yields the perturbation field in the postshock region.
The problem thus reduces to solving the wave equation (1hl) as a function of and for . The spatial domain is bounded by the leading (first) reflected acoustic wave propagating backwards, , and the shock wave front traveling upwards, . One of the boundary conditions stems from the assumption of isolated shock, according to which no reflected sonic waves can reach the shock from behind. This is consistent with the largeradius approximation, , and the linear perturbation analysis that neglects the secondorder interaction of these sonic waves with the entropic or rotational perturbations downstream. The upstream flow distorting the shock front is what ultimately determines the amplitude of the shock oscillations and the postshock inhomogeneities.
The boundary conditions at the shock are readily obtained from the linearized RankineHugoniot relations (1a)(1c) along with the conservation of momentum transverse to the shock front, namely
(1hma)  
(1hmb)  
(1hmc)  
(1hmd) 
where is the time derivative of the dimensionless shock front deformation. As sketched in Figure 3, the transverse velocity induced behind the shock, included in (1hmd), generates a counterpressure effect to balance the shock ripple. The characteristic shock reaction time would determine the shock evolution, which in turn is affected by the energy absorbed in the nuclear dissociation process. This effect is quantified in the energy equation (1hmc) via functions
(1hmn) 
and
(1hmo) 
that allow us to differentiate adiabatic shock waves from reacting shock waves such as shocks with nuclear dissociation or detonations. Functions (1hmn) and (1hmo) are equivalent to those provided in [28], with the latter being here corrected by a factor .
From (1hma)(1hmd), one can obtain one of the two relations for the shock boundary condition involving and , while the other can be obtained using the material derivative of the longitudinal velocity perturbation behind the shock, namely
(1hmpa)  
(1hmpb)  
where
(1hmpq) 
are the factors accompanying the pressure perturbation.
The initial condition for the shock perturbations can be obtained from the requirement that the initial shock is planar, i.e., . Consequently, the initial pressure and streamwise velocity perturbations obey condition , thus yielding , for the initial pressure perturbation right behind the shock.
3 Results
3.1 Temporal evolution of the shock front
In order to study the transient response of the shock to upstream perturbation, the transformation
(1hmpr) 
is employed [35, 38]. The initial condition, , corresponds to the moment when the initial unperturbed shock first meets the incident vorticity perturbations. The const condition represents a planar surface moving in the postshock gas along the direction, from the weak discontinuity at to the reacting shock front at . In terms of these variables, equation (1hl) for sound waves reads
(1hmps) 
The boundary conditions at the shock front reduces to
(1hmpt) 
and
(1hmpu) 
where
(1hmpv) 
is the characteristic shock oscillation frequency induced by the incident shear wave.
The Laplace transform is conveniently employed to reduce the above system of partial differential equations to an algebraic system. That is, the integral
(1hmpw) 
applied to the functions defining the shock boundary conditions yields an algebraic system of function . The Laplace transform of can be computed using the isolated boundary condition, namely . From this, one can obtain the Laplace transform of the pressure perturbation at the shock:
(1hmpx) 
where
(1hmpy) 
is the factor accounting for the periodic excitation amplitude, and
(1hmpz) 
The solution of equation (1hmps) for the pressure field can be expressed as a combination of the Bessel functions [44], as shown in [28] in this particular context. It is however illustrative to use the inverse of the Laplace transform to the function (1hmpx), which yields
(1hmpaa) 
as the temporal evolution of the shock pressure perturbations, with
(1hmpab) 
being the auxiliary function. The corresponding is readily given by the variable change .
Selfinduced stable oscillations are found to depend on the slope of the RankineHugoniot curve [45, 46, 47, 39, 40, 48] and the corresponding spontaneous acoustic radiation is determined by the condition . A parametrical study reveals that it would occur only for , i.e., when the net positive energy release increases with the shock intensity. In this case, as the shock will oscillate only with the excitement frequency coming from upstream perturbations, , thereby yielding an asymptotic response qualitatively similar to the one found for adiabatic shock waves [35]. Nonetheless, the endothermic contribution may have a qualitative impact on the transient evolution towards the longtime dynamics. Although transient evolution always decays in time with for , when the function
(1hmpac) 
is lower than unity (or in [39]), the initial degree of damping is significantly modified. That is, for , corresponding highly endothermic shocks, the oscillations associated to the transient response are effectively shorten, while the contrary occurs for , the latter case corresponding to the regular shock family of solutions. The delimiting curve is computed in Figure 4 for as a function of the shock strength and the endothermic parameter .
Irrespective of the transient behavior, the longtime reaction of the shock pressure to monofrequency perturbations is
(1hmpad) 
where
(1hmpae) 
for , and
(1hmpaf) 
for . As in previous studies of the interaction of shocks with vorticity perturbations [49, 50, 51, 35, 28], in the long wavelength regime (), the sonic disturbances immediately behind the shock consists of two orthogonal contributions , and . In this regime, the sonic waves decay exponentially as they move away from the shock. In contrast to this, in the short wavelength regime (), the solution is represented by constantamplitude sonic waves. The critical value corresponds to the case when stable acoustic waves travel parallel to the shock surface in the shock reference frame.
Likewise, the integration of (1hmpb) with respect to yields the temporal evolution of the shock oscillation amplitude:
(1hmpag)  
with the variable . The associated longtime function of the shock oscillations is
(1hmpah) 
where the coefficients are obtained from the pressure fluctuations at the shock:
(1hmpai) 
A direct comparison of the longtime response, provided by (1hmpah), and the exact temporal evolution, given by (1hmpag), is computed in Figure 5 as a function of for a shock strength with two different dissociation sensitivities (a) and (b). The figure is qualitatively similar to figure 5 in [28], where different input parameters had been selected in this occasion. The upstream shearwave angle is , which yields dimensionless oscillation frequencies greater than unity in both cases: and for the adiabatic and endothermic case, respectively. The transient solution is found to achieve the asymptotic regime in a relatively short period of time, with the panel on the left doing it seemingly faster, in agreement with Figure 4. As expected, the amplitude of the longtime oscillations is found to be smaller in the endothermic case.
Figure 6 shows the amplitude of the oscillations as a function of the shock strength , shearwave angle , and for three different values of the dissociation degree (a), (b), and (c). In agreement with Figure 5, the amplitude is found to be generally smaller for endothermic shocks. The limit is also computed in Figure 6 as a function of and for different endothermic intensities. The zones on the left of these dashed curves correspond to pressure radiating conditions (highfrequency regime), and the zones on the right refer to nonradiating conditions (lowfrequency regime), with the area of the latter being reduced with the increase of nuclear dissociation.
3.2 Postshock flow variables
The downstream flow perturbations are conveniently split decomposed those of acoustic type, which travel at the speed of sound relative to the fluid particles, and those of entropicrotational nature, which move with the fluid particles [27].
The exact temporal evolution of the pressure field downstream is readily obtained through Bessel functions, as derived in [28], or by direct integration of the wave equation provided that . The former strategy has been employed in Figure 7 to compute the pressure field in the domain . Computations reveal two welldistinguished regimes: acoustically radiating and nonradiating conditions. The acoustic radiation condition is then determined by , a condition that depends on the upstream shear wave, since depends on the relative properties of the perturbation field ahead of the shock. Small values of represent the interaction with upstream vortices highly stretched in the streamwise direction , while the opposite is true for . In the latter low modenumber scenario (), the problem reduces to the onedimensional interaction of the shock with radial perturbation waves.
The sonic waves traveling in the postshock region are functions of , where the frequency and the wavenumber are obtained from the shock oscillation frequency and the adiabatic dispersion relation :
(1hmpaj) 
It is readily seen that can be either positive or negative. The case represents acoustic waves traveling in the direction of the postshock flow, while the waves moving in the opposite direction towards the shock have . The solution corresponding to shock oscillation frequency represents the standing acoustic waves that separate the solution traveling to the left from the one traveling to the right in the reference frame comoving with the postshock fluid. At large distances from the shock in the downstream region (far larger than the wavelength of the perturbations), the asymptotic pressure and the isentropic density perturbations are given by
(1hmpak) 
where is the amplitude of the pressure perturbations. The isentropic temperature variations induced by the acoustic radiation are just .
The acoustic contribution of the velocity perturbations are readily obtained through the Euler momentum equations
(1hmpal) 
which can be used to write the longtime response
(1hmpam)  
(1hmpan) 
where
(1hmpao) 
are the associated amplitudes.
In absence of dissipative effects, vorticity disturbances downstream remain frozen to the fluid particles, with the amplitude being determined by the vorticity produced at the shock, namely
(1hmpap) 
where
(1hmpaq) 
accounts for the shockcompression of the vortices, a onedimensional effect, and
(1hmpar) 
indicates the contribution of shock deformations, which is a twodimensional effect.
The linear perturbations in the velocity field satisfy the equation
(1hmpas) 
with the breakdown of irrotationalacoustic and steadyrotational perturbations obeying, separately,
(1hmpat) 
With the acoustic field being given by (1hmpam) and (1hmpan), the spatial distribution of the rotationalvelocity perturbations is now derived to provide the complete velocity field. The solenoidal part is calculated by tracking the vorticity left behind by the oscillating shock front from , yielding
(1hmpaua)  
(1hmpaub) 
with the particular solution being provided by
(1hmpauav)  
and with the auxiliary integration function being defined in (1hmpab).
The function (1hmpauav) is qualitatively similar to that shown in [35] and [38] for adiabatic and reacting shocks waves, respectively. Details of its derivation, omitted here for the sake of conciseness, involve the use of the inverse Laplace transform technique. The asymptotic rotational contribution of the velocity field is also written as a piecewise function of the shock oscillation frequency, with the longitudinal
(1hmpauaw) 
and transverse
(1hmpauax) 
contributions being characterized by their amplitudes
(1hmpauay) 
and
(1hmpauaz) 
for the streamwise and crosswise components in the longwavelength regime, respectively, and
(1hmpauba) 
for the corresponding shortwavelength regime amplitudes.
The density variations due to the entropy waves is obtained from RankineHugoniot relations (1hma)(1hmc) by subtracting the contribution of sonic waves:
(1hmpaubb) 
where is the density perturbation amplitude relative to the pressure at the shock. The temperature variation corresponding to the entropy waves is given by , where the temperature is normalized with the base flow temperature.
Figure 8 shows the spatial distribution of the rotational velocity field and entropic density perturbations from the origin for , and . The upper panels (a)(c) display the functions , , and as a function of , respectively, and the lower panel (d) shows the vector velocity field superposed to the density isocontours in the plane . The exact solution for the rotationalvelocity contribution (solid line) is found to approach the asymptotic solution (reddashed line) in a fairly short distance, while the entropicdensity function exhibits a longer transient period towards the longtime solution.
The amplitude of the asymptotic rotational and acoustic velocity perturbations, as well as the corresponding entropic and acoustic density perturbations, have been conveniently expressed in terms of the shock pressure perturbation amplitude. Likewise, the asymptotic shock oscillation amplitude has been written as a function of in equation (1hmpai). It is then immediate to obtain the value of any perturbation variable with the aid of Figure 6. Further computations, as the longtime amplitudes for the velocity and density perturbations, have been computed in Figures A1A3 of [28]. As found in Figure 6, there exists a peak in the perturbation amplitude near the critical frequency and the effect of endothermicity, along with the associated amplitude change, is to stretch the peak in the frequency domain. This effect occurs for any perturbation variable: pressure, velocity or density. The acoustic contribution is found to provide a negligible contribution in comparison to the rotational or entropic counterpart.
Density perturbations generated by the interaction of upstream asphericities with the shock are found to play a pivotal role in driving postshock turbulence by buoyancy effects, which translates into a reduction of the critical neutrino luminosity necessary for producing explosion [9, 52]. In [28], the effect of vorticity waves with isotropic distribution of orientations was considered. This is perhaps not an accurate representation. The convective motion is characterized by a dominant eddy with a specific size that undergo stretching due to the accelerated collapse, which increases the value of the shear angle .
Anticipating that the entropic contribution dominates the postshock density perturbations, the correlation between the upstream turbulent Mach number and the entropicdensity perturbations is used to evaluate the contribution of the shockgenerated density fluctuations through the factor
(1hmpaubc) 
which is computed in Figure 9 as a function of the upstream Mach number and the shearwave angle . The upper region, corresponding to the highlyelongated vorticity limit (), is found to yield stronger density perturbations than those produced for (), the lower limit. Modification of the endothermic contribution does not change the qualitative picture significantly, but it increases the figures. That is, despite pressure perturbations and shockripple deviations decrease with the endothermicity associated to the nuclear dissociation mechanism, entropicdensity perturbations grow with the factor .
4 Interaction with highlystretched vortices
Natural convective cells are typically formed when temperature gradients are counteraligned with the gravity field, conditions that are met in massive stars. When, in addition, they are advected by a highlyaccelerating radiallyconverging flow, they elongate in the radial direction. In this scenario, the upstream perturbation field is strongly anisotropic, so that the canonical shockturbulence interaction approach is no longer applicable. It is, however, possible to construct a representative problem setup with the formulation previously presented, as depicted in Figure 10.
In the planarshock limit, i.e., when the characteristic length of the upstream perturbation field is much smaller than the shock radius, the formulation of the problem shown in previous section can be particularized to the case or, equivalently, . In this slender limit, the acousticallyinduced shock oscillation period is much shorter than the characteristic residence time of the convective cells crossing out the shock, thereby providing the following upstream modulation
(1hmpaubd) 
as the transverse contribution becomes negligible.
The corresponding boundary conditions at the shock front are appropriately adapted to yield
(1hmpaube) 
and
(1hmpaubf) 
which have been conveniently written as a function of the variables and [34]. After some straightforward manipulation, it is obtained the Laplace Transform of the pressure perturbations, namely
(1hmpaubg) 
provided that , and with the factor in (1hmpy) being now
(1hmpaubh) 
The initially planar shape of the shock front is distorted as a result of the interaction with the upstream modulated velocity field. For , the shock encounters an always positive velocity perturbation aligned with the shock propagation, which pulls the shock front upwards. The opposite would apply for , where an always negative velocity input pushes the front backwards. Then, the tangential velocity perturbation induced by the shock corrugation, with the associated mass flux, tends to restore the shock shape. This two countereffects are responsible of the shock oscillations in the initial stage. Eventually, the shock approaches a steadystate regime when the two opposed effects balance, something that cannot occur when the upstream nonuniform flow is made of vortices with . The corresponding longtime functions take the forms and for the shock pressure and ripple variables, respectively, with the associated amplitudes being
(1hmpaubi) 
The evolution of the shock pressure and shock ripple perturbations, when the front travels through a stripelike pattern of longitudinal perturbations, is computed in Figure 11 for , and for and . The transient behaviour is readily obtained by just taking the limit in (1hmpaa) and (1hmpag), respectively. It is observed that endothermic nuclear dissociation exhibits a stabilizing mechanism, since the longtime amplitudes are smaller and the oscillations decay faster for than for . This is in consonance with Figure 4, where the case and lies on the region , while regular adiabatic shocks with lie on the zone .
In regard to the farfield perturbations, it is immediate to see that the acoustic contribution is negligible, as pressure perturbations decay exponentially with the distance from the shock for . What remains constant, in absence of diffusive effects, is the entropyvorticity perturbation generated across the shock, namely
(1hmpaubj) 
Likewise, by direct inspection of (1hmpaub) and anticipating that tangential velocity generated behind the shock when constant is of acoustic type, and then evanescent, it is found that
(1hmpaubk) 
for the longitudinal and transverse components of the rotational velocity field.
5 Conclusions
When a shock wave encounters a vorticity wave on its way, their interaction results in a deformation of the shock. The latter induces pressure changes that results in the radiation of acoustic waves downstream. Along with the sonic mode, postshock perturbations include vortical and entropic disturbances that are convected by the fluid particles downstream. The distinguished feature of nucleardissociating shocks is the endothermic contribution, which depends on the shock intensity relative to the flow stream. Then, in a likely nonuniform context, perturbations ahead of the shock as the shearpattern considered in this work, the amount of nuclei that are dissociated is affected by the shock perturbation, which in turn affects the energy balance across the shock, and ultimately the amplitude of the perturbations downstream. This effect can be easily studied when considering the interaction of the shock with intermediate vortical scales, i.e., those whose characteristic length is sufficiently small for the shock to be considered a planar front, yet sufficiently large for the shock to be a seen as a discontinuity front.
This problem is particularly relevant in the context of corecollapse supernova explosions, where the shock after bounce swallows the convective structures generated upstream within the supersonic inwardstraveling mass. The shockperturbation interaction modifies and creates additional perturbations downstream, which affects the critical conditions for the supernova explosion. When the upstream flow is assumed to be dominantly isotropic, theory predicts that the injection of nonradial motion and the buoyancydriven convection triggered by entropy waves reduce the critical neutrino luminosity by 1224 per cent, for typical problem parameters [9, 25, 28]. When convective cells are not sufficiently turbulized by the inwards acceleration of the fluid particles, the isotropic assumption may be inaccurate. Motivated by this fact, in this work the emphasis is placed in bidimensional structures made of isodensity vortices. The temporal evolution of the oscillating shock and the longtime asymptotic expressions have been derived analytically as a function of the dominant governing parameters: the shock strength , the nuclear dissociation degree , and the incident shearwave angle . Likewise, the exact and asymptotic spatial distribution of the perturbations in the shocked gas have been derived explicitly with use made of the Laplace Transform technique. The limit of high elongated vortices, corresponding to , is also evaluated in terms of closedform expressions.
The effect of the endothermic nuclear dissociation is found to diminish the amplitude of the shock oscillations and to reduce the acoustic radiation. Contrarily, the entropic density disturbances grow with the factor . In some conditions, for sufficiently endothermic shock, the transient evolution is significantly shorten. The distribution of the downstream entropic and vortical perturbations have been provided analytically. These perturbations serve as initial conditions of the postshock regime, as density variances become buoyant in a temporal scale that exceeds the shock influence. This phenomenon will be studied in more depth in a future work.
Acknowledgements
This work is supported by the Ministry of Science, MEC (ENE201565852C21R) and Fundación Iberdrola España (BINVhBbhOeJQ), Spain (for C. Huete), by MES RK statetargeted program BR05236454, MES RK grant No. 346, NU ORAU grant SST 2015021 and NU grant No. 090118FD5348 (for E. Abdikamalov).
References
References
 [1] Woosley SE, Heger A, and Weaver TA. 2002 The evolution and explosion of massive stars. Rev. Mod. Phys. 74 1015–1071.
 [2] Janka HT. 2012 Explosion mechanisms of corecollapse supernovae. Annu. Rev. Nucl. Part. Sci. 62 407.
 [3] Burrows A. 2013 Colloquium: Perspectives on corecollapse supernova theory. Rev. Modern Phys. 85 245.
 [4] Janka HT, Melson T, and Summa A. 2016 Physics of corecollapse supernovae in three dimensions: a sneak preview. Annu. Rev. Nucl. Part. Sci. 66 341–375.
 [5] Lentz EJ, Mezzacappa A, Messer OB, Hix WR, and Bruenn SW. 2012 Interplay of neutrino opacities in corecollapse supernova simulations. Astrophys. J. 760 94.
 [6] Kotake K, Takiwaki T, Fischer T, Nakamura K, and MartínezPinedo G. 2018 Impact of Neutrino Opacities on Corecollapse Supernova Simulations. Astrophys. J. 853 170.
 [7] Burrows A, Dessart L, Livne E, Ott CD, and Murphy J. 2007 Simulations of magnetically driven supernova and hypernova explosions in the context of rapid rotation. Astrophys. J. 664 416.
 [8] Mösta P, Richers S, Ott CD, Haas R, Piro AL, Boydstun K, Abdikamalov E, Reisswig C, and Schnetter E. 2014 Magnetorotational corecollapse supernovae in three dimensions. Astrophys. J. Lett. 785 L29.
 [9] Müller B. 2016 The Status of MultiDimensional CoreCollapse Supernova Models. Publ. Astron. Soc. Aust 33 E048.
 [10] Radice D, Abdikamalov E, Ott CD. et al. 2018 Turbulence in corecollapse supernovae J. Phys. G 45 053003.
 [11] Murphy J, Dolence JC, and Burrows A. 2013 The dominance of neutrinodriven convection in corecollapse supernovae. Astrophys. J. 771 52.
 [12] Blondin JM, Mezzacappa A, and DeMarino C. 2003 Stability of standing accretion shocks, with an eye toward corecollapse supernovae. Astrophys. J. 584 971–980.
 [13] Foglizzo T, Scheck L, and Janka HT. 2006 Neutrinodriven convection versus advection in corecollapse supernovae. Astrophys. J. 652 1436–1450.
 [14] Foglizzo T, Galletti P, Scheck L, and Janka HT. 2007, Instability of a stalled accretion shock: evidence for the advectiveacoustic cycle. Astrophys. J. 654 1006–1021.
 [15] Fernández R. 2015 Threedimensional simulations of SASIand convectiondominated corecollapse supernovae. Mon. Notices Royal Astron. Soc. 452, 2071–2086.
 [16] Hanke F, Müller B, Wongwathanarat A, Marek A, and Janka HT. 2013 SASI activity in threedimensional neutrinohydrodynamics simulations of supernova cores. Astrophys. J. 770 66.
 [17] Couch SM, and Ott CD. 2013 Revival of the stalled corecollapse supernova shock triggered by precollapse asphericity in the progenitor star. Astrophys. J. Lett. 778 L7.
 [18] Couch, S. M., & Ott, C. D. 2015 The Role of Turbulence in Neutrinodriven Corecollapse Supernova Explosions. Astrophys. J., 799, 5.
 [19] Couch SM, Chatzopoulos E, Arnett WD, and Timmes FX. 2015 The threedimensional evolution to core collapse of a massive star. Astrophys. J. Lett. 808 L21 .
 [20] Müller B, and Janka, HT. 2015 NonRadial Instabilities and Progenitor Asphericities in Core–Collapse Supernovae. Mon. Notices Royal Astron. Soc. 448 2141–2174.
 [21] Müller, B., Viallet, M., Heger, A., & Janka, H.T. 2016, The Last Minutes of Oxygen Shell Burning in a Massive Star. Astrophys. J., 833, 124.
 [22] Collins, C., Müller, B., & Heger, A. 2018, Properties of convective oxygen and silicon burning shells in supernova progenitors. Mon. Notices Royal Astron. Soc., 473, 1695.
 [23] Nagakura, H., Takahashi, K., & Yamamoto, Y. 2018, On the importance of progenitor asymmetry to shock revival in corecollapse supernovae, Mon. Notices Royal Astron. Soc. 483 208–222.
 [24] Takahashi, K., Iwakami, W., Yamamoto, Y., & Yamada, S. 2016, Links between the Shock Instability in Corecollapse Supernovae and Asymmetric Accretions of Envelopes. Astrophys. J., 831, 75.
 [25] Abdikamalov E, Zhaksylykov A, Radice D, and Berdibek S. 2016 Shock–turbulence interaction in core–collapse supernovae. Mon. Notices Royal Astron. Soc. 461 3864–3876.
 [26] Abdikamalov E, Huete C, Nussupbekov A, and Berdibek S. 2018 Turbulence Generation by Shock–Acoustic–Wave Interaction in Core–Collapse Supernovae. Particles 1 7.
 [27] Kovasznay LSG. 1953 Turbulence in supersonic flow. J. Aero. Sci. 20 657–682.
 [28] Huete C, Abdikamalov E, and Radice D. 2018 The impact of vorticity waves on the shock dynamics in corecollapse supernovae. Mon. Notices Royal Astron. Soc. 475 3305–3323.
 [29] Müller B, Melson T, Heger A, and Janka HT. 2017 Supernova simulations from a 3D progenitor model – Impact of perturbations and evolution of explosion properties. Mon. Notices Royal Astron. Soc. 472 491–513.
 [30] Wouchuk JG. 2001 Growth rate of the linear Richtmyer–Meshkov instability when a shock is reflected. Phys. Rev. E 63 056303.
 [31] Wouchuk JG. 2001 Growth rate of the Richtmyer–Meshkov instability when a rarefaction is reflected. Phys. Plasmas 8 2890–2907.
 [32] Campos FC, Wouchuk JG. 2014 Analytical asymptotic velocities in linear Richtmyer–Meshkov–like flows. Phys. Rev. E 90 053007.
 [33] Jackson TL, Kapila AK, and Hussaini MY, 1990 Convection of a pattern of vorticity through a reacting shock wave. Phys. Fluids A 2 1260.
 [34] Velikovich AL, Wouchuk JG, Huete Ruiz de Lira C, Metzler N, Zalesak S, and Schmitt AJ. 2007 Shock front distortion and Richtmyer–Meshkovtype growth caused by a small preshock nonuniformity. Phys. Plasmas 14 072706.
 [35] Wouchuk JG, Huete Ruiz de Lira C, and Velikovich AL. 2009 Analytical linear theory for the interaction of a planar shock wave with an isotropic turbulent vorticity field. Phys. Rev. E 79 066315.
 [36] Huete Ruiz de Lira C. 2010 Turbulence generation by a shock wave interacting with a random density inhomogeneity field. Phys. Scripta T142 014022.
 [37] Huete C, Sánchez AL, and Williams FA. 2013 Theory of interactions of thin strong detonations with turbulent gases. Phys. Fluids 25 076105.
 [38] Huete C, Jin T, MartÃnezRuiz D, and Luo K. 2017 Interaction of a planar reacting shock wave with an isotropic turbulent vorticity field. Phys. Rev. E 96 053104.
 [39] Bates JW. 2004 Initialvalueproblem solution for isolated rippled shock fronts in arbitrary fluid media. Phys. Rev. E 69 056313.
 [40] Bates JW. 2007 Instability of isolated planar shock waves. Phys. Fluids 19 094102.
 [41] Bruenn, S. W. 1985 Stellar core collapse  Numerical model and infall epoch. Astrophys. J. S. 58 771
 [42] Fernández R, and Thompson C. 2009 Stability of a spherical accretion shock with nuclear dissociation. Astrophys. J. 697 1827–1841.
 [43] Fernández R, and Thompson C. 2009 Dynamics of a spherical accretion shock with neutrino heating and alphaparticle recombination. Astrophys. J. 703 1464–1485.
 [44] Zaidel PM. 1960 Shock wave from a slightly curved piston. J. Appl. Math. Mech 24 316–327.
 [45] D’yakov, SP. 1954 The stability of shockwaves: Investigation of the problem of stability of shock waves in arbitrary media. Zh. Eksp. Teor. Fiz. 27 288–295.
 [46] Kontorovich VM. 1957 On the stability of shock waves. Zh. Eksp. Teor. Fiz. 33 1525–1526.
 [47] Wouchuk JG, and López Cabada J. 2004 Spontaneous acoustic emission of a corrugated shock wave in the presence of a reflecting surface. Phys. Rev. E 70 046303.
 [48] Clavin P, and Williams FA. 2012 Analytical studies of the dynamics of gaseous detonations. Phil. Trans. R. Soc. A. 13 597624.
 [49] Ribner HS. 1954 Convection of a pattern of vorticity through a shock wave. NACA Report 1164.
 [50] Lee S, Lele LK, and Moin P. 1997 Interaction of isotropic turbulence with shock waves: effect of shock strength. J. Fluid Mech. 340 225–247.
 [51] Mahesh K, Lele SK, and Moin P. 1997 The influence of entropy fluctuations on the interaction of turbulence with a shock wave. J. Fluid Mech. 334 353–379.
 [52] Mabanta, Q. A., & Murphy, J. W. 2018 How Turbulence Enables Corecollapse Supernova Explosions, Astrophys. J., 856, 22