Thermo-chemical convection is the primary process that controls the large-scale dynamics of the mantle of the Earth and terrestrial planets. Its numerical simulation is one the principal tools for exploiting the constraints posed by geological and geochemical surface observations performed by planetary spacecrafts. In the present work, the authors discuss the modeling of active compositional fields in the framework of solid-state mantle convection using the cylindrical/spherical code Gaia. They compare an Eulerian method based on double-diffusive convection against a Lagrangian, particle-based method. Through a series of increasingly complex benchmark tests, the authors show the superiority of the particle method when it comes to model the advection of compositional interfaces with sharp density and viscosity contrasts. They finally apply this technique to simulate the Rayleigh-Taylor overturn of the Mars’ and Mercury’s primordial magma oceans.