N-Component Free Energy Lattice Boltzmann Method with Reduction Consistency and Global Momentum Conservation

本文提出了一种新颖的 N 组分自由能格子玻尔兹曼方法，该方法确保严格满足归约一致性并将全局动量守恒精度提升至机器精度，在从静态液滴到复杂微流体应用的各种多相流模拟中均展现出高精度。

要点

想象你是一位厨师，试图在计算机上模拟一种复杂的汤。这碗汤不仅仅包含水和盐；它含有数十种不同的成分——油、醋、香料、香草——它们彼此之间并不相容。有些成分倾向于聚拢在一起，有些则希望彼此分离，而且它们各自拥有不同的“粘性”（粘度）和“排斥力”（表面张力）水平。

长期以来，计算机模拟一次只能处理两三种这样的成分。如果你尝试加入第四种，模拟就会陷入混乱。它可能会凭空捏造出一种新的成分，仅仅因为数学计算变得杂乱无章；或者整锅汤可能会开始在虚拟的厨房桌面上自行漂移，违背物理定律。

本文介绍了一种新的、更智能的方法来模拟这些“多成分”流体，所采用的技术称为格子玻尔兹曼方法（LBM）。可以将 LBM 想象成一个由微小方格组成的网格，流体粒子在这些方格之间跳跃。作者建立了一套新的规则来规范这些粒子的跳跃方式，确保以下两个关键事项得以实现：

1. “无幽灵成分”规则（还原一致性）

问题所在： 在以往的模拟中，如果你有一锅包含四种成分的汤，但只实际加入了三种，计算机可能会突然“幻觉”出第四种成分凭空出现。这就像你用面粉、糖和鸡蛋烤蛋糕，面糊却突然在没有添加任何可可粉的情况下变成了巧克力味。这会毁掉整个模拟。

解决方案： 作者制定了一条严格规则：如果某种成分不存在，它就绝不可能出现。他们通过在数学公式中加入一个“校正因子”来实现这一点。想象一下俱乐部门口的保安在核对宾客名单。如果名单上写着“禁止巧克力”，那么无论其他成分如何“舞动”，保安都会确保没有任何巧克力分子进入派对。这确保了：如果你移除第四种成分，四种流体的模拟行为将与三种流体的模拟完全一致。

2. “无漂移锅”规则（动量守恒）

问题所在： 在旧方法中，维持油水分离的力（表面张力）的计算方式存在微小的“泄漏”。这就像有一台微小而隐形的风扇对着你的汤锅吹。随着时间的推移，整锅汤会缓慢地在桌面上滑动，尽管无人触碰。这导致模拟结果不准确。

解决方案： 作者重新设计了这些力的数学计算，确保每一个方向的推力都被另一个方向的拉力完美平衡。这就像一场拔河比赛，绳子处于完美的中心位置；无论双方队伍如何用力，绳子都不会向左或向右漂移。这使得流体能够精确地停留在它应该所在的位置，精度达到计算机所能实现的极限。

他们进行的测试（“品尝测试”）

为了证明他们的新配方有效，他们运行了多项模拟：

• 液透镜： 他们将不同流体的液滴相互滴落，观察是否形成了正确的角度（就像油浮在水面上那样）。他们的模型与理论角度完美吻合。
• 贾纳斯液滴： 他们模拟了一种具有两个不同“面”的液滴（就像一枚有正反面的硬币）。旧方法会导致这些液滴发生漂移；而他们的新方法能让它们保持完全静止，直到应当移动的时刻。
• 分层流动： 他们模拟了六种不同流体层在管道中流动，每一层具有不同的厚度（粘度）。流动情况与数学预测完全一致。
• 相分离： 他们观察流体随时间分离的过程（就像瓶子里的油和醋分离）。他们的模型准确预测了分离发生的速度，符合现实世界的物理定律。

他们展示的实际应用

该论文证明，这种新方法能够处理涉及多种流体的复杂现实场景：

• 图案化液体表面： 他们模拟了一滴液体在覆盖着交替条纹的不同润滑流体的表面上移动。液滴会在条纹边缘“卡住”（钉扎），然后向前跳跃，这种行为用旧工具很难模拟。
• 微流控乳液： 他们模拟了一台能够制造“液滴中包含液滴”的微型机器（就像用液体制成的俄罗斯套娃）。他们的方法成功制造了一个包含流体 B 液滴的流体 A 液滴，而流体 B 液滴中又包含了一个流体 C 液滴。

核心结论

作者构建了一个稳健、无“幽灵”且无“漂移”的模拟器，适用于任意数量成分的流体。这使得科学家能够以前所未有的精度和稳定性研究复杂系统——例如蛋白质在细胞内的分离过程，或如何设计更优的药物输送液滴。他们不仅修正了数学公式，更使得在不导致计算机混淆的情况下，模拟流体那种混乱、多层的现实成为可能。

技术摘要

技术摘要：具有还原一致性与全局动量守恒的 N 组分自由能格子玻尔兹曼方法

问题陈述
涉及三种或更多不混溶组分的多组分流体流动，在从药物递送、微流体到生物相分离的各类应用中至关重要。尽管针对二元或三元系统的数值技术已经存在，但能够处理任意数量（$N$）不同相的方法却十分有限。主要挑战在于构建一个还原一致的宏观方程组：当一种流体缺失时，$N$组分系统必须精确复现$N-1$组分系统的行为。缺乏这一性质，缺失的流体组分可能会在三相点处自发成核，严重损害模拟的精度和稳定性。此外，现有的格子玻尔兹曼方法（LBM）通常采用无法在全局范围内守恒动量的表面张力力离散化方案，导致非物理的全域漂移。

方法论
作者构建了$N$组分流体的连续介质模型，并开发了相应的自由能 LBM。该模型求解不可压缩连续性方程和纳维 - 斯托克斯方程，并与$N-1$个 Cahn-Hilliard 方程耦合，以追踪流体体积分数（$C_i$）。

1. 通过通量修正实现还原一致性：
   不同于以往依赖特定退化形式扩散迁移率（这限制了扩散物理过程）的方法，该方法保持了自由的迁移率参数。为了确保还原一致性，作者在 Cahn-Hilliard 方程中引入了源项。具体而言，他们通过在化学势中添加修正项（$\phi_j$和$\xi_j$）来重构扩散通量$J_i$。

   • 项$\phi_j$旨在抵消错误的体通量，否则这些通量会导致由其他组分形成的三相点处出现缺失流体的成核。
   • 项$\xi_j$在展宽参数为正时，为自由能提供可选的稳定性，防止自由能在物理范围$C_i \in [0,1]$之外无界。
   • 这些系数基于表面张力矩阵解析计算得出，且不会影响模拟期间的计算性能。

2. 动量守恒的表面张力力：
   表面张力力$F_s = \sum \mu_i \nabla C_i$的标准数值实现，在使用有限差分模板时往往无法精确满足链式法则，从而破坏动量守恒。作者通过将力表示为应力张量的散度，推导出了动量守恒的离散化方案。为了避免计算完整应力张量散度的计算成本，他们修改了化学力公式为：
   $$F_s = \nabla E_B + \sum_{i,j; i \neq j} \kappa_{ij} \nabla^2 C_j \nabla C_i$$
   这种形式确保了力在域上的积分为零（达到机器精度），从而消除了非物理的域漂移。

3. 格子玻尔兹曼实现：
   该方法利用 D2Q9 模板和双弛豫时间（TRT）碰撞算子来处理粘度对比。流体速度和压力通过带有强迫项的标准 LBM 方程更新，而相场则使用一组包含修正化学势的独立分布函数（$g_{i,k}$）进行演化。

关键结果与基准测试
该方法经过一系列静态和动态基准测试的验证，与理论预测表现出极好的一致性：

• 液透镜与 Janus 液滴： 在四流体液透镜设置中，该方法正确捕捉了三相点处的 Neumann 角，与解析预测的偏差在$1^\circ$以内。至关重要的是，还原一致性项（$\phi_j$）防止了缺失流体（组分 5 和 6）的成核，并稳定了那些原本会因局部三相点触发不稳定的模拟。对于 Janus 液滴，动量守恒力消除了非守恒公式中观察到的非物理漂移。
• 相分离： 该模型成功复现了已建立的流体动力学（$L \propto t^{2/3}$）和扩散（$L \propto t^{1/3}$）标度律，用于描述域粗化。作者证明，先前通过浓度依赖迁移率强制执行还原一致性的方法未能恢复正确的扩散标度，而他们的通量修正方法则取得了成功。
• 分层泊肃叶流动： 具有不同粘度的六层流动模拟显示，模拟的速度剖面与从纳维 - 斯托克斯方程推导出的解析解之间具有极好的一致性。
• 应用：
  • 图案化液体表面（PaLS）： 该方法模拟了液滴在具有交替润滑剂斑块表面上的运动，捕捉了粘滑运动以及在变化体积力作用下从钉扎状态到运动状态的转变。
  • 乳液液滴生成： 模拟了用于生成三重乳液液滴的微流体设计，成功模拟了多种不混溶相的封装和 pinch-off（颈缩断裂）过程。

意义与主张
本文声称提出了首个能够模拟任意数量不混溶组分流体系统，同时严格强制执行还原一致性和全局动量守恒的自由能 LBM。

• 迁移率的独立性： 一个关键优势在于，还原一致性的强制执行独立于扩散迁移率。这使得能够精确控制扩散动力学并恢复已建立的相分离标度律，而先前的方法无法同时实现这两点。
• 消除漂移： 动量守恒的离散化方案解决了先前 LBM 方案中存在的显著速度漂移问题，确保了长时间模拟的物理保真度。
• 适用性： 作者将该方法定位为支持生物医学系统（如生物分子凝聚物）、软物质（如 DNA 纳米星）和微流体设备实验工作的工具。他们指出，该方法最适合于$N \lesssim 10$的问题，其中需要独立选择界面张力，例如最近涉及多达九个不混溶相的 DNA 纳米星实验。

作者总结道，虽然该方法在测试范围内表现稳健，但未来的扩展应解决与固体边界的润湿相互作用以及支持大密度对比的问题。