Figure 1: Cooling curves for material of solar abundance assuming collisional ionization equilibrium (solid) or non-equilibrium ionization (dashed). The former was calculated using the MEKAL thermal plasma code (Mewe et al.1985; Kaastra 1992) distributed in XSPEC (v11.2.0), while the latter was taken from data in Sutherland & Dopita (1993). These curves are normalized so that the net cooling rate per unit volume, , where is the electron number density and is the total number density of all of the ions. | |
Open with DEXTER |
Figure 2: The slope, , of the respective cooling curves in Fig. 1 (CIE - solid; NEI - dashed). | |
Open with DEXTER |
Figure 3: Schematic of the grid set up for implusive shock generation. Supersonic flow with density and flow speed collides with a subsonic CDL with density and flow speed . The upstream boundary condition is inflow, while the downstream boundary condition is outflow. All simulations are 1-D. | |
Open with DEXTER |
Figure 4: Time-space diagrams of the density evolution of a 1-dimensional Mach 1.4 radiative shock with different cooling exponents . Supersonic flow enters the grid from the top, and the cooled postshock gas flows off the grid at the bottom. Lighter shades indicate lower densities. Distances are marked in units of (the value of this is different in each panel - see Eq. (4) and Table 2) while time is shown in units of ( ). The value of is -2.5, -2.0, and -1.5 in the top, middle, and bottom panels respectively. | |
Open with DEXTER |
Figure 5: Time-space diagrams of the density evolution of a 1-dimensional radiative shock with cooling exponents and different Mach numbers. M is varied from 1.4, 2, 3 and 5 from the top to bottom panels. Note that the axis and density scaling changes from panel to panel, unlike in Fig. 4 where the scaling is kept constant. | |
Open with DEXTER |
Figure 6: Power spectra (mean square amplitude) for the simulations shown in Fig. 5 where . M varies from 1.4, 2, 3, 5 from the top to bottom panels respectively. | |
Open with DEXTER |
Figure 7: Time-space diagrams of the density evolution of a 1-dimensional radiative shock with M=3 and cooling exponent . In the top panel the flow is initialized to the steady state solution. In the middle panel the shock is impulsively generated. In the bottom panel the downstream boundary condition is reflecting. | |
Open with DEXTER |
Figure 8: The shock position as a function of time for simulations with M=3 and . The set up with an isolated steady-state shock as the initial condition is shown at the top, the set up with impulsive shock generation is shown in the middle, and the set up with a reflecting downstream boundary ( i.e. "flow into a wall'') is shown at the bottom. In the latter case the growth in width of the cold dense layer has been removed to aid comparison with the other scenarios, and each dataset is offset with respect to the others. | |
Open with DEXTER |
Figure 9: Time-space diagrams of the density evolution of a 1-dimensional radiative shock with cooling exponent and M=5 ( top) and M=10 ( bottom). | |
Open with DEXTER |
Figure 10: The time-evolution of the shock position as a function of when M=1.4 and . | |
Open with DEXTER |
Figure 11: The shock position as a function of time when M=1.4 and . The ratio of the temperature of the CDL to the pre-shock temperature, , is noted in each panel. | |
Open with DEXTER |
Figure 12: The value of , the critical value for damping of the overstability, as a function of the reflection coefficient, R, for four different values of . The lines are only meant to guide the eye. | |
Open with DEXTER |