求教基于LBM的大密度单组分两相流模拟

sky.jun

我最近的工作基于Peng Yuan的《Equations of state in a lattice Boltzmann model》，采用Carnahan-Starling状态方程，当密度比达到10^3左右时，平板两相分离模拟的气相初始密度无法下去，只能通过迭代计算达到并保持平衡稳定，不会出现发散扰动。比如，临界温度为0.09433306467862，当温度T=0.055时，平板两相分离模拟结果为2.14E-04：0.41407，密度比为1935，初始气相密度只能设到2E-03，无法直接设成分离后平衡稳定值2.14E-04；对于液滴两相分离模拟，由于表面张力的存在，当温度T=0.053时，密度比才达到10^3，气相初始密度可直接设到3.2E-04，液相初始密度设为0.42，但是液滴半径有上限限制。

我想求教一下，有没有哪位DX采用过类似模型做大密度比单组分两相流模拟，在平板两相分离模拟中有没有遇到类似的气相初始密度下限的问题？谢谢各位了