<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Jonathan,<div class=""><br class=""></div><div class="">What I think Tomek is saying (and what I normally do for the Gamma EOS) is define Abar so that it takes into account the contributions from _both_ the ions and the electrons. </div><div class=""><br class=""></div><div class="">So what he is saying is that you should set eos_singleSpeciesA to abar/(1+zbar), in your notation. </div><div class=""><br class=""></div><div class="">Best,</div><div class=""><br class=""></div><div class="">John</div><div class=""><br class=""></div><div class=""><blockquote type="cite" class=""><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-color: rgb(204, 204, 204); border-left-style: solid; padding-left: 1ex;"><div text="#000000" bgcolor="#FFFFFF" class="">On Sep 16, 2015, at 11:14 AM, Jonathan Slavin <<a href="mailto:jslavin@cfa.harvard.edu" class="">jslavin@cfa.harvard.edu</a>> wrote:</div></blockquote></div></div></blockquote><div><blockquote type="cite" class=""><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div dir="ltr" class=""><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Yes, that's exactly what I said. How does one get FLASH to use zbar? I just tried again after running setup with the added argument</div><div class="gmail_default" style="">-with-unit=physics/Eos/EosMain/Gamma<br class=""></div><div class="gmail_default" style="">but I still get that d*T/P * kboltz/amu_cgs = 1.28 (= abar) instead of 0.612 (= abar/(1 + zbar)) like I would expect.</div><div class="gmail_default" style=""><br class=""></div><div class="gmail_default" style="">Jon</div></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Wed, Sep 16, 2015 at 10:23 AM, Tomasz Plewa <span dir="ltr" class=""><<a href="mailto:tplewa@fsu.edu" target="_blank" class="">tplewa@fsu.edu</a>></span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF" class="">
Hi Jon -<br class="">
<br class="">
It is frequently assumed electron pressure is negligible compared to
ions, which is true for a nonionized gas. Obviously, the correct
formula should take into account contributions from both ions and
electrons. For a fully ionized gas, <br class="">
<br class="">
\mu = \frac{\sum_j n_j A_j }{\sum_j n_j ( 1 + Z_j) }<br class="">
<br class="">
Tomek<br class="">
--<div class=""><div class="h5"><br class="">
<div class="">On 09/16/15 10:00, Jonathan Slavin
wrote:<br class="">
</div>
<blockquote type="cite" class="">
<div dir="ltr" class="">
<div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Hi all,</div>
<div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br class="">
</div>
<div class="gmail_default" style="font-family:arial,helvetica,sans-serif">I'm doing
simulations aimed at understanding aspects of supernova
remnant evolution. Since the SNR is expanding the general
interstellar medium, the proper mean mass per atom (i.e.
neutrals and ions) is about 1.28 amu, given 90% H and 10% He
by number. As I understand it, this is set by the runtime
parameter eos_singleSpeciesA, which corresponds to abar. If I
want to assume constant, uniform, complete ionization then, as
I understand it, I need to set zbar via the runtime parameter
eos_singleSpeciesZ, to a value of ~1.091 = 1.2/1.1 = (n(H+) +
2*n(He++))/(n(H) + n(He)). At least that's what I thought
from the docs where at one point it's stated that n(e) =
zbar*n(ion). However, when I look at
density*temperature/pressure*kB/amu, I get abar. So it seems
that zbar is not used. What I expect is that</div>
<div class="gmail_default" style="font-family:arial,helvetica,sans-serif">P = (1 +
zbar)/abar*kB*density*temperature</div>
<div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Where P is the
total pressure (electron pressure + atom pressure). Is this
not true? Do I need to specify some extra unit for zbar to be
used (maybe -with-unit=physics/Eos/EosMain/Gamma)? Am I
misunderstanding the meaning of zbar? Any help would be
appreciated.</div>
<div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br class="">
</div>
<div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Thanks,</div>
<div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Jon</div>
<div class=""><br class="">
</div>
-- <br class="">
<div class="">
<div dir="ltr" class="">
<div class="">
<div dir="ltr" class="">________________________________________________________<br class="">
Jonathan D. Slavin Harvard-Smithsonian
CfA<br class="">
<a href="mailto:jslavin@cfa.harvard.edu" target="_blank" class="">jslavin@cfa.harvard.edu</a>
60 Garden Street, MS 83<br class="">
phone: <a href="tel:%28617%29%20496-7981" value="+16174967981" target="_blank" class="">(617) 496-7981</a> Cambridge, MA 02138-1516<br class="">
cell: <a href="tel:%28781%29%20363-0035" value="+17813630035" target="_blank" class="">(781) 363-0035</a> USA<br class="">
________________________________________________________<br class="">
<br class="">
</div>
</div>
</div>
</div>
</div>
</blockquote>
<br class="">
</div></div></div>
</blockquote></div><br class=""><br clear="all" class=""><div class=""><br class=""></div>-- <br class=""><div class="gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class="">________________________________________________________<br class="">Jonathan D. Slavin Harvard-Smithsonian CfA<br class=""><a href="mailto:jslavin@cfa.harvard.edu" target="_blank" class="">jslavin@cfa.harvard.edu</a> 60 Garden Street, MS 83<br class="">phone: (617) 496-7981 Cambridge, MA 02138-1516<br class="">cell: (781) 363-0035 USA<br class="">________________________________________________________<br class=""><br class=""></div></div></div></div>
</div>
</div></blockquote></div><br class=""></div></body></html>