Numerical solution of a non-linear conservation law applicable to the interior dynamics of partially molten planets

Abstract

The energy balance of a partially molten rocky planet can be expressed as a non-linear diffusion equation using mixing length theory to quantify heat transport by both convection and mixing of the melt and solid phases. In this formulation the effective or eddy diffusivity depends on the entropy gradient, $\partial S/\partial r$, as well as entropy. First we present a simplified model with semi-analytical solutions, highlighting the large dynamic range of $\partial S/\partial r$, around 12 orders of magnitude, for physically-relevant parameters. It also elucidates the thermal structure of a magma ocean during the earliest stage of crystal formation. This motivates the development of a simple, stable numerical scheme able to capture the large dynamic range of $\partial S/\partial r$ and provide a flexible and robust method for time-integrating the energy equation. We then consider a full model including energy fluxes associated with convection, mixing, gravitational separation, and conduction that all depend on the thermophysical properties of the melt and solid phases. This model is discretised and evolved by applying the finite volume method (FVM), allowing for extended precision calculations and using $\partial S/\partial r$ as the solution variable. The FVM is well-suited to this problem since it is naturally energy conserving, flexible, and intuitive to incorporate arbitrary non-linear fluxes that rely on lookup data. Special attention is given to the numerically challenging scenario in which crystals first form in the centre of a magma ocean. Our computational framework is immediately applicable to modelling high melt fraction phenomena in Earth and planetary science research. Furthermore, it provides a template for solving similar non-linear diffusion equations arising in other disciplines, particularly for non-linear functional forms of the diffusion coefficient.

References in corpus (5)

• Scaling in thermal convection: A unifying theory
• Finite volume schemes for diffusion equations: introduction to and review of modern methods
• A Possible Carbon-rich Interior in Super-Earth 55 Cancri e
• Formation of Terrestrial Planets from Protoplanets under a Realistic Accretion Condition
• Differentiation of Vesta: Implications for a shallow magma ocean

Cited by in corpus (18)

• Atmosphere Origins for Exoplanet Sub-Neptunes
• Retention of water in terrestrial magma oceans and carbon-rich early atmospheres
• Linking the evolution of terrestrial interiors and an early outgassed atmosphere to astrophysical observations
• Vertically resolved magma ocean-protoatmosphere evolution: H${}_2$, H${}_2$O, CO${}_2$, CH${}_4$, CO, O${}_2$, and N${}_2$ as primary absorbers
• What factors affect the duration and outgassing of the terrestrial magma ocean?
• The influence of bulk composition on long-term interior-atmosphere evolution of terrestrial exoplanets
• Super-Earths and Earth-like Exoplanets
• Direct imaging of molten protoplanets in nearby young stellar associations
• Magma ocean evolution at arbitrary redox state
• Interior dynamics of super-Earth 55 Cancri e
• Thermal Evolution and magnetic history of rocky planets
• Convective shutdown in the atmospheres of lava worlds
• The evolution of catastrophically evaporating rocky planets
• Self-limited tidal heating and prolonged magma oceans in the L 98-59 system
• Assessing the Interior Structure of Terrestrial Exoplanets with Implications for Habitability
• Volatile-rich evolution of molten super-Earth L 98-59 d
• Evolution of Mercury's Earliest Atmosphere
• Mantle Differentiation, Mixing and Interior-Exterior Exchange