Critical heat flux in a fuel rod: a boiling crisis with solve_ivp

Afterwards you can simulate heat flow through a fuel rod's pellet, gap, and cladding with solve_ivp and follow a boiling crisis with terminal events.