[Mesa-users] Double Roche lobe overflow

[Ken Shen]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20171211/3b8fdeab/attachment.html>