多孔介质渗流问题

vfengchuiguo

花了一上午时间把开版以来的帖子看了遍扫了扫盲避免问一些让大家笑掉牙的问题。但是•••心中疑惑还是没有解决，我觉得还是写下来博大家一笑好了。
1.LBM实验设计及收敛性：
最近刚看完Lattice Boltzmann Modeling：An introduction for geoscientists and engineers。
作者主要是用LBGK模型。书中有一例子，单组分单相流（D2Q9模型）过程中为了保证收敛性，Poiseuille Flow的模拟参数确定过程时这样的：
① 格子速度umax要小于0.1 lu/ts，这里取0.1
② 一般松弛时间τ=1，进而运动粘度υ=1/3（τ-1/2）=1/6 lu^2/ts,.
③ 由公式umax=3/2u确定平均流速u，进而根据Re=u 2a/υ得到特征长度2a。
现在我要设计一个多孔介质渗流模拟实验，可能要进行③->②->①的反过程，不知道这样做对不对：
Step1由于多孔介质图像已经确定，所以格子区域也确定了，继而得到特征长度。（顺便问一下，多孔介质的特征长度是一般怎么取。取宽度是不行的吧？）。
Step2 然后给定一个雷诺数Re。
Step3 根据step1和Step2，u/υ的比值确定了，但是多孔介质渗流没有找到类似泊肃叶流的umax与u的关系式，只能限制υ取1/6了。 从而得到参考速度u，判断u是否太大，如果太大稍微调节一下υ。
2.关于边界的问题，比较菜见谅。
如果用速度边界的话，进出口速度需要一样么（为了保证质量守恒？）
我怎么“无意”还看到了有人进口用压力边界，出口速度边界？错的？
3.两相流的问题
多孔介质求渗透率有两种方式实现，一种是给定压力边界，一种是施加定常的体力。
以施加体力为例，利用达西定律求某相相对渗透率的时候，是否意味着我要施加定常体力模拟一次该相单相流时候的结果得到流速，再施加同样的体力模拟一次多相流得到该相流速及饱和度，用流速之比得到相渗。若要求得相渗曲线，我就需要施加不同体力针对不同相得到不同的点。额，好麻烦，各位大牛们都怎么求呢。

总之由于只是纸上谈兵，好多是自己臆想的，还请赐教。下面会进一步学习MRT，阅读Pan的文章。。。苦恼ing

小七工作室

你是做那个方面扥?有具体模型?

vfengchuiguo

多孔介质，因为项目涉及到页岩油气，可能要用LB方法，主要还在学习与调研。 首先多孔介质的基本参数得求一下，相渗、毛管力曲线等等。然后再做一些多尺度下渗流机理的的工作，模型还没有吧 阁下有什么高见？

luo@odu.edu

The following paper may be of help to you:

C. Pan, L.-S. Luo, and C.T. Miller.
An evaluation of lattice Boltzmann schemes for porous medium flow simulation.
Computers & Fluids 35(8/9):898-909 (September, 2006).