CFD冲蚀模拟的一些理论

在CFD中计算颗粒对固体壁面的冲蚀往往采用冲蚀模型（Erosion Model）。

1 冲蚀速率（Erosion Rate）

冲蚀速率定义为壁面材料在单位时间单位面积上损失的质量（单位：\(kg/(m^2 \cdot s)\)）。通过计算每一个颗粒对壁面的累积损伤来计算冲蚀速率。
\[ E_f =\frac {1}{A_f} \sum_{\pi(f)}\dot m_\pi e_r \]

式中：

• \(A_f\)为网格单元面积
• \(\dot m_\pi\)为冲击壁面的颗粒质量流量
• \(e_r\)为冲蚀率

冲蚀速率取决于颗粒流动状态（以何种方式冲击壁面）以及选择的冲蚀率。

2 STAR CCM+中的冲蚀率模型

STAR CCM+内置了四种常用的冲蚀率：

• Ahlert Correlation
• DNV Correlation
• Neilson-Gilchrist Correlation
• Oka Correlation
  这四种模型定义了由于颗粒直接冲击壁面造成的冲击损伤，还可以利用Archard Correlation则定义了由于颗粒切向或以低角度冲击壁面造成的磨蚀。

2.1 Ahlert Correlation

Ahlert Correlation[1]定义冲蚀率为：
\[E_r=KF_sf(\alpha)(\frac{u_{rel}}{u_{ref}})^n\]
式中：

• \(K\)为与材料相关的常数
• \(F_s\)为颗粒形状因子。对于完全球形颗粒取0.2，对于半圆形颗粒取0.5，不规则颗粒取1
• \(f(\alpha)\)为颗粒冲击角函数。冲击函数被分为两部分：冲击角低于转变角度\(\alpha_0\)时表示为冲击角\(\alpha\)的多项式函数（角度单位为弧度），高于转变角度时，遵循三角函数关系：
  \[ f(\alpha)=xcos^2\alpha sin \omega \alpha + ysin^2 \alpha + z \]
  式中常数\(\omega\)、\(x\)及\(y\)均为用户自定义，Z值为程序内部计算得到以满足曲线的连续性。
• \(u_{rel}\)为粒子与壁面间的相对速度，\(|u_{ref}|, u_{ref}=u_{particle}-u_{wall}\)
• \(u_{ref}\)为颗粒参考速度常数
• {n}为指数

STAR CCM+中Ahlert correlation默认参数为液体驱动半圆形砂颗粒侵蚀铝，参数采用文献[2]。

2.2 DNV Correlection

DNV Correlection[3]对于冲蚀率表示为：
\[ e_r =K f(\theta)(\frac{u_{rel}}{u_{ref}})^n \]
式中：

• \(K\)为材料相关阐述，对于普通钢材，其默认参数为2.0e-9
• \(f(\theta)\)为粒子冲