wd_nova_burstΒΆ
This test case checks the evolution of a nova outburst for one cycle.
This test case has 2 parts. Click to see a larger version of a plot.
- Part 1 (
inlist_setup
) loads1.1M_lgTc_7.7.mod
, a prebuilt carbon-oxygen white dwarf from the make_o_ne_wd test suite, sets the topical depth to 30, evolves for a few steps to an age of 1e7 year and terminates. - Part 2 (
inlist_wd_nova_burst
) continues the evolution with the accretion of a hydrogen-helium mixture at a rate of 1e-9 Msun/yr. Eventually hydrogen burning in the accerted envelope causes the luminosity to exceed 1e4 Lsun and an alert is written to the terminat via therun_star_extras.f90
:
num_bursts 1
405 8.277879 4.365E+05 8.624940 8.625037 1.115473 0.000025 0.000000 0.000000 0.309197 329.049014 1497 36
-5.0945E+00 7.679057 -1.661720 -2.464781 7.731723 -9.000000 1.115449 0.000000 0.400729 0.019380 65.113954 10
2.4766E+04 7.926520 4.191187 4.972103 -6.862099 -4.289137 1.115442 0.038204 0.474628 1.000000 0.327E-01 max increase
rel_E_err -1.9543755810854924D-08
log_rel_run_E_err -6.0755826271591920
and when the luminosity falls below 1e3 Lsun a nova cycle is considered finished and the run terminates:
585 7.835616 3.884E+05 1.624521 1.624523 1.115450 0.000001 0.000000 0.000000 0.332277 329.045159 1617 48
-1.4608E+00 7.679051 -2.147858 -16.859525 0.457865 -9.000000 1.115449 0.000000 0.400729 0.019454 65.113548 5
2.4768E+04 7.926494 3.015873 -3.748867 -67.588316 -3.897373 1.115442 0.038204 0.474628 1.000000 -0.133E-06 Lnuc_H
rel_E_err -5.3661942520223188D-10
log_rel_run_E_err -5.9590177286527863
have finished burst
pgstar commands used for the plots above:
&pgstar
file_white_on_black_flag = .true. ! white_on_black flags -- true means white foreground color on black background
!file_device = 'png' ! png
!file_extension = 'png'
file_device = 'vcps' ! postscript
file_extension = 'ps'
pgstar_interval = 1
Grid8_win_flag = .true.
Grid8_win_width = 12
Grid8_title = 'wd_nova_burst'
Summary_Burn_xaxis_name = 'logxq'
Summary_Burn_xaxis_reversed = .true.
Summary_Burn_xmin = -14 ! -101d0 ! only used if /= -101d0
Summary_Burn_xmax = -3 ! -101d0 ! only used if /= -101d0
Summary_Burn_title = ' '
Abundance_title = ' '
Abundance_xaxis_name = 'logxq'
Abundance_xaxis_reversed = .true.
Abundance_xmin = -14 ! -101d0 ! only used if /= -101d0
Abundance_xmax = -3 ! -101d0 ! only used if /= -101d0
Abundance_txt_scale = 0.8
Abundance_legend_max_cnt = 0
HR_title = ' '
HR_txt_scale = 0.8
HR_logT_min = 4.5
HR_logT_max = 6.0
HR_logL_min = -1.0
HR_logL_max = 5.0
TRho_title = ' '
TRho_txt_scale = 0.8
TRho_logT_min = -101d0
TRho_logT_max = -101d0
TRho_logRho_min = -101d0
TRho_logRho_max = -101d0
TRho_Profile_title = ' '
TRho_Profile_txt_scale = 0.7
TRho_Profile_xmin = -8.0
TRho_Profile_xmax = 9.0
TRho_Profile_ymin = 5.0
TRho_Profile_ymax = 9.0
Text_Summary1_name(5,2) = 'time_step_sec'
Text_Summary1_name(6,2) = 'max_tau_conv'
Grid8_file_flag = .true.
Grid8_file_dir = 'pgstar_out'
Grid8_file_prefix = 'grid8_'
Grid8_file_interval = 10000
Grid8_file_width = -1
Grid8_file_aspect_ratio = -1
/ ! end of pgstar namelist
Last-Updated: 08Jul2021 (MESA 094ff71) by fxt.