[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