[mesa-users] accretion_powered_irradiation

Chen WC chenwc at pku.edu.cn
Sat Jun 11 13:22:23 EDT 2016


Hi, Bill,
Thank you for your response.

My inlist_project file is as follows:
******
&binary_job

   inlist_names(1) = 'inlist1'
   inlist_names(2) = 'inlist2'

   evolve_both_stars = .false.

/ ! end of binary_job namelist

&binary_controls

   m1 = 3.0d0  ! donor mass in Msun
   m2 = 1.4d0 ! companion mass in Msun
   initial_period_in_days =1.0d0

   fr = 0.05
   fr_limit = 1.0d-2

   limit_retention_by_mdot_edd = .true.
   accretion_powered_irradiation = .true.
   col_depth_for_eps_extra = 100.0
   accretor_radius_for_irrad =1.0d6
!   max_F_irr = 5.0d11

   use_other_jdot_mb = .true.
!   use_other_jdot_ml = .true.
   use_other_mdot_edd = .true.
!   keep_mb_on = .true.
   max_tries_to_achieve = 50
   mdot_scheme = 'Ritter'
/ ! end of binary_controls name list
*****

Some output results in the terminal are as follows:

******
hydro_newton_step failed to converge
       retry cnt, step, log10(dt/yr), retry_factor     1  2737       -0.788149        0.500000
                   1                2737     2.216895761E+08     2.117224534E+00    -5.713074589E+00     6.755361544E-03        photos1/x700
                                               lg F_irr        2737    1.2698970004336019D+01
                   1                2738     2.216895765E+08     2.117223777E+00    -5.713074589E+00     6.749303810E-03        photos1/x700
                                               lg F_irr        2738    1.2698970004336019D+01
                   1                2739     2.216895770E+08     2.117222869E+00    -5.713074589E+00     6.752663357E-03        photos1/x700
                                               lg F_irr        2739    1.2698970004336019D+01
                   1                2740     2.216895775E+08     2.117221779E+00    -5.713074589E+00     6.763619824E-03        photos1/x700
                                               lg F_irr        2740    1.2698970004336019D+01
                   1                2741     2.216895782E+08     2.117220471E+00    -5.713074589E+00     6.776681402E-03        photos1/x700
                                               lg F_irr        2741    1.2698970004336019D+01
                   1                2742     2.216895790E+08     2.117218902E+00    -5.713074589E+00     6.785790596E-03        photos1/x700
                                               lg F_irr        2742    1.2698970004336019D+01
                   1                2742     2.216895790E+08     2.117218823E+00    -5.691885290E+00     6.706529261E-03        photos1/x700
                                               lg F_irr        2742    1.2698970004336019D+01
                   1                2742     2.216895790E+08     2.117218823E+00    -5.691885290E+00     6.776681402E-03        photos1/x700
hydro_newton_step returned ierr
       retry cnt, step, log10(dt/yr), retry_factor     1  2742       -0.693273        0.500000
                   1                2742     2.216895786E+08     2.117219686E+00    -5.713074589E+00     6.770370486E-03        photos1/x700
                                               lg F_irr        2742    1.2698970004336019D+01
                   1                2743     2.216895791E+08     2.117218745E+00    -5.713074589E+00     6.770294906E-03        photos1/x700
                                               lg F_irr        2743    1.2698970004336019D+01
                   1                2744     2.216895797E+08     2.117217615E+00    -5.713074589E+00     6.774228347E-03        photos1/x700
                                               lg F_irr        2744    1.2698970004336019D+01
                   1                2745     2.216895804E+08     2.117216259E+00    -5.713074589E+00     6.779083382E-03        photos1/x700
                                               lg F_irr        2745    1.2698970004336019D+01
                   1                2746     2.216895812E+08     2.117214632E+00    -5.713074589E+00     6.796953754E-03        photos1/x700
                                               lg F_irr        2746    1.2698970004336019D+01
                   1                2746     2.216895812E+08     2.117214551E+00    -5.691885290E+00     6.718226481E-03        photos1/x700
                                               lg F_irr        2746    1.2698970004336019D+01
                   1                2746     2.216895812E+08     2.117214593E+00    -5.702907945E+00     6.758997650E-03        photos1/x700
                                               lg F_irr        2746    1.2698970004336019D+01
                   1                2747     2.216895822E+08     2.117212595E+00    -5.702907945E+00     6.746944676E-03        photos1/x700
                                               lg F_irr        2747    1.2698970004336019D+01
                   1                2748     2.216895835E+08     2.117210196E+00    -5.702907945E+00     6.725341516E-03        photos1/x700
                                               lg F_irr        2748    1.2698970004336019D+01
                   1                2749     2.216895849E+08     2.117207318E+00    -5.702907945E+00     6.695434808E-03        photos1/x700
                                               lg F_irr        2749    1.2698970004336019D+01
                   1                2749     2.216895849E+08     2.117207455E+00    -5.724097244E+00     6.778111874E-03        photos1/x700
                                               lg F_irr        2749    1.2698970004336019D+01
                   1                2749     2.216895849E+08     2.117207390E+00    -5.713849448E+00     6.739682584E-03        photos1/x700
                                               lg F_irr        2749    1.2698970004336019D+01
                star               model                 age                mass             lg_mdot           (r-rl)/rl          last photo
                   1                2750     2.216895867E+08     2.117204022E+00    -5.713849448E+00     6.714233214E-03        photos1/x700
                                               lg F_irr        2750    1.2698970004336019D+01
       2750   7.343212  1.214E+04   0.969037   0.969037   2.117204   2.117204   0.289812   0.011725   0.280162   0.601416  -3.186960    904    567
   0.241207   1.727500   0.311526 -36.435469  -0.214368  -5.713849   0.000000   0.690597   0.001600   0.020000   0.378468   0.058189      4     49
 2.2169E+08  17.035797   1.912756 -14.750470   2.908839  -9.265046   0.000000   0.000076   0.002099  1.959E-02  2.012E-02  0.630E-04  max increase

save photos1/x750 for model 2750

***********





Cheers,

Wen-Cong




在 2016年6月10日,下午6:20,Bill Wolf <wmwolf at physics.ucsb.edu> 写道:

> Hi Chen,
> 
> Full inlists and model files would be helpful for us to reproduce this error.
> 
> My guess, though, is that you’re having periodic retries and/or backups. Often I’ll have a series of steps where the tilmestep is increasing as much as it can (max_increase) and then once it has gotten to some critical timestep size, the hydro solver struggles to converge and forces a retry, bringing the timesteps back down.
> 
> If this is a case, I have no simple solution for you, but you might try to use the hydro debugging tools to see what the physical situation is that is causing problems. See my guide over here: http://wmwolf.github.io/projects/mesa_debugging/ for an example of how to do this.
> 
> Bill Wolf
> 
> ________________________ 
> William Wolf 
> wmwolf at physics.ucsb.edu 
> 
> UCSB Department of Physics 
> 
> On June 10, 2016 at 9:59:52 AM, Chen WC (chenwc at pku.edu.cn) wrote:
> 
>> Dear All,
>> I use the “binary_donor_only_implicit_mdot” to calculate the evolution of an intermediate-mass X-ray binary including a 1.4 Msun NS and a 3.0 Msun donor star (with a solar composition),
>> and the initial orbital period is 1.0 day.
>> 
>> In calculation, I consider the irradiation process by X-ray luminosity of the NS. In the list_project:
>> ****
>> accretion_powered_irradiation = .true.
>> col_depth_for_eps_extra = 100.0 ! g/cm^2
>> ****
>> Other input parameters such as accretor_radius_for_irrd, max_F_irr were adopted to be default.
>> 
>> 
>> However, after the mass of the donor star decrease to be ~ 2.3-2.4 Msun, the calculation is very difficult to continue.
>> The timestep decrease to be <0.1 yr, and the terminal output  for timestep control displays “max increase”.
>> 
>> I would appreciate if someone can help me to solve this problem.
>> 
>> 
>> Cheers,
>> 
>> 
>> Chen Wen-Cong
>> 
>> 
>> ------------------------------------------------------------------------------ 
>> What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic 
>> patterns at an interface-level. Reveals which users, apps, and protocols are 
>> consuming the most bandwidth. Provides multi-vendor support for NetFlow, 
>> J-Flow, sFlow and other flows. Make informed decisions using capacity 
>> planning reports. https://ad.doubleclick.net/ddm/clk/305295220;132659582;e_______________________________________________ 
>> mesa-users mailing list 
>> mesa-users at lists.sourceforge.net 
>> https://lists.sourceforge.net/lists/listinfo/mesa-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160611/c55dbd8d/attachment.html>


More information about the Mesa-users mailing list