jianmu
/
lammps
forked from lijiext/lammps
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@11550 f3b2605a-c512-4ea7-a41b-209d697bcdaa

This commit is contained in:
sjplimp 2014-02-12 22:52:36 +00:00
parent 75e0153477
commit 3791069665
13 changed files with 3349 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

290
examples/ASPHERE/box/log.box.1Feb14 Normal file
View File

@ -0,0 +1,290 @@
LAMMPS (1 Feb 2014)
# SRD diffusion demo - rigid box particles
units lj
atom_style sphere
atom_modify map array first big
dimension 2
# read in clusters of rigid bodies
fix molprop all property/atom mol
read_data data.box fix molprop NULL Molecules
orthogonal box = (-19.1447 -19.1447 -0.5) to (19.1447 19.1447 0.5)
4 by 2 by 1 MPI processor grid
reading atoms ...
420 atoms
set type 1 mass 1.0
420 settings made for mass
group big type 1
420 atoms in group big
velocity big create 1.44 87287 loop geom
# equilibrate big particles
pair_style soft 1.12
pair_coeff 1 1 0.0
pair_coeff 2 2 0.0 0.0
pair_coeff 1 2 0.0 0.0
variable prefactor equal ramp(0,60)
fix soft all adapt 1 pair soft a * * v_prefactor
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
fix 1 big rigid molecule
30 rigid bodies with 420 atoms
fix 2 all enforce2d
#dump 1 all atom 10 dump.box.equil
thermo 100
run 1000
Memory usage per processor = 0.443012 Mbytes
Step Temp E_pair E_mol TotEng Press
0 0 0 0 0 -8.770036
100 0 3.4203218 0 3.4203218 0.49351517
200 0 5.3315601 0 5.3315601 0.26965614
300 0 7.2094878 0 7.2094878 0.43326479
400 0 9.2008598 0 9.2008598 0.27649318
500 0 11.562353 0 11.562353 0.60730421
600 0 13.731564 0 13.731564 0.27272312
700 0 15.850142 0 15.850142 0.0032862218
800 0 18.395444 0 18.395444 1.0272658
900 0 20.3322 0 20.3322 -0.16729762
1000 0 22.600716 0 22.600716 -0.13023705
Loop time of 0.0534293 on 8 procs for 1000 steps with 420 atoms
Pair time (%) = 0.00388658 (7.27425)
Neigh time (%) = 0.0012759 (2.38801)
Comm time (%) = 0.0151907 (28.4315)
Outpt time (%) = 0.000330567 (0.618701)
Other time (%) = 0.0327455 (61.2876)
Nlocal: 52.5 ave 67 max 35 min
Histogram: 1 1 1 0 1 0 0 1 1 2
Nghost: 30.625 ave 41 max 17 min
Histogram: 1 0 0 1 1 1 2 0 1 1
Neighs: 70.25 ave 90 max 44 min
Histogram: 2 0 0 1 1 0 1 0 0 3
Total # of neighbors = 562
Ave neighs/atom = 1.3381
Neighbor list builds = 158
Dangerous builds = 0
#undump 1
unfix soft
unfix 1
unfix 2
# add small particles as hi density lattice
region plane block INF INF INF INF -0.001 0.001 units box
lattice sq 85.0
Lattice spacing in x,y,z = 0.108465 0.108465 0.108465
create_atoms 2 region plane
Created 124609 atoms
set type 2 mass 0.1
124609 settings made for mass
group small type 2
124609 atoms in group small
velocity small create 1.0 593849 loop geom
# delete overlaps
# must set 1-2 cutoff to non-zero value
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 1.0
delete_atoms overlap 1.0 small big
Deleted 51030 atoms, new total = 73999
# SRD run
reset_timestep 0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
communicate multi group big vel yes
neigh_modify include big
# no pairwise interactions with small particles
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 0.0
# use fix SRD to push small particles out from inside big ones
# if comment out, big particles won't see SRD particles
timestep 0.001
fix 1 big rigid molecule
30 rigid bodies with 420 atoms
fix 2 small srd 20 big 1.0 0.25 49894 shift yes 54979 search 0.2 collision slip inside ignore overlap yes
fix 3 all enforce2d
# diagnostics
compute tbig big temp/sphere
variable pebig equal pe*atoms/count(big)
variable ebig equal etotal*atoms/count(big)
thermo_style custom step temp f_2[8] etotal v_pebig v_ebig press f_2[1] f_2[2] f_2[3] f_2[4] f_2[5] f_2[6] f_2[7] f_2[8] f_2[9] f_2[10] f_2[11] f_2[12]
thermo_modify temp tbig
WARNING: Temperature for thermo pressure is not for group all (../thermo.cpp:439)
thermo 1000
#dump 1 all atom 1000 dump.box
#dump 1 all image 1000 image.*.jpg type type zoom 1.6
#dump_modify 1 pad 6 adiam 1 1 adiam 2 0.2
run 100000
SRD info:
SRD/big particles = 73579 420
big particle diameter max/min = 1 1
SRD temperature & lamda = 1 0.0632456
SRD max distance & max velocity = 0.252982 12.6491
SRD grid counts: 153 153 1
SRD grid size: request, actual (xyz) = 0.25, 0.250257 0.250257 1
SRD per actual grid cell = 4.05573
SRD viscosity = 0.229297
big/SRD mass density ratio = 0.196614
WARNING: Fix srd grid size > 1/4 of big particle diameter (../fix_srd.cpp:2875)
# of rescaled SRD velocities = 0
ave/max small velocity = 4.19668 7.72762
ave/max big velocity = 2.12757 5.83591
Memory usage per processor = 14.2707 Mbytes
Step Temp 2[8] TotEng pebig ebig Press 2[1] 2[2] 2[3] 2[4] 2[5] 2[6] 2[7] 2[8] 2[9] 2[10] 2[11] 2[12]
0 27.381112 0 0.25599398 42.234595 45.103092 0.74047865 0 0 0 0 0 0 0 0 0 0 0 0
1000 3.4640099 16291 0.24158286 42.201127 42.564024 10.312083 10282 38 38 0 957 90440 23409 16291 0.91054529 0 69 0
2000 1.4820124 16228 0.24022248 42.169082 42.324341 9.5969069 10066 30 30 0 1726 90440 23409 16228 0.91140743 0 69 0
3000 1.3797075 16198 0.24017061 42.170661 42.315202 11.093574 10151 43 43 0 2410 90440 23409 16198 0.91909014 0 69 0
4000 1.3403741 16287 0.24014493 42.170257 42.310677 10.101336 10172 38 38 0 3094 90440 23409 16287 0.92752333 0 69 0
5000 1.0444432 16295 0.2399965 42.175108 42.284525 9.5498999 10375 26 26 0 3799 90440 23409 16295 0.93506273 0 69 0
6000 0.78295352 16299 0.2398933 42.184319 42.266343 9.535916 10090 37 37 0 4563 90440 23409 16299 0.93041308 0 69 0
7000 1.1675613 16333 0.24006565 42.174393 42.296709 11.594065 10426 37 37 0 5266 90440 23409 16333 0.94031182 0 69 0
8000 1.1343357 16302 0.24001364 42.168711 42.287546 9.8091173 9979 39 39 0 5993 90440 23409 16302 0.94116994 0 69 0
9000 1.4165164 16306 0.24010969 42.156072 42.304469 10.736953 10191 38 39 1 6683 90440 23409 16306 0.93957787 0 69 0
10000 0.63657628 16306 0.23958347 42.145066 42.211755 6.3923231 10115 37 37 0 7409 90440 23409 16306 0.9420945 0 69 0
11000 0.9421837 16286 0.23976832 42.145619 42.244324 9.9611184 10079 45 45 0 8122 90440 23409 16286 0.94240902 0 69 0
12000 0.9735159 16295 0.23980202 42.148274 42.250262 11.862106 10345 36 36 0 8854 90440 23409 16295 0.94994179 0 69 0
13000 0.93356977 16392 0.23970059 42.134589 42.232391 11.29915 10321 41 41 0 9607 90440 23409 16392 0.94747955 0 69 0
14000 1.0596151 16306 0.2397451 42.129226 42.240233 12.305503 10075 37 37 0 10323 90440 23409 16306 0.94857373 0 69 0
15000 1.4275311 16285 0.23986502 42.11181 42.261361 8.6785032 9963 43 43 0 11054 90440 23409 16285 0.95130035 0 69 0
16000 1.1798179 16367 0.23973717 42.115236 42.238836 9.7201163 10138 38 38 0 11785 90440 23409 16367 0.94838949 0 69 0
17000 1.2809316 16368 0.23978706 42.113433 42.247626 7.4568393 10239 31 31 0 12545 90440 23409 16368 0.94818482 0 69 0
18000 1.0552359 16355 0.23962676 42.108834 42.219382 11.024507 10270 39 39 0 13350 90440 23409 16355 0.94825339 0 69 0
19000 1.1373204 16343 0.23964456 42.103371 42.222519 9.1505749 10245 37 37 0 14198 90440 23409 16343 0.94776971 0 69 0
20000 0.89607862 16295 0.23962161 42.1246 42.218475 8.377591 10173 48 48 0 14985 90440 23409 16295 0.95836326 0 69 0
21000 0.93823225 16323 0.23958542 42.113808 42.212099 9.9128627 10398 39 39 0 15722 90440 23409 16323 0.94732512 0 69 0
22000 1.0132051 16270 0.23958637 42.106121 42.212267 12.521464 10414 37 37 1 16538 90440 23409 16270 0.94611956 0 69 0
23000 1.2399208 16420 0.2396696 42.097034 42.226931 9.1607507 10169 32 32 0 17263 90440 23409 16420 0.95141771 0 69 0
24000 1.2359103 16395 0.23957354 42.080529 42.210006 10.051289 10225 34 34 1 18057 90440 23409 16395 0.95487046 0 69 0
25000 0.96280337 16322 0.23944374 42.086271 42.187136 11.818136 10183 40 40 0 18818 90440 23409 16322 0.95691179 0 69 0
26000 0.83139747 16344 0.23940306 42.09287 42.179968 10.844701 10090 38 38 0 19555 90440 23409 16344 0.95332651 0 69 0
27000 1.107274 16270 0.23949948 42.080957 42.196957 14.435937 10028 46 46 0 20355 90440 23409 16270 0.95321506 0 69 0
28000 1.0892583 16399 0.23938721 42.063064 42.177177 9.3768225 10011 38 38 0 21195 90440 23409 16399 0.95935372 0 69 0
29000 1.1099534 16320 0.23940855 42.064655 42.180936 8.8044985 9919 34 34 0 22020 90440 23409 16320 0.95465302 0 69 0
30000 1.0119272 16302 0.23934223 42.06324 42.169252 12.13721 10160 47 47 0 22750 90440 23409 16302 0.94922806 0 69 0
31000 0.9567755 16335 0.23916982 42.038642 42.138875 8.3903033 10007 30 30 0 23497 90440 23409 16335 0.95223124 0 69 0
32000 0.86048362 16291 0.23910468 42.037252 42.127398 11.646647 10347 41 41 0 24266 90440 23409 16291 0.95484469 0 69 0
33000 1.1129645 16284 0.23928206 42.042053 42.15865 9.0389494 9909 33 33 0 25073 90440 23409 16284 0.96246044 0 69 0
34000 1.0741874 16303 0.23920774 42.033022 42.145556 12.710332 10125 35 35 0 25845 90440 23409 16303 0.95343962 0 69 0
35000 0.87606515 16293 0.23916564 42.046361 42.138139 13.784203 10044 42 42 0 26569 90440 23409 16293 0.95827172 0 69 0
36000 0.8268686 16315 0.23912916 42.045087 42.131712 13.474499 10240 43 44 0 27324 90440 23409 16315 0.9539089 0 69 0
37000 1.1732368 16336 0.23930076 42.039035 42.161945 8.7392785 9996 34 34 0 28094 90440 23409 16336 0.95672234 0 69 0
38000 1.089263 16323 0.2391967 42.029498 42.143611 9.464664 10104 39 39 0 28856 90440 23409 16323 0.95362271 0 69 0
39000 1.06444 16335 0.23919427 42.03167 42.143183 10.036973 9961 36 36 0 29607 90440 23409 16335 0.95464483 0 69 0
40000 1.1390869 16376 0.23927054 42.037288 42.156621 9.8133871 10206 38 38 0 30409 90440 23409 16376 0.95191201 0 69 0
41000 0.76020035 16347 0.23906101 42.040064 42.119704 10.826738 10406 41 41 0 31195 90440 23409 16347 0.94685002 0 69 0
42000 1.3087439 16311 0.23945151 42.051398 42.188505 11.375199 10219 45 45 0 31938 90440 23409 16311 0.95326292 0 69 0
43000 0.98306175 16313 0.23923796 42.047893 42.15088 10.536429 10269 47 47 0 32660 90440 23409 16313 0.94791465 0 69 0
44000 1.1067947 16299 0.23929754 42.045427 42.161377 10.366386 10166 33 33 0 33471 90440 23409 16299 0.9520033 0 69 0
45000 1.1754605 16269 0.2394174 42.059352 42.182496 7.898603 9985 37 37 0 34245 90440 23409 16269 0.95500544 0 69 0
46000 0.93860976 16272 0.23926541 42.057387 42.155717 11.566786 10108 35 35 0 35050 90440 23409 16272 0.95183186 0 69 0
47000 1.3576761 16315 0.23942 42.040721 42.182954 10.037574 9902 33 33 0 35807 90440 23409 16315 0.95525561 0 69 0
48000 1.3465813 16308 0.23938501 42.035719 42.176789 8.9584586 9830 33 33 0 36625 90440 23409 16308 0.95044669 0 69 0
49000 0.93823974 16340 0.23917669 42.041793 42.140085 10.474409 9930 35 35 0 37393 90440 23409 16340 0.94869206 0 69 0
50000 1.0643183 16332 0.23923523 42.0389 42.1504 8.9372951 9886 43 43 0 38182 90440 23409 16332 0.94865149 0 69 0
51000 0.87518396 16310 0.23914944 42.043598 42.135284 11.284535 9695 40 40 0 38995 90440 23409 16310 0.9560088 0 69 0
52000 1.1910982 16336 0.23929054 42.035364 42.160145 12.437325 9961 44 44 0 39787 90440 23409 16336 0.95735399 0 69 0
53000 1.0023949 16384 0.23923896 42.046043 42.151056 10.386358 9859 36 36 0 40569 90440 23409 16384 0.95693245 0 69 0
54000 0.98715708 16261 0.23921241 42.042962 42.146379 9.9941499 10027 35 35 0 41343 90440 23409 16261 0.951038 0 69 0
55000 1.0004435 16299 0.23914811 42.030241 42.135049 12.345206 9830 47 48 0 42056 90440 23409 16299 0.95562154 0 69 0
56000 1.0479458 16296 0.23914874 42.025377 42.135162 11.587952 9801 39 39 0 42915 90440 23409 16296 0.94868699 0 69 0
57000 0.94369119 16319 0.23902076 42.013749 42.112612 10.195934 9764 46 46 0 43652 90440 23409 16319 0.95438171 0 69 0
58000 1.1549569 16347 0.23924582 42.031269 42.152265 9.511812 9910 41 41 0 44366 90440 23409 16347 0.9484367 0 69 0
59000 0.81453047 16402 0.23895926 42.016445 42.101776 10.359564 9967 46 46 0 45115 90440 23409 16402 0.95256147 0 69 0
60000 0.9253868 16360 0.23897085 42.006873 42.103818 12.308731 9839 44 44 0 45872 90440 23409 16360 0.95753176 0 69 0
61000 0.95173279 16298 0.23908799 42.024752 42.124457 10.985428 9848 45 45 0 46621 90440 23409 16298 0.95937218 0 69 0
62000 1.0168886 16273 0.23910318 42.020603 42.127134 13.493586 9842 41 41 0 47356 90440 23409 16273 0.95720763 0 69 0
63000 0.59853847 16306 0.23888991 42.026854 42.089558 11.654045 10030 47 47 0 48152 90440 23409 16306 0.95164201 0 69 0
64000 0.88165231 16295 0.23904152 42.023907 42.11627 8.6017877 9808 31 31 0 48944 90440 23409 16295 0.95791125 0 69 0
65000 0.83089386 16306 0.23895929 42.014736 42.101782 11.57787 10042 51 51 0 49729 90440 23409 16306 0.95848466 0 69 0
66000 1.1172469 16324 0.23919188 42.025716 42.142761 11.423941 9802 55 55 0 50498 90440 23409 16324 0.94666759 0 69 0
67000 0.97470955 16356 0.23903681 42.013327 42.11544 11.345125 10043 40 40 0 51251 90440 23409 16356 0.95264876 0 69 0
68000 0.90592443 16274 0.23906764 42.025966 42.120872 11.050097 9782 44 44 0 52009 90440 23409 16274 0.95334971 0 69 0
69000 1.0820724 16352 0.23920813 42.032265 42.145625 11.892518 9803 43 43 0 52800 90440 23409 16352 0.95290952 0 69 0
70000 1.0413192 16301 0.23929529 42.051891 42.160982 13.193444 9847 48 48 0 53595 90440 23409 16301 0.95613844 0 69 0
71000 0.95433723 16430 0.23919345 42.04306 42.143038 12.116035 9898 46 46 0 54350 90440 23409 16430 0.95532726 0 69 0
72000 0.80923962 16318 0.23909626 42.041137 42.125915 10.240888 9800 37 37 0 55133 90440 23409 16318 0.95196829 0 69 0
73000 0.93773964 16330 0.23920083 42.046099 42.144338 13.100955 9780 37 37 0 55931 90440 23409 16330 0.95075325 0 69 0
74000 1.0547823 16296 0.23913003 42.021364 42.131865 12.215247 10308 45 45 0 56705 90440 23409 16296 0.95515593 0 69 0
75000 0.74939607 16287 0.23910386 42.048745 42.127253 10.270685 10092 43 43 0 57457 90440 23409 16287 0.95365363 0 69 0
76000 0.95221448 16326 0.23921923 42.047825 42.147581 10.33549 10191 47 47 0 58254 90440 23409 16326 0.95389896 0 69 0
77000 0.94854568 16249 0.23925669 42.054809 42.15418 10.317964 10055 35 35 1 59026 90440 23409 16249 0.95716276 0 69 0
78000 1.2663624 16282 0.23949043 42.062695 42.195362 10.496265 9893 41 41 0 59780 90440 23409 16282 0.94847998 0 69 0
79000 0.98933199 16264 0.2392933 42.056987 42.160631 12.876564 9946 41 41 0 60518 90440 23409 16264 0.95354593 0 69 0
80000 0.8212992 16341 0.23927469 42.07131 42.157351 8.5608431 9838 40 40 0 61296 90440 23409 16341 0.95234534 0 69 0
81000 0.98231967 16339 0.23928757 42.056712 42.159622 11.491156 10104 37 38 0 62061 90440 23409 16339 0.94886686 0 69 0
82000 0.99158973 16299 0.23925843 42.050607 42.154488 9.4538952 10047 37 37 0 62852 90440 23409 16299 0.95656275 0 69 0
83000 1.2341517 16344 0.239348 42.040977 42.170269 10.524371 9784 37 37 0 63608 90440 23409 16344 0.94780456 0 69 0
84000 1.1753084 16288 0.23928447 42.035947 42.159074 9.1748883 9831 37 37 0 64314 90440 23409 16288 0.94217045 0 69 0
85000 0.84912682 16310 0.23915288 42.046935 42.135891 11.43337 9890 42 42 0 64992 90440 23409 16310 0.9423407 0 69 0
86000 0.88939653 16275 0.23909757 42.032971 42.126146 12.34417 9961 46 46 0 65677 90440 23409 16275 0.93887016 0 69 0
87000 1.1303842 16248 0.23922093 42.029458 42.147879 11.934641 9945 47 47 0 66417 90440 23409 16248 0.93840688 0 69 0
88000 1.2294052 16326 0.23916023 42.00839 42.137185 9.7233305 9838 35 36 0 67149 90440 23409 16326 0.94945034 0 69 0
89000 1.0455638 16345 0.23918575 42.032147 42.141682 8.6325214 9871 35 35 0 67855 90440 23409 16345 0.94891424 0 69 0
90000 0.79047982 16288 0.23897011 42.020876 42.103688 10.271997 10053 45 45 0 68551 90440 23409 16288 0.94963234 0 69 0
91000 0.85451256 16320 0.23898107 42.016098 42.105619 12.567859 10053 42 42 0 69246 90440 23409 16320 0.9506585 0 69 0
92000 1.1082033 16363 0.23917252 42.023254 42.139351 11.730554 9858 46 46 0 70046 90440 23409 16363 0.95352701 0 69 0
93000 0.90368051 16312 0.23897927 42.010631 42.105303 10.70167 9793 37 37 0 70824 90440 23409 16312 0.95992717 0 69 0
94000 0.95749706 16350 0.23909592 42.025545 42.125855 10.174353 9970 42 42 0 71646 90440 23409 16350 0.95082103 0 69 0
95000 0.90723392 16298 0.23898239 42.010809 42.105853 11.561525 9938 33 33 0 72414 90440 23409 16298 0.94768554 0 69 0
96000 1.0249944 16288 0.23910408 42.019912 42.127292 10.755062 9945 38 38 0 73184 90440 23409 16288 0.94637131 0 69 0
97000 0.96263604 16339 0.23911151 42.027754 42.128601 9.9303923 9816 36 36 0 73853 90440 23409 16339 0.94851816 0 69 0
98000 0.93105892 16301 0.23897506 42.00702 42.10456 12.088691 9933 39 39 0 74604 90440 23409 16301 0.94785265 0 69 0
99000 1.2448977 16286 0.23920972 42.015487 42.145905 11.565661 10001 51 51 0 75328 90440 23409 16286 0.9483854 0 69 0
100000 0.95492648 16337 0.23919603 42.043453 42.143492 11.223512 9793 32 32 0 76045 90440 23409 16337 0.95143273 0 69 0
Loop time of 51.3335 on 8 procs for 100000 steps with 73999 atoms
Pair time (%) = 0.586967 (1.14344)
Neigh time (%) = 3.95145 (7.69761)
Comm time (%) = 2.97199 (5.78957)
Outpt time (%) = 0.0085291 (0.0166151)
Other time (%) = 43.8145 (85.3528)
Nlocal: 9249.88 ave 10642 max 7533 min
Histogram: 1 0 0 1 1 1 2 0 1 1
Nghost: 57.375 ave 89 max 34 min
Histogram: 1 2 0 1 1 0 2 0 0 1
Neighs: 335.875 ave 616 max 209 min
Histogram: 2 1 0 3 1 0 0 0 0 1
Total # of neighbors = 2687
Ave neighs/atom = 0.0363113
Neighbor list builds = 5000
Dangerous builds = 0
Please see the log.cite file for references relevant to this simulation

291
examples/ASPHERE/box/log.box.mp.1Feb14 Normal file
View File

@ -0,0 +1,291 @@
LAMMPS (1 Feb 2014)
# SRD viscosity demo - rigid box particles
units lj
atom_style sphere
atom_modify map array first big
dimension 2
# read in clusters of rigid bodies
fix molprop all property/atom mol
read_data data.box fix molprop NULL Molecules
orthogonal box = (-19.1447 -19.1447 -0.5) to (19.1447 19.1447 0.5)
4 by 2 by 1 MPI processor grid
reading atoms ...
420 atoms
set type 1 mass 1.0
420 settings made for mass
group big type 1
420 atoms in group big
velocity big create 1.44 87287 loop geom
# equilibrate big particles
pair_style soft 1.12
pair_coeff 1 1 0.0
pair_coeff 2 2 0.0 0.0
pair_coeff 1 2 0.0 0.0
variable prefactor equal ramp(0,60)
fix soft all adapt 1 pair soft a * * v_prefactor
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
fix 1 big rigid molecule
30 rigid bodies with 420 atoms
fix 2 all enforce2d
#dump 1 all atom 10 dump.box.equil
thermo 100
run 1000
Memory usage per processor = 0.443012 Mbytes
Step Temp E_pair E_mol TotEng Press
0 0 0 0 0 -8.770036
100 0 3.4203218 0 3.4203218 0.49351517
200 0 5.3315601 0 5.3315601 0.26965614
300 0 7.2094878 0 7.2094878 0.43326479
400 0 9.2008598 0 9.2008598 0.27649318
500 0 11.562353 0 11.562353 0.60730421
600 0 13.731564 0 13.731564 0.27272312
700 0 15.850142 0 15.850142 0.0032862218
800 0 18.395444 0 18.395444 1.0272658
900 0 20.3322 0 20.3322 -0.16729762
1000 0 22.600716 0 22.600716 -0.13023705
Loop time of 0.0540034 on 8 procs for 1000 steps with 420 atoms
Pair time (%) = 0.00388333 (7.1909)
Neigh time (%) = 0.0012731 (2.35743)
Comm time (%) = 0.0160792 (29.7743)
Outpt time (%) = 0.000274509 (0.508318)
Other time (%) = 0.0324934 (60.169)
Nlocal: 52.5 ave 67 max 35 min
Histogram: 1 1 1 0 1 0 0 1 1 2
Nghost: 30.625 ave 41 max 17 min
Histogram: 1 0 0 1 1 1 2 0 1 1
Neighs: 70.25 ave 90 max 44 min
Histogram: 2 0 0 1 1 0 1 0 0 3
Total # of neighbors = 562
Ave neighs/atom = 1.3381
Neighbor list builds = 158
Dangerous builds = 0
#undump 1
unfix soft
unfix 1
unfix 2
# add small particles as hi density lattice
region plane block INF INF INF INF -0.001 0.001 units box
lattice sq 85.0
Lattice spacing in x,y,z = 0.108465 0.108465 0.108465
create_atoms 2 region plane
Created 124609 atoms
set type 2 mass 0.1
124609 settings made for mass
group small type 2
124609 atoms in group small
velocity small create 1.0 593849 loop geom
# delete overlaps
# must set 1-2 cutoff to non-zero value
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 1.0
delete_atoms overlap 1.0 small big
Deleted 51030 atoms, new total = 73999
# SRD run
reset_timestep 0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
communicate multi group big vel yes
neigh_modify include big
# no pairwise interactions with small particles
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 0.0
# use fix SRD to push small particles out from inside big ones
# if comment out, big particles won't see SRD particles
timestep 0.001
fix 1 big rigid molecule
30 rigid bodies with 420 atoms
fix 2 small srd 20 big 1.0 0.25 49894 shift yes 54979 search 0.2 collision slip inside ignore overlap yes
fix 3 small viscosity 10 x y 50
fix 4 all enforce2d
# diagnostics
compute tbig big temp/sphere
variable pebig equal pe*atoms/count(big)
variable ebig equal etotal*atoms/count(big)
thermo_style custom step temp f_2[8] etotal v_pebig v_ebig press f_2[1] f_2[2] f_2[3] f_2[4] f_2[5] f_2[6] f_2[7] f_2[8] f_2[9] f_2[10] f_2[11] f_2[12]
thermo_modify temp tbig
WARNING: Temperature for thermo pressure is not for group all (../thermo.cpp:439)
thermo 1000
#dump 1 all atom 1000 dump.box.mp
#dump 1 all image 1000 image.*.jpg type type zoom 1.6
#dump_modify 1 pad 6 adiam 1 1 adiam 2 0.2
run 100000
SRD info:
SRD/big particles = 73579 420
big particle diameter max/min = 1 1
SRD temperature & lamda = 1 0.0632456
SRD max distance & max velocity = 0.252982 12.6491
SRD grid counts: 153 153 1
SRD grid size: request, actual (xyz) = 0.25, 0.250257 0.250257 1
SRD per actual grid cell = 4.05573
SRD viscosity = 0.229297
big/SRD mass density ratio = 0.196614
WARNING: Fix srd grid size > 1/4 of big particle diameter (../fix_srd.cpp:2875)
# of rescaled SRD velocities = 0
ave/max small velocity = 4.19668 7.72762
ave/max big velocity = 2.12757 5.83591
Memory usage per processor = 14.2707 Mbytes
Step Temp 2[8] TotEng pebig ebig Press 2[1] 2[2] 2[3] 2[4] 2[5] 2[6] 2[7] 2[8] 2[9] 2[10] 2[11] 2[12]
0 27.381112 0 0.25599398 42.234595 45.103092 0.74047865 0 0 0 0 0 0 0 0 0 0 0 0
1000 3.0049838 16229 0.24126663 42.1935 42.508307 10.275827 10158 39 39 0 970 90440 23409 16229 0.90908897 0 5 0
2000 2.1454941 16282 0.24068985 42.18192 42.406686 8.3551805 10114 39 39 0 1766 90440 23409 16282 0.91411077 0 5 0
3000 1.7194781 16297 0.24043991 42.182513 42.362649 9.0500757 10293 35 35 0 2454 90440 23409 16297 0.91613392 0 5 0
4000 1.9700398 16189 0.24063155 42.190029 42.396415 11.380017 10043 27 28 0 3083 90440 23409 16189 0.92535621 0 5 0
5000 2.2262815 16301 0.24071918 42.178625 42.411854 10.059361 10169 29 29 0 3730 90440 23409 16301 0.92861813 0 5 0
6000 2.2291507 16232 0.24063671 42.163793 42.397323 7.276142 9899 40 40 0 4445 90440 23409 16232 0.93275834 0 5 0
7000 2.7841984 16212 0.2407729 42.129641 42.421319 10.838996 10144 39 39 0 5148 90440 23409 16212 0.93329636 0 5 0
8000 3.1799279 16294 0.24102189 42.132053 42.465189 9.3204165 10068 31 31 0 5897 90440 23409 16294 0.93229671 0 5 0
9000 3.6101913 16308 0.24126974 42.130646 42.508857 9.9117415 10113 30 30 0 6716 90440 23409 16308 0.93107227 0 5 0
10000 4.1022731 16319 0.24161934 42.140689 42.570451 9.8799373 10205 34 34 0 7480 90440 23409 16319 0.92489062 0 5 0
11000 4.2790298 16269 0.2417129 42.138657 42.586936 10.14836 10002 37 37 2 8283 90440 23409 16269 0.92331323 0 10 0
12000 4.4299614 16323 0.24179577 42.137445 42.601537 11.786696 10068 37 37 1 9085 90440 23409 16323 0.93507414 0 10 0
13000 4.1451629 16263 0.24166347 42.143972 42.578227 7.5907603 9931 27 27 0 9967 90440 23409 16263 0.93269278 0 10 0
14000 5.0478037 16299 0.24222245 42.147894 42.676711 9.616107 10010 39 39 0 10856 90440 23409 16299 0.93132197 0 10 0
15000 4.9761523 16318 0.24208522 42.131222 42.652533 10.028283 10181 45 45 0 11744 90440 23409 16318 0.92914494 0 10 0
16000 5.0281695 16300 0.24209471 42.127445 42.654206 13.025657 10410 38 38 0 12657 90440 23409 16300 0.9322983 0 10 0
17000 6.0553511 16294 0.24278108 42.140765 42.775135 8.3971756 9898 31 31 0 13606 90440 23409 16294 0.93207864 0 10 0
18000 5.8870865 16349 0.24261355 42.128876 42.745619 9.4579631 9872 33 33 0 14510 90440 23409 16349 0.92911772 0 10 0
19000 6.19767 16260 0.2425369 42.082834 42.732114 10.497489 10068 37 37 0 15478 90440 23409 16260 0.9258668 0 10 0
20000 5.0588231 16330 0.24187554 42.085618 42.61559 9.3582625 9935 35 35 0 16470 90440 23409 16330 0.92639501 0 10 0
21000 6.4272357 16337 0.24285051 42.11404 42.78737 11.339789 9973 38 38 0 17363 90440 23409 16337 0.91830732 0 10 0
22000 6.3660195 16301 0.24257121 42.071243 42.738159 10.340724 9855 41 41 0 18300 90440 23409 16301 0.92735973 0 13 0
23000 5.9141967 16265 0.24245321 42.097787 42.71737 10.888473 9874 46 46 0 19294 90440 23409 16265 0.92958577 0 13 0
24000 7.3326392 16305 0.24333743 42.104977 42.873159 8.159204 9831 32 32 0 20337 90440 23409 16305 0.92498949 0 13 0
25000 7.8712498 16317 0.24365329 42.104202 42.928809 13.3861 10107 51 52 0 21343 90440 23409 16317 0.92093733 0 13 0
26000 7.6909446 16302 0.24354679 42.104327 42.910045 9.3626418 10101 45 45 0 22390 90440 23409 16302 0.91881953 0 13 0
27000 7.7525711 16296 0.24366541 42.11877 42.930944 9.3545743 9869 31 31 0 23458 90440 23409 16296 0.92276633 0 13 0
28000 7.8688635 16287 0.24363999 42.102109 42.926466 9.6143934 10023 37 37 0 24515 90440 23409 16287 0.92561111 0 13 0
29000 8.1116441 16298 0.24396496 42.13393 42.983722 9.8457785 10018 37 37 0 25571 90440 23409 16298 0.92896968 0 13 0
30000 7.5432306 16233 0.24356273 42.12261 42.912853 10.484741 9978 33 33 0 26714 90440 23409 16233 0.92084331 0 13 0
31000 7.9849131 16251 0.24395097 42.144741 42.981256 8.728219 9926 36 36 0 27794 90440 23409 16251 0.91985968 0 13 0
32000 8.8200827 16255 0.24461412 42.174088 43.098096 9.4477224 10118 31 31 0 28911 90440 23409 16255 0.92764378 0 13 0
33000 9.0462578 16266 0.24476502 42.176979 43.124682 9.0256529 9948 41 41 0 30010 90440 23409 16266 0.92496675 0 13 0
34000 11.115911 16273 0.24591951 42.163566 43.32809 10.89006 10013 35 35 1 31155 90440 23409 16273 0.9209011 0 59 0
35000 9.5287387 16348 0.24494148 42.157524 43.155772 9.1589486 10058 33 33 0 32262 90440 23409 16348 0.91852697 0 59 0
36000 8.6525285 16271 0.24454224 42.178976 43.085431 9.6115865 9936 39 39 0 33447 90440 23409 16271 0.91738637 0 59 0
37000 9.3244405 16237 0.24490339 42.172216 43.149062 11.766374 10057 49 49 0 34660 90440 23409 16237 0.92209291 0 59 0
38000 9.2952308 16272 0.24496591 42.186292 43.160078 9.5173055 10039 32 32 1 35916 90440 23409 16272 0.91915422 0 59 0
39000 10.12307 16312 0.24550971 42.195376 43.255888 10.678168 9989 43 43 1 37116 90440 23409 16312 0.92499805 0 59 0
40000 11.051307 16266 0.24606279 42.195578 43.353334 10.078073 9867 37 37 0 38363 90440 23409 16266 0.9229677 0 59 0
41000 10.400692 16313 0.24566631 42.193882 43.283479 9.2181368 10040 39 39 0 39612 90440 23409 16313 0.92195363 0 59 0
42000 9.0646248 16277 0.24487452 42.194348 43.143975 8.2531891 10173 37 37 0 40821 90440 23409 16277 0.91941681 0 59 0
43000 9.8646504 16258 0.24532194 42.189366 43.222806 10.34072 10037 37 37 0 42068 90440 23409 16258 0.91264623 0 59 0
44000 10.441846 16296 0.24553053 42.165649 43.259556 12.147264 10299 31 31 0 43277 90440 23409 16296 0.91599514 0 59 0
45000 10.458302 16261 0.24544296 42.148496 43.244128 8.8999271 10208 26 26 0 44492 90440 23409 16261 0.91716968 0 59 0
46000 10.117248 16269 0.24511999 42.127322 43.187224 9.0053892 10132 36 36 0 45691 90440 23409 16269 0.91327538 0 59 0
47000 9.6094606 16247 0.24494761 42.150148 43.156854 11.844069 9936 43 43 0 46930 90440 23409 16247 0.91513796 0 59 0
48000 12.544005 16250 0.24661869 42.137144 43.451278 14.420938 9999 43 43 0 48233 90440 23409 16250 0.91642267 0 59 0
49000 11.948168 16257 0.24628573 42.140902 43.392614 8.653163 10093 38 38 0 49517 90440 23409 16257 0.90604463 0 59 0
50000 12.315254 16233 0.2466687 42.16992 43.460089 9.1125366 10050 39 39 0 50806 90440 23409 16233 0.90490961 0 59 0
51000 11.364234 16265 0.24609762 42.168932 43.359471 9.8694192 9860 37 37 0 52074 90440 23409 16265 0.909623 0 59 0
52000 10.627787 16231 0.24578998 42.191882 43.305269 10.869196 9877 33 33 0 53420 90440 23409 16231 0.90981484 0 59 0
53000 11.809423 16261 0.24656286 42.204263 43.441441 13.38992 9828 39 39 0 54766 90440 23409 16261 0.90710704 0 59 0
54000 10.63283 16251 0.24554697 42.148537 43.262453 10.181102 10098 29 29 0 56085 90440 23409 16251 0.90800902 0 59 0
55000 10.803075 16254 0.24562954 42.145251 43.277002 7.9856721 9964 26 26 0 57451 90440 23409 16254 0.90786888 0 59 0
56000 10.844803 16274 0.2458735 42.183861 43.319983 9.9407901 10031 36 36 0 58695 90440 23409 16274 0.89847818 0 59 0
57000 10.623373 16218 0.24574768 42.184891 43.297816 12.015798 9792 43 43 0 60012 90440 23409 16218 0.90944412 0 59 0
58000 10.679964 16262 0.24583434 42.194231 43.313084 11.651504 9776 45 45 0 61295 90440 23409 16262 0.90511687 0 59 0
59000 10.35931 16278 0.24567475 42.199706 43.284967 10.599471 10069 30 30 1 62506 90440 23409 16278 0.90064789 0 59 0
60000 12.034078 16208 0.24677506 42.218115 43.478828 10.511331 9936 31 31 0 63825 90440 23409 16208 0.90250225 0 59 0
61000 11.915856 16208 0.24666669 42.211407 43.459734 11.982448 10227 34 34 0 65189 90440 23409 16208 0.90006711 0 59 0
62000 11.388719 16269 0.24635856 42.212341 43.405445 8.4867265 9991 29 29 0 66502 90440 23409 16269 0.90418592 0 59 0
63000 10.513272 16243 0.24593174 42.228854 43.330244 9.2843187 10133 39 39 0 67836 90440 23409 16243 0.91282932 0 59 0
64000 11.465761 16203 0.24638806 42.209467 43.410642 8.8736184 10023 38 39 0 69216 90440 23409 16203 0.90742452 0 59 0
65000 9.9798452 16292 0.24547666 42.204557 43.250064 12.475112 9766 36 36 0 70559 90440 23409 16292 0.90921789 0 59 0
66000 12.184131 16177 0.2468175 42.209872 43.486305 8.4442171 10049 33 33 0 71851 90440 23409 16177 0.90771921 0 59 0
67000 11.845593 16250 0.24645726 42.181869 43.422835 8.2359478 9678 33 33 0 73241 90440 23409 16250 0.9069606 0 59 0
68000 10.616202 16220 0.2456631 42.17074 43.282913 7.6844013 9954 28 28 0 74707 90440 23409 16220 0.90830051 0 59 0
69000 11.708481 16174 0.246425 42.190549 43.417152 9.6601201 10019 34 34 0 76143 90440 23409 16174 0.90853993 0 59 0
70000 13.149115 16175 0.24730865 42.195313 43.57284 7.8721271 9794 31 31 0 77520 90440 23409 16175 0.90460472 0 59 0
71000 12.420743 16210 0.24676949 42.176626 43.477846 10.370102 9886 36 36 0 78901 90440 23409 16210 0.90567986 0 59 0
72000 13.424842 16222 0.24726211 42.158228 43.56464 10.837461 9892 36 36 0 80221 90440 23409 16222 0.90896224 0 59 0
73000 12.983985 16177 0.2470467 42.16646 43.526687 7.3246122 9975 31 31 0 81578 90440 23409 16177 0.90172487 0 59 0
74000 11.873207 16217 0.24644792 42.177329 43.421189 11.257297 9937 40 40 0 83004 90440 23409 16217 0.9069376 0 59 0
75000 11.719883 16219 0.24631358 42.169723 43.397521 8.2731229 10281 28 28 0 84503 90440 23409 16219 0.90087971 0 59 0
76000 11.648725 16220 0.24633978 42.181793 43.402136 11.364684 9961 31 31 0 85910 90440 23409 16220 0.90445518 0 59 0
77000 11.315159 16169 0.24618794 42.189986 43.375384 10.639734 9853 35 35 0 87272 90440 23409 16169 0.90967681 0 59 0
78000 11.026923 16262 0.2459278 42.17435 43.329551 9.3461335 10068 36 36 0 88765 90440 23409 16262 0.90250223 0 59 0
79000 11.923316 16205 0.24657893 42.195162 43.444272 8.3200696 9914 29 29 1 90167 90440 23409 16205 0.89792511 0 59 0
80000 12.28723 16233 0.246793 42.194756 43.481989 9.8594414 9898 41 41 0 91543 90440 23409 16233 0.89920018 0 59 0
81000 12.657271 16252 0.24700788 42.193847 43.519847 9.5478351 10045 38 38 0 92888 90440 23409 16252 0.90032559 0 59 0
82000 10.415047 16249 0.24563975 42.1877 43.2788 8.4502786 9729 30 30 0 94355 90440 23409 16249 0.90516702 0 59 0
83000 13.038087 16237 0.24718282 42.184776 43.550671 10.813181 9834 48 48 0 95792 90440 23409 16237 0.90160106 0 59 0
84000 11.994728 16253 0.24664911 42.200046 43.456637 10.907779 9632 48 48 0 97203 90440 23409 16253 0.90075335 0 59 0
85000 11.560331 16253 0.24626448 42.177788 43.38887 9.8019244 9637 27 27 0 98581 90440 23409 16253 0.90204694 0 59 0
86000 12.294346 16214 0.2466302 42.165325 43.453305 8.4731319 9776 30 30 1 100022 90440 23409 16214 0.8976977 0 59 0
87000 11.439691 16227 0.24593937 42.133145 43.331589 11.138197 9408 35 35 0 101446 90440 23409 16227 0.90678956 0 59 0
88000 12.269451 16218 0.24642733 42.132191 43.417562 9.2022903 9627 28 28 0 102882 90440 23409 16218 0.89354177 0 59 0
89000 11.630085 16281 0.24612243 42.145453 43.363843 9.4906584 9713 45 45 0 104390 90440 23409 16281 0.90416008 0 59 0
90000 12.068036 16226 0.24636746 42.142743 43.407014 13.55402 9721 30 30 1 105786 90440 23409 16226 0.90693102 0 59 0
91000 11.423374 16180 0.24593187 42.133533 43.330267 11.12182 9490 33 33 0 107319 90440 23409 16180 0.89908597 0 59 0
92000 11.759782 16253 0.24594059 42.099827 43.331804 8.5785997 9808 36 37 0 108822 90440 23409 16253 0.90323757 0 59 0
93000 12.860286 16203 0.24670999 42.120094 43.467362 9.7352723 9342 32 32 0 110320 90440 23409 16203 0.89837993 0 59 0
94000 12.65067 16201 0.24650827 42.106515 43.431823 11.726981 9708 37 37 0 111721 90440 23409 16201 0.90153973 0 59 0
95000 12.938398 16301 0.2467405 42.117288 43.472739 13.229664 9490 42 42 0 113106 90440 23409 16301 0.89222858 0 59 0
96000 13.125 16243 0.24672164 42.094416 43.469416 11.492955 9666 43 43 0 114441 90440 23409 16243 0.89369219 0 59 0
97000 14.779601 16270 0.24775001 42.102263 43.650602 9.7132418 9517 31 31 0 115831 90440 23409 16270 0.88854614 0 59 0
98000 15.427308 16254 0.24808383 42.093224 43.709418 11.821283 9510 36 37 0 117263 90440 23409 16254 0.88982996 0 59 0
99000 16.212104 16248 0.24846852 42.078784 43.777195 10.786344 9807 37 37 0 118659 90440 23409 16248 0.89136464 0 59 0
100000 15.896 16229 0.24841013 42.101613 43.766908 10.836017 9640 28 28 0 120081 90440 23409 16229 0.90018614 0 59 0
Loop time of 53.3125 on 8 procs for 100000 steps with 73999 atoms
Pair time (%) = 0.559711 (1.04987)
Neigh time (%) = 3.94754 (7.40453)
Comm time (%) = 3.46588 (6.50106)
Outpt time (%) = 0.00787139 (0.0147646)
Other time (%) = 45.3315 (85.0298)
Nlocal: 9249.88 ave 10477 max 8359 min
Histogram: 1 2 1 1 0 1 0 0 0 2
Nghost: 54 ave 65 max 45 min
Histogram: 2 1 0 0 1 1 0 2 0 1
Neighs: 316.875 ave 387 max 170 min
Histogram: 1 0 0 0 1 0 0 4 0 2
Total # of neighbors = 2535
Ave neighs/atom = 0.0342572
Neighbor list builds = 5000
Dangerous builds = 0
Please see the log.cite file for references relevant to this simulation

286
examples/ASPHERE/dimer/log.dimer.1Feb14 Normal file
View File

@ -0,0 +1,286 @@
LAMMPS (1 Feb 2014)
# SRD diffusion demo - dimer particles
units lj
atom_style sphere
atom_modify map array first big
dimension 2
# read in clusters of rigid bodies
fix molprop all property/atom mol
read_data data.dimer fix molprop NULL Molecules
orthogonal box = (-9.34165 -9.34165 -0.5) to (9.34165 9.34165 0.5)
4 by 2 by 1 MPI processor grid
reading atoms ...
200 atoms
set type 1 mass 1.0
200 settings made for mass
group big type 1
200 atoms in group big
velocity big create 1.44 87287 loop geom
# equilibrate big particles
pair_style soft 1.12
pair_coeff 1 1 0.0
pair_coeff 2 2 0.0 0.0
pair_coeff 1 2 0.0 0.0
variable prefactor equal ramp(0,60)
fix soft all adapt 1 pair soft a * * v_prefactor
fix 1 big rigid molecule
100 rigid bodies with 200 atoms
fix 2 all enforce2d
#dump 1 all atom 10 dump.dimer.equil
thermo 100
run 1000
Memory usage per processor = 3.1029 Mbytes
Step Temp E_pair E_mol TotEng Press
0 3.4028231 0 0 0.83369167 -0.55065517
100 9.5167872 2.392105 0 4.7237178 2.8319556
200 13.564465 3.0352634 0 6.3585572 3.6388732
300 13.133051 4.3835112 0 7.6011086 4.060051
400 14.584346 5.5141254 0 9.0872901 4.4231056
500 15.331515 6.6554832 0 10.411704 4.847642
600 16.953755 7.4794317 0 11.633102 5.2020696
700 16.503256 8.8572339 0 12.900532 5.6694613
800 17.006131 9.877343 0 14.043845 5.5483363
900 17.305927 11.081594 0 15.321546 6.2908201
1000 18.122491 12.126462 0 16.566472 5.7536055
Loop time of 0.185949 on 8 procs for 1000 steps with 200 atoms
Pair time (%) = 0.00187448 (1.00806)
Neigh time (%) = 0.000561714 (0.30208)
Comm time (%) = 0.0652371 (35.0833)
Outpt time (%) = 0.00209856 (1.12857)
Other time (%) = 0.116177 (62.4779)
Nlocal: 25 ave 30 max 20 min
Histogram: 1 0 1 1 0 0 3 1 0 1
Nghost: 28.375 ave 36 max 22 min
Histogram: 1 0 1 2 1 0 2 0 0 1
Neighs: 33.125 ave 45 max 22 min
Histogram: 1 0 2 1 0 1 1 0 1 1
Total # of neighbors = 265
Ave neighs/atom = 1.325
Neighbor list builds = 99
Dangerous builds = 96
#undump 1
unfix soft
unfix 1
unfix 2
# add small particles as hi density lattice
region plane block INF INF INF INF -0.001 0.001 units box
lattice sq 85.0
Lattice spacing in x,y,z = 0.108465 0.108465 0.108465
create_atoms 2 region plane
Created 29929 atoms
set type 2 mass 0.1
29929 settings made for mass
group small type 2
29929 atoms in group small
velocity small create 1.0 593849 loop geom
# delete overlaps
# must set 1-2 cutoff to non-zero value
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 0.5
delete_atoms overlap 0.5 small big
Deleted 12759 atoms, new total = 17370
# SRD run
reset_timestep 0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
communicate multi group big vel yes
neigh_modify include big
# no pairwise interactions with small particles
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 0.0
# use fix SRD to push small particles out from inside big ones
# if comment out, big particles won't see SRD particles
timestep 0.001
fix 1 big rigid molecule
100 rigid bodies with 200 atoms
fix 2 small srd 20 big 1.0 0.25 49894 shift yes 54979 search 0.2 collision slip inside ignore overlap yes
fix 3 all enforce2d
# diagnostics
compute tbig big temp/sphere
variable pebig equal pe*atoms/count(big)
variable ebig equal etotal*atoms/count(big)
thermo_style custom step temp f_2[8] etotal v_pebig v_ebig press f_2[1] f_2[2] f_2[3] f_2[4] f_2[5] f_2[6] f_2[7] f_2[8] f_2[9] f_2[10] f_2[11] f_2[12]
thermo_modify temp tbig
WARNING: Temperature for thermo pressure is not for group all (../thermo.cpp:439)
thermo 1000
#dump 1 all atom 1000 dump.dimer
#dump 1 all image 1000 image.*.jpg type type zoom 1.6
#dump_modify 1 pad 6 adiam 1 1 adiam 2 0.2
run 100000
SRD info:
SRD/big particles = 17170 200
big particle diameter max/min = 1 1
SRD temperature & lamda = 1 0.0632456
SRD max distance & max velocity = 0.252982 12.6491
SRD grid counts: 75 75 1
SRD grid size: request, actual (xyz) = 0.25, 0.249111 0.249111 1
SRD per actual grid cell = 5.5499
SRD viscosity = 0.235551
big/SRD mass density ratio = 0.142367
# of rescaled SRD velocities = 0
ave/max small velocity = 4.20347 7.72323
ave/max big velocity = 2.64047 6.60911
Memory usage per processor = 7.0395 Mbytes
Step Temp 2[8] TotEng pebig ebig Press 2[1] 2[2] 2[3] 2[4] 2[5] 2[6] 2[7] 2[8] 2[9] 2[10] 2[11] 2[12]
0 6.9744766 0 0.30175295 21.011259 26.207244 2.7219331 0 0 0 0 0 0 0 0 0 0 0 0
1000 1.1941672 3369 0.24718395 20.578272 21.467926 9.8944086 13658 56 56 2 664 40320 5625 3369 1.1358997 0 4 0
2000 1.3247909 3390 0.24717547 20.48022 21.467189 6.5341926 13389 48 48 3 1988 40320 5625 3390 1.1814916 0 144 0
3000 1.386904 3355 0.24769415 20.478993 21.512237 6.1207029 13507 45 45 0 2871 40320 5625 3355 1.1657697 0 144 0
4000 1.1567347 3376 0.24573733 20.48052 21.342287 3.5540965 13475 46 46 0 3731 40320 5625 3376 1.1485984 0 144 0
5000 1.2565155 3376 0.24620928 20.447172 21.383276 2.5808161 13267 55 55 0 4559 40320 5625 3376 1.1778096 0 144 0
6000 1.1489985 3375 0.24520477 20.440031 21.296035 7.8204262 13362 59 59 0 5354 40320 5625 3375 1.1530517 0 144 0
7000 1.1608248 3399 0.24587989 20.489854 21.354669 7.9622779 13333 59 59 2 6254 40320 5625 3399 1.1809325 0 144 0
8000 1.1749441 3374 0.24700359 20.576929 21.452262 6.7498405 13495 64 64 0 7134 40320 5625 3374 1.1445035 0 144 0
9000 1.1285085 3405 0.24536183 20.468936 21.309675 5.8959937 13333 60 60 0 8008 40320 5625 3405 1.162836 0 144 0
10000 1.1462675 3425 0.24564758 20.480523 21.334492 7.0047543 13487 45 45 0 8967 40320 5625 3425 1.1881074 0 144 0
11000 1.2854766 3388 0.246088 20.415062 21.372742 7.8508136 13491 60 60 1 9833 40320 5625 3388 1.1698514 0 144 0
12000 1.2926536 3360 0.24602616 20.404345 21.367372 7.9286141 13142 64 64 1 11029 40320 5625 3360 1.1658039 0 161 0
13000 1.2289767 3412 0.24595657 20.44574 21.361328 2.2809592 13162 56 56 1 12018 40320 5625 3412 1.1665531 0 161 0
14000 1.3518235 3380 0.24629055 20.383226 21.390334 8.0062608 13310 59 59 7 13146 40320 5625 3380 1.1733296 0 161 0
15000 1.2971143 3355 0.24671134 20.46053 21.42688 9.7871945 12967 60 60 3 14084 40320 5625 3355 1.1703643 0 161 0
16000 1.4160948 3384 0.2478594 20.471598 21.526589 7.3141215 13123 56 56 1 14929 40320 5625 3384 1.1666782 0 161 0
17000 1.3172038 3391 0.247377 20.503375 21.484692 6.1541908 12957 54 54 0 15778 40320 5625 3391 1.1639931 0 161 0
18000 0.98748293 3369 0.2443499 20.486114 21.221789 4.0981146 13355 50 50 1 17031 40320 5625 3369 1.2083414 0 192 0
19000 1.4619261 3367 0.24867817 20.508564 21.597699 11.68547 13519 58 58 0 18334 40320 5625 3367 1.1667505 0 192 0
20000 1.44087 3391 0.24747923 20.420123 21.493571 4.0409936 13594 67 67 0 19303 40320 5625 3391 1.1869387 0 192 0
21000 1.1155365 3377 0.24474418 20.424957 21.256032 8.4985692 13423 52 52 0 20245 40320 5625 3377 1.1726608 0 192 0
22000 1.2045785 3382 0.24466159 20.351448 21.248859 7.0518674 13271 62 62 0 21128 40320 5625 3382 1.1658557 0 192 0
23000 1.2634616 3378 0.24571326 20.398918 21.340197 6.8382109 13246 61 61 1 22067 40320 5625 3378 1.170196 0 192 0
24000 1.365723 3375 0.24749053 20.477089 21.494552 5.9584152 13017 46 46 1 23017 40320 5625 3375 1.1841192 0 192 0
25000 1.1421045 3367 0.24497142 20.4249 21.275768 10.381661 13049 66 66 1 23907 40320 5625 3367 1.1856767 0 192 0
26000 1.4008939 3340 0.24678624 20.389719 21.433385 7.3799658 12971 53 53 5 24850 40320 5625 3340 1.1713233 0 192 0
27000 1.1835022 3379 0.24522732 20.416284 21.297993 7.2251705 12717 46 46 1 25723 40320 5625 3379 1.1798333 0 192 0
28000 1.0496668 3398 0.24461816 20.463086 21.245087 6.5485338 13385 60 60 0 26990 40320 5625 3398 1.1765205 0 192 0
29000 1.1779927 3453 0.24501083 20.401586 21.27919 4.8540048 13157 77 77 2 28023 40320 5625 3453 1.1741595 0 192 0
30000 1.2277637 3375 0.24554215 20.410652 21.325336 10.087705 12939 72 72 1 28855 40320 5625 3375 1.1627414 0 192 0
31000 1.1575066 3397 0.24504213 20.419566 21.281909 6.4849648 12977 66 66 0 29745 40320 5625 3397 1.1575018 0 192 0
32000 1.1752733 3381 0.2451073 20.41199 21.287569 7.2784509 12971 48 48 0 30561 40320 5625 3381 1.1790395 0 192 0
33000 1.1743895 3370 0.24485514 20.390748 21.265669 9.084299 12810 64 64 0 31342 40320 5625 3370 1.1429061 0 192 0
34000 1.2975876 3377 0.24608262 20.405573 21.372276 10.624304 12610 60 60 1 32140 40320 5625 3377 1.1685718 0 192 0
35000 1.7701252 3356 0.2497356 20.370794 21.689537 7.107766 12983 45 45 8 33862 40320 5625 3356 1.1720964 0 242 0
36000 1.2711521 3393 0.24675748 20.483878 21.430887 8.7087273 13175 60 60 0 34920 40320 5625 3393 1.1796212 0 242 0
37000 1.1140649 3373 0.24450996 20.405711 21.23569 8.0025138 12964 50 50 1 35770 40320 5625 3373 1.1772234 0 242 0
38000 1.3128798 3379 0.24583813 20.372946 21.351041 5.0554655 12939 56 56 1 36679 40320 5625 3379 1.1702949 0 242 0
39000 1.2093703 3365 0.24536578 20.409037 21.310018 5.6243735 12873 62 62 3 37511 40320 5625 3365 1.1751064 0 242 0
40000 1.2095402 3365 0.24485106 20.364207 21.265315 3.9565189 12876 54 54 1 38503 40320 5625 3365 1.1637516 0 242 0
41000 1.1146238 3350 0.24476111 20.427107 21.257502 6.5033046 12656 53 53 0 39313 40320 5625 3350 1.1553197 0 242 0
42000 1.2302245 3370 0.24452048 20.320087 21.236604 10.7363 12676 52 52 1 40108 40320 5625 3370 1.1492163 0 242 0
43000 1.2749689 3404 0.24511103 20.338041 21.287893 9.18604 13057 72 72 1 41014 40320 5625 3404 1.1572511 0 242 0
44000 1.1989211 3385 0.24450633 20.342179 21.235375 7.6714085 13086 52 53 1 41840 40320 5625 3385 1.1530041 0 242 0
45000 1.6465972 3370 0.249994 20.485264 21.711979 7.4908607 12894 58 58 7 43545 40320 5625 3370 1.1978571 0 361 0
46000 1.3082316 3412 0.24665637 20.447473 21.422106 8.9900563 13136 53 53 4 44510 40320 5625 3412 1.1686573 0 361 0
47000 1.2163258 3355 0.24493724 20.366636 21.272799 6.9309438 13020 56 56 0 45414 40320 5625 3355 1.1620943 0 361 0
48000 1.349966 3362 0.24683156 20.431596 21.437321 8.8527164 13044 64 64 1 46260 40320 5625 3362 1.1683452 0 361 0
49000 1.3903769 3347 0.24726403 20.439051 21.474881 5.9250552 13103 58 59 2 49229 40320 5625 3347 1.1924538 0 578 0
50000 1.2737918 3349 0.24563843 20.384723 21.333698 7.2419874 12892 41 41 1 50225 40320 5625 3349 1.1803233 0 578 0
51000 1.2229799 3355 0.24504936 20.371417 21.282537 8.5556146 12912 63 63 0 51096 40320 5625 3355 1.1333843 0 578 0
52000 1.2326547 3386 0.24596384 20.443632 21.361959 6.1757404 12790 66 66 1 52114 40320 5625 3386 1.1685608 0 578 0
53000 1.3695441 3389 0.24696147 20.428294 21.448604 8.7464732 12997 60 60 1 53115 40320 5625 3389 1.1631228 0 578 0
54000 1.0711204 3400 0.24322154 20.325806 21.12379 8.7228921 13005 56 56 0 53913 40320 5625 3400 1.1250986 0 578 0
55000 1.3412565 3360 0.24617922 20.381429 21.380665 6.9489559 12960 46 46 0 55249 40320 5625 3360 1.176467 0 578 0
56000 1.3910761 3355 0.2467942 20.397725 21.434077 6.6762734 12966 65 65 2 56653 40320 5625 3355 1.1715559 0 578 0
57000 1.1177307 3388 0.24424738 20.380176 21.212885 5.5073089 13153 62 62 1 57550 40320 5625 3388 1.1576307 0 578 0
58000 1.3136279 3365 0.24661231 20.439627 21.418279 9.1390284 12799 52 52 1 58366 40320 5625 3365 1.159026 0 578 0
59000 1.2667621 3374 0.24556104 20.383238 21.326976 3.4246418 13127 54 54 1 59170 40320 5625 3374 1.1461442 0 578 0
60000 1.1903121 3326 0.24492127 20.384629 21.271412 6.9439977 13078 47 47 2 60105 40320 5625 3326 1.1665886 0 578 0
61000 1.209313 3362 0.2455341 20.423699 21.324637 9.1212002 13225 61 61 4 61358 40320 5625 3362 1.1725894 0 578 0
62000 1.1720074 3366 0.24476301 20.384522 21.257667 7.7592884 13187 53 54 2 62257 40320 5625 3366 1.1599665 0 578 0
63000 1.2155508 3386 0.24486685 20.361101 21.266686 5.7864126 13090 42 42 0 63137 40320 5625 3386 1.1724234 0 578 0
64000 1.1967326 3379 0.2451834 20.402613 21.294179 4.8815049 13033 48 48 4 64152 40320 5625 3379 1.1818756 0 578 0
65000 1.1732378 3386 0.24500733 20.404824 21.278887 11.381026 12998 59 59 0 65071 40320 5625 3386 1.1625543 0 578 0
66000 1.1960739 3385 0.24407516 20.306852 21.197927 12.17249 12994 61 61 0 65941 40320 5625 3385 1.1826516 0 578 0
67000 1.3468331 3366 0.24657491 20.41164 21.415031 7.8115623 12894 62 62 1 66788 40320 5625 3366 1.1750424 0 578 0
68000 1.0742924 3369 0.24347736 20.345661 21.146009 9.5240252 13065 68 68 2 67665 40320 5625 3369 1.1561575 0 578 0
69000 1.1411254 3380 0.2450544 20.432836 21.282975 8.0392218 13134 49 49 1 68471 40320 5625 3380 1.1405485 0 578 0
70000 1.1585312 3388 0.24435206 20.358871 21.221977 7.5303626 13113 59 59 0 69320 40320 5625 3388 1.1726416 0 578 0
71000 1.3418857 3382 0.24598802 20.364355 21.36406 9.2840877 13265 68 68 0 70522 40320 5625 3382 1.1696767 0 578 0
72000 1.1547857 3388 0.24328952 20.26938 21.129695 7.9643121 13270 50 50 1 71384 40320 5625 3388 1.1647353 0 578 0
73000 1.1114283 3385 0.24397276 20.36102 21.189034 13.168935 13110 62 62 0 72335 40320 5625 3385 1.1625358 0 578 0
74000 1.6120075 3370 0.25009739 20.520013 21.720959 12.065782 12902 51 51 9 73743 40320 5625 3370 1.1873738 0 578 0
75000 1.3271087 3353 0.24600453 20.376798 21.365493 7.5902315 13076 51 51 1 74746 40320 5625 3353 1.1706259 0 578 0
76000 1.1908626 3386 0.2451628 20.405197 21.292389 6.8084086 13124 58 58 1 75702 40320 5625 3386 1.1765809 0 578 0
77000 1.042536 3391 0.24320847 20.345966 21.122655 6.4883579 12846 57 57 0 76930 40320 5625 3391 1.1769829 0 578 0
78000 1.3791243 3355 0.24628683 20.362563 21.390011 9.2397346 12897 56 56 12 78068 40320 5625 3355 1.1905167 0 578 0
79000 1.2058397 3391 0.2453076 20.406614 21.304965 9.7949819 13020 51 51 1 79065 40320 5625 3391 1.1777257 0 578 0
80000 1.0634737 3383 0.24440921 20.434652 21.22694 7.5296946 12973 61 61 2 80029 40320 5625 3383 1.1648017 0 578 0
81000 1.3768966 3373 0.24730623 20.452758 21.478546 7.1508584 13297 56 56 1 81172 40320 5625 3373 1.1740603 0 578 0
82000 1.3348382 3375 0.24680945 20.440946 21.4354 9.3251946 12822 56 56 2 82447 40320 5625 3375 1.1573749 0 578 0
83000 1.2042501 3393 0.24484972 20.368032 21.265198 7.5283729 12990 51 51 0 83344 40320 5625 3393 1.1833162 0 578 0
84000 1.1985223 3405 0.24537841 20.418216 21.311115 7.848856 13189 60 60 0 84147 40320 5625 3405 1.173911 0 578 0
85000 1.1811593 3391 0.24418657 20.32764 21.207604 8.3169438 13303 51 51 0 85022 40320 5625 3391 1.1586027 0 578 0
86000 1.2595689 3385 0.24556774 20.389179 21.327558 8.3808674 12966 56 56 1 85962 40320 5625 3385 1.1617838 0 578 0
87000 1.1419322 3356 0.24399697 20.340397 21.191137 6.8911116 13166 63 63 0 86753 40320 5625 3356 1.1729274 0 578 0
88000 1.275219 3380 0.24586129 20.403015 21.353053 6.838563 12963 56 56 0 87729 40320 5625 3380 1.1511968 0 578 0
89000 1.2794858 3366 0.24507826 20.33183 21.285047 9.3370019 12889 64 64 1 88575 40320 5625 3366 1.1660311 0 578 0
90000 1.0549593 3402 0.24289618 20.309589 21.095533 4.9944605 13076 59 59 1 89560 40320 5625 3402 1.1744335 0 578 0
91000 1.214812 3402 0.24442304 20.323106 21.228141 7.5681019 13229 44 44 1 90440 40320 5625 3402 1.1696689 0 578 0
92000 1.2019714 3374 0.24469176 20.356011 21.25148 7.6728432 13331 64 64 2 91243 40320 5625 3374 1.1678246 0 578 0
93000 1.0964004 3364 0.24364709 20.343931 21.16075 5.9650235 13086 55 55 0 92088 40320 5625 3364 1.1570838 0 578 0
94000 1.1722133 3393 0.24332541 20.259513 21.132812 7.1743779 13206 66 66 0 92858 40320 5625 3393 1.154328 0 578 0
95000 1.2413711 3373 0.24549874 20.396744 21.321566 8.131371 12967 61 61 1 93789 40320 5625 3373 1.1635681 0 578 0
96000 1.1702409 3363 0.24528289 20.430989 21.302819 11.399295 12776 63 63 0 94571 40320 5625 3363 1.1367365 0 578 0
97000 1.0696537 3368 0.242673 20.279258 21.07615 8.537175 13054 55 55 0 95403 40320 5625 3368 1.1648141 0 578 0
98000 1.1635952 3369 0.24486659 20.399785 21.266664 8.6714443 12750 70 70 0 96296 40320 5625 3369 1.1686272 0 578 0
99000 1.1286255 3397 0.24343772 20.30174 21.142566 6.7153212 12810 44 45 1 97021 40320 5625 3397 1.137818 0 578 0
100000 1.2847818 3367 0.24509282 20.329149 21.286312 8.7486629 13260 64 64 0 97855 40320 5625 3367 1.1597138 0 578 0
Loop time of 20.8136 on 8 procs for 100000 steps with 17370 atoms
Pair time (%) = 0.322913 (1.55145)
Neigh time (%) = 1.11753 (5.36921)
Comm time (%) = 1.72335 (8.27991)
Outpt time (%) = 0.00594518 (0.0285638)
Other time (%) = 17.6439 (84.7709)
Nlocal: 2171.25 ave 2428 max 1747 min
Histogram: 1 0 0 0 1 1 2 1 0 2
Nghost: 61.25 ave 69 max 54 min
Histogram: 1 1 0 1 1 1 0 2 0 1
Neighs: 175.625 ave 225 max 132 min
Histogram: 1 1 1 1 1 0 1 0 0 2
Total # of neighbors = 1405
Ave neighs/atom = 0.0808866
Neighbor list builds = 5156
Dangerous builds = 1
Please see the log.cite file for references relevant to this simulation

237
examples/ASPHERE/dimer/log.dimer.mp.1Feb14 Normal file
View File

@ -0,0 +1,237 @@
LAMMPS (1 Feb 2014)
# SRD viscosity demo - dimer particles
units lj
atom_style sphere
atom_modify map array first big
dimension 2
# read in clusters of rigid bodies
fix molprop all property/atom mol
read_data data.dimer fix molprop NULL Molecules
orthogonal box = (-9.34165 -9.34165 -0.5) to (9.34165 9.34165 0.5)
4 by 2 by 1 MPI processor grid
reading atoms ...
200 atoms
set type 1 mass 1.0
200 settings made for mass
group big type 1
200 atoms in group big
velocity big create 1.44 87287 loop geom
# equilibrate big particles
pair_style soft 1.12
pair_coeff 1 1 0.0
pair_coeff 2 2 0.0 0.0
pair_coeff 1 2 0.0 0.0
variable prefactor equal ramp(0,60)
fix soft all adapt 1 pair soft a * * v_prefactor
fix 1 big rigid molecule
100 rigid bodies with 200 atoms
fix 2 all enforce2d
#dump 1 all atom 10 dump.dimer.equil
thermo 100
run 1000
Memory usage per processor = 3.1029 Mbytes
Step Temp E_pair E_mol TotEng Press
0 3.4028231 0 0 0.83369167 -0.55065517
100 9.5167872 2.392105 0 4.7237178 2.8319556
200 13.564465 3.0352634 0 6.3585572 3.6388732
300 13.133051 4.3835112 0 7.6011086 4.060051
400 14.584346 5.5141254 0 9.0872901 4.4231056
500 15.331515 6.6554832 0 10.411704 4.847642
600 16.953755 7.4794317 0 11.633102 5.2020696
700 16.503256 8.8572339 0 12.900532 5.6694613
800 17.006131 9.877343 0 14.043845 5.5483363
900 17.305927 11.081594 0 15.321546 6.2908201
1000 18.122491 12.126462 0 16.566472 5.7536055
Loop time of 0.0603173 on 8 procs for 1000 steps with 200 atoms
Pair time (%) = 0.00182396 (3.02395)
Neigh time (%) = 0.000557959 (0.92504)
Comm time (%) = 0.0127766 (21.1823)
Outpt time (%) = 0.000280827 (0.465583)
Other time (%) = 0.0448779 (74.4031)
Nlocal: 25 ave 30 max 20 min
Histogram: 1 0 1 1 0 0 3 1 0 1
Nghost: 28.375 ave 36 max 22 min
Histogram: 1 0 1 2 1 0 2 0 0 1
Neighs: 33.125 ave 45 max 22 min
Histogram: 1 0 2 1 0 1 1 0 1 1
Total # of neighbors = 265
Ave neighs/atom = 1.325
Neighbor list builds = 99
Dangerous builds = 96
#undump 1
unfix soft
unfix 1
unfix 2
# add small particles as hi density lattice
region plane block INF INF INF INF -0.001 0.001 units box
lattice sq 85.0
Lattice spacing in x,y,z = 0.108465 0.108465 0.108465
create_atoms 2 region plane
Created 29929 atoms
set type 2 mass 0.1
29929 settings made for mass
group small type 2
29929 atoms in group small
velocity small create 1.0 593849 loop geom
# delete overlaps
# must set 1-2 cutoff to non-zero value
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 0.5
delete_atoms overlap 0.5 small big
Deleted 12759 atoms, new total = 17370
# SRD run
reset_timestep 0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
communicate multi group big vel yes
neigh_modify include big
# no pairwise interactions with small particles
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 0.0
# use fix SRD to push small particles out from inside big ones
# if comment out, big particles won't see SRD particles
timestep 0.001
fix 1 big rigid molecule
100 rigid bodies with 200 atoms
fix 2 small srd 20 big 1.0 0.25 49894 shift yes 54979 search 0.2 collision slip inside ignore overlap yes
fix 3 small viscosity 10 x y 50
fix 4 all enforce2d
# diagnostics
compute tbig big temp/sphere
variable pebig equal pe*atoms/count(big)
variable ebig equal etotal*atoms/count(big)
thermo_style custom step temp f_2[8] etotal v_pebig v_ebig press f_2[1] f_2[2] f_2[3] f_2[4] f_2[5] f_2[6] f_2[7] f_2[8] f_2[9] f_2[10] f_2[11] f_2[12]
thermo_modify temp tbig
WARNING: Temperature for thermo pressure is not for group all (../thermo.cpp:439)
thermo 1000
#dump 1 all atom 500 dump.dimer.mp
#dump 1 all image 500 image.*.jpg type type zoom 1.6
#dump_modify 1 pad 6 adiam 1 1 adiam 2 0.2
run 50000
SRD info:
SRD/big particles = 17170 200
big particle diameter max/min = 1 1
SRD temperature & lamda = 1 0.0632456
SRD max distance & max velocity = 0.252982 12.6491
SRD grid counts: 75 75 1
SRD grid size: request, actual (xyz) = 0.25, 0.249111 0.249111 1
SRD per actual grid cell = 5.5499
SRD viscosity = 0.235551
big/SRD mass density ratio = 0.142367
# of rescaled SRD velocities = 0
ave/max small velocity = 4.20347 7.72323
ave/max big velocity = 2.64047 6.60911
Memory usage per processor = 7.0395 Mbytes
Step Temp 2[8] TotEng pebig ebig Press 2[1] 2[2] 2[3] 2[4] 2[5] 2[6] 2[7] 2[8] 2[9] 2[10] 2[11] 2[12]
0 6.9744766 0 0.30175295 21.011259 26.207244 2.7219331 0 0 0 0 0 0 0 0 0 0 0 0
1000 1.906187 3319 0.25279446 20.53509 21.955199 9.3471734 14560 77 77 9 6051 40320 5625 3319 1.2073941 0 1062 0
2000 1.6009081 3420 0.25060769 20.572601 21.765278 12.835365 13521 63 63 1 7515 40320 5625 3420 1.1621445 0 1062 0
3000 1.8288965 3342 0.25195284 20.519576 21.882104 7.1720562 13550 62 62 1 8561 40320 5625 3342 1.1411375 0 1062 0
4000 1.8940522 3383 0.25185767 20.462769 21.873838 7.3011842 13210 55 55 0 9536 40320 5625 3383 1.1541274 0 1062 0
5000 1.7848964 3377 0.25063675 20.438054 21.767801 5.7235852 13258 60 60 1 10628 40320 5625 3377 1.1401366 0 1062 0
6000 2.0436442 3393 0.2528503 20.437534 21.960048 4.3266887 13369 65 65 1 11693 40320 5625 3393 1.1363335 0 1062 0
7000 2.0248598 3377 0.25304247 20.468218 21.976738 9.1535003 13126 69 69 2 12819 40320 5625 3377 1.1556137 0 1062 0
8000 1.9627172 3397 0.25201977 20.425693 21.887917 7.9563872 13216 58 58 3 14113 40320 5625 3397 1.144738 0 1062 0
9000 2.4805409 3385 0.25688024 20.462046 22.310049 6.0178529 13031 42 42 1 15329 40320 5625 3385 1.1667813 0 1062 0
10000 2.534712 3377 0.25802283 20.520922 22.409283 5.8494568 12804 49 50 1 16835 40320 5625 3377 1.1559636 0 1062 0
11000 2.4295377 3368 0.25671703 20.485869 22.295874 7.0797587 13097 56 56 2 18180 40320 5625 3368 1.1634145 0 1062 0
12000 2.5699389 3377 0.25784825 20.479516 22.39412 10.440503 13263 63 63 5 19848 40320 5625 3377 1.1714081 0 1062 0
13000 2.5416625 3359 0.25893582 20.595038 22.488576 6.823793 12926 57 57 1 21307 40320 5625 3359 1.1478982 0 1062 0
14000 2.9400702 3353 0.26140738 20.512878 22.703231 10.592709 12848 67 67 0 22782 40320 5625 3353 1.1642018 0 1062 0
15000 2.7315042 3364 0.26050572 20.589951 22.624922 7.6633608 12640 49 49 0 24273 40320 5625 3364 1.1532394 0 1062 0
16000 2.7176527 3389 0.25840703 20.418 22.442651 8.0150175 13019 59 59 1 25970 40320 5625 3389 1.1555937 0 1062 0
17000 2.883266 3355 0.26040537 20.468173 22.616206 6.053208 12916 54 54 1 27532 40320 5625 3355 1.1363169 0 1062 0
18000 2.8072496 3369 0.26013039 20.500924 22.592325 5.7179969 12989 45 45 0 29011 40320 5625 3369 1.1341799 0 1062 0
19000 2.8007066 3372 0.2602809 20.51887 22.605396 8.9063895 12817 59 59 1 30504 40320 5625 3372 1.1253001 0 1062 0
20000 3.0147455 3352 0.26233709 20.537991 22.783977 10.01786 13099 63 63 2 32100 40320 5625 3352 1.1263738 0 1062 0
21000 2.5502874 3354 0.25762968 20.475174 22.375138 7.4525361 12917 55 55 3 33904 40320 5625 3354 1.1432482 0 1062 0
22000 2.5643661 3357 0.25754499 20.45733 22.367782 9.5567227 12893 58 58 0 35500 40320 5625 3357 1.1309581 0 1062 0
23000 2.8451307 3387 0.26090715 20.540164 22.659786 7.8308779 12735 65 65 2 37080 40320 5625 3387 1.1666359 0 1062 0
24000 2.6241046 3341 0.25916876 20.553849 22.508807 9.7905794 12815 54 54 1 38864 40320 5625 3341 1.1422596 0 1062 0
25000 2.6054468 3374 0.25897446 20.550874 22.491932 9.4099539 13094 59 59 4 40505 40320 5625 3374 1.1887482 0 1062 0
26000 3.1370949 3356 0.26301193 20.50545 22.842586 9.930938 12785 59 59 0 42126 40320 5625 3356 1.1539158 0 1062 0
27000 2.757373 3361 0.26029623 20.552484 22.606727 7.4336931 12687 59 59 2 43728 40320 5625 3361 1.1408006 0 1062 0
28000 2.9939441 3393 0.26250984 20.568491 22.79898 9.201853 12802 66 66 0 45174 40320 5625 3393 1.1429688 0 1062 0
29000 3.1611971 3376 0.26252235 20.444974 22.800066 8.2146139 13042 48 48 0 46696 40320 5625 3376 1.1370708 0 1062 0
30000 2.2624796 3374 0.25475303 20.439753 22.125301 9.3304907 13249 53 53 1 48164 40320 5625 3374 1.1062325 0 1062 0
31000 2.6864602 3393 0.2607928 20.648442 22.649855 7.6108593 13028 58 58 0 49556 40320 5625 3393 1.1154125 0 1062 0
32000 2.8404087 3377 0.26037369 20.497351 22.613455 6.9017135 13055 63 63 1 51004 40320 5625 3377 1.1320087 0 1062 0
33000 2.8467378 3365 0.26095982 20.543541 22.664361 7.194747 13306 58 58 0 52481 40320 5625 3365 1.1312548 0 1062 0
34000 3.0833296 3382 0.26259198 20.509033 22.806114 7.44833 12988 59 60 1 54049 40320 5625 3382 1.1447344 0 1062 0
35000 3.1068536 3382 0.26441995 20.650266 22.964872 9.3804156 13139 54 54 0 55749 40320 5625 3382 1.1272861 0 1062 0
36000 2.6998442 3369 0.25864735 20.452138 22.463522 6.3327985 13317 52 52 0 57348 40320 5625 3369 1.1312386 0 1062 0
37000 2.600171 3376 0.25882039 20.541424 22.478551 9.3909788 12841 57 57 0 58880 40320 5625 3376 1.1319429 0 1062 0
38000 2.3962821 3367 0.25728672 20.560122 22.345352 7.4364329 13039 66 66 0 60627 40320 5625 3367 1.1468141 0 1062 0
39000 2.8548804 3422 0.26044489 20.492753 22.619638 6.0832552 13249 50 50 0 62327 40320 5625 3422 1.1442726 0 1062 0
40000 2.7721124 3375 0.25927374 20.452701 22.517924 10.368128 13197 55 55 1 63903 40320 5625 3375 1.140807 0 1062 0
41000 3.0463159 3378 0.26123216 20.418508 22.688013 10.179554 13030 64 64 0 65629 40320 5625 3378 1.1277894 0 1062 0
42000 2.9625288 3388 0.26205988 20.552817 22.759901 8.8326646 12958 64 64 0 67242 40320 5625 3388 1.159201 0 1062 0
43000 2.7524867 3347 0.2592403 20.464417 22.51502 10.30221 12815 65 65 2 68734 40320 5625 3347 1.1444123 0 1062 0
44000 2.8779981 3354 0.26009923 20.44551 22.589618 12.127509 12743 67 67 3 70309 40320 5625 3354 1.1467915 0 1062 0
45000 2.7833999 3365 0.25968205 20.479753 22.553386 8.6942365 12909 46 47 0 71846 40320 5625 3365 1.1480348 0 1062 0
46000 2.6564678 3370 0.25935897 20.546258 22.525327 9.2618375 12923 61 61 1 73491 40320 5625 3370 1.1528284 0 1062 0
47000 2.821625 3414 0.25999655 20.47859 22.5807 7.5359036 12861 55 55 2 75201 40320 5625 3414 1.1395745 0 1062 0
48000 2.9301623 3351 0.26139433 20.519127 22.702098 9.5177842 13037 49 49 2 76819 40320 5625 3351 1.1472505 0 1062 0
49000 2.528824 3390 0.25834586 20.553364 22.437338 8.953785 13127 49 49 0 78450 40320 5625 3390 1.1495903 0 1062 0
50000 3.0054471 3375 0.261046 20.432787 22.671845 7.9749959 13089 56 56 2 80717 40320 5625 3375 1.1574724 0 1062 0
Loop time of 10.5055 on 8 procs for 50000 steps with 17370 atoms
Pair time (%) = 0.156617 (1.49081)
Neigh time (%) = 0.556925 (5.30126)
Comm time (%) = 0.902042 (8.58637)
Outpt time (%) = 0.00305247 (0.0290559)
Other time (%) = 8.88688 (84.5925)
Nlocal: 2171.25 ave 2567 max 1849 min
Histogram: 2 0 0 1 1 2 1 0 0 1
Nghost: 62.625 ave 71 max 57 min
Histogram: 2 0 1 0 2 2 0 0 0 1
Neighs: 168.25 ave 232 max 120 min
Histogram: 1 0 1 2 2 1 0 0 0 1
Total # of neighbors = 1346
Ave neighs/atom = 0.0774899
Neighbor list builds = 2599
Dangerous builds = 17
Please see the log.cite file for references relevant to this simulation

190
examples/ASPHERE/ellipsoid/log.ellipsoid.1Feb14 Normal file
View File

@ -0,0 +1,190 @@
LAMMPS (1 Feb 2014)
# SRD diffusion demo - ellipsoids
units lj
atom_style ellipsoid
atom_modify first big
dimension 2
# create big ellipsoidal particles
lattice sq 0.14
Lattice spacing in x,y,z = 2.67261 2.67261 2.67261
region box block 0 10 0 10 -0.5 0.5
create_box 2 box
Created orthogonal box = (0 0 -1.33631) to (26.7261 26.7261 1.33631)
4 by 2 by 1 MPI processor grid
create_atoms 1 region box
Created 100 atoms
set type 1 mass 1.0
100 settings made for mass
set type 1 shape 3.0 1.0 1.0
100 settings made for shape
group big type 1
100 atoms in group big
set group big quat/random 29898
100 settings made for quat/random
velocity big create 1.44 87287 loop geom
# equilibrate big particles
pair_style gayberne 1.0 3.0 1.0 4.0
pair_coeff 1 1 1.0 1.0 1 1 1 1 1 1
pair_coeff 1 2 1.0 1.0 1 1 1 1 1 1 0.0
pair_coeff 2 2 1.0 1.0 1 1 1 1 1 1 0.0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
fix 1 big nve/asphere
fix 2 all enforce2d
compute rot big temp/asphere
#dump 1 all custom 10 dump.ellipsoid.equil id type x y z # quatw quati quatj quatk
thermo_style custom step temp c_rot epair etotal press
thermo 100
run 1000
Memory usage per processor = 3.31932 Mbytes
Step Temp rot E_pair TotEng Press
0 1.44 0.95677852 2.7038078 4.1294078 10.518912
100 2.5524145 2.801098 -0.37027046 2.1566199 0.84703874
200 2.6266386 2.7938164 -0.35322565 2.2471465 1.004886
300 2.9987557 2.9499545 -0.58917376 2.3795944 0.73081788
400 2.8557446 2.8208128 -0.39904801 2.4281391 0.91808964
500 2.4399047 2.8255746 -0.40056447 2.0149412 1.0538908
600 2.854258 2.9166789 -0.53424483 2.2914706 0.8117508
700 2.9593679 2.8231211 -0.40051714 2.5292571 1.1630889
800 2.7632971 2.9060854 -0.52075339 2.2149107 0.77106814
900 2.9905601 2.8869672 -0.49099457 2.4696599 0.69616725
1000 2.8470146 2.9004954 -0.51281252 2.305732 0.68820531
Loop time of 0.0485955 on 8 procs for 1000 steps with 100 atoms
Pair time (%) = 0.0201517 (41.4682)
Neigh time (%) = 0.000492364 (1.01319)
Comm time (%) = 0.0191883 (39.4858)
Outpt time (%) = 0.000273198 (0.562188)
Other time (%) = 0.00848994 (17.4706)
Nlocal: 12.5 ave 16 max 8 min
Histogram: 1 1 0 0 0 3 0 0 1 2
Nghost: 34.125 ave 42 max 28 min
Histogram: 1 0 2 2 0 1 0 1 0 1
Neighs: 46.125 ave 65 max 20 min
Histogram: 1 0 1 0 1 2 0 0 1 2
Total # of neighbors = 369
Ave neighs/atom = 3.69
Neighbor list builds = 174
Dangerous builds = 0
#undump 1
unfix 1
unfix 2
# add small particles as hi density lattice
region plane block INF INF INF INF -0.001 0.001 units box
lattice sq 120.0
Lattice spacing in x,y,z = 0.0912871 0.0912871 0.0912871
create_atoms 2 region plane
Created 85849 atoms
set type 2 mass 0.01
85849 settings made for mass
group small type 2
85849 atoms in group small
velocity small create 1.0 593849 loop geom
# delete overlaps
# must set 1-2 cutoff to non-zero value
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 2.0
neigh_modify one 10000
delete_atoms overlap 1.6 small big
Deleted 71060 atoms, new total = 14889
# SRD run
reset_timestep 0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
communicate multi group big vel yes
neigh_modify include big
# no pairwise interactions with small particles
pair_style gayberne 1.0 3.0 1.0 4.0
pair_coeff 1 1 1.0 1.0 1 1 1 1 1 1
pair_coeff 1 2 1.0 1.0 1 1 1 1 1 1 0.0
pair_coeff 2 2 1.0 1.0 1 1 1 1 1 1 0.0
# use fix SRD to push small particles out from inside big ones
# if comment out, big particles won't see SRD particles
timestep 0.0005
fix 1 big nve/asphere
fix 2 small srd 20 big 1.0 0.25 49894 shift yes 54979 collision noslip search 0.2 inside ignore exact no bounce 50
fix 3 all enforce2d
# diagnostics
compute tbig big temp/asphere
variable pebig equal pe*atoms/count(big)
variable ebig equal etotal*atoms/count(big)
thermo_style custom step temp c_rot f_2[9] etotal v_pebig v_ebig press f_2[1] f_2[2] f_2[3] f_2[4] f_2[5] f_2[6] f_2[7] f_2[8] f_2[9] f_2[10] f_2[11] f_2[12]
thermo_modify temp tbig
WARNING: Temperature for thermo pressure is not for group all (../thermo.cpp:439)
thermo 1000
#dump 1 all custom 1000 dump.ellipsoid id type x y z # quatw quati quatj quatk
#dump 1 all image 1000 image.*.jpg type type zoom 1.6
#dump_modify 1 pad 6 adiam 1 1 adiam 2 0.2
run 100000
SRD info:
SRD/big particles = 14789 100
big particle diameter max/min = 3 1
SRD temperature & lamda = 1 0.1
SRD max distance & max velocity = 0.4 40
SRD grid counts: 107 107 1
SRD grid size: request, actual (xyz) = 0.25, 0.249777 0.249777 2.67261
SRD per actual grid cell = 1.92757
SRD viscosity = 0.688101
big/SRD mass density ratio = 1.37367
# of rescaled SRD velocities = 0
ave/max small velocity = 13.3093 24.3359
ave/max big velocity = 2.08028 5.05671
Memory usage per processor = 10.9077 Mbytes
Step Temp rot 2[9] TotEng pebig ebig Press 2[1] 2[2] 2[3] 2[4] 2[5] 2[6] 2[7] 2[8] 2[9] 2[10] 2[11] 2[12]
0 2.9004954 2.9004954 0 0.025582146 -0.51281252 3.8089257 0.89865242 0 0 0 0 0 0 0 0 0 0 0 0
1000 3.3106106 3.3106106 1.0635866 0.029956176 -0.47263474 4.4601751 1.8111247 14000 35 160 0 224946 98040 11449 4182 1.0635866 1047 50 0
2000 2.2814487 2.2814487 1.2237286 0.017972589 -0.72341981 2.6759388 1.3454977 15318 47 276 0 528340 98040 11449 4136 1.2237286 2383 50 0
3000 4.2476485 4.2476485 1.2992612 0.038526442 -0.59279435 5.7362019 1.6896403 16506 63 410 0 830594 98040 11449 4172 1.2992612 3834 50 0
4000 2.4798681 2.4798681 1.3453264 0.022837826 -0.29467953 3.4003239 2.0765724 17201 51 403 0 1185573 98040 11449 4233 1.3453264 5456 50 0
5000 3.0852629 3.0852629 1.3915904 0.02760985 -0.48621111 4.1108305 1.2676357 17540 61 591 0 1594856 98040 11449 4306 1.3915904 7419 50 0
6000 6.7201122 6.7201122 1.3649223 0.064627842 -0.39052776 9.6224394 2.5524159 17997 59 564 0 2010250 98040 11449 4213 1.3649223 9280 50 0
7000 2.720343 2.720343 1.3909677 0.024282471 -0.43789405 3.6154171 1.5288376 17922 64 342 0 2392388 98040 11449 4228 1.3909677 10947 50 0
8000 1.9181626 1.9181626 1.3514483 0.014305579 -0.72810462 2.1299577 0.89423917 18188 42 260 0 2751626 98040 11449 4236 1.3514483 12556 50 0
9000 2.9090186 2.9090186 1.3688847 0.02558926 -0.52445278 3.8099849 1.1756318 18536 63 440 0 3112508 98040 11449 4227 1.3688847 14164 50 0
10000 1.9647403 1.9647403 1.3979066 0.015375037 -0.63827389 2.2891892 0.65597214 18817 60 466 0 3501352 98040 11449 4257 1.3979066 16016 50 0
11000 1.9989051 1.9989051 1.3872043 0.016191243 -0.56765443 2.4107142 0.80582429 18756 58 767 0 3856028 98040 11449 4254 1.3872043 17611 50 0
12000 1.7538755 1.7538755 1.4594883 0.013930543 -0.53915598 2.0741185 1.0563875 18602 58 335 0 4248666 98040 11449 4251 1.4594883 19350 50 0
13000 13.730871 13.730871 1.4864545 0.13447802 -0.43656501 20.022433 3.6701166 19144 74 681 0 4708523 98040 11449 4104 1.4864545 21589 50 0
14000 5.1078782 5.1078782 1.3637807 0.049721182 -0.20775172 7.4029868 2.559935 18494 65 713 0 5139143 98040 11449 4167 1.3637807 23555 50 0
15000 2.5843962 2.5843962 1.3450247 0.022294373 -0.53134114 3.3194092 1.1069147 18489 52 260 0 5489252 98040 11449 4262 1.3450247 25209 50 0
ERROR: Lost atoms: original 14889 current 14819 (../thermo.cpp:392)

247
examples/ASPHERE/ellipsoid/log.ellipsoid.mp.1Feb14 Normal file
View File

@ -0,0 +1,247 @@
LAMMPS (1 Feb 2014)
# SRD viscosity demo - ellipsoids
units lj
atom_style ellipsoid
atom_modify first big
dimension 2
# create big ellipsoidal particles
lattice sq 0.14
Lattice spacing in x,y,z = 2.67261 2.67261 2.67261
region box block 0 10 0 10 -0.5 0.5
create_box 2 box
Created orthogonal box = (0 0 -1.33631) to (26.7261 26.7261 1.33631)
4 by 2 by 1 MPI processor grid
create_atoms 1 region box
Created 100 atoms
set type 1 mass 1.0
100 settings made for mass
set type 1 shape 3.0 1.0 1.0
100 settings made for shape
group big type 1
100 atoms in group big
set group big quat/random 29898
100 settings made for quat/random
velocity big create 1.44 87287 loop geom
# equilibrate big particles
pair_style gayberne 1.0 3.0 1.0 4.0
pair_coeff 1 1 1.0 1.0 1 1 1 1 1 1
pair_coeff 1 2 1.0 1.0 1 1 1 1 1 1 0.0
pair_coeff 2 2 1.0 1.0 1 1 1 1 1 1 0.0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
fix 1 big nve/asphere
fix 2 all enforce2d
compute rot big temp/asphere
#dump 1 all custom 10 dump.ellipsoid.equil id type x y z # quatw quati quatj quatk
thermo_style custom step temp c_rot epair etotal press
thermo 100
run 1000
Memory usage per processor = 3.31932 Mbytes
Step Temp rot E_pair TotEng Press
0 1.44 0.95677852 2.7038078 4.1294078 10.518912
100 2.5524145 2.801098 -0.37027046 2.1566199 0.84703874
200 2.6266386 2.7938164 -0.35322565 2.2471465 1.004886
300 2.9987557 2.9499545 -0.58917376 2.3795944 0.73081788
400 2.8557446 2.8208128 -0.39904801 2.4281391 0.91808964
500 2.4399047 2.8255746 -0.40056447 2.0149412 1.0538908
600 2.854258 2.9166789 -0.53424483 2.2914706 0.8117508
700 2.9593679 2.8231211 -0.40051714 2.5292571 1.1630889
800 2.7632971 2.9060854 -0.52075339 2.2149107 0.77106814
900 2.9905601 2.8869672 -0.49099457 2.4696599 0.69616725
1000 2.8470146 2.9004954 -0.51281252 2.305732 0.68820531
Loop time of 0.0581853 on 8 procs for 1000 steps with 100 atoms
Pair time (%) = 0.0202803 (34.8546)
Neigh time (%) = 0.000480682 (0.826122)
Comm time (%) = 0.0264942 (45.5342)
Outpt time (%) = 0.000326395 (0.560958)
Other time (%) = 0.0106037 (18.2241)
Nlocal: 12.5 ave 16 max 8 min
Histogram: 1 1 0 0 0 3 0 0 1 2
Nghost: 34.125 ave 42 max 28 min
Histogram: 1 0 2 2 0 1 0 1 0 1
Neighs: 46.125 ave 65 max 20 min
Histogram: 1 0 1 0 1 2 0 0 1 2
Total # of neighbors = 369
Ave neighs/atom = 3.69
Neighbor list builds = 174
Dangerous builds = 0
#undump 1
unfix 1
unfix 2
# add small particles as hi density lattice
region plane block INF INF INF INF -0.001 0.001 units box
lattice sq 120.0
Lattice spacing in x,y,z = 0.0912871 0.0912871 0.0912871
create_atoms 2 region plane
Created 85849 atoms
set type 2 mass 0.01
85849 settings made for mass
group small type 2
85849 atoms in group small
velocity small create 1.0 593849 loop geom
# delete overlaps
# must set 1-2 cutoff to non-zero value
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 2.0
neigh_modify one 10000
delete_atoms overlap 1.6 small big
Deleted 71060 atoms, new total = 14889
# SRD run
reset_timestep 0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
communicate multi group big vel yes
neigh_modify include big
# no pairwise interactions with small particles
pair_style gayberne 1.0 3.0 1.0 4.0
pair_coeff 1 1 1.0 1.0 1 1 1 1 1 1
pair_coeff 1 2 1.0 1.0 1 1 1 1 1 1 0.0
pair_coeff 2 2 1.0 1.0 1 1 1 1 1 1 0.0
# use fix SRD to push small particles out from inside big ones
# if comment out, big particles won't see SRD particles
timestep 0.0005
fix 1 big nve/asphere
fix 2 small srd 20 big 1.0 0.25 49894 shift yes 54979 collision noslip search 0.2 inside ignore exact no bounce 50
fix 3 small viscosity 20 x y 50
fix 4 all enforce2d
# diagnostics
compute tbig big temp/asphere
variable pebig equal pe*atoms/count(big)
variable ebig equal etotal*atoms/count(big)
thermo_style custom step temp c_rot f_2[9] etotal v_pebig v_ebig press f_2[1] f_2[2] f_2[3] f_2[4] f_2[5] f_2[6] f_2[7] f_2[8] f_2[9] f_2[10] f_2[11] f_2[12]
thermo_modify temp tbig
WARNING: Temperature for thermo pressure is not for group all (../thermo.cpp:439)
thermo 1000
#dump 1 all custom 500 dump.ellipsoid.mp id type x y z # quatw quati quatj quatk
#dump 1 all image 500 image.*.jpg type type zoom 1.6
#dump_modify 1 pad 6 adiam 1 1 adiam 2 0.2
run 50000
SRD info:
SRD/big particles = 14789 100
big particle diameter max/min = 3 1
SRD temperature & lamda = 1 0.1
SRD max distance & max velocity = 0.4 40
SRD grid counts: 107 107 1
SRD grid size: request, actual (xyz) = 0.25, 0.249777 0.249777 2.67261
SRD per actual grid cell = 1.92757
SRD viscosity = 0.688101
big/SRD mass density ratio = 1.37367
# of rescaled SRD velocities = 0
ave/max small velocity = 13.3093 24.3359
ave/max big velocity = 2.08028 5.05671
Memory usage per processor = 10.9077 Mbytes
Step Temp rot 2[9] TotEng pebig ebig Press 2[1] 2[2] 2[3] 2[4] 2[5] 2[6] 2[7] 2[8] 2[9] 2[10] 2[11] 2[12]
0 2.9004954 2.9004954 0 0.025582146 -0.51281252 3.8089257 0.89865242 0 0 0 0 0 0 0 0 0 0 0 0
1000 3.1663129 3.1663129 1.0612788 0.028825811 -0.42593131 4.291875 1.8917479 14347 45 319 0 247730 98040 11449 4128 1.0612788 1160 50 0
2000 2.3854865 2.3854865 1.2711805 0.019550917 -0.6434388 2.910936 1.4342405 15199 41 151 0 526468 98040 11449 4191 1.2711805 2321 50 0
3000 2.6796789 2.6796789 1.2494843 0.023573743 -0.48282696 3.5098946 1.6090734 16441 52 230 0 827117 98040 11449 4203 1.2494843 3659 50 0
4000 2.7945661 2.7945661 1.3021617 0.024256529 -0.5523488 3.6115547 1.2702377 16935 52 474 0 1127934 98040 11449 4215 1.3021617 4966 50 0
5000 2.7061959 2.7061959 1.338483 0.023736966 -0.498035 3.5341969 1.2801942 17578 49 277 0 1435089 98040 11449 4261 1.338483 6351 50 0
6000 3.2698703 3.2698703 1.3285507 0.029840303 -0.42918398 4.4429227 1.4528039 17882 63 228 0 1776449 98040 11449 4252 1.3285507 7935 50 0
7000 3.9340641 3.9340641 1.3823898 0.035368008 -0.59581277 5.2659427 1.1321242 18069 49 462 0 2090062 98040 11449 4265 1.3823898 9343 50 0
8000 3.6037728 3.6037728 1.4150583 0.032400676 -0.54548483 4.8241367 1.2091317 18275 56 369 0 2486575 98040 11449 4229 1.4150583 11121 50 0
9000 105.15275 105.15275 1.3630934 1.0504642 -0.27397996 156.40361 23.278069 18139 107 732 0 2861341 98040 11449 4217 1.3630934 12781 50 0
10000 3.5224307 3.5224307 1.3897056 0.0324883 -0.4112387 4.837183 1.6093705 19100 61 347 0 3248400 98040 11449 4204 1.3897056 14424 50 0
11000 4.1949959 4.1949959 1.421272 0.03853305 -0.51335814 5.7371858 1.3586332 18670 52 590 0 3722041 98040 11449 4263 1.421272 16373 50 0
12000 3.1692742 3.1692742 1.3999982 0.028463927 -0.48422456 4.237994 1.1713046 18787 61 323 0 4117517 98040 11449 4202 1.3999982 18125 50 0
13000 3.5868757 3.5868757 1.3193021 0.032298734 -0.53548629 4.8089585 1.0664078 18572 67 224 0 4417636 98040 11449 4277 1.3193021 19546 50 0
14000 3.708681 3.708681 1.317797 0.033355183 -0.55968141 4.9662533 1.3959088 18284 67 207 0 4675359 98040 11449 4253 1.317797 20631 50 0
15000 3.7410393 3.7410393 1.4017441 0.033377852 -0.60452027 4.9696283 1.2506767 18385 51 189 0 4956150 98040 11449 4275 1.4017441 21763 50 0
16000 5.3298522 5.3298522 1.3456075 0.049673767 -0.54555257 7.3959272 1.841432 18693 63 549 0 5301393 98040 11449 4196 1.3456075 23216 50 0
17000 3.7425275 3.7425275 1.321554 0.036608704 -0.12569607 5.4506699 2.3767421 18377 68 526 0 5637038 98040 11449 4261 1.321554 24829 50 0
18000 6.1904009 6.1904009 1.3905006 0.05857936 -0.50181651 8.7218809 1.7936734 18422 63 522 0 6049781 98040 11449 4280 1.3905006 26646 50 0
19000 3.1770374 3.1770374 1.3694174 0.028096405 -0.55051202 4.1832737 0.98036773 18342 54 343 0 6383253 98040 11449 4239 1.3694174 28038 50 0
20000 3.5793112 3.5793112 1.3160294 0.031413769 -0.65597761 4.677196 1.0713128 18739 60 272 0 6692090 98040 11449 4262 1.3160294 29389 50 0
21000 3.2300481 3.2300481 1.3407923 0.027955254 -0.65051393 4.1622578 0.89783071 18582 55 456 0 6999361 98040 11449 4237 1.3407923 30653 50 0
22000 3.5879886 3.5879886 1.3790201 0.031692391 -0.6274228 4.7186802 0.97752405 18768 63 552 0 7350833 98040 11449 4308 1.3790201 32120 50 0
23000 3.5666316 3.5666316 1.4123826 0.031275136 -0.6577261 4.656555 0.99950405 18864 49 231 0 7741276 98040 11449 4265 1.4123826 33939 50 0
24000 3.7306353 3.7306353 1.3854831 0.032426295 -0.7306956 4.827951 0.93575015 18488 63 552 0 8096740 98040 11449 4238 1.3854831 35503 50 0
25000 4.2971901 4.2971901 1.4207289 0.038358698 -0.69158669 5.7112265 1.0308127 18922 58 430 0 8518804 98040 11449 4254 1.4207289 37429 50 0
26000 9.514164 9.514164 1.4209186 0.09200552 -0.47740255 13.698702 2.4102159 18934 63 570 0 8931705 98040 11449 4233 1.4209186 39132 50 0
27000 3.9153535 3.9153535 1.3814961 0.036650525 -0.37698007 5.4568967 1.4531475 19153 65 553 0 9277922 98040 11449 4299 1.3814961 40660 50 0
28000 3.6560048 3.6560048 1.3520603 0.033100525 -0.51910991 4.9283372 1.3474705 19529 59 447 0 9642607 98040 11449 4278 1.3520603 42335 50 0
29000 14.835118 14.835118 1.3981622 0.14535835 -0.4619207 21.642405 3.6917556 19265 52 102 0 9937354 98040 11449 4224 1.3981622 43664 50 0
30000 5.3595251 5.3595251 1.3636194 0.050893845 -0.40810785 7.5775845 1.7023989 19065 43 131 0 10261222 98040 11449 4170 1.3636194 45193 50 0
31000 5.1821698 5.1821698 1.3180039 0.049074014 -0.41480304 7.30663 1.6943582 19142 58 282 0 10584027 98040 11449 4228 1.3180039 46636 50 0
32000 4.3056534 4.3056534 1.3508467 0.040670309 -0.36002116 6.0554024 1.7875405 19065 36 284 0 10933981 98040 11449 4286 1.3508467 48240 50 0
33000 3.8277999 3.8277999 1.3770405 0.034411971 -0.57982354 5.1235984 1.1046562 18838 71 729 0 11329901 98040 11449 4257 1.3770405 50046 50 0
34000 4.6562056 4.6562056 1.3938501 0.043777831 -0.41966518 6.5180812 1.3815841 18914 61 269 0 11758952 98040 11449 4221 1.3938501 51773 50 0
35000 3.4251492 3.4251492 1.3813031 0.031876928 -0.35731644 4.7461558 1.4152845 18838 60 283 0 12159425 98040 11449 4238 1.3813031 53477 50 0
36000 4.2565688 4.2565688 1.3474871 0.039214915 -0.5035788 5.8387088 1.1222014 18765 56 240 0 12568711 98040 11449 4249 1.3474871 55170 50 0
37000 5.3745505 5.3745505 1.4215553 0.053606226 -0.02664916 7.981431 2.8011643 19103 58 504 0 12941065 98040 11449 4270 1.4215553 56743 50 0
38000 4.0346121 4.0346121 1.4354472 0.036753623 -0.53932498 5.472247 1.1475792 18829 62 623 0 13351480 98040 11449 4226 1.4354472 58443 50 0
39000 4.0471971 4.0471971 1.417888 0.037884509 -0.3896991 5.6406246 1.3824278 19047 70 579 0 13752101 98040 11449 4219 1.417888 60208 50 0
40000 4.1611575 4.1611575 1.3712048 0.039505792 -0.31810741 5.8820173 1.6615204 18811 65 378 0 14143802 98040 11449 4313 1.3712048 61793 50 0
41000 4.852673 4.852673 1.3719985 0.044831688 -0.5554927 6.67499 1.2744313 18899 53 414 0 14505662 98040 11449 4274 1.3719985 63358 50 0
42000 4.8338875 4.8338875 1.3595624 0.045465888 -0.43307636 6.7694161 1.6297026 18708 46 360 0 14839866 98040 11449 4224 1.3595624 64895 50 0
43000 5.0243115 5.0243115 1.4371431 0.047111836 -0.47174292 7.0144813 1.5899963 18704 55 291 0 15195383 98040 11449 4209 1.4371431 66440 50 0
44000 7.2748224 7.2748224 1.3693754 0.069498669 -0.49182853 10.347657 1.9007074 18751 73 423 0 15572639 98040 11449 4253 1.3693754 68148 50 0
45000 5.4888002 5.4888002 1.3501677 0.050885254 -0.60200685 7.5763055 1.2841374 18662 51 314 0 15972088 98040 11449 4227 1.3501677 69998 50 0
46000 3.2144774 3.2144774 1.3342547 0.028194354 -0.59171403 4.1978573 0.9459541 18948 50 534 0 16308430 98040 11449 4253 1.3342547 71549 50 0
47000 4.0640905 4.0640905 1.3473081 0.036715084 -0.58898597 5.4665088 1.2043249 19188 57 236 0 16670321 98040 11449 4239 1.3473081 73216 50 0
48000 3.7843815 3.7843815 1.3413467 0.034360184 -0.52284075 5.1158877 1.1195984 19477 57 283 0 17000588 98040 11449 4307 1.3413467 74662 50 0
49000 4.3042567 4.3042567 1.354248 0.039635615 -0.51199566 5.9013468 1.410884 19628 68 456 0 17319627 98040 11449 4277 1.354248 75983 50 0
50000 4.015575 4.015575 1.3058373 0.036258162 -0.58472888 5.3984778 0.84500807 19257 63 416 0 17675269 98040 11449 4254 1.3058373 77468 50 0
Loop time of 11.9502 on 8 procs for 50000 steps with 14889 atoms
Pair time (%) = 1.09132 (9.13222)
Neigh time (%) = 0.959464 (8.02883)
Comm time (%) = 1.19057 (9.96277)
Outpt time (%) = 0.00310364 (0.0259714)
Other time (%) = 8.70577 (72.8502)
Nlocal: 1861.12 ave 1979 max 1606 min
Histogram: 1 0 0 0 1 1 1 0 0 4
Nghost: 35 ave 39 max 32 min
Histogram: 1 0 2 0 3 1 0 0 0 1
Neighs: 46 ave 61 max 31 min
Histogram: 1 1 0 1 1 0 2 0 1 1
Total # of neighbors = 368
Ave neighs/atom = 0.0247162
Neighbor list builds = 2624
Dangerous builds = 0
Please see the log.cite file for references relevant to this simulation

213
examples/ASPHERE/line/log.line.1Feb14 Normal file
View File

@ -0,0 +1,213 @@
LAMMPS (1 Feb 2014)
# Aspherical shear demo - 2d line box and triangle mixture, implicit solvent
units lj
atom_style line
dimension 2
read_data data.line
orthogonal box = (-22.5539 -22.5539 -0.5) to (22.5539 22.5539 0.5)
4 by 2 by 1 MPI processor grid
reading atoms ...
350 atoms
350 lines
velocity all create 1.44 320984 loop geom
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
neigh_modify exclude molecule all
pair_style line/lj 2.5
pair_coeff * * 1.0 0.25
fix 2 all rigid molecule langevin 2.0 2.0 1.0 492983
100 rigid bodies with 350 atoms
fix 3 all deform 1 x scale 0.3 y scale 0.3
fix 4 all enforce2d
compute 10 all property/atom end1x end1y end2x end2y
#dump 1 all custom 500 dump1.atom id type x y z ix iy iz
#dump 2 all custom 500 dump1.line id type # c_10[1] c_10[2] c_10[3] c_10[4]
timestep 0.004
compute 1 all erotate/asphere
compute 2 all ke
compute 3 all pe
variable toteng equal (c_1+c_2+c_3)/atoms
thermo 1000
thermo_style custom step temp f_2 pe ke c_1 c_2 c_3 v_toteng
run 10000
Memory usage per processor = 2.6072 Mbytes
Step Temp 2 PotEng KinEng 1 2 3 toteng
0 0 1.2780105 0 -0 0.037823677 0.50989511 0 0.54771879
1000 0 1.9896906 -0.13333756 -0 0.12630626 0.72641827 -0.13333756 0.71938697
2000 0 2.0408541 -0.24906647 -0 0.13199238 0.74265938 -0.24906647 0.62558529
3000 0 1.9921323 -0.39759798 -0 0.11671173 0.73705927 -0.39759798 0.45617302
4000 0 2.1392159 -0.36475197 -0 0.11923802 0.7975688 -0.36475197 0.55205485
5000 0 2.139715 -0.52582091 -0 0.15846417 0.75855653 -0.52582091 0.39119979
6000 0 2.1313904 -0.65532027 -0 0.11090422 0.80254883 -0.65532027 0.25813278
7000 0 1.9416614 -0.81322598 -0 0.11722471 0.71491587 -0.81322598 0.018914608
8000 0 1.9388183 -1.0581149 -0 0.10142762 0.72949452 -1.0581149 -0.22719275
9000 0 2.2830265 -1.583347 -0 0.14583927 0.83260066 -1.583347 -0.60490709
10000 0 4.1416666 -2.6034045 -0 0.18839177 1.5866082 -2.6034045 -0.82840455
Loop time of 2.84713 on 8 procs for 10000 steps with 350 atoms
Pair time (%) = 1.26607 (44.4682)
Neigh time (%) = 0.0163046 (0.572669)
Comm time (%) = 0.786148 (27.612)
Outpt time (%) = 0.000368953 (0.0129588)
Other time (%) = 0.778241 (27.3342)
Nlocal: 43.75 ave 51 max 39 min
Histogram: 1 2 0 2 0 1 1 0 0 1
Nghost: 170.25 ave 180 max 160 min
Histogram: 2 1 0 0 0 0 2 0 1 2
Neighs: 963.125 ave 1209 max 767 min
Histogram: 2 0 2 0 1 0 1 0 1 1
Total # of neighbors = 7705
Ave neighs/atom = 22.0143
Neighbor list builds = 987
Dangerous builds = 0
#undump 1
#undump 2
unfix 3
change_box all triclinic
triclinic box = (-6.76616 -6.76616 -0.5) to (6.76616 6.76616 0.5) with tilt (0 0 0)
#dump 1 all custom 500 dump2.atom id type x y z ix iy iz
#dump 2 all custom 500 dump2.line id type # c_10[1] c_10[2] c_10[3] c_10[4]
fix 3 all deform 1 xy erate 0.01 units box
run 100000
Memory usage per processor = 2.75978 Mbytes
Step Temp 2 PotEng KinEng 1 2 3 toteng
10000 0 4.1416666 -2.6024092 -0 0.18839177 1.5866082 -2.6024092 -0.82740923
11000 0 1.8408319 -3.1031477 -0 0.12073234 0.66819561 -3.1031477 -2.3142198
12000 0 2.0793172 -3.0329681 -0 0.086927592 0.80420833 -3.0329681 -2.1418322
13000 0 2.2022136 -2.99776 -0 0.14309291 0.8007129 -2.99776 -2.0539542
14000 0 1.9510757 -3.094649 -0 0.09482969 0.74134559 -3.094649 -2.2584737
15000 0 1.9874689 -3.1431753 -0 0.1083061 0.7434663 -3.1431753 -2.2914029
16000 0 1.8484778 -2.9491537 -0 0.079102883 0.71310191 -2.9491537 -2.1569489
17000 0 2.1978438 -2.9675694 -0 0.11677634 0.82515673 -2.9675694 -2.0256363
18000 0 2.0293397 -2.9860257 -0 0.1287845 0.7409325 -2.9860257 -2.1163087
19000 0 2.0077219 -3.005622 -0 0.12697603 0.7334762 -3.005622 -2.1451698
20000 0 2.1806369 -3.0622132 -0 0.11066657 0.82389212 -3.0622132 -2.1276545
21000 0 1.8156509 -3.1031481 -0 0.10227614 0.67585994 -3.1031481 -2.325012
22000 0 2.1028516 -3.0861182 -0 0.098877162 0.80234493 -3.0861182 -2.1848961
23000 0 1.8994891 -3.0110243 -0 0.10961187 0.70445488 -3.0110243 -2.1969576
24000 0 1.9305389 -3.0057136 -0 0.11735151 0.7100223 -3.0057136 -2.1783398
25000 0 1.9553918 -3.0848948 -0 0.13217467 0.70585039 -3.0848948 -2.2468697
26000 0 1.8903754 -2.9543658 -0 0.099925113 0.71023579 -2.9543658 -2.1442049
27000 0 2.2624684 -3.2416154 -0 0.11398815 0.85564117 -3.2416154 -2.2719861
28000 0 2.0335234 -3.1795174 -0 0.10291986 0.76859015 -3.1795174 -2.3080074
29000 0 1.7056403 -3.1198739 -0 0.076174496 0.65481419 -3.1198739 -2.3888853
30000 0 2.1203465 -3.0863113 -0 0.11355683 0.79516311 -3.0863113 -2.1775914
31000 0 1.8446708 -3.2764357 -0 0.11006455 0.68050865 -3.2764357 -2.4858625
32000 0 1.7947968 -3.081031 -0 0.087540776 0.68165784 -3.081031 -2.3118324
33000 0 2.2576228 -3.1042222 -0 0.15666855 0.81088407 -3.1042222 -2.1366696
34000 0 1.8522306 -3.0037311 -0 0.08276626 0.71104684 -3.0037311 -2.209918
35000 0 2.0611686 -2.9877406 -0 0.10822003 0.77513794 -2.9877406 -2.1043826
36000 0 2.0739798 -2.981184 -0 0.11784198 0.77100651 -2.981184 -2.0923355
37000 0 2.2120215 -2.8961216 -0 0.13172943 0.8162798 -2.8961216 -1.9481124
38000 0 2.1097357 -2.9746777 -0 0.11881736 0.78535507 -2.9746777 -2.0705053
39000 0 1.8928127 -2.9482169 -0 0.10915752 0.70204792 -2.9482169 -2.1370115
40000 0 1.8515483 -3.003524 -0 0.096981256 0.69653943 -3.003524 -2.2100033
41000 0 1.8569065 -3.0834675 -0 0.10562739 0.6901897 -3.0834675 -2.2876504
42000 0 2.2638459 -3.0577224 -0 0.10896442 0.86125524 -3.0577224 -2.0875027
43000 0 2.0992151 -3.1219185 -0 0.13280434 0.76685929 -3.1219185 -2.2222548
44000 0 2.3479986 -3.0702445 -0 0.14438131 0.86190379 -3.0702445 -2.0639594
45000 0 2.0295235 -2.9619688 -0 0.13726365 0.73253213 -2.9619688 -2.092173
46000 0 2.0597998 -3.1282569 -0 0.10139093 0.7813804 -3.1282569 -2.2454856
47000 0 1.8173541 -3.1122647 -0 0.10380346 0.67506259 -3.1122647 -2.3333986
48000 0 1.9650208 -3.22153 -0 0.10115952 0.74099226 -3.22153 -2.3793782
49000 0 1.8836303 -3.0781944 -0 0.10561306 0.70165705 -3.0781944 -2.2709243
50000 0 1.7799289 -3.1090208 -0 0.10522105 0.65760561 -3.1090208 -2.3461941
51000 0 1.7270244 -2.894789 -0 0.079316066 0.66083727 -2.894789 -2.1546357
52000 0 2.1036512 -3.0708266 -0 0.11762021 0.7839446 -3.0708266 -2.1692618
53000 0 2.1749106 -3.043193 -0 0.11908953 0.81301499 -3.043193 -2.1110885
54000 0 1.7245018 -3.1528646 -0 0.11118993 0.62788226 -3.1528646 -2.4137924
55000 0 1.7599209 -3.0543405 -0 0.098996756 0.65525507 -3.0543405 -2.3000886
56000 0 1.9085953 -3.1072383 -0 0.10931083 0.70865859 -3.1072383 -2.2892689
57000 0 1.9682028 -3.131335 -0 0.094973379 0.74854212 -3.131335 -2.2878195
58000 0 1.7813545 -3.167275 -0 0.085523136 0.6779145 -3.167275 -2.4038374
59000 0 2.0324497 -3.1103534 -0 0.13343285 0.73761703 -3.1103534 -2.2393035
60000 0 1.697349 -3.152831 -0 0.086035815 0.64139945 -3.152831 -2.4253957
61000 0 2.061537 -3.0730484 -0 0.12753143 0.75598441 -3.0730484 -2.1895325
62000 0 1.8186587 -3.1590894 -0 0.096701457 0.68272371 -3.1590894 -2.3796643
63000 0 1.9692063 -3.1053354 -0 0.094864849 0.7490807 -3.1053354 -2.2613899
64000 0 2.0113722 -3.0898117 -0 0.12640828 0.73560838 -3.0898117 -2.227795
65000 0 1.9350085 -3.2290712 -0 0.11850297 0.71078638 -3.2290712 -2.3997819
66000 0 2.0569306 -3.0489252 -0 0.12731012 0.75423156 -3.0489252 -2.1673835
67000 0 2.0119071 -3.10894 -0 0.099142151 0.76310375 -3.10894 -2.2466941
68000 0 2.0615321 -3.167552 -0 0.11066625 0.77284751 -3.167552 -2.2840383
69000 0 1.8731077 -3.1864825 -0 0.12982742 0.672933 -3.1864825 -2.3837221
70000 0 2.109088 -3.1787326 -0 0.13107619 0.77281866 -3.1787326 -2.2748377
71000 0 1.9764202 -3.0731116 -0 0.12035585 0.72668137 -3.0731116 -2.2260744
72000 0 1.6460591 -3.0865952 -0 0.09486718 0.61058673 -3.0865952 -2.3811413
73000 0 1.9756864 -3.0431747 -0 0.1041775 0.74254526 -3.0431747 -2.196452
74000 0 2.0553658 -3.1290715 -0 0.12267262 0.75819842 -3.1290715 -2.2482004
75000 0 1.550997 -2.9813403 -0 0.094350599 0.57036238 -2.9813403 -2.3166273
76000 0 1.9824851 -3.0057088 -0 0.11621875 0.73341774 -3.0057088 -2.1560723
77000 0 1.8451505 -3.0063772 -0 0.12602143 0.66475735 -3.0063772 -2.2155984
78000 0 1.934755 -3.0603773 -0 0.11800375 0.71117696 -3.0603773 -2.2311966
79000 0 2.1982093 -3.095224 -0 0.12840585 0.81368383 -3.095224 -2.1531343
80000 0 1.7547828 -3.0620364 -0 0.092552269 0.65949751 -3.0620364 -2.3099866
81000 0 2.1064588 -3.1631358 -0 0.11415036 0.78861768 -3.1631358 -2.2603677
82000 0 2.0039288 -3.1190532 -0 0.13374009 0.72508654 -3.1190532 -2.2602266
83000 0 1.7749465 -3.2172309 -0 0.11017601 0.65051537 -3.2172309 -2.4565396
84000 0 1.7865147 -3.1433093 -0 0.086859958 0.67878922 -3.1433093 -2.3776602
85000 0 1.4541199 -3.0123872 -0 0.096301993 0.52689225 -3.0123872 -2.3891929
86000 0 2.0345567 -2.9577061 -0 0.11756878 0.7543841 -2.9577061 -2.0857532
87000 0 1.9655671 -3.0235755 -0 0.13068174 0.71170417 -3.0235755 -2.1811896
88000 0 2.1238733 -3.0972095 -0 0.15660293 0.75362848 -3.0972095 -2.1869781
89000 0 1.93413 -3.0018414 -0 0.11736651 0.71154634 -3.0018414 -2.1729285
90000 0 2.0106701 -3.0778989 -0 0.1150105 0.74670527 -3.0778989 -2.2161831
91000 0 1.9700224 -3.0072749 -0 0.1200151 0.72428023 -3.0072749 -2.1629796
92000 0 2.3062898 -3.1245149 -0 0.16045862 0.82795128 -3.1245149 -2.136105
93000 0 1.9977984 -3.1072212 -0 0.098381856 0.75781746 -3.1072212 -2.2510219
94000 0 1.9684844 -3.0923519 -0 0.11228362 0.73135255 -3.0923519 -2.2487157
95000 0 1.8913269 -3.1044766 -0 0.088944817 0.72162386 -3.1044766 -2.293908
96000 0 2.139258 -2.9032081 -0 0.12122357 0.7956013 -2.9032081 -1.9863833
97000 0 2.1549839 -3.0073823 -0 0.14177897 0.78178555 -3.0073823 -2.0838178
98000 0 2.001084 -3.0411655 -0 0.16549603 0.69211139 -3.0411655 -2.1835581
99000 0 2.269108 -3.0749812 -0 0.11370514 0.85876972 -3.0749812 -2.1025064
100000 0 2.0270482 -3.0802101 -0 0.098599302 0.77013563 -3.0802101 -2.2114752
101000 0 1.9412796 -3.1543623 -0 0.11953004 0.71244692 -3.1543623 -2.3223853
102000 0 2.0146461 -3.0140006 -0 0.10152606 0.7618937 -3.0140006 -2.1505808
103000 0 1.7377282 -3.1862597 -0 0.082111131 0.66262952 -3.1862597 -2.441519
104000 0 1.7973897 -3.1055088 -0 0.093370304 0.67693958 -3.1055088 -2.3351989
105000 0 2.2615194 -2.9636424 -0 0.12584102 0.84338157 -2.9636424 -1.9944198
106000 0 1.8974182 -2.9505576 -0 0.085928679 0.72725057 -2.9505576 -2.1373783
107000 0 2.0691161 -3.0804349 -0 0.12411017 0.76265389 -3.0804349 -2.1936708
108000 0 2.0457472 -2.981702 -0 0.09308074 0.78366806 -2.981702 -2.1049532
109000 0 1.6610604 -3.1463569 -0 0.066318676 0.64556436 -3.1463569 -2.4344738
110000 0 1.9606721 -3.080013 -0 0.13164067 0.70864736 -3.080013 -2.2397249
Loop time of 62.2251 on 8 procs for 100000 steps with 350 atoms
Pair time (%) = 43.4946 (69.8987)
Neigh time (%) = 0.395421 (0.635469)
Comm time (%) = 10.3551 (16.6414)
Outpt time (%) = 0.00358662 (0.00576394)
Other time (%) = 7.97644 (12.8187)
Nlocal: 43.75 ave 51 max 33 min
Histogram: 1 0 1 0 0 2 0 1 2 1
Nghost: 168.375 ave 180 max 155 min
Histogram: 1 1 0 0 1 0 2 2 0 1
Neighs: 971 ave 1278 max 631 min
Histogram: 1 1 0 1 0 1 2 1 0 1
Total # of neighbors = 7768
Ave neighs/atom = 22.1943
Neighbor list builds = 7621
Dangerous builds = 0

244
examples/ASPHERE/line/log.line.srd.1Feb14 Normal file
View File

@ -0,0 +1,244 @@
LAMMPS (1 Feb 2014)
# Aspherical shear demo - 2d line boxes, solvated by SRD particles
units lj
atom_style line
atom_modify first big
dimension 2
read_data data.line.srd
orthogonal box = (-28.7968 -28.7968 -0.5) to (28.7968 28.7968 0.5)
4 by 2 by 1 MPI processor grid
reading atoms ...
400 atoms
400 lines
# add small particles as hi density lattice
lattice sq 0.4
Lattice spacing in x,y,z = 1.58114 1.58114 1.58114
region plane block INF INF INF INF -0.001 0.001
lattice sq 10.0
Lattice spacing in x,y,z = 0.316228 0.316228 0.316228
create_atoms 2 region plane
Created 33489 atoms
group big type 1
400 atoms in group big
group small type 2
33489 atoms in group small
set group small mass 0.01
33489 settings made for mass
# delete overlaps
# must set 1-2 cutoff to non-zero value
pair_style lj/cut 1.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0
delete_atoms overlap 1.5 small big
Deleted 13605 atoms, new total = 20284
# SRD run
reset_timestep 0
velocity small create 1.44 87287 loop geom
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
neigh_modify exclude molecule big include big
communicate multi group big vel yes
neigh_modify include big
# no pairwise interactions with small particles
pair_style line/lj 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 0.0
# use fix SRD to push small particles out from inside big ones
# if comment out, big particles won't see SRD particles
timestep 0.001
fix 1 big rigid molecule
100 rigid bodies with 400 atoms
fix 2 small srd 20 big 1.0 0.25 49894 search 0.2 cubic warn 0.0001 shift yes 49829 overlap yes collision noslip
fix 3 all deform 1 x scale 0.35 y scale 0.35
fix 4 all enforce2d
# diagnostics
compute tsmall small temp/deform
compute tbig big temp
variable pebig equal pe*atoms/count(big)
variable ebig equal etotal*atoms/count(big)
compute 1 big erotate/asphere
compute 2 all ke
compute 3 all pe
variable toteng equal (c_1+c_2+c_3)/atoms
thermo 1000
thermo_style custom step temp c_tsmall f_2[9] c_1 etotal v_pebig v_ebig press
thermo_modify temp tbig
WARNING: Temperature for thermo pressure is not for group all (../thermo.cpp:439)
compute 10 big property/atom end1x end1y end2x end2y
#dump 1 all custom 500 dump1.atom.srd id type x y z ix iy iz
#dump 2 all custom 500 dump1.line.srd id type # c_10[1] c_10[2] c_10[3] c_10[4]
run 10000
WARNING: Using fix srd with box deformation but no SRD thermostat (../fix_srd.cpp:385)
SRD info:
SRD/big particles = 19884 400
big particle diameter max/min = 1.99882 0.503306
SRD temperature & lamda = 1 0.2
SRD max distance & max velocity = 0.8 40
SRD grid counts: 230 230 1
SRD grid size: request, actual (xyz) = 0.25, 0.250407 0.250407 1
SRD per actual grid cell = 0.444963
SRD viscosity = 4.2356
big/SRD mass density ratio = 14.0918
WARNING: SRD bin size for fix srd differs from user request (../fix_srd.cpp:2853)
WARNING: Fix srd grid size > 1/4 of big particle diameter (../fix_srd.cpp:2875)
# of rescaled SRD velocities = 0
ave/max small velocity = 15.906 29.1054
ave/max big velocity = 0 0
WARNING: Using compute temp/deform with inconsistent fix deform remap option (../compute_temp_deform.cpp:76)
Memory usage per processor = 7.79007 Mbytes
Step Temp tsmall 2[9] 1 TotEng pebig ebig Press
0 0 1.4528554 0 0 0 0 0 0
1000 0 1.1122612 1.1071958 0.00058011072 0 0 0 0.32625408
2000 0 1.0254475 1.0231236 0.00072347646 -2.3144253e-06 -0.00011736451 -0.00011736451 0.44526308
3000 0 1.0089214 0.99966408 0.00061511441 -0.00023298189 -0.011814512 -0.011814512 0.87208025
4000 0 1.0066185 0.99391102 0.00079065587 -0.0021630644 -0.10968899 -0.10968899 1.3901494
5000 0 1.0109003 1.0009124 0.00071588733 -0.0040023493 -0.20295913 -0.20295913 1.2736445
6000 0 1.0119642 1.0105049 0.00080741361 -0.0067908165 -0.34436231 -0.34436231 2.5492857
7000 0 1.0048989 0.98912274 0.00098851463 -0.010933999 -0.55446309 -0.55446309 3.5707613
8000 0 1.0021597 0.99377939 0.0008405671 -0.015690841 -0.79568252 -0.79568252 6.685381
9000 0 1.0197546 1.00114 0.0012645557 -0.021124602 -1.0712286 -1.0712286 18.907993
10000 0 1.0401079 1.0184189 0.0036111452 0.0099395852 0.50403637 0.50403637 96.909257
Loop time of 3.93391 on 8 procs for 10000 steps with 20284 atoms
Pair time (%) = 0.110261 (2.80284)
Neigh time (%) = 0.32853 (8.35124)
Comm time (%) = 0.199613 (5.07416)
Outpt time (%) = 0.00108692 (0.0276295)
Other time (%) = 3.29442 (83.7441)
Nlocal: 2535.5 ave 2850 max 2226 min
Histogram: 1 0 1 1 1 1 1 0 1 1
Nghost: 109.25 ave 117 max 102 min
Histogram: 1 0 2 0 1 1 1 1 0 1
Neighs: 526 ave 665 max 407 min
Histogram: 1 2 0 1 0 1 1 1 0 1
Total # of neighbors = 4208
Ave neighs/atom = 0.207454
Neighbor list builds = 522
Dangerous builds = 0
#undump 1
#undump 2
unfix 3
change_box all triclinic
triclinic box = (-10.0789 -10.0789 -0.5) to (10.0789 10.0789 0.5) with tilt (0 0 0)
fix 2 small srd 20 big 1.0 0.25 49894 search 0.2 cubic warn 0.0001 shift yes 49829 overlap yes collision noslip tstat yes
#dump 1 all custom 500 dump2.atom.srd id type x y z ix iy iz
#dump 2 all custom 500 dump2.line.srd id type # c_10[1] c_10[2] c_10[3] c_10[4]
fix 3 all deform 1 xy erate 0.05 units box remap v
run 40000
SRD info:
SRD/big particles = 19884 400
big particle diameter max/min = 1.99882 0.503306
SRD temperature & lamda = 1 0.2
SRD max distance & max velocity = 0.8 40
SRD grid counts: 81 81 1
SRD grid size: request, actual (xyz) = 0.25, 0.248861 0.248861 1
SRD per actual grid cell = -11.3331
SRD viscosity = -1920.36
big/SRD mass density ratio = -0.546468
WARNING: SRD bin size for fix srd differs from user request (../fix_srd.cpp:2853)
WARNING: Fix srd grid size > 1/4 of big particle diameter (../fix_srd.cpp:2875)
WARNING: Fix srd viscosity < 0.0 due to low SRD density (../fix_srd.cpp:2877)
# of rescaled SRD velocities = 1
ave/max small velocity = 12.724 40
ave/max big velocity = 1.54523 5.36901
Memory usage per processor = 4.23847 Mbytes
Step Temp tsmall 2[9] 1 TotEng pebig ebig Press
10000 0 1.0277413 0 0.0036111452 0.010135973 0.51399517 0.51399517 60.794865
11000 0 1.0024742 1 0.00040482851 -0.031643325 -1.604633 -1.604633 55.531632
12000 0 1.0036177 1 0.00036416993 -0.030844063 -1.5641024 -1.5641024 56.694788
13000 0 1.0043067 1 0.00068862721 -0.030673787 -1.5554677 -1.5554677 56.852537
14000 0 1.0042533 1 0.00050997053 -0.035967153 -1.8238943 -1.8238943 51.254109
15000 0 1.003956 1 0.00048458218 -0.038060147 -1.9300301 -1.9300301 50.346943
16000 0 1.0059929 1 0.0004740426 -0.037433402 -1.8982478 -1.8982478 49.82532
17000 0 1.0042401 1 0.00039855238 -0.034314446 -1.7400855 -1.7400855 48.430648
18000 0 1.0025952 1 0.00046207703 -0.034166102 -1.732563 -1.732563 53.883613
19000 0 1.0025881 1 0.00034818943 -0.037341367 -1.8935807 -1.8935807 48.90662
20000 0 1.0009932 1 0.00031244041 -0.035274728 -1.7887814 -1.7887814 47.665935
21000 0 1.0025603 1 0.0005514826 -0.034350221 -1.7418997 -1.7418997 63.290704
22000 0 1.0038575 1 0.00034453716 -0.03576596 -1.8136918 -1.8136918 49.805726
23000 0 1.0030321 1 0.00049535709 -0.035873083 -1.819124 -1.819124 48.973244
24000 0 1.0048593 1 0.00034225992 -0.03322158 -1.6846663 -1.6846663 52.557383
25000 0 1.0024423 1 0.00059157362 -0.030152319 -1.5290241 -1.5290241 56.733821
26000 0 1.0022703 1 0.00047512976 -0.029563064 -1.499143 -1.499143 56.271943
27000 0 1.0024537 1 0.00054158319 -0.02957039 -1.4995145 -1.4995145 55.246787
28000 0 1.0023325 1 0.00051129428 -0.035115211 -1.7806924 -1.7806924 49.718172
29000 0 1.002865 1 0.00056728135 -0.035484703 -1.7994293 -1.7994293 48.387401
30000 0 1.0022693 1 0.00040284402 -0.033844303 -1.7162446 -1.7162446 50.486676
31000 0 1.0037106 1 0.00056291948 -0.036529699 -1.852421 -1.852421 52.011803
32000 0 1.0022542 1 0.00041877447 -0.036903458 -1.8713743 -1.8713743 52.712289
33000 0 1.0021364 1 0.00040871451 -0.033616728 -1.7047043 -1.7047043 55.199593
34000 0 1.003779 1 0.00040667148 -0.031745412 -1.6098098 -1.6098098 51.055261
35000 0 1.0015547 1 0.00042509365 -0.033603064 -1.7040114 -1.7040114 53.887389
36000 0 1.0024849 1 0.00040589851 -0.032151258 -1.6303903 -1.6303903 55.673615
37000 0 1.0038204 1 0.00040542177 -0.033643242 -1.7060488 -1.7060488 51.834694
38000 0 1.0031777 1 0.00040237126 -0.034711811 -1.7602359 -1.7602359 56.275778
39000 0 1.002777 1 0.00051021239 -0.034375104 -1.7431615 -1.7431615 48.461145
40000 0 1.0053762 1 0.00041566465 -0.036016384 -1.8263908 -1.8263908 48.937456
41000 0 1.0039078 1 0.00049726673 -0.036481888 -1.8499966 -1.8499966 49.290465
42000 0 1.0033408 1 0.00045236191 -0.034254524 -1.7370469 -1.7370469 52.146603
43000 0 1.0039893 1 0.00048362351 -0.032519474 -1.6490625 -1.6490625 53.320599
44000 0 1.0036667 1 0.00046834006 -0.031099875 -1.5770747 -1.5770747 55.907531
45000 0 1.0035765 1 0.0005195031 -0.029559708 -1.4989728 -1.4989728 56.191165
46000 0 1.004514 1 0.00050134348 -0.028852909 -1.463131 -1.463131 57.547696
47000 0 1.003864 1 0.00038963389 -0.031339229 -1.5892123 -1.5892123 55.775764
48000 0 1.003145 1 0.00035436452 -0.032390682 -1.6425315 -1.6425315 55.233354
49000 0 1.0042615 1 0.00053161075 -0.029960468 -1.5192953 -1.5192953 54.611685
50000 0 1.0037096 1 0.00031191389 -0.02942011 -1.4918938 -1.4918938 58.653762
Loop time of 27.1275 on 8 procs for 40000 steps with 20284 atoms
Pair time (%) = 1.54277 (5.68709)
Neigh time (%) = 1.26864 (4.67658)
Comm time (%) = 1.19871 (4.4188)
Outpt time (%) = 0.00421953 (0.0155544)
Other time (%) = 23.1132 (85.202)
Nlocal: 2535.5 ave 2670 max 2406 min
Histogram: 1 1 1 0 0 2 1 1 0 1
Nghost: 114.5 ave 124 max 103 min
Histogram: 1 0 1 0 2 1 0 1 0 2
Neighs: 521.5 ave 692 max 442 min
Histogram: 3 1 1 1 0 0 0 1 0 1
Total # of neighbors = 4172
Ave neighs/atom = 0.205679
Neighbor list builds = 2002
Dangerous builds = 0
Please see the log.cite file for references relevant to this simulation

288
examples/ASPHERE/poly/log.poly.1Feb14 Normal file
View File

@ -0,0 +1,288 @@
LAMMPS (1 Feb 2014)
# SRD diffusion demo - poydisperse spheres
units lj
atom_style sphere
atom_modify first big
dimension 2
# create big particles with 3 different types and diameters
lattice sq 0.3
Lattice spacing in x,y,z = 1.82574 1.82574 1.82574
region box block 0 10 0 10 -0.5 0.5
create_box 4 box
Created orthogonal box = (0 0 -0.912871) to (18.2574 18.2574 0.912871)
2 by 4 by 1 MPI processor grid
create_atoms 1 region box
Created 100 atoms
group big type 1
100 atoms in group big
set group big type/fraction 2 0.33 394895
35 settings made for type/fraction
set group big type/fraction 3 0.5 989894
57 settings made for type/fraction
group big type 2 3
100 atoms in group big
set type 1*3 mass 1.0
100 settings made for mass
velocity big create 1.44 87287 loop geom
# equilibrate big particles, repulsive only to prevent aggregation
pair_style lj/cut 1.12
pair_coeff 1 1 1.0 1.0 1.12
pair_coeff 2 2 1.0 2.0 2.24
pair_coeff 3 3 1.0 1.5 1.68
pair_coeff 4 4 0.0 1.0 0.0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
fix 1 big nve
fix 2 all enforce2d
#dump 1 all atom 10 dump.poly.equil
run 1000
Memory usage per processor = 2.7472 Mbytes
Step Temp E_pair E_mol TotEng Press
0 1.44 -0.16013916 0 1.2654608 1.1298975
1000 1.3367862 -0.30816328 0 1.0152551 1.5440006
Loop time of 0.0221665 on 8 procs for 1000 steps with 100 atoms
Pair time (%) = 0.000323534 (1.45957)
Neigh time (%) = 0.000274181 (1.23692)
Comm time (%) = 0.0146933 (66.2864)
Outpt time (%) = 2.864e-05 (0.129204)
Other time (%) = 0.00684676 (30.8879)
Nlocal: 12.5 ave 18 max 10 min
Histogram: 1 3 1 1 0 1 0 0 0 1
Nghost: 27.75 ave 31 max 24 min
Histogram: 1 0 1 0 2 2 0 0 0 2
Neighs: 19 ave 30 max 16 min
Histogram: 4 1 2 0 0 0 0 0 0 1
Total # of neighbors = 152
Ave neighs/atom = 1.52
Neighbor list builds = 115
Dangerous builds = 0
#undump 1
unfix 1
unfix 2
# add small particles as hi density lattice
region plane block INF INF INF INF -0.001 0.001 units box
lattice sq 250.0
Lattice spacing in x,y,z = 0.0632456 0.0632456 0.0632456
create_atoms 4 region plane
Created 83521 atoms
set type 4 mass 0.1
83521 settings made for mass
group small type 4
83521 atoms in group small
velocity small create 1.0 593849 loop geom
# delete overlaps
# must set *-4 cutoffs to non-zero values
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 1.0 2.0
pair_coeff 3 3 1.0 1.5
pair_coeff 1 4 0.0 1.0 0.5
pair_coeff 2 4 0.0 1.0 1.0
pair_coeff 3 4 0.0 1.0 0.75
delete_atoms overlap 1.0 small big
Deleted 63410 atoms, new total = 20211
# SRD run
reset_timestep 0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
communicate multi group big vel yes
neigh_modify include big
# no pairwise interactions with small particles
pair_style lj/cut 1.12
pair_coeff 1 1 1.0 1.0 1.12
pair_coeff 2 2 1.0 2.0 2.24
pair_coeff 3 3 1.0 1.5 1.68
pair_coeff 4 4 0.0 1.0 0.0
# use fix SRD to push small particles out from inside big ones
# if comment out, big particles won't see SRD particles
timestep 0.001
fix 1 big nve
fix 2 small srd 20 big 1.0 0.25 49894 shift yes 54979 search 0.2 inside ignore
fix 3 all enforce2d
# diagnostics
compute tbig big temp/sphere
variable pebig equal pe*atoms/count(big)
variable ebig equal etotal*atoms/count(big)
thermo_style custom step temp f_2[8] etotal v_pebig v_ebig press f_2[1] f_2[2] f_2[3] f_2[4] f_2[5] f_2[6] f_2[7] f_2[8] f_2[9] f_2[10] f_2[11]
thermo_modify temp tbig
WARNING: Temperature for thermo pressure is not for group all (../thermo.cpp:439)
thermo 1000
#dump 1 all atom 1000 dump.poly
#dump 1 all image 1000 image.*.jpg type type zoom 1.6
#dump_modify 1 pad 6 adiam 1 1 adiam 2 2.0 adiam 3 1.5 adiam 4 0.1
run 100000
SRD info:
SRD/big particles = 20111 100
big particle diameter max/min = 1 1
SRD temperature & lamda = 1 0.0632456
SRD max distance & max velocity = 0.252982 12.6491
SRD grid counts: 73 73 1
SRD grid size: request, actual (xyz) = 0.25, 0.250102 0.250102 1.82574
SRD per actual grid cell = 4.93717
SRD viscosity = 0.23322
big/SRD mass density ratio = 0.161311
WARNING: Fix srd grid size > 1/4 of big particle diameter (../fix_srd.cpp:2875)
# of rescaled SRD velocities = 0
ave/max small velocity = 4.19344 7.74495
ave/max big velocity = 1.43991 3.5724
Memory usage per processor = 11.8633 Mbytes
Step Temp 2[8] TotEng pebig ebig Press 2[1] 2[2] 2[3] 2[4] 2[5] 2[6] 2[7] 2[8] 2[9] 2[10] 2[11]
0 0.88820023 0 0.0050232797 -0.30816328 1.0152551 1.5440006 0 0 0 0 0 0 0 0 0 0 0
1000 0.84607633 3934 0.0065245151 0.058016014 1.3186697 4.2649214 7470 26 26 0 665 36480 5329 3934 0.95810956 0 2
2000 0.86586785 3944 0.0045898954 -0.36247933 0.92766376 2.1563866 7786 30 30 0 1374 36480 5329 3944 1.0799777 0 49
3000 0.78968677 4001 0.00417319 -0.33318986 0.84344343 1.9302042 7957 34 34 0 2129 36480 5329 4001 1.1247747 0 49
4000 0.7110667 4005 0.0031642256 -0.41996775 0.63952164 1.9679689 8188 24 24 0 2925 36480 5329 4005 1.1345941 0 49
5000 0.80379425 4052 0.0038912991 -0.41118297 0.78647047 1.6851034 8326 30 30 0 3794 36480 5329 4052 1.1310774 0 49
6000 0.96927432 4048 0.0054656969 -0.33954674 1.104672 1.8976158 8380 23 23 0 4643 36480 5329 4048 1.1861264 0 49
7000 0.75821666 4050 0.0036201309 -0.39807817 0.73166465 1.6087451 8450 51 51 0 5699 36480 5329 4050 1.1664119 0 49
8000 0.8712338 4035 0.0047243103 -0.343308 0.95483036 1.4707918 8458 42 42 0 6974 36480 5329 4035 1.1947807 0 91
9000 0.87995969 4048 0.0042778416 -0.44654537 0.86459457 1.5567027 8379 36 36 0 7965 36480 5329 4048 1.1675384 0 91
10000 0.83009098 4053 0.0042178709 -0.38436167 0.85247389 1.9119485 8373 34 34 0 8952 36480 5329 4053 1.1933986 0 91
11000 0.73595909 4015 0.0037439098 -0.33989742 0.75668162 1.6495307 8333 22 22 0 9983 36480 5329 4015 1.1959519 0 91
12000 0.84043443 4034 0.0041823741 -0.40694768 0.84529962 1.5115578 8566 31 31 0 11023 36480 5329 4034 1.2064988 0 91
13000 0.64253331 4044 0.0034868675 -0.25264383 0.7047308 1.739124 8528 31 31 0 12101 36480 5329 4044 1.1912336 0 91
14000 0.99330171 4042 0.0055047339 -0.36745778 1.1125618 1.8215191 8548 36 36 0 13098 36480 5329 4042 1.1819124 0 91
15000 0.90379859 4049 0.0049047446 -0.35536197 0.99129793 1.7159536 8459 37 37 0 14130 36480 5329 4049 1.1806762 0 91
16000 0.87446498 4034 0.0053318245 -0.22533777 1.0776151 2.0239395 8690 33 33 0 15176 36480 5329 4034 1.1912858 0 91
17000 0.95733995 4035 0.0053759665 -0.33989994 1.0865366 1.6243881 8518 27 27 0 16147 36480 5329 4035 1.2074487 0 91
18000 1.03543 4032 0.0052946469 -0.47268961 1.0701011 1.5444562 8464 28 28 0 17230 36480 5329 4032 1.1959755 0 91
19000 0.75997582 4051 0.0034371936 -0.43767278 0.6946912 1.4209137 8436 39 39 0 18276 36480 5329 4051 1.182461 0 91
20000 0.89615234 4037 0.0048910458 -0.34673773 0.98852926 1.5843879 8433 28 28 0 19275 36480 5329 4037 1.1850673 0 91
21000 0.85849578 4063 0.0044978314 -0.37010201 0.9090567 1.5445809 8504 40 40 0 20278 36480 5329 4063 1.1931031 0 91
22000 0.79050471 4066 0.0038260812 -0.40456274 0.77328927 1.6214277 8463 44 44 0 21287 36480 5329 4066 1.1948824 0 91
23000 0.89563532 4026 0.0051030186 -0.30312555 1.0313711 1.4848223 8548 35 35 0 22422 36480 5329 4026 1.1975081 0 91
24000 0.83528717 4036 0.0043136275 -0.37275064 0.87182725 1.5504059 8527 34 34 0 23404 36480 5329 4036 1.1755463 0 91
25000 0.64179358 4022 0.0035185741 -0.24513342 0.71113901 1.7238889 8452 47 47 0 24522 36480 5329 4022 1.2066339 0 91
26000 0.86710098 4028 0.0046099591 -0.36026164 0.93171883 1.4255297 8734 33 33 0 25524 36480 5329 4028 1.2001233 0 91
27000 0.8327003 4054 0.0045323379 -0.32469262 0.91603082 1.2654649 8468 34 34 0 26520 36480 5329 4054 1.1814901 0 91
28000 0.84396759 4049 0.0042805934 -0.39236096 0.86515074 1.2929473 8536 37 37 0 27657 36480 5329 4049 1.1897647 0 91
29000 0.77639299 4052 0.0038143413 -0.38590904 0.77091652 1.3627099 8637 36 36 0 28773 36480 5329 4052 1.1858017 0 91
30000 0.69994805 4046 0.0034430868 -0.34704032 0.69588228 1.5351247 8704 46 46 0 29818 36480 5329 4046 1.1724555 0 91
31000 0.73907898 4039 0.003034759 -0.48787254 0.61335514 1.2291035 8486 29 29 0 30848 36480 5329 4039 1.2063852 0 91
32000 0.84649298 4052 0.0044326201 -0.36539769 0.89587685 1.4591662 8513 28 28 0 31917 36480 5329 4052 1.1917883 0 91
33000 0.7188891 4031 0.0033440675 -0.39527528 0.67586949 1.2495703 8800 37 37 0 32938 36480 5329 4031 1.191087 0 91
34000 0.93588921 4054 0.0054112843 -0.30080024 1.0936747 1.3650172 8585 33 33 0 33913 36480 5329 4054 1.1840415 0 91
35000 0.86446961 4039 0.0053251497 -0.21179371 1.076266 1.3791947 8607 22 22 0 34865 36480 5329 4039 1.1916342 0 91
36000 0.90510006 4011 0.0054168442 -0.2538007 1.0947984 1.5747653 8479 31 31 0 35842 36480 5329 4011 1.1760305 0 91
37000 0.95332175 4058 0.0055020495 -0.30843018 1.1120192 1.4326924 8512 29 29 0 36751 36480 5329 4058 1.180738 0 91
38000 0.71413841 4043 0.0038380378 -0.2883604 0.77570582 1.4283837 8611 39 39 0 37689 36480 5329 4043 1.1693806 0 91
39000 0.81161335 4064 0.0044323475 -0.31348213 0.89582176 1.4378305 8697 28 28 0 38696 36480 5329 4064 1.1987209 0 91
40000 0.8294437 4062 0.0051262098 -0.19981285 1.0360583 1.6010501 8559 31 31 0 39743 36480 5329 4062 1.2092371 0 91
41000 0.79325147 4079 0.0045273868 -0.26691453 0.91503015 1.4080428 8775 53 53 0 40806 36480 5329 4079 1.1989619 0 91
42000 0.95693891 4027 0.005623893 -0.28919396 1.136645 1.5334192 8509 34 34 0 41793 36480 5329 4027 1.1931035 0 91
43000 0.66458435 4029 0.0035902162 -0.26461209 0.72561859 1.4253918 8559 44 44 0 42768 36480 5329 4029 1.2111853 0 91
44000 0.84442459 4022 0.0050980388 -0.22782802 1.0303646 1.4587515 8522 32 32 0 43812 36480 5329 4022 1.1943268 0 91
45000 0.96505326 4051 0.0057148128 -0.28290855 1.1550208 1.3014361 8496 32 32 0 44803 36480 5329 4051 1.1912109 0 91
46000 0.75345782 4064 0.0040245838 -0.30924351 0.81340864 1.7179797 8568 29 29 0 45795 36480 5329 4064 1.1975591 0 91
47000 0.87187806 4017 0.004406241 -0.40855293 0.89054538 1.3105257 8586 36 36 0 46822 36480 5329 4017 1.1834816 0 91
48000 0.83729235 4078 0.0045115864 -0.33572888 0.91183673 1.3856733 8513 38 38 0 47866 36480 5329 4078 1.1996374 0 91
49000 0.65112065 4048 0.0030445536 -0.35483503 0.61533473 1.2659451 8705 32 32 0 48907 36480 5329 4048 1.1907782 0 91
50000 0.82730473 4066 0.0045492584 -0.31323342 0.91945063 1.1803982 8562 41 41 0 49965 36480 5329 4066 1.1696702 0 91
51000 0.74252279 4048 0.0046741422 -0.16166808 0.94469088 1.5539445 8564 35 35 0 50948 36480 5329 4048 1.1748481 0 91
52000 0.85671386 4051 0.0058575179 -0.092640701 1.183863 1.6307333 8609 19 19 0 51894 36480 5329 4051 1.1770671 0 91
53000 0.7233951 4068 0.0043421658 -0.20026358 0.87759512 1.453924 8582 28 28 0 52893 36480 5329 4068 1.2027667 0 91
54000 0.72117217 4051 0.0040247209 -0.26111018 0.81343635 1.3496518 8679 25 25 0 53831 36480 5329 4051 1.1823115 0 91
55000 0.87940919 4044 0.0047057065 -0.35924936 0.95107033 1.240057 8559 34 34 0 54813 36480 5329 4044 1.1815144 0 91
56000 0.83192698 4044 0.0046298732 -0.30382753 0.93574367 1.547229 8441 37 37 0 55740 36480 5329 4044 1.1761232 0 91
57000 0.91043729 4058 0.0050029377 -0.34540782 1.0111437 1.4951603 8550 32 32 0 56678 36480 5329 4058 1.2024504 0 91
58000 0.97648852 4036 0.0053292287 -0.37787747 1.0770904 1.4058278 8554 31 31 0 57683 36480 5329 4036 1.1746239 0 91
59000 0.72239185 4018 0.0046615138 -0.1342253 0.94213856 1.8337414 8716 35 35 0 58640 36480 5329 4018 1.1697661 0 91
60000 0.73800253 4032 0.0037590001 -0.33989225 0.75973152 1.399426 8561 32 32 0 59593 36480 5329 4032 1.1958194 0 91
61000 0.72778603 4060 0.0042599892 -0.22341476 0.86098642 1.3778202 8541 31 31 0 60573 36480 5329 4060 1.1980875 0 91
62000 0.80350669 4044 0.0042316822 -0.34195968 0.85526529 1.4379455 8487 37 37 0 61608 36480 5329 4044 1.1942365 0 91
63000 0.90084181 4066 0.0043040273 -0.47236734 0.86988695 1.4672845 8504 32 32 0 62670 36480 5329 4066 1.1829834 0 91
64000 0.87220486 4059 0.0047840275 -0.33268544 0.9668998 1.5398216 8588 34 34 0 63605 36480 5329 4059 1.1884975 0 91
65000 0.82469157 4057 0.004405933 -0.33830733 0.89048312 1.5284565 8783 32 32 0 64693 36480 5329 4057 1.1952722 0 91
66000 0.88692288 4081 0.0047264561 -0.36625105 0.95526405 1.3151847 8745 36 36 0 65755 36480 5329 4081 1.1834684 0 91
67000 0.83322374 4031 0.0047732147 -0.27678896 0.96471442 1.5085732 8415 34 34 0 66796 36480 5329 4031 1.1969 0 91
68000 0.75375865 4020 0.0044379413 -0.22614806 0.89695233 1.6351001 8573 30 30 0 67876 36480 5329 4020 1.1911308 0 91
69000 0.91169697 4055 0.0050371842 -0.34036318 1.0180653 1.2495539 8625 20 20 0 68979 36480 5329 4055 1.1860768 0 91
70000 0.80502435 4075 0.0042565586 -0.33919322 0.86029306 1.3194042 8571 36 36 0 69947 36480 5329 4075 1.18994 0 91
71000 0.8928873 4055 0.0050586026 -0.3080079 1.0223942 1.3949364 8698 32 32 0 71037 36480 5329 4055 1.2154821 0 91
72000 0.95397384 4050 0.0052866702 -0.35293209 1.0684889 1.2950893 8613 35 35 0 72090 36480 5329 4050 1.1790804 0 91
73000 0.78742879 4053 0.004378612 -0.28830762 0.88496127 1.5168042 8624 30 30 0 73174 36480 5329 4053 1.2008591 0 91
74000 0.79498726 4046 0.0042999917 -0.31545969 0.86907133 1.3751647 8562 22 22 0 74250 36480 5329 4046 1.1921334 0 91
75000 0.88929145 4051 0.0048995187 -0.33480254 0.99024172 1.6589606 8393 33 33 0 75296 36480 5329 4051 1.1753137 0 91
76000 0.98366685 4045 0.0049217494 -0.47092883 0.99473478 1.5271346 8555 29 29 0 76301 36480 5329 4045 1.2031327 0 91
77000 0.72618119 4051 0.0035028638 -0.37404618 0.7079638 1.5178344 8565 37 37 0 77315 36480 5329 4051 1.1993308 0 91
78000 0.80445682 4028 0.003992176 -0.39178196 0.8068587 1.3682816 8624 32 32 0 78299 36480 5329 4028 1.2065456 0 91
79000 1.091703 4052 0.0063793569 -0.33730561 1.2893318 1.5717928 8677 39 39 0 79249 36480 5329 4052 1.1725194 0 91
80000 0.77212865 4052 0.0041275152 -0.31625959 0.8342121 1.5583982 8526 33 33 0 80284 36480 5329 4052 1.1981751 0 91
81000 0.80732682 4061 0.004618782 -0.26941493 0.93350203 1.5517327 8444 36 36 0 81331 36480 5329 4061 1.2095831 0 91
82000 0.80166933 4057 0.0042732734 -0.33081602 0.86367129 1.3937505 8476 29 29 0 82335 36480 5329 4057 1.2047807 0 91
83000 0.76980598 4049 0.0041331379 -0.31166241 0.83534851 1.1411979 8561 34 34 0 83375 36480 5329 4049 1.1810941 0 91
84000 0.85529188 4076 0.0046148117 -0.34168531 0.93269959 1.3284464 8664 27 27 0 84510 36480 5329 4076 1.192734 0 91
85000 0.96167157 4034 0.0053261317 -0.35642617 1.0764645 1.3195071 8557 36 36 0 85517 36480 5329 4034 1.2000123 0 91
86000 0.88130862 4047 0.0046646382 -0.37037981 0.94277003 1.3595079 8634 36 36 0 86578 36480 5329 4047 1.1654584 0 91
87000 0.68799517 4037 0.0042184867 -0.17251446 0.85259834 1.5087731 8363 40 40 0 87606 36480 5329 4037 1.1833599 0 91
88000 0.83191621 4041 0.0048307277 -0.26321678 0.97633838 1.4008912 8597 39 39 0 88490 36480 5329 4041 1.1607031 0 91
89000 0.71456519 4036 0.0039378285 -0.26882761 0.79587453 1.2251569 8397 36 36 0 89436 36480 5329 4036 1.1991362 0 91
90000 0.75198251 4048 0.0043377839 -0.24374445 0.87670949 1.4896261 8551 29 29 0 90453 36480 5329 4048 1.1899998 0 91
91000 0.83310145 4032 0.0049902932 -0.232733 1.0085882 1.4140378 8583 48 48 0 91563 36480 5329 4032 1.1953029 0 91
92000 0.8819465 4050 0.0058087956 -0.14008461 1.1740157 1.753218 8610 46 46 0 92618 36480 5329 4050 1.1972268 0 91
93000 0.76721608 4048 0.0041923848 -0.29582907 0.8473229 1.4126521 8671 32 32 0 93690 36480 5329 4048 1.2027935 0 91
94000 0.8942606 4050 0.0051973822 -0.28200539 1.0504429 1.5743434 8751 21 21 0 94750 36480 5329 4050 1.1804594 0 91
95000 0.73609582 4036 0.0035091666 -0.38754512 0.70923765 1.28494 8660 31 31 0 95713 36480 5329 4036 1.1883125 0 91
96000 0.79430429 4078 0.0043400156 -0.30635284 0.87716056 1.4351927 8600 38 38 0 96695 36480 5329 4078 1.1704833 0 91
97000 0.74817621 4044 0.0038163094 -0.34346826 0.77131429 1.255104 8377 36 36 0 97683 36480 5329 4044 1.1756023 0 91
98000 0.78160255 4048 0.0039385343 -0.36857064 0.79601716 1.5320263 8350 25 25 0 98631 36480 5329 4048 1.181196 0 91
99000 0.804448 4053 0.004503827 -0.28835905 0.91026848 1.7569548 8506 35 35 0 99687 36480 5329 4053 1.1950898 0 91
100000 0.74975993 4032 0.004263269 -0.255493 0.86164929 1.5983872 8621 33 33 0 100701 36480 5329 4032 1.1910217 0 91
Loop time of 13.4746 on 8 procs for 100000 steps with 20211 atoms
Pair time (%) = 0.0563478 (0.418178)
Neigh time (%) = 0.794382 (5.89541)
Comm time (%) = 1.23773 (9.18569)
Outpt time (%) = 0.00554895 (0.0411809)
Other time (%) = 11.3806 (84.4595)
Nlocal: 2526.38 ave 2725 max 2337 min
Histogram: 2 0 0 1 1 1 1 0 1 1
Nghost: 23.875 ave 27 max 21 min
Histogram: 1 1 0 2 0 2 0 0 0 2
Neighs: 17 ave 25 max 12 min
Histogram: 3 1 0 1 1 0 0 0 0 2
Total # of neighbors = 136
Ave neighs/atom = 0.00672901
Neighbor list builds = 5008
Dangerous builds = 0
Please see the log.cite file for references relevant to this simulation

239
examples/ASPHERE/poly/log.poly.mp.1Feb14 Normal file
View File

@ -0,0 +1,239 @@
LAMMPS (1 Feb 2014)
# SRD viscosity demo - poydisperse spheres
units lj
atom_style sphere
atom_modify first big
dimension 2
# create big particles with 3 different types and diameters
lattice sq 0.3
Lattice spacing in x,y,z = 1.82574 1.82574 1.82574
region box block 0 10 0 10 -0.5 0.5
create_box 4 box
Created orthogonal box = (0 0 -0.912871) to (18.2574 18.2574 0.912871)
2 by 4 by 1 MPI processor grid
create_atoms 1 region box
Created 100 atoms
group big type 1
100 atoms in group big
set group big type/fraction 2 0.33 394895
35 settings made for type/fraction
set group big type/fraction 3 0.5 989894
57 settings made for type/fraction
group big type 2 3
100 atoms in group big
set type 1*3 mass 1.0
100 settings made for mass
velocity big create 1.44 87287 loop geom
# equilibrate big particles, repulsive only to prevent aggregation
pair_style lj/cut 1.12
pair_coeff 1 1 1.0 1.0 1.12
pair_coeff 2 2 1.0 2.0 2.24
pair_coeff 3 3 1.0 1.5 1.68
pair_coeff 4 4 0.0 1.0 0.0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
fix 1 big nve
fix 2 all enforce2d
#dump 1 all atom 10 dump.poly.equil
run 1000
Memory usage per processor = 2.7472 Mbytes
Step Temp E_pair E_mol TotEng Press
0 1.44 -0.16013916 0 1.2654608 1.1298975
1000 1.3367862 -0.30816328 0 1.0152551 1.5440006
Loop time of 0.0173425 on 8 procs for 1000 steps with 100 atoms
Pair time (%) = 0.00033614 (1.93825)
Neigh time (%) = 0.000266045 (1.53407)
Comm time (%) = 0.0115004 (66.3132)
Outpt time (%) = 2.78056e-05 (0.160332)
Other time (%) = 0.00521213 (30.0541)
Nlocal: 12.5 ave 18 max 10 min
Histogram: 1 3 1 1 0 1 0 0 0 1
Nghost: 27.75 ave 31 max 24 min
Histogram: 1 0 1 0 2 2 0 0 0 2
Neighs: 19 ave 30 max 16 min
Histogram: 4 1 2 0 0 0 0 0 0 1
Total # of neighbors = 152
Ave neighs/atom = 1.52
Neighbor list builds = 115
Dangerous builds = 0
#undump 1
unfix 1
unfix 2
# add small particles as hi density lattice
region plane block INF INF INF INF -0.001 0.001 units box
lattice sq 250.0
Lattice spacing in x,y,z = 0.0632456 0.0632456 0.0632456
create_atoms 4 region plane
Created 83521 atoms
set type 4 mass 0.1
83521 settings made for mass
group small type 4
83521 atoms in group small
velocity small create 1.0 593849 loop geom
# delete overlaps
# must set *-4 cutoffs to non-zero values
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 1.0 2.0
pair_coeff 3 3 1.0 1.5
pair_coeff 1 4 0.0 1.0 0.5
pair_coeff 2 4 0.0 1.0 1.0
pair_coeff 3 4 0.0 1.0 0.75
delete_atoms overlap 1.0 small big
Deleted 63410 atoms, new total = 20211
# SRD run
reset_timestep 0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
communicate multi group big vel yes
neigh_modify include big
# no pairwise interactions with small particles
pair_style lj/cut 1.12
pair_coeff 1 1 1.0 1.0 1.12
pair_coeff 2 2 1.0 2.0 2.24
pair_coeff 3 3 1.0 1.5 1.68
pair_coeff 4 4 0.0 1.0 0.0
# use fix SRD to push small particles out from inside big ones
# if comment out, big particles won't see SRD particles
timestep 0.001
fix 1 big nve
fix 2 small srd 20 big 1.0 0.25 49894 shift yes 54979 search 0.2 inside ignore
fix 3 small viscosity 10 x y 50
fix 4 all enforce2d
# diagnostics
compute tbig big temp/sphere
variable pebig equal pe*atoms/count(big)
variable ebig equal etotal*atoms/count(big)
thermo_style custom step temp f_2[8] etotal v_pebig v_ebig press f_2[1] f_2[2] f_2[3] f_2[4] f_2[5] f_2[6] f_2[7] f_2[8] f_2[9] f_2[10] f_2[11]
thermo_modify temp tbig
WARNING: Temperature for thermo pressure is not for group all (../thermo.cpp:439)
thermo 1000
#dump 1 all atom 500 dump.poly.mp
#dump 1 all image 500 image.*.jpg type type zoom 1.6
#dump_modify 1 pad 6 adiam 1 1 adiam 2 2.0 adiam 3 1.5 adiam 4 0.1
run 50000
SRD info:
SRD/big particles = 20111 100
big particle diameter max/min = 1 1
SRD temperature & lamda = 1 0.0632456
SRD max distance & max velocity = 0.252982 12.6491
SRD grid counts: 73 73 1
SRD grid size: request, actual (xyz) = 0.25, 0.250102 0.250102 1.82574
SRD per actual grid cell = 4.93717
SRD viscosity = 0.23322
big/SRD mass density ratio = 0.161311
WARNING: Fix srd grid size > 1/4 of big particle diameter (../fix_srd.cpp:2875)
# of rescaled SRD velocities = 0
ave/max small velocity = 4.19344 7.74495
ave/max big velocity = 1.43991 3.5724
Memory usage per processor = 11.8633 Mbytes
Step Temp 2[8] TotEng pebig ebig Press 2[1] 2[2] 2[3] 2[4] 2[5] 2[6] 2[7] 2[8] 2[9] 2[10] 2[11]
0 0.88820023 0 0.0050232797 -0.30816328 1.0152551 1.5440006 0 0 0 0 0 0 0 0 0 0 0
1000 0.76223471 3964 0.0058318994 0.042955471 1.1786852 3.9611295 7412 29 29 0 611 36480 5329 3964 0.97424113 0 1
2000 0.83394503 3951 0.0050432506 -0.22328672 1.0192914 2.7492771 7829 30 30 0 1225 36480 5329 3951 1.0832295 0 3
3000 0.86856945 3991 0.0048205728 -0.31988251 0.97428597 2.4105317 7961 34 34 0 2000 36480 5329 3991 1.1243219 0 6
4000 0.88106888 4009 0.0066731082 0.035909261 1.3487019 3.025346 8002 33 33 0 2895 36480 5329 4009 1.1359575 0 6
5000 1.0764198 4036 0.0070136315 -0.18634051 1.4175251 2.7617336 8117 36 36 0 3808 36480 5329 4036 1.1147743 0 6
6000 1.0433139 4004 0.0059532709 -0.35132214 1.2032156 2.3597383 8004 23 23 0 4845 36480 5329 4004 1.1192873 0 6
7000 1.1551053 4026 0.0074766715 -0.20999681 1.5111101 2.280806 8269 36 36 0 5873 36480 5329 4026 1.1704244 0 6
8000 1.2222713 4008 0.0078175708 -0.24117508 1.5800092 1.967726 8218 33 33 0 7077 36480 5329 4008 1.1496529 0 6
9000 1.2240884 4012 0.0081583857 -0.17500037 1.6488913 2.3684675 8502 27 27 0 8359 36480 5329 4012 1.1664455 0 6
10000 1.2071315 4035 0.0073838443 -0.30627724 1.4923488 2.124366 8376 33 33 0 9778 36480 5329 4035 1.1942304 0 6
11000 1.1694272 4036 0.0077950192 -0.16699527 1.5754513 2.4203067 8623 40 40 0 11210 36480 5329 4036 1.1802896 0 42
12000 1.0911561 4044 0.0074117515 -0.12783343 1.4979891 2.4859591 8490 30 30 0 12647 36480 5329 4044 1.171701 0 42
13000 1.3535072 4023 0.0091503262 -0.16735325 1.8493724 2.7718999 8320 33 33 0 14082 36480 5329 4023 1.1841757 0 42
14000 1.4191753 4043 0.0095368758 -0.18707328 1.927498 2.8442438 8422 25 25 0 15462 36480 5329 4043 1.1567854 0 42
15000 1.703769 4053 0.010892787 -0.33707457 2.2015412 2.4160876 8123 34 34 0 16985 36480 5329 4053 1.1875794 0 42
16000 1.4270161 4032 0.010142553 -0.076342616 2.0499114 2.9568071 8360 34 34 0 18487 36480 5329 4032 1.164586 0 42
17000 1.2267528 4013 0.0080642511 -0.19799584 1.6298658 2.5084784 8477 27 27 0 20101 36480 5329 4013 1.1529604 0 42
18000 1.4228586 4050 0.01009573 -0.079611446 2.0404479 2.8013833 8519 33 33 0 21633 36480 5329 4050 1.1549291 0 42
19000 1.3041379 4040 0.0082901927 -0.26763462 1.6755308 1.9403854 8483 30 30 0 23194 36480 5329 4040 1.1592747 0 42
20000 1.3837406 4030 0.0092446911 -0.19332904 1.8684445 2.4550531 8425 25 25 0 24761 36480 5329 4030 1.1757648 0 42
21000 1.7882552 4031 0.012076007 -0.22381856 2.4406817 2.4206092 8581 40 40 0 26371 36480 5329 4031 1.1591574 0 42
22000 1.5079803 4015 0.010373742 -0.15025377 2.0966369 2.6321039 8383 33 33 0 27989 36480 5329 4015 1.1495522 0 42
23000 1.3835348 4044 0.0087922321 -0.28446888 1.776998 2.2729548 8441 30 30 0 29575 36480 5329 4044 1.1616377 0 42
24000 1.5854836 4044 0.0099174248 -0.35795983 2.0044107 2.2250073 8462 35 35 0 31298 36480 5329 4044 1.1796851 0 42
25000 1.5314412 4062 0.010880793 -0.082730301 2.1991171 2.6858232 8473 28 28 0 33079 36480 5329 4062 1.1726593 0 54
26000 1.2896316 4070 0.0090198232 -0.098554631 1.8229965 2.3437649 8453 40 40 0 34816 36480 5329 4070 1.1764482 0 54
27000 1.6956611 4035 0.012064024 -0.088275024 2.43826 2.8138661 8380 32 32 0 36453 36480 5329 4035 1.1622504 0 54
28000 1.3560979 4043 0.0088193796 -0.23810107 1.7824848 2.1799536 8611 30 30 0 38202 36480 5329 4043 1.1681248 0 54
29000 1.5837646 4029 0.010790526 -0.17893612 2.1808731 2.3116753 8480 26 26 0 39846 36480 5329 4029 1.1774511 0 54
30000 1.6549448 4035 0.01137773 -0.16631482 2.2995529 2.6236656 8637 30 30 0 41461 36480 5329 4035 1.1538405 0 54
31000 1.445315 4050 0.01013447 -0.1052417 2.0482777 2.4440312 8504 38 38 0 43099 36480 5329 4050 1.1690375 0 54
32000 1.5572315 4037 0.010828477 -0.13173155 2.1885434 2.5911192 8453 30 30 0 44760 36480 5329 4037 1.1612147 0 54
33000 2.018109 4068 0.013467519 -0.28506208 2.7219203 2.5784978 8701 40 40 0 46462 36480 5329 4068 1.1727186 0 54
34000 1.2455134 4066 0.0083640288 -0.16536113 1.6904539 2.3054829 8428 32 32 0 48055 36480 5329 4066 1.1712016 0 54
35000 1.6248596 4012 0.011213824 -0.15461501 2.2664259 2.8697537 8453 31 31 0 49658 36480 5329 4012 1.1693962 0 54
36000 1.5485873 4039 0.01112737 -0.058442367 2.2489527 3.0386676 8529 42 42 0 51330 36480 5329 4039 1.1750818 0 54
37000 1.5561543 4024 0.011313246 -0.032149816 2.2865202 2.7583779 8484 45 45 0 52981 36480 5329 4024 1.1613247 0 54
38000 1.5174762 4043 0.010120957 -0.21549294 2.0455467 2.2489655 8528 37 37 0 54637 36480 5329 4043 1.1692017 0 54
39000 1.5936544 4045 0.010946409 -0.16216642 2.2123787 2.7579117 8460 35 35 0 56452 36480 5329 4045 1.1670575 0 54
40000 1.7884514 4051 0.011913021 -0.2570519 2.4077408 2.494169 8539 37 37 0 58163 36480 5329 4051 1.1744075 0 54
41000 1.4173685 4042 0.0096196931 -0.16764288 1.9442362 1.9949409 8511 42 42 0 59967 36480 5329 4042 1.1747318 0 54
42000 1.4946304 4055 0.0098254549 -0.24117659 1.9858227 1.9958514 8395 39 39 0 61697 36480 5329 4055 1.1745062 0 54
43000 1.6486882 4041 0.011190255 -0.19488288 2.2616625 2.2173122 8473 28 28 0 63447 36480 5329 4041 1.1631975 0 54
44000 1.7473914 4055 0.011743128 -0.23020951 2.3734036 2.3913609 8421 30 30 0 65095 36480 5329 4055 1.1567659 0 54
45000 1.7826681 4027 0.011993913 -0.23208563 2.4240898 2.4582759 8539 35 35 0 66749 36480 5329 4027 1.1513285 0 54
46000 1.3988187 4049 0.0091738733 -0.23010838 1.8541315 2.3866226 8447 38 38 0 68459 36480 5329 4049 1.1728215 0 54
47000 1.8402398 4064 0.01271321 -0.17249031 2.569467 2.7748374 8524 33 33 0 70139 36480 5329 4064 1.1555276 0 54
48000 1.5634287 4042 0.011047794 -0.096639164 2.2328696 2.7173437 8561 33 33 0 71923 36480 5329 4042 1.177462 0 54
49000 1.789207 4058 0.012453079 -0.14902655 2.5168918 2.9310418 8545 38 38 0 73721 36480 5329 4058 1.1581934 0 54
50000 1.3475703 4010 0.0092484369 -0.13867816 1.8692016 2.1619389 8664 32 32 0 75439 36480 5329 4010 1.1589691 0 54
Loop time of 6.84 on 8 procs for 50000 steps with 20211 atoms
Pair time (%) = 0.0295508 (0.432029)
Neigh time (%) = 0.400359 (5.8532)
Comm time (%) = 0.622996 (9.10813)
Outpt time (%) = 0.00241029 (0.0352382)
Other time (%) = 5.78468 (84.5714)
Nlocal: 2526.38 ave 2651 max 2342 min
Histogram: 1 0 0 2 1 0 0 1 1 2
Nghost: 24.5 ave 27 max 23 min
Histogram: 1 0 4 0 0 2 0 0 0 1
Neighs: 17.375 ave 24 max 13 min
Histogram: 2 2 0 1 1 0 0 0 0 2
Total # of neighbors = 139
Ave neighs/atom = 0.00687744
Neighbor list builds = 2502
Dangerous builds = 0
Please see the log.cite file for references relevant to this simulation

290
examples/ASPHERE/star/log.star.1Feb14 Normal file
View File

@ -0,0 +1,290 @@
LAMMPS (1 Feb 2014)
# SRD diffusion demo - rigid star particles
units lj
atom_style sphere
atom_modify map array first big
dimension 2
# read in clusters of rigid bodies
fix molprop all property/atom mol
read_data data.star fix molprop NULL Molecules
orthogonal box = (-13.2934 -13.2934 -0.5) to (13.2934 13.2934 0.5)
4 by 2 by 1 MPI processor grid
reading atoms ...
270 atoms
set type 1 mass 1.0
270 settings made for mass
group big type 1
270 atoms in group big
velocity big create 1.44 87287 loop geom
# equilibrate big particles
pair_style soft 1.12
pair_coeff 1 1 0.0
pair_coeff 2 2 0.0 0.0
pair_coeff 1 2 0.0 0.0
variable prefactor equal ramp(0,60)
fix soft all adapt 1 pair soft a * * v_prefactor
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
fix 1 big rigid molecule
30 rigid bodies with 270 atoms
fix 2 all enforce2d
#dump 1 all atom 10 dump.star.equil
thermo 100
run 1000
Memory usage per processor = 2.64859 Mbytes
Step Temp E_pair E_mol TotEng Press
0 0 0 0 0 0.24721146
100 0 3.291475 0 3.291475 0.69158089
200 0 4.6176595 0 4.6176595 0.71737945
300 0 7.0921814 0 7.0921814 1.2163497
400 0 8.3666709 0 8.3666709 0.50645957
500 0 10.630838 0 10.630838 0.72764487
600 0 12.42157 0 12.42157 1.0130599
700 0 14.349074 0 14.349074 0.74795509
800 0 16.464746 0 16.464746 1.054549
900 0 18.253108 0 18.253108 0.51551753
1000 0 20.699563 0 20.699563 1.8084662
Loop time of 0.040709 on 8 procs for 1000 steps with 270 atoms
Pair time (%) = 0.00249004 (6.11668)
Neigh time (%) = 0.000995547 (2.44552)
Comm time (%) = 0.0116902 (28.7165)
Outpt time (%) = 0.000232756 (0.571755)
Other time (%) = 0.0253005 (62.1496)
Nlocal: 33.75 ave 42 max 22 min
Histogram: 1 0 1 0 0 2 0 1 1 2
Nghost: 23.375 ave 30 max 14 min
Histogram: 1 0 1 0 0 2 1 0 1 2
Neighs: 52.25 ave 69 max 27 min
Histogram: 1 0 1 0 1 0 1 1 1 2
Total # of neighbors = 418
Ave neighs/atom = 1.54815
Neighbor list builds = 176
Dangerous builds = 0
#undump 1
unfix soft
unfix 1
unfix 2
# add small particles as hi density lattice
region plane block INF INF INF INF -0.001 0.001 units box
lattice sq 85.0
Lattice spacing in x,y,z = 0.108465 0.108465 0.108465
create_atoms 2 region plane
Created 60025 atoms
set type 2 mass 0.1
60025 settings made for mass
group small type 2
60025 atoms in group small
velocity small create 1.0 593849 loop geom
# delete overlaps
# must set 1-2 cutoff to non-zero value
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 0.5
delete_atoms overlap 0.5 small big
Deleted 16305 atoms, new total = 43990
# SRD run
reset_timestep 0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
communicate multi group big vel yes
neigh_modify include big
# no pairwise interactions with small particles
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 0.0
# use fix SRD to push small particles out from inside big ones
# if comment out, big particles won't see SRD particles
timestep 0.001
fix 1 big rigid molecule
30 rigid bodies with 270 atoms
fix 2 small srd 20 big 1.0 0.25 49894 shift yes 54979 search 0.2 collision slip inside ignore overlap yes
fix 3 all enforce2d
# diagnostics
compute tbig big temp/sphere
variable pebig equal pe*atoms/count(big)
variable ebig equal etotal*atoms/count(big)
thermo_style custom step temp f_2[8] etotal v_pebig v_ebig press f_2[1] f_2[2] f_2[3] f_2[4] f_2[5] f_2[6] f_2[7] f_2[8] f_2[9] f_2[10] f_2[11] f_2[12]
thermo_modify temp tbig
WARNING: Temperature for thermo pressure is not for group all (../thermo.cpp:439)
thermo 1000
#dump 1 all atom 1000 dump.star
#dump 1 all image 1000 image.*.jpg type type zoom 1.6
#dump_modify 1 pad 6 adiam 1 1 adiam 2 0.2
run 100000
SRD info:
SRD/big particles = 43720 270
big particle diameter max/min = 1 1
SRD temperature & lamda = 1 0.0632456
SRD max distance & max velocity = 0.252982 12.6491
SRD grid counts: 106 106 1
SRD grid size: request, actual (xyz) = 0.25, 0.250819 0.250819 1
SRD per actual grid cell = 5.55866
SRD viscosity = 0.235582
big/SRD mass density ratio = 0.144099
WARNING: Fix srd grid size > 1/4 of big particle diameter (../fix_srd.cpp:2875)
# of rescaled SRD velocities = 0
ave/max small velocity = 4.19085 7.72582
ave/max big velocity = 2.20262 5.4168
Memory usage per processor = 7.83759 Mbytes
Step Temp 2[8] TotEng pebig ebig Press 2[1] 2[2] 2[3] 2[4] 2[5] 2[6] 2[7] 2[8] 2[9] 2[10] 2[11] 2[12]
0 18.299191 0 0.25067776 37.859815 40.841906 5.7680841 0 0 0 0 0 0 0 0 0 0 0 0
1000 1.3155168 8201 0.22979575 37.225304 37.439685 8.9588433 14577 56 56 0 687 57680 11236 8201 1.0283377 0 5 0
2000 0.82290815 8187 0.22918912 37.206746 37.340849 6.0607727 14326 56 56 0 1524 57680 11236 8187 1.0127464 0 5 0
3000 0.99928391 8168 0.22921842 37.182777 37.345623 10.268998 14389 68 68 0 2274 57680 11236 8168 1.023497 0 5 0
4000 1.0348458 8200 0.22920437 37.174693 37.343335 9.8954805 14335 51 51 0 3047 57680 11236 8200 1.0156984 0 5 0
5000 1.1496765 8199 0.22939095 37.186379 37.373733 8.0143223 14444 74 74 0 3868 57680 11236 8199 1.0104858 0 5 0
6000 0.91767381 8225 0.2292466 37.200668 37.350215 8.6189298 14131 61 61 0 4594 57680 11236 8225 1.0147272 0 5 0
7000 0.74806752 8203 0.22888634 37.169611 37.291518 9.8613421 14121 53 53 1 5322 57680 11236 8203 1.0198289 0 19 0
8000 1.1177867 8213 0.22924801 37.168286 37.350444 7.911691 14380 62 62 0 6024 57680 11236 8213 1.0066799 0 19 0
9000 1.1258411 8209 0.22939703 37.191253 37.374723 6.7377179 14350 51 51 0 6821 57680 11236 8209 1.0081225 0 19 0
10000 1.0597344 8182 0.22919436 37.169006 37.341703 9.7769156 14372 48 48 0 7543 57680 11236 8182 1.0099373 0 19 0
11000 1.3337199 8159 0.22943046 37.162823 37.38017 6.0897666 14359 58 58 0 8245 57680 11236 8159 1.0046221 0 19 0
12000 1.0280975 8184 0.22894286 37.133186 37.300728 2.555317 14162 51 51 0 8935 57680 11236 8184 1.0010306 0 19 0
13000 1.2318533 8190 0.22921487 37.144299 37.345045 6.1739989 14285 64 64 0 9634 57680 11236 8190 1.0138112 0 19 0
14000 1.1851991 8215 0.22915708 37.142485 37.335629 7.8380175 14151 66 66 0 10407 57680 11236 8215 1.0008011 0 19 0
15000 0.89853218 8214 0.22897294 37.159201 37.305628 4.0743953 14077 61 61 0 11069 57680 11236 8214 1.0041341 0 19 0
16000 1.0295974 8176 0.22891205 37.127921 37.295708 8.7527117 14000 48 48 0 11832 57680 11236 8176 1.0089653 0 23 0
17000 1.1793664 8219 0.22906268 37.128056 37.320249 5.2583304 14268 57 57 0 12490 57680 11236 8219 1.0157079 0 23 0
18000 0.90486919 8180 0.22891008 37.147927 37.295387 4.8498125 13939 64 64 0 13154 57680 11236 8180 0.99700921 0 23 0
19000 0.81486319 8219 0.22881911 37.147773 37.280566 1.6338041 13934 59 59 1 13861 57680 11236 8219 1.0099189 0 23 0
20000 1.0146287 8180 0.22896254 37.138587 37.303934 2.1682593 14057 45 45 0 14593 57680 11236 8180 1.0009466 0 23 0
21000 1.3031463 8164 0.22920646 37.131309 37.343674 6.570906 13944 65 65 0 15323 57680 11236 8164 0.98955239 0 23 0
22000 1.1258162 8174 0.22915084 37.151147 37.334613 9.0817519 14052 56 56 0 16042 57680 11236 8174 1.0083228 0 23 0
23000 1.1421222 8217 0.22906722 37.134865 37.320989 8.0529401 14143 66 66 0 16797 57680 11236 8217 1.0111161 0 23 0
24000 0.77402445 8168 0.22883056 37.156294 37.282431 6.8609899 14258 62 62 0 17559 57680 11236 8168 1.0080535 0 23 0
25000 0.94006632 8215 0.22879533 37.123496 37.276692 3.8570839 14446 48 48 0 18267 57680 11236 8215 1.0087899 0 23 0
26000 1.2808563 8192 0.22933112 37.155253 37.363985 5.3591484 13960 54 54 1 19004 57680 11236 8192 1.0101851 0 79 0
27000 1.0266075 8171 0.22898652 37.140542 37.307841 2.9961422 14049 56 56 1 19683 57680 11236 8171 0.99950217 0 79 0
28000 0.98185838 8199 0.22867129 37.096475 37.256481 12.080529 14106 61 62 0 20417 57680 11236 8199 1.0026548 0 79 0
29000 0.96921352 8192 0.22861802 37.089856 37.247802 8.257635 13891 63 64 0 21169 57680 11236 8192 1.0128041 0 79 0
30000 1.0107665 8223 0.22866616 37.090929 37.255646 9.4280788 14021 61 61 1 21863 57680 11236 8223 1.0034994 0 79 0
31000 1.1859435 8195 0.22877222 37.079661 37.272926 7.3863353 14093 78 78 0 22480 57680 11236 8195 0.9945582 0 79 0
32000 1.2951374 8172 0.22897044 37.094161 37.30522 6.4322363 13792 52 52 0 23151 57680 11236 8172 1.0026092 0 79 0
33000 0.82899919 8176 0.22823894 37.050944 37.186041 0.60210493 13797 62 62 0 23903 57680 11236 8176 1.0051075 0 79 0
34000 0.96685714 8165 0.22828901 37.036637 37.194199 8.1409866 13837 50 50 0 24591 57680 11236 8165 1.0157895 0 79 0
35000 0.95645753 8206 0.22833008 37.045022 37.200889 8.9956195 14130 63 63 0 25340 57680 11236 8206 1.0030424 0 79 0
36000 1.2933379 8197 0.22879201 37.065383 37.27615 8.3938627 13978 80 80 0 26085 57680 11236 8197 1.0094497 0 79 0
37000 1.1511434 8176 0.22861016 37.058929 37.246523 4.9481995 13856 49 49 1 26759 57680 11236 8176 1.0063722 0 79 0
38000 0.93310274 8181 0.22837532 37.056199 37.20826 5.9575015 13949 51 51 0 27445 57680 11236 8181 1.0152607 0 79 0
39000 1.0352467 8215 0.22861437 37.078501 37.247208 9.8519532 14109 69 69 0 28245 57680 11236 8215 1.005523 0 79 0
40000 1.2696427 8193 0.22885162 37.078957 37.285862 8.4680587 13931 54 55 0 28960 57680 11236 8193 1.0014074 0 79 0
41000 1.3270531 8193 0.22909573 37.109374 37.325634 4.3731598 14115 62 62 0 29668 57680 11236 8193 1.0022493 0 79 0
42000 1.3650738 8183 0.22911257 37.105921 37.328377 6.6345101 13659 52 52 0 30400 57680 11236 8183 1.0145834 0 79 0
43000 1.0780084 8209 0.22855316 37.06156 37.237235 5.5147087 14062 52 52 0 31092 57680 11236 8209 1.0031365 0 79 0
44000 0.96146724 8193 0.22846546 37.066263 37.222946 9.7967735 13699 56 56 0 31771 57680 11236 8193 1.0072364 0 79 0
45000 1.1733618 8200 0.2286481 37.061489 37.252703 9.8058881 13967 51 52 0 32532 57680 11236 8200 1.0068663 0 79 0
46000 0.94784989 8211 0.22845052 37.066048 37.220513 5.8834756 13944 34 35 0 33213 57680 11236 8211 1.0080225 0 79 0
47000 1.2354771 8207 0.22859598 37.042876 37.244213 6.1458407 13851 56 56 0 33886 57680 11236 8207 1.0151781 0 79 0
48000 1.0834097 8164 0.22835968 37.029157 37.205712 6.7378814 13944 56 56 0 34578 57680 11236 8164 1.0128575 0 79 0
49000 0.9386786 8208 0.22837254 37.054838 37.207808 5.0809383 13618 53 53 0 35328 57680 11236 8208 1.0045622 0 79 0
50000 1.047176 8154 0.22859126 37.072791 37.243442 4.8016981 13847 55 55 0 36059 57680 11236 8154 1.0005378 0 79 0
51000 1.0943417 8222 0.22844191 37.040773 37.21911 7.9359541 13897 61 61 0 36711 57680 11236 8222 1.0036343 0 79 0
52000 1.2811259 8200 0.22871149 37.054255 37.263031 4.9958687 13744 56 56 0 37402 57680 11236 8200 1.0049188 0 79 0
53000 1.0997162 8218 0.22853173 37.05453 37.233743 4.6772882 13967 65 65 0 38082 57680 11236 8218 1.0054533 0 79 0
54000 0.97625718 8185 0.2283994 37.05309 37.212184 5.9026406 13741 69 69 0 38761 57680 11236 8185 0.99393606 0 79 0
55000 0.97919465 8192 0.22835234 37.044944 37.204516 6.8419755 13990 58 58 0 39453 57680 11236 8192 1.0103052 0 79 0
56000 0.86300824 8181 0.22842164 37.075169 37.215808 6.8288348 13751 57 57 0 40118 57680 11236 8181 0.99683099 0 79 0
57000 0.89429559 8181 0.22851743 37.085677 37.231414 9.378853 13746 72 72 0 40730 57680 11236 8181 1.0086882 0 79 0
58000 1.0853441 8172 0.22857719 37.064279 37.24115 7.4580035 13869 58 58 1 41462 57680 11236 8172 1.0092141 0 79 0
59000 1.0149632 8199 0.22868187 37.092803 37.258205 8.2913359 13720 54 54 0 42159 57680 11236 8199 1.0082738 0 79 0
60000 0.99214038 8241 0.22871146 37.101344 37.263026 7.9905647 13776 67 67 0 42858 57680 11236 8241 1.0027929 0 79 0
61000 0.95721926 8211 0.2285241 37.07651 37.232501 11.837258 13907 62 62 0 43565 57680 11236 8211 1.0201117 0 79 0
62000 1.1457294 8195 0.22866514 37.068769 37.25548 9.7934944 13796 40 40 1 44292 57680 11236 8195 1.005192 0 79 0
63000 1.1096177 8193 0.22863037 37.068989 37.249815 8.0809047 13984 54 55 0 44986 57680 11236 8193 1.0076438 0 79 0
64000 1.1643495 8168 0.22860016 37.055147 37.244893 4.7376475 13857 50 50 1 45657 57680 11236 8168 1.002011 0 79 0
65000 1.2172585 8201 0.22845803 37.023368 37.221736 8.9744616 13610 53 53 0 46295 57680 11236 8201 1.0033999 0 79 0
66000 0.99924611 8189 0.22807344 36.996237 37.159077 8.5854786 13849 49 49 0 46944 57680 11236 8189 1.015528 0 79 0
67000 1.2543542 8165 0.22831075 36.993327 37.19774 5.0081411 13906 57 58 1 47666 57680 11236 8165 1.0189753 0 79 0
68000 1.0360874 8205 0.22812149 36.99806 37.166904 2.2431831 13650 52 52 0 48363 57680 11236 8205 1.0076864 0 79 0
69000 0.69755552 8218 0.22766249 36.978446 37.092121 4.2769033 13654 51 51 0 49044 57680 11236 8218 1.009863 0 79 0
70000 0.90921994 8167 0.22796236 36.992809 37.140978 9.3555738 13927 52 52 0 49793 57680 11236 8167 1.0015807 0 79 0
71000 0.95608837 8177 0.22797832 36.987772 37.143579 8.0740708 13612 63 63 0 50464 57680 11236 8177 1.0097491 0 79 0
72000 1.0374855 8192 0.22803745 36.98414 37.153212 4.4789755 13561 52 52 0 51127 57680 11236 8192 0.99496724 0 79 0
73000 0.90636338 8214 0.22783374 36.97232 37.120024 9.3364171 13598 53 53 0 51799 57680 11236 8214 1.0078219 0 79 0
74000 1.1869638 8179 0.22834625 37.010094 37.203525 4.8292929 13771 51 51 0 52500 57680 11236 8179 0.98865277 0 79 0
75000 1.0177351 8188 0.22826506 37.024443 37.190296 5.1808505 13508 55 55 0 53183 57680 11236 8188 0.99903558 0 79 0
76000 1.062496 8230 0.22820987 37.008158 37.181305 8.1032524 13769 58 58 0 53864 57680 11236 8230 0.99397021 0 79 0
77000 0.9762027 8131 0.22831271 37.038975 37.19806 10.309798 13787 47 47 0 54553 57680 11236 8131 1.0066289 0 79 0
78000 1.0688533 8160 0.22855677 37.06364 37.237824 7.4800503 13224 58 58 0 55248 57680 11236 8160 1.0073477 0 79 0
79000 1.2598513 8151 0.22844942 37.015024 37.220333 7.4724924 13859 63 63 0 55911 57680 11236 8151 1.0065168 0 79 0
80000 1.1044575 8176 0.22834306 37.023019 37.203005 3.0142586 13579 62 62 0 56611 57680 11236 8176 1.0008704 0 79 0
81000 1.2756491 8205 0.22853748 37.026797 37.23468 6.434309 13702 52 52 0 57313 57680 11236 8205 1.0058366 0 79 0
82000 1.1276926 8193 0.22857049 37.056287 37.240059 11.825248 13621 63 63 0 58011 57680 11236 8193 1.0113537 0 79 0
83000 0.96941691 8173 0.22824759 37.029471 37.18745 2.9152799 13478 48 49 0 58712 57680 11236 8173 1.0168819 0 79 0
84000 1.0770325 8186 0.22821292 37.006285 37.181802 5.0626072 13733 60 61 0 59433 57680 11236 8186 1.0177757 0 79 0
85000 1.1959489 8175 0.2284266 37.021719 37.216615 5.3844747 13964 55 55 0 60226 57680 11236 8175 1.007371 0 79 0
86000 1.035779 8204 0.22845952 37.053185 37.221979 7.5428558 13713 64 64 0 60982 57680 11236 8204 1.0110801 0 79 0
87000 1.1878573 8204 0.22836412 37.012859 37.206436 3.1124447 13655 47 47 0 61668 57680 11236 8204 1.0034968 0 79 0
88000 0.98402924 8206 0.22815704 37.012337 37.172698 10.613422 13947 56 56 0 62325 57680 11236 8206 1.0081468 0 79 0
89000 0.83763529 8199 0.22807071 37.022127 37.158631 5.4628281 13966 50 50 0 63012 57680 11236 8199 0.99956322 0 79 0
90000 1.1976672 8187 0.22849599 37.032745 37.227921 5.5583762 13903 43 43 0 63697 57680 11236 8187 1.0047573 0 79 0
91000 0.94863742 8188 0.22818061 37.021945 37.176537 5.157692 13859 46 46 0 64382 57680 11236 8188 1.0041894 0 79 0
92000 1.2038427 8172 0.22846103 37.026042 37.222224 6.685059 13675 57 57 0 65081 57680 11236 8172 1.0036906 0 79 0
93000 0.95835618 8210 0.22828548 37.037446 37.193622 8.0038227 13826 57 58 0 65761 57680 11236 8210 1.0076646 0 79 0
94000 1.2451267 8197 0.22856138 37.035665 37.238575 3.6305666 13918 62 62 0 66462 57680 11236 8197 1.0110913 0 79 0
95000 0.82128549 8206 0.22809174 37.028218 37.162057 8.7024273 13817 44 44 0 67185 57680 11236 8206 1.0032224 0 79 0
96000 0.83473128 8199 0.22820935 37.04519 37.18122 8.1258155 14167 54 54 0 67912 57680 11236 8199 1.0035561 0 79 0
97000 0.98924493 8196 0.22847944 37.064014 37.225225 8.9964529 14039 48 51 0 68608 57680 11236 8196 1.0007515 0 79 0
98000 1.1749668 8211 0.22879851 37.085733 37.27721 7.0283734 13850 58 59 0 69305 57680 11236 8211 0.99525651 0 79 0
99000 1.3269705 8163 0.22890919 37.078994 37.295241 9.7962294 13885 63 63 0 69994 57680 11236 8163 1.011015 0 107 0
100000 1.0020436 8171 0.22867221 37.093336 37.256632 4.9455304 14102 55 56 0 70707 57680 11236 8171 1.0072865 0 107 0
Loop time of 33.9872 on 8 procs for 100000 steps with 43990 atoms
Pair time (%) = 0.379981 (1.11801)
Neigh time (%) = 2.01801 (5.93757)
Comm time (%) = 2.30255 (6.77476)
Outpt time (%) = 0.00694308 (0.0204285)
Other time (%) = 29.2797 (86.1492)
Nlocal: 5498.75 ave 6445 max 4431 min
Histogram: 1 0 0 2 1 1 0 1 1 1
Nghost: 62.5 ave 70 max 57 min
Histogram: 3 0 0 1 1 0 1 0 1 1
Neighs: 207 ave 365 max 92 min
Histogram: 1 1 1 1 2 0 1 0 0 1
Total # of neighbors = 1656
Ave neighs/atom = 0.0376449
Neighbor list builds = 5000
Dangerous builds = 0
Please see the log.cite file for references relevant to this simulation

291
examples/ASPHERE/star/log.star.mp.1Feb14 Normal file
View File

@ -0,0 +1,291 @@
LAMMPS (1 Feb 2014)
# SRD viscosity demo - rigid star particles
units lj
atom_style sphere
atom_modify map array first big
dimension 2
# read in clusters of rigid bodies
fix molprop all property/atom mol
read_data data.star fix molprop NULL Molecules
orthogonal box = (-13.2934 -13.2934 -0.5) to (13.2934 13.2934 0.5)
4 by 2 by 1 MPI processor grid
reading atoms ...
270 atoms
set type 1 mass 1.0
270 settings made for mass
group big type 1
270 atoms in group big
velocity big create 1.44 87287 loop geom
# equilibrate big particles
pair_style soft 1.12
pair_coeff 1 1 0.0
pair_coeff 2 2 0.0 0.0
pair_coeff 1 2 0.0 0.0
variable prefactor equal ramp(0,60)
fix soft all adapt 1 pair soft a * * v_prefactor
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
fix 1 big rigid molecule
30 rigid bodies with 270 atoms
fix 2 all enforce2d
#dump 1 all atom 10 dump.star.equil
thermo 100
run 1000
Memory usage per processor = 2.64859 Mbytes
Step Temp E_pair E_mol TotEng Press
0 0 0 0 0 0.24721146
100 0 3.291475 0 3.291475 0.69158089
200 0 4.6176595 0 4.6176595 0.71737945
300 0 7.0921814 0 7.0921814 1.2163497
400 0 8.3666709 0 8.3666709 0.50645957
500 0 10.630838 0 10.630838 0.72764487
600 0 12.42157 0 12.42157 1.0130599
700 0 14.349074 0 14.349074 0.74795509
800 0 16.464746 0 16.464746 1.054549
900 0 18.253108 0 18.253108 0.51551753
1000 0 20.699563 0 20.699563 1.8084662
Loop time of 0.0484946 on 8 procs for 1000 steps with 270 atoms
Pair time (%) = 0.00250834 (5.17242)
Neigh time (%) = 0.000999629 (2.06132)
Comm time (%) = 0.0154877 (31.9369)
Outpt time (%) = 0.000273734 (0.564464)
Other time (%) = 0.0292252 (60.2649)
Nlocal: 33.75 ave 42 max 22 min
Histogram: 1 0 1 0 0 2 0 1 1 2
Nghost: 23.375 ave 30 max 14 min
Histogram: 1 0 1 0 0 2 1 0 1 2
Neighs: 52.25 ave 69 max 27 min
Histogram: 1 0 1 0 1 0 1 1 1 2
Total # of neighbors = 418
Ave neighs/atom = 1.54815
Neighbor list builds = 176
Dangerous builds = 0
#undump 1
unfix soft
unfix 1
unfix 2
# add small particles as hi density lattice
region plane block INF INF INF INF -0.001 0.001 units box
lattice sq 85.0
Lattice spacing in x,y,z = 0.108465 0.108465 0.108465
create_atoms 2 region plane
Created 60025 atoms
set type 2 mass 0.1
60025 settings made for mass
group small type 2
60025 atoms in group small
velocity small create 1.0 593849 loop geom
# delete overlaps
# must set 1-2 cutoff to non-zero value
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 0.5
delete_atoms overlap 0.5 small big
Deleted 16305 atoms, new total = 43990
# SRD run
reset_timestep 0
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
communicate multi group big vel yes
neigh_modify include big
# no pairwise interactions with small particles
pair_style lj/cut 2.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 0.0
# use fix SRD to push small particles out from inside big ones
# if comment out, big particles won't see SRD particles
timestep 0.001
fix 1 big rigid molecule
30 rigid bodies with 270 atoms
fix 2 small srd 20 big 1.0 0.25 49894 shift yes 54979 search 0.2 collision slip inside ignore overlap yes
fix 3 small viscosity 10 x y 50
fix 4 all enforce2d
# diagnostics
compute tbig big temp/sphere
variable pebig equal pe*atoms/count(big)
variable ebig equal etotal*atoms/count(big)
thermo_style custom step temp f_2[8] etotal v_pebig v_ebig press f_2[1] f_2[2] f_2[3] f_2[4] f_2[5] f_2[6] f_2[7] f_2[8] f_2[9] f_2[10] f_2[11] f_2[12]
thermo_modify temp tbig
WARNING: Temperature for thermo pressure is not for group all (../thermo.cpp:439)
thermo 1000
#dump 1 all atom 1000 dump.star.mp
#dump 1 all image 1000 image.*.jpg type type zoom 1.6
#dump_modify 1 pad 6 adiam 1 1 adiam 2 0.2
run 100000
SRD info:
SRD/big particles = 43720 270
big particle diameter max/min = 1 1
SRD temperature & lamda = 1 0.0632456
SRD max distance & max velocity = 0.252982 12.6491
SRD grid counts: 106 106 1
SRD grid size: request, actual (xyz) = 0.25, 0.250819 0.250819 1
SRD per actual grid cell = 5.55866
SRD viscosity = 0.235582
big/SRD mass density ratio = 0.144099
WARNING: Fix srd grid size > 1/4 of big particle diameter (../fix_srd.cpp:2875)
# of rescaled SRD velocities = 0
ave/max small velocity = 4.19085 7.72582
ave/max big velocity = 2.20262 5.4168
Memory usage per processor = 7.83759 Mbytes
Step Temp 2[8] TotEng pebig ebig Press 2[1] 2[2] 2[3] 2[4] 2[5] 2[6] 2[7] 2[8] 2[9] 2[10] 2[11] 2[12]
0 18.299191 0 0.25067776 37.859815 40.841906 5.7680841 0 0 0 0 0 0 0 0 0 0 0 0
1000 1.3197888 8196 0.2296966 37.208454 37.42353 7.4694369 14388 59 59 0 707 57680 11236 8196 1.0259754 0 5 0
2000 1.5789784 8192 0.22993088 37.204387 37.461702 7.6252927 14352 56 56 0 1521 57680 11236 8192 1.0168955 0 7 0
3000 1.286186 8153 0.2296947 37.213621 37.423222 8.8624394 14533 73 73 0 2299 57680 11236 8153 1.0109574 0 7 0
4000 2.0862082 8194 0.23044604 37.20566 37.545635 7.8064365 14228 48 48 0 3068 57680 11236 8194 1.0046224 0 7 0
5000 1.892401 8172 0.23019301 37.196017 37.504409 7.9041579 14290 55 55 0 3866 57680 11236 8172 1.005673 0 7 0
6000 2.1063411 8138 0.23041636 37.197543 37.540799 4.8468885 13803 47 47 0 4638 57680 11236 8138 1.0083722 0 7 0
7000 2.3549938 8174 0.23071159 37.205123 37.5889 6.5703154 13773 58 58 0 5398 57680 11236 8174 1.0039365 0 7 0
8000 3.0660305 8216 0.23155654 37.226914 37.726563 11.977404 13898 65 65 0 6180 57680 11236 8216 0.99586836 0 7 0
9000 3.1510761 8203 0.23143756 37.19367 37.707178 7.2978028 13815 54 54 0 6985 57680 11236 8203 0.99794556 0 7 0
10000 2.7472522 8196 0.23096076 37.181796 37.629496 11.336568 14012 59 59 0 7812 57680 11236 8196 0.99210203 0 11 0
11000 2.1776587 8186 0.23038876 37.181424 37.536302 9.2106614 14063 64 64 0 8627 57680 11236 8186 1.003993 0 11 0
12000 3.0146734 8196 0.2312123 37.179198 37.670478 5.9303705 14077 60 60 0 9436 57680 11236 8196 0.98496016 0 23 0
13000 3.0029829 8179 0.23150759 37.229214 37.718589 7.3329148 14116 74 74 0 10339 57680 11236 8179 1.0039911 0 23 0
14000 3.2942972 8191 0.23166868 37.207985 37.744833 9.6428846 14019 56 57 0 11290 57680 11236 8191 0.99410574 0 23 0
15000 3.8480961 8213 0.23228155 37.217589 37.844686 2.1025549 13746 68 68 0 12175 57680 11236 8213 1.0011325 0 23 0
16000 3.2573043 8195 0.2316418 37.209634 37.740454 8.1255708 13682 55 55 0 13104 57680 11236 8195 1.0061399 0 23 0
17000 3.3579037 8186 0.2315443 37.177356 37.72457 9.8606039 13833 54 54 0 14034 57680 11236 8186 1.004156 0 23 0
18000 4.7903055 8194 0.23307181 37.192799 37.973441 3.8090185 13914 55 55 1 15012 57680 11236 8194 0.98633014 0 23 0
19000 4.4696839 8183 0.23293212 37.222288 37.950681 7.1931333 13499 54 54 0 15960 57680 11236 8183 0.98359656 0 23 0
20000 5.0239019 8126 0.23355887 37.234085 38.052794 5.1504203 13600 58 58 0 16972 57680 11236 8126 0.98555755 0 23 0
21000 4.6595948 8155 0.23275722 37.162845 37.922186 9.0375076 13965 73 74 0 17966 57680 11236 8155 0.99574348 0 23 0
22000 4.043444 8124 0.23247599 37.217434 37.876366 3.9265975 13584 53 53 2 18989 57680 11236 8124 0.9943952 0 23 0
23000 4.9981028 8127 0.23314549 37.170939 37.985445 4.2352191 13589 47 47 0 19958 57680 11236 8127 0.99276741 0 47 0
24000 4.7542655 8145 0.23295434 37.179533 37.954302 7.365415 13753 51 51 0 20955 57680 11236 8145 0.99024306 0 47 0
25000 4.2987054 8177 0.23252979 37.184601 37.885131 10.735771 13536 58 58 0 21994 57680 11236 8177 0.99215977 0 47 0
26000 4.4136261 8170 0.23331629 37.294015 38.013272 14.25925 13842 58 58 0 22979 57680 11236 8170 0.99283174 0 47 0
27000 4.0546007 8173 0.23283443 37.274016 37.934766 6.22969 13877 51 51 0 23997 57680 11236 8173 0.98496361 0 47 0
28000 4.5919614 8151 0.23364384 37.318319 38.066638 4.9917488 13836 60 61 1 25071 57680 11236 8151 0.99390956 0 47 0
29000 4.8025535 8189 0.23313849 37.201667 37.984305 4.7895172 13978 46 47 0 26121 57680 11236 8189 0.99018295 0 47 0
30000 4.2835226 8139 0.23247583 37.178284 37.87634 8.994878 13988 61 61 0 27141 57680 11236 8139 0.98735471 0 47 0
31000 5.4049428 8133 0.23373767 37.201121 38.081927 7.6419291 13771 58 58 0 28207 57680 11236 8133 0.98354883 0 47 0
32000 5.0299584 8177 0.23336211 37.201041 38.020738 8.1206741 14036 55 55 0 29303 57680 11236 8177 0.98239013 0 47 0
33000 4.2993816 8158 0.23320229 37.29406 37.9947 10.745167 13583 51 51 0 30344 57680 11236 8158 0.98952524 0 47 0
34000 5.8334257 8179 0.23458355 37.269109 38.219741 11.094625 13921 62 62 0 31432 57680 11236 8179 0.98822406 0 47 0
35000 6.481318 8168 0.23501133 37.233224 38.289439 7.9397445 13942 56 56 0 32518 57680 11236 8168 0.98830425 0 47 0
36000 4.967891 8186 0.23366927 37.261199 38.070782 8.3979019 13952 51 51 0 33626 57680 11236 8186 0.99579665 0 47 0
37000 5.3458808 8145 0.2337339 37.210132 38.081312 5.4177865 13659 55 55 0 34715 57680 11236 8145 1.0009582 0 47 0
38000 6.3597955 8159 0.23479285 37.217431 38.253842 6.1408411 13754 55 56 0 35754 57680 11236 8159 0.98853252 0 47 0
39000 6.1217173 8195 0.23447066 37.203737 38.20135 9.6016075 13756 70 70 0 36807 57680 11236 8195 0.98652299 0 47 0
40000 5.0218059 8159 0.23370899 37.258886 38.077254 6.6460004 13603 56 56 1 37933 57680 11236 8159 0.98674585 0 47 0
41000 5.5873722 8138 0.23429215 37.26173 38.172265 1.4183057 13634 51 52 1 39050 57680 11236 8138 0.98833769 0 47 0
42000 5.5502013 8170 0.23405343 37.228895 38.133372 7.8338637 13913 72 72 0 40105 57680 11236 8170 0.99124609 0 47 0
43000 4.6854143 8161 0.23326283 37.241013 38.004562 7.857085 13759 48 48 0 41187 57680 11236 8161 0.99649848 0 47 0
44000 6.3172352 8180 0.2348893 37.240081 38.269557 12.59288 13810 55 55 0 42292 57680 11236 8180 0.99164669 0 47 0
45000 6.1213322 8166 0.23449878 37.208381 38.205931 6.3760889 13854 48 49 0 43406 57680 11236 8166 0.98981447 0 47 0
46000 5.0015874 8195 0.23354809 37.235965 38.051038 10.992836 13962 62 62 0 44527 57680 11236 8195 0.99128184 0 47 0
47000 5.0499615 8190 0.23376151 37.262854 38.085811 8.7985849 14105 64 64 0 45653 57680 11236 8190 0.99020272 0 47 0
48000 5.202255 8170 0.23400688 37.278013 38.125788 4.7389097 13754 51 51 0 46705 57680 11236 8170 0.98792999 0 47 0
49000 5.6144741 8167 0.23429848 37.258346 38.173297 6.3247626 13953 51 51 1 47846 57680 11236 8167 0.99270232 0 47 0
50000 6.3027505 8145 0.2347122 37.213588 38.240703 6.4424496 13843 59 59 0 49061 57680 11236 8145 0.99810926 0 47 0
51000 5.2050769 8154 0.23390178 37.26043 38.108665 6.6143765 13999 51 51 0 50075 57680 11236 8154 0.99525977 0 47 0
52000 5.4749884 8175 0.23418567 37.262697 38.154917 8.7256782 13916 51 51 0 51249 57680 11236 8175 0.99508289 0 47 0
53000 5.1367735 8158 0.23393299 37.276645 38.113748 9.3248422 13910 41 41 0 52381 57680 11236 8158 0.98956839 0 47 0
54000 4.8680584 8211 0.23402832 37.335968 38.129281 4.8979582 14214 63 64 0 53439 57680 11236 8211 0.98838908 0 47 0
55000 4.6223848 8167 0.23347097 37.285196 38.038473 6.1898955 14162 46 46 0 54504 57680 11236 8167 0.97935238 0 47 0
56000 4.9075747 8168 0.23409393 37.340218 38.139971 12.20686 13980 60 60 0 55540 57680 11236 8168 0.99885966 0 47 0
57000 4.8736147 8143 0.23354696 37.256636 38.050854 3.9348529 13857 59 59 0 56680 57680 11236 8143 0.98434054 0 47 0
58000 5.6134201 8162 0.23393276 37.198932 38.113711 3.7493067 13866 55 55 2 57767 57680 11236 8162 0.99041155 0 47 0
59000 5.440371 8197 0.23353392 37.162151 38.04873 4.1297786 13748 67 67 0 58894 57680 11236 8197 0.98824623 0 47 0
60000 5.2754075 8148 0.23359083 37.198307 38.058003 7.1427242 13781 55 55 0 60028 57680 11236 8148 0.9831827 0 47 0
61000 5.8228233 8153 0.23438274 37.23812 38.187025 8.8839625 14019 48 48 0 61204 57680 11236 8153 0.98566731 0 47 0
62000 5.7623285 8129 0.23480056 37.316052 38.255098 9.992341 13791 67 68 0 62385 57680 11236 8129 0.97338459 0 47 0
63000 6.095315 8178 0.23455792 37.222257 38.215567 10.722535 14210 54 54 0 63448 57680 11236 8178 0.98135973 0 47 0
64000 5.2494942 8189 0.23380158 37.236866 38.092339 4.4010329 13869 51 51 0 64538 57680 11236 8189 0.98071569 0 47 0
65000 5.6534858 8151 0.23405328 37.212039 38.133348 8.9144953 14065 53 53 0 65661 57680 11236 8151 0.99589337 0 47 0
66000 5.3999138 8165 0.23453061 37.33113 38.211116 7.3102657 13994 61 62 0 66785 57680 11236 8165 0.98715854 0 47 0
67000 6.3619419 8183 0.23459518 37.184875 38.221636 4.0425684 13859 49 49 0 67944 57680 11236 8183 0.98479883 0 47 0
68000 6.2170116 8157 0.23470153 37.225822 38.238965 9.3090925 13598 50 50 0 69121 57680 11236 8157 0.98683618 0 47 0
69000 6.8545121 8179 0.23542272 37.239433 38.356465 9.4175179 13911 59 59 0 70316 57680 11236 8179 0.99828763 0 47 0
70000 5.1476232 8154 0.23389653 37.268937 38.107809 7.6438157 13882 55 55 0 71380 57680 11236 8154 0.99088995 0 47 0
71000 5.9284504 8161 0.23499148 37.320087 38.286205 9.1242845 13968 60 60 0 72478 57680 11236 8161 0.98627289 0 47 0
72000 5.8044136 8180 0.23473312 37.298207 38.244111 6.6084298 14066 56 56 0 73653 57680 11236 8180 0.97959812 0 47 0
73000 5.3690655 8148 0.23424676 37.289911 38.16487 8.7025002 14433 55 55 1 74740 57680 11236 8148 0.97949605 0 47 0
74000 4.8144157 8192 0.23361625 37.277572 38.062144 7.6685178 14425 58 58 1 75918 57680 11236 8192 0.98949387 0 47 0
75000 4.9424688 8143 0.23362551 37.258213 38.063652 5.5882668 13938 42 42 0 77080 57680 11236 8143 0.99602516 0 47 0
76000 6.3017203 8178 0.23500247 37.261048 38.287995 5.2760401 13952 47 47 0 78240 57680 11236 8178 0.98484281 0 47 0
77000 5.3086333 8208 0.23389726 37.242817 38.107928 4.034731 14052 47 47 0 79331 57680 11236 8208 0.9896033 0 47 0
78000 6.1152265 8157 0.23546115 37.366171 38.362726 11.521489 13980 58 58 0 80446 57680 11236 8157 0.98562154 0 47 0
79000 6.5034669 8173 0.23569944 37.341726 38.40155 5.0939756 14030 74 74 0 81663 57680 11236 8173 0.98837711 0 47 0
80000 6.1688245 8182 0.23495003 37.274162 38.279451 6.3337593 14292 50 50 0 82872 57680 11236 8182 0.99316806 0 47 0
81000 6.3818913 8189 0.23516767 37.274899 38.314911 9.0232294 14259 56 56 0 84072 57680 11236 8189 0.98467136 0 47 0
82000 6.502454 8171 0.23498158 37.224932 38.284591 7.6685939 13970 63 63 0 85190 57680 11236 8171 0.99511647 0 47 0
83000 7.2245962 8142 0.23564338 37.215075 38.392416 6.0363083 13789 61 61 0 86416 57680 11236 8142 0.99223565 0 47 0
84000 6.3529888 8192 0.23477435 37.215526 38.250828 7.2109855 14062 52 52 0 87555 57680 11236 8192 0.98623425 0 47 0
85000 6.5185235 8177 0.23516364 37.251976 38.314254 5.266764 14325 52 53 0 88732 57680 11236 8177 0.97895516 0 47 0
86000 5.742959 8170 0.23432816 37.242243 38.178133 6.0834874 13922 48 48 0 89836 57680 11236 8170 0.99203022 0 47 0
87000 6.548676 8177 0.23495538 37.213132 38.280324 6.2502744 14147 48 48 0 91035 57680 11236 8177 0.98829721 0 47 0
88000 5.9500587 8193 0.23446634 37.231006 38.200645 6.2196893 13960 44 44 0 92206 57680 11236 8193 0.98956808 0 47 0
89000 5.0591447 8174 0.23386574 37.278339 38.102793 11.164066 14071 54 54 0 93437 57680 11236 8174 0.99827748 0 47 0
90000 6.4930795 8181 0.23544182 37.301445 38.359576 6.6167028 14183 63 63 1 94640 57680 11236 8181 0.99559871 0 47 0
91000 6.489104 8187 0.23523152 37.267829 38.325313 6.4932786 14223 51 51 1 95718 57680 11236 8187 0.99090796 0 47 0
92000 6.3543839 8184 0.23502776 37.256587 38.292116 7.1553807 14265 66 66 1 96903 57680 11236 8184 0.98979354 0 47 0
93000 5.3818591 8195 0.23393122 37.236417 38.11346 7.9445503 14327 70 70 0 98138 57680 11236 8195 1.0030252 0 47 0
94000 5.6818206 8171 0.2343696 37.258958 38.184885 5.7311019 13943 63 63 0 99306 57680 11236 8171 0.99691591 0 47 0
95000 4.8295204 8160 0.23371941 37.291918 38.078951 4.5703653 13682 48 48 1 100491 57680 11236 8160 1.0078235 0 47 0
96000 6.7415726 8191 0.23539566 37.253429 38.352056 11.837168 13783 53 53 1 101761 57680 11236 8191 1.0010125 0 49 0
97000 6.5298825 8160 0.2354153 37.291126 38.355255 7.6012271 13688 65 65 0 103010 57680 11236 8160 0.99694671 0 49 0
98000 6.2804965 8175 0.23496877 37.259015 38.282504 5.4636503 14041 55 55 0 104259 57680 11236 8175 1.0042653 0 49 0
99000 5.7684681 8173 0.2344998 37.26605 38.206096 1.7807753 14082 58 59 0 105418 57680 11236 8173 0.99124491 0 49 0
100000 6.1371898 8170 0.23495644 37.28036 38.280495 7.7665232 13908 53 53 0 106633 57680 11236 8170 0.98281247 0 49 0
Loop time of 34.1892 on 8 procs for 100000 steps with 43990 atoms
Pair time (%) = 0.361362 (1.05695)
Neigh time (%) = 2.01143 (5.88324)
Comm time (%) = 2.27139 (6.64359)
Outpt time (%) = 0.00694269 (0.0203067)
Other time (%) = 29.5381 (86.3959)
Nlocal: 5498.75 ave 6478 max 4726 min
Histogram: 2 0 1 0 1 2 1 0 0 1
Nghost: 52.875 ave 59 max 44 min
Histogram: 1 0 1 0 1 0 1 0 2 2
Neighs: 186.25 ave 293 max 81 min
Histogram: 1 0 1 3 0 0 1 0 0 2
Total # of neighbors = 1490
Ave neighs/atom = 0.0338713
Neighbor list builds = 5000
Dangerous builds = 0
Please see the log.cite file for references relevant to this simulation

243
examples/ASPHERE/tri/log.tri.srd.1Feb14 Normal file
View File

@ -0,0 +1,243 @@
LAMMPS (1 Feb 2014)
# Aspherical shear demo - 3d triangle boxes, solvated by SRD particles
units lj
atom_style tri
atom_modify first big
read_data data.tri.srd
orthogonal box = (-8.43734 -8.43734 -8.43734) to (8.43734 8.43734 8.43734)
2 by 2 by 2 MPI processor grid
reading atoms ...
1500 atoms
1500 triangles
# add small particles as hi density lattice
lattice sc 0.4
Lattice spacing in x,y,z = 1.35721 1.35721 1.35721
region box block INF INF INF INF INF INF
lattice sc 20.0
Lattice spacing in x,y,z = 0.368403 0.368403 0.368403
create_atoms 2 region box
Created 91125 atoms
group big type 1
1500 atoms in group big
group small type 2
91125 atoms in group small
set group small mass 0.01
91125 settings made for mass
# delete overlaps
# must set 1-2 cutoff to non-zero value
pair_style lj/cut 1.5
pair_coeff 1 1 1.0 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0
delete_atoms overlap 1.5 small big
Deleted 76354 atoms, new total = 16271
# SRD run
reset_timestep 0
velocity small create 1.44 87287 loop geom
neighbor 0.3 bin
neigh_modify delay 0 every 1 check yes
neigh_modify exclude molecule big include big
communicate multi group big vel yes
neigh_modify include big
# no pairwise interactions with small particles
pair_style tri/lj 3.5
pair_coeff 1 1 0.1 1.0
pair_coeff 2 2 0.0 1.0 0.0
pair_coeff 1 2 0.0 1.0 0.0
# use fix SRD to push small particles out from inside big ones
# if comment out, big particles won't see SRD particles
timestep 0.001
fix 1 big rigid molecule #langevin 1.0 1.0 0.1 12398
125 rigid bodies with 1500 atoms
fix 2 small srd 20 big 1.0 0.25 49894 search 0.2 cubic warn 0.0001 shift yes 49829 overlap yes collision noslip
fix 3 all deform 1 x scale 0.6 y scale 0.6 z scale 0.6
# diagnostics
compute tsmall small temp/deform
compute tbig big temp
variable pebig equal pe*atoms/count(big)
variable ebig equal etotal*atoms/count(big)
compute 1 big erotate/asphere
compute 2 all ke
compute 3 all pe
variable toteng equal (c_1+c_2+c_3)/atoms
thermo 1000
thermo_style custom step f_1 c_tsmall f_2[9] temp press
thermo_modify temp tbig
WARNING: Temperature for thermo pressure is not for group all (../thermo.cpp:439)
compute 10 all property/atom corner1x corner1y corner1z corner2x corner2y corner2z corner3x corner3y corner3z
#dump 1 all custom 500 dump1.atom.srd id type x y z ix iy iz
#dump 2 all custom 500 dump1.tri.srd id type # c_10[1] c_10[2] c_10[3] c_10[4] c_10[5] c_10[6] # c_10[7] c_10[8] c_10[9]
run 10000
WARNING: Using fix srd with box deformation but no SRD thermostat (../fix_srd.cpp:385)
SRD info:
SRD/big particles = 14771 1500
big particle diameter max/min = 1.46014 0.436602
SRD temperature & lamda = 1 0.2
SRD max distance & max velocity = 0.8 40
SRD grid counts: 67 67 67
SRD grid size: request, actual (xyz) = 0.25, 0.251861 0.251861 0.251861
SRD per actual grid cell = 0.0647662
SRD viscosity = -1.09837
big/SRD mass density ratio = 24.668
WARNING: SRD bin size for fix srd differs from user request (../fix_srd.cpp:2853)
WARNING: Fix srd grid size > 1/4 of big particle diameter (../fix_srd.cpp:2875)
WARNING: Fix srd viscosity < 0.0 due to low SRD density (../fix_srd.cpp:2877)
# of rescaled SRD velocities = 0
ave/max small velocity = 19.9708 35.1504
ave/max big velocity = 0 0
WARNING: Using compute temp/deform with inconsistent fix deform remap option (../compute_temp_deform.cpp:76)
Memory usage per processor = 116.135 Mbytes
Step 1 tsmall 2[9] Temp Press
0 0 1.4405441 0 0 -0.15917996
1000 1.0535509 1.1241378 1.1224038 0 0.15526438
2000 2.4635987 1.0240667 0.94231519 0 0.011069846
3000 3.1093028 1.0070585 1.0176028 0 0.34124888
4000 2.9524179 1.0101774 1.0311733 0 -0.22118101
5000 2.9175438 1.0109144 1.0284445 0 0.1839507
6000 3.3200928 0.99108983 0.95968219 0 1.0229339
7000 3.3020156 0.99661428 0.98000621 0 2.3333851
8000 4.1500999 0.99129645 0.99723707 0 4.3054414
9000 6.537028 1.0074533 1.0251098 0 18.782913
10000 16.233245 1.0145766 1.0219787 0 147.83787
Loop time of 228.558 on 8 procs for 10000 steps with 16271 atoms
Pair time (%) = 145.018 (63.4492)
Neigh time (%) = 32.6283 (14.2757)
Comm time (%) = 43.3283 (18.9572)
Outpt time (%) = 0.000848889 (0.00037141)
Other time (%) = 7.58254 (3.31755)
Nlocal: 2033.88 ave 2092 max 1902 min
Histogram: 1 0 0 0 1 1 0 1 1 3
Nghost: 2805 ave 2855 max 2751 min
Histogram: 1 0 2 1 0 0 0 2 1 1
Neighs: 30026.5 ave 38700 max 24367 min
Histogram: 1 0 2 1 3 0 0 0 0 1
Total # of neighbors = 240212
Ave neighs/atom = 14.7632
Neighbor list builds = 501
Dangerous builds = 0
#undump 1
#undump 2
unfix 3
change_box all triclinic
triclinic box = (-5.0624 -5.0624 -5.0624) to (5.0624 5.0624 5.0624) with tilt (0 0 0)
fix 2 small srd 20 big 1.0 0.25 49894 search 0.2 cubic warn 0.0001 shift yes 49829 overlap yes collision noslip tstat yes
#dump 1 all custom 500 dump2.atom.srd id type x y z ix iy iz
#dump 2 all custom 500 dump2.tri.srd id type # c_10[1] c_10[2] c_10[3] c_10[4] c_10[5] c_10[6] # c_10[7] c_10[8] c_10[9]
fix 3 all deform 1 xy erate 0.05 units box remap v
run 40000
SRD info:
SRD/big particles = 14771 1500
big particle diameter max/min = 1.46014 0.436602
SRD temperature & lamda = 1 0.2
SRD max distance & max velocity = 0.8 40
SRD grid counts: 40 40 40
SRD grid size: request, actual (xyz) = 0.25, 0.25312 0.25312 0.25312
SRD per actual grid cell = -1.93929
SRD viscosity = -0.36972
big/SRD mass density ratio = -0.836253
WARNING: SRD bin size for fix srd differs from user request (../fix_srd.cpp:2853)
WARNING: Fix srd grid size > 1/4 of big particle diameter (../fix_srd.cpp:2875)
WARNING: Fix srd viscosity < 0.0 due to low SRD density (../fix_srd.cpp:2877)
# of rescaled SRD velocities = 1
ave/max small velocity = 16.0639 40
ave/max big velocity = 2.05735 6.73052
Memory usage per processor = 65.2471 Mbytes
Step 1 tsmall 2[9] Temp Press
10000 16.233245 1.0144194 0 0 146.27032
11000 17.301043 1.0619209 1 0 84.226323
12000 13.481177 1.063945 1 0 80.42834
13000 12.616615 1.0536125 1 0 68.365053
14000 12.16592 1.0581344 1 0 59.709941
15000 10.811557 1.0650453 1 0 59.869798
16000 9.9303081 1.0641012 1 0 59.289126
17000 10.452639 1.0589904 1 0 52.680235
18000 9.2488947 1.0556713 1 0 51.044123
19000 8.7179788 1.0566791 1 0 50.927924
20000 7.8641299 1.0668891 1 0 53.961014
21000 8.4333975 1.0628003 1 0 54.677068
22000 8.554623 1.0591872 1 0 51.784737
23000 7.895462 1.0598512 1 0 53.26885
24000 7.7014869 1.055569 1 0 54.236941
25000 7.508138 1.0581539 1 0 53.297709
26000 8.0707381 1.0586451 1 0 51.411341
27000 8.224223 1.0585509 1 0 52.775784
28000 8.8720461 1.0648837 1 0 49.681447
29000 8.503895 1.0538278 1 0 52.478556
30000 8.2272425 1.0625583 1 0 52.795709
31000 8.4626161 1.0599557 1 0 52.905343
32000 7.7410527 1.0587681 1 0 57.464472
33000 7.2765481 1.0677531 1 0 59.46728
34000 7.036016 1.0640728 1 0 60.636828
35000 8.497203 1.0577458 1 0 55.257312
36000 8.9218761 1.0585038 1 0 55.668201
37000 8.0469115 1.0561193 1 0 58.711565
38000 8.6456966 1.0663818 1 0 53.691658
39000 8.566667 1.0616928 1 0 52.866468
40000 7.7081259 1.0577046 1 0 54.105829
41000 7.4051952 1.0570277 1 0 56.094997
42000 7.331432 1.0632261 1 0 55.661067
43000 7.5587053 1.0602077 1 0 54.30762
44000 7.4521497 1.0594391 1 0 52.993548
45000 7.6234559 1.0641992 1 0 53.352202
46000 8.0659693 1.0609367 1 0 50.302682
47000 7.5227882 1.0563543 1 0 50.43551
48000 8.0798684 1.0628968 1 0 48.419201
49000 7.130784 1.0594143 1 0 50.656081
50000 7.3219268 1.0608236 1 0 54.808358
Loop time of 1651.84 on 8 procs for 40000 steps with 16271 atoms
Pair time (%) = 1232.76 (74.6296)
Neigh time (%) = 171.893 (10.4061)
Comm time (%) = 193.625 (11.7218)
Outpt time (%) = 0.00338289 (0.000204795)
Other time (%) = 53.5578 (3.24231)
Nlocal: 2033.88 ave 2229 max 1953 min
Histogram: 2 1 2 2 0 0 0 0 0 1
Nghost: 2757.88 ave 2810 max 2692 min
Histogram: 1 1 1 0 1 0 0 1 1 2
Neighs: 30098.6 ave 37369 max 25592 min
Histogram: 1 3 0 1 0 1 1 0 0 1
Total # of neighbors = 240789
Ave neighs/atom = 14.7987
Neighbor list builds = 2007
Dangerous builds = 0
Please see the log.cite file for references relevant to this simulation