《Chinese Journal of Aeronautics》：A GPU-accelerated implicit and consistent Multiphase Lattice Boltzmann Flux Solver for incompressible interface dynamics
编辑推荐：
为解决传统多相格子玻尔兹曼通量求解器（MLBFS）在压力-速度耦合及界面捕捉上面临的数值挑战，研究人员开展了基于一致性MLBFS（CMLBFS）的双时间步长界面捕捉方法研究。该研究通过重构宏观变量、修正动量方程对流项及引入双时间步长压力修正，显著提升了数值稳定性与界面分辨率，为高密度比、高雷诺数多相流模拟提供了高效工具。
在计算流体力学领域，多相流界面捕捉一直是个“精细活”与“力气活”并存的难题。无论是航空航天领域的燃油晃动，还是微流控芯片中的液滴操控，科学家们总希望能像用高清摄像机拍摄一样，精确捕捉到流体界面的每一丝变化。然而，现实往往很骨感：传统的数值方法要么因为复杂的压力泊松方程求解而计算量巨大，要么在模拟高密度比、高雷诺数流动时出现数值“崩溃”。特别是基于Cahn-Hilliard（C-H）相场模型的格子玻尔兹曼方法（LBM），虽然能优雅地描述界面演化，却常常受限于均匀网格和显式格式，难以胜任复杂工程计算。
正是在这样的背景下，一项发表于《Chinese Journal of Aeronautics》的研究提出了一种名为“一致性多相格子玻尔兹曼通量求解器（CMLBFS）”的新方法，并引入了双时间步长策略，旨在攻克上述瓶颈。该研究不仅重新推导了更严谨的宏观方程，还通过巧妙的“分治”策略处理对流项，并利用“伪时间”迭代解决了压力场的精度问题，最终实现了对复杂多相流动的高效、高分辨率模拟。
关键技术方法
本研究构建了基于C-H相场模型的一致性多相格子玻尔兹曼通量求解器（CMLBFS）框架。针对传统MLBFS的不足，研究采用了三项关键技术：1）一致性通量重构：利用LBM分布函数重构宏观通量，并统一用于C-H方程和Navier-Stokes（N-S）方程，修正了动量方程中对流项带来的误差；2）对流项分裂离散：对C-H方程中的对流项采用“分治”策略，其中序参量采用三阶WENO（加权本质无振荡）格式重构以保证界面锐利，速度则直接采用宏观重构速度；3）双时间步长压力修正：引入伪时间维度，在物理时间步内通过隐式BDPLUR（块数据并行下上松弛）算法迭代求解压力演化方程，确保压力-速度耦合精度，摆脱了传统ACM（人工压缩性方法）的局限性。
主要研究结果
1. 一致性离散框架的构建与验证
研究人员首先从多相LBM出发，通过Chapman-Enskog展开恢复了包含压力演化方程、动量方程和C-H方程的宏观方程组。与传统MLBFS不同，CMLBFS将LBM视为纯粹的“通量重构器”，并将重构出的宏观变量（如密度、速度）一致性地用于所有方程的离散。为了验证该框架的准确性，研究团队模拟了经典的二维液滴变形问题。结果显示，在不同雷诺数（Re）条件下，CMLBFS计算得到的液滴长轴和短轴随时间的变化曲线与文献结果高度吻合，且显著优于未进行一致性修正的原始MLBFS。这表明，通过一致性重构修正动量方程中的对流项，有效降低了因C-H方程与质量守恒方程不完全等价而引入的误差。
2. 双时间步长策略提升压力场精度
针对显式时间推进难以满足真实流动中压力-速度耦合关系的问题，研究提出了双时间步长方法。在每个物理时间步内，压力演化方程仅在伪时间内进行迭代求解直至收敛，而不随物理时间显式推进。通过对比单时间步长（显式）与双时间步长（隐式）在稳态和非稳态流动中的表现，研究发现双时间步长CMLBFS能获得更接近不可压缩N-S方程理论解的压力场。特别是在模拟高密度比（1000:1）的Rayleigh-Taylor不稳定性时，该方法成功捕捉到了清晰的界面失稳和涡旋结构，而传统方法则容易出现数值振荡。
3. 复杂多相流动的鲁棒性测试
为了检验算法的极限性能，研究团队挑战了多个高难度算例。在三维液滴碰撞案例中，CMLBFS精确再现了液滴碰撞后的融合与反弹过程，界面厚度在整个模拟过程中保持稳定，未出现非物理的扩散或收缩。在带障碍物的液滴撞击流动中，即使网格是非均匀的，该方法也表现出了良好的网格适应性，准确预测了液滴的飞溅模式。这些结果证明了该方法的鲁棒性，能够处理包含大密度比、高雷诺数和复杂几何边界的真实工程问题。
研究结论与意义
该研究成功开发了一种基于一致性离散和双时间步长策略的多相格子玻尔兹曼通量求解器（CMLBFS）。该方法的创新之处在于：1）通过一致性通量重构，消除了传统方法中不同方程离散格式不一致带来的数值扩散，提升了界面捕捉的精度；2）引入双时间步长压力修正，解决了显式格式下压力-速度耦合弱的问题，获得了更物理的压力场；3）保持了LBM无需求解压力泊松方程、易于并行计算的优势。
这项研究的重要意义在于，它为模拟航空航天、能源化工等领域中常见的复杂多相流动提供了一种高精度、高效率的数值工具。相比于传统的基于N-S方程的求解器，CMLBFS计算复杂度更低，更易于在GPU等高性能计算平台上实现加速，有望在未来的工程设计与分析中发挥关键作用。
