<div dir="ltr"><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">Thanks for your reply. I'm still not clear though on whether there is a problem with doing split hydro and unsplit diffusion. I would think that it's not a problem since the usual practice is use operator splitting such that the hydro evolution is done separately from the diffusion, but I don't know that for a fact. Can anyone comment?</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 class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 5, 2015 at 5:25 PM, Ernesto Zurbriggen <span dir="ltr"><<a href="mailto:ernesto@mail.oac.uncor.edu" target="_blank">ernesto@mail.oac.uncor.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Jon,<br>
integrate using Hydro/HydroMain/split or Hydro/HydroMain/unsplit is very different. That is if you put the shortcut "+uhd" in the setup line or not (in the hydro case). Things I told you are for the unsplit case and using Diffuse/DiffuseMain/Unsplit.<br>
<br>
While for the case Hydro/HydroMain/split and Diffuse/DiffuseMain/split (maybe also with Diffuse/DiffuseMain/unsplit) the parameter "addThermalFlux" is unkown for the code.<br>
<br>
On the other hand, every time you modify something in your flash.par or Config file the best would be run the setup line again.<div class="HOEnZb"><div class="h5"><br>
<br>
<br>
<br>
Quoting "Slavin, Jonathan" <<a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>>:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Ernesto,<br>
<br>
I just tried the run again with the parameter addThermalFlux set to False<br>
and got in the log file:<br>
WARNING: Ignored Parameters :<br>
These parameters were found in the flash.par file, but they were<br>
not declared in any Config file for the simulation!<br>
...<br>
addThermalFlux<br>
<br>
so it seems that the code doesn't use addThermalFlux. Could the parameter<br>
have some other name? I didn't rerun setup or build, but I wouldn't think<br>
that would be necessary for a runtime parameter change. Do you think it's<br>
a problem to use split PPM and unsplit diffusion? I think the diffusion is<br>
done separately from the hydro via operator splitting.<br>
<br>
Regards,<br>
Jon<br>
<br>
On Mon, Oct 5, 2015 at 11:53 AM, Ernesto Zurbriggen <<br>
<a href="mailto:ernesto@mail.oac.uncor.edu" target="_blank">ernesto@mail.oac.uncor.edu</a>> wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Jon,<br>
as I understand from what you told us, you should put<br>
addThermalFlux=.false.<br>
Because of you are already considering thermal conduction in the Diffuse<br>
unit, so you need the code doesn't enter the subroutine<br>
hy_uhd_addThermalFluxes.F90. In case of keeping addThermalFlux=.true., you<br>
are considering thermal conduction two times. You need to turn off one of<br>
them.<br>
<br>
<br>
I think of the code has a bug putting addThermalFlux=.true. by default.<br>
That would be all.<br>
<br>
Best!<br>
<br>
<br>
<br>
Quoting "Slavin, Jonathan" <<a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>>:<br>
<br>
Hi Ernesto,<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
I don't specify the runtime parameter addThermalFlux. Instead I set<br>
diff_useEleCond, useConductivity, useDiffuse, and useDiffuseTherm to<br>
.true.<br>
and specify<br>
REQUIRES<br>
physics/materialProperties/Conductivity/ConductivityMain/PowerLaw<br>
REQUESTS physics/Diffuse/DiffuseMain/Unsplit<br>
Should I change one of those?<br>
<br>
I note that in the log file that<br>
hy_fluxrespresentation = fluxes<br>
I'm not sure exactly what that means. My goal here is to verfy that<br>
thermal conduction is operating as I expect by comparing evaporation of a<br>
spherical cloud in 2D cylindrical symmetry with analytical results derived<br>
for steady flow including saturation of thermal conduction (Dalton &<br>
Balbus<br>
1993).<br>
<br>
Regards,<br>
Jon<br>
<br>
On Sat, Oct 3, 2015 at 1:00 PM, <<a href="mailto:flash-users-request@flash.uchicago.edu" target="_blank">flash-users-request@flash.uchicago.edu</a>><br>
wrote:<br>
<br>
From: Ernesto Zurbriggen <<a href="mailto:ernesto@mail.oac.uncor.edu" target="_blank">ernesto@mail.oac.uncor.edu</a>><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
To: <a href="mailto:flash-users@flash.uchicago.edu" target="_blank">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<br>
can see for a single species run Ye=Zbar/Abar too. And I think that is<br>
not<br>
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<br>
Diffuse/Split or Diffuse/Unplit units, I think the runtime parameter<br>
"addThermalFlux" must be false. On the contrary, you might be taking into<br>
account thermal conduction twice, within the Diffuse unit and Hydro unit<br>
via hy_uhd_addThermalFluxes.F90<br>
<br>
Best!<br>
<br>
<br>
</blockquote>
<br>
<br>
<br>
--<br>
________________________________________________________<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: <a href="tel:%28617%29%20496-7981" value="+16174967981" target="_blank">(617) 496-7981</a> Cambridge, MA 02138-1516<br>
cell: <a href="tel:%28781%29%20363-0035" value="+17813630035" target="_blank">(781) 363-0035</a> USA<br>
________________________________________________________<br>
<br>
<br>
</blockquote>
<br>
<br>
----------------------------------------------------------------<br>
This message was sent using IMP, the Internet Messaging Program.<br>
<br>
<br>
</blockquote>
<br>
<br>
--<br>
________________________________________________________<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: <a href="tel:%28617%29%20496-7981" value="+16174967981" target="_blank">(617) 496-7981</a> Cambridge, MA 02138-1516<br>
cell: <a href="tel:%28781%29%20363-0035" value="+17813630035" target="_blank">(781) 363-0035</a> USA<br>
________________________________________________________<br>
<br>
</blockquote>
<br>
<br>
<br>
----------------------------------------------------------------<br>
This message was sent using IMP, the Internet Messaging Program.<br>
<br>
</div></div></blockquote></div><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>