| 查看: 758 | 回复: 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这是导热测定过程,所得结果温度波动较大,因为我没加热流,觉得温度应该是至少相等的 |
» 猜你喜欢
323求调剂
已经有6人回复
一志愿北京化工大学 070300 学硕 336分 求调剂
已经有4人回复
352求调剂
已经有3人回复
一志愿东华大学化学070300,求调剂
已经有8人回复
277材料科学与工程080500求调剂
已经有7人回复
317求调剂
已经有18人回复
293求调剂
已经有5人回复
280分求调剂 一志愿085802
已经有7人回复
0854电子信息求调剂
已经有3人回复
263求调剂
已经有4人回复













回复此楼