<div dir="ltr">Oh I see, thank you for your suggestions! I will try firstly establish my own IONMIX files. <div>Thank you again!</div><div><div>Best </div></div><div>Siyi YU</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 21, 2017 at 12:52 PM, Carlo Graziani <span dir="ltr"><<a href="mailto:carlo@oddjob.uchicago.edu" target="_blank">carlo@oddjob.uchicago.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Klaus is, of course, correct. However, note that the ionization of Hydrogen<br>
is a straightforward analytic solution of the Saha equations, requiring no<br>
laborious computation. If you Google "Saha equation Hydrogen" you'll get<br>
260,000+ hits, many of the top of which are informative. Many thermal physics<br>
and stellar astrophysics texts also supply the solution.<br>
<br>
Which is all by way of saying, if you don't want to mess around with the<br>
internals of FLASH's EOS (and I wouldn't blame you if you didn't), it<br>
should not be too much trouble to make your own IONMIX-style table for<br>
Hydrogen, using the analytic solution.<br>
<br>
Carlo<div class="HOEnZb"><div class="h5"><br>
<br>
On 06/21/2017 02:29 PM, Klaus Weide wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
On Wed, 21 Jun 2017, Siyi Yu wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hey Klaus,<br>
I am also trying to implement hytrogen ionization in my simulation. I am<br>
wondering, instead of using IOMIX documents, is it possible to implement a<br>
self-written function that calculates electron densities from temperature<br>
please?<br>
</blockquote>
<br>
Hi Siyi,<br>
<br>
In principle this is possible. Even in the best case, it will probably<br>
require some significant (re)writing of code. How much, and how difficult,<br>
will depend a lot on what you want to do. Some initial questions to narrow<br>
down your requirements:<br>
<br>
* Do you want to use this in 1T or in 3T (or RadFLAH?) simulations?<br>
<br>
* Do you Nele / ZBAR depend on only (density,temperature,compositi<wbr>on),<br>
or on additional independent variables? (maybe including time /<br>
history?)<br>
<br>
* Do you just want one Eos implementation for all of your material(s),<br>
or do you want an implementation that models just the behavior of the hydrogen<br>
fraction in a multi-material, multi-type situation and can be combined<br>
with other implementations for other materials?<br>
<br>
* Do you want an Eos that is thermally ideal or something more<br>
complicated?<br>
<br>
* Do you need the derived quantities that Eos can potentially<br>
provide - see Eos.h - and if yes, which ones, or would you be satisfied<br>
(at least initially) with just getting the main variables<br>
like EOS_PRES, EOS_EINT, EOS_TEMP returned correctly by Eos?<br>
<br>
Klaus<br>
<br>
</blockquote>
<br>
<br></div></div><span class="HOEnZb"><font color="#888888">
-- <br>
Carlo Graziani <a href="tel:%28773%29%20702-7973" value="+17737027973" target="_blank">(773) 702-7973</a> (Voice)<br>
University of Chicago Astronomy <a href="tel:%28773%29%20702-6645" value="+17737026645" target="_blank">(773) 702-6645</a> (FAX)<br>
5640 South Ellis Avenue -----------------------------<wbr>--------<br>
ERC 549 | The less a statesman amounts to, the<br>
Chicago, IL 60637 | more he loves the flag.<br>
<a href="mailto:carlo@oddjob.uchicago.edu" target="_blank">carlo@oddjob.uchicago.edu</a> | --- Kin Hubbard<br>
<br>
<br>
</font></span></blockquote></div><br></div>