24小时热门版块排行榜

返回列表

200909140505

铜虫 (小有名气)

应助: 3 (幼儿园)
金币: 30.7
散金: 11
红花: 2
帖子: 75
在线: 52.9小时
虫号: 1175337
注册: 2010-12-22
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] fix heat comand计算晶格热导

最近在学习用nemd计算晶格热导，运行了lammps官网上关于计算Ar晶格热导的例子，以下是log.lammps文件。
Step Temp Jx Jy Jz k11 k22 k33
   0 72.521833 2.498502e-06 -2.0047631e-07 -6.7129653e-07 0.0088667428 5.708616e-05 0.00064007888
2000    66.11113 1.4391576e-06 1.6829237e-06 -1.4599868e-06 0.98583388 1.3459874 0.20604396
4000 68.656729 1.5105571e-06 1.9964751e-06 -2.6599568e-06 0.24943423 0.48354981 1.6447711
6000    65.51534 -2.6166856e-07 1.4559872e-06 7.0150247e-07 0.16570196 0.58961854 2.3785217
8000 68.351915 -1.3303494e-06 4.3049495e-06 2.2655487e-06    0.130139 0.52400031 1.6302988
10000 68.442156 -3.7559096e-06 -1.2824083e-06 2.6004241e-08 0.49872705 1.0163064 1.532751
（1）后六列分别表示三个方向的热流和热导率吧？
（2）另外一个输出文件内容如下，这个就不太看得懂了，求解释。
# Time-correlated data for fix JJ
# Timestep Number-of-time-windows
# Index TimeDelta Ncount c_flux[1]*c_flux[1] c_flux[2]*c_flux[2] c_flux[3]*c_flux[3]
2000 200
1 0 201 0.000485669 0.000744015 0.000526979
2 10 200 0.000454311 0.000712065 0.000493873
3 20 199 0.000387833 0.000648309 0.00041936
4 30 198 0.00033294 0.000606816 0.000353151
5 40 197 0.000302059 0.000591577 0.000320354
6 50 196 0.000277342 0.000570019 0.00031651
7 60 195 0.000243456 0.000525581 0.000324793
8 70 194 0.000203369 0.000473704 0.000328459
9 80 193 0.000170657 0.000432944 0.000315548
10 90 192 0.000151629 0.000403438 0.000286584
11 100 191 0.000139125 0.000378841 0.000257655
12 110 190 0.000125903 0.000357865 0.000244263
13 120 189 0.000113193 0.000342032 0.000247021
14 130 188 0.000106422 0.000326952 0.000257794
15 140 187 0.00010766 0.000304277 0.000267782
16 150 186 0.000112893 0.000272026 0.000268346
17 160 185 0.000116312 0.000237737 0.000254589
18 170 184 0.000116059 0.000211744 0.000230521
19 180 183 0.000113949 0.000196721 0.00020776

回复此楼