| 查看: 817 | 回复: 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这是导热测定过程，所得结果温度波动较大，因为我没加热流，觉得温度应该是至少相等的 |
回复此楼» 猜你喜欢
申博求助
已经有4人回复
-
27年申博自荐-能源/环境方向
已经有3人回复
-
闲聊
已经有3人回复
-
亲密关系
已经有13人回复
-
27秋申博自荐
已经有4人回复