<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Hi Dongwook,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">I was using HLLD. However my problem was solved when I set the initial time step after the second explosion to dtinit, which is a very small value. So I think it was related to the large pressure gradients present when the second explosion is initiated. As far as I can tell, the problem was between fine resolution pixels, not at a fine-coarse boundary.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Jon</div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 8, 2017 at 7:05 PM, Dongwook Lee <span dir="ltr"><<a href="mailto:dongwook@flash.uchicago.edu" target="_blank">dongwook@flash.uchicago.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">Hi Jon,<div><br></div><div>Do you see large B fluctuations aligned with the fine-coarse block boundaries; and they get larger over time?</div><div><br></div><div>Just to add one quick suggestion: what Riemann solver you've been using? If HLLD, can you try to switch to HLLC or HLL and see if the fluctuations still persist?</div><div><br></div><div>Cheers,</div><div>Dongwook</div></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">On Fri, May 26, 2017 at 1:38 PM, Slavin, Jonathan <span dir="ltr"><<a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.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 style="font-family:arial,helvetica,sans-serif;font-size:small">To answer Sean and Ernesto's questions:</div><div style="font-family:arial,helvetica,sans-serif;font-size:small">I am using the USM solver and AMR. I have now determined that div B is still 0 after the second explosion, which is not too surprising since, for this run, I have not changed the magnetic field and include div B cleaning (killdivb = .true.). So the problem is not with non-zero div B but with the large B fluctuations that are initiated. I noticed that there are small scale fluctuations in the region where the second explosion is initiated before the explosion, though the field is weak ~ 0.01 muG. The fluctuations after the explosion go from ~ -1000 muG - +1000 muG. I do expect amplification of B field fluctuations at the shock, but not to that extent.</div><div style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div style="font-family:arial,helvetica,sans-serif;font-size:small">The way I'm initiating the second explosion is (almost) the same way I initiate the first one - I set the pressure of parcels within a given radius to a value such that the sum of the energy in those parcels is equal to the value I want for the explosion (1.E51 ergs). I set the density to a value so as to have the mass total to 8 solar masses. So all the energy is thermal at first. I'll see about calling the EOS unit, since I haven't done that. Any additional info on how to do that would be appreciated.</div><div style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div style="font-family:arial,helvetica,sans-serif;font-size:small">Thanks,</div><div style="font-family:arial,helvetica,sans-serif;font-size:small">Jon</div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 26, 2017 at 1:00 PM, <span dir="ltr"><<a href="mailto:flash-users-request@flash.uchicago.edu" target="_blank">flash-users-request@flash.uch<wbr>icago.edu</a>></span> wrote:<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span>---------- Forwarded message ----------<br>From: ERNESTO ZURBRIGGEN <<a href="mailto:ezurbriggen@unc.edu.ar" target="_blank">ezurbriggen@unc.edu.ar</a>><br>To: <a href="mailto:flash-users@flash.uchicago.edu" target="_blank">flash-users@flash.uchicago<wbr>.edu</a><br>Cc: <br>Bcc: <br>Date: Fri, 26 May 2017 10:56:45 -0300<br>Subject: Re: [FLASH-USERS] ensuring div B = 0 when adjusting evolution<br><div dir="ltr"><div>Hi Jon! </div><div><span style="font-family:arial,helvetica,sans-serif"><br></span></div><div><div><font face="arial, helvetica, sans-serif">Have you tried to set off a less intense second supernova explosion? Have you observed the same problems in that case? The second explotion might be much too intense. </font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">On the other hand, how do you set the later explotion? I mean, applying the explotion, are you consistently modifying the thermodynamical variables? For example, if you instantaneously modify the density and the temperature, then you also should call the Eos unit to keep the consistence. I have experimented situations in which just modifying velocities and keeping the thermodynamics unaltered, I also had to call the Eos unit in order to maintain consistence</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Some runtime parameter that might help being on are 'shockDetect' and 'shockLowerCFL' (this last one I think is just in realease 4.4). </font></div></div><div><br></div><div><br></div>Best!</div></span><div class="gmail_extra"><div><br></div><span>-- <br><div class="m_-5390688968651369799m_7673830170894663220m_8547132186880594962gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><font face="monospace, monospace"><i><font color="#000000">Ernesto Zurbriggen</font></i><i> </i></font><div><font face="monospace, monospace"><i><br></i></font><div><font face="monospace, monospace"><i><font color="#444444">Instituto de Astronomía Teórica y Experimental </font><font color="#274e13">(IATE)</font><font color="#444444">. </font></i></font></div><div><font face="monospace, monospace"><i><font color="#444444">Observatorio Astronómico de Córdoba </font><font color="#274e13">(OAC)</font><font color="#444444">, </font></i></font><i style="color:rgb(68,68,68);font-family:monospace,monospace;font-size:12.8px">Universidad Nacional de Córdoba </i><i style="font-family:monospace,monospace;font-size:12.8px"><font color="#274e13">(UNC)</font></i><i style="color:rgb(68,68,68);font-family:monospace,monospace;font-size:12.8px">. </i></div></div><div><i style="color:rgb(68,68,68);font-family:monospace,monospace;font-size:12.8px">Teléfono: +54 0351 4331064-5, interno 222. </i></div><div><font color="#444444"><font face="monospace, monospace"><i>Córdoba, Argentina.</i></font> </font><br></div></div></div></div></div></div></div></div>
</span></div><div><div class="m_-5390688968651369799h5">
<br><br>---------- Forwarded message ----------<br>From: "Sean M. Couch" <<a href="mailto:couch@pa.msu.edu" target="_blank">couch@pa.msu.edu</a>><br>To: "Slavin, Jonathan" <<a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>>, flash-users <<a href="mailto:flash-users@flash.uchicago.edu" target="_blank">flash-users@flash.uchicago.ed<wbr>u</a>><br>Cc: <br>Bcc: <br>Date: Fri, 26 May 2017 14:07:22 +0000<br>Subject: Re: [FLASH-USERS] ensuring div B = 0 when adjusting evolution<br>
<div style="word-wrap:break-word">
<div id="m_-5390688968651369799m_7673830170894663220m_5462410415856234389bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px;color:rgba(0,0,0,1.0);margin:0px;line-height:auto">
Hi Jon,</div>
<div id="m_-5390688968651369799m_7673830170894663220m_5462410415856234389bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px;color:rgba(0,0,0,1.0);margin:0px;line-height:auto">
<br>
</div>
<div id="m_-5390688968651369799m_7673830170894663220m_5462410415856234389bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px;color:rgba(0,0,0,1.0);margin:0px;line-height:auto">
Can you give a little more info? Are you using the USM solver? Are you using AMR? Have you actually checked that divB>0 in the output data? A log file from a representative run would be useful.</div>
<div id="m_-5390688968651369799m_7673830170894663220m_5462410415856234389bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px;color:rgba(0,0,0,1.0);margin:0px;line-height:auto">
<br>
</div>
<div id="m_-5390688968651369799m_7673830170894663220m_5462410415856234389bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px;color:rgba(0,0,0,1.0);margin:0px;line-height:auto">
In my experience, the USM solver in cylindrical R-Z coordinates with AMR can be….touchy. But it should work and maintain divB=0! (See, e.g., <a href="http://adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=2013ApJ...773..136J&link_type=EJOURNAL" target="_blank">http://adsabs.harvard.ed<wbr>u/cgi-bin/nph-data_query?bibco<wbr>de=2013ApJ...773..136J&link_ty<wbr>pe=EJOURNAL</a>).</div>
<div id="m_-5390688968651369799m_7673830170894663220m_5462410415856234389bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px;color:rgba(0,0,0,1.0);margin:0px;line-height:auto">
<br>
</div>
<div id="m_-5390688968651369799m_7673830170894663220m_5462410415856234389bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px;color:rgba(0,0,0,1.0);margin:0px;line-height:auto">
Cheers,</div>
<div id="m_-5390688968651369799m_7673830170894663220m_5462410415856234389bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px;color:rgba(0,0,0,1.0);margin:0px;line-height:auto">
Sean</div>
<br>
<div class="m_-5390688968651369799m_7673830170894663220m_5462410415856234389bloop_sign" id="m_-5390688968651369799m_7673830170894663220m_5462410415856234389bloop_sign_1495807385531248128">
<div style="font-family:Helvetica;font-size:14px;word-wrap:break-word">
<br>
------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>----------------<br>
Sean M. Couch</div>
<div style="font-family:Helvetica;font-size:14px;word-wrap:break-word">
Assistant Professor<br>
Department of Physics and Astronomy<br>
Department of Computational Mathematics, Science, and Engineering<br>
National Superconducting Cyclotron Laboratory/Facility for Rare Isotope Beams<br>
Michigan State University<br>
567 Wilson Rd, 3250 BPS<br>
East Lansing, MI 48824<br>
<a href="tel:(517)%20884-5035" value="+15178845035" target="_blank">(517) 884-5035</a> —— <a href="mailto:couch@pa.msu.edu" target="_blank">couch@pa.msu.edu</a> —— <a href="http://www.pa.msu.edu/~couch" target="_blank">www.pa.msu.edu/~couch</a></div>
<div style="font-family:Helvetica;font-size:14px;word-wrap:break-word">
<br>
<div><br>
<br>
</div>
</div>
</div>
<br>
<p class="m_-5390688968651369799m_7673830170894663220m_5462410415856234389airmail_on">On May 25, 2017 at 4:41:45 PM, Slavin, Jonathan (<a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>) wrote:</p>
<blockquote type="cite" class="m_-5390688968651369799m_7673830170894663220m_5462410415856234389clean_bq"><span>
<div>
<div></div>
<div>
<div dir="ltr">
<div style="font-family:arial,helvetica,sans-serif;font-size:small">
Hi all,</div>
<div style="font-family:arial,helvetica,sans-serif;font-size:small">
<br>
</div>
<div style="font-family:arial,helvetica,sans-serif;font-size:small">
I'm running MHD simulations where I set off a second supernova explosion within a pre-existing remnant. I'm using Simulation_adjustEvolution for this. However, I'm running into a problem with the magnetic field just after initiating the second explosion.
I'm getting very large pixel-to-pixel variations in Bx and By at the edge of the new expanding blast wave. I'm thinking that it could be because of non-zero div B in the region in which the explosion is generated. So my question is, does anyone have a suggestion
for div B cleaning at the point that the explosion is started, i.e. within Simulation_adjustEvolution, to prevent the problems I'm having? I don't really expect the B field to be dynamically important inside the remnant (beta >> 1), so accuracy of the B field
is probably not important there. I should add that I'm doing these calculations in 2D cylindrical symmetry (R-Z). Thanks in advance for any help.</div>
<div style="font-family:arial,helvetica,sans-serif;font-size:small">
<br>
</div>
<div style="font-family:arial,helvetica,sans-serif;font-size:small">
Regards,</div>
<div style="font-family:arial,helvetica,sans-serif;font-size:small">
Jon</div>
</div></div></div></span></blockquote></div><br></div></div>______________________________<wbr>_________________<br>
flash-users mailing list<br>
<a href="mailto:flash-users@flash.uchicago.edu" target="_blank">flash-users@flash.uchicago.edu</a><br>
<a href="http://flash.uchicago.edu/mailman/listinfo/flash-users" rel="noreferrer" target="_blank">http://flash.uchicago.edu/mail<wbr>man/listinfo/flash-users</a><br></blockquote></div><span>-- <br><div class="m_-5390688968651369799m_7673830170894663220gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">______________________________<wbr>__________________________<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:(617)%20496-7981" value="+16174967981" target="_blank">(617) 496-7981</a> Cambridge, MA 02138-1516<br>cell: <a href="tel:(781)%20363-0035" value="+17813630035" target="_blank">(781) 363-0035</a> USA<br>______________________________<wbr>__________________________<br><br></div></div></div></div>
</span></div></div>
</blockquote></div><br><br clear="all"><div><br></div></div></div><span class="HOEnZb"><font color="#888888">-- <br><div class="m_-5390688968651369799gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif"><br>==============================<wbr>===========</font></div><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif">Dongwook Lee, Ph.D., Assistant Professor</font></div><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif">Applied Mathematics and Statistics<br></font></div><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif">University of California, Santa Cruz</font></div><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif">Baskin Engineering, Room 353C</font></div><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif">1156 High Street, Santa Cruz, CA 95064</font></div><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif"><a href="https://users.soe.ucsc.edu/~dongwook/" target="_blank">https://users.soe.ucsc.edu/~<wbr>dongwook/</a></font></div><div><br></div></div></div>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" 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: (617) 496-7981 Cambridge, MA 02138-1516<br>cell: (781) 363-0035 USA<br>________________________________________________________<br><br></div></div></div></div>
</div></div>