24小时热门版块排行榜    

查看: 740  |  回复: 0

tang_wei

金虫 (正式写手)

[求助] 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这是导热测定过程,所得结果温度波动较大,因为我没加热流,觉得温度应该是至少相等的
回复此楼

» 猜你喜欢

已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tang_wei 的主题更新
信息提示
请填处理意见