| 查看: 1039 | 回复: 3 | ||
诗意栖居木虫 (小有名气)
|
[求助]
狭缝中汽液相变的分子动力学模拟
|
用fortran语言编程，壁面是铂原子，液相是氩原子，模拟狭缝中的汽液相变问题，氩原子之间是LJ力，壁面原子（铂原子）与氩的力也是LJ势能，壁面是沿X方向的，在Y，Z方向上是周期边界，氩原子竟然能跑到壁面外面，请求指导 这是模拟的模型 [ Last edited by ghcacj on 2011-11-23 at 10:38 ] |
回复此楼» 猜你喜欢
293求调剂
已经有12人回复
-
070300化学319求调剂
已经有4人回复
-
268求调剂
已经有9人回复
-
301求调剂
已经有4人回复
-
0703化学336分求调剂
已经有5人回复
-
301求调剂
已经有4人回复
-
材料专硕306英一数二
已经有9人回复
-
304求调剂
已经有8人回复
-
材料工程专硕调剂
已经有5人回复
-
344求调剂
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 相变动力学（Kinetics of phase transitions）
已经有256人回复
- 求助有关材料相变动力学研究的问题
已经有6人回复
- DPD这种耗散粒子动力学应用于蛋白质自组装研究
已经有3人回复
- 【原创】分子对接，同源模建，分子动力学软件介绍帖
已经有35人回复
- 【求助】相变模拟
已经有6人回复
2楼2011-11-23 12:22:55
【答案】应助回帖
★ ★ ★
诗意栖居(金币+3): 在程序运行二十部左右，有十几个粒子跑出来，但以后的过程中飞出的粒子数就不再增加了，如果像你说的那样，是L-J势不够，应该会越来越多的粒子跑出来吧 2011-11-24 09:00:57
ghcacj(金币+3): 谢谢 2011-11-25 10:39:31
诗意栖居(金币+3): 在程序运行二十部左右，有十几个粒子跑出来，但以后的过程中飞出的粒子数就不再增加了，如果像你说的那样，是L-J势不够，应该会越来越多的粒子跑出来吧 2011-11-24 09:00:57
ghcacj(金币+3): 谢谢 2011-11-25 10:39:31
|
LZ： 我的分析如下： 一：如果氩原子和壁面只是LJ相互作用的话是不够的，实际上还应考虑碰撞相互作用（效果就是反弹氩原子）。因为LJ势相对来说偏软（作用比较弱），故可能使氩原子穿透壁面。 可能解决方案： 一：简单处理的话把可以加入一人为的弹性散射程序，当氩原子和壁面接触时，作弹性散射考虑。 二：复杂的解决方案：短程相互作用引入范德瓦尔斯相互作用。 看情况而定吧  谢谢 |
3楼2011-11-23 19:41:31
4楼2011-11-24 13:04:46
