SPARTA (23 Dec 2017) # flow around a sphere on a uniform grid # particles are created initially and input at boundaries at stream velocity variable x index 10 variable y index 10 variable z index 10 variable n equal 10*$x*$y*$z variable n equal 10*256*$y*$z variable n equal 10*256*160*$z variable n equal 10*256*160*160 variable fnum equal 7.33e+21/v_n seed 12345 dimension 3 global nrho 7.03E18 global fnum ${fnum} global fnum 111846923828125 timestep 1.E-5 global gridcut 0.1 global surfmax 300 boundary o o o create_box -5.0 5.0 -5.0 5.0 -5.0 5.0 Created orthogonal box = (-5 -5 -5) to (5 5 5) create_grid $x $y $z create_grid 256 $y $z create_grid 256 160 $z create_grid 256 160 160 WARNING: Could not acquire nearby ghost cells b/c grid partition is not clumped (../grid.cpp:381) Created 6553600 child grid cells parent cells = 1 CPU time = 0.0334005 secs create/ghost percent = 58.207 41.793 balance_grid rcb cell Balance grid migrated 6508442 cells CPU time = 0.436929 secs reassign/sort/migrate/ghost percent = 39.9011 0.997168 31.6419 27.4598 species ar.species Ar mixture air Ar frac 1.0 mixture air group species mixture air Ar vstream 2500.0 0.0 0.0 temp 300.00 fix in emit/face air all collide vss air ar.vss read_surf data.sphere 602 points 1200 triangles -1 1 xlo xhi -1 1 ylo yhi -1 1 zlo zhi 0.100631 min triangle edge length 0.00439601 min triangle area 6696 = cells with surfs 16924 = total surfs in all grid cells 6 = max surfs in one grid cell 1.6101 = min surf-size/cell-size ratio 0 0 = number of pushed cells 6696 0 = cells overlapping surfs, overlap cells with unmarked corner pts 6523080 23824 6696 = cells outside/inside/overlapping surfs 6696 = surf cells with 1,2,etc splits 995.862 995.862 = cell-wise and global flow volume CPU time = 0.210944 secs read/sort/check/surf2grid/ghost/inout/particle percent = 2.72849 13.0035 1.21321 34.7916 48.2632 33.4064 0.0110835 surf_collide 1 diffuse 300.0 1.0 surf_modify all collide 1 create_particles air n 0 Created 62593655 particles CPU time = 0.150747 secs fix 3 balance 100 1.1 rcb part collide_modify vremax 100 yes #compute 1 grid all all nrho mass u v w #fix 1 ave/grid 10 100 1000 c_1 ave running #dump 1 grid all 1000 tmp_flow.* id xc yc zc f_1 #compute 2 surf all all press ke #fix 2 ave/surf all 10 100 1000 c_2 ave running #dump 2 surf all 1000 tmp_surf.* id f_2 compute g grid all all n compute max reduce max c_g[1] stats_style step cpu np nattempt ncoll nscoll c_max stats 100 # equilibrate flow for 1000 steps # then benchmark for 1000 steps run 1000 WARNING: One or more fix inflow faces oppose streaming velocity (../fix_emit_face.cpp:195) Memory usage per proc in Mbytes: particles (ave,min,max) = 45.375 42.1875 45.5625 grid (ave,min,max) = 11.0089 9.79147 11.6665 surf (ave,min,max) = 0.0778974 0.0778961 0.0778999 total (ave,min,max) = 56.809 54.2791 57.6541 Step CPU Np Natt Ncoll Nscoll c_max 0 0 62593655 0 0 0 11 100 22.776733 62586614 25941 21275 7532 179 200 44.994143 62726151 32200 26683 7703 177 300 66.631362 63244830 34491 28826 7816 202 400 89.368828 63905935 35260 29322 8069 195 500 111.17448 64399917 35457 29512 8087 199 600 133.45732 64657894 35467 29545 7854 187 700 155.88805 64759556 35532 29703 7962 185 800 178.28469 64800911 35672 29730 8042 196 900 200.55186 64816018 35574 29645 7779 194 1000 222.8265 64825943 35529 29633 8023 182 Loop time of 222.829 on 144 procs for 1000 steps with 64825943 particles Particle moves = 64141283307 (64.1B) Cells touched = 108717723977 (109B) Particle comms = 939281984 (939M) Boundary collides = 0 (0K) Boundary exits = 179974359 (180M) SurfColl checks = 468610055 (469M) SurfColl occurs = 7781516 (7.78M) Surf reactions = 0 (0K) Collide attempts = 946050206 (946M) Collide occurs = 719369321 (719M) Reactions = 0 (0K) Particles stuck = 0 Particle-moves/CPUsec/proc: 1.99896e+06 Particle-moves/step: 6.41413e+07 Cell-touches/particle/step: 1.69497 Particle comm iterations/step: 1 Particle fraction communicated: 0.014644 Particle fraction colliding with boundary: 0 Particle fraction exiting boundary: 0.00280591 Surface-checks/particle/step: 0.0073059 Surface-collisions/particle/step: 0.000121318 Surf-reactions/particle/step: 0 Collision-attempts/particle/step: 0.0147495 Collisions/particle/step: 0.0112154 Reactions/particle/step: 0 Move time (%) = 139.413 (62.565) Coll time (%) = 29.6698 (13.315) Sort time (%) = 42.9816 (19.2891) Comm time (%) = 7.55419 (3.39013) Modfy time (%) = 2.70078 (1.21204) Outpt time (%) = 0.504752 (0.22652) Other time (%) = 0.00502221 (0.00225384) Particles: 450180 ave 481514 max 442926 min Histogram: 84 15 11 8 9 5 1 3 4 4 Cells: 45511.1 ave 73160 max 32227 min Histogram: 4 10 29 93 2 2 0 0 0 4 GhostCell: 19550.2 ave 28270 max 10178 min Histogram: 7 4 13 17 25 23 25 15 10 5 EmptyCell: 0 ave 0 max 0 min Histogram: 144 0 0 0 0 0 0 0 0 0 collide_modify vremax 100 no run 1000 WARNING: One or more fix inflow faces oppose streaming velocity (../fix_emit_face.cpp:195) Memory usage per proc in Mbytes: particles (ave,min,max) = 48.1992 47.25 54 grid (ave,min,max) = 11.5046 9.79147 16.9165 surf (ave,min,max) = 0.0778974 0.0778961 0.0778999 total (ave,min,max) = 60.1289 57.4706 68.1775 Step CPU Np Natt Ncoll Nscoll c_max 1000 0 64825943 35529 29633 0 182 1100 22.299004 64830352 35670 29791 8052 191 1200 44.59277 64831337 35666 29622 7843 189 1300 66.862579 64831376 35551 29711 7945 180 1400 89.530279 64833884 35557 29730 7819 185 1500 112.19785 64835886 35659 29757 7895 183 1600 134.74355 64836702 35504 29500 8048 188 1700 157.27141 64831681 35215 29394 8011 192 1800 179.85426 64828840 35528 29781 7945 187 1900 202.4263 64825082 35442 29691 7979 186 2000 224.99076 64825414 35699 29805 7733 183 Loop time of 224.993 on 144 procs for 1000 steps with 64825414 particles Particle moves = 65012848201 (65B) Cells touched = 110119138190 (110B) Particle comms = 928585687 (929M) Boundary collides = 0 (0K) Boundary exits = 182201979 (182M) SurfColl checks = 488614999 (489M) SurfColl occurs = 7987178 (7.99M) Surf reactions = 0 (0K) Collide attempts = 995660366 (996M) Collide occurs = 757266000 (757M) Reactions = 0 (0K) Particles stuck = 0 Particle-moves/CPUsec/proc: 2.00663e+06 Particle-moves/step: 6.50128e+07 Cell-touches/particle/step: 1.69381 Particle comm iterations/step: 1 Particle fraction communicated: 0.0142831 Particle fraction colliding with boundary: 0 Particle fraction exiting boundary: 0.00280255 Surface-checks/particle/step: 0.00751567 Surface-collisions/particle/step: 0.000122855 Surf-reactions/particle/step: 0 Collision-attempts/particle/step: 0.0153148 Collisions/particle/step: 0.0116479 Reactions/particle/step: 0 Move time (%) = 142.866 (63.4981) Coll time (%) = 30.078 (13.3684) Sort time (%) = 42.97 (19.0984) Comm time (%) = 7.87038 (3.49806) Modfy time (%) = 0.67816 (0.301414) Outpt time (%) = 0.525432 (0.233533) Other time (%) = 0.00469419 (0.00208637) Particles: 450176 ave 481428 max 442808 min Histogram: 87 12 11 6 8 8 3 2 3 4 Cells: 45511.1 ave 73160 max 32227 min Histogram: 4 10 29 93 2 2 0 0 0 4 GhostCell: 19550.2 ave 28270 max 10178 min Histogram: 7 4 13 17 25 23 25 15 10 5 EmptyCell: 0 ave 0 max 0 min Histogram: 144 0 0 0 0 0 0 0 0 0