Cooling Examples¶
Cooling Rate¶
This example is generated in examples/cooling_rate with the Grackle cooling.
The plot are the same than in [Smith2016] and uses a solar metallicity.
The code has two different modes 1D or 2D. In the first mode, the code generates particles at different temperatures but same density and computes the cooling rate associated to theses conditions.
In the second mode, the code generates a grid of particles at different density and temperatures and then computes the cooling. I suppose that in [Smith2016], they use the non equilibrium mode of Grackle to generate this graph. It requires a careful initialization of the element fractions which is not easy due to the non monotonic behavior of the cooling function. For simplicity, only Grackle with the equilibrium mode is shown in this documentation.
