[mesa-users] Problem with stellar winds during mass transfer

Pablo Marchant pamarca at gmail.com
Wed Mar 6 16:25:46 EST 2013


Cole, find attached a working directory with the binary routines i'm using
right now. Important details are that the inlist_binary_rlo must specify
"doing_binary_rlo = .true." and each star inlist must set
"use_other_adjust_mdot = .true.". Though I'm surprised about what you say,
I think in principle any added material should add no angular momentum in
the current version. You should explicitly check for the total angular
momentum of the system, you should be losing a small amount due to this (as
far as I know, in most cases orbital angular momentum will be far greater
than the spin angular momentum of each star).

I'm planning on getting to implement rotation processes in binaries soon,
but I have first another couple of non-rotating things to check on my
models.

Cheers!


On Wed, Mar 6, 2013 at 10:18 PM, Cole Adams <adams164 at gmail.com> wrote:

> I was not aware that the binary rotation routines were not implemented,
> since what I have been running appears to account for angular momentum
> transfer when both stars have rotation enabled. I was able to use Bill's
> workaround by turning off velocity (not surface velocity though) and it
> seemed to have worked while still maintaining the rotational physics.
>
> I managed a workaround for my current problem by re-arranging the code in
> winds.f but it still doesn't work completely. I would be glad to test some
> of your code for you.
>
> Cole
>
>
> On Wed, Mar 6, 2013 at 2:36 PM, Pablo Marchant <pamarca at gmail.com> wrote:
>
>> Cole, it works exactly as you say. I'be been working on modifications to
>> the binary routines and my current sources properly separate the two things
>> (and add some missing physics), I could share them with you if you're
>> interested in doing some testing.
>>
>> Regarding the topic of rotation though, there is no support for that in
>> the current version of the binary routines. Of course, you can have your
>> two stars rotating, but they won't transfer angular momentum between them.
>> In this sense I really don't understand what you say about spinning up the
>> accretor, have you implemented this yourself?
>>
>> Also, from what you say I guess you managed to go past the problem you
>> previously mentioned on this mailing list. Is that so? I haven't had much
>> more time to check that.
>>
>> Cheers!
>>
>>
>> On Wed, Mar 6, 2013 at 9:25 PM, Cole Adams <adams164 at gmail.com> wrote:
>>
>>> Hi,
>>>
>>> I've been running binary systems with rotation, and I've noticed a
>>> problem with the way that the winds and mass transfer function. During mass
>>> transfer, the wind and accretion components of mdot appear to not be
>>> treated separately. This leads to a rotational boost of mass loss not
>>> behaving well during mass transfer, since the accretion effectively
>>> overrides any wind loss effects. This then leads to the secondary star
>>> being spun up to speeds that are not reasonable (omega/omega_crit > 1)
>>> since the secondary star can't eject any mass or angular momentum. Looking
>>> at the code in winds.f it seems that if accretion is greater than the wind
>>> loss then wind loss is completely ignored. Is there a way to work around
>>> this?
>>>
>>> Cole
>>>
>>>
>>> ------------------------------------------------------------------------------
>>> Symantec Endpoint Protection 12 positioned as A LEADER in The Forrester
>>> Wave(TM): Endpoint Security, Q1 2013 and "remains a good choice" in the
>>> endpoint security space. For insight on selecting the right partner to
>>> tackle endpoint security challenges, access the full report.
>>> http://p.sf.net/sfu/symantec-dev2dev
>>> _______________________________________________
>>> 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
>>
>
>


-- 
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/20130306/1bd2499d/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: binary.tar.gz
Type: application/x-gzip
Size: 4037038 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20130306/1bd2499d/attachment.bin>


More information about the Mesa-users mailing list