[mesa-users] Rotation : relax_initial_omega_div_omega_crit + preserving omega in several inlists
hector.mr at pitt.edu
Thu Jun 16 23:48:31 EDT 2016
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
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
/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
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!
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Mesa-users