[FLASH-USERS] Sedov explosion in 3D reaches max simTime very fast without actually simulating the explosion
Ryan Farber
rjfarber at umich.edu
Wed Jul 9 18:56:50 EDT 2025
Hi Sumit,
The parameter names are zmin and zmax, not z_min and z_max. Maybe that’s one issue.
Anyway, from what you shared there’s just a single block (no refinement occurred). The question is why. It’d be helpful if you shared your.log file.
If you’re just trying to run the default Sedov problem in 3D with only changes in flash.par I’m a bit surprised you have this issue. Though as you wrote, setting lrefine_min can be helpful if you aren’t initially resolving the features you’re trying to refine on.
Best wishes,
Ryan
Sent from my iPhone
> On 9. Jul 2025, at 15:43, Sumit Banerjee <sumitb at g.clemson.edu> wrote:
>
>
> 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)
> _______________________________________________
> flash-users mailing list
> flash-users at flash.rochester.edu
>
> For list info, including unsubscribe:
> https://flash.rochester.edu/mailman/listinfo/flash-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20250709/2f227174/attachment-0001.htm>
More information about the flash-users
mailing list