| 查看: 4136 | 回复: 14 | ||
leigp木虫 (正式写手)
|
[求助]
fix heat计算石墨烯热导率，固定两端原子是否合适？ 已有1人参与
|
|
|
大家好，我刚接触lammps，想采用fix heat方法计算石墨烯纳米带的热导率。我看很多文献上为了防止石墨烯纳米带的整体移动通常将纳米带两端的原子固定。但是我采用该方法进行计算后，发现在弛豫过程中，虽然体系整体的温度保持在400K，但是由于两端原子固定的原因导致温度分布出现中间高两端低的现象，如下图： 当体系达到平衡后，在体系的左端加入热量，右端抽取热量。由于固定原子的存在导致与其相邻的原子速度相对较低，因此如果抽取的热量过高时则会导致原子速度降低至零而出现错误。下面是我的in文件，希望了解的大神能帮我看一下，是否是我in文件中设置有问题，在此谢过了！！ # MD simulation of Graphene thermal conductivity # Initialization units metal variable T equal 400 variable V equal vol variable dt equal 0.0005 variable p equal 10000 # correlation length variable s equal 1 # sample interval variable d equal $p*$s # dump interval variable width equal 14.6544 variable kB equal 1.3806504e-23 # [J/K] Boltzmann variable eV2J equal 1.602e-19 variable A2m equal 1.0e-10 variable ps2s equal 1.0e-12 variable convert equal ${eV2J}*${eV2J}/${ps2s}/${A2m} dimension 3 newton on boundary p p p atom_style full neighbor 0.3 bin neigh_modify check yes lattice fcc 1 region box block 0.000 14.6544 0.000 100.00 0.000 50.000 units lattice y方向为热流施加方向，为防止两端原子的相互作用将y方向盒子长度设置为10nm,并将石墨烯纳米带放置于盒子中部 read_data graphene.data ################### set group ############################ region up block INF INF INF 26.0 INF INF units lattice region down block INF INF 81 INF INF INF units lattice region steady union 2 up down group steady region steady region remain block INF INF 26.0 81 INF INF units lattice group remain region remain region hot block INF INF 26.0 36.0 INF INF units lattice group hot region hot region cold block INF INF 71 81 INF INF units lattice group cold region cold velocity remain create $T 458 mom yes dist gaussian units lattice velocity steady set 0.0 0.0 0.0 # Potential ********************************************************* pair_style airebo 2.0 pair_coeff * * CH.airebo C fix nve remain nve fix temp remain temp/berendsen $T $T 0.0005 fix 1 steady setforce 0.0 0.0 0.0 compute mytemp remain temp compute hottemp hot temp compute coldtemp cold temp compute ke remain ke/atom variable temp atom c_ke*${eV2J}/(1.5*${kB}) fix temp_profile1 remain ave/spatial 1 1000000 1000000 y lower 1.0 v_temp file 11.profile units lattice thermo_style custom step c_mytemp c_hottemp c_coldtemp etotal vol thermo_modify lost warn thermo 500 timestep ${dt} dump 1 all atom 1000 dump.lammps run 2000000 unfix temp fix hot remain heat 1 1 region hot fix cold remain heat 1 -1 region cold run 2000000 reset_timestep 0 fix temp_profile2 remain ave/spatial 1 2000000 2000000 y lower 1.0 v_temp file temp.profile units lattice dump 2 all atom 1000 graphene.lammps # Run run 10000000 # every 1 set group steady vx 0 vy 0 vz 0  未施加温度梯度时的温度分布 |
回复此楼» 猜你喜欢
271求调剂
已经有30人回复
-
材料类284调剂
已经有50人回复
-
291分调剂
已经有9人回复
-
调剂求收留
已经有34人回复
-
291 求调剂
已经有38人回复
-
22408 312求调剂
已经有17人回复
-
一志愿华中农业071010，320求调剂
已经有6人回复
-
290调剂生物0860
已经有41人回复
-
291求调剂
已经有3人回复
-
211本科材料化工求调剂
已经有23人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 石墨烯的德拜温度
已经有8人回复
- 求助：lammps非周期边界
已经有22人回复
9楼2019-11-20 09:17:56
w76990480
荣誉版主 (文坛精英)
super-emuch
- 应助: 47 (小学生)
- 贵宾: 13.729
- 金币: 109757
- 散金: 29386
- 红花: 337
- 沙发: 50
- 帖子: 23155
- 在线: 3166.8小时
- 虫号: 942733
- 注册: 2010-01-15
- 专业: 固体力学
- 管辖: 文献求助
2楼2014-12-02 10:46:22
leigp
木虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 1811.8
- 散金: 115
- 红花: 5
- 帖子: 638
- 在线: 225小时
- 虫号: 758557
- 注册: 2009-04-27
- 性别: GG
- 专业: 分离过程
3楼2014-12-02 13:38:27
lihuan12593
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 26.8
- 散金: 180
- 红花: 4
- 帖子: 98
- 在线: 58.2小时
- 虫号: 1785524
- 注册: 2012-04-29
- 性别: GG
- 专业: 核技术及其应用
4楼2015-11-07 11:12:42
对马拉二胡
至尊木虫 (职业作家)
- 应助: 27 (小学生)
- 金币: 26097.2
- 散金: 244
- 红花: 13
- 帖子: 4095
- 在线: 939.3小时
- 虫号: 2471734
- 注册: 2013-05-19
- 性别: GG
- 专业: 无机非金属材料
|
层主好，最近一直被fix heat 命令困扰着。 用fix heat 命令 fix 1 hot heat 1 4.3 region hot fix 2 cold heat 1 -4.3 region cold ERROR: Fix heat kinetic energy went negative (../fix_heat.cpp:201) Last command: run 100000 然后我改变了fix 2 中抽出热流值（-4.3<efulx<-2.1）还是出现以上错误 如果eflux设置>-2.1，体系温度又上升太高，请问这该怎么办呢？ 还有输出和抽出的热流值是否一定要一致呢？ |
5楼2017-06-13 16:51:44
对马拉二胡
至尊木虫 (职业作家)
- 应助: 27 (小学生)
- 金币: 26097.2
- 散金: 244
- 红花: 13
- 帖子: 4095
- 在线: 939.3小时
- 虫号: 2471734
- 注册: 2013-05-19
- 性别: GG
- 专业: 无机非金属材料
|
楼主好，最近一直被fix heat 命令困扰着。 用fix heat 命令 fix 1 hot heat 1 4.3 region hot fix 2 cold heat 1 -4.3 region cold ERROR: Fix heat kinetic energy went negative (../fix_heat.cpp:201) Last command: run 100000 然后我改变了fix 2 中抽出热流值（-4.3<efulx<-2.1）还是出现以上错误 如果eflux设置>-2.1，体系温度又上升太高，请问这该怎么办呢？ 还有输出和抽出的热流值是否一定要一致呢？ |
» 本帖已获得的红花（最新10朵）
wHy~!~sm6楼2017-06-13 16:51:58
7楼2018-12-12 19:29:10
8楼2019-04-09 21:47:02
10楼2019-11-29 09:08:25
