Initializing AMReX (26.10-139-g01130277b04d)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 OMP initialized with 1 OMP threads AMReX (26.10-139-g01130277b04d) initialized MFIX git describe: 26.10-45-g517829837eb0 AMReX-Hydro git hash: 154572aa6cf22c35a69729c74a0276909b96d5e3 Number of levels: 1 RunType is Standard Auto-generating level-set parameters: eb_refinement = 1 levelset_pad = 2 eb_pad = 2 ================================================================================ Warning: Overwriting amr.check_file with value from mfix.check_file ================================================================================ ================================================================================ Warning: Overwriting amr.plt_regtest with value from mfix.plt_regtest ================================================================================ ================================================================================ Warning: Overwriting amr.plot_int with value from output.plot.int ================================================================================ ================================================================================ Warning: Overwriting amr.plot_file with value from output.plot.file ================================================================================ Simulation is transient -------------------------------------------------------------------------------- Info: Max effective viscosity factor not specified. Using default 1e+06 -------------------------------------------------------------------------------- Particle-particle conduction is disabled. Particle-fluid-particle conduction is disabled. Load balance setup: Type: DualGrid Strategy: KnapSack Weighting: ParticleCount Drag model: BVK2 Convection model not specified. Setting default: mfix.convection.model = RanzMarshall Building cylinder geometry. Internal Flow: 1 Radius: 0.0039 Height: -1 Direction: 0 Rotation angle: 0 Rotation axis: 0 Center: 0.004, 0.004, 0.004 Done making the fluid eb levels ... On level 0 updated BC areas: BC: 0 area: 4.776607274e-05 BC: 1 area: 4.776607274e-05 ================================================================================ Warning: Setting particle max grid size to fluid max grid size ================================================================================ In MakeBaseGrids: BA HAS 16 GRIDS Auto generating particles ... Total number of generated particles: 215,120 Time spent in initializing particles 0.805695796 **************************************** Coverage report: Grids Cells regular: 4 44960 covered: 0 12800 cut: 12 6240 **************************************** Bounding box for levelset imposed walls: Lo: 1e-15 -0.008 -0.008 Hi: 0.016 0.016 0.016 Clean up auto-generated particles. IC Region Generated Removed Remaining **************** **************** **************** **************** full-domain 0 0 0 bed 215,120 58,061 157,059 ---------------- ---------------- ---------------- ---------------- Total 215,120 58,061 157,059 **************** **************** **************** **************** Setting original sum_vol to 2.998926818e-07 Initial projection: Before projection: max(abs(u/v/w)) = 0.015 0 0 max(abs(gpx/gpy/gpz)) = 0 0 0 Min and Max of epf 0.2712117833 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 44.6112468 MLMG: Initial rhs = 22.3056234 MLMG: Initial residual (resid0) = 22.3056234 MLMG: Final Iter. 12 resid, resid/bnorm = 1.475768396e-10, 6.616127108e-12 MLMG: Timers: Solve = 0.149024345 Iter = 0.11377801 Bottom = 0.003637103 After projection: max(abs(u/v/w)) = 0.02199251829 0.001131161343 0.001189525565 Doing initial pressure iterations with dt = 0.0002 In initial_iterations: iter = 0 Before projection: max(abs(u/v/w)) = 0.014211585 0.004882459189 0.004977002289 max(abs(gpx/gpy/gpz)) = 0 0 0 Min and Max of epf 0.2712117833 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 148.002174 MLMG: Initial rhs = 74.00108702 MLMG: Initial residual (resid0) = 74.00108702 MLMG: Final Iter. 12 resid, resid/bnorm = 6.681746545e-10, 9.029254588e-12 MLMG: Timers: Solve = 0.147665793 Iter = 0.113768103 Bottom = 0.003814108 After projection: max(abs(u/v/w)) = 0.02182752578 0.002007162363 0.002102075339 In initial_iterations: iter = 1 Before projection: max(abs(u/v/w)) = 0.01715471639 0.004871918643 0.004962549299 max(abs(gpx/gpy/gpz)) = 7.863898951 0.6377384995 0.6836572841 Min and Max of epf 0.2712117833 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 147.5196537 MLMG: Initial rhs = 73.75982684 MLMG: Initial residual (resid0) = 73.75982684 MLMG: Final Iter. 12 resid, resid/bnorm = 5.329256481e-10, 7.225147765e-12 MLMG: Timers: Solve = 0.147526141 Iter = 0.114058145 Bottom = 0.003841185 After projection: max(abs(u/v/w)) = 0.02177113833 0.002773963131 0.002889364751 In initial_iterations: iter = 2 Before projection: max(abs(u/v/w)) = 0.01774292217 0.004861723965 0.00494825502 max(abs(gpx/gpy/gpz)) = 15.69527858 1.377175469 1.446111575 Min and Max of epf 0.2712117833 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 147.0385839 MLMG: Initial rhs = 73.51929195 MLMG: Initial residual (resid0) = 73.51929195 MLMG: Final Iter. 12 resid, resid/bnorm = 2.257256049e-10, 3.070290789e-12 MLMG: Timers: Solve = 0.14780659 Iter = 0.114135454 Bottom = 0.003879743 After projection: max(abs(u/v/w)) = 0.0217799692 0.003183185454 0.003316757802 Enclosed domain volume is 2.998926818e-07 Time spent in init 3.904140455 Writing plotfile BENCH05-Size0008_wide-Chemistry_plt00000 at time 0 Regridding at step 0 Updating particle ebfactory Regridding level-set on lev = 0 Also regridding refined level-set On level 0 updated BC areas: BC: 0 area: 4.776607274e-05 BC: 1 area: 4.776607274e-05 ============ NEW TIME STEP ============ Step 1: from old_time 0 to new time 1.807175116e-05 with dt = 1.807175116e-05 Before projection (with small dt modification): max(abs(u/v/w)) = 0.01505648261 0.004738505028 0.00475849595 max(abs(gpx/gpy/gpz)) = 23.4662823 2.128787563 2.205468023 Min and Max of epf 0.2712117833 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 182.2931523 MLMG: Initial rhs = 91.14657615 MLMG: Initial residual (resid0) = 91.14657615 MLMG: Final Iter. 12 resid, resid/bnorm = 7.251361733e-10, 7.955714893e-12 MLMG: Timers: Solve = 0.147834456 Iter = 0.113782161 Bottom = 0.003657836 After projection (with small dt modification): max(abs(u/v/w)) = 0.02281090062 0.002039600208 0.002315378864 Evolving particles on level: 0 ... with fluid dt 1.807175116e-05 done. Time per fluid step 0.811805211 Time per 38 particle steps 0.649161134 Coupling time per step 0.067655891 Time per chemistry step 0.028219972 Time per step 1.61473481 Regridding at step 1 On level 0 updated BC areas: BC: 0 area: 4.776607274e-05 BC: 1 area: 4.776607274e-05 ============ NEW TIME STEP ============ Step 2: from old_time 1.807175116e-05 to new time 3.604288715e-05 with dt = 1.797113599e-05 Before projection: max(abs(u/v/w)) = 0.02681226687 0.005836315522 0.005691195537 max(abs(gpx/gpy/gpz)) = 107.9543545 8.374418462 9.273648694 Min and Max of epf 0.2712174435 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 318.5347437 MLMG: Initial rhs = 159.2673719 MLMG: Initial residual (resid0) = 159.2673719 MLMG: Final Iter. 12 resid, resid/bnorm = 1.239833236e-09, 7.784602845e-12 MLMG: Timers: Solve = 0.147050942 Iter = 0.113619768 Bottom = 0.003829542 After projection: max(abs(u/v/w)) = 0.02524924737 0.003281428983 0.003695600231 Evolving particles on level: 0 ... with fluid dt 1.797113599e-05 done. Time per fluid step 0.828344415 Time per 38 particle steps 0.644713252 Coupling time per step 0.065142767 Time per chemistry step 0.028311869 Time per step 1.576408875 Regridding at step 2 On level 0 updated BC areas: BC: 0 area: 4.776607274e-05 BC: 1 area: 4.776607274e-05 ============ NEW TIME STEP ============ Step 3: from old_time 3.604288715e-05 to new time 5.391431677e-05 with dt = 1.787142962e-05 Before projection: max(abs(u/v/w)) = 0.04753113255 0.006581978554 0.006336075521 max(abs(gpx/gpy/gpz)) = 192.2816315 14.05904283 15.17102024 Min and Max of epf 0.2712339391 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 452.0718766 MLMG: Initial rhs = 226.0359383 MLMG: Initial residual (resid0) = 226.0359383 MLMG: Final Iter. 12 resid, resid/bnorm = 1.424811547e-09, 6.303473498e-12 MLMG: Timers: Solve = 0.148395278 Iter = 0.114483003 Bottom = 0.003702879 After projection: max(abs(u/v/w)) = 0.02768826314 0.004167676031 0.004629295839 Evolving particles on level: 0 ... with fluid dt 1.787142962e-05 done. Time per fluid step 0.828045125 Time per 37 particle steps 0.63399807 Coupling time per step 0.065361456 Time per chemistry step 0.02777389 Time per step 1.565035983 Regridding at step 3 On level 0 updated BC areas: BC: 0 area: 4.776607274e-05 BC: 1 area: 4.776607274e-05 ============ NEW TIME STEP ============ Step 4: from old_time 5.391431677e-05 to new time 7.16863061e-05 with dt = 1.777198933e-05 Before projection: max(abs(u/v/w)) = 0.06759605467 0.007374890746 0.006990073112 max(abs(gpx/gpy/gpz)) = 275.5045065 21.1020983 21.45035686 Min and Max of epf 0.2712609676 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 581.8752264 MLMG: Initial rhs = 290.9376132 MLMG: Initial residual (resid0) = 290.9376132 MLMG: Final Iter. 12 resid, resid/bnorm = 1.067606559e-09, 3.669537764e-12 MLMG: Timers: Solve = 0.148002014 Iter = 0.114581164 Bottom = 0.003710244 After projection: max(abs(u/v/w)) = 0.03007612074 0.004888809133 0.00537161067 Evolving particles on level: 0 ... with fluid dt 1.777198933e-05 done. Time per fluid step 0.826240639 Time per 37 particle steps 0.633209892 Coupling time per step 0.06595344 Time per chemistry step 0.027791993 Time per step 1.562897392 Regridding at step 4 On level 0 updated BC areas: BC: 0 area: 4.776607274e-05 BC: 1 area: 4.776607274e-05 ============ NEW TIME STEP ============ Step 5: from old_time 7.16863061e-05 to new time 8.935580204e-05 with dt = 1.766949594e-05 Before projection: max(abs(u/v/w)) = 0.0870940195 0.008118231492 0.007598569142 max(abs(gpx/gpy/gpz)) = 357.7096866 28.59431563 28.9476102 Min and Max of epf 0.2712982317 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 707.8484696 MLMG: Initial rhs = 353.9242348 MLMG: Initial residual (resid0) = 353.9242348 MLMG: Final Iter. 12 resid, resid/bnorm = 1.194579768e-09, 3.375241509e-12 MLMG: Timers: Solve = 0.147673514 Iter = 0.113814124 Bottom = 0.003696619 After projection: max(abs(u/v/w)) = 0.0324114816 0.005388470709 0.006029773269 Evolving particles on level: 0 ... with fluid dt 1.766949594e-05 done. Time per fluid step 0.825033971 Time per 37 particle steps 0.631608861 Coupling time per step 0.064500033 Time per chemistry step 0.027896054 Time per step 1.558842092 Regridding at step 5 On level 0 updated BC areas: BC: 0 area: 4.776607274e-05 BC: 1 area: 4.776607274e-05 ============ NEW TIME STEP ============ Step 6: from old_time 8.935580204e-05 to new time 0.0001069226982 with dt = 1.756689616e-05 Before projection: max(abs(u/v/w)) = 0.106063641 0.008833618049 0.008747788134 max(abs(gpx/gpy/gpz)) = 438.9039291 36.39268723 37.40026379 Min and Max of epf 0.2713454319 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 830.035934 MLMG: Initial rhs = 415.017967 MLMG: Initial residual (resid0) = 415.017967 MLMG: Final Iter. 12 resid, resid/bnorm = 1.421181395e-09, 3.424385226e-12 MLMG: Timers: Solve = 0.148202752 Iter = 0.114429168 Bottom = 0.003591215 After projection: max(abs(u/v/w)) = 0.0346951635 0.005790911664 0.006514597912 Evolving particles on level: 0 ... with fluid dt 1.756689616e-05 done. Time per fluid step 0.825801451 Time per 37 particle steps 0.632723463 Coupling time per step 0.065779767 Time per chemistry step 0.028275922 Time per step 1.562530228 Regridding at step 6 On level 0 updated BC areas: BC: 0 area: 4.776607274e-05 BC: 1 area: 4.776607274e-05 ============ NEW TIME STEP ============ Step 7: from old_time 0.0001069226982 to new time 0.000124387851 with dt = 1.746515278e-05 Before projection: max(abs(u/v/w)) = 0.124528642 0.009978864224 0.0107785246 max(abs(gpx/gpy/gpz)) = 519.0669236 44.4249925 46.31407092 Min and Max of epf 0.2714022764 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 948.2954258 MLMG: Initial rhs = 474.1477129 MLMG: Initial residual (resid0) = 474.1477129 MLMG: Final Iter. 12 resid, resid/bnorm = 2.815490063e-09, 5.938001991e-12 MLMG: Timers: Solve = 0.150236221 Iter = 0.115687441 Bottom = 0.00368973 After projection: max(abs(u/v/w)) = 0.03692823426 0.00608817008 0.0068872866 Evolving particles on level: 0 ... with fluid dt 1.746515278e-05 done. Time per fluid step 0.827382743 Time per 36 particle steps 0.62343415 Coupling time per step 0.064569859 Time per chemistry step 0.027851169 Time per step 1.552957852 Regridding at step 7 On level 0 updated BC areas: BC: 0 area: 4.776607274e-05 BC: 1 area: 4.776607274e-05 ============ NEW TIME STEP ============ Step 8: from old_time 0.000124387851 to new time 0.000141752461 with dt = 1.736461007e-05 Before projection: max(abs(u/v/w)) = 0.1425041963 0.0117948001 0.01288046407 max(abs(gpx/gpy/gpz)) = 598.1791447 52.62228142 55.63630627 Min and Max of epf 0.2714684852 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 1062.678076 MLMG: Initial rhs = 531.3390379 MLMG: Initial residual (resid0) = 531.3390379 MLMG: Final Iter. 12 resid, resid/bnorm = 2.1388189e-09, 4.025337397e-12 MLMG: Timers: Solve = 0.149648553 Iter = 0.115116864 Bottom = 0.003581956 After projection: max(abs(u/v/w)) = 0.03911162421 0.006321537108 0.007191503641 Evolving particles on level: 0 ... with fluid dt 1.736461007e-05 done. Time per fluid step 0.833199966 Time per 36 particle steps 0.61945249 Coupling time per step 0.065899692 Time per chemistry step 0.028116956 Time per step 1.555962887 Regridding at step 8 On level 0 updated BC areas: BC: 0 area: 4.776607274e-05 BC: 1 area: 4.776607274e-05 ============ NEW TIME STEP ============ Step 9: from old_time 0.000141752461 to new time 0.0001590181054 with dt = 1.726564433e-05 Before projection: max(abs(u/v/w)) = 0.160004503 0.01363416457 0.01504175748 max(abs(gpx/gpy/gpz)) = 676.2218055 60.94418239 65.31987066 Min and Max of epf 0.271543779 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 1173.276694 MLMG: Initial rhs = 586.6383472 MLMG: Initial residual (resid0) = 586.6383472 MLMG: Final Iter. 12 resid, resid/bnorm = 2.22436114e-09, 3.791707702e-12 MLMG: Timers: Solve = 0.148073799 Iter = 0.114177038 Bottom = 0.0036223 After projection: max(abs(u/v/w)) = 0.04124602932 0.006513875686 0.007451450232 Evolving particles on level: 0 ... with fluid dt 1.726564433e-05 done. Time per fluid step 0.82627895 Time per 36 particle steps 0.620027413 Coupling time per step 0.064550595 Time per chemistry step 0.028311755 Time per step 1.549530577 Regridding at step 9 On level 0 updated BC areas: BC: 0 area: 4.776607274e-05 BC: 1 area: 4.776607274e-05 ============ NEW TIME STEP ============ Step 10: from old_time 0.0001590181054 to new time 0.0001761866629 with dt = 1.716855748e-05 Before projection: max(abs(u/v/w)) = 0.1770420254 0.01548962913 0.01725173592 max(abs(gpx/gpy/gpz)) = 753.1789677 69.35410345 75.32084892 Min and Max of epf 0.2716278916 1 Nodal Projection: >> Before projection: * On lev 0 max(abs(rhs)) = 1280.194078 MLMG: Initial rhs = 640.0970389 MLMG: Initial residual (resid0) = 640.0970389 MLMG: Final Iter. 12 resid, resid/bnorm = 3.436030127e-09, 5.367983163e-12 MLMG: Timers: Solve = 0.147151002 Iter = 0.11378565 Bottom = 0.003822899 After projection: max(abs(u/v/w)) = 0.04333221318 0.006679636746 0.00768215459 Evolving particles on level: 0 ... with fluid dt 1.716855748e-05 done. Time per fluid step 0.824788277 Time per 36 particle steps 0.635654821 Coupling time per step 0.076508471 Time per chemistry step 0.028305593 Time per step 1.564447304 ************************************************************************************************************** Species mass balance for interval: 0.000000 to 0.000176 Species mass(t+dt) mass(t) production mass(in) mass(out) net accu flux CH4 2.1642e-09 2.1642e-09 0.0000e+00 8.9481e-13 9.0856e-13 -4.1359e-25 -1.3748e-14 H2 4.3284e-09 4.3284e-09 0.0000e+00 1.7896e-12 1.8171e-12 -8.2718e-25 -2.7497e-14 CO 2.1642e-09 2.1642e-09 0.0000e+00 8.9481e-13 9.0856e-13 -4.1359e-25 -1.3748e-14 H2O 4.3284e-09 4.3284e-09 0.0000e+00 1.7896e-12 1.8171e-12 -8.2718e-25 -2.7497e-14 CO2 2.1642e-09 2.1642e-09 0.0000e+00 8.9481e-13 9.0856e-13 -4.1359e-25 -1.3748e-14 O2 2.1642e-09 2.1642e-09 0.0000e+00 8.9481e-13 9.0856e-13 -4.1359e-25 -1.3748e-14 N2 4.3284e-09 4.3284e-09 0.0000e+00 1.7896e-12 1.8171e-12 -8.2718e-25 -2.7497e-14 Fe2O3 4.9342e-05 4.9342e-05 0.0000e+00 0.0000e+00 0.0000e+00 6.7763e-21 0.0000e+00 FeO 4.9342e-05 4.9342e-05 0.0000e+00 0.0000e+00 0.0000e+00 6.7763e-21 0.0000e+00 inert 6.5789e-05 6.5789e-05 0.0000e+00 0.0000e+00 0.0000e+00 -1.3553e-20 0.0000e+00 -------------------------------------------------------------------------------------------------------------- Total 1.6449e-04 1.6449e-04 0.0000e+00 8.9481e-12 9.0856e-12 -4.1359e-24 -1.3748e-13 total number of particles = 157,059 |(total acc) - (net flux)| = 1.3748e-13 ************************************************************************************************************** Writing plotfile BENCH05-Size0008_wide-Chemistry_plt00010 at time 0.0001761866629 Time spent in main (after init) 15.70591534 Time spent in main 19.6100558 AMReX (26.10-139-g01130277b04d) finalized