[Mesa-users] Double Roche lobe overflow

[Ken Shen]

Thanks all!  Jared, I hadn't thought about just hard-wiring the binary_id
in the binary_ptr call inside the stars' other_hooks.  That should work!

Pablo (& Tassos), that scheme looks similar to what I'm envisioning, but
somewhat different; I'm thinking about small nova-like envelopes getting
ejected from the surface, so they shouldn't lead to the merger of the
system once they go past L2.


Thanks again,
Ken







On Mon, Dec 11, 2017 at 2:04 PM, Pablo Marchant <pamarca at gmail.com> wrote:

> Ninja'd by Tassos it seems
>
> On Mon, Dec 11, 2017 at 4:03 PM, Pablo Marchant <pamarca at gmail.com> wrote:
>
>> Also, there is an implementation of a simple contact prescription I
>> developed, which is explained in
>>
>> http://adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=201
>> 6A%26A...588A..50M&db_key=AST&link_type=ABSTRACT&high=51ac7f347015896
>>
>> you can use this by setting mdot_scheme = "contact" in your
>> inlist_project. Not quite what you say, but I think under certain
>> conditions its a more physically correct choice.
>>
>> Cheers
>>
>> On Mon, Dec 11, 2017 at 3:55 PM, Jared Brooks via Mesa-users <
>> mesa-users at lists.mesastar.org> wrote:
>>
>>> Hey Ken,
>>>
>>> Check out my other_wind_hook on mesastar.org (
>>> http://cococubed.asu.edu/mesa_market/add-ons.html). It sounds like it
>>> should fit your needs. I tried my best to document it well, so hopefully it
>>> should make sense without further explanation here. Please let me know if
>>> it is confusing and/or needs more explanation.
>>>
>>> -Jared
>>>
>>> On Mon, Dec 11, 2017 at 1:50 PM Ken Shen <kenshen at astro.berkeley.edu>
>>> wrote:
>>>
>>>> Hi all, hope I'm not duplicating a previous question.  I'm interested
>>>> in using MESA's binary capabilities to model systems where both components
>>>> might overflow their Roche lobes at the same time.  The basic prescription
>>>> would be to evaluate both stars' overflow values ("rl_relative_gap" in
>>>> binary_data), calculate Mdots as necessary via Ritter or Kolb, and transfer
>>>> from one star to the other or from both stars out of the system, depending
>>>> on which one(s) are undergoing RLOF.  Mass leaving the system would carry
>>>> the specific angular momentum of the binary.
>>>>
>>>> I've looked into the various other_ binary hooks, but none seems to be
>>>> quite right, especially for changing the nominal accretor's Mdot as it
>>>> expands to fill its Roche lobe.  For single stars, there is an
>>>> other_adjust_mdot, which would work except it doesn't know about the Roche
>>>> radius, and there is no binary equivalent.  Maybe there's a way to pass in
>>>> binary data into other_adjust_mdot?
>>>>
>>>> Any tips?
>>>>
>>>>
>>>> Thanks,
>>>> Ken
>>>>
>>>>
>>
>>
>> --
>> Pablo Marchant Campos
>> M.Sc on Astrophysics, Universidad Católica de Chile
>> PhD on Astrophysics, Argelander-Institut für Astronomie, Universität Bonn
>>
>
>
>
> --
> Pablo Marchant Campos
> M.Sc on Astrophysics, Universidad Católica de Chile
> PhD on Astrophysics, Argelander-Institut für Astronomie, Universität Bonn
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20171211/4873ca5d/attachment.html>