| 查看: 814 | 回复: 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这是导热测定过程，所得结果温度波动较大，因为我没加热流，觉得温度应该是至少相等的 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
20