music-mesa-tables 0.2.1

MESA tables for post-processing of MUSIC simulations

music-mesa-tables

Fast implementation of the 'mesa' equation of state as used in the MUSIC project.

MUSIC is a stellar convection code. See e.g. Viallet et al (2016), Goffrey et al (2017).

The so-called 'mesa' equation of state implemented in MUSIC relies on tables from the MESA project (version 15140), reinterpolated in (density, internal energy) space. In the regimes typically simulated in MUSIC, those tables come from the work of Rogers and Nayfonov (2002). The reinterpolated tables are bundled in this package.

This package is intended for post-processing of data from MUSIC simulations. To build stellar evolution models, refer to the MESA project.

---

This offers a Python API. Here is a simple example:

import music_mesa_tables as mmt
import numpy as np

eos = mmt.CstCompoEos(metallicity=0.02, he_frac=0.28)

density = np.array([1.05, 15.7, 134.9])
internal_energy = np.array([1e12, 1e15, 3e15])
state = mmt.CstCompoState(eos, density, internal_energy)

temperature = 10 ** state.compute(mmt.StateVar.LogTemperature)
# array([2.21745558e+03, 5.00852231e+06, 1.48317986e+07])