[mesa-users] deuterium in solar calibration test and internal energy

[Aaron Dotter]

Hi Masanobu,

1) deuterium = h2 is included in add_pp_extras and therefore also in
pp_and_cno_extras.net.  Maybe you're just writing out the information to
the profile or history file?  You can add it to one or both, see the
*_columns.list in mesa/star/defaults.

2) C_V = (dE/dT)_V does not mean the same as E = C_v * T.  You would have
to integrate C_V * dT to get U.


I hope this helps.

Cheers,
Aaron



On Fri, Nov 21, 2014 at 6:59 PM, Masanobu Kunitomo <
kunitomo at geo.titech.ac.jp> wrote:

> Dear all,
>
> Today I would like to ask two questions: (1) whether it is possible to
> include deuterium in the solar calibration test and (2) the definition of
> the internal energy.
>
>
> (1)
> I am interested in  the evolution of 1Msun star from Pre-Main Sequence to
> MS phase. In my project, I use the input parameters such as mixing length
> alpha, f_overshoot, and initial composition (X, Y, Z) determined by solar
> calibration test (in mesa/star/test_suite/solar_calibration) with
> “create_pre_main_sequence_model".
>
>
> I found that deuterium is not included in the solar calibration test even
> if I set “pp_and_cno_extras.net”, “create_pre_main_sequence_model =
> .true.” and “set_uniform_initial_composition” with "initial_h2 = 2d-5". I
> got the only initial h1 fraction by the test. Since I am interested in the
> PMS evolution in which deuterium burning is important, in my project I add
> deuterium and subtract its value from h1 mass fraction determined by the
> test.
> However, this causes the disagreement with the solar values (e.g., R, L,
> Rcz etc.) at the solar age, although they are well reproduced in the solar
> calibration test. (In particular, Y_surf is obviously changed due to the
> addition of deuterium which is changed to 3He.) I guess this is because the
> composition is changed by the addition of deuterium and the enhancement of
> 3He.
>
> Therefore, I would like to include deuterium in the solar calibration
> test. Is it possible?
>
>
> (2)
> I think the specific internal energy, E, is given by eos module as
> described in Sec. 4.2 of paper I. I found that the internal energy disagree
> with the energy estimated by cvT (up to about 20%), where cv is the
> specific heat at constant volume. I thought that the specific internal
> energy should be the same value as cvT. Is this wrong?
>
>
> I appreciate any comments. Thanks in advance!
>
> Best wishes,
> Masanobu
>
> --
> Masanobu Kunitomo
> Ph.D. student / JSPS Research Fellow
> Tokyo Institute of Technology
> Dept. of Earth and Planetary Sciences
> Email: kunitomo at geo.titech.ac.jp
>
>
>
>
> ------------------------------------------------------------------------------
> Download BIRT iHub F-Type - The Free Enterprise-Grade BIRT Server
> from Actuate! Instantly Supercharge Your Business Reports and Dashboards
> with Interactivity, Sharing, Native Excel Exports, App Integration & more
> Get technology previously reserved for billion-dollar corporations, FREE
>
> http://pubads.g.doubleclick.net/gampad/clk?id=157005751&iu=/4140/ostg.clktrk
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20141121/6770d565/attachment.html>