Energy Conservation in the Thin Layer Approximation: I. The Spherical Classic Case for Supernovae Remnants ()
1. Introduction
The thin layer approximation assumes that the mass ejected in the explosion of a supernova (SN) resides in a thin layer. This approximation is usually applied in the late stage of the explosion in order to explain the supernova remnant (SNR), see [1] [2] [3]. The physical quantity which is conserved in the previous approaches is the momentum, equal to the swept mass multiplied by the velocity at a given radius of expansion
equated to these quantities at a radius r. Some natural questions therefore arise:
· Can we model the expansion of an SNR when the energy is conserved rather than the momentum?
· Can we model the energy conservation when the density of the interstellar medium (ISM) decreases with the distance from the point of the explosion?
In order to answer the above questions, Section 2 reviews the standard laws of conservation, Section 3 introduces the conservation of energy and Section 4 applies the derived equations of motion to 4 SNRs.
2. Laws of Conservation
We summarise four laws of conservation useful to model some astrophysical phenomena in which the temperature and the pressure are absent. The first law is the conservation of momentum in spherical coordinates in the framework of the thin layer approximation. The Newton’s second law for an expanding sphere in the framework of the thin shell approximation along a solid angle
is
(1)
where r is the advancing radius,
is the density assumed to be constant, v the velocity and P the internal pressure, see formula (10.27) in [4]. Let us assume
(cold model) and the above equation in two different points of expansion becomes
(2)
where
and
are the swept masses at
and r, while
and v are the velocities of the thin layer at
and r. This first law has been widely used to model the SNRs, see [5] - [10]. This conservation law can be expressed as
a differential equation of the first order by inserting
:
(3)
In the case where the ISM has constant density, the analytical solution for the trajectory is
(4)
and the velocity is
(5)
where
and
are the position and the velocity when
. The second law is the conservation of energy which will be introduced in details in the next section. An example is given by the energy conserving phase in the interstellar bubbles, see [4]. The third law of conservation is given by the conservation of momentum flux which is the rate of transfer of momentum through a unit area
(6)
where
is the density at position x,
is the area at position x and
is the velocity at position x, see Formula A27 in [11]. This law is useful to model the radiogalaxies where there is a continuous flow of matter from the central region to the periphery, see [12]. The fourth law of conservation is given by the conservation of energy flux which is the rate of transfer of energy through a unit area
(7)
where
is the density at position x,
is the area at position x and
is the velocity at position x, see Formula A28 in [11]. This law is useful to model the astrophysical jets, see [13].
3. Energy Conservation
The conservation of kinetic energy in spherical coordinates within the framework of the thin layer approximation when the thermal effects are negligible is
(8)
where
and
are the swept masses at
and r, while
and v are the velocities of the thin layer at
and r. The above conservation law, when written as a differential equation, is
(9)
The velocity as a function of the momentary radius is
(10)
In the following, the case of constant density as well as 7 profiles of decreasing density will be considered.
3.1. Medium with Constant Density
When the ISM is considered to have constant density, the analytical solution for the trajectory when the energy is conserved is
(11)
which has the asymptotic behaviour
,
(12)
The velocity as function of the radius is
(13)
and the velocity as a function of time is
(14)
where
and
are the position and the velocity when
.
3.2. Constant Density and Back Reaction
The radiative losses per unit length are assumed to be proportional to the flux of momentum
(15)
where
is a constant and
is density in the thin advancing layer which is
. Inserting in the above equation the velocity to first order as given by Equation (13) the radiative losses,
, are
(16)
The sum of the radiative losses between
and r is given by the following integral, L,
(17)
The conservation of energy in presence of the back reaction due to the radiative losses is
(18)
The analytical solution for the velocity to second order,
, is
(19)
The inclusion of back reaction allows the evaluation of the SRS’s maximum length
, which can be derived imposing to zero the above velocity.
(20)
3.3. Medium with an Hyperbolic Profile of Density
We assume that the medium around the SN scales with the piecewise dependence
(21)
where
is the density at
and
is the radius after which the density starts to decrease. The mass swept,
, in the interval
is
The total mass swept,
, in the interval
is
The application of energy conservation gives the velocity as a function of the radius:
(22)
Separation of variables followed by integration gives
(23)
In this equation it is not possible to extract the radius as a function of time, and therefore a numerical procedure is adopted in order to derive the trajectory.
3.4. Medium with an Inverse Square Profile for the Density
We now assume that the medium around the SN scales with the piecewise dependence (which avoids a pole at
)
(24)
where
is the density at
and
is the radius after which the density starts to decrease.
The total mass swept,
, in the interval
is
Applying the conservation of energy, the velocity as a function of the radius is
(25)
The trajectory, i.e. the radius as a function of time, is
(26)
which has the asymptotic behavior,
,
(27)
The velocity as a function of time is
(28)
3.5. Medium with a Power Law Profile for the Density
We now assume that the medium around the SN scales as
(29)
where
is the density at
,
is the radius after which the density starts to decrease and
.
The total mass swept,
, in the interval
is
The application of energy conservation gives the differential equation
(30)
The velocity as a function of the radius is
(31)
There is no analytical solution for the trajectory, and therefore we have implemented a numerical procedure. The first approximation for the trajectory is obtained by a series solution of Equation (30) to fourth order,
(32)
The second approximation for the trajectory is found by first deriving an asymptotic expansion of Equation (31), namely
(33)
Then, the asymptotic approximate trajectory turns out to be
(34)
3.6. Medium with an Exponential Profile for the Density
We assume that the medium around the SN scales with the piecewise dependence
(35)
where
is the density at
and
is the radius after which the density starts to decrease. The total mass swept,
, in the interval
is
The application of energy conservation gives the differential equation
(36)
The velocity as a function of the radius is
(37)
where
(38)
and
(39)
There is no analytical solution for the trajectory, and therefore we present a series solution of Equation (36) to fourth order:
(40)
3.7. Medium with a Gaussian Profile for the Density
We assume that the medium around the SN scales with the piecewise dependence
(41)
where
is the density at
and
is the radius after which the density starts to decrease. The total mass swept,
, in the interval
is
(42)
where
is the error function, defined by
(43)
See [14].
The differential equation when the energy is conserved is
(44)
In the absence of an analytical solution for this differential equation, we present an approximation using the fourth order Taylor series:
(45)
3.8. Autogravitating Medium
We assume that the medium around the SN scales with the piecewise dependence
(46)
where
is the density at
,
is the radius after which the density starts to decrease and sech is the hyperbolic secant ( [15] [16] [17] [18]).
The total mass swept,
, in the interval
is
(47)
where the polylog operator is defined by
(48)
and
is a Dirichlet series. The differential equation when the energy is conserved is
(49)
where
(50)
The velocity as a function of the radius is
(51)
where
(52)
In the absence of an analytical solution for this differential equation, we present the approximation arising from the fourth order Taylor series:
(53)
3.9. Medium with an NFW Profile
We assume that the medium around the SN scales with the Navarro-Frenk-White (NFW) distribution as follows:
(54)
where
is the density at
, and
is the radius after which the density starts to decrease, see [19]. The total mass swept,
, in the interval [0, r] is
(55)
The differential equation when the energy is conserved for an NFW profile is
(56)
where
(57)
The velocity as a function of the radius is
(58)
where
(59)
This differential equation does not have an analytical solution, so we present the approximation arising from the fourth order Taylor series:
(60)
4. Astrophysical Applications
We now test the reliability of the numerical and approximate solutions on four SNRs: Tycho, see [20], Cas A, see [21], Cygnus loop, see [22], and SN 1006, see [23]. The three astronomically measurable parameters are the time since the explosion in years, t, the actual observed radius in pc, r, and the present velocity of expansion in km∙s−1, see Table 1. The astrophysical units are pc for length and yr for time. With these units, the initial velocity is
. In all the models here considered, the initial velocity,
, is constant in the time interval
.
The goodness of the model is evaluated through the percentage error
of the radius, which is
(61)
where
is the radius of the SNR as given by the astronomical observations and
is the radius suggested by the model. In an analogous way, we can define the percentage error of the velocity. Another useful astrophysical variable is the predicted decrease in the theoretical velocity in 10 years,
.
4.1. Constant Density
The numerical results for the medium with constant density are presented in Table 2.
4.2. Power Law Densities
The results for a medium with an hyperbolic density are presented in Table 3, those for the medium with an inverse square profile of density are presented in Table 4, and those for the medium with an inverse power law profile of density are presented in Table 5.
Table 1. Observed astronomical parameters of the SNRs.
Table 2. Theoretical parameters of the SNRs for the equation of motion in the case of conservation of energy with constant density, see Section 3.1.
Table 3. Theoretical parameters of the SNRs for the equation of motion in the case of conservation of energy with an hyperbolic profile of density, see Section 3.3.
Table 4. Theoretical parameters of the SNRs for the equation of motion in the case of conservation of energy with an inverse square profile of density, see Section 3.4.
Table 5. Theoretical parameters of the SNRs for the equation of motion in the case of conservation of energy with a power law profile of density when
, see Section 3.5.
In the case of a density which decreases with a power law profile we have already pointed out the absence of an analytical solution. As a consequence, Figure 1 presents the asymptotic approximate trajectory as given by (34) for Tycho in the full range of time. Figure 2 presents the Taylor approximation of the trajectory as given by (32) in the restricted range of time.
Figure 1. Numerical solution (full red line) and asymptotic approximate solution (blue dashed line) for the inverse power law with Parameters as in Table 5 for Tycho.
Figure 2. Numerical solution (full red line) and Taylor approximation (blue dashed line) for the inverse power law with. Parameters as in Table 5 for Tycho.
4.3. Presence of an Exponential
The astrophysical parameters for an exponential profile of density are presented in Table 6 and the fit of the trajectory with a Taylor expansion, see Equation (40), is presented in Figure 3.
The astrophysical parameters for a Gaussian profile of density are presented in Table 7 and the fit of the trajectory with a Taylor expansion, see Equation (45), is presented in Figure 4.
4.4. Autogravitating Medium
The astrophysical parameters for an autogravitating medium are presented in Table 8 and the fit of the trajectory with a Taylor expansion, see Equation (53), is presented in Figure 5.
4.5. NFW Profile
The astrophysical parameters for an NFW profile of density are presented in Table 9 and the fit of the trajectory with a Taylor expansion, see Equation (60), is presented in Figure 6.
Table 6. Theoretical parameters of the SNRs for the equation of motion in the case of conservation of energy with an exponential profile of density, see Section 3.6.
Table 7. Theoretical parameters of the SNRs for the equation of motion in the case of conservation of energy with a Gaussian profile of density, see Section 3.7.
Table 8. Theoretical parameters of the SNRs for the equation of motion in the case of conservation of energy with an autogravitating profile of density, see Section 3.8.
Table 9. Theoretical parameters of the SNRs for the equation of motion in the case of conservation of energy with an NFW profile of density, see Section 3.9.
Figure 3. Numerical solution (full red line) and Taylor approximation (blue dashed line) for the exponential profile. Parameters as in Table 6 for Tycho.
Figure 4. Numerical solution (full red line) and Taylor approximation (blue dashed line) for the Gaussian profile. Parameters as in Table 7 for Tycho.
Figure 5. Numerical solution (full red line) and Taylor approximation (blue dashed line) for the autogravitating profile. Parameters as in Table 8 for Tycho.
Figure 6. Numerical solution (full red line) and Taylor approximation (blue dashed line) for an NFW profile. Parameters as in Table 9 for Tycho.
5. Conclusion
The thin layer approximation in the framework of the conservation of energy is an alternative to the use of the conservation of momentum in order to find the equation of motion for a supernova remnant (SNR). In the case where the interstellar medium (ISM) has a constant density, it is possible to find the trajectory in an analytical form, see Equation (11). The case of energy conservation in a medium with variable density was also explored but an analytical trajectory was found only
Table 10. Synoptical parameters of the best model for SNRs with different density profiles.
in the case of a medium characterized by an inverse square decrease of density, see Equation (26). The other profiles of density require a numerical integration in order to find the trajectory. A Taylor series can provide the trajectory for a short interval of time: see Figure 2 for a power law, Figure 3 for an exponential law, Figure 4 for a Gaussian law, Figure 5 for an autogravitating medium and Figure 6 for a Navarro-Frenk-White (NFW) density profile. As an astrophysical target we have chosen to reproduce 4 standard SNRs. The match between the observed and simulated radius as well as that between the observed velocity and the simulated velocity has been analysed in terms of the percentage error, see Tables 2-9. Table 10 presents in column 2 the best model for the SNRs here analysed. The solution for the velocity to first order allows the insertion of the back reaction, i.e. the radiative losses, in the equation for the energy conservation, see Equation (18), and as a consequence the velocity corrected to second order, see Equation (19). The radiative losses allow evaluating the length at which the advancing velocity of the SNR is zero.