韦伯数 — CAE用语解说

分类: 术语集 | 2026-01-15

CAE visualization for weber number - technical simulation diagram

韦伯数

什么是韦伯数

🧑🎓

韦伯数与表面张力有关系吗？与雷诺数和弗劳德数有不同的应用场景吗？

韦伯数的理论基础

韦伯数的定义和物理意义

🧑🎓

教科书中说韦伯数是"惯性力与表面张力的比"，但具体在什么现象中很重要呢？

🎓

在涉及界面变形的流动中，韦伯数起着决定性作用，如液滴分裂和气泡合并。例如，在柴油发动机的燃料喷射中，高压喷射的燃料液柱被空气阻力分裂成细小液滴（雾状），这种分裂模式由韦伯数控制。其公式为

$$ We = \frac{\rho u^2 L}{\sigma} $$

其中，ρ是密度，u是特征速度，L是特征长度（如液滴直径），σ是表面张力系数。

🧑🎓

特征长度L应该怎样选择？用液滴的初始直径，还是分裂后的卫星液滴直径？

🎓

好问题。通常选择能代表问题尺度的长度。讨论喷流的一次分裂时，使用喷嘴口径或液柱直径；处理二次分裂（液滴的进一步细分化）时，使用母液滴的直径。在实验中，通常用激光衍射法测得的Sauter平均直径（SMD, D32）作为L。例如，汽油喷射的研究中，喷嘴直径0.2mm、喷射速度100m/s、表面张力0.025N/m，则

$$ We = \frac{750 \times 100^2 \times 0.0002}{0.025} \approx 60,000 $$

这个非常大的值预示会发生剧烈的分裂。

🧑🎓

We很大时，表面张力的效应可以忽略吗？但液滴能形成意味着表面张力在某处起作用，对吧？

🎓

完全正确。在局部和微观上，表面张力始终在起作用。宏观上We很大，但分裂的最后阶段产生的微小液滴，其相对速度u变小，局部We接近1。此时表面张力的作用是将液滴集合成球形。换句话说，现象是多尺度的，需要分开考虑全局We和局部We。全局We描述分裂过程，局部We描述液滴稳定化。

韦伯数的数值计算方法

CFD中的界面追踪和We的作用

🧑🎓

在CFD中模拟液滴分裂时，韦伯数如何影响求解器的设置？

🎓

虽然韦伯数不是直接的设置参数，但作为无次元数，它主导解的行为，因此成为确定网格分辨率和时间步长的基准。例如，使用VOF法或Level Set法追踪界面时，为了精确计算表面张力项，需要足够细的网格来捕捉界面曲率。毛细数

$$ Ca = \frac{\mu u}{\sigma} = \frac{We}{Re} $$

与界面附近网格大小Δx遵循以下经验规则：

$$ \Delta x < \frac{\sigma}{\rho u^2} \cdot f(We) $$

在Ansys Fluent的VOF模型中，推荐使用"Continuum Surface Force"法，该法考虑相邻单元间的梯度进行曲率计算。

🧑🎓

听说表面张力的计算容易不稳定甚至发散。We很大和We很小时，哪个在数值上更困难？

🎓

数值上，We非常小的领域（表面张力主导）更困难。因为当表面张力占主导时，规定界面形状的压力跳跃

$$ \Delta p = \sigma \kappa $$

（κ是曲率）变得非常灵敏，界面的任何数值"噪声"都会导致大的压力波动，从而导致发散。为了稳定计算，很多求解器（如OpenFOAM的`interFoam`）采用"曲率平滑"或"括号法"来处理表面张力项。相反，We非常大的领域中，表面张力可以实际忽略，普通流体求解器可以处理，但代价是生成的微小尺度涡旋和复杂界面变形需要极高的网格分辨率。

韦伯数的实务应用

喷雾仿真的设置流程

🧑🎓

实际开始进行燃料喷射仿真时，有考虑We的工作流程吗？

🎓

有的，推荐的步骤如下。
1. **事前估计**: 根据喷嘴直径、喷射速度、物性值计算全局We。如果是10^4～10^5量级，说明一次分裂很剧烈，界面追踪法（VOF等）的计算成本可能巨大。
2. **模型选择**: 全局We很大时，只在喷嘴附近用VOF详细计算一次分裂，下游的二次分裂及以后使用拉格朗日粒子法（Discrete Phase Model）替代，采用混合方法（如Fluent的VOF-to-DPM）更现实。
3. **无次元数标度**: 当实际尺度的计算困难时，通过调整模型尺度或流体物性使We和Re一致，进行缩尺模型仿真。

🧑🎓

结果验证时怎样用韦伯数？只是简单比较液滴大小与实验数据？

🎓

比较液滴粒径分布（PSD）很重要，但更根本的做法是用无次元数验证标度规律。例如，在一次分裂区域，将测得的Sauter平均直径D32除以喷嘴直径d得到无次元化的

$$ \frac{D32}{d} $$

，以此为纵轴，全局We为横轴作图。许多实验数据遵循

$$ \frac{D32}{d} \propto We^{-0.5} $$

到

$$ We^{-0.6} $$

的相关性。检查仿真结果是否落在这条趋势线上。同时，用高速摄像头视频与仿真结果的界面形状（如液柱破裂长度等）进行定量比较，以We为参数。

韦伯数的软件比较

各求解器的界面处理和无次元数

🧑🎓

Ansys Fluent、STAR-CCM+、OpenFOAM在处理与韦伯数相关的界面上有差异吗？

🎓

基础方程没有大的差异，但表面张力项的离散化方法和稳定化技术各有特色。
**Ansys Fluent**: VOF法为主流。表面张力模型除了"CSF"法外，还有更高精度的"CSS"法（从相邻单元重构曲率）可选。在中等～高We的喷雾问题中，"VOF-to-DPM"功能很强大。
**Siemens STAR-CCM+**: 拥有"VOF Wave"和"Eulerian Multiphase"模型。其"Segregated Flow"求解器有表面张力项的隐式处理选项，对We较小（表面张力主导）的流动，计算稳定性更高。
**OpenFOAM**: 标准求解器是`interFoam`。表面张力项的计算包括"界面锐化保持法"和"曲率计算平滑化"的平衡，可细致设置。用户可在源代码级改动算法，常用于针对特定We范围的定制研究。

🧑🎓

COMSOL Multiphysics这样的基于FEM的软件呢？界面追踪的方法似乎不同。

🎓

完全同意。COMSOL主要采用"移动网格（ALE）法"和"Level Set法"。尤其是Level Set法，将界面表示为隐函数的等值面，能自然处理界面的拓扑变化（合并、分裂），曲率计算也很连续。因此，在We中等、界面形状复杂变形的问题（如微流体装置中的液滴生成）上较强。但FEM的计算成本比FVM求解器高，不太适合We非常高的乱流喷雾等极为动态的现象。COMSOL在物理场界面中直接设置"表面张力"，无次元数用于结果评估。

韦伯数的故障排除

界面仿真的收敛问题

🧑🎓

液滴在计算域内移动的仿真中，液滴形状变得崎岖，甚至计算发散。We约为5。原因和对策是什么？

🎓

We=5属于表面张力影响强的区域。最常见的原因是"寄生流"。这是界面曲率计算的误差与表面张力项离散化的不匹配导致界面附近产生非物理微小涡的现象。对策如下。
1. **网格改进**: 界面至少用5～10个单元分解。结构网格为佳。
2. **压力-速度耦合求解**: 相比SIMPLE法，PISO法或COUPLED法对这类问题稳定性更高（Fluent的情况）。
3. **表面张力模型调整**: 设置用于曲率计算的界面平滑宽度。太窄对噪声敏感，太宽界面形状会变钝。
4. **减小时间步长**: 根据毛细时间常数

$$ t_c = \sqrt{\frac{\rho L^3}{\sigma}} $$

设置时间步长为其1/10以下。

🧑🎓

相反，高速喷流的VOF计算We超过1000时，界面异常扩散，细液滴无法再现。

🎓

这是"数值界面扩散"。高We时界面变形剧烈，普通VOF输运格式导致界面模糊。对策有两个。
1. **使用界面锐化保持法**: Fluent的"Geo-Reconstruct"、OpenFOAM的`MULES`算法等特殊对流格式维持界面锐度。
2. **自适应网格加密（AMR）**: 动态细分界面所在单元及其周围，如STAR-CCM+的"Adaptive Mesh Refinement"或Fluent的"Dynamic Mesh Adaptation"。这样既保持界面分辨率，又控制总单元数。
根本上，这个We区间用VOF捕捉整个液滴有极限。现实的做法是，用VOF捕捉一次分裂的大局行为，微小液滴用粒子法建模，采用混合方法。