This repository accompanies the paper An entropy method for geodynamic modeling of phase transitions: Capturing sharp and broad transitions in a multi-phase assemblage by Juliane Dannberg, Rene Gassmoeller, Ranpeng Li, Carolina Lithgow-Bertelloni, Lars Stixrude This repository contains ASPECT data, and model setups. To run the models in this repo follow the steps: The global spherical convection models of the publication presented in Section 3.3 were run using the following dependencies: ----------------------------------------------------------------------------- -- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion. -- . version 2.4.0-pre (main, 8dcb250) -- . using deal.II 9.4.0-pre -- . with 64 bit indices and vectorization level 2 (256 bits) -- . using Trilinos 12.18.1 -- . using p4est 2.3.2 ----------------------------------------------------------------------------- This configuration was also used to compute the convection_box_opx.prm model in the compare_phase_function_and_table folder (Table model in Section 3.2). Running these models requires plugins that are located in the aspect/benchmarks/entropy_adiabat/plugins/ folder and that need to be compiled. Navigate into this directory and follow the steps: cmake -D Aspect_DIR=PATH_TO_ASPECT (replace PATH_TO_ASPECT with the directory where you compiled ASPECT). make Now the models in this repository can be started. You can start them with the ASPECT executable in your build folder. The models that reproduce the setup of Christensen & Yuen (1985) presented in Section 3.1 require a certain ASPECT branch because the adiabatic profile is computed in a slightly different way in Christensen & Yuen (1985) compared to how ASPECT usually does it. The easiest way to get the correct branch is to follow these git commands in your main ASPECT directory: git remote add gassmoeller https://github.com/gassmoeller/aspect.git (This is only necessary once and gives an error if the remote already exists, if so nothing happens) git remote update gassmoeller git checkout gassmoeller/solve_entropy_equation Recompile ASPECT as usual, with cmake . and make -jX (replace X by number of cores) Now the models in the other directories can be started. You can start them with the executable in your build folder. They were run using the following dependencies: ----------------------------------------------------------------------------- -- . using deal.II 9.2.0 -- . with 32 bit indices and vectorization level 2 (256 bits) -- . using Trilinos 12.18.1 -- . using p4est 2.2.0 ----------------------------------------------------------------------------- This configuration was also used to compute the convection_box_phase_function.prm model in the compare_phase_function_and_table folder (corresponding to the three phase function models in Section 3.2). The models correspond to the manuscript sections as follows: Section 3.1: The effect of phase transition thickness: christensen_yuen_1985_Ra_width Section 3.2: The effect of latent heat: compare_phase_function_and_table Section 3.3: Spherical convection models: spherical_convection_models The directories contain parameter files that should run individually for one of the setups of the publication. Some also contain shell scripts that will start all models of a model series. Some also contain python scripts to generate the figures. The models presented in Appendix A are part of ASPECT's main repository and are located in the folder aspect/benchmarks/entropy_adiabat/. entropy_adiabat.prm corresponds to A1 Advection along an adiabat entropy_half_space.prm corresponds to A2 Conductive cooling of a half-space entropy_conduction.prm corresponds to A3 Conductive cooling of an adiabat