<div dir="ltr"><div class="gmail_extra"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">​Hi Ernesto,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">I don't specify the runtime parameter addThermalFlux.  Instead I set diff_useEleCond, useConductivity, useDiffuse, and useDiffuseTherm to .true. and specify </div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">​<span style="font-family:arial,helvetica,sans-serif">REQUIRES physics/materialProperties/Conductivity/ConductivityMain/PowerLaw</span></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">REQUESTS physics/Diffuse/DiffuseMain/Unsplit</div><div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">​Should I change one of those?​</div></div><div><br></div><div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">​I note that in the log file that</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">hy_fluxrespresentation = fluxes</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">I'm not sure exactly what that means.  My goal here is to verfy that ​thermal conduction is operating as I expect by comparing evaporation of a spherical cloud in 2D cylindrical symmetry with analytical results derived for steady flow including saturation of thermal conduction (Dalton & Balbus 1993).</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Regards,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Jon</div></div><br><div class="gmail_quote">On Sat, Oct 3, 2015 at 1:00 PM,  <span dir="ltr"><<a href="mailto:flash-users-request@flash.uchicago.edu" target="_blank">flash-users-request@flash.uchicago.edu</a>></span> wrote:<br><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 id=":14o" class="" style="overflow:hidden">From: Ernesto Zurbriggen <<a href="mailto:ernesto@mail.oac.uncor.edu">ernesto@mail.oac.uncor.edu</a>><br>To: <a href="mailto:flash-users@flash.uchicago.edu">flash-users@flash.uchicago.edu</a><br>Cc: <br>Date: Fri, 02 Oct 2015 17:45:40 -0300<br>Subject: Re: [FLASH-USERS] value of Ye for single species runs<br>Jon,<br>
looking at the end of the subroutine Eos/EosMain/Eos_getAbarZbar.F90 you can see for a single species run Ye=Zbar/Abar too. And I think that is not a stub implementation for the single-species gamma-law EOS.<br>
<br>
On the other hand, if you turn on the thermal conduction managed by Diffuse/Split or Diffuse/Unplit units, I think the runtime parameter "addThermalFlux" must be false. On the contrary, you might be taking into account thermal conduction twice, within the Diffuse unit and Hydro unit via hy_uhd_addThermalFluxes.F90<br>
<br>
Best!<br></div></blockquote></div><br><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr">________________________________________________________<br>Jonathan D. Slavin                 Harvard-Smithsonian CfA<br><a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>       60 Garden Street, MS 83<br>phone: (617) 496-7981       Cambridge, MA 02138-1516<br>cell: (781) 363-0035             USA<br>________________________________________________________<br><br></div></div></div></div>
</div></div>