[Mesa-users] D_omega_flag, am_nu_rot_flag, and a possible memory leak

[agilkis at alumni.technion.ac.il]

Hi mesa-users,


I am looking into core rotation rates of massive stars and the influence of the Spruit-Tayler (ST) dynamo on them. Following a question about rotational mixing during the MESA Summer School, Frank suggested to me using D_omega_flag to get smoother profiles. This seemed to do the job at first, but further runs seemed to raise a couple of issues:

1) I noticed that when using D_omega_flag, the core rotation rate becomes faster. I think this is because for the ST dynamo mixing of chemical elements is calculated separately from the angular momentum transport, and using D_omega_flag (with use_D_omega_for_am_nu_rot set to its default value of .true.) results in an underestimation of the core-envelope coupling by magnetic fields (according to Spruit 2002, mixing of elements is less efficient than the angular momentum transport).

2) To fix the above, I tried to use am_nu_rot_flag as well as D_omega_flag. This gives a reasonable core rotation rate, similar to the case with no smoothing. However, I noticed that the memory usage in this case kept increasing throughout various runs, reaching up to 10GB in some cases. Without am_nu_rot_flag, the runs used between 1GB and 2GB.


Attached is a figure of the rotation profile at core helium depletion for no smoothing (inlist0), using D_omega_flag (inlist1), and using both D_omega_flag and am_nu_rot_flag (inlist2).


Is the memory leak a known issue? If it is, has it been fixed in a more advanced version of MESA? What would be recommended to do about the smoothing of the profile? I think the current options (version 10398) are (A) to not do smoothing; (B) smooth only with D_omega_flag and set use_D_omega_for_am_nu_rot to .false.


Thanks!

Avishai
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181212/d86c9c3d/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: M15V100deplHe4Omega.png
Type: image/png
Size: 27187 bytes
Desc: M15V100deplHe4Omega.png
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181212/d86c9c3d/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: M15V100HR.png
Type: image/png
Size: 34697 bytes
Desc: M15V100HR.png
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181212/d86c9c3d/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: M15V100vrot.png
Type: image/png
Size: 29218 bytes
Desc: M15V100vrot.png
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181212/d86c9c3d/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist0
Type: application/octet-stream
Size: 6412 bytes
Desc: inlist0
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181212/d86c9c3d/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist1
Type: application/octet-stream
Size: 6410 bytes
Desc: inlist1
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181212/d86c9c3d/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist2
Type: application/octet-stream
Size: 6408 bytes
Desc: inlist2
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181212/d86c9c3d/attachment-0002.obj>