两相格子动力学方案中压力计算的数值研究：显式与隐式方法的比较

《Computers & Fluids》：A numerical study of pressure calculation in a two-phase lattice kinetic scheme: A comparison between explicit and implicit methods

　　

在本研究中，研究人员（研究人员）对改进的格子动力学方案（improved lattice kinetic scheme, LKS）中一种特征显式压力计算方法进行了数值评估。该显式方案利用压力演化方程，并结合内迭代以增强计算稳定性。为进行对比分析，研究人员引入了一种隐式方法，通过简化标记与网格（simplified marker and cell, SMAC）方法求解泊松方程（Poisson equation）来计算压力。主要目标是阐明改进的LKS中内迭代次数与声波误差（acoustic error）之间的关系，并通过与隐式方法的比较分析加以验证。研究人员在相同条件下对三个典型两相流问题进行了基准模拟：(i) 静止液滴（stationary droplet），(ii) 单上升气泡（single rising bubble），以及(iii) 瑞利-泰勒不稳定性（Rayleigh–Taylor instability）。研究结果表明，改进的LKS中内迭代次数不仅影响数值稳定性，还影响声学尺度上的压力振荡。具体而言，声波的频率依赖于内迭代次数。特别是在重力等外力作用下，使用较多迭代次数会导致声波频率降低且衰减更快。隐式方案成功抑制了此类声波尺度效应，使压力演化更为平滑；然而，在计算效率方面，显式方法显著优于隐式方法。这些发现阐明了在两相模拟中数值稳定性、声波误差与计算成本之间的权衡关系。

**论文解读：两相改进格子动力学方案中显式与隐式压力计算方法的数值对比研究**

**研究背景与问题**
不可混溶两相流在化工反应器气液相互作用、传热冷却技术、海啸期间空气卷吸、推进系统空化以及生物医学中微泡输运等工程领域广泛存在。这些流动常涉及界面变形、相互作用和不稳定性等复杂物理过程。准确理解与定量预测两相流行为对工程系统的安全、性能优化和精确设计至关重要。数值模拟是研究两相流的有效手段，计算流体动力学（CFD）中的压力-速度耦合方法（如投影法、SMAC法、SIMPLE法）已较为成熟，但传统方法需求解压力泊松方程，该过程涉及稀疏矩阵处理与线性求解器收敛问题，在高分辨率、大规模图形处理器（GPU）模拟中成为关键瓶颈。

格子玻尔兹曼方法（LBM）因其显式格式、无需求解压力泊松方程、算法简单且适合并行计算而受到关注。格子动力学方案（LKS）作为LBM的一种变体，具有单松弛时间（SRT）模型的简洁性。改进的LKS通过引入超粘性和声波误差抑制，成功模拟了大密度比两相流。其关键特征在于利用压力演化方程显式计算压力，并通过内迭代和加速参数调整各相声速以增强稳定性。然而，LKS本质上是弱可压缩方法，声波（acoustic wave）误差不可避免，尤其在涉及外力时。以往研究尚未数值探究内迭代策略与声波行为之间的具体关系。

**研究目的与意义**
本研究旨在通过数值分析，阐明改进的LKS中显式压力计算方法的内迭代次数与声波误差之间的关联，并与隐式方法（基于SMAC法求解泊松方程）进行对比。研究意义在于：一方面揭示显式方案的声波误差特性，为可靠应用提供依据；另一方面对比两类方法的计算效率，为大规模两相流模拟选择合适策略提供参考。论文发表在《Computers》。

**主要关键技术方法**
研究人员采用了以下关键技术方法：
- 改进的LKS框架：基于D2Q9速度模型，采用保守Allen-Cahn（AC）方程（一种相场模型）追踪两相界面，通过压力演化方程显式计算压力，并引入加速参数ω和内迭代次数n调整声速。
- 隐式SMAC方法：将SMAC方法引入改进的LKS，通过求解压力泊松方程获得压力，作为对比基准。
- 基准测试：选取三个典型两相流问题——静止液滴、单上升气泡、瑞利-泰勒不稳定性，在相同条件下用显式和隐式方案模拟，对比压力场、速度场、界面变形及速度散度。
- 离散傅里叶变换：定量分析压力时间序列的声波频率特性。

**研究结果**
**3.1 静止液滴**
通过模拟静止液滴验证表面张力计算精度。显式方法（n=2, ω_max=400）与隐式SMAC方法均能准确复现拉普拉斯压力（Laplace pressure），误差分别为1.69%和1.78%，差异在2%以内。但计算时间差异显著：显式方法每时间步0.003秒，隐式方法（红黑SOR）需0.54秒，显式方法效率高约180倍。此外，对比了压力泊松方程迭代求解器的收敛速度：LBM求解器收敛最慢，红黑SOR最快。

**3.2 单上升气泡**
模拟密度比10的上升气泡，分析内迭代次数n的影响。设置显式方案case1（n=2, ω_max=5）和case2（n=5, ω_max=2），满足nω_max=ρ_L。界面变形、气泡质心垂直位置Y_c和平均上升速度u_c在三种情况下几乎重叠，表明声波对宏观流动影响甚微。但压力场差异显著：显式方法因有限声速产生声波振荡，隐式方法通过泊松方程消除声波。压力时间序列的离散傅里叶变换显示，case1（n=2）峰值频率高于case2（n=5），表明内迭代次数增加可抑制高频声波生成。声波随时间逐渐衰减，最终与隐式方法结果一致。速度散度方面，显式方法误差分布全域，隐式方法误差集中在界面附近。

**3.3 瑞利-泰勒不稳定性**
模拟低密度比（case1: ρ_L/ρ_G=3）和高密度比（case2: ρ_L/ρ_G=1000）两种情况。界面演化：低密度比下出现开尔文-亥姆霍兹不稳定性（Kelvin–Helmholtz instability）导致涡旋与液滴破碎，高密度比下落滴直线下落，与文献一致。显式与隐式方法在整体界面位置上吻合良好，但小液滴运动存在细微差异。压力场：由于上下壁面边界，显式方法中声波被反复反射，持续振荡而不衰减；隐式方法无此现象。离散傅里叶变换显示显式方法有明确高频峰值，隐式方法高频分量被抑制。速度散度特征与上升气泡类似。此外，内迭代次数增加（如n=5）在超粘性系数λ=1时仍不稳定，表明高雷诺数不稳定需其他稳定策略。

**讨论与结论**
讨论部分强调两种方案各有优劣：显式方案计算效率极高，适合大规模GPU模拟；但不可避免地产生声波误差，内迭代次数会影响声波频率和衰减速率，需综合考虑稳定性与误差控制。声波对宏观流动和界面追踪影响较小，但对小气泡或液滴运动影响更明显，导致与隐式方案结果差异。隐式方案通过求解泊松方程消除声波，压力演化更平滑，但计算成本高。

研究结论翻译如下：
在本研究中，研究人员对改进的LKS用于两相流[51,52]中的压力计算方法进行了数值研究。改进的LKS的主要优势在于使用压力演化方程，避免了求解计算昂贵的压力泊松方程。研究人员重点数值分析了该方案特征稳定化技术（即使用加速参数和压力计算的内迭代）的影响。为便于研究，研究人员将基于SMAC方法[10]的隐式压力求解器引入改进的LKS框架，提供了与传统求解压力泊松方程方法相反的对比方法。通过显式与隐式压力计算方案的对比分析，研究人员成功地从传统隐式求解器的视角凸显了显式方法的优势与局限。研究人员证明，显式和隐式方案各有优缺点。显式方案最显著的优势在于其卓越的计算速度（即效率）。近年来，GPU作为计算平台不仅在流体模拟中，也在其他领域变得越来越普遍。因此，高分辨率模拟中泊松方程的收敛性差问题日益突出。随着基于GPU的大规模模拟持续扩展，隐式方案的计算成本将变得更为显著。研究人员的发现证实了显式方案非常适合大规模模拟。然而，显式方案不可避免地会产生虚假的声波，这是与隐式方案相比的一个显著缺点。虽然Inamuro等人[51,52]报道了压力计算中使用的内迭代次数影响数值稳定性，但研究人员进一步揭示了它还会影响声波的频率和阻尼率。因此，内迭代次数应谨慎选择，不仅要考虑稳定性，还要考虑声波误差的控制。研究人员还报告了声波的存在对流动速度和界面追踪结果影响可忽略不计。这是因为声波传播的时间尺度与流动演化和界面运动的时间尺度显著不同。然而，在涉及小气泡或液滴运动的情况下，影响变得更明显，导致与隐式方案存在差异。由于在显式方案中声波误差不可避免，因此在采用显式方法时严格评估此类误差的潜在影响至关重要。