[Mesa-users] Pitfalls of eosPTEH?

[Zhang, Zhaoxi M. (Michael)]

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<mailto: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<mailto: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<mailto:mesa-developers at lists.mesastar.org>; mesa-users at lists.mesastar.org<mailto: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<mailto: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<mailto:mesa-users at lists.mesastar.org>
https://lists.mesastar.org/mailman/listinfo/mesa-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181228/e08bbdc1/attachment.html>