[Mesa-users] [Mesa-developers] Pitfalls of eosPTEH?

[Pablo Marchant]

Think you want to turn off ELM and DT2 as well to restore the old EOS
completely.

Cheers

On Fri, Dec 28, 2018 at 4:51 PM Zhang, Zhaoxi M. (Michael) <
mz at astro.caltech.edu> wrote:

> Hi Bill,
>
>
> From controls.defaults, it seems that these values are the default as of
> now:
>
>
> use_eosPTEH_for_low_density = .true.
> use_eosPTEH_for_high_Z = .true.
> Z_for_all_PTEH = 0.040d0
> Z_for_any_PTEH = 0.039d0
> logRho_min_for_all_OPAL = -9.1d0
> logRho_min_for_any_OPAL = -9.2d0 ! -9.6d0
>
>
> Does setting use_eosPTEH_for_low_density and use_eosPTEH_for_high_Z to
> .false. restore the old (pre-10382) behavior, where I could disable PTEH
> completely?  It doesn't seem so, because my models are still crashing.  If
> not, how do I disable PTEH completely and get the old EOS back?
>
>
>
> Cheers,
>
>
> Michael
> ------------------------------
> *From:* Bill Paxton <paxton at kitp.ucsb.edu>
> *Sent:* Friday, December 28, 2018 1:01:53 PM
> *To:* Zhang, Zhaoxi M. (Michael)
> *Cc:* mesa-developers at lists.mesastar.org; mesa-users at lists.mesastar.org
> *Subject:* Re: [Mesa-users] Pitfalls of eosPTEH?
>
>
> On Dec 28, 2018, at 12:48 PM, Zhang, Zhaoxi M. (Michael) <
> mz at astro.caltech.edu> wrote:
>
> Hi Bill,
>
> I am reusing the old models.  They still work in version 10204 by setting
> use_eosPTEH_with_eosDT = .false. and use_PTEH_for_high_Z = .false. , but
> since version 10382, use_eosPTEH_with_eosDT is gone as an option and the
> old models no longer work.  Do you know how to restore the pre-10382
> behavior?
>
>
>          !### use_eosPTEH_for_low_density
>          !### use_eosPTEH_for_high_Z
>          !### Z_for_all_PTEH
>          !### Z_for_any_PTEH
>          !### logRho_min_for_all_OPAL
>          !### logRho_min_for_any_OPAL
>
>          ! overwrite the control of same name in the `Eos_General_Info`
> structure
>
>          ! if use_eosPTEH_for_low_density is true,
>          ! for logT < 3.7, use PTEH for logRho < about 7.5
>          ! for higher logT, use PTEH for logRho < logRho_min_for_any_OPAL
>          ! and do blend of OPAL and PTEH for logRho <
> logRho_min_for_any_OPAL.
>
>          ! if use_eosPTEH_for_high_Z is true.
>          ! use PTEH in place of HELM for Z above OPAL/SCVH range
>          ! Z >= Z_for_all_PTEH, then use pure PTEH with no blend
>          ! Z < Z_for_any_PTEH, then do not use PTEH
>          ! blend with PTEH for intermediate Z.
>          ! if use_eosDT2, blend with DT2.
>          ! else blend with OPAL/SCVH data.
>
>       use_eosPTEH_for_low_density = .true.
>       use_eosPTEH_for_high_Z = .true.
>       Z_for_all_PTEH = 0.040d0
>       Z_for_any_PTEH = 0.039d0
>       logRho_min_for_all_OPAL = -9.1d0
>       logRho_min_for_any_OPAL = -9.2d0 ! -9.6d0
>
>
> I can also try to generate new models, but the white dwarf models in
> star_data/white_dwarf_models have been deleted in the newest version, so I
> don't have an easy starting point.  Do you plan to regenerate them?
>
>
> no.
>
>
>
>
> Thanks for your help,
>
> Michael
> ------------------------------
> *From:* Bill Paxton <paxton at kitp.ucsb.edu>
> *Sent:* Friday, December 28, 2018 7:35:49 AM
> *To:* Zhang, Zhaoxi M. (Michael)
> *Cc:* mesa-developers at lists.mesastar.org; mesa-users at lists.mesastar.org
> *Subject:* Re: [Mesa-users] Pitfalls of eosPTEH?
>
> Are you reusing old models or building new ones?
>
> For the cases you mention, the differences between the older eos and PTEH
> could be large, so it might not be possible to load a model created with
> the one and expect it to work with the other.
>
> -b
>
>
>
>
> On Dec 28, 2018, at 1:16 AM, Zhang, Zhaoxi M. (Michael) <
> mz at astro.caltech.edu> wrote:
>
> Hello MESA users & developers,
>
> I'm trying to reproduce the simulations I did for my paper (
> https://arxiv.org/abs/1812.08793) with the latest MESA code, but they
> fail immediately with a lot of retries and backups.  I managed to determine
> that version number 10,204--specifically, the introduction of eosPTEH--is
> what broke my simulations.  Turning off PTEH by setting
> use_eosPTEH_with_eosDT = .false. and use_PTEH_for_high_Z = .false. fixes
> the problem.
>
> Is there a reason why PTEH isn't good to use in my scenario, namely iron
> enriched C/O white dwarfs with hot inflated envelopes?  Reading the latest
> release notes, it seems like low density and high metallicity environments
> are precisely what PTEH is meant for, so I have no idea why MESA can't even
> find a first model.
>
> Thanks for any help,
>
>
> Michael
> _______________________________________________
> 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/20181228/d9334818/attachment.html>