高速列车气动分析：为何0.31马赫场景下不可压缩流体模型仍适用？

1次阅读

共计 1513 个字符，预计需要花费 4 分钟才能阅读完成。

在高速列车设计中，气动性能分析至关重要。当列车速度达到 0.31 马赫（约 380km/h）时，许多工程师会疑惑：是否需要采用可压缩流体模型？本文将系统性地论证在此速度级下，不可压缩流体模型仍然适用的理论基础，并通过实际案例分析其工程适用性。

不可压缩流体模型的基本假设是流体密度变化可忽略不计。判断这一假设是否成立，主要考察两个无量纲参数：

雷诺数（Re）：
$$Re = \frac{\rho v L}{\mu}$$
其中 $\rho$ 为空气密度（~1.225kg/m³），$v$ 为车速（106m/s @380km/h），$L$ 为特征长度（取车头长度 10m），$\mu$ 为动力粘度（1.789×10⁻⁵Pa·s）。计算得 Re≈7.3×10⁷，表明流动处于高度湍流状态。
马赫数（Ma）：
$$Ma = \frac{v}{c}$$
其中 $c$ 为声速（340m/s）。0.31 马赫时，根据等熵流动关系，密度变化仅为：
$$\frac{\Delta\rho}{\rho} \approx \frac{Ma^2}{4} = 2.4\%$$
这一变化量在工程误差允许范围内（通常 <5% 可视为不可压缩）。

我们采用相同网格（2000 万单元）对 CR400AF 列车模型进行对比计算：

压力分布差异 ：
在车头驻点区域，可压缩模型最大压力系数 Cp=1.02，不可压缩模型 Cp=1.00，相对误差 2%；
车顶过渡区压力梯度偏差 <1.5kPa

计算效率对比 ：

| 模型类型       | 计算耗时（h）| 内存占用（GB）|
|----------------|--------------|----------------|
| 可压缩（SST）| 38.6         | 124            |
| 不可压缩（k-ε）| 12.4         | 78             |

通过用户自定义函数实现马赫数动态校正：

#include "udf.h"
DEFINE_ADJUST(mach_correction, domain)
{
  Thread *t;
  face_t f;
  real Ma = RP_Get_Real("mach-number"); // 获取全局马赫数

  if (Ma > 0.3) {real beta = sqrt(1.0 - pow(Ma,2)); // 普朗特 - 格劳厄特修正

    begin_f_loop(f,t)
    {F_Cp(f,t) *= beta; // 修正压力系数
      F_Shear_Stress(f,t) /= beta; // 修正剪切应力
    }
    end_f_loop(f,t)
  }
}

（注：需在 Materials 面板中将空气密度设为 constant）

不同外形列车的临界马赫数存在差异：