[mesa-users] Full ionization

[Francis Timmes]

in addition to neglecting ionization, your P_gas neglects the effects 
of degeneracy (~ chemical potential / kT) and relativity (~ kT/m_e c^2) 
on the free electrons along with ignoring interactions between the 
components of the plasma (i.e., coulomb corrections).

examine the 1 and 10 msun main-sequence thermodynamic tracks in the rho-T plane.
armed with the information above, consider re-assessing the explanation of 
ionization causing all the delta behavior.

fxt



> On Oct 27, 2014, at 6:49 AM, Ehsan Moravveji <e.moravveji at gmail.com> wrote:
> 
> Hello MESA users,
> 
> On my continuing plan to prepare some classroom exercises based on MESA outputs, I came across an interesting feature that would like to share with you, and ask for your comments.
> 
> The idea is to check the assumption of P_ideal = P_gas + P_rad, where P_gas means pressure from fully-ionized ideal gas: P_gas=R_gas*Rho*T/mu.
> 
> Then, I calculate the above P_ideal inside two models: 1Msun and 10Msun. I also read the total pressure returned from EOS and stored in the output profile file as P_mesa.
> 
> Then, I would like to check for the difference between P_mesa and P_tot, so I define delta to be a relative difference between the two:
>      delta = (P_mesa - P_tot) / P_mesa * 100 [%]
> delta essentially tests how well the stellar interior can be considered as a fully ionized mixture of ideal gas + radiation.
> Now, please take a look at the three attached plots. 
> 
> <ans03-profile-1Msun.png>
> Here, for a 1Msun model, delta is large and also blows up close to the surface. To me, the reason is partial ionization at cooler temperatures, and therefore the logical departure from the assumption of the gas being ideal fully ionized.
> 
> The second figure is quite similar to the first but for a 10 Msun model:
> <ans03-profile-10Msun.png>
> Here, MESA is doing better, since at higher temperatures, the full ionization is relatively better fulfilled.
> 
> Then, i look at the above “delta” quantity that I defined, only in the core and on the surface of the 10 Msun model, and trace it during the main sequence. 
> The x-axis gives the drop in core hydrogen mass fraction.
> Here, I find an interesting feature in the core (bottom panel): delta is somehow cyclic! 
> 
> <ans03-hist-10Msun.png>
> How should I understand this?
> I set use_lnPgas = .true., so MESA is using the interpolated EOS tables. Is what I’m seeing representative of cyclic interpolation goodness?
> Sometimes delta is close to 0.0% and sometimes it goes up to ~0.30 %.
> 
> I use v.6794, and the inlist is attached. To avoid writing a long email, I cut out some details, so let me know if something is unclear.
> 
> Cheers
> Ehsan.
> 
> <inlist_exercise>
> 
> 
> 
> 
> ------------------------------------------------------------------------------
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users