[Mesa-users] dt limit for binary system

[Amedeo Romagnolo]

Hello Josiah,

thanks a lot for your reply.
Just a small question: but in this way one wouldn't be able to see either
when the CE dissipates or how the stellar companion behaves inside the CE
(accreting, etc...), right?

Cheers,

Amedeo

On Wed, 2 Sep 2020 at 02:17, Josiah Schwab <jwschwab at ucsc.edu> wrote:

> Hi Amedeo,
>
> is there no way one could take pre-CE simulations, simulate a simplified
>> version of CE on his own and then start MESA again once the CE is supposed
>> to dissipate?
>>
>
> This is certainly something one could do, but there is no built-in or
> "cookbook" way of doing it, so you'll have to do the hard work yourself.
>
> One possible procedure might look something like:
>
> A) Stop the binary run just before the CE and save the models.
>
> B) Take the model that will lose its envelope and put it in a single-star
> run.  Remove its outer layers rapidly using controls like
> http://mesa.sourceforge.net/star_job_defaults.html#remove_surface_at_cell_k
> or
> http://mesa.sourceforge.net/star_job_defaults.html#relax_mass
> Once things have settled down, save a model.
>
> C) Calculate what you think the post-CE orbital period should be, possibly
> using information from the models in at step A.
>
> D) Start a new binary run with one of the stellar models from step A and
> the stripped model from step B, and start them with the initial period you
> calculated in step C.
>
> You'll likely need to tailor that approach based on what sorts of
> scientific questions you're looking to answer.
>
> Hope that helps,
> Josiah
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200903/ce2c0d37/attachment.htm>