<div dir="ltr"><div dir="ltr">Hi Maksim! <div><br></div><div>It may happen your perturbation is too strong. It generates a ~10^4 pressure jump which is pretty much. </div><div>To support that you will need a strong setup. </div><div><br></div><div>You could try either beginning with a much less energetic pulse, let's say a ~5 pressure jump, </div><div>or: increase your spatial resolution; use a more diffusive riermannSolver and slopeLimiter; </div><div>setting smaller cfl values, etc. </div><div><br></div><div>Best for you! </div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El vie, 12 feb 2021 a las 1:56, Maksim Kozlov (<<a href="mailto:maksim.kozlov@nu.edu.kz" target="_blank">maksim.kozlov@nu.edu.kz</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div id="gmail-m_7138675922743617335gmail-m_-2663961290291028641gmail-:mw"><div dir="ltr">Dear Friends,<div>I would really appreciate your help with the following issue:</div><div>I
 am running a SBlast simulation in 3D - hydrodynamic expansion of initial perturbation with two fluids (3D graph of initial perturbation of internal energy is attached - energyZ). Simulation starts to execute 
perfectly - initial ball-shaped perturbation starts to expand but 
between 16th and 17th time step (please see attached slices of distributions of density, pressure and internal energy) something
 happens - in a few cells values of density, pressure and internal 
energy  become crazy ~10^38 and code crashes after that with Floating 
point exception - erroneous arithmetic operation. Any ideas on what can 
cause this problem or how can I track it down would be greatly appreciated!<br></div><div>Thank you very much!</div></div></div></div>
</blockquote></div>
<br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div dir="ltr" style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;vertical-align:baseline;color:inherit"><font face="monospace, monospace" style="background-color:rgb(234,209,220)"><font color="#000000"><span style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;margin:0px;padding:0px;border:0px;vertical-align:baseline;color:inherit"><div style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;margin:0px;padding:0px;border:0px;vertical-align:baseline;color:inherit"><i style="font-family:"Courier New",monospace;background-color:rgb(255,255,255);color:rgb(68,68,68);line-height:normal">Dr. Ernesto Zurbriggen, FAPESP posdoctoral fellow</i><i style="color:rgb(68,68,68)"><br></i></div></span></font></font><font face="monospace, monospace" style="background-color:rgb(234,209,220)"><font color="#000000"><span style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;margin:0px;padding:0px;border:0px;vertical-align:baseline;color:inherit"><div style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;margin:0px;padding:0px;border:0px;vertical-align:baseline;color:inherit"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;vertical-align:baseline;color:inherit;background-color:rgb(255,255,255)"><i style="color:rgb(68,68,68)"></i></span></div></span></font></font></div><div dir="ltr" style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;vertical-align:baseline;color:inherit"><div style="margin:0px;padding:0px;border:0px;font-stretch:inherit;line-height:inherit;font-family:Calibri,Arial,Helvetica,sans-serif;vertical-align:baseline;color:rgb(0,0,0)"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:"Courier New",monospace;vertical-align:baseline;color:inherit"><i style="color:rgb(68,68,68);line-height:normal"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;vertical-align:baseline;color:inherit"><i style="font-family:monospace,monospace;background-color:rgb(234,209,220)"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;vertical-align:baseline;color:inherit"><div style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;vertical-align:baseline;color:inherit"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:"Courier New",monospace;vertical-align:baseline;color:inherit;background-color:rgb(255,255,255)"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:inherit;vertical-align:baseline;color:inherit">Email</span></span><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:"Courier New",monospace;vertical-align:baseline;color:inherit;background-color:rgb(255,255,255)">:</span><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:inherit;vertical-align:baseline;color:inherit;background-color:rgb(255,255,255)"><i style="line-height:normal"> <a href="mailto:ernesto.zurbriggen@craam.mackenzie.br" target="_blank">ernesto.zurbriggen@craam.mackenzie.br</a></i></span><i><br></i></div><div style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;font-family:inherit;vertical-align:baseline;color:inherit"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:inherit;vertical-align:baseline;color:inherit;background-color:rgb(255,255,255)"><i style="line-height:normal">Phone: +55 11-2114-8783</i></span><i><br></i></div></span></i></span><br></i></span></div><div style="margin:0px;padding:0px;border:0px;font-stretch:inherit;line-height:inherit;font-family:Calibri,Arial,Helvetica,sans-serif;vertical-align:baseline;color:rgb(0,0,0)"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:"Courier New",monospace;vertical-align:baseline;color:inherit"><i style="color:rgb(68,68,68);line-height:normal">Universidade Presbiteriana Mackenzie, </i></span><i style="color:rgb(68,68,68);font-family:monospace,monospace;background-color:rgb(234,209,220)"><br></i></div><div style="margin:0px;padding:0px;border:0px;font-stretch:inherit;line-height:inherit;font-family:Calibri,Arial,Helvetica,sans-serif;vertical-align:baseline;color:rgb(0,0,0)"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:"Courier New",monospace;vertical-align:baseline;color:inherit"><i style="color:rgb(68,68,68);line-height:normal">Escola de Engenharia, </i></span><i style="color:rgb(68,68,68);font-family:monospace,monospace;background-color:rgb(234,209,220)"><br></i></div><div style="margin:0px;padding:0px;border:0px;font-stretch:inherit;line-height:inherit;font-family:Calibri,Arial,Helvetica,sans-serif;vertical-align:baseline;color:rgb(0,0,0)"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:"Courier New",monospace;vertical-align:baseline;color:inherit"><i style="color:rgb(68,68,68);line-height:normal">Centro de Rádio Astronomia e Astrofísica Mackenzie (CRAAM),</i></span><i style="color:rgb(68,68,68);font-family:monospace,monospace;background-color:rgb(234,209,220)"><br></i></div><div style="margin:0px;padding:0px;border:0px;font-stretch:inherit;line-height:inherit;font-family:Calibri,Arial,Helvetica,sans-serif;vertical-align:baseline;color:rgb(0,0,0)"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:normal;font-family:"Courier New",monospace;vertical-align:baseline;color:inherit"><i style="color:rgb(68,68,68);line-height:normal">Rua da Consolação 896, 01302-907, SP, São Paulo, Brasil.</i></span></div></div></div></div></div>