[Mesa-users] entropy relaxation skipped

[Bradley Munson]

Hi Akaash,

It's hard to diagnose the problem when we do not have the necessary files to reproduce it. I think you are interpreting the output correctly (the entropy relaxation struggles for a bit and then ultimately gets "skipped" without taking a successful step).

Have you tried playing around with the following options? The defaults for both are 1d-9, but the following values typically work for me with similar looking entropy profiles:

timescale_for_relax_entropy = 1d-12
max_dt_for_relax_entropy = 5d-13

​Try adjusting those options and see if your model at least takes some steps in the right direction.

Also, what is the resolution of your input entropy file like? Could it be too sparse?

Best,

Brad
________________________________
From: Mesa-users <mesa-users-bounces at lists.mesastar.org> on behalf of Akaash Conhye via Mesa-users <mesa-users at lists.mesastar.org>
Sent: Wednesday, December 16, 2020 12:13 PM
To: Francis Timmes <fxt44 at mac.com>
Cc: mesa-users <mesa-users at lists.mesastar.org>
Subject: Re: [Mesa-users] entropy relaxation skipped

Hi Francis,

Removing this max_model_number line allows for proper time evolution to occur i.e. the white dwarf goes about its cooling. This line is something I found in the relax_composition_j_entropy template on the test suite. It is a means to stop the simulation right when relaxation finishes, at least as far as I know. The entropy relaxation process is just skipped over after a few retries and backups when the line is removed. This is precisely what I find quite strange.

Thanks,

Akaash
________________________________
From: Francis Timmes <fxt44 at mac.com>
Sent: 16 December 2020 14:01
To: Conhye, Akaash (UG - Physics) <ac01283 at surrey.ac.uk>
Cc: Frank Timmes <fxt44 at mac.com>; mesa-users <mesa-users at lists.mesastar.org>
Subject: Re: [Mesa-users] entropy relaxation skipped

the run is terminating as instructed:

> stop because model_number >= max_model_number


try setting max_model_number to a larger value in your inlist(s).

fxt




> On Dec 16, 2020, at 6:13 AM, Akaash Conhye via Mesa-users <mesa-users at lists.mesastar.org> wrote:
>
> Hi all,
>
> I'm trying to relax composition and entropy for a white dwarf model. The profiles come from hydrodynamic simulations with a Helmholtz EOS without coulomb corrections.
>
> composition relaxation occurs just fine. However, my problem is entropy relaxation. In the terminal output file I've attached, MESA tries a few times with entropy relaxation and gives up instead of the usual timescale reduction to the minimum timestep. Does anyone know why this might be occurring?
>
> Also attached is a comparison of the entropy profiles pre-relaxation (default) post relaxation (relax363) and the target profile (from hydro simulation)
>
> Thanks,
>
> Akaash
>
> <relax363_s_and_mu_no_initial_dt_max_dt_s_1d0_mlt_cox_mesh_coeff_0.5_timescale_default_eddington_t_tau_iterated.out><sprofile_default_relax363_hydro_loglog.png>_______________________________________________
> mesa-users at lists.mesastar.org
> https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.mesastar.org%2Fmailman%2Flistinfo%2Fmesa-users&data=04%7C01%7Cac01283%40surrey.ac.uk%7C726c5d348f2b4d3d31c208d8a1cb26c2%7C6b902693107440aa9e21d89446a2ebb5%7C0%7C0%7C637437241188199657%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=NfuXykXuLCuGFCq%2BEbFJT%2FCthWojEgoE3s8L12ThyJ0%3D&reserved=0<https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.mesastar.org%2Fmailman%2Flistinfo%2Fmesa-users&data=04%7C01%7Cbmunso2%40lsu.edu%7C89c769534502464b8b3608d8a1ee6321%7C2d4dad3f50ae47d983a09ae2b1f466f8%7C0%7C0%7C637437392525074864%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=lmnGm%2Fyeq%2B30V7%2FA2Uq1K0TZwS13C3z5%2FFEA%2B2VfqaI%3D&reserved=0>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20201216/b6de5a80/attachment.htm>