[Mesa-users] Irradiating a Planet

[Evan Bauer]

Thanks Brian. The plot is very helpful. I think I agree with your interpretation, so I guess the ionization stuff from MESA 4 isn’t that relevant in this case. MESA doesn’t remove material from the model unless you tell it how, so this looks like a case where it’s just trying to expand to infinity because that’s the best it can do with unbound surface material without some prescription for mass loss.

Cheers,
Evan


> On Dec 2, 2019, at 4:07 PM, Brian Jackson <bjackson at boisestate.edu> wrote:
> 
> Evan,
> Thanks for responding.
> 
> What I am seeing does not seem to be the same as Figure 56 -- see attached
> screenshot. Figure 56 shows that both energies are negative when the
> mass coordinate q is equal to the total mass (i.e., at the top of the
> atmosphere). In my model, however, the specific energy is negative
> throughout much of the model, but then near the surface (near q = 1),
> the energy becomes positive. When that happens in the time evolution,
> the model slows way down and takes smaller and smaller time steps. I
> also notice that the model planet's radius begins to increase without
> any obvious limit, eventually exceeding the size of the Sun (which
> seems strange for a 29 Earth mass planet with about a 10% gaseous
> atmosphere).
> 
> So am I wrong in thinking the planet has become unbound?
> 
> In case it's helpful, I've attached the final inlist, pgstar inlist,
> and loaded model. Apologies if I missed anything.
> 
> Thanks again,
> Brian
> 
> 
> On Mon, Dec 2, 2019 at 11:09 AM Evan Bauer <ebauer at physics.ucsb.edu> wrote:
>> 
>> Hi Brian,
>> 
>> My guess is what you’re seeing is that the ionization and molecular dissociation energies are included in the definition of energy from the EOS. See section 8.4 and figure 56 of MESA 4. Does that seem consistent with what you see in your model?
>> 
>> Cheers,
>> Evan
>> 
>> 
>>> On Nov 29, 2019, at 3:41 PM, Brian Jackson <bjackson at boisestate.edu> wrote:
>>> 
>>> All,
>>> Happy Post-Thanksgiving.
>>> 
>>> If I can trouble you all once more, I've made some modifications to my
>>> calculations based on your input. However, when I track the
>>> total_energy parameter during the evolution, I find that, as I
>>> irradiate my planet, near the top of the atmosphere, the total_energy
>>> becomes positive. If I'm interpreting correctly, that means the
>>> planet's atmosphere is becoming unbound.
>>> 
>>> Could someone please confirm this interpretation or correct it?
>>> 
>>> Thanks,
>>> Brian
>>> 
>>> On Wed, Nov 20, 2019 at 7:38 PM Josiah Schwab via Mesa-users
>>> <mesa-users at lists.mesastar.org> wrote:
>>>> 
>>>> Hi Brian,
>>>> 
>>>>> After creating that model (attached as
>>>>> planet_setEntropy_29ME_9percent_10kB.mod), I then attempted to evolve
>>>>> that model for several Gyrs under irradiation (as indicated in the
>>>>> attached inlist). However, the calculation stalls out as the time-step
>>>>> gets smaller and smaller.
>>>> 
>>>> Thanks for including the inlists/model.
>>>> 
>>>>> Could you suggest any way to keep that calculation from stalling? I am
>>>>> using Mac OS 10.13.1, MESA 11701, and the most recent version of the
>>>>> SDK (August 30 2019).
>>>> 
>>>> In the terminal output, the dt_limit is "varcontrol".  That means that
>>>> the timestep is limited by zone-to-zone changes in temperature, density,
>>>> etc.  In the T-Rho diagram, things look pretty static. But if you add
>>>> 
>>>>   Profile_Panels1_win_flag = .true.
>>>>   Profile_Panels1_xmin = 8.78e-5
>>>> 
>>>> to your pgstar inlist, you'll be able to see heat from the irradiated
>>>> outer material moving deeper (in a Lagrangian sense) into the star.
>>>> 
>>>> You have to limit the x-axis to easily see this because this thermal
>>>> wave is quite thin in mass.  Following that requires short timesteps.
>>>> 
>>>> I just wanted to push back a bit on the statement that the calculation
>>>> is stalled.  Something is happening in the model and MESA is
>>>> timestepping at the timescale to follow that.  I recognize that this is
>>>> not the Gyr timescale that you're interested in.
>>>> 
>>>> 
>>>> I don't know what your science case needs, but I would probably try
>>>> starting from a higher entropy model.
>>>> 
>>>> Josiah
>>>> _______________________________________________
>>>> mesa-users at lists.mesastar.org
>>>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>>> 
>>> _______________________________________________
>>> mesa-users at lists.mesastar.org
>>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>> 
>> 
> <Screen Shot 2019-12-02 at 5.04.51 PM.png><inlist_pgstar><planet_irradiated_29ME_9percent_10kB_logflux_9.524.mod><inlist>_______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
>