[mesa-users] Rotation : relax_initial_omega_div_omega_crit + preserving omega in several inlists

[Bill Paxton]

Hi Hector,

Let me take a look.  I need your inlist and starting model.  

Bill



On Jun 16, 2016, at 8:48 PM, Héctor MR wrote:

>   Dear all,
> 
>   I have a question regarding the &star_job command relax_initial_omega_div_omega_crit. [Disclosure: this is the first time I am implementing rotation in MESA]
> 
>   I am analyzing rotation on an accreting, 1Msun-WD with radius R=8.25E-3*Rsun, so that the break-up frequency sqrt(G*M/R**3) is ~0.83; however, when I use new_omega_div_omega_crit = 0.1, I get
> 
> 
>   num_steps_to_relax_rotation         100
>         relax to omega: wanted-current, current, wanted           1    1.7485786361779884D-05    0.0000000000000000D+00    1.7485786361779884D-05
>         relax to omega: wanted-current, current, wanted           2    1.7330581050599486D-05    1.7485788969702041D-07    1.7505438940296507D-05
>         relax to omega: wanted-current, current, wanted           3    1.7177934014914814D-05    3.5010756599501660D-07    1.7528041580909830D-05
>         relax to omega: wanted-current, current, wanted           4    1.7016693789573586D-05    5.2583873499324927D-07    1.7542532524566834D-05
>         relax to omega: wanted-current, current, wanted           5    1.6850228116883350D-05    7.0169673253464221D-07    1.7551924849417992D-05
>         relax to omega: wanted-current, current, wanted           6    1.6680900334687578D-05    8.7758847614349309D-07    1.7558488810831069D-05
>         relax to omega: wanted-current, current, wanted           7    1.6513051889555891D-05    1.0534860275206089D-06    1.7566537917076500D-05
> 
> 
>   Rather than, let's say, 0.083, as I was expecting. I have done some "research" with "grep -r" in $MESA_DIR/star/private:
> 
>   /relax.f90: integer function relax_omega_check_model
> 
>   /do_one_utils: integer function do_check_limits
> 
>   /star_utils.f90: real(dp) function omega_crit(s, k), subroutine set_surf_avg_rotation_info(s)
> 
>   /hydro_vars.f90: subroutine set_rmid_and_Amid        !definition of rmid
> 
> 
>   And have not found anything weird. Is there anyone who can tell me what I am either doing wrong, obviating or misuderstanding? This sums everything up:
> 
>       new_rotation_flag = .true.
>       change_initial_rotation_flag = .true.
> 
>       !set_initial_omega_div_omega_crit = .true.
>       new_omega_div_omega_crit = 0.1
> 
>       relax_initial_omega_div_omega_crit = .true.
>       num_steps_to_relax_rotation = 100
> 
> 
>   Finally, my inlists are split in two pieces. How can I preserve the last value of omega_div_omega_crit in the first one for the initial angular velocity in the second one? Via x_ctrl(1)? And if so, how?
> 
> 
>   Thank you so much!
> 
> 
>   --
>   Héctor
> ------------------------------------------------------------------------------
> What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic
> patterns at an interface-level. Reveals which users, apps, and protocols are 
> consuming the most bandwidth. Provides multi-vendor support for NetFlow, 
> J-Flow, sFlow and other flows. Make informed decisions using capacity planning
> reports. http://sdm.link/zohomanageengine_______________________________________________
> 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/20160622/ed5e7796/attachment.html>