Numerical approaches to solid-state convection in planetary bodies