24小时热门版块排行榜

返回列表

tang_wei

金虫 (正式写手)

应助: 24 (小学生)
金币: 1354.3
红花: 28
帖子: 428
在线: 81.7小时
虫号: 4067192
注册: 2015-09-12
专业: 高分子物理与高分子物理化

[求助] fix heat 导热率测定的问题

我是用fix heat 方法，发现即使我所加的热流为0的话，体系每一层的温度也不是想的的，波动较大
velocity          all scale ${ct}
unfix             1
fix                1 all  nve
#fix             add    all heat 1  0.03  region hot
#fix             subtract all heat 1 -0.03  region cold

compute          Thot  all temp/region hot
compute          Tcold all temp/region cold

variable          Tdelta equal c_Thot-c_Tcold

thermo_style    custom step temp c_Thot c_Tcold v_Tdelta
thermo             1000
run             200000

compute          ke all ke/atom
variable          temp atom c_ke/(1.5*1.986857e-3)

fix             3 all ave/spatial 10 100 1000 y lower 0.02 v_temp file tmp.profile units reduced
fix             ave all ave/time 10 100 1000 v_Tdelta file Tdelta.profile
thermo_style    custom step temp c_Thot c_Tcold v_Tdelta f_ave

run             100000这是导热测定过程，所得结果温度波动较大，因为我没加热流，觉得温度应该是至少相等的

回复此楼

» 猜你喜欢

1楼 2015-11-20 09:30:42

已阅关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 tang_wei 的主题更新

返回列表