.. _16M_predictive_mix:

******************
16M_predictive_mix
******************

.. tags:: star, predictive-mixing, convective-core, massive-star, main-sequence

This test suite example re-creates the 16-solar mass main-sequence
evolution with the inclusion of predictive mixing (using the Ledoux
criterion), as detailed in Section 2 of the MESA IV instrument paper
(Paxton et al 2018).

This test case has two parts. Click to see a larger view of a plot.

* Part 1 (``inlist_start``) creates a 16 |Msun| pre-main-sequence model and evolves it for 10 time steps.

* Part 2 (``inlist_16M_predictive_mix``) continues the evolution until core hydrogen depletion (mass fraction h1_center < 1e-6).

A Kippenhahn diagram shows the evolution of a retreating convective core on the main sequence, the blue region between model numbers 240 and 385.
This can be compared to Figure 3 in the MESA IV instrument paper, and the convective pre-mixing Kippenhahn diagram in :ref:`16M_conv_premix`.

.. image:: ../../../star/test_suite/16M_predictive_mix/docs/kipp_00000385.svg
   :width: 100%


pgstar commands used for the plots above:

.. code-block:: console

 &pgstar

   file_white_on_black_flag = .true. ! white_on_black flags -- true means white foreground color on black background
   file_device = 'png'            ! png

   !file_device = 'vcps'          ! postscript

    pgstar_interval = 10
    file_digits = 8


   Kipp_win_flag = .true.

   ! window properties
   Kipp_win_width = 12
   Kipp_win_aspect_ratio = 0.75
   Kipp_txt_scale = 0.9
   Kipp_title = ''

   ! y axis limits
   Kipp_mass_max = 16.0
   Kipp_mass_min = 0
   Kipp_show_mass_boundaries = .true.

   ! x axis limits
   Kipp_xaxis_name = 'model_number'
   Kipp_xmax = -101              ! maximum step number.  negative means use default.
   Kipp_xmin = 0         ! minimum step number.  negative means use default.

   Kipp_show_mixing = .true.
   Kipp_show_burn = .true.
   Kipp_show_luminosities = .true.

   ! file output
   Kipp_file_flag = .true.
   Kipp_file_dir = 'kipp_png'
   Kipp_file_prefix = 'kipp_'
   Kipp_file_interval = 10     ! output when mod(model_number,file_interval)==0
   Kipp_file_width = 12        ! (inches) negative means use same value as for window
   Kipp_file_aspect_ratio = -1 ! negative means use same value as for window

 / ! end of pgstar namelist


Last-Updated: 28May2021 (MESA ebecc10) by fxt