CFD仿真VOF界面捕捉稳定性研究：数值扩散与表面张力伪速度的影响及优化

1. 界面稳定性问题分析

在VOF方法中，界面稳定性是多相流模拟的核心问题之一。界面捕捉的稳定性直接影响模拟结果的物理准确性，尤其在处理复杂界面形变、表面张力驱动流动以及高密度比流体时，数值误差可能导致非物理现象。以下是界面稳定性中两个主要问题的分析：数值扩散和表面张力伪速度。

 1.1 数值扩散

成因：

数值扩散是由于VOF方法中体积分数的对流项离散化误差引起的，尤其在使用低阶格式（如一阶迎风格式）时更为显著。

在界面捕捉过程中，界面形状会随着流动逐渐扩散，导致界面模糊，失去物理意义。

当网格分辨率不足时，界面重构（如PLIC方法）也可能引入额外的扩散误差。

影响：

界面厚度增加：数值扩散会导致界面从一个清晰的薄层变成模糊的过渡区域。

物理现象失真：例如在液滴动力学中，数值扩散可能导致液滴体积损失或形状畸变。

多相流模拟精度下降：尤其在小尺度问题（如微流体）中，数值扩散会显著影响界面动力学行为。

  1.2 表面张力“伪速度”（spurious currents）

成因：

表面张力的计算通常采用CSF（Continuum      Surface Force）模型[1]，该模型将表面张力分布到界面附近的网格中。如图 1所示，界面的离散近似充当物理平滑界面上的扰动，在其附近产生虚假的毛细管流。其毛细管流的量级与网格尺度有关，见式。