1SPARTA (7 Jul 2014) 2# advect particles via VSS collisional flow on a uniform grid 3# particles reflect off global box boundaries 4 5variable x index 10 6variable y index 10 7variable z index 10 8 9variable lx equal $x*1.0e-5 10variable lx equal 20*1.0e-5 11variable ly equal $y*1.0e-5 12variable ly equal 20*1.0e-5 13variable lz equal $z*1.0e-5 14variable lz equal 25*1.0e-5 15 16variable n equal 10*$x*$y*$z 17variable n equal 10*20*$y*$z 18variable n equal 10*20*20*$z 19variable n equal 10*20*20*25 20 21seed 12345 22dimension 3 23global gridcut 1.0e-5 24 25boundary rr rr rr 26 27create_box 0 ${lx} 0 ${ly} 0 ${lz} 28create_box 0 0.0002 0 ${ly} 0 ${lz} 29create_box 0 0.0002 0 0.0002 0 ${lz} 30create_box 0 0.0002 0 0.0002 0 0.00025 31Created orthogonal box = (0 0 0) to (0.0002 0.0002 0.00025) 32create_grid $x $y $z 33create_grid 20 $y $z 34create_grid 20 20 $z 35create_grid 20 20 25 36Created 10000 child grid cells 37 parent cells = 1 38 CPU time = 0.00143218 secs 39 create/ghost percent = 91.127 8.87298 40 41balance_grid rcb part 42Balance grid migrated 0 cells 43 CPU time = 0.0143762 secs 44 reassign/sort/migrate/ghost percent = 5.37663 0.243789 0.494212 93.8854 45 46species ar.species Ar 47mixture air Ar vstream 0.0 0.0 0.0 temp 273.15 48 49global nrho 7.07043E22 50global fnum 7.07043E6 51 52collide vss air ar.vss 53 54create_particles air n $n 55create_particles air n 100000 56Created 100000 particles 57 CPU time = 0.0186059 secs 58 59stats 10 60compute temp temp 61stats_style step cpu np nattempt ncoll c_temp 62 63# first equilibrate with large timestep to unsort particles 64# then benchmark with normal timestep 65 66timestep 7.00E-8 67run 30 68Memory usage per proc in Mbytes: 69 particles (ave,min,max) = 10.2997 10.2997 10.2997 70 grid (ave,min,max) = 1.83113 1.83113 1.83113 71 surf (ave,min,max) = 0 0 0 72 total (ave,min,max) = 12.1308 12.1308 12.1308 73Step CPU Np Natt Ncoll temp 74 0 0 100000 0 0 274.57169 75 10 0.36405206 100000 96414 70599 274.57169 76 20 0.75538611 100000 98884 70648 274.57169 77 30 1.1479352 100000 100568 70919 274.57169 78Loop time of 1.14796 on 1 procs for 30 steps with 100000 particles 79 80Particle moves = 3000000 (3M) 81Cells touched = 14457156 (14.5M) 82Particle comms = 0 (0K) 83Boundary collides = 561056 (0.561M) 84Boundary exits = 0 (0K) 85SurfColl checks = 0 (0K) 86SurfColl occurs = 0 (0K) 87Collide attempts = 2907303 (2.91M) 88Collide occurs = 2109908 (2.11M) 89Particles stuck = 0 90 91Particle-moves/CPUsec/proc: 2.61332e+06 92Particle-moves/step: 100000 93Cell-touches/particle/step: 4.81905 94Particle comm iterations/step: 1 95Particle fraction communicated: 0 96Particle fraction colliding with boundary: 0.187019 97Particle fraction exiting boundary: 0 98Surface-checks/particle/step: 0 99Surface-collisions/particle/step: 0 100Collision-attempts/particle/step: 0.969101 101Collisions/particle/step: 0.703303 102 103Move time (%) = 0.601222 (52.3729) 104Coll time (%) = 0.510326 (44.4548) 105Sort time (%) = 0.0343685 (2.99386) 106Comm time (%) = 6.7234e-05 (0.0058568) 107Outpt time (%) = 0.00194478 (0.169411) 108Other time (%) = 3.6478e-05 (0.00317763) 109 110Particles: 100000 ave 100000 max 100000 min 111Histogram: 1 0 0 0 0 0 0 0 0 0 112Cells: 10000 ave 10000 max 10000 min 113Histogram: 1 0 0 0 0 0 0 0 0 0 114GhostCell: 0 ave 0 max 0 min 115Histogram: 1 0 0 0 0 0 0 0 0 0 116EmptyCell: 0 ave 0 max 0 min 117Histogram: 1 0 0 0 0 0 0 0 0 0 118timestep 7.00E-9 119run 100 120Memory usage per proc in Mbytes: 121 particles (ave,min,max) = 10.2997 10.2997 10.2997 122 grid (ave,min,max) = 1.83113 1.83113 1.83113 123 surf (ave,min,max) = 0 0 0 124 total (ave,min,max) = 12.1308 12.1308 12.1308 125Step CPU Np Natt Ncoll temp 126 30 0 100000 100568 70919 274.57169 127 40 0.075831175 100000 8986 7092 274.57169 128 50 0.15269303 100000 9153 7098 274.57169 129 60 0.22981811 100000 9335 7013 274.57169 130 70 0.30708814 100000 9401 7084 274.57169 131 80 0.38458395 100000 9507 7049 274.57169 132 90 0.46223307 100000 9473 7003 274.57169 133 100 0.53989601 100000 9574 7051 274.57169 134 110 0.61755514 100000 9601 7057 274.57169 135 120 0.69535398 100000 9680 7035 274.57169 136 130 0.77528715 100000 9791 7073 274.57169 137Loop time of 0.775309 on 1 procs for 100 steps with 100000 particles 138 139Particle moves = 10000000 (10M) 140Cells touched = 13818853 (13.8M) 141Particle comms = 0 (0K) 142Boundary collides = 186694 (0.187M) 143Boundary exits = 0 (0K) 144SurfColl checks = 0 (0K) 145SurfColl occurs = 0 (0K) 146Collide attempts = 937935 (0.938M) 147Collide occurs = 702813 (0.703M) 148Particles stuck = 0 149 150Particle-moves/CPUsec/proc: 1.28981e+07 151Particle-moves/step: 100000 152Cell-touches/particle/step: 1.38189 153Particle comm iterations/step: 1 154Particle fraction communicated: 0 155Particle fraction colliding with boundary: 0.0186694 156Particle fraction exiting boundary: 0 157Surface-checks/particle/step: 0 158Surface-collisions/particle/step: 0 159Collision-attempts/particle/step: 0.0937935 160Collisions/particle/step: 0.0702813 161 162Move time (%) = 0.342805 (44.2153) 163Coll time (%) = 0.313939 (40.4921) 164Sort time (%) = 0.111967 (14.4416) 165Comm time (%) = 0.000215769 (0.02783) 166Outpt time (%) = 0.0063014 (0.81276) 167Other time (%) = 8.13007e-05 (0.0104862) 168 169Particles: 100000 ave 100000 max 100000 min 170Histogram: 1 0 0 0 0 0 0 0 0 0 171Cells: 10000 ave 10000 max 10000 min 172Histogram: 1 0 0 0 0 0 0 0 0 0 173GhostCell: 0 ave 0 max 0 min 174Histogram: 1 0 0 0 0 0 0 0 0 0 175EmptyCell: 0 ave 0 max 0 min 176Histogram: 1 0 0 0 0 0 0 0 0 0 177 178