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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 专家会诊 » 关于poiseuille流
131/2返回列表12下一页
查看: 2177  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

cuiliu1007

银虫 (小有名气)

[求助] 关于poiseuille流

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

» 猜你喜欢

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

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

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的回帖

paladinlee

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
cuiliu1007: 金币+5, ★★★很有帮助 2013-05-14 17:46:56
引用回帖:
4楼: Originally posted by cuiliu1007 at 2013-05-13 16:27:48
请问，采用将外力持续加在整个管道中的粒子上的方式后，边界应该怎么处理呢？我的边界是放入由粒子组成的粒子墙，采用边界粒子和反弹边界联用的方式。可是我得到的poiseuille曲线的幂指数在1.5到3.0的范围内，幂指 ...

一般外力持续加在管道所有粒子上后，由于没有压力差，管道两边都采用周期边界条件即可。流体与壁面的采用碰撞反弹模型或简单地LJ作用模型都可以。
赞一下 回复此楼
我是胖人，不是粗人
5楼2013-05-14 02:35:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuiliu1007

银虫 (小有名气)
引用回帖:
5楼: Originally posted by paladinlee at 2013-05-14 02:35:15
一般外力持续加在管道所有粒子上后，由于没有压力差，管道两边都采用周期边界条件即可。流体与壁面的采用碰撞反弹模型或简单地LJ作用模型都可以。...

在我的模型里，是将外力持续加在管道中的所有粒子上，外力沿x正方向的，x方向采用的是周期边界，y和z方向设置的是由墙粒子组成的截面为原型的管道，墙粒子与溶液粒子之间采用的是WCA模型。但是我得到的poiseuille曲线的幂指数就不是2.0，而是在1.5~3.0的范围内。我不知道我的模型是哪里出了问题？体系中的墙粒子和溶液粒子采用nose-hoover方法分别控温，是对peculiar速度控温的。
赞一下 回复此楼
6楼2013-05-14 17:45:31
已阅   回复此楼   关注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的回帖

cuiliu1007

银虫 (小有名气)
引用回帖:
7楼: Originally posted by paladinlee at 2013-05-15 01:08:23
WCA和LJ差不多，你的浸润性系数应该是默认为1吧？
“poiseuille曲线的幂指数就不是2.0，而是在1.5~3.0的范围内”，是波动吗？你算了多久？一般平衡态统计的结果还是要对时间平均的。
还有就是恒温的方法没问题， ...

WCA的作用参数是设为1.0的。我的结果中速度曲线的幂指数不是在1.5~3.0之间波动，基本是随外力的增加单调上升的，就是说外力小时1.5，随外力增大幂指数会增加到3.0.我的体系都是在NVT系棕下计算的，先让它在平衡MD下（即不加外力场）松弛100万步，之后再非平衡MD下松弛100万步（加外力场），之后在随后的600万步里取数据，所得的速度曲线是几千个时间点的数据的统计平均。
不太明白“力加在所有流体粒子上对流体控温有影响”是指什么？我在控温的时候已经去除了流体的stream velocity。
另外，“加外力在NVE下模拟“是不是就是不控温？
赞一下 回复此楼
8楼2013-05-15 18:06:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paladinlee

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★
cuiliu1007: 金币+9, ★★★★★最佳答案 2013-05-16 19:31:46
引用回帖:
8楼: Originally posted by cuiliu1007 at 2013-05-15 18:06:10
WCA的作用参数是设为1.0的。我的结果中速度曲线的幂指数不是在1.5~3.0之间波动，基本是随外力的增加单调上升的，就是说外力小时1.5，随外力增大幂指数会增加到3.0.我的体系都是在NVT系棕下计算的，先让它在平衡MD下 ...

一般管内流体的控温流程是，前期不受外力的整个体系NVT平衡后。加上外力，NVT计算。这里一般对与管道接触的部分流体进行控温，也就是流体与壁面热交换。实际流体的温度也是通过这种热交换存在的。像我建议你参考的那文献，就可以归属到这一类里面（在壁面另一方向设置虚拟的恒温区）。
还有就是NVE模拟了，前期NVT稳定后，体系温度稳定。在NVE下模拟，这时控制的是volume 和energy，temperature不需要控制。
赞一下 回复此楼
我是胖人，不是粗人
9楼2013-05-16 00:18:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuiliu1007

银虫 (小有名气)
引用回帖:
9楼: Originally posted by paladinlee at 2013-05-16 00:18:52
一般管内流体的控温流程是，前期不受外力的整个体系NVT平衡后。加上外力，NVT计算。这里一般对与管道接触的部分流体进行控温，也就是流体与壁面热交换。实际流体的温度也是通过这种热交换存在的。像我建议你参考的 ...

好的，我试试NVE的方式。还有两个问题：
1.我是不是可以将“对与管道接触的部分流体进行控温”调节一下，使墙粒子的厚度增加，即多几层墙粒子，之后我直接对墙粒子控温而不控与管道接触的部分流体呢？因为对流体控温就涉及到peculiar速度的计算。
2.“一般管内流体的控温流程是，前期不受外力的整个体系NVT平衡后。加上外力，NVT计算。”这段话是不是说不加外力之前NVT，加外力之后，与管道接触的部分流体NVT，其它流体NVE？

还有，你推荐的文献我已经找到了，正在研读，谢谢！
非常感谢你的耐心解答!
赞一下 回复此楼
10楼2013-05-16 19:31:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cuiliu1007 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22408 366分，本科211，一志愿西工大 +3 Rubt 2026-04-09 3/150 2026-04-10 09:38 by liuhuiying09
[考研] 化工调剂求导师收留！一志愿失利，踏实肯干，有植物提取科研经历 +15 yzyzx 2026-04-09 15/750 2026-04-10 09:36 by grey499
[考研] 求调剂 +5 archer.. 2026-04-09 7/350 2026-04-09 22:18 by lbsjt
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH（等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 计算机11408，286分求调剂 +9 木子念晞 2026-04-05 9/450 2026-04-09 20:04 by vgtyfty
[考研] 267求调剂 +5 再忙也要吃饭啊 2026-04-09 5/250 2026-04-09 18:47 by stone_128
[硕博家园] 新一代电子信息294求调剂 不挑学校 +5 Ytyt11 2026-04-09 6/300 2026-04-09 14:40 by Ytyt11
[考研] 311求调剂 +6 surte 2026-04-08 13/650 2026-04-09 14:00 by surte
[考研] 生物学308分求调剂（一志愿华东师大） +13 相信必会光芒万� 2026-04-06 16/800 2026-04-09 13:54 by 徐良白眉大侠
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 327求调剂 +10 Xxjc1107. 2026-04-06 11/550 2026-04-09 01:21 by lature00
[考研] 软件工程求调剂22软工296分求调剂，接受跨调 +4 yangchen2017 2026-04-08 5/250 2026-04-08 21:56 by 土木硕士招生
[考研] 277、学硕，求调剂 数一104， +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
[考研] 081700，311，求调剂 +17 冬十三 2026-04-04 18/900 2026-04-07 12:50 by Sammy2
[考研] 285求调剂 +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by 蓝云思雨
[考研] 调剂 一志愿吉林大学357分 +5 .Starry. 2026-04-04 5/250 2026-04-06 09:28 by cql1109
[考研] 生物与医药086000调剂一志愿西北农林320分 +3 美美女士 2026-04-03 3/150 2026-04-05 21:55 by 学员8dgXkO
[考研] 308求调剂 +4 maverick^_^ 2026-04-03 4/200 2026-04-05 19:08 by 蓝云思雨
[考研] 考研生物学考A区211，初试322，科目生化和生物综合，求调剂 +6 。。。54 2026-04-03 6/300 2026-04-05 14:54 by JOKER0401
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-03 5/250 2026-04-03 15:06 by arrow8852
信息提示
关闭
请填处理意见
关闭