Shan-Chen模型

7楼 · 发表于 2011-4-20 21:01:57

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？注册

x

drho = 0.001;
delta_rho = -drho*(1-2.0*rand(lx));
for i=1:9
fIn(i,1:lx,1:ly) = tNS(i).*(1.0 + delta_rho);
gIn(i,1:lx,1:ly) = tNS(i).*(1.0 - delta_rho);
end

关于fIn和gIn里密度的一加一减是什么意思？

nanosmile · 发表于 2011-4-23 12:59:46

初始化过程
随机生成delta_rho，范围在[-0.001,0.001]
然后通过f=tNS(i)来初试化分布函数
一加一减只是增加一些随机波动
这样处理之后各点的密度仍为1

nanosmile · 发表于 2011-4-23 13:06:59

楼主用的模型是那篇文献里的？

7楼 · 发表于 2011-4-23 15:17:09

应该是模拟一个多孔骨架出来的，但是没有看到空隙率或者其他的类似初始条件。
在网上找的，第三个
http://www.lbmethod.org/doku.php?id=numerics:matlab_samples

nanosmile · 发表于 2011-4-23 19:19:53

这个程序是模拟两相流体的
初始化是在一个正方形区域内，两流体1,2均匀分布，
就是你问题里的命令。流体1,2的平均密度都为1，每个点的密度在1附近波动，
幅值为10e-3
没有其它力作用，只是通过1,2间的分子排斥力使两流体分离。
程序里好像没加边界条件，可以加个反弹边界试一下。
matlab我不太懂，大概是这个意思。好像是这篇文章里的模型。

7楼 · 发表于 2011-4-23 22:31:11

非常感谢，以及无私的分享文献！

		自动登录	找回密码
密码			注册

Shan-Chen模型

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复 1# 7楼的帖子

回复 1# 7楼的帖子

那初始条件是什么？

回复 4# 7楼的帖子

回复 5# nanosmile 的帖子

站点信息

Shan-Chen模型

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复 1# 7楼 的帖子

回复 1# 7楼 的帖子

那初始条件是什么？

回复 4# 7楼 的帖子

回复 5# nanosmile 的帖子

站点信息

回复 1# 7楼的帖子

回复 1# 7楼的帖子

回复 4# 7楼的帖子