<div dir="ltr">Hi Binh,<div><br></div><div>I think the density inside the rarefied region is going below the value of the runtime parameter "smlrho", which defaults to 1e-10. You can try lowering that in your flash.par. <br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr" style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px">*************************************************************************</div><div dir="ltr" style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px">Adam Reyes<div><img src="https://flash.rochester.edu/site/common/images/email_signature/FLASH.jpg"><br></div><div>Code Group Leader, Flash Center for Computational Science <span> </span><br></div><div>Research Scientist, Dept. of Physics and Astronomy</div><div>University of Rochester</div><div>River Campus: Bausch and Lomb Hall, 369 <span> </span><br></div><div>500 Wilson Blvd. PO Box 270171, Rochester, NY 14627</div><div dir="ltr">Email<span> </span><a href="mailto:adam.reyes@rochester.edu" target="_blank">adam.reyes@rochester.edu</a></div><div>Web<span> </span><a href="https://flash.rochester.edu/" target="_blank">https://flash.rochester.edu</a><br></div><div> (he / him / his)</div><div><img src="https://flash.rochester.edu/site/common/images/email_signature/FLASH-pride-sml.png"><br></div>*************************************************************************</div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jul 21, 2023 at 3:53 PM Yên Bình <<a href="mailto:yenbinhpy308@gmail.com">yenbinhpy308@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear Adam Reyes, <br></div><div><br></div><div>I can change the explosion energy to 10^51 erg but when I try to minimize the initial radius to below 10^14, the simulation cracks and sometimes produces weird results like in the photo.</div><div><br></div><div>It does not seem to be related to other runtime parameters, it thinks it is because of the grid ? I try to set nxb and nyb to higher value like 200 but it output this error:</div><div><i>Abort(1) on node 0 (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0 <br></i></div><div><i><br></i></div><div>How can I set the initial radius to 10^4 like a star core before an explosion? <br></div><div><br></div><div>Thank you for your kind support,</div><div><br></div><div>Binh<br></div><div><i></i></div><div><br></div><div><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><img src="cid:ii_lkcn3jls1" alt="image.png" width="513" height="492"><br><br><br><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:0px 1px 20px 0px"><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:0px 1px 0px 0px"><p style="margin:1px"><br></p></td></tr></tbody></table></td></tr><tr><td style="padding:0px 1px 25px 0px;font-family:Verdana,Geneva,sans-serif;font-size:13px;line-height:17px;white-space:nowrap"><p style="line-height:17px;font-weight:700;color:rgb(0,0,0);margin:1px">Nguyen Thi Yen Binh</p><p style="line-height:17px;color:rgb(136,136,136);margin:1px"><br></p><p style="line-height:17px;color:rgb(136,136,136);margin:1px">Summer Student, Academia Sinica Institute of Astronomy and Astrophysics<br></p><p style="line-height:17px;color:rgb(136,136,136);margin:1px">Space Science and Satellite Technology, University of Science and Technology of Ha Noi</p></td></tr><tr><td style="padding:0px 1px 0px 0px"><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:25px 1px 25px 0px;border-top-width:2px;border-top-style:solid;border-top-color:rgb(123,146,152);border-bottom-width:2px;border-bottom-style:solid;border-bottom-color:rgb(123,146,152)"><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:0px 1px 0px 0px"><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:1px 5px 1px 0px;vertical-align:middle" valign="middle"><p style="margin:1px"><img src="https://signaturehound.com/api/v1/png/email/default/fa0000.png" alt="" style="display: block; border: 0px; max-width: 18px;" width="18" height="18"></p></td><td style="line-height:17px;padding:1px 0px;vertical-align:middle"><p style="margin:1px"><a href="mailto:binhnty.bi12-060@st.usth.edu.vn" style="font-family:Verdana,Geneva,sans-serif;font-size:13px;line-height:17px;white-space:nowrap;color:rgb(136,136,136)" target="_blank">binhnty.bi12-060@st.usth.edu.vn</a></p></td></tr><tr><td style="padding:1px 5px 1px 0px;vertical-align:middle" valign="middle"><p style="margin:1px"><img src="https://signaturehound.com/api/v1/png/telephone/default/fa0000.png" alt="" style="display: block; border: 0px; max-width: 18px;" width="18" height="18"></p></td><td style="line-height:17px;padding:1px 0px;vertical-align:middle"><p style="margin:1px"><a href="tel:+84367533986" style="font-family:Verdana,Geneva,sans-serif;font-size:13px;line-height:17px;white-space:nowrap;color:rgb(136,136,136)" target="_blank">+84367533986</a></p></td></tr><tr><td style="padding:1px 5px 1px 0px;vertical-align:middle" valign="middle"><p style="margin:1px"><img src="https://signaturehound.com/api/v1/png/map/default/fa0000.png" alt="" style="display: block; border: 0px; max-width: 18px;" width="18" height="18"></p></td><td style="line-height:17px;padding:1px 0px;vertical-align:middle"><p style="margin:1px"><span style="font-family:Verdana,Geneva,sans-serif;font-size:13px;line-height:17px;white-space:nowrap;color:rgb(136,136,136)">Tuy Hoa, Phu Yen, Viet Nam</span></p></td></tr></tbody></table></td><td style="padding:0px 1px 0px 0px" width="15"></td></tr></tbody></table></td></tr></tbody></table></td></tr><tr><td style="padding:0px 1px 0px 0px"><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="font-size:0px;line-height:0px;padding:25px 1px 0px 0px" width="30"><p style="margin:1px"><a href="https://urldefense.com/v3/__https://www.facebook.com/wbinhw__;!!CGUSO5OYRnA7CQ!eHzY5sAvt9JVcQA5cJKZc6ry3jGnKBHVycUBNVDfU0KfAwSYw4QTVhFLbkxC7sqZjy9eKhHdhKa2njTzIDK0EINBdi0$" target="_blank"><img src="https://signaturehound.com/api/v1/png/facebook/default/1d62b0.png" alt="" style="display: block; border: 0px; max-width: 30px;" width="30" height="30"></a></p></td><td style="padding:0px 0px 1px" width="3"></td><td style="font-size:0px;line-height:0px;padding:25px 1px 0px 0px" width="30"><p style="margin:1px"><a href="https://urldefense.com/v3/__https://www.instagram.com/yenbinne/__;!!CGUSO5OYRnA7CQ!eHzY5sAvt9JVcQA5cJKZc6ry3jGnKBHVycUBNVDfU0KfAwSYw4QTVhFLbkxC7sqZjy9eKhHdhKa2njTzIDK0HytT4xk$" target="_blank"><img src="https://signaturehound.com/api/v1/png/instagram/default/1d62b0.png" alt="" style="display: block; border: 0px; max-width: 30px;" width="30" height="30"></a></p></td><td style="padding:0px 0px 1px" width="3"></td><td style="font-size:0px;line-height:0px;padding:25px 1px 0px 0px" width="30"><p style="margin:1px"><a href="https://urldefense.com/v3/__https://github.com/pea9__;!!CGUSO5OYRnA7CQ!eHzY5sAvt9JVcQA5cJKZc6ry3jGnKBHVycUBNVDfU0KfAwSYw4QTVhFLbkxC7sqZjy9eKhHdhKa2njTzIDK0Q3WW6Mc$" target="_blank"><img src="https://signaturehound.com/api/v1/png/github/default/1d62b0.png" alt="" style="display: block; border: 0px; max-width: 30px;" width="30" height="30"></a></p></td><td style="padding:0px 0px 1px" width="3"></td></tr></tbody></table></td></tr><tr><td style="padding:0px 1px 0px 0px"><br></td></tr></tbody></table><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:0px 1px 20px 0px"><br></td></tr><tr><td style="padding:0px 1px 25px 0px;font-family:Verdana,Geneva,sans-serif;font-size:13px;line-height:17px;white-space:nowrap"><br></td></tr><tr><td style="padding:0px 1px 0px 0px"><br></td></tr><tr><td style="padding:0px 1px 0px 0px"><br></td></tr><tr><td style="padding:0px 1px 0px 0px"><br></td></tr></tbody></table></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Vào Th 4, 19 thg 7, 2023 vào lúc 16:43 Adam Reyes <<a href="mailto:adam.reyes@rochester.edu" target="_blank">adam.reyes@rochester.edu</a>> đã viết:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi Binh,<div><br></div><div>You probably want to play with the runtime parameters "sim_rInit" and "sim_expEnergy" in your flash.par. These control the energy of the explosion and the initial radius in which it is deposited for the initial conditions. You'll want to scale these to your much larger domain. </div><div><br></div><div>The FLASH setup of the problem including the many runtime parameters as well as the dimensionless analytic solution are described <a href="https://flash.rochester.edu/site/flashcode/user_support/flash4_ug_4p7/node187.html#SECTION010114000000000000000" target="_blank">here</a> in the user guide.<br><br clear="all"><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><div dir="ltr" style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px">*************************************************************************</div><div dir="ltr" style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px">Adam Reyes<div><img src="https://flash.rochester.edu/site/common/images/email_signature/FLASH.jpg"><br></div><div>Code Group Leader, Flash Center for Computational Science <span> </span><br></div><div>Research Scientist, Dept. of Physics and Astronomy</div><div>University of Rochester</div><div>River Campus: Bausch and Lomb Hall, 369 <span> </span><br></div><div>500 Wilson Blvd. PO Box 270171, Rochester, NY 14627</div><div dir="ltr">Email<span> </span><a href="mailto:adam.reyes@rochester.edu" target="_blank">adam.reyes@rochester.edu</a></div><div>Web<span> </span><a href="https://flash.rochester.edu/" target="_blank">https://flash.rochester.edu</a><br></div><div> (he / him / his)</div><div><img src="https://flash.rochester.edu/site/common/images/email_signature/FLASH-pride-sml.png"><br></div>*************************************************************************</div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jul 19, 2023 at 7:44 AM Yên Bình <<a href="mailto:yenbinhpy308@gmail.com" target="_blank">yenbinhpy308@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear FLASH Users,<div><br></div><div>I have a problem with setting up the code for Sedov. I try to set up the simulation box to 10^16 x 10^16:</div><i>sim_xctr = 5.0e15<br></i><div><i>sim_yctr = 5.0e15</i> </div><div><i>xmin = 0.<br>xmax = 1.0e16<br>ymin = 0.<br>ymax = 1.0e16</i><br></div><div><br></div><div>But when I try to run the simulation, the plot file does not have any sign of an explosion (just a blank). I think it's because the box is too big so I need to set the right value of the running time:</div><div><br></div><div><i># checkpoint file output parameters<br>checkpointFileIntervalTime = 0.01<br>checkpointFileIntervalStep = 0<br>checkpointFileNumber = 0<br><br># plotfile output parameters<br>plotfileIntervalTime = 1<br>plotfileIntervalStep = 0<br>plotfileNumber = 0<br><br>nend = 10000<br>tmax = 1.0e8</i><br></div><div><br></div><div><div>But even if I set the time very high (10^8), nothing happens.</div><div><br></div><div>I really appreciate any help to solve this problem. Thank you in advance.</div><div><br></div><div>Binh.</div><div><br clear="all"><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><br><img src="cid:ii_lk9aj1ys0" alt="image.png" width="542" height="502"><br><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:0px 1px 20px 0px"><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:0px 1px 0px 0px"><p style="margin:1px"><br></p></td></tr></tbody></table></td></tr><tr><td style="padding:0px 1px 25px 0px;font-family:Verdana,Geneva,sans-serif;font-size:13px;line-height:17px;white-space:nowrap"><p style="line-height:17px;font-weight:700;color:rgb(0,0,0);margin:1px">Nguyen Thi Yen Binh</p><p style="line-height:17px;color:rgb(136,136,136);margin:1px">College Student</p><p style="line-height:17px;color:rgb(136,136,136);margin:1px">Space Science and Satellite Technology</p><p style="line-height:17px;color:rgb(136,136,136);margin:1px">University of Science and Technology of Ha Noi</p></td></tr><tr><td style="padding:0px 1px 0px 0px"><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:25px 1px 25px 0px;border-top-width:2px;border-top-style:solid;border-top-color:rgb(123,146,152);border-bottom-width:2px;border-bottom-style:solid;border-bottom-color:rgb(123,146,152)"><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:0px 1px 0px 0px"><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:1px 5px 1px 0px;vertical-align:middle" valign="middle"><p style="margin:1px"><img src="https://signaturehound.com/api/v1/png/email/default/fa0000.png" alt="" style="display: block; border: 0px; max-width: 18px;" width="18" height="18"></p></td><td style="line-height:17px;padding:1px 0px;vertical-align:middle"><p style="margin:1px"><a href="mailto:binhnty.bi12-060@st.usth.edu.vn" style="font-family:Verdana,Geneva,sans-serif;font-size:13px;line-height:17px;white-space:nowrap;color:rgb(136,136,136)" target="_blank">binhnty.bi12-060@st.usth.edu.vn</a></p></td></tr><tr><td style="padding:1px 5px 1px 0px;vertical-align:middle" valign="middle"><p style="margin:1px"><img src="https://signaturehound.com/api/v1/png/telephone/default/fa0000.png" alt="" style="display: block; border: 0px; max-width: 18px;" width="18" height="18"></p></td><td style="line-height:17px;padding:1px 0px;vertical-align:middle"><p style="margin:1px"><a href="tel:+84367533986" style="font-family:Verdana,Geneva,sans-serif;font-size:13px;line-height:17px;white-space:nowrap;color:rgb(136,136,136)" target="_blank">+84367533986</a></p></td></tr><tr><td style="padding:1px 5px 1px 0px;vertical-align:middle" valign="middle"><p style="margin:1px"><img src="https://signaturehound.com/api/v1/png/map/default/fa0000.png" alt="" style="display: block; border: 0px; max-width: 18px;" width="18" height="18"></p></td><td style="line-height:17px;padding:1px 0px;vertical-align:middle"><p style="margin:1px"><span style="font-family:Verdana,Geneva,sans-serif;font-size:13px;line-height:17px;white-space:nowrap;color:rgb(136,136,136)">Tuy Hoa, Phu Yen, Viet Nam</span></p></td></tr></tbody></table></td><td style="padding:0px 1px 0px 0px" width="15"></td></tr></tbody></table></td></tr></tbody></table></td></tr><tr><td style="padding:0px 1px 0px 0px"><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="font-size:0px;line-height:0px;padding:25px 1px 0px 0px" width="30"><p style="margin:1px"><a href="https://urldefense.com/v3/__https://www.facebook.com/wbinhw__;!!CGUSO5OYRnA7CQ!eHzY5sAvt9JVcQA5cJKZc6ry3jGnKBHVycUBNVDfU0KfAwSYw4QTVhFLbkxC7sqZjy9eKhHdhKa2njTzIDK0EINBdi0$" target="_blank"><img src="https://signaturehound.com/api/v1/png/facebook/default/1d62b0.png" alt="" style="display: block; border: 0px; max-width: 30px;" width="30" height="30"></a></p></td><td style="padding:0px 0px 1px" width="3"></td><td style="font-size:0px;line-height:0px;padding:25px 1px 0px 0px" width="30"><p style="margin:1px"><a href="https://urldefense.com/v3/__https://www.instagram.com/yenbinne/__;!!CGUSO5OYRnA7CQ!eHzY5sAvt9JVcQA5cJKZc6ry3jGnKBHVycUBNVDfU0KfAwSYw4QTVhFLbkxC7sqZjy9eKhHdhKa2njTzIDK0HytT4xk$" target="_blank"><img src="https://signaturehound.com/api/v1/png/instagram/default/1d62b0.png" alt="" style="display: block; border: 0px; max-width: 30px;" width="30" height="30"></a></p></td><td style="padding:0px 0px 1px" width="3"></td><td style="font-size:0px;line-height:0px;padding:25px 1px 0px 0px" width="30"><p style="margin:1px"><a href="https://urldefense.com/v3/__https://github.com/pea9__;!!CGUSO5OYRnA7CQ!eHzY5sAvt9JVcQA5cJKZc6ry3jGnKBHVycUBNVDfU0KfAwSYw4QTVhFLbkxC7sqZjy9eKhHdhKa2njTzIDK0Q3WW6Mc$" target="_blank"><img src="https://signaturehound.com/api/v1/png/github/default/1d62b0.png" alt="" style="display: block; border: 0px; max-width: 30px;" width="30" height="30"></a></p></td><td style="padding:0px 0px 1px" width="3"></td></tr></tbody></table></td></tr><tr><td style="padding:0px 1px 0px 0px"><br></td></tr></tbody></table><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:0px 1px 20px 0px"><br></td></tr><tr><td style="padding:0px 1px 25px 0px;font-family:Verdana,Geneva,sans-serif;font-size:13px;line-height:17px;white-space:nowrap"><br></td></tr><tr><td style="padding:0px 1px 0px 0px"><br></td></tr><tr><td style="padding:0px 1px 0px 0px"><br></td></tr><tr><td style="padding:0px 1px 0px 0px"><br></td></tr></tbody></table></div></div></div></div></div></div>
_______________________________________________<br>
flash-users mailing list<br>
<a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a><br>
<br>
For list info, including unsubscribe:<br>
<a href="https://flash.rochester.edu/mailman/listinfo/flash-users" rel="noreferrer" target="_blank">https://flash.rochester.edu/mailman/listinfo/flash-users</a><br>
</blockquote></div>
</blockquote></div>
</blockquote></div>