1.5M_with_diffusion

The test checks the functionality of element diffusion. The test vehicle is a 1.5 \({\rm M}_\odot\) solar metallicity model.

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

  • Part 1 (inlist_to_ZAMS) creates the pre-main-sequence model and evolves it to the main sequence. In the mass fraction plot, note the roughness of the 12C mass fraction profile (orange curve) from convection and that the surface hydrogen profile is flat (yellow curve).

../_images/abund_zams_0000181.svg

  • Part 2 (1.5M_with_diffusion) continues the evolution until the central hydrogen mass fraction drops below 0.01. In the mass fraction plot, note the 12C mass fraction profile (orange) has been smoothed by element diffusion and that the surface is nearly pure hydrogen (yellow upwards spike) from heavier elements diffusing inwards (e.g., orange downwards spike in 16O).

../_images/abund_h_depletion_0000430.svg

pgstar commands used for the plots above, e.g.:

&pgstar

   pgstar_interval = 50

   ! device
   file_white_on_black_flag = .true.
   file_device = 'vcps' ! postscript
   file_extension = 'ps'

   Abundance_win_flag = .false. ! switch to .true. to visualize in real-time
   Abundance_file_flag = .true.
   Abundance_file_dir = 'pgstar_out2'
   Abundance_file_prefix = 'abund_h_depletion_'
   Abundance_win_width = 9
   Abundance_win_aspect_ratio = 1.0d0

   Grid2_win_flag = .false. ! switch to .true. to visualize in real-time
   Grid2_file_flag = .false.

   Grid2_win_width = 15
   file_digits = 7
   Grid2_file_dir = 'pgstar_out2'
   Grid2_file_prefix = 'grid'
   Grid2_file_interval = 5
   Grid2_file_width = 15

/ ! end of pgstar namelist

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

Last-Run: 17Jul2024 (MESA 0cc374fa) by pmocz on C916PXT6XW in 423 seconds using 8 threads.