[mesa-users] Keplerian angular momentum for accreted material

[Robert Farmer]

Hi

How about trying the relax mass (on a non rotating model) then once you've
changed the mass then run the relax omega to spin the star up? ie split
into two separate steps.

In general, for the relax routines, if you get into convergence problems
try making the change per step smaller either by taking more steps or
having a smaller delta_ per step (depending on which relax options you are
using)

Rob

On Sun, Nov 13, 2016 at 11:01 AM, Héctor MR <hector.mr at pitt.edu> wrote:

>   Hola Pablo,
>
>   Sure! I will give you more details.
>
>   I am trying to implement accretion+rotation in a WD, for which I am
> relaxing the initial omega_div_omega_crit (plus mass_change).
>
>   If I use the control use_accreted_material_j = .true., I will have
> serious trouble to relax the initial mass change, which I have not managed
> to do this way. However, if I implement it in run_star_extras.f (via s%
> use_accreted_material_j = .true. in extras_finish_step plus my desired
> value for j), things seem to work. This, I do not understand.
>
>   If I do not relax the initial mass change, I will get "stop set_omega"
> plus the model number after relax_initial_omega_div_omega_crit, terminal
> output that I have located in adjust_mass.f90.
>
>   Finally, if I manage to relax both of them, my model will get stuck in
> small time steps <1 s ("failed in mixing angular momentum"), and then the
> usual "stopping because of convergence problems dt < min_timestep_limit".
>
>   I wish I could provide you with some plots, but I honestly do not know
> how to tackle debugging during the relaxation processes, and all the bugs I
> am reporting are difficult to diagnose for me. I am attaching the inlist I
> am using plus the model, the nuclear network, run_star_extras and the
> history+profile files hoping that you can find a solution (please, let me
> know if I forgot to include something). I would be quite grateful if you
> could lend me a helping hand with these issues.
>
>   Thanks!
>
>
>   --
>   H.
>
> ------------------------------------------------------------
> ------------------
> Developer Access Program for Intel Xeon Phi Processors
> Access to Intel Xeon Phi processor-based developer platforms.
> With one year of Intel Parallel Studio XE.
> Training and support from Colfax.
> Order your platform today. http://sdm.link/xeonphi
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20161116/72f1e59f/attachment.html>