24小时热门版块排行榜    

查看: 1616  |  回复: 3

yuelingxiao

新虫 (小有名气)


[交流] DPD 关于粒子数密度的取值

甚是迷惑,请大侠指点
在C.A.Marsh等人所写【Dissipative particle dynamics: The equilibrium for finite time steps】推导出:在时间步长确定的情况下,为使模拟结果稳定,存在一个密度上限。如图‘密度上限’所示。
在Rbert D.Groot等人所写【Dissipative particle dynamics: Bridging the gap between atomisticand mesoscopic simulation】中,又得出密度的一个下限(density>2)。如图‘密度下限’中所示。
不知道这两个密度是不是同一个值?
另外也不清楚在第二个文献中的密度是如何计算的?density=N/L^3,N:粒子总数, L模拟空间边长。其中L是不是应该是一个无量纲的呢,肯定不是真实值(若取为真实值的话一个立方米的空间内才3个粒子,不符合),那它到底是和谁的无量纲呢?

密度上限



密度下限

[ Last edited by ghcacj on 2011-11-24 at 21:25 ]
回复此楼

» 收录本帖的淘帖专辑推荐

和毕设有关的 Immiscible的淘贴

» 本帖@通知

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
yuelingxiao(金币+1):谢谢参与
ghcacj(金币+2): 谢谢 2011-11-25 10:32:16
引用回帖:
1楼: Originally posted by yuelingxiao at 2011-11-24 19:42:54:
甚是迷惑,请大侠指点
在C.A.Marsh等人所写【Dissipative particle dynamics: The equilibrium for finite time steps】推导出:在时间步长确定的情况下,为使模拟结果稳定,存在一个密度上限。如图‘密度上限’ ...

我没读过你提到的第一篇文章。如果你不是做方法的,那没有必要纠缠于所谓的密度上下限,采用常用的3$r_0^{-3}$即可。

DPD中的基本单位为$r_0$, $t_0$ 和$m_0$,这在Groot的文章中有提及。
2楼2011-11-25 03:04:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuelingxiao

新虫 (小有名气)


引用回帖:
2楼: Originally posted by yahoohoo at 2011-11-25 03:04:30:
我没读过你提到的第一篇文章。如果你不是做方法的,那没有必要纠缠于所谓的密度上下限,采用常用的3$r_0^{-3}$即可。

DPD中的基本单位为$r_0$, $t_0$ 和$m_0$,这在Groot的文章中有提及。

谢谢您的解答,还是不能完全清晰。
这个基本单位$r_0$和$t_0$指的是什么?
$m_0$该指DPD粒子的质量吧?

Groot在哪篇文章中有说明?【Dissipative particle dynamics: Bridging the gap between atomisticand mesoscopic simulation】这篇里没有呀?

期待回音
3楼2011-11-27 12:32:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuelingxiao(金币+4): 2011-11-28 09:55:18
ghcacj(金币+1): 谢谢 2011-11-28 12:02:20
引用回帖:
3楼: Originally posted by yuelingxiao at 2011-11-27 12:32:17:
谢谢您的解答,还是不能完全清晰。
这个基本单位$r_0$和$t_0$指的是什么?
$m_0$该指DPD粒子的质量吧?

Groot在哪篇文章中有说明?【Dissipative particle dynamics: Bridging the gap between atomistican ...

是这篇文章,多读几遍就理解了。
4楼2011-11-27 17:11:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yuelingxiao 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:38 by k06BKNblrGRH
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 Aj1rhIDL5ixY 2026-09-14 5/250 2026-09-17 05:14 by k06BKNblrGRH
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 5/250 2026-09-17 05:14 by k06BKNblrGRH
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 04:52 by k06BKNblrGRH
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:26 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-17 04:16 by k06BKNblrGRH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-17 04:04 by k06BKNblrGRH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 6/300 2026-09-17 02:38 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 6/300 2026-09-17 01:50 by Tql5LQhh5rLK
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-17 01:26 by Tql5LQhh5rLK
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-17 01:16 by Tql5LQhh5rLK
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +8 QUjhNVAcOSff 2026-09-13 11/550 2026-09-16 18:37 by rRvxEj3zorCy
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 16:51 by i3KBRD19ERkG
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 08:51 by SjQobnC04j06
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-16 02:38 by DgNGHc3h5tPl
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见