<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
<div>Jon,</div>
<div>I forgot to say;</div>
<div>V=r/t, and Rho = r^(-n)* t^(-3) for the power-law ejecta. You need to have a non-zero starting time.</div>
<div>Given the time, V is just proportional to r.</div>
<div>The built-in BC won't work for SNR simulations.</div>
<div>Hope this help.</div>
<div>- Yueh</div>
<div><br>
</div>
<div><br>
</div>
<div>Jon,</div>
<div><br>
</div>
<div>The free expansion ejecta must have a velocity profile V=R/T.</div>
<div>Your inner BC is NOT placed at zero, which is good. Inside the inner BC, i.e. the inn ner ghostzones, you still need to set up the velocity and density profile so that the flow satisfies the free expansion. If the free expansion is not satisfied, the flow
will first clash at the inner boundary, and the numbers (pressure and density) may eventually run 'wild.'</div>
<div><br>
</div>
<div>1-D spherical POLAR symmetry: don't use theta near 0 (the pole). Try theta (th) at theta=pi/2 (equator) and a reasonably small delta theta (dth). </div>
<div><br>
</div>
<div>I think some trials and errors on the BC (and IC) would help in this cases.</div>
<div><br>
</div>
<div>-Yueh</div>
<div><br>
</div>
</div>
<div id="signature">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Arial,Helvetica,sans-serif">
<font face="Calibri" size="3" style="font-size:12pt"><span style="line-height:22.399999618530273px"><br>
</span></font></div>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> flash-users-bounces@flash.uchicago.edu <flash-users-bounces@flash.uchicago.edu> on behalf of Slavin, Jonathan <jslavin@cfa.harvard.edu><br>
<b>Sent:</b> Thursday, March 29, 2018 9:29 PM<br>
<b>To:</b> flash-users<br>
<b>Subject:</b> [FLASH-USERS] noisiness in run with high ekin/eint</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div class="x_gmail_default" style="font-family:tahoma,sans-serif; font-size:small">
Hi,</div>
<div class="x_gmail_default" style="font-family:tahoma,sans-serif; font-size:small">
<br>
</div>
<div class="x_gmail_default" style="font-family:tahoma,sans-serif; font-size:small">
I've been doing runs in 1D spherical symmetry of supernova remnants starting with fast cool ejecta. This is one of those situations in which the kinetic energy far exceeds the thermal energy (at least in the early stages). </div>
<div class="x_gmail_default" style="font-family:tahoma,sans-serif; font-size:small">
<br>
</div>
<div class="x_gmail_default" style="font-family:tahoma,sans-serif; font-size:small">
I know about eintSwitch which is supposed to help with numerical issues that arise in such situations, but I've been finding that setting it lower actually results in less noisiness, though the noisiness that remains is, if nothing else, unsightly. What happens
is that the pressure develops fairly large fluctuations that seem to be purely numerical since there doesn't seem to be any other source for them. The pressure is small in those regions, relative to the kinetic energy, but still the noise seems to affect
the density. </div>
<div class="x_gmail_default" style="font-family:tahoma,sans-serif; font-size:small">
<br>
</div>
<div class="x_gmail_default" style="font-family:tahoma,sans-serif; font-size:small">
So I'm wondering if others have found parameter settings or approaches that work well in this sort of situation. If so, I'd like to hear about them. I should mention that I've been using split PPM. I did try to use unsplit hydro but it failed with the pressure
dropping to the smallp value very quickly. I may try to investigate that some more.</div>
<div class="x_gmail_default" style="font-family:tahoma,sans-serif; font-size:small">
<br>
</div>
<div class="x_gmail_default" style="font-family:tahoma,sans-serif; font-size:small">
Regards,</div>
<div class="x_gmail_default" style="font-family:tahoma,sans-serif; font-size:small">
Jon</div>
<div><br>
</div>
-- <br>
<div class="x_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>
</body>
</html>