User:Tohline/ProjectsUnderway/Core Collapse Supernovae
From VisTrailsWiki
Contents 
A Template for Gravitational Wave Signals from CoreCollapse Supernovae
 Tiled Menu  Tables of Content  Banner Video  Tohline Home Page  
29 August 2014: The initial content of this chapter has been drawn from a LaTeX document dated 23 December 2009 that J. E. Tohline created while interacting with Cody Arceneaux (at the time, an LSU undergraduate physics major) and Sarah Caudell (at the time, an LSU physics graduate student). I've decided to build on this project foundation that was laid in the 2008  2009 time frame.
Introduction
Our principal objective is to help the gravitationalwave community more fully understand the underlying physics that is fundamentally responsible for the characteristic features that are expected to arise in the signals that are detected from corecollapse supernovae. Building on early, simplified models of nonspherical core collapse presented by Thuan & Ostriker (1974) and by Tohilne (1984), our discussion is intended to supplement and complement published works from the past couple of decades that have focused on analyzing results from largescale, multidimensional hydrodynamic (or magnetohydrodynamic and fully relativistic) models, such as the the seminal work of Zwerger and Mueller (1997) and reviews by Ott (2009), Fryer and New (2011) and Logue et al. (2012).As an example, Figure 2 from Ott (2009) is reprinted here, on the right. Its published caption reads, in part: "Note the generic shape of the waveforms, exhibiting one pronounced spike at core bounce and a subsequent ring down. Very Rapid precollapse rotation … results in a significant slowdown of core bounce, leading to a loweramplitude and lowerfrequency GW burst." As early as 1991, Mönchmeyer et al. (1991) identified four phases of a collapse evolution:
 Phase I: Contraction and core flattening. "This phase … is characterized by a growing overall flattening of the density stratification inside and outside [of an] inner core (IC)."
 Phase II: Rapid contraction of the IC. "Both the IC and the outer core get strongly deformed during this phase. … The collapse of the IC is finally decelerated and stopped … by either the action of centrifugal forces … or the rise of [the effective adiabatic index] …"
 Phase III: Core bounce. "Depending on the stiffness of the [equation of state] and the amount of rotational and kinetic energy the IC overshoots its (rotational) equilibrium position" then bounces back.
 Phase IV: Shock propagation and postbounce oscillations. "… the … IC oscillates [about its equilibrium configuration] with various volume and surface modes. The amplitudes and frequencies of these modes strongly depend on the kinetic energy of the IC at bounce, the stiffness of the [equation of state], and the central and averaged density of the IC …. The ICoscillations are damped … [on a] time scale [that] strongly depends on the conditions of the shock heated material just outside the IC. After a model dependent number of oscillations, the IC approaches rotational equilibrium …."
Here, we will develop analytic and semianalytic models that encapsulate the physical phenomena that give rise to and that define the relevant scales of each of these evolutionary phases, then sew the models together to demonstrate how all of the principal features in the corecollapse gravitationalwave signal can be fairly readily understood. We begin with a discussion of the freefall collapse of rotationally flattened spheroids.
FreeFall Collapse
Nod to LyndenBell's Early Contributions
I want to begin this section by paying tribute to Donald LyndenBell who, in 1962 (Mathematical Proceedings of the Cambridge Philosophical Society, vol. 58, pp. 709711), was the first to appreciate the relatively simple behavior that should be exhibited by the freefall collapse of a uniformly rotating, uniformdensity spheroid. In an article less than two pages in length, LyndenBell first noted that the governing dynamical equations (written in cylindrical coordinates) take the form,






where the dots denote differentiation with respect to time, is a constant and the two timedependent coefficients, and , come from the gravitational potential that, according to Lyttleton (1953), has the form,



Then LyndenBell deduced that, (a) "the result of the motion is merely a change of scales"; (b) the collapsing system "remains uniform [in density], and the boundary remains spheroidal"; and (c) "the collapse … will be through a series of [uniformly rotating] spheroids." Two years later in a separate article, LyndenBell (1964, ApJ, 139, 1195) presented results from the numerical integration of this governing set of dynamical equations (see his Figure 1, reprinted to the right, here). The various publications by other authors who also have modeled the freefall collapse of rotating or nonrotating spheroids in various contexts (see our discussion that follows) have not always acknowledged LyndenBell's pioneering analysis of this problem.
Nonrotating, Spherically Symmetric Collapse
When describing the freefall (pressurefree) collapse from rest of a uniformdensity sphere of mass, , and initial radius, — hence, initial density — it is convenient to use the Lagrangian radial coordinate, , which tracks the radius of the sphere at any time, . The relevant equation of motion is (see Lin, Mestel & Shu 1965 or our accompanying discussion),






which integrates to give,
where,






(Detailed steps through this derivation, as well as an indepth discussion of the freefall collapse problem, are provided in an accompanying chapter of this H_Book.) A set of dimensionless expressions drawn from this analytic solution that will prove useful to our discussion of gravitationalwave signals from corecollapse supernovae is provided in the following table. All lengths have been normalized to , and times have been normalized to the freefall time, . In the lower half of the table, these analytic functions have been evaluated at six different times during the freefall collapse.
Analytic Solution to Spherical FreeFall Collapse 


 
















































Note that the second timederivative of the moment of inertia, , goes through zero when , that is, at time,
Nonrotating, OblateSpheroidal Collapse
The analogous collapse from rest of a nonrotating, uniformdensity, oblate spheroid with equatorial radius, , and polar radius, , is governed by the equations,






where, to within an additive constant,



We should clarify and emphasize that this expression for the timedependent gravitational potential has been written in terms of the (timevarying) eccentricity of the spheroid, , as measured in the meridional plane. Specifically,






and, as is derived in our accompanying discussion of the properties of homogeneous ellipsoids,






Hence,



and the pair of governing dynamical equations become,






where we have used the relation that is valid for uniformdensity, oblate spheoids,
Related Discussions
© 2014  2021 by Joel E. Tohline 