Traditionally, 1D models based on scaling laws have been used to parameterized convective heat transfer rocks in the interior of terrestrial planets like Earth, Mars, Mercury and Venus tackle computational bottleneck high-fidelity forward runs 2D or 3D. However, these are limited amount physics they can model (e.g. depth dependent material properties) predict only mean quantities such as mantle...