MESAHub / MESAHub/mesa

High values for the thermodynamic consistency metric dsp in blending involving FreeEOS, and as well in OPAL/SCVH tables

Open
#822 0 comments 0 reactions 2 assignees Claimed by @Debraheem View on GitHub
eos
Dominant language
Fortran
Stars
249
Forks
82
Avg merge
1d 3h
Merged PRs (30d)
3

Description

Hi,

Context: In my scientific problem, I am evolving the envelope of very massive stars as they expand due to the injection of additional heat. During the course of this evolution, as the outer part of the envelope expands and cools, the model began closing the blending zone at the lower temperature limit of FreeEOS. This led to abnormal behavior in the reported values of temperature and density.

Upon closer inspection, I found that this behavior correlates with unusually high values of the thermodynamic consistency metric `dsp`, which is defined as:

$$
\text{dsp} = -\rho \left( \frac{ds}{d\ln \rho} \right) \bigg/ \left( \frac{dp}{dT} \right) - 1
$$

This metric quantifies the consistency of the thermodynamic derivatives. Large deviations from zero may indicate that the EOS is no longer self-consistent in that regime.




My `eos` controls for this run are:
```fortran
use_simple_Skye_blends = .true.

logT_min_for_any_Skye = 7.9
logT_min_for_all_Skye = 8.0

!logQ_max_FreeEOS_hi = 5.6d0
!logQ_max_FreeEOS_lo = 4.0d0


logT_max_FreeEOS_hi = 8.2d0 ! FreeEOS tabulations go up to here
logT_max_FreeEOS_lo = 8.1d0

logT_min_FreeEOS_hi = 3.0d0
logT_min_FreeEOS_lo = 2.7d0
logQ_min_FreeEOS_hi = -8.0d0
logQ_min_FreeEOS_lo = -10.0d0

logT_low_all_HELM = 2.2d0 ! HELM for lgT <= this
logT_low_all_SCVH = 2.3d0 ! SCVH for lgT >= this
```

For the default EOS blending limits one still find inconsistent values for dsp in similar regions:




Similarly if one sets the lower T blending for freeEOS higher as:
```fortran
logT_min_FreeEOS_hi = 3.6d0
logT_min_FreeEOS_lo = 3.5d0

```




Or if one removes the OPAL/SCVH tables:
```fortran
use_OPAL_SCVH = .false.
```




Given this, I am unable to set arrangement of the EOS that won't produce any weird behaviour while crossing that region of low T and low Rho. I'm opening this issue to better understand why `dsp` becomes large in this part of the envelope and whether this behavior is expected at the edge of FreeEOS's valid temperature range.

Contributor guide

Open the contributing guide

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.