plotting local chi ("cooling strength")

I was a bit skeptical that calculating cooling lengths from 1D radiative shock simulations would directly translate to these 3D jet simulations, so I implemented a field variable to calculate the cooling length on the fly for each cell. Doing it this way results in chi's that are very different. Here is a pseudocolor plot of chi from what I have been calling weak cooling.

I had estimated a chi of approximately 0.4 using my 1D radiative shock simulations, but as you can see here chi can locally be much different. However, you sort of have to take the ambient with a grain of salt. You get extremely low chi's in those regions because nothing is moving, and the material has not been shocked yet. Here is a lineout of chi along the center of the jet.

You can see that the ambient is very low, and then chi does some interesting things as it goes through the bow shock, and then it becomes constant once you enter the jet beam. For most of the jet chi > 1. Here is a zoom in of the chi lineout at the head of the jet.

In this specific region chi does fall between 0.1 and 1, so perhaps a chi of 0.4 is a good approximation after all.


UPDATE

The above plots are not the original post; they are more accurate now.

Also, since plotting hte cooling length of the ambient is kind of strange, I also decided to plot cooling time. Now the ambient looks exactly the same as the jet beam…

Attachments (4)

Download all attachments as: .zip

Comments

No comments.