[FLASH-USERS] Cooling module and shocks

[Michiel Bustraan]

Dear FLASH Users,


I've been using the cooling module in my supernova simulations, and have been seeing some unusual results.


My simulations involve a shock moving radially outward into a circumstellar medium.

My runs without cooling have good results, with the structure of the shock matching those in the literature.

The attached plots are the density, temperature, and radial velocity from a run without cooling.


I've been looking into the cooling module and implementing a fairly simple cooling function.

While the cooling seems to be exactly as I would expect (I plot the variable solnData(COOL_VAR,i,j,k) = abs(dedt) against my cooling function),

the impact the cooling has on the shock doesn't seem right.

The cooling causes an extreme narrowing of the shell between the forward and reverse shock, while the reverse shock almost disappears into the shell.


I wanted to ask if anyone else encountered something like this using the cooling module.


Also, just for assurance, I wanted to check what the cooling module actually calculates.

Am I correct that the variable dedt, defined by the cooling function in cool_deriv, in the Cool.F90 file, should be the cooling rate (Lambda [erg cm^3 s^-1]) multiplied by the number density squared?

And the variable EINT_VAR is the internal energy per mass [erg g^-1]? Meaning it has to be multiplied by mass density to be the energy density.

?

I know the question is a bit vague, but any advice would be appreciated.


Thank you,

Michiel Bustraan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20170731/2a10bf58/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: density_1d_20.png
Type: image/png
Size: 25614 bytes
Desc: density_1d_20.png
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20170731/2a10bf58/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: temp_1d_20.png
Type: image/png
Size: 22640 bytes
Desc: temp_1d_20.png
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20170731/2a10bf58/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: velx_1d_20.png
Type: image/png
Size: 21305 bytes
Desc: velx_1d_20.png
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20170731/2a10bf58/attachment-0002.png>