[Mesa-users] Dust formation

[Aaron Dotter]

Hi Roberto,

I am dealing with gas having temperatures lower than 10^3 K, in some
> situations lower than 10^2 K and densities of the order of 10^-13gcm-3.
> According to the MESA paper this should below the thresholds tabulated. In
> this case how does the code handle the opacities?
>

The code decides what to do in cases like this from the logical flag
clip_to_kap_table_boundaries, which is set to .true. in
kap/public/kap_def.f90.

If you have a temperature that is below the table limit, then you'll get an
opacity value that is consistent with the lowest temperature in the table.
Likewise for density.

If you were to set this flag to false (and recompile!), then whenever you
go outside of table the call to kap would return an error and the run would
stop.

Depending on your case, you might look at using the Freedman et al.
opacities. These have a much lower temperature limit (logT_min = 1.88) but
also note that these are H-free.  This is discussed in MESA 2.

Cheers,
Aaron
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181116/d823c695/attachment.html>