
<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class=""><br class=""></div><div class="">Hi Doris, Hi Aaron,</div><div class=""><br class=""></div><div class="">from my experience with this, it might not help much anyway to reset smlrho. The actual problem is that the PPM solver crashes, because it is not very robust (I’ve had a lot of these non-convergence problems in the past when I was still using PPM). You can try to set the CFL number smaller and/or nriem higher to push it further, but eventually, you’d need a more diffusive solver (e.g., HLLE) or (actually better) a stable solver that keeps density and pressure positive by construction (<a href="http://adsabs.harvard.edu/abs/2011JCoPh.230.3331W" class="">http://adsabs.harvard.edu/abs/2011JCoPh.230.3331W</a>). The latter is not in the public version, but if you’d like to collaborate, I can send the solver to you. It is an MHD solver, but of course, you can use it with B=0 and I also have a stripped pure-hydro version of it.</div><div class=""><br class=""></div><div class="">Also, since this is a collapse, eventually you will need to create a sink particle in the centre to cut off the collapse in a controlled way in order to follow what’s happening on larger scales (see sink particle technique: <a href="http://adsabs.harvard.edu/abs/2010ApJ...713..269F" class="">http://adsabs.harvard.edu/abs/2010ApJ...713..269F</a>). The sink particle method is already in the public version, so feel free to switch it on and see if you can get further. The way you’d use the sink particles is explained in the user guide. Depending on what your science goal with this collapse simulation is, this is certainly something to consider in any case.</div><div class=""><br class=""></div><div class="">I hope this helps.</div><div class=""><br class=""></div><div class="">Kind regards,</div><div class="">Christoph</div><br class=""><div class=""><blockquote type="cite" class=""><div class="">On 24 Mar 2016, at 11:26, Jung Lin Doris Lee <<a href="mailto:dorislee@berkeley.edu" class="">dorislee@berkeley.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><p class=""><span class="">Flash community,</span></p><p class="">We have a (hopefully simple) question about how to best change the numerical floors for state quantities like density and pressure. For context, we are trying to reproduce the cloud collapse problem of Foster & Chevalier (<a href="http://adsabs.harvard.edu/abs/1993ApJ...416..303F" class=""><span class="">ADS</span></a>), which require densities of 10^-19~10^-27g/cm^3. The setup has a spherical core in the central ~1/2 of the domain with a radial profile in density. Outside the sphere is a tenuous and hotter medium (a reduced density but increased temperature so that it remains in pressure balance with the edge of the sphere). <br class=""><span class=""></span></p><p class="">For example, when we initialize the data, we receive the following warning:<br class=""><span class=""></span></p><p class=""><span class="">WARNING after gc filling: min. unk(DENS_VAR)=0.7259099999999999388370E-20 PE=0 block=10 type=1</span></p><p class=""><span class="">And the initial outputs show that the density has become entirely uniform after which the program crashes because of non-convergence in the Riemann subroutine. Grep-ing around, we see that there are quantities like “smallp” or “smlrho”, whose default values are 10^-10g/cm^3, indeed larger than the values we are initializing. We tried setting a value for smlrho in flash.par, but it did not appear to do anything.</span></p><p class=""><span class="">Can these values be reset in flash.par versus going into the source code? We are glancing through the user guide and example problems and can not find an example of how to do this. </span></p><div class=""><span class=""></span><br class="webkit-block-placeholder"></div><p class=""><span class="">Best,</span></p><p class=""><span class="">Doris and Aaron</span></p></div>

</div></blockquote></div><br class=""></body></html>