mesa-users] dlnP/dlnMu ?

[Bill Paxton]

Hi,

Since mesa/eos returns both mu and free_e = 1/mu_e, and their  
partials wrt T and rho, is it possible to remove the full ionization  
assumption?

--Bill


BTW: I've copied this to mesa-users just to get in the habit of  
leaving a trace of such conservations that will undoubtedly come up  
again in the future.  Please humor me on this and include mesa-users  
in your messages too.


On May 4, 2008, at 11:02 AM, Frank Timmes wrote:

> hi ken,
>
> ok. assuming full ionization,
>
> mu_i = 1 /(sum Y_i) = Abar
>
> mu_e = 1 /(sum Y_i Z_i) = 1/Y_e
>
> mu = ( 1/mu_i + 1/mu_e )**(-1) = 1/(sum Y_i (Z_i + 1)) = Abar/(Zbar  
> + 1)
>
> d(mu) = dAbar/(Zbar + 1) + Abar/(Zbar + 1)**2 d(Zbar)
>
> the routines return derivatives with respect
> to composition like dP/d(Abar) and dP/d(Zbar).
>
> given the relation between mu, Abar and Zbar
> can you fill in the rest?
>
> fxt
>
>
>
>> I mean the mean mass per particle, like, 1/mu = 1/mu_i + 1/mu_e.
>
>
>