[Mesa-users] How to extract values

Fri Jul 22 17:51:20 UTC 2022

hi harish,

> ... simulate upto log Tc = 9.5.
> 

> ... recreate the results of paper "The outcomes of carbon–oxygen white dwarfs accreting CO-rich material " https://doi.org/10.1093/mnras/sty3176 in which they achieved that temperature ranges.

you are doing their case 1, with mdot < 2.45e-5 msun/yr.
their figure 2 shows log(Tc) ~8.7, about the same as you find, and does not reach log(T) ~ 9.5.

once carbon burning gets going in earnest, the model is finished (timesteps get small) 
as mesa currently cannot follow the deflagration front as it propagates toward the surface. 

fxt

> On Jul 22, 2022, at 9:27 AM, harish sorrout via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> Thank you Rob, 
> 
>  Now I have added the extra profile columns in run_star_extras.f90. By using mesa 12778, I have accreted the carbon and oxygen on to  CO white dwarf with accretion rate 2.3d-6. Although white dwarf accreted 1.369 solar mass and central temperature log 8.69 and then simulation terminate with following errors:
> num_newton_iterations dt_limit        1.1642857143        17        14
>                                       log eps nuc ratio         368         764         164    5.0531126467963694D-09   -5.0531126467963694D-09    1.
> 5160424637814009D+01    1.5160424814208914D+01
>   timestep_hold > model_number, so no timestep increase         165         164
>                                  WARNING: rel_run_E_err         164    4.1014717634964032D-02
>  retry     165
>  retry     165
>  retry     165
>  retry     165
>  retry     165
>  retry     165
>  retry     165
> num_newton_iterations dt_limit        1.2357142857        18        14
>                                       log eps nuc ratio         362         775         165   -1.6962520136639194D-12    1.6962520136639194D-12    8.
> 3615846268716751D+00    8.3615846268390168D+00
>   timestep_hold > model_number, so no timestep increase         166         165
>                                  WARNING: rel_run_E_err         165    4.1014909100627141D-02
>  png/trho_000165.png/png
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>  retry     166
>                                                      dt    8.4894335288525781D-16
>                                      min_timestep_limit    1.0000000000000001D-15
> 
>  stopping because of problems dt < min_timestep_limit
> 
>                                  WARNING: rel_run_E_err         166    4.1014909100627141D-02
> 
> terminated evolution: cannot find acceptable model
> 
> termination code: min_timestep_limit
> DATE: 2022-07-22
> TIME: 21:46:20
> 
> ******************************************************************
> failed to create lgT_8.45.mod when running inlist_to_c_ignition_header
> ******************************************************************
> 
> SIr, I know it is difficult to further simulate upto log Tc = 9.5. but I want to recreate the results of paper "The outcomes of carbon–oxygen white dwarfs accreting CO-rich material " https://doi.org/10.1093/mnras/sty3176 in which they achieved that temperature ranges. Also I have checked their published inlist. 
> So, I am requesting that you and all mesa users please suggest what parameters I should use in my inlist and what physics I am missing  to achieve that temperature range.  Please do me a favour.
> 
> Thankyou 
>  
> with regards 
> Harish Kumar
> 
> On Tue, Jul 19, 2022 at 8:11 PM Rob Farmer <robert.j.farmer37 at gmail.com> wrote:
> > By doing so, I did not get any errors while compiling. Sir Is there anything wrong with the above command. 
> Sir I could not understand your words  "You want to pass the string "burn_c" not a variable called burn_c".
> Please enlighten me.
> 
> Yes your still passing a character called burn_c and not the string 'burn_c' (with the quotes)
> 
> E_cburn(nz) = star_get_profile_output(s,'burn_c',k)
> 
> Rob
> 
> 
> On Fri, 15 Jul 2022 at 17:38, harish sorrout via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> ---------- Forwarded message ---------
> From: harish sorrout <harishsorrout at gmail.com>
> Date: Fri, 15 Jul 2022, 21:02
> Subject: Re: [Mesa-users] How to extract values
> To: Rob Farmer <robert.j.farmer37 at gmail.com>
> 
> 
> Sir 
> As you told to use 
>             real(dp) function star_get_profile_output(s,name,k) to get the value of 'burn_c' .  So for this I need to put 'burn_c' at the place of 'name' argument. 
> So i select character (len=maxlen_profile_column_name) :: burn_c.  And use
>   E_cburn(nz) = star_get_profile_output(s,burn_c,k) where real(dp) :: E_cburn(nz). and 
> then used burn_timescales = (s% cp(k) * s% T) / E_cburn(k)
> 
> By doing so, I did not get any errors while compiling. Sir Is there anything wrong with the above command. 
> Sir I could not understand your words  "You want to pass the string "burn_c" not a variable called burn_c".
> Please enlighten me.
> 
> Harish Kumar
> 
> 
> 
> 
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 
> <inlist_to_lgT_8.45><0.85Mwd_0.5crelax.mod>_______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 