《Journal of Computational Physics》：A rotated-Riemann discrete-equations method for two-phase interfaces with phase transition
编辑推荐：
本文提出了一种基于鲁棒离散方程界面捕捉法（RDEMIC）的旋转黎曼求解器，用于带相变的两相流模拟，并采用HLLC与HLLP近似方法。在多维设置中，若界面与网格方向不一致，直接在网格方向上应用相变黎曼求解器会导致非物理解。为恢复旋转不变性，研究人员在界面对齐坐标
本文提出了一种基于鲁棒离散方程界面捕捉法（RDEMIC）的旋转黎曼求解器，用于带相变的两相流模拟，并采用HLLC与HLLP近似方法。在多维设置中，若界面与网格方向不一致，直接在网格方向上应用相变黎曼求解器会导致非物理解。为恢复旋转不变性，研究人员在界面对齐坐标系中求解界面黎曼问题，结合界面法向与切向解，以保持与网格对齐离散格式的一致性。相变黎曼求解器基于由动力学关系（kinetic relation）获得的界面质量通量进行构建。两种不同复杂度的动力学关系被耦合至求解器，可提供稳健的预测与快速的收敛性。涉及蒸发与冷凝的数值测试，从一维基准问题到激波-界面相互作用，均表明该旋转黎曼求解器能在存在质量输运的情况下准确捕捉界面动力学，并消除多维设置中非物理的网格对齐伪影。
论文解读：旋转黎曼离散方程法在相变两相界面模拟中的应用
一、研究背景与问题
可压缩多相流动的数值模拟是计算流体力学的重要研究方向，其核心挑战之一在于对相界面的准确捕捉与追踪。黎曼求解器作为Godunov型有限体积法的基石，通过在单元界面处求解局部黎曼问题，为双曲守恒律方程提供物理一致的数值通量。经典HLL求解器仅考虑两波结构，忽略了中间接触波，导致剪切波与接触间断的过度抹平。HLLC求解器通过恢复接触波，成为Euler方程与Navier–Stokes方程的完备黎曼求解器，并被广泛用于两相流动模拟。
当考虑相变效应时，界面处发生质量交换，接触波不再严格对应相界面。Fechter等人提出的HLLCP求解器引入了与相变相关的第四波，HLLP求解器则进一步简化，在相界面的蒸发或冷凝侧移除接触间断。这些求解器通过动力学关系（kinetic relation）将界面质量通量与界面两侧的流体状态耦合，从而闭合方程组。
然而，在有限体积法中，黎曼问题通常在计算单元界面处求解，求解结果天然与网格对齐。对于无相变的两相流系统，Euler方程的旋转不变性保证了网格对齐求解的合理性。但一旦引入相变，界面黎曼问题与界面无关区域的黎曼问题本质不同，网格对齐求解将导致非物理解。具体而言，当界面相对网格存在倾角时，一维流动配置在多维网格上的求解结果依赖于倾角，无法恢复正确的界面解。这一网格对齐伪影严重限制了相变两相流模拟的可靠性与普适性。
二、研究内容与意义
针对上述问题，本研究在鲁棒离散方程界面捕捉法（RDEMIC）框架内，将HLLCP/HLLP相变黎曼求解器与旋转黎曼求解器相结合，提出了一种与网格方向无关的相变两相流数值格式。该方法的核心理念是将界面黎曼问题变换到界面对齐坐标系中求解，再通过变换张量将通量与交换项投影回网格对齐方向，从而恢复系统的旋转不变性。数值测试涵盖蒸发与冷凝，从一维基准问题到激波-界面相互作用，均验证了该方法的准确性与鲁棒性。
本研究的意义在于：首次将旋转黎曼求解器思想引入相变两相流模拟，系统解决了相变界面求解对网格方向的依赖问题，为多维相变流动的高保真数值模拟提供了可靠工具。该方法在ALPACA框架中实现，结合块状多分辨率（MR）离散与自适应局部时间步进，具备面向高性能计算（HPC）系统的可扩展性。
三、主要技术方法
研究人员基于以下关键技术方法开展研究：
控制方程与离散化：采用两相（液相与气相）的离散方程系统，控制方程在向量形式下包含守恒量、通量张量与相互作用张量。相界面通过体积分数场表示，界面处通过相互作用张量耦合两相，其中包含质量通量向量与界面速度。采用Noble–Abel刚性气体（NASG）状态方程闭合方程组，该方程在经典刚性气体模型基础上引入排斥效应，并包含能量偏移参数以对齐两相内能在共同参考态的值。
空间离散与时间积分：在笛卡尔网格上进行空间离散，采用RDEMIC方法逐维计算数值通量与交换项。RDEMIC基于离散方程法（DEM），通过多个成对黎曼问题构造离散通量与非守恒项，包括单相与两相相互作用。时间积分采用三阶强稳定性保持Runge–Kutta格式，时间步长基于CFL条件确定。界面重构采用多维THINC（双曲正切界面捕捉）格式及其保正扩展，通过双曲正切函数重构体积分数，有效限制扩散界面的厚度。
HLLC/HLLCP/HLLP黎曼求解器：标准HLLC求解器通过Rankine–Hugoniot跳跃条件与接触波处的压力速度连续性确定中间状态。HLLCP求解器在相界面上引入第四波，HLLP求解器则简化接触波结构。相变求解器通过界面质量通量与动力学关系耦合，采用割线法迭代求解非线性方程组，以相对残差小于10?6为收敛标准。
动力学关系：采用两种不同复杂度的动力学关系。Hertz–Knudsen关系仅依赖于气液两相的热力学状态，包含蒸发与冷凝系数。Schrage关系在冷凝项中引入依赖于蒸气漂移速度的修正因子，考虑了蒸气相对界面的速度效应，增强了动力学关系与跳跃条件之间的耦合。
旋转黎曼求解器：构造界面法向单位向量与两个切向向量构成正交基，通过变换张量将网格对齐坐标系中的状态与通量变换到界面对齐坐标系，求解后再变换回网格对齐系统并投影到面法向方向，获得面通量。
四、研究结果
数值模型建立：研究人员扩展了原有的多相流系统，在相互作用张量中加入质量输运项。在无相变情况下，相互作用张量恢复为原始形式。通过引入体积平均状态与通量向量，交换项可解释为跨界面的拉格朗日通量。质量通量定义为密度与相速度和界面速度之差的乘积，且在两相间相等以满足质量守恒。
界面速度估计：由于相变发生时相速度与界面速度不同，研究人员基于质量通量定义与两相间质量守恒，推导了界面速度的估计公式。界面法向向量由符号距离函数的梯度计算。
HLLCP/HLLP求解器的构建：相变黎曼问题的波扇结构比标准黎曼问题更复杂，蒸发时在蒸气侧出现额外的中间蒸气区域，冷凝时在液体侧出现冷凝液体区域。HLLCP求解器解析全部四波结构，HLLP求解器在蒸发或冷凝侧移除接触波，仅保留相界面两侧各一个中间状态。通过相界面跳跃条件、质量通量与动力学关系闭合方程组，HLLCP还需引入密度比近似条件。
迭代求解策略：由于所有状态均以质量通量表示，可从初始值出发迭代更新。采用割线法加速收敛，典型情况下两到三次迭代即可达到10?6的相对残差，最多不超过五次。HLLP与HLLCP在所有测试案例中结果几乎不可区分，因此采用算法相对简单的HLLP。
守恒性质分析：HLLCP/HLLP构造的界面交换项满足跳跃条件，保证跨界面的总质量与动量守恒。能量跳跃条件包含相变潜热，界面能量传递包含非守恒贡献，反映蒸发与冷凝的热力学效应。无相变时该贡献消失，恢复原始严格守恒的界面相互作用。
旋转求解器的构造与验证：通过将界面黎曼问题变换到界面对齐坐标系求解，再变换回网格对齐系统，消除了网格对齐伪影。数值测试表明，旋转黎曼求解器在不同界面倾角下均能恢复正确的界面解，而网格对齐求解器的结果依赖于倾角。
数值测试：从一维基准问题到激波-界面相互作用，涵盖蒸发与冷凝过程，旋转黎曼求解器均准确捕捉了存在质量输运时的界面动力学。
五、讨论与结论
讨论部分指出，HLLCP与HLLP本身并不构成相变模型，相变模型完全编码在动力学关系中，因此这些求解器可与任何能表达为动力学关系形式的相变模型耦合。Hertz–Knudsen关系的主要局限在于忽略了蒸气相对界面的速度，导致蒸发质量通量被低估。Schrage关系通过引入依赖于蒸气漂移速度的修正因子解决了这一问题，但也增加了动力学关系与跳跃条件之间的耦合强度，验证了迭代求解策略处理更复杂动力学关系的能力。研究还指出存在更先进的动力学关系，如基于Knudsen层理论的模型，但详细讨论与比较超出了本研究范围。
研究结论表明，旋转黎曼求解器成功恢复了相变两相流系统在存在质量输运时的旋转不变性，消除了多维设置中非物理的网格对齐伪影。HLLP与HLLCP在测试中结果几乎不可区分，HLLP因算法相对简单而被采用。两种动力学关系均能提供稳健的预测与快速的收敛性，验证了该数值框架的灵活性与通用性。该方法在ALPACA框架中实现，结合多分辨率离散与自适应局部时间步进，为面向HPC系统的大规模相变两相流模拟提供了可靠工具。
