Numerical integration was performed in Python, using the RK45 method with the solve_ivp routine having an accuracy of 10−15 for vanishing the dimensionless pressure at the surface. 5.2.5 Cracking and convective instability We wanted to explore the incidence of convective instability with cracking...