当前位置: 首页 > 分子模拟 >计算热导率的时候,为什么体系温度一直特别高?

计算热导率的时候,为什么体系温度一直特别高?

作者 举重若轻I
来源: 小木虫 1050 21 举报帖子
+关注

各位大侠,请问一下我在计算二氧化铪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 返回小木虫查看更多

今日热帖
  • 精华评论
猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓