Initializing AMReX (26.11-5-g1d9bd90c3583)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 OMP initialized with 1 OMP threads AMReX (26.11-5-g1d9bd90c3583) initialized MFIX git describe: 26.11-11-gbf4ab648cdd0 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.82074871 **************************************** 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.147901241 Iter = 0.113046659 Bottom = 0.003786977 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.146817867 Iter = 0.11353343 Bottom = 0.003626505 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.146729222 Iter = 0.113231659 Bottom = 0.003791784 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.146702609 Iter = 0.113449389 Bottom = 0.004032919 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.912654526 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.253280199e-10, 7.957819707e-12 MLMG: Timers: Solve = 0.147298462 Iter = 0.114001092 Bottom = 0.003638137 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.80906247 Time per 38 particle steps 0.654589654 Coupling time per step 0.069320456 Time per chemistry step 0.064798422 Time per step 1.653757651 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.238944503e-09, 7.779022709e-12 MLMG: Timers: Solve = 0.146844403 Iter = 0.113546136 Bottom = 0.003610943 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.82533078 Time per 38 particle steps 0.650203151 Coupling time per step 0.064783543 Time per chemistry step 0.065318515 Time per step 1.616124311 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.430390972e-09, 6.328157298e-12 MLMG: Timers: Solve = 0.146745682 Iter = 0.113407665 Bottom = 0.003714388 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.826340699 Time per 37 particle steps 0.636801773 Coupling time per step 0.063867436 Time per chemistry step 0.066076863 Time per step 1.603059774 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.065678545e-09, 3.662910868e-12 MLMG: Timers: Solve = 0.147408131 Iter = 0.113695365 Bottom = 0.003808721 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.826140659 Time per 37 particle steps 0.636976557 Coupling time per step 0.064522567 Time per chemistry step 0.065373735 Time per step 1.602613642 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.198314337e-09, 3.385793395e-12 MLMG: Timers: Solve = 0.147437874 Iter = 0.113898636 Bottom = 0.004023921 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.826990501 Time per 37 particle steps 0.637403809 Coupling time per step 0.064725421 Time per chemistry step 0.066127358 Time per step 1.604682749 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.415784823e-09, 3.411382001e-12 MLMG: Timers: Solve = 0.146836094 Iter = 0.113526832 Bottom = 0.003778409 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.827620418 Time per 37 particle steps 0.635712416 Coupling time per step 0.064282215 Time per chemistry step 0.066086877 Time per step 1.603421102 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.802867272e-09, 5.911379925e-12 MLMG: Timers: Solve = 0.147283961 Iter = 0.11391899 Bottom = 0.003884016 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.828598544 Time per 36 particle steps 0.624881971 Coupling time per step 0.063903693 Time per chemistry step 0.066010323 Time per step 1.593661263 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.236156149e-09, 4.208529751e-12 MLMG: Timers: Solve = 0.147465044 Iter = 0.11405364 Bottom = 0.003595983 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.84444537 Time per 36 particle steps 0.624067309 Coupling time per step 0.06477171 Time per chemistry step 0.065499769 Time per step 1.6085814 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.089734608e-09, 3.562219581e-12 MLMG: Timers: Solve = 0.147383384 Iter = 0.113955993 Bottom = 0.003578526 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.826234545 Time per 36 particle steps 0.6246303 Coupling time per step 0.063888947 Time per chemistry step 0.066105531 Time per step 1.590706038 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.769907053e-09, 5.889586772e-12 MLMG: Timers: Solve = 0.147502782 Iter = 0.113979328 Bottom = 0.003648775 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.827840819 Time per 36 particle steps 0.621550796 Coupling time per step 0.06567387 Time per chemistry step 0.065569084 Time per step 1.590324248 ************************************************************************************************************** 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) 16.10633124 Time spent in main 20.01898577 AMReX (26.11-5-g1d9bd90c3583) finalized