[Mesa-users] Dust formation
trobolo dinni
trobolo.trobolo.dinni5 at gmail.com
Mon Nov 19 20:19:10 EST 2018
Hi Aaron,
thank you.
Cheers,
Roberto
On Mon, 19 Nov 2018 at 20:32, Aaron Dotter <aaron.dotter at gmail.com> wrote:
> Hi Roberto,
>
> Are you aware of any other opacity tables dealing with low T and Rho and
>> including H?
>>
>
> No, I'm not. If anyone who reads this knows of a source, I hope they'll
> share it with us.
>
> Aaron
>
>>
>>
>> On Sat, 17 Nov 2018 at 03:59, Aaron Dotter <aaron.dotter at gmail.com>
>> wrote:
>>
>>> 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/20181120/0fd15c47/attachment.html>
More information about the Mesa-users
mailing list