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*64*$y*$z variable n equal 10*64*40*$z variable n equal 10*64*40*40 variable fnum equal 7.33e+21/v_n seed 12345 dimension 3 global nrho 7.03E18 global fnum ${fnum} global fnum 7.158203125e+15 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 64 $y $z create_grid 64 40 $z create_grid 64 40 40 WARNING: Could not acquire nearby ghost cells b/c grid partition is not clumped (../grid.cpp:381) Created 102400 child grid cells parent cells = 1 CPU time = 0.0152829 secs create/ghost percent = 48.7278 51.2722 balance_grid rcb cell Balance grid migrated 101600 cells CPU time = 0.0794759 secs reassign/sort/migrate/ghost percent = 59.1815 2.98459 21.9574 15.8766 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 424 = cells with surfs 4096 = total surfs in all grid cells 19 = max surfs in one grid cell 0.402525 = min surf-size/cell-size ratio 0 0 = number of pushed cells 424 0 = cells overlapping surfs, overlap cells with unmarked corner pts 101744 232 424 = cells outside/inside/overlapping surfs 424 = surf cells with 1,2,etc splits 995.862 995.862 = cell-wise and global flow volume CPU time = 0.062232 secs read/sort/check/surf2grid/ghost/inout/particle percent = 29.5767 4.49927 12.6929 42.6032 10.6279 18.6993 0.0758563 surf_collide 1 diffuse 300.0 1.0 surf_modify all collide 1 create_particles air n 0 Created 978025 particles CPU time = 0.0172389 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) = 1.6875 1.6875 1.6875 grid (ave,min,max) = 1.66647 1.66647 1.66647 surf (ave,min,max) = 0.0779014 0.0778999 0.0779037 total (ave,min,max) = 3.43797 3.43797 3.43797 Step CPU Np Natt Ncoll Nscoll c_max 0 0 978025 0 0 0 11 100 0.33666801 977328 342 304 118 79 200 0.70951009 979581 443 381 118 90 300 1.0519819 987941 467 419 134 89 400 1.4266579 998461 502 434 112 92 500 1.7702639 1007835 511 434 117 93 600 2.1274149 1011643 554 486 136 93 700 2.4850969 1015079 527 457 141 105 800 2.8374901 1015193 513 450 139 96 900 3.188709 1015271 512 460 117 97 1000 3.5468431 1015732 526 459 115 94 Loop time of 3.54699 on 128 procs for 1000 steps with 1015732 particles Particle moves = 1003319638 (1B) Cells touched = 1175336269 (1.18B) Particle comms = 17152084 (17.2M) Boundary collides = 0 (0K) Boundary exits = 2808260 (2.81M) SurfColl checks = 70434937 (70.4M) SurfColl occurs = 120704 (0.121M) Surf reactions = 0 (0K) Collide attempts = 14729081 (14.7M) Collide occurs = 11232168 (11.2M) Reactions = 0 (0K) Particles stuck = 0 Particle-moves/CPUsec/proc: 2.20988e+06 Particle-moves/step: 1.00332e+06 Cell-touches/particle/step: 1.17145 Particle comm iterations/step: 1 Particle fraction communicated: 0.0170953 Particle fraction colliding with boundary: 0 Particle fraction exiting boundary: 0.00279897 Surface-checks/particle/step: 0.0702019 Surface-collisions/particle/step: 0.000120305 Surf-reactions/particle/step: 0 Collision-attempts/particle/step: 0.0146803 Collisions/particle/step: 0.011195 Reactions/particle/step: 0 Move time (%) = 2.50511 (70.6265) Coll time (%) = 0.328274 (9.25501) Sort time (%) = 0.33379 (9.41051) Comm time (%) = 0.260512 (7.34459) Modfy time (%) = 0.110112 (3.10438) Outpt time (%) = 0.00786634 (0.221775) Other time (%) = 0.0013205 (0.0372287) Particles: 7935.41 ave 8656 max 7571 min Histogram: 5 31 28 27 14 9 8 3 1 2 Cells: 800 ave 1222 max 582 min Histogram: 4 0 32 86 2 0 0 0 0 4 GhostCell: 736.031 ave 1219 max 389 min Histogram: 4 20 13 29 20 8 19 10 4 1 EmptyCell: 0 ave 0 max 0 min Histogram: 128 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) = 1.6875 1.6875 1.6875 grid (ave,min,max) = 1.66647 1.66647 1.66647 surf (ave,min,max) = 0.0779014 0.0778999 0.0779037 total (ave,min,max) = 3.43797 3.43631 3.44119 Step CPU Np Natt Ncoll Nscoll c_max 1000 0 1015732 526 459 0 94 1100 0.35891795 1015794 520 461 138 102 1200 0.71506286 1014711 504 445 135 84 1300 1.0711319 1015716 506 446 127 98 1400 1.428072 1015103 515 470 114 100 1500 1.782191 1016091 504 442 120 93 1600 2.134984 1015713 513 456 105 89 1700 2.487741 1014905 484 420 145 107 1800 2.8402328 1015096 495 438 135 91 1900 3.1934829 1015790 496 444 133 101 2000 3.5680079 1015244 514 457 134 88 Loop time of 3.56817 on 128 procs for 1000 steps with 1015244 particles Particle moves = 1018368727 (1.02B) Cells touched = 1192633088 (1.19B) Particle comms = 17427115 (17.4M) Boundary collides = 0 (0K) Boundary exits = 2847768 (2.85M) SurfColl checks = 74416905 (74.4M) SurfColl occurs = 123918 (0.124M) Surf reactions = 0 (0K) Collide attempts = 15578584 (15.6M) Collide occurs = 11893854 (11.9M) Reactions = 0 (0K) Particles stuck = 0 Particle-moves/CPUsec/proc: 2.22972e+06 Particle-moves/step: 1.01837e+06 Cell-touches/particle/step: 1.17112 Particle comm iterations/step: 1 Particle fraction communicated: 0.0171128 Particle fraction colliding with boundary: 0 Particle fraction exiting boundary: 0.0027964 Surface-checks/particle/step: 0.0730746 Surface-collisions/particle/step: 0.000121683 Surf-reactions/particle/step: 0 Collision-attempts/particle/step: 0.0152976 Collisions/particle/step: 0.0116793 Reactions/particle/step: 0 Move time (%) = 2.55583 (71.6287) Coll time (%) = 0.342365 (9.59498) Sort time (%) = 0.33898 (9.50011) Comm time (%) = 0.259699 (7.27821) Modfy time (%) = 0.0617859 (1.73159) Outpt time (%) = 0.00819612 (0.229701) Other time (%) = 0.00131108 (0.0367438) Particles: 7931.59 ave 7955 max 7909 min Histogram: 2 2 1 18 33 55 15 0 0 2 Cells: 800 ave 1184 max 575 min Histogram: 4 3 25 30 60 2 0 0 1 3 GhostCell: 704.977 ave 1138 max 300 min Histogram: 3 2 18 32 11 30 10 12 5 5 EmptyCell: 0 ave 0 max 0 min Histogram: 128 0 0 0 0 0 0 0 0 0