24小时热门版块排行榜    

查看: 1600  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金中了 +8 laoda193707 2026-08-06 8/400 2026-08-07 07:37 by emeraldine
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ANUsWpSnuWQD 2026-08-06 4/200 2026-08-07 04:46 by nSM3ys5fefM3
[基金申请] 听说今天filecode变了 +19 布布和一二 2026-08-06 34/1700 2026-08-07 03:13 by 虫友是什么虫
[基金申请] filecode +10 等待解的谜 2026-08-06 14/700 2026-08-06 21:16 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +3 archvillain 2026-08-06 7/350 2026-08-06 21:05 by archvillain
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 6/300 2026-08-06 20:19 by 苏知砚
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[教师之家] 咨询面上基金 +4 李长云 2026-07-31 7/350 2026-08-06 11:20 by 李长云
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +22 低垂的野花 2026-07-31 30/1500 2026-08-05 18:03 by 低垂的野花
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[论文投稿] 十年后又回来了,论文投稿求助 +3 哈哈114477 2026-08-01 3/150 2026-08-04 15:40 by tegsgjy20
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见