24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1486  |  回复: 3

wang06x

木虫 (小有名气)

[求助] 热传导边界问题

最近看论坛上一些热传导例子,帮助很大。首先感谢前辈们的无私奉献精神。
同时我有些问题很困惑,希望大家给予解答。十分感激!
1.我发现例子中3个方向都是周期性边界条件,那么请问G-K公式、MP法 以及fix heat法这三种方法可以使用非周期边界条件么?如 f p p
2.我自己尝试了把条件改为f p p可是  出现了丢失原子(lost atoms)的情况,请问这是为什么?看相关帖子说可能是有覆盖的原子,应该用minimize 命令,可是我用了之后,虽然不像原来的那样一开始就出现lost atoms ,但是跑了几步之后仍然出现了。
3.还有的人说可以用fix  wall/reflect 命令,我分别尝试了命令:fix 1 all wall/refiect xlo -150  xhi 150 units box(盒子x方向尺寸是-100-100.单位是metal) 或者fix 2 all wall/reflect  xlo EDGE xhi EDGE 。可是还是会出现lost atoms.是我的命令用错了么?
4.论坛上例子中NEMD法貌似都是先给定热流密度的,那么请问有没有使用分子动力学模拟给定温度梯度的热传导方面相关理论文献啊?
问题有点多。水平有限,真的想不出来。希望大家给予解答,十分感谢!
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐,对您同样有帮助:

低调
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

咫尺天涯1991

新虫 (初入文坛)

楼主解决了吗?我也遇到原子丢失的情况了
2楼2016-07-15 09:51:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

第一,用自由边界条件是可能丢失原子的,这就类似于初中学过的升华现象。所以,用nemd时最好用固定边界条件,即固定两端一小部分原子。用nemd方法时也可以用周期边界条件,此时热流会从热源左右两端流出,从热汇左右两端流入。即使用了这两种边界条件,系统在高温下还是有可能不稳定,因为这是非平衡模拟,就是不那么稳定。第二,可以先用速度标度法或动量交换法产生热流,再测量温度梯度,也可以用热浴,如郎之万热浴或者nose-hoover chain热浴控制热源和热汇的温度。它们都大同小异。第三,如果用emd方法,即格林-久保公式,在输运方向必须用周期边界条件,在其它方向的边界条件视模拟体系而定。以上所说都是热导率模拟的最基本知识,文献有几百篇,其中重要的有几十篇,楼主可自行查阅,不可看看lammps的例子就拉到。我看过lammps作者关于热导率模拟的一个ppt,感觉他并不是很懂。

[ 发自手机版 http://muchong.com/3g ]
3楼2016-07-15 18:30:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

第一,用自由边界条件是可能丢失原子的,这就类似于初中学过的升华现象。所以,用nemd时最好用固定边界条件,即固定两端一小部分原子。用nemd方法时也可以用周期边界条件,此时热流会从热源左右两端流出,从热汇左右两端流入。即使用了这两种边界条件,系统在高温下还是有可能不稳定,因为这是非平衡模拟,就是不那么稳定。第二,可以先用速度标度法或动量交换法产生热流,再测量温度梯度,也可以用热浴,如郎之万热浴或者nose-hoover chain热浴控制热源和热汇的温度。它们都大同小异。第三,如果用emd方法,即格林-久保公式,在输运方向必须用周期边界条件,在其它方向的边界条件视模拟体系而定。以上所说都是热导率模拟的最基本知识,文献有几百篇,其中重要的有几十篇,楼主可自行查阅,不可看看lammps的例子就拉到。我看过lammps作者关于热导率模拟的一个ppt,感觉他并不是很懂。

[ 发自手机版 http://muchong.com/3g ]
4楼2016-07-15 18:31:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wang06x 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854调剂 +9 长弓傲 2026-04-12 11/550 2026-04-12 22:13 by paopaotu326
[考研] 电子信息270求调剂 +18 terminal469 2026-04-07 18/900 2026-04-12 16:23 by ajpv风雷
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 求调剂,一志愿材料科学与工程985,365分, +8 材化李可 2026-04-11 10/500 2026-04-12 08:42 by 852137818
[考研] 求调剂 +18 MAX怅惘 2026-04-09 20/1000 2026-04-11 23:31 by zhen~
[考研] 271求调剂 +20 2261744733 2026-04-11 22/1100 2026-04-11 23:14 by labixiaoqiao
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 070300化学279求调剂 +19 哈哈哈^_^ 2026-04-08 20/1000 2026-04-11 20:43 by stoner78
[考研] 调剂 +6 青灯不负 2026-04-09 6/300 2026-04-11 20:35 by dongdian1
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 一志愿华东师范生物学326分,求调剂 +8 刘墨墨 2026-04-09 8/400 2026-04-10 12:00 by pengliang8036
[考研] 0702物理学学硕299求调剂 +6 祁柒连 2026-04-06 6/300 2026-04-10 11:10 by Roomoo
[考研] 江苏大学 工科调剂 捡漏 +3 Evan_Liu 2026-04-09 5/250 2026-04-10 10:22 by Evan_Liu
[考研] 材料调剂 +5 hzhahg 2026-04-06 5/250 2026-04-10 10:10 by may_新宇
[考研] 初试分332,一志愿报考西北工业大学, +11 故人?? 2026-04-09 11/550 2026-04-09 21:54 by JineShine
[论文投稿] 求助文献原文 10+3 18500821399 2026-04-08 3/150 2026-04-09 16:56 by 北京莱茵润色
[考研] 283电子信息求调剂 +4 三石WL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
[考研] 材料调剂 +13 汉123456 2026-04-07 14/700 2026-04-07 22:53 by 来看流星雨10
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
信息提示
请填处理意见