SPARTA (23 Dec 2017) # advect particles via free molecular 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 64*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*64*$y*$z variable n equal 10*64*80*$z variable n equal 10*64*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.00064 0 ${ly} 0 ${lz} create_box 0 0.00064 0 0.0008 0 ${lz} create_box 0 0.00064 0 0.0008 0 0.0008 Created orthogonal box = (0 0 0) to (0.00064 0.0008 0.0008) create_grid $x $y $z create_grid 64 $y $z create_grid 64 80 $z create_grid 64 80 80 WARNING: Could not acquire nearby ghost cells b/c grid partition is not clumped (../grid.cpp:381) Created 409600 child grid cells parent cells = 1 CPU time = 0.0207412 secs create/ghost percent = 40.9483 59.0517 balance_grid rcb part Balance grid migrated 408800 cells CPU time = 0.287677 secs reassign/sort/migrate/ghost percent = 76.5497 0.0566051 10.0237 13.3699 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 create_particles air n $n create_particles air n 4096000 Created 4096000 particles CPU time = 0.01543 secs stats 100 compute temp temp stats_style step cpu np nattempt ncoll c_temp # 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) = 1.6875 1.6875 1.6875 grid (ave,min,max) = 1.51388 1.51388 1.51388 surf (ave,min,max) = 0 0 0 total (ave,min,max) = 3.20138 3.20138 3.20138 Step CPU Np Natt Ncoll c_temp 0 0 4096000 0 0 273.25033 30 0.36512303 4096000 0 0 273.25033 Loop time of 0.365287 on 512 procs for 30 steps with 4096000 particles Particle moves = 122880000 (123M) Cells touched = 609595083 (610M) Particle comms = 42183657 (42.2M) Boundary collides = 6651988 (6.65M) Boundary exits = 0 (0K) SurfColl checks = 0 (0K) SurfColl occurs = 0 (0K) Surf reactions = 0 (0K) Collide attempts = 0 (0K) Collide occurs = 0 (0K) Reactions = 0 (0K) Particles stuck = 0 Particle-moves/CPUsec/proc: 657017 Particle-moves/step: 4.096e+06 Cell-touches/particle/step: 4.9609 Particle comm iterations/step: 3.43333 Particle fraction communicated: 0.343291 Particle fraction colliding with boundary: 0.054134 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 Collisions/particle/step: 0 Reactions/particle/step: 0 Move time (%) = 0.24587 (67.3086) Coll time (%) = 0 (0) Sort time (%) = 0 (0) Comm time (%) = 0.118671 (32.4871) Modfy time (%) = 0 (0) Outpt time (%) = 0.000615905 (0.168608) Other time (%) = 0.000130375 (0.0356911) Particles: 8000 ave 8248 max 7740 min Histogram: 4 19 38 69 108 114 76 60 18 6 Cells: 800 ave 800 max 800 min Histogram: 512 0 0 0 0 0 0 0 0 0 GhostCell: 1234.88 ave 1552 max 640 min Histogram: 8 12 42 24 61 50 115 25 50 125 EmptyCell: 594.881 ave 912 max 144 min Histogram: 5 29 11 83 48 111 57 64 56 48 timestep 7.00E-9 run $t run 100 Memory usage per proc in Mbytes: particles (ave,min,max) = 1.6875 1.6875 1.6875 grid (ave,min,max) = 1.51388 1.51388 1.51388 surf (ave,min,max) = 0 0 0 total (ave,min,max) = 3.20138 3.20138 3.20138 Step CPU Np Natt Ncoll c_temp 30 0 4096000 0 0 273.25033 100 0.16746902 4096000 0 0 273.25033 130 0.23947811 4096000 0 0 273.25033 Loop time of 0.23963 on 512 procs for 100 steps with 4096000 particles Particle moves = 409600000 (410M) Cells touched = 571096905 (571M) Particle comms = 15323685 (15.3M) Boundary collides = 2215578 (2.22M) Boundary exits = 0 (0K) SurfColl checks = 0 (0K) SurfColl occurs = 0 (0K) Surf reactions = 0 (0K) Collide attempts = 0 (0K) Collide occurs = 0 (0K) Reactions = 0 (0K) Particles stuck = 0 Particle-moves/CPUsec/proc: 3.33848e+06 Particle-moves/step: 4.096e+06 Cell-touches/particle/step: 1.39428 Particle comm iterations/step: 1 Particle fraction communicated: 0.0374113 Particle fraction colliding with boundary: 0.00540913 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 Collisions/particle/step: 0 Reactions/particle/step: 0 Move time (%) = 0.166639 (69.5403) Coll time (%) = 0 (0) Sort time (%) = 0 (0) Comm time (%) = 0.0714085 (29.7995) Modfy time (%) = 0 (0) Outpt time (%) = 0.00132524 (0.553036) Other time (%) = 0.00025681 (0.107169) Particles: 8000 ave 8303 max 7735 min Histogram: 6 26 50 93 117 114 61 39 2 4 Cells: 800 ave 800 max 800 min Histogram: 512 0 0 0 0 0 0 0 0 0 GhostCell: 1234.88 ave 1552 max 640 min Histogram: 8 12 42 24 61 50 115 25 50 125 EmptyCell: 594.881 ave 912 max 144 min Histogram: 5 29 11 83 48 111 57 64 56 48