应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 热导率计算问题
【调剂】北京石油化工学院2024年16个专业接受调剂
91/1返回列表
查看: 2910  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

qujinfeng

金虫 (小有名气)

[求助] 热导率计算问题已有1人参与

我是计算一个石墨烯纳米带300k下的热导率   下面是我的in文件   不知道为什么体系就是达不到平衡  整体温度一直在285到315之间震荡   局部温度震荡范围更大   请问我该怎么做   是什么问题   怎么样才能让体系弛豫达到平衡状态呢  求指教！  附件是计算的log文件  
boundary  s s s
units  metal
dimension  3
atom_style  atomic
atom_modify map array
neighbor 2.0 nsq
read_data   data.graphene1550
pair_style  airebo 3.0 1 0
pair_coeff * * CH.airebo  C
region  box  block    0.0   112.311    0.0  36.900    0.0  0.5   units box
region    begin      block   7.21        8.63      INF INF  INF  INF    units box
region    hot        block   8.63        17.152    INF INF  INF  INF    units box
region    a          block   17.152      25.674    INF INF  INF  INF    units box
region    b          block   25.674      34.195    INF INF  INF  INF    units box
region    c          block   34.195      42.717    INF INF  INF  INF    units box
region    d          block   42.717      51.239    INF INF  INF  INF    units box
region    e          block   51.239      59.76     INF INF  INF  INF    units box
region    e1         block   59.76       68.282    INF INF  INF  INF    units box
region    d1         block   68.282      76.804    INF INF  INF  INF    units box
region    c1         block   76.804      85.325    INF INF  INF  INF    units box
region    b1         block   85.325      93.847    INF INF  INF  INF    units box
region    a1         block   93.847      102.369   INF INF  INF  INF    units box
region    cold       block   102.369     110.891   INF INF  INF  INF    units box
region    end        block   110.891    112.311    INF INF  INF  INF    units box
region    boundary   union 2  begin end  units box
group   begin    region  begin   
group   hot      region  hot     
group   a        region  a      
group   b        region  b      
group   c        region  c      
group   d        region  d      
group   e        region  e      
group   e1       region  e1      
group   d1       region  d1      
group   c1       region  c1      
group   b1       region  b1      
group   a1       region  a1      
group   cold     region  cold   
group   end      region  end     
group   boundary  region boundary
group   1  subtract all boundary
velocity  1  create 300 7986875  rot yes dist gaussian units box
velocity begin  set 0.0  0.0  0.0  units box
velocity end    set 0.0  0.0  0.0  units box

fix 1  begin setforce 0.0 0.0 0.0
fix 2  end   setforce 0.0 0.0 0.0
fix a 1 nvt temp 300 300 0.025   
compute   ke 1 ke/atom                     
variable  temp  atom  c_ke*1.6e+4/(1.5*1.38)

fix       temp_profile    1   ave/spatial    1   100000  100000  x (0,0,0) 2.10202  v_temp  file  temp.profile  units  box
compute   c_temp   c   temp
compute   c1_temp  c1  temp
compute   1_temp   1   temp
variable  delta_temp   equal   c_c_temp-c_c1_temp                                         
fix       delta_out  1  ave/time   1  100000   100000   v_delta_temp  file  delta_temp.dat
thermo_style   custom   step c_1_temp  etotal  pe  ke   c_c_temp   c_c1_temp
#thermo_modify   lost  warn                                                  
timestep  0.0005                                                            
thermo    1000                                                               
dump   1  all  xyz  100  g.xyz         
unfix  a                                                                                                                  
fix 3  1   nve                                                                                                      
fix   hot   all  heat  1  1   region hot                                                                 
fix   cold  all  heat  1  -1  region  cold                                                            
variable   thermal_conductivity  equal   1*1.6e-19*40.448*1e-10/(36.9*3.4*1e-20*f_delta_out*1e-12)  
fix   thermal_conductivity_out  all   ave/time   100000  1   100000   v_thermal_conductivity  file thermal_conductivity.dat
                                                                                                                           
#run  5000000
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com
  • 附件 1 : log.lammps
    • 2012-06-28 20:29:51, 30.64 K

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-06-28 20:30:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

维特的日记

木虫 (小有名气)

【答案】应助回帖

我感觉是你nvt系综的damp值设置有问题
另外你计算热导率公式里的36.9*3.4*10e-20是那个区域的面积 有问题把
赞一下 回复此楼
2楼2014-08-20 22:23:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛犊小子

铁杆木虫 (著名写手)
帮顶…………
回复此楼
3楼2014-08-22 10:46:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛犊小子

铁杆木虫 (著名写手)
那个40.448是什么呀
赞一下 回复此楼
4楼2014-08-27 09:35:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tang_wei

金虫 (正式写手)
你好，请问一下你的问题解决了吗？我也是用fix heat 方法，发现即使我所加的热流为0的话，体系每一层的温度也不是想的的，波动较大
赞一下(1人) 回复此楼
5楼2015-11-20 09:24:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青春逍遥啊

新虫 (小有名气)
楼主您好，很抱歉我无法回答您的问题，我是刚刚入门做热导率的，我想问您variable  temp  atom  c_ke*1.6e+4/(1.5*1.38)是什么意思，我看手册没弄明白，如果您有时间，麻烦您可以教我。谢谢。
赞一下 回复此楼
6楼2015-12-02 19:57:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

计算化学

新虫 (初入文坛)
引用回帖:
5楼: Originally posted by tang_wei at 2015-11-20 09:24:56
你好，请问一下你的问题解决了吗？我也是用fix heat 方法，发现即使我所加的热流为0的话，体系每一层的温度也不是想的的，波动较大

我的也是，就是两端设置了冷热区能看到，但中间部分就像一条直线，且波动特大，请教一下，你知道原因吗
赞一下 回复此楼
7楼2018-05-24 09:38:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青春逍遥啊

新虫 (小有名气)
引用回帖:
7楼: Originally posted by 计算化学 at 2018-05-24 09:38:47
我的也是，就是两端设置了冷热区能看到，但中间部分就像一条直线，且波动特大，请教一下，你知道原因吗...

可以用其他的方法计算，比如mp之类的呢

发自小木虫Android客户端
赞一下 回复此楼
8楼2018-06-12 18:46:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

机智康

新虫 (初入文坛)
问题解决了吗？

发自小木虫Android客户端
回复此楼
9楼2020-05-07 13:35:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qujinfeng 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] SCI 评审意见回复 模板 15+3 Kilig0317 2024-04-18 4/200 2024-04-19 08:46 by bobvan
[基金申请] 国自然青年基金只有一篇能第一标注，其他都是第二标注能结题吗？ 10+4 lice_1987 2024-04-18 6/300 2024-04-19 08:29 by kudofaye
[基金申请] 院士建议：对勇于提出解决卡脖子问题新解决方案的学者加大资助力度 +4 zju2000 2024-04-16 7/350 2024-04-19 07:55 by newfuzzy1
[论文投稿] The total environment格式 25+4 lyyyyds0 2024-04-13 9/450 2024-04-18 09:33 by 我为人人研
[博后之家] 博后进站年龄可以超过35岁？ +8 suesong0818 2024-04-14 8/400 2024-04-18 08:38 by charles-c
[考研] 0854计算机技术316分调剂，求求导师 捞我一下，有学就上 +6 zhushijie218 2024-04-16 7/350 2024-04-17 23:33 by petro
[考研] 332求调剂 +3 木叶下1999 2024-04-16 5/250 2024-04-17 00:56 by angeliar
[基金申请] 请问教育部人文社科难度大吗 +9 锦衣卫寒战 2024-04-15 14/700 2024-04-16 23:10 by 锦衣卫寒战
[考研] 材料专硕329调剂遗留难民 +9 Kaylawander 2024-04-13 9/450 2024-04-15 19:21 by mumin1990
[考研] 328求调剂 +3 send rgsc 2024-04-14 7/350 2024-04-15 18:17 by zw_muchong
[论文投稿] Journal of The Science of Food and Agriculture投稿遇到格式问题 5+3 tongtongco 2024-04-12 3/150 2024-04-15 09:19 by shuoyi
[考研] 化学、材料类最后调剂机会！！！ +3 加油努力就好 2024-04-14 10/500 2024-04-15 09:05 by 任pen
[考研] 274求调剂 +5 四个大字，你滴� 2024-04-13 6/300 2024-04-15 00:01 by 四个大字，你滴�
[考研] 280求调剂 +4 黑皮冰棒 2024-04-12 7/350 2024-04-14 23:06 by 永字号
[考研] 338求调剂 +3 18280338551 2024-04-14 5/250 2024-04-14 10:03 by tcni
[考研] 一志愿中国科学院大学生物与医药专业，总分332，求调剂。 +4 2114020553 2024-04-13 4/200 2024-04-13 17:21 by lincunhui
[考研] 086000生医没考数学304分，一志愿西农求调剂 +3 王有田有学上 2024-04-13 5/250 2024-04-13 16:18 by coco1981
[论文投稿] material research express还没交开源版面费，已经线上出版了 10+3 lzjslyz 2024-04-12 3/150 2024-04-13 11:09 by kingzgh
[考研] 化工专硕求调剂 +9 青梅duoduo 2024-04-12 11/550 2024-04-13 09:33 by lincunhui
[硕博家园] 生物炭采购 +3 锦鲤附体@ 2024-04-12 3/150 2024-04-12 22:10 by Danny614
信息提示
关闭
请填处理意见
关闭