[mesa-users] Turning on rotation near the ZAMS

[Richard Townsend]

Ah, it’s just occurred to me that v_surf is the suface *radial* velocity, not the surface rotational velocity. I shall go and don my dunce’s hat right away!

cheers,

Rich

> On Aug 11, 2015, at 1:28 AM, Richard Townsend <townsend at astro.wisc.edu> wrote:
> 
> Hmmm, I *can* reproduce your plot! 
> 
> So, the model itself is fine — but the v_surf value written to the history file appears to be bogus.
> 
> Bill, any insights?
> 
> cheers,
> 
> Rich
> 
>> On Aug 11, 2015, at 1:09 AM, Robert Farmer <rjfarmer at asu.edu <mailto:rjfarmer at asu.edu>> wrote:
>> 
>> So this is my final plot of omega_div_omega_crit and v_rot for the inlist provided, a few steps before the end, and mine seems fine. The outer layer is >100km/s and ~0.25. The omega_div_omega_crit=0.1 also looks similar. 
>> 
>> 
>> Rob 
>> 
>> On Mon, Aug 10, 2015 at 4:44 PM, Richard Townsend <townsend at astro.wisc.edu <mailto:townsend at astro.wisc.edu>> wrote:
>> On the order of 30-50 km/s, I imagine. And in the inlist I attached, I use omega_crit = 0.25, which should be in excess of 100 km/s. Similar result….
>> 
>> cheers,
>> 
>> Rich
>> 
>>> On Aug 11, 2015, at 12:43 AM, Robert Farmer <rjfarmer at asu.edu <mailto:rjfarmer at asu.edu>> wrote:
>>> 
>>> What would 0.1 Omega/Omega_crit in km/s be for a your star? I.e how far off are we?
>>> 
>>> Rob
>>> 
>>> On Mon, Aug 10, 2015 at 4:34 PM, Richard Townsend <townsend at astro.wisc.edu <mailto:townsend at astro.wisc.edu>> wrote:
>>> Hi folks —
>>> 
>>> I’m trying to turn rotation on near the ZAMS, by adding the following flags to the &star section of my inlist:
>>> 
>>>       near_zams_relax_omega_div_omega_crit = .true.
>>>       new_omega_div_omega_crit = 0.1
>>> 
>>> I then ensure the star evolves as a solid body by adding the following to the &controls section:
>>> 
>>>       set_uniform_am_nu_non_rot = .true.
>>> 
>>> Unfortunately, I’m not getting the results I expect — during the ZAMS/post-ZAMS evolution, the star remains very slowly rotating, with v_surf_km_s not rising above a km/s or so.
>>> 
>>> I’m clearly doing something wrong — can anyone tell me what? Input files attached.
>>> 
>>> Many thanks,
>>> 
>>> Rich
>>> 
>>> 
>>> 
>>> 
>>> 
>>> 
>>> ------------------------------------------------------------------------------
>>> 
>>> _______________________________________________
>>> mesa-users mailing list
>>> mesa-users at lists.sourceforge.net <mailto:mesa-users at lists.sourceforge.net>
>>> https://lists.sourceforge.net/lists/listinfo/mesa-users <https://lists.sourceforge.net/lists/listinfo/mesa-users>
>> 
>> <profile_panels1_001125.png>
> 
> ------------------------------------------------------------------------------
> _______________________________________________
> 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/20150811/17807a61/attachment.html>