<div dir="ltr"><div>Hi Jonathan, Marissa, Tomek :</div><div><br></div> Allow me to amend my prior statement -- for the point-blast case, the code is defining the pressure of the point-like initial condition as<div>
<p class="gmail-p1"><span class="gmail-s1"> sim_pExp </span><span class="gmail-s2">=</span><span class="gmail-s1"> (sim_gamma</span><span class="gmail-s2">-</span><span class="gmail-s3">1</span><span class="gmail-s1">.) </span><span class="gmail-s2">*</span><span class="gmail-s1"> sim_expEnergy </span><span class="gmail-s2">/</span><span class="gmail-s1"> sim_</span><span class="gmail-s4">vctr</span></p><p class="gmail-p1"><span class="gmail-s4">in Simulation_init.F90. The parameter sim_vctr is in turn defined as the volume, area, or length within a region of sim_rInit in 3D, 2D, and 1D, respectively. The total internal energy <img alt="\int e_{\rm int} dV" title="\int e_{\rm int} dV" class="va_li" src="https://latex.codecogs.com/gif.latex?%5Cdpi%7B300%7D%5Cinline%09%5Cint%09e%5F{%5Crm%09int}%09dV" id="l0.7804726829727564" style="display: inline; vertical-align: -5.134px;" height="19" width="58"> = <img alt="\int p / [(\gamma - 1) \rho] dV" title="\int p / [(\gamma - 1) \rho] dV" class="va_li" src="https://latex.codecogs.com/gif.latex?%5Cdpi%7B300%7D%5Cinline%09%5Cint%09p%09/%09[(%5Cgamma%09-%091)%09%5Crho]%09dV" id="l0.2844168348743963" height="19" style="display: inline; vertical-align: -5.134px;" width="120"> should identically return a total energy of sim_expEnergy, <i>for an initial density of 1, on a continuous domain</i>. Because of discretization, however, this will not be rigorously true, but you should be able to achieve as close an approximation as you want by increasing the effective resolution of the initialization -- that is, either increasing the sim_rInit parameter or adding more resolution.</span></p><p class="gmail-p1"> All of the previous caveats for point-like initialization still apply.</p><p class="gmail-p1"> Best wishes,</p><p class="gmail-p1"> Bob</p></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jul 26, 2016 at 11:25 AM, Marissa Adams <span dir="ltr"><<a href="mailto:madams@pas.rochester.edu" target="_blank">madams@pas.rochester.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">We understand that. The total energy is less than the input value. This is the problem.<br></div></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">On Tue, Jul 26, 2016 at 11:23 AM, Tomasz Plewa <span dir="ltr"><<a href="mailto:tplewa@fsu.edu" target="_blank">tplewa@fsu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
<tt>The relevant data is the total energy shown in the first line of
the dat file. That should stay constant with relative accuracy
better than 1e-10.<br>
<br>
The actual initial value of the total will vary as one changes
problem parameters.<br>
<br>
Tomek<br>
--</tt><div><div><br>
<div>On 07/26/16 10:30, Slavin, Jonathan
wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<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>
</div>
<div class="gmail_default" style="font-family:arial,helvetica,sans-serif">I've been
running supernova remnant evolution simulations. The
initialization is done with code taken from the supplied Sedov
simulation with certain small changes. One issue that I've
run into is that the total energy is not very accurate. That
is, though the energy is conserved, the total energy,
initially input as thermal energy, is below that which is
given in the sim_expEnergy parameter by ~5%. Has anyone else
run into this? I don't think that it has to do with any of
the changes that I've made to Simulation_initBlock or
Simulation_init. I'm wondering if increasing the
sim_nsubzones parameter might help. Any advice would be
appreciated.</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,</div>
<div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Jon</div>
<div><br>
</div>
-- <br>
<div data-smartmail="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: <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>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<br>
</div></div></div>
</blockquote></div><br><br clear="all"><br></div></div><span class="">-- <br><div data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><span><font size="1"><span style="color:rgb(0,0,0)"><span style="font-family:monospace,monospace"><div>Marissa Adams <span><font size="1"><span style="color:rgb(0,0,0)"><span style="font-family:monospace,monospace">E-mail:</span></span><span style="color:rgb(0,0,0)"><span style="font-family:monospace,monospace"><a> </a><a href="mailto:madams@pas.rochester.edu" target="_blank">madams@pas.rochester.edu</a></span></span></font></span><br>Graduate Student <span><font size="1"><span style="color:rgb(0,0,0)"><span style="font-family:monospace,monospace">Ph: <a href="tel:%28585%29%20402-5779" value="+15854025779" target="_blank">(585) 402-5779</a></span></span></font></span></div><div>Department of Physics & Astronomy <span><font size="1"><span style="color:rgb(0,0,0)"><span style="font-family:monospace,monospace">Website: <a href="http://www.pas.rochester.edu/%7Emadams/" target="_blank">http://www.pas.rochester.edu/~madams/</a></span></span></font></span></div>University of Rochester<br>478 Bausch & Lomb Hall<br>P.O. Box 270171<br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Dr. Robert Fisher<br>Associate Professor / Graduate Program Director<br>University of Massachusetts/Dartmouth<br>Department of Physics<br>285 Old Westport Road<br>North Dartmouth, Massachusetts 02747<br><a href="mailto:robert.fisher@umassd.edu" target="_blank">robert.fisher@umassd.edu</a><br><a href="http://www.novastella.org" target="_blank">http://www.novastella.org</a></div></div>
</span></span></font></span></div></div>
</div></div></div></div></div></div></span></div></blockquote></div></div>