[FLASH-USERS] Sedov explosion in 3D reaches max simTime very fast without actually simulating the explosion

Sumit Banerjee sumitb at g.clemson.edu
Wed Jul 9 18:42:34 EDT 2025


Flash users,

I hope everything is fine.

I am trying to run the Sedov explosion in 3D and I changed the flash.par by
adding the parameters, z_min, z_max, zl_boundary_type and zr_boundary_type.
However, when I run the simulation, it reaches the max simTime real quick
and the output files do not show any explosion. Below is the output of the
run..

 Grid_init: resolution based on runtime params:
>   lrefine          dx              dy              dz
>       1             0.125           0.125           0.125
>       2             0.062           0.062           0.062
>       3             0.031           0.031           0.031
>       4             0.016           0.016           0.016
>       5             0.008           0.008           0.008
>       6             0.004           0.004           0.004
>  MaterialProperties initialized
>  Cosmology initialized
>  Source terms initialized
>  sim_rhoAmbient is   1.0000000000000000
>   iteration, no. not moved =            0           0
>  refined: total leaf blocks =            1
>  refined: total blocks =            1
>  [amr_morton_process]: Initializing surr_blks using standard orrery
> implementation
>   Finished with Grid_initDomain, no restart
>  Ready to call Hydro_init
>  Hydro initialized
>  Gravity initialized
>  Initial dt verified
>  *** Wrote checkpoint file to sedov_hdf5_chk_0000 ****
>  *** Wrote plotfile to sedov_hdf5_plt_cnt_0000 ****
>  Initial plotfile written
>  Driver init all done
>        n          t         dt  (         x,          y,          z) |
>  dt_hydro  dt_Diff
>        1 1.0000E-10 2.0000E-10  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>        2 3.0000E-10 4.0000E-10  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>        3 7.0000E-10 8.0000E-10  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>        4 1.5000E-09 1.6000E-09  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>        5 3.1000E-09 3.2000E-09  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>        6 6.3000E-09 6.4000E-09  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>        7 1.2700E-08 1.2800E-08  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>        8 2.5500E-08 2.5600E-08  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>        9 5.1100E-08 5.1200E-08  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       10 1.0230E-07 1.0240E-07  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       11 2.0470E-07 2.0480E-07  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       12 4.0950E-07 4.0960E-07  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       13 8.1910E-07 8.1920E-07  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       14 1.6383E-06 1.6384E-06  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       15 3.2767E-06 3.2768E-06  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       16 6.5535E-06 6.5536E-06  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       17 1.3107E-05 1.3107E-05  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       18 2.6214E-05 2.6214E-05  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       19 5.2429E-05 5.2429E-05  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       20 1.0486E-04 1.0486E-04  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       21 2.0972E-04 2.0972E-04  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       22 4.1943E-04 4.1943E-04  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       23 8.3886E-04 8.3886E-04  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       24 1.6777E-03 1.6777E-03  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       25 3.3554E-03 3.3554E-03  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       26 6.7109E-03 6.7109E-03  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       27 1.3422E-02 1.3422E-02  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       28 2.6844E-02 2.6844E-02  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       29 5.3687E-02 5.3687E-02  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>       30 1.0737E-01 1.0737E-01  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>  *** Wrote checkpoint file to sedov_hdf5_chk_0001 ****
>       31 2.1475E-01 2.1475E-01  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>  *** Wrote checkpoint file to sedov_hdf5_chk_0002 ****
>       32 4.2950E-01 4.2950E-01  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>  *** Wrote checkpoint file to sedov_hdf5_chk_0003 ****
>       33 8.5899E-01 8.5899E-01  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>  *** Wrote checkpoint file to sedov_hdf5_chk_0004 ****
>       34 1.7180E+00 1.7180E+00  ( 6.250E-02,  6.250E-02,  6.250E-02) |
>  2.673E+01 6.250E+27
>  *** Wrote checkpoint file to sedov_hdf5_chk_0005 ****
>  exiting: reached max SimTime
>  *** Wrote plotfile to sedov_forced_hdf5_plt_cnt_0000 ****


I tried reading through the past issues and found out a suggestion to
increase the parameters lrefine_min or lrefine_min_init but they aren't
there in the flash.par file under the
source/Simulation/SimulationMain/Sedov/flash.par file. I tried checking a
few par files written for 3D runs but couldn't find anything useful yet.

I tried running the Sod problem in 3D and it runs successfully with the
similar adjustments to the par file as above so I have no idea what's going
on here.

Could somebody explain what is going on and how I can get over this
problem?

*Clear Skies,*
*Sumit (Sam)*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20250709/603e8fb2/attachment.htm>


More information about the flash-users mailing list