计算热导率的时候,为什么体系温度一直特别高?
各位大侠,请问一下我在计算二氧化铪hfo2热导率的时候,采用的常热流密度法,为什么最后统计的温度特别高,吓死人。
以下是In文件,请各位大侠指点迷津。
图片横坐标沿厚度方向取的点,纵坐标是温度。为什么会那么高?
units metal
boundary p p p
dimension 3
atom_style charge
neighbor 2.0 nsq
read_data hfo2-2-3-10.data
group Hf type 1
group O type 2
pair_style comb
pair_coeff * * ffield.comb Hf O
timestep 0.0004
min_style cg
minimize 1.0e-12 1.0e-12 1000 1000
reset_timestep 0
variable temp equal temp
variable ke equal ke
variable pe equal pe
variable press equal press
variable vol equal vol
variable etotal equal etotal
compute tot_temp all temp
compute myKE all ke/atom
variable temp0 atom 16000*c_myKE/(1.5*1.38)
fix info all ave/time 1 10000 10000 c_tot_temp v_etotal file cryTinfo-2-3-10.data
dump 1 all custom 5000 cryTmovie2-6-14.lammpstrj id element x y z
thermo 1000
velocity all create 150 56789 rot yes mom yes dist gaussian units box
fix NVE all nve
fix ALL all langevin 300 300 10 1983
run 1000000
fix heat_swap all thermal/conductivity 100 z 100
fix energy all ave/time 100 100 10000 f_heat_swap file cryTenergy2-3-10.data
compute c1 all chunk/atom bin/1d z lower 1 units box
fix 1 all ave/chunk 1 1000000 1000000 c1 v_temp0 file cryTfile2-3-10.data
run 1000000
%YB1RS3884$G3NRT`S@7YTX.png 返回小木虫查看更多
,