24小时热门版块排行榜

返回列表

congmingbing

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 2089.9
散金: 106
帖子: 333
在线: 412.3小时
虫号: 446335
注册: 2007-10-29
专业: 凝聚态物性I:结构、力学和

[交流] 【求助】lammps计算heat/flux出错已有1人参与

在lammps手册上看到了compute heat/flux command，用提供的输入文件执行，出错了，错误信息是
ERROR: Illegal compute heat/flux command，不知道这个错误怎么解决，感觉像是版本的问题。我用的lammps是在主页上下载的lammps-15Jan10.tar.gz，这个应该是最新版了吧。
刚开始学lammps,还请高手指点！万分感谢。
下面是in文件。
＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃

atom_style    atomic
units real
dimension 3
boundary p p p
lattice fcc  5.376  orient x 1 0 0 orient y 0 1 0 orient z 0 0 1
region    box block 0 4 0 4 0 4
create_box 1 box
create_atoms 1 box
mass 1 39.948
pair_style lj/cut 13.0
pair_coeff * * 0.2381 3.405
group every region box
velocity all create 70 102486 mom yes rot yes dist gaussian
timestep 4.0
thermo          10

# ------------- Equilibration and thermalisation ----------------

fix NPT all npt 70 70 10 xyz 0.0 0.0 100.0 drag 0.2
run 8000
unfix          NPT

# --------------- Equilibration in nve -----------------

fix NVE all nve
run 8000

# -------------- Flux calculation in nve ---------------

reset_timestep  0
compute          myKE all ke/atom
compute          myPE all pe/atom
compute          myStress all stress/atom virial
compute flux all heat/flux myKE myPE myStress
log      flux.log
variable       J equal c_flux[1]/vol
thermo_style custom step temp v_J
run          100000

回复此楼