Initializing AMReX (26.10-131-g5a9a0cc7f8d8)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 OMP initialized with 1 OMP threads AMReX (26.10-131-g5a9a0cc7f8d8) initialized MFIX git describe: 26.10-44-gd5d5db38d15a 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.793556972 **************************************** 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.148277319 Iter = 0.113442131 Bottom = 0.003802509 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.147183729 Iter = 0.113889928 Bottom = 0.00387582 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.147180125 Iter = 0.113685023 Bottom = 0.003861002 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.146462153 Iter = 0.113177454 Bottom = 0.003705515 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.883819324 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.147141409 Iter = 0.113766667 Bottom = 0.003740536 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.809356611 Time per 38 particle steps 0.65298213 Coupling time per step 0.070182208 Time per chemistry step 0.028456835 Time per step 1.61562421 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.146959949 Iter = 0.113533036 Bottom = 0.003973541 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.828369081 Time per 38 particle steps 0.648876836 Coupling time per step 0.064338242 Time per chemistry step 0.028497782 Time per step 1.580034609 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.147402977 Iter = 0.113995912 Bottom = 0.003723176 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.826335051 Time per 37 particle steps 0.635320881 Coupling time per step 0.066098397 Time per chemistry step 0.02840568 Time per step 1.565481479 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.147511696 Iter = 0.114112557 Bottom = 0.003674026 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.822929498 Time per 37 particle steps 0.6556551 Coupling time per step 0.064708702 Time per chemistry step 0.028337706 Time per step 1.581194894 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.149323973 Iter = 0.115644828 Bottom = 0.004339289 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.836538986 Time per 37 particle steps 0.637276334 Coupling time per step 0.064809335 Time per chemistry step 0.028005813 Time per step 1.576281067 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.146990329 Iter = 0.113520095 Bottom = 0.00374314 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.825060407 Time per 37 particle steps 0.638408114 Coupling time per step 0.064930401 Time per chemistry step 0.02838747 Time per step 1.566344024 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.14701659 Iter = 0.113477249 Bottom = 0.003626053 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.825018505 Time per 36 particle steps 0.623094407 Coupling time per step 0.064288615 Time per chemistry step 0.028507724 Time per step 1.550780093 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.147141196 Iter = 0.113800798 Bottom = 0.00392972 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.82405389 Time per 36 particle steps 0.622518183 Coupling time per step 0.066713341 Time per chemistry step 0.027927908 Time per step 1.550823802 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.147210629 Iter = 0.113840213 Bottom = 0.003833658 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.824714551 Time per 36 particle steps 0.622045499 Coupling time per step 0.064921723 Time per chemistry step 0.027978262 Time per step 1.549138966 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.147730845 Iter = 0.113794236 Bottom = 0.003776801 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.823997867 Time per 36 particle steps 0.622058144 Coupling time per step 0.064647175 Time per chemistry step 0.028156473 Time per step 1.548527715 ************************************************************************************************************** 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.72622973 Time spent in main 19.61004905 AMReX (26.10-131-g5a9a0cc7f8d8) finalized