This repository provides a script to be used with dedalus to model Rayleigh-Bénard convection using a fully compressible model. The equation of state (EOS) is of the Murnaghan type and applies well to planetary interiors. The Prandtl number is assumed infinite, as relevant to solid planetary mantles.
Details of the model and results can be found in the paper
Yanick Ricard, Thierry Alboussière, Stéphane Labrosse, Jezabel Curbelo, and Fabien Dubuffet. Fully compressible convection for planetary mantles Geophys. J. Int. (2022), 30, 932-956. https://doi.org/10.1093/gji/ggac102