  Over the past decades, global geodynamical models have been used to
investigate the thermal evolution of terrestrial planets. With the increase of
computational power and improvement of numerical techniques, these models have
become more complex, and simulations are now able to use a high resolution 3D
spherical shell geometry and to account for strongly varying viscosity, as
appropriate for mantle materials. In this study we review global 3D geodynamic
models that have been used to study the thermal evolution and interior dynamics
of Mars. We discuss how these models can be combined with local and global
observations to constrain the planet's thermal history. In particular, we use
the recent InSight estimates of the crustal thickness, upper mantle structure,
and core size to show how these constraints can be combined with 3D geodynamic
models to improve our understanding of the interior dynamics, present-day
thermal state and temperature variations in the interior of Mars.
