[Mesa-users] Doubt on overshoot parameters
Pedro Nuno David Gomes
pedro.david.gomes at ist.utl.pt
Sun Dec 9 08:04:34 EST 2018
Dear MESA users,
I have a doubt on the overshoot parameters I was hoping you could help
me clarify. I'm sure this is a very simple doubt, but I can't wrap my
head around it.
I have an inlist of a previous MESA release (release 7456) with the
following controls regarding the overshooting,
overshoot_f_below_nonburn = 0.018
overshoot_f_above_burn_h = 0.018
overshoot_f_above_burn_he = 0.018
From what I understood, using these controls one establishes that
overshoot occurs at a distance 0.018*H from the boundaries of convective
regions, with H being the scale height. Right?
Now, I was trying to run the same controls in the latest MESA release
but I got an error: "Fortran runtime error: Cannot match namelist object
name overshoot_f_below_nonburn".
So, the previous controls changed with posterior releases. I'm having
some doubts understanding how to implement these controls in the current
release. At first I simply thought about replacing them directly with
(what I think are) the updated controls,
overshoot_f_above_burn_h_core = 0.018
overshoot_f_above_burn_h_shell = 0.018
overshoot_f_below_nonburn_shell = 0.018
(I'm not considering helium burning regions because I'm only interested
in analysing the lower Red Giant Branch). But when I run MESA I get an
error: "ERROR: when using overshooting above, must set f0_above > 0".
From what I read on the controls defaults, namely
"The switch from convective mixing to overshooting happens at a
distance f0*Hp into the convection zone
(...) where Hp is the pressure scale height at that location.",
I also have to define f0. So, given this explanation, what I thought I
had to do was replace f with f0,
overshoot_f0_above_burn_h_core = 0.018
overshoot_f0_above_burn_h_shell = 0.018
overshoot_f0_below_nonburn_shell = 0.018
And in this case MESA runs. However, I'm not sure if this is the way -
from what I read after, namely
"for example, if you want overshooting of 0.2 scale heights beyond
the normal edge, you might want
to back up 0.05 scale heights to get the diffusion coeff from
near the edge and then go out by 0.25
scale heights from there to reach 0.2 Hp beyond the old boundary.
In the inlist this would mean setting
the “f0” to 0.05 and the “f” to 0.25."
I have to define in the inlist controls with f and f0 simultaneoulsy;
so, from my understanding, this would mean writting, for instance,
overshoot_f_above_burn_h_core = 0.018
overshoot_f_above_burn_h_shell = 0.018
overshoot_f_below_nonburn_shell = 0.018
overshoot_f0_above_burn_h_core = 0.001
overshoot_f0_above_burn_h_shell = 0.001
overshoot_f0_below_nonburn_shell = 0.001
Is any of these the right way, or am I missing something? If this is the
case, could you help me understand how to replace the previous controls
with these values?
Thank you for your attention. Regards,
Pedro
More information about the Mesa-users
mailing list