[mesa-users] terminate because of overflowing initial model

[Pablo Marchant]

Ying, I happen to have the implementation for that wind prescription, as it
was used for some of the models in the last MESA paper (paxton et al. 2015
in case you haven't seen it) that reproduce the results of Wellstein &
Langer 1999. You can download the required inlists and run_star_extras.f
from

http://mesastar.org/results/mesa-iii-instrument-paper-inlists

The files you need are in the folder section2_binaries/6_massive_binary of
the compressed file. Just copy the inlists into your work directory and the
run_star_extras.f into the src folder of your work directory, and do a
./clean and a ./mk to recompile. Please go through the run_star_extras.f to
understand what is being done, the code essentially switches to the
prescription you want when X<0.4 at the surface, and uses the wind scheme
defined in the inlist otherwise.

And regarding the error, it's just what Rob says, your system is too
compact and one model is above its Roche lobe radius. You can remove the
terminating condition by setting

terminate_if_initial_overflow = .false.

in the &binary_controls section of inlist_project, but this will likely
result in convergence problems as you would obtain a very large mass
transfer rate, and even if it works, results might not make much sense. So,
I don't recommend you do this, just start with a longer orbital period.

On Sun, Aug 23, 2015 at 11:45 PM, Robert Farmer <rjfarmer at asu.edu> wrote:

> Hi
> You get this error becuase the first time step your star1 will be roche
> lobe overflowing. Try increasing the intial period such that the first step
> wont transfer mass.
>
> Rob
>
> On Sun, Aug 23, 2015 at 2:29 PM, Ying Qin <Ying.Qin at unige.ch> wrote:
>
>> Hi mesa users,
>>
>> I want to computer the evolution of close Wolf-Rayet binaries(The WR
>> companion is approximated as a point mass), including tidal angular
>> momentum exchange, different rotation of the Wolf-Rayet star, internal
>> magnetic fields, stellar wind mass loss, and mass transfer.
>>
>> 1. About the tidal angular momentum exchange, I add the following lines
>> in inlist_project
>>    do_tidal_sync = .true.
>>    do_initial_orbit_sync_1 = .true.
>>
>>
>> 2. I also added the rotation and magnetic fields in inlist1, please see
>> them in  inlist1.
>>
>>
>> 3. The mass loss of WR star is calculated according to Hamann et al. 1995
>> and and Hanann et al.1982.
>>
>>           log(Mdot/Msun yr^-1) =  eta * (-11.95 + 1.5logL/Lsun ),  where
>> logL/L_sun > 4.5(Hamann et al. 1995)
>>
>>           log(Mdot/Msun yr^-1) =  eta * (-35.8 + 6.8logL/Lsun ),  where
>> logL/L_sun < 4.5(Hamann et al. 1982)
>>
>> I know other winds are used in mesa,  but can I bother you to tell me how
>> to use this wind mass loss.
>>
>> 4. About mass transfer, I chose mdot_scheme = “Kolb” in inlist_project.
>>
>> I tried to use 'Dutch wind’ , not the wind I listed above,  to run mesa,
>> but it stopped with error:
>> termination code: Terminate because of overflowing initial model
>>
>>
>> In the attachment are my inlist1 and inlist_project, any comments are
>> kindly welcome, and thank you in advance!
>> If possible, could you please explain them clearly, because I am learning
>> to use mesa.
>>
>> Ying
>>
>>
>>
>> ------------------------------------------------------------------------------
>>
>> _______________________________________________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>>
>>
>
>
> ------------------------------------------------------------------------------
>
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>


-- 
Pablo Marchant Campos
M.Sc on Astrophysics, Universidad Católica de Chile
PhD student, Argelander-Institut für Astronomie
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150823/b9cb84f1/attachment.html>