<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-GB" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">Hello MESA users,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I have come across a bit of a conundrum. I am currently trying to replicate a test that a collaborator of mine used with their own code.<o:p></o:p></p>
<p class="MsoNormal">I am modifying the temperature gradients in overshooting regions, comparing it’s inclusion in a solar-like star of 1Msun to a model of the same star with no modification.<o:p></o:p></p>
<p class="MsoNormal">Both stars are evolved from a zero age main-sequence model. The aim is to compare the profiles of these two models once they terminate at a set age, modifying the mixing_length_alpha of the modified model so that it reaches solar structure
at the terminal age. <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The difficulty is that the changes in the structural parameters I expect to see are very small. In order for the linear interpolation to be accurate I have been advised to get the compared mesh points to have a relative radial difference
of 1x10-5 Rsun. My concern is that since MESA is lagrangian in nature I will be unable to achieve this degree of relative accuracy. I was wondering if any of you smart people know anyway to control the mesh so that this kind of radial comparison is possible?
Or if it’s just realistically unfeasible? Please feel free to correct me if I have made any horrendous assumptions.
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">My main models have been run with mesh_delta_coeff = 1d-1 but I’ve also run models with 5d-1 and 1d0 to decrease simulation run time. mesh_min_dlnR = 1d-14.<o:p></o:p></p>
<p class="MsoNormal">All other mesh grid controls are default. I am using Gold tolerances with:<o:p></o:p></p>
<p class="MsoNormal">gold_tol_residual_norm1 = 1d-11<o:p></o:p></p>
<p class="MsoNormal">gold_tol_max_residual1 = 1d-9<o:p></o:p></p>
<p class="MsoNormal">gold_iter_for_resid_tol2 = 10<o:p></o:p></p>
<p class="MsoNormal">gold_tol_residual_norm2 = 1d-10<o:p></o:p></p>
<p class="MsoNormal">gold_tol_max_residual2 = 1d-8<o:p></o:p></p>
<p class="MsoNormal">gold_iter_for_resid_tol3 = 15<o:p></o:p></p>
<p class="MsoNormal">gold_tol_residual_norm3 = 1d-6<o:p></o:p></p>
<p class="MsoNormal">gold_tol_max_residual3 = 1d-4<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">For my 1d-1 mesh_delta_coeff model the average radial shift of a cell in the modified model compared to its nearest cell in the base model is 5.361718481731894e-05 rsun, while the relative radial shift is 2.2747013522397857e-04 rsun.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I am using MacOS Catalina 0.15.7 (19H524) , MESA Star version 12115, and Mesa SDK x86_64-osx-macos-19.11.2.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Cheers,<o:p></o:p></p>
<p class="MsoNormal">Josh<o:p></o:p></p>
<div style="mso-element:para-border-div;border:none;border-bottom:solid windowtext 1.5pt;padding:0cm 0cm 1.0pt 0cm">
<p class="MsoNormal" style="border:none;padding:0cm"><o:p> </o:p></p>
</div>
<p class="MsoNormal"><span style="color:black">Joshua Clarke</span><span style="font-size:12.0pt;color:black"><o:p></o:p></span></p>
<p class="MsoNormal" style="caret-color: rgb(0, 0, 0)"><span style="color:black">PhD Student</span><span style="font-size:12.0pt;color:black"><o:p></o:p></span></p>
<p class="MsoNormal" style="caret-color: rgb(0, 0, 0)"><span style="color:black">Department of Physics & Astronomy, University of Exeter,</span><span style="font-size:12.0pt;color:black"><o:p></o:p></span></p>
<div style="mso-element:para-border-div;border:none;border-bottom:solid windowtext 1.5pt;padding:0cm 0cm 1.0pt 0cm">
<p class="MsoNormal" style="border:none;padding:0cm"><span style="color:black">Stocker Road, Exeter, EX4 4QL, UK</span><span style="font-size:12.0pt;color:black"><o:p></o:p></span></p>
</div>
<p class="MsoNormal" style="caret-color: rgb(0, 0, 0)"><span style="color:black"> </span><span style="font-size:12.0pt;color:black"><o:p></o:p></span></p>
<p class="MsoNormal" style="caret-color: rgb(0, 0, 0)"><span style="color:black"> </span><span style="font-size:12.0pt;color:black"><o:p></o:p></span></p>
<p class="MsoNormal" style="caret-color: rgb(0, 0, 0)"><span style="color:black"> </span><span style="font-size:12.0pt;color:black"><o:p></o:p></span></p>
<p class="MsoNormal" style="caret-color: rgb(0, 0, 0)"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>