[求助]为什么我用LBM算出来的速度那么小？

didu

LBE looks simple. But it is not. If u have not done some test cases and compared with some exact/computational solutions, it is hard to find the problems. I met some people who are/were using LBE. They could not know how to achieve real physics.
It';s really worth to do some tests, trust me.

didu

Also, it is impossible that the code is unstable when tau is around 1 and initial conditions is correct.That also implied that your code has problems.
So check u initial condition and debug your code with a channel flow, which has periodic boundaries at inlet and outlet.

飘风

LBM中速度是无量纲的,其参考速度也就是粒子的迁移速度,接近于声速,对于不可压缩流动而言,无量纲速度应该远小于1
TAO是松驰系数,无量纲粘性表示为TAO-0.5)/3,TAO接近0.5时粘性很小,可能导致不稳定.

caiqd

在LBM中，速度1相对于“光速”，所以，计算出的速度非常小是正常的。需要选择一个基准速度点，即把所有的速度乘一个适当的系数。

lien

各位达人，我找到了问题的所在。主要是我在initialize的时候，局部的赋予了初始密度和速度。
我把整个求解域的所有cell的密度初始化为1，速度初始化为0，然后在stream和collision的过程中添加一个inlet，如z=5处，设置宏观速度为0.06。
结果看起来还像那么回事，放进去的粒子可以显示出涡街现象，计算到万步之后还是会出现计算的不稳定，出现负的密度值，这个可能是正常的么？
在调试的过程中，发现如果不初始化所有节点的密度为某个值，而只是初始化部分节点，结果是很快就会出现负的密度值。对于单相流，这个可以理解，比如气体，因为气体首先肯定就充满了整个求解域。但是对于两相流，比如两种液体a和b，最开始的时候，a在求解域内的某个区域A内不存在，即在此区域内a的密度为0，在A内b的密度为1，A之外的其它区域a的密度为1，b的密度为0。这样的情形，在现实中存在，但是在lbm中，我这样初始化的话，很快就出现了负密度值的情形（两相流lbm参考的是Diffusion in a multi-component Lattice Boltzmann Equation, Xiaowen Shan etc.）。请问大家在用lbm计算两相流的时候，是如何作的？？