24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1487  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 一志愿华中农微生物,288分,三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 求调剂,一志愿材料科学与工程985,365分, +8 材化李可 2026-04-11 10/500 2026-04-12 08:42 by 852137818
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 一志愿厦大0856,306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 293求调剂 +8 勇远库爱314 2026-04-06 8/400 2026-04-11 20:25 by 蓝云思雨
[考研] 346,工科求调剂 +3 moser233 2026-04-09 3/150 2026-04-11 10:04 by zhq0425
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
[考研] 一志愿京区985,085401,与本科专业一致,电子信息工程, +4 阳光开朗的男孩 2026-04-10 4/200 2026-04-10 18:27 by shenrf
[考研] 一志愿沪9,326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 江苏大学 工科调剂 捡漏 +3 Evan_Liu 2026-04-09 5/250 2026-04-10 10:22 by Evan_Liu
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 332,085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文 +13 Emmy~ 2026-04-09 13/650 2026-04-09 14:47 by only周
[考研] 308求调剂 +17 墨墨漠 2026-04-06 17/850 2026-04-09 09:25 by 壹往無前
[考研] 材料与化工专硕306分找合适调剂 +27 沧海轻舟e 2026-04-06 28/1400 2026-04-08 22:06 by wdyheheeh
[考研] 求调剂 +15 熊二想上岸 2026-04-06 15/750 2026-04-08 04:53 by 无际的草原
[考研] 材料工程专业日语生求调剂 +9 111623 2026-04-07 9/450 2026-04-07 23:31 by 一只好果子?
信息提示
请填处理意见