Pulsed Jet Simulations and Cooling Length

I'm setting up a diagnostic variable, so that I can see what the cooling length is everywhere on the grid in my pulsed jet simulations.

I have… (cooling length) = (velocity)*(internal energy) / (cooling strength)

Where cooling strength is in units of energy/time.

This should give the value to me in computational units. My simulation has 30 cells per cu at the highest level. Say I want at least 10 cells per cooling length, so if my cooling length drops below 0.33 then I'm under-resolved.

The only thing I"m not sure about is if it's correct to use the velocity, because it will be 0 in the ambient. So the cooling length would be zero everywhere, but maybe that's ok and we just ignore the ambient?

Below is an early snapshot of the hydro run with the Raga set up. The hydro case should have the strongest cooling since magnetic fields tend to puff up cooling regions. You can see that the cooling length starts at 0 because it's in the ambient. Using this criteria the region directly behind the bow shock, and behind that first internal shock appear to be under-resolved.

Attachments (1)

Download all attachments as: .zip

Comments

No comments.