Version 2 (modified by 12 years ago) ( diff ) | ,
---|

# 1D Radiative Shocks

These simulations follow Ch. 4 of Delamarter '01 and his treatment of the 1D steady radiative shock problem. The purpose of this problem is to check that the cooling source terms are being handled correctly. Here are some useful links:

## Equations

The shock jump equations for a stationary shock are used to solve for the initial post-shock values. The post-shock velocity v_{2} can be written as:

where v_{1} is the ambient velocity, and M is the ambient mach number. Remember that the mach number M = v_{1}/c where c is the ambient sound speed, and where P_{1} is ambient pressure, and is ambient density. The post-shock density and pressure ( and P_{2} respectively) can be found by using mass flux and momentum flux conservation across the shock:

These post-shock values become the boundary conditions for the fluid equations in the cooling region:

where

, v, and P represent the density, velocity, and pressure in the cooling region as functions of x, and is the cooling rate.

## Initial Parameters

n_{1} = 60 particles/cc

v_{1} = 10^{7} cm/s

T_{1} = 10^{4} K

For analytic cooling,

where n and T and the number density and temperature in the post-shock region respectively.

For these simulations,

= 2

^{-34} erg*cm^{3}/s/K^{2}

cell length = 2.5 x 10^{15} cm

problem domain = 400 cells

final time ~ 4000 years

## Results

A run with no cooling was done to make sure the post-shock values were correct. With no cooling (adiabatic), the hydrodynamic quantities should jump discontinuously at the shock, and then remain constant while the shock itself remains stationary. Only the first 100 cells of the simulation are shown here.

Adiabatic | |

Pressure Movie | Temperature Movie |
---|---|

Density Movie | Velocity Movie |

Then, a run with the aforementioned analytic cooling parameters was done. The hydrodynamic quantities will change discontinuously at the shock as before, but then change continuously according to the cooling rate. The ambient temperature is set to be the floor temperature, so that no cooling occurs at or below this temperature. Once the post-shock temperature reaches this floor temperature, the hydrodynamic quantities will remain constant. This is called the quiescent region. Only the first 100 cells of the simulation are shown here.

Analytic Cooling | |

Pressure Movie | Temperature Movie |
---|---|

Density Movie | Velocity Movie |

For Dalgarno-McCray (DM) cooling, some of the initial parameters had to be changed. The initial velocity used in the previous simulations would give rise to a post-shock temperature in the unstable region of the DM curve. This led to the radiative instability, and thus unsteady shocks. To avoid this, the initial velocity was lowered to 80 km/s. Also, the cooling rate from the DM curve is much higher than that of the analytic form. In order to see the cooling region, the grid had to be shrunk, resulting in a cell length of 10^{13} cm.

DM Cooling |

Temperature Movie |
---|

For Non-equilibrium (NEQ) cooling, the length scale had to be adjusted. The shock velocity is the same as the DM case, but the cell length is now 7.75e11 cm.

NEQ Cooling |

Temperature Movie |
---|

Z Cooling has the same initial parameters and scaling as the NEQ case.

Z Cooling |

Temperature Movie |
---|

Here is a fun plot that shows the temperature profiles for DM, NEQ, and Z cooling. This gives you a good idea of the temperature ranges that each type of cooling works in. Also, you can see the differences in the cooling length which is why I would normally change the length scale depending on the type of cooling.

## Including MHD

Making these simulations work with MHD was no simple task. Having a magnetic field adds another jump condition, and it makes the equations within the cooling region more complicated. Here are what the equations look like now:

### Attachments (19)

- ad_press.gif (87.2 KB ) - added by 12 years ago.
- ad_press0000.png (11.6 KB ) - added by 12 years ago.
- ad_rho.gif (97.0 KB ) - added by 12 years ago.
- ad_rho0000.png (13.9 KB ) - added by 12 years ago.
- ad_temp.gif (91.8 KB ) - added by 12 years ago.
- ad_temp0000.png (12.6 KB ) - added by 12 years ago.
- ad_vel.gif (95.8 KB ) - added by 12 years ago.
- ad_vel0000.png (14.0 KB ) - added by 12 years ago.
- an_press.gif (85.1 KB ) - added by 12 years ago.
- an_press0000.png (11.5 KB ) - added by 12 years ago.
- an_rho.gif (97.8 KB ) - added by 12 years ago.
- an_rho0000.png (14.1 KB ) - added by 12 years ago.
- an_temp.gif (94.3 KB ) - added by 12 years ago.
- an_temp0000.png (13.1 KB ) - added by 12 years ago.
- an_vel.gif (97.2 KB ) - added by 12 years ago.
- an_vel0000.png (14.3 KB ) - added by 12 years ago.
- DM_temp.gif (1.2 MB ) - added by 12 years ago.
- DM_temp0010.png (19.3 KB ) - added by 12 years ago.
- allcool_temp.png (35.6 KB ) - added by 12 years ago.

**Note:**See TracWiki for help on using the wiki.