应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 专家会诊 » 关于poiseuille流
51/1返回列表
查看: 2047  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cuiliu1007

银虫 (小有名气)

[求助] 关于poiseuille流

请问，用分子动力学方法模拟高分子溶液的poiseuille流，体系中的外力如何加？如果模拟模型是一个长为L，直径为R的圆形管道，那外力是持续的加在整个管道中的粒子上，还是只加在端头的粒子上？是否需要等体系达平衡在收集数据？
期待您的回答！
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-05-09 16:05:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paladinlee

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★
jiaoyixiong: 金币+3, 鼓励交流 2013-05-15 08:25:55
引用回帖:
6楼: Originally posted by cuiliu1007 at 2013-05-14 17:45:31
在我的模型里，是将外力持续加在管道中的所有粒子上，外力沿x正方向的，x方向采用的是周期边界，y和z方向设置的是由墙粒子组成的截面为原型的管道，墙粒子与溶液粒子之间采用的是WCA模型。但是我得到的poiseuille曲 ...

WCA和LJ差不多，你的浸润性系数应该是默认为1吧？
“poiseuille曲线的幂指数就不是2.0，而是在1.5~3.0的范围内”，是波动吗？你算了多久？一般平衡态统计的结果还是要对时间平均的。
还有就是恒温的方法没问题，但是可能你的力加在所有流体粒子上会对流体控温影响，可以参考Soong, C. Y., T. H. Yen, et al. (2007). "Molecular dynamics simulation of nanochannel flows with effects of wall lattice-fluid interactions." Physical Review E 76(036303): 1-14.对部分流体粒子恒温。或者不加外力时NVT弛豫，加外力在NVE下模拟
赞一下 回复此楼
我是胖人，不是粗人
7楼2013-05-15 01:08:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★
感谢参与，应助指数 +1
cuiliu1007: 金币+1, ★有帮助 2013-05-13 16:23:04
去参考lammps    的example/flow  里的例子
赞一下 回复此楼
2楼2013-05-09 20:22:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paladinlee

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
cuiliu1007: 金币+5, ★★★很有帮助 2013-05-13 16:23:16
jiaoyixiong: 金币+2, 鼓励交流 2013-05-13 16:56:12
两种方法都可以，前一种（外力是持续的加在整个管道中的粒子上）是整个管道流向压力恒定，后一种是有压力差。具体看你想要模拟什么结果。
数据需要体系达到稳态平衡后再统计。
赞一下 回复此楼
我是胖人，不是粗人
3楼2013-05-09 22:45:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuiliu1007

银虫 (小有名气)
引用回帖:
3楼: Originally posted by paladinlee at 2013-05-09 22:45:17
两种方法都可以，前一种（外力是持续的加在整个管道中的粒子上）是整个管道流向压力恒定，后一种是有压力差。具体看你想要模拟什么结果。
数据需要体系达到稳态平衡后再统计。

请问，采用将外力持续加在整个管道中的粒子上的方式后，边界应该怎么处理呢？我的边界是放入由粒子组成的粒子墙，采用边界粒子和反弹边界联用的方式。可是我得到的poiseuille曲线的幂指数在1.5到3.0的范围内，幂指数的大小依赖外力的大小，外力大幂指数高。这是不是说我体系中的处理是不合适的，结果不对？我应该怎么处理才对呢？
赞一下 回复此楼
4楼2013-05-13 16:27:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭