[Mesa-users] Fwd: Rotation in MESA r22.11.1

Thu Nov 21 01:35:32 UTC 2024

Thanks, Pablo.

I looked at the inlist you suggested:

$MESA_DIR/star/test_suite/20M_z2m2_high_rotation/vim
inlist_to_end_core_he_burn
and copied the following into my inlist_ms in the &controls section,
leaving everything else the same as I had before (in particular,
omega/omega_crit = 0.6):

  am_nu_visc_factor = 0
  am_D_mix_factor = 0.0333333333333333d00
  D_DSI_factor = 0
  D_SH_factor = 1
  D_SSI_factor = 1
  D_ES_factor = 1
  D_GSF_factor = 1
  D_ST_factor = 1

  smooth_D_ST = 5
  smooth_nu_ST = 5
  smooth_D_omega = 5
  smooth_am_nu_rot = 5

However, it did not work. The track has a very strange shape, with no blue
hook. Additionally, I get a bunch of warning messages along the way, like:

...
        405   7.990808   3654.444   3.769356   3.769356   4.800000
3.625954   0.000000   0.008917   0.499213  -1.351599   1379      0
 3.2136E+00   7.990808   2.236979  -2.888566   2.602095 -99.000000
1.174046   0.986427   0.000169   0.013769   0.094427      6
 1.8059E+08   3.581852   3.679937  -4.504155 -99.000000  -8.185846
0.000000   0.000174   0.001460   0.013573  0.000E+00    varcontrol

        406   7.994394   3633.284   3.762294   3.762294   4.800000
3.625954   0.000000   0.008917   0.499237  -1.341105   1395      0
 3.1576E+00   7.994394   2.228715  -2.714662   2.595033 -99.000000
1.174046   0.986426   0.000169   0.013769   0.094351      6
 1.8059E+08   3.591556   3.653320  -4.339422 -99.000000  -8.160579
0.000000   0.000174   0.001460   0.013574  0.000E+00    varcontrol

 retry:  adjust_correction failed in eval_equations -- give up in solver
  407
 retry:  adjust_correction failed in eval_equations -- give up in solver
  407
 retry:  adjust_correction failed in eval_equations -- give up in solver
  407
 retry:  adjust_correction failed in eval_equations -- give up in solver
  407
 retry:  adjust_correction failed in eval_equations -- give up in solver
  407
 retry:  adjust_correction failed in eval_equations -- give up in solver
  407
 retry:  adjust_correction failed in eval_equations -- give up in solver
  407
 retry:  adjust_correction failed in eval_equations -- give up in solver
  407
 retry: max residual jumped -- give up in solver     407
 retry: max residual jumped -- give up in solver     407
 retry: max residual jumped -- give up in solver     407
 retry: max residual jumped -- give up in solver     407
 retry: max residual jumped -- give up in solver     407
 retry:  logT < hydro_mtx_min_allowed_logT          88     407
 retry:  logRho > hydro_mtx_max_allowed_logRho          62     407
 retry:  logRho > hydro_mtx_max_allowed_logRho          64     407
 retry:  logT < hydro_mtx_min_allowed_logT         162     407
        407   7.994394   3996.459   3.762291   3.762291   4.800000
3.625954   0.000000   0.008917   0.499237  -1.341105   1395     17
-2.0169E+00   7.994394   2.186812  -2.714674   2.595030 -99.000000
1.174046   0.986426   0.000169   0.013769   0.094351     22
 1.8059E+08   3.591555   3.735019  -4.339434 -99.000000  -8.316486
0.000000   0.000174   0.001460   0.013574  0.000E+00         retry

        408   7.994394   3974.201   3.762304   3.762304   4.800000
3.625954   0.000000   0.008917   0.499237  -1.341105   1372     17
-2.0169E+00   7.994394   2.188980  -2.714673   2.595043 -99.000000
1.174046   0.986426   0.000169   0.013769   0.094351      4
 1.8059E+08   3.591555   3.729653  -4.339433 -99.000000  -8.310856
0.000000   0.000174   0.001460   0.013574  0.000E+00          hold
...

from which MESA seems to recover, and then it ends normally. But the shape
is wrong. Evidently I'm still doing things wrong. I looked at the
documentation again, as well as the controls.defaults and star_job.defaults
files, but couldn't figure it out.

Willie

On Mon, Nov 18, 2024 at 6:24 AM Pablo Marchant <pamarca at gmail.com> wrote:

> Hi Guillermo,
>
> looking at your inlist, you do include rotation, which will give you the
> associated centrifugal effect, but you don't have something specified for
> rotational mixing, which I presume dominates the luminosity boost in the
> PARSEC tracks. By default no rotational mixing effects are included, so you
> need to turn them on. You can take a look at the section commented
> "rotation" in the inlist:
>
> $MESA_DIR/star/test_suite/20M_z2m2_high_rotation/vim
> inlist_to_end_core_he_burn
>
> as an example.
>
> Cheers
>
> On Sun, Nov 17, 2024 at 9:22 PM Torres, Guillermo via Mesa-users <
> mesa-users at lists.mesastar.org> wrote:
>
>> Dear all,
>>
>> I am new to MESA, and I'm trying to compute main-sequence models with
>> r22.11.1 for a range of degrees of rotation (omega/omega_crit), for a star
>> with 4.8 MSun. However, the tracks don't seem to change much as I increase
>> the rotation, compared to tracks from the PARSEC v2.0 series computed for
>> the same mass and composition (see plot attached). In particular, I get
>> very little change in luminosity. So I figure I must be doing something
>> wrong. I looked at previous threads on rotation, and tried several things I
>> saw there, but with no luck. I wonder if someone can tell me what I need to
>> change in my inlist to make this work. I attach my main-sequence inlist
>> (which loads a previously computed pre-main-sequence model for the same
>> mass and composition, and no rotation), as well as my PMS inlist.
>>
>> Thanks,
>> Willie
>> _______________________________________________
>> mesa-users at lists.mesastar.org
>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>
>>
>
> --
> Pablo Marchant Campos
> M.Sc on Astrophysics, Universidad Católica de Chile
> PhD on Astrophysics, Argelander-Institut für Astronomie, Universität Bonn
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20241120/c37e8165/attachment.htm>