SPARTA (23 Dec 2017) # advect particles via VSS collisional flow on a uniform grid # particles reflect off global box boundaries variable x index 10 variable y index 10 variable z index 10 variable lx equal $x*1.0e-5 variable lx equal 128*1.0e-5 variable ly equal $y*1.0e-5 variable ly equal 80*1.0e-5 variable lz equal $z*1.0e-5 variable lz equal 80*1.0e-5 variable n equal 10*$x*$y*$z variable n equal 10*128*$y*$z variable n equal 10*128*80*$z variable n equal 10*128*80*80 seed 12345 dimension 3 global gridcut 1.0e-5 boundary rr rr rr create_box 0 ${lx} 0 ${ly} 0 ${lz} create_box 0 0.00128 0 ${ly} 0 ${lz} create_box 0 0.00128 0 0.0008 0 ${lz} create_box 0 0.00128 0 0.0008 0 0.0008 Created orthogonal box = (0 0 0) to (0.00128 0.0008 0.0008) create_grid $x $y $z create_grid 128 $y $z create_grid 128 80 $z create_grid 128 80 80 WARNING: Could not acquire nearby ghost cells b/c grid partition is not clumped (../grid.cpp:381) Created 819200 child grid cells parent cells = 1 CPU time = 0.00547194 secs create/ghost percent = 68.7508 31.2492 balance_grid rcb part Balance grid migrated 806400 cells CPU time = 0.078418 secs reassign/sort/migrate/ghost percent = 26.328 0.487065 21.893 51.292 species ar.species Ar mixture air Ar vstream 0.0 0.0 0.0 temp 273.15 global nrho 7.07043E22 global fnum 7.07043E6 collide vss air ar.vss create_particles air n $n create_particles air n 8192000 Created 8192000 particles CPU time = 0.02302 secs stats 100 compute temp temp stats_style step cpu np nattempt ncoll c_temp # first equilibrate with large timestep to unsort particles # then benchmark with normal timestep timestep 7.00E-8 run 30 Memory usage per proc in Mbytes: particles (ave,min,max) = 6.75 6.75 6.75 grid (ave,min,max) = 2.43673 1.51388 2.45138 surf (ave,min,max) = 0 0 0 total (ave,min,max) = 9.18673 8.26388 9.20138 Step CPU Np Natt Ncoll c_temp 0 0 8192000 0 0 273.10066 30 1.677686 8192000 8219584 5785021 273.10066 Loop time of 1.67771 on 128 procs for 30 steps with 8192000 particles Particle moves = 245760000 (246M) Cells touched = 1219466825 (1.22B) Particle comms = 40614921 (40.6M) Boundary collides = 10740751 (10.7M) Boundary exits = 0 (0K) SurfColl checks = 0 (0K) SurfColl occurs = 0 (0K) Surf reactions = 0 (0K) Collide attempts = 237957411 (238M) Collide occurs = 172770718 (173M) Reactions = 0 (0K) Particles stuck = 0 Particle-moves/CPUsec/proc: 1.14442e+06 Particle-moves/step: 8.192e+06 Cell-touches/particle/step: 4.96202 Particle comm iterations/step: 3.03333 Particle fraction communicated: 0.165263 Particle fraction colliding with boundary: 0.0437042 Particle fraction exiting boundary: 0 Surface-checks/particle/step: 0 Surface-collisions/particle/step: 0 Surf-reactions/particle/step: 0 Collision-attempts/particle/step: 0.968251 Collisions/particle/step: 0.703006 Reactions/particle/step: 0 Move time (%) = 0.657252 (39.1756) Coll time (%) = 0.612886 (36.5311) Sort time (%) = 0.126322 (7.52943) Comm time (%) = 0.277378 (16.5331) Modfy time (%) = 0 (0) Outpt time (%) = 0.00383741 (0.228729) Other time (%) = 3.33805e-05 (0.00198965) Particles: 64000 ave 64612 max 63342 min Histogram: 2 6 10 11 25 34 18 14 5 3 Cells: 6400 ave 6400 max 6400 min Histogram: 128 0 0 0 0 0 0 0 0 0 GhostCell: 3692.47 ave 5120 max 1538 min Histogram: 2 4 8 20 14 24 8 24 0 24 EmptyCell: 1636.48 ave 2808 max 441 min Histogram: 3 17 0 37 1 25 17 13 10 5 timestep 7.00E-9 run $t run 100 Memory usage per proc in Mbytes: particles (ave,min,max) = 6.75 6.75 6.75 grid (ave,min,max) = 2.43673 1.51388 2.45138 surf (ave,min,max) = 0 0 0 total (ave,min,max) = 9.18673 8.26388 9.20138 Step CPU Np Natt Ncoll c_temp 30 0 8192000 8219584 5785021 273.10066 100 1.2982361 8192000 785641 578437 273.10066 130 1.858108 8192000 796410 578743 273.10066 Loop time of 1.85813 on 128 procs for 100 steps with 8192000 particles Particle moves = 819200000 (819M) Cells touched = 1142931857 (1.14B) Particle comms = 14067862 (14.1M) Boundary collides = 3580661 (3.58M) Boundary exits = 0 (0K) SurfColl checks = 0 (0K) SurfColl occurs = 0 (0K) Surf reactions = 0 (0K) Collide attempts = 76743469 (76.7M) Collide occurs = 57489265 (57.5M) Reactions = 0 (0K) Particles stuck = 0 Particle-moves/CPUsec/proc: 3.44433e+06 Particle-moves/step: 8.192e+06 Cell-touches/particle/step: 1.39518 Particle comm iterations/step: 1 Particle fraction communicated: 0.0171727 Particle fraction colliding with boundary: 0.00437092 Particle fraction exiting boundary: 0 Surface-checks/particle/step: 0 Surface-collisions/particle/step: 0 Surf-reactions/particle/step: 0 Collision-attempts/particle/step: 0.093681 Collisions/particle/step: 0.0701773 Reactions/particle/step: 0 Move time (%) = 0.956096 (51.4548) Coll time (%) = 0.400042 (21.5293) Sort time (%) = 0.419065 (22.5531) Comm time (%) = 0.0755734 (4.06718) Modfy time (%) = 0 (0) Outpt time (%) = 0.00724021 (0.389651) Other time (%) = 0.000111945 (0.00602461) Particles: 64000 ave 64768 max 63430 min Histogram: 9 11 12 23 25 26 11 6 3 2 Cells: 6400 ave 6400 max 6400 min Histogram: 128 0 0 0 0 0 0 0 0 0 GhostCell: 3692.47 ave 5120 max 1538 min Histogram: 2 4 8 20 14 24 8 24 0 24 EmptyCell: 1636.48 ave 2808 max 441 min Histogram: 3 17 0 37 1 25 17 13 10 5