[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