24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1534  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +3 光伏爱好者 2026-05-17 4/200 2026-05-19 01:19 by 风起于浮萍
[基金申请] 评审有感 +6 popular289 2026-05-18 11/550 2026-05-18 22:40 by iaeyuan
[基金申请] 今年审到国自然15份,谈谈感受 +16 国自然国社科中 2026-05-17 16/800 2026-05-18 14:58 by gy116024
[考博] 找博士生导师 +4 小代想上岸 2026-05-15 5/250 2026-05-17 17:44 by 脱颖而出
[考博] 2026博士还有哪些学校有名额 +6 小王求读研 2026-05-15 7/350 2026-05-17 16:54 by 知音湖畔
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +13 dd921ww 2026-05-12 16/800 2026-05-17 12:38 by 于轩
[考博] 26/27博士推荐 +3 1木头人13949 2026-05-13 3/150 2026-05-17 09:41 by YuY66
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:06 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:25 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:45 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-16 17:45 by j6b2pdz07o
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
请填处理意见