沙尘环境下高速列车气动特性分析

文恒 于梦阁 盛旭高 殷硕

摘要：  为研究沙尘环境下高速列车明线运行时的气动特性，基于剪切应力传输模型SST k w双方程湍流模型和拉格朗日离散相模型，与无沙环境下的高速列车气动特性进行比较，计算分析不同沙粒浓度、不同车速下的高速列车气动特性。计算结果表明：沙尘环境下，当车速一定时，列车整车气动阻力、头车气动阻力、尾车气动阻力均随沙粒浓度增加而逐渐增大，且与沙粒浓度近似呈线性关系;对于气动升力，当车速一定时，头车气动升力绝对值随沙粒浓度的增加而增大，尾车气动升力随车速的增加而降低。该研究成果可为高速列车在沙尘环境中的运行安全提供理论参考。

关键词：  高速列车; 沙粒浓度; 气动特性; 气动阻力; 气动升力

中图分类号： U271.91; U270.1+1  文献标识码： A

近几年，随着科学技术的不断进步，国家加快铁路网建设，驱动区域经济地协调发展。兰新高铁穿越我国西部高寒风沙区域，由于兰新高铁线路的特殊性，列车途经的百里风区、三十里风区是内陆大风天气频发的地区之一，时常会引发大风灾害性气象[1 3] ，强风地区大多缺少植被覆盖，地表裸露，而风速又往往远远大于起沙风速，因此地表的沙粒在大风的作用下撞击动车组，使得高速列车的气动特性明显变差[4] 。沙尘环境属于多相流中的气 固两相流问题，多采用欧拉 欧拉模型和欧拉 拉格朗日模型进行模拟计算。欧拉 欧拉模型主要描述两相的运动，用于高浓度离散相的问题;欧拉 拉格朗日模型适用于离散相的体积分数在10%～12%以下的问题[5] ，而且关注离散相的运动轨迹。C.Paz等人[6] 采用欧拉 拉格朗日方法研究了高速列车在沙尘环境中的运行安全性;熊红兵等人[7] 研究了沙尘暴环境下高速列车运行时的气动特性;李田等人[8] 采用欧拉 欧拉方法研究了不同沙尘暴环境下高速列车的动力学性能;倪守隆[9] 采用欧拉 欧拉方法研究确定了高速列车在沙尘暴环境下运行的安全域;高琛光[10] 采用欧拉 拉格朗日方法研究了不同风速、不同风向、头尾中间车设备舱的流动特性。王洪涛等人[11] 通过实验研究了近地面的沙粒体积分数随高度的变化，实验结果表明，即使是风速很高，但在接近地面处，沙粒相所占的体积分数也只是在10-4 量级，在很靠近沙床表面时，沙子的浓度随风速变化不大，沙粒在空气中所占的体积分数在10-4 量级以下，远小于10%，因此本文选用欧拉 拉格朗日模型进行研究。经过以上分析，目前国内外关于高速列车在沙尘环境下的空气动力学研究正处于起步阶段，以往的研究所采用的沙粒浓度也相对较低，本文应用欧拉 拉格朗日方法建立沙尘环境下高速列车空气动力学计算模型，开展沙尘环境下高速列车的空气动力学计算，研究车速、沙粒浓度对高速列车空气动力学性能的影响。

1 计算模型

一般情况下，高速列车明线运行时，其外部流场为不可压缩的定常流[12 13] ，工程中的流场计算多是采用基于雷诺时均方法的湍流模型，因此，高速列车绕流流场控制方程为定常不可压缩雷诺时均Navier-Stokes方程，质量以及动量守恒方程分别为[14]

j xj =0 （1）

ρ  j i  xj =-  xi +  ij xj +  xj  -ρu′i u′j  -Mi  （2）

式中， j为气流速度在j方向的时间平均值;xj为j方向的坐标，ρ为气体密度;ui为i方向气流的速度;xi为i方向的坐标;p为压强;u′i 为i方向气体湍流脉动速度;ui 为气流速度在i方向的时间平均值;uj为j方向气流的速度;u′j 为j方向气体湍流脉动速度; 为气体压力的时间平均值; ij 为黏性应力张量的时间平均值;-ρu′i u′j  为雷诺应力;Mi为相间动量交换第i个分量。

本文研究高速列车在沙尘环境中运行所引起的空气流动属于湍流流动，进行数值计算时湍流模型多种多样，本文选取的剪切应力传输模型（shear stress transfer，SST）广泛应用于在高速列车的空气动力学计算中[15 18] ，具有比较高的精度，其表达式为

κ i xi =  xj Γκ κ xj +Gκ-Yκ （3）

ρ  ω i xi =  xj Γω ω xj +Gω-Yω+Dω （4）

式中，Γκ和Γω、Yκ和Yω分别为κ和ω的扩散率和湍流耗散项;Gκ为κ的湍流生成项;Gω为ω的湍流生成项;Dω为交叉扩散项。

应用欧拉 拉格朗日方法建模时，假设沙粒是球形的，且不会发生变形，只考虑空气阻力和重力的作用，沙粒之间的相互作用不予考虑，在欧拉 拉格朗日参考框架下，沙粒运动方程为

duΡi dt = 18μ ρΡ d2Ρ  cReΡ 24 ui -uΡi + gi ρΡ -ρ ρΡ  （5）

式中，uΡi 为沙粒速度第i个分量;t为时间;μ为动力黏度;ρΡ为沙粒密度;dΡ为沙粒粒径;c为沙粒曳力系数;ReΡ为相对雷诺数;gi為重力加速度第i个分量。

2 数值模型及边界条件设定

2.1 计算区域及边界条件

依照中国铁路某高速列车（china railway high-speed，CRH）的几何外形，建立高速列车三维几何模型，采用头车 中间车 尾车构成的3节车编组模型，其中头车26.5 m，中间车25 m，尾车26.5 m，整车78 m，宽

3.38 m，高3.7 m。对列车进行简化，保留挡风板，忽略转向架、门把手、受电弓等细部特征[19] 。本文CFD计算是模拟风洞吹风的方式，在计算域流场中，假设列车静止不动，在列车前侧计算域入口设置大小等于速，方向与车速相反的来流风速，沙尘入射面选为计算区域前端入口。连续相边界条件为：列车前侧计算区域设置为速度入口，列车后侧计算区域设置为压力出口，计算区域底面及列车表面采用无滑移壁面条件。离散相边界条件为：地面边界条件设置为trap，列车壁面边界条件设置为reflect，入口和出口边界条件均设置为escape。沙尘环境计算区域及边界条件如图1所示，列车鼻尖距离计算域左侧为110 m，距离计算域右侧为290 m。

[7] Xiong H B， Wen-Guang Y U， Chen D W， et al. Numerical study on the aerodynamic performance and safe running of high-speed trains in sandstorms[J]. Journal of Zhejiang University SCIENCE A， 2011， 12（12）： 971 978.

[8] 李田， 張继业， 张卫华. 沙尘暴环境下的高速列车运行安全分析[J]. 四川大学学报： 工程科学版， 2012， 44（S2）： 5 8.

[9] 倪守隆. 风雨与沙尘暴环境下列车运行安全性研究[D]. 大连： 大连交通大学， 2015.

[10] 高琛光. 风沙环境下高速列车设备舱流动特性数值模拟研究[D]. 兰州： 兰州交通大学， 2016.

[11] 王洪涛， 董治宝， 张晓航. 风沙流中沙粒浓度分布的实验研究[J]. 地球科学进展， 2004， 19（5）： 732 735.

[12] 黄成涛， 王立新. 风雨对飞机飞行安全性的影响[J]. 航空学报， 2010， 31（4）： 694 700.

[13] 李陶. 风沙环境下高速列车气动效应研究[D]. 兰州： 兰州交通大学， 2018.

[14] Wu Z L， Cao Y H. Numerical simulation of flow over an airfoil in heavy rain via a two-way coupled eulerian-Lagrangian approach[J]. International Journal of Multiphase Flow， 2015， 69： 81 92.

[15] Muoz-Paniagua J， García J. Aerodynamic drag optimization of a high-speed train[J]. Journal of Wind Engineering and Industrial Aerodynamics， 2020， 204： 104215 104220.

[16] 于梦阁， 李海庆， 刘加利， 等. 强风雨环境下高速列车空气动力学性能研究[J]. 机械工程学报， 2020， 56（4）： 185 192.

[17] 于梦阁， 李田， 张骞， 等. 强降雨环境下高速列车空气动力学性能[J]. 交通运输工程学报， 2019， 19（5）： 96 105.

[18] 周鹏， 李田， 张继业， 等. 真空管道超级列车激波簇结构研究[J]. 机械工程学报， 2020， 56（2）： 86 97.

[19] 李人宪， 刘杰， 戚振宕， 等. 明线会车压力波幅变化规律研究[J]. 机械工程学报， 2011， 47（4）： 125 130.

[20] 李凯崇， 蒋富强， 薛春晓， 等. 兰新铁路十三间房段的戈壁风沙流特征分析[J]. 铁道工程学报， 2010， 27（3）： 15 18.

Analysis of Aerodynamic Characteristics of High-Speed Trains in Sand and Dust Environment

WEN Heng， YU Mengge， SHENG Xugao， YIN Shuo

（School of Electromechanic Engineering， Qingdao University， Qingdao 266071， China）

Abstract：  In order to study the aerodynamic characteristics of high-speed trains in the open environment under sand and dust environment， based on the SST k w double equation turbulence model and Lagrangian discrete phase model， the aerodynamic characteristics of high-speed trains in a non-sand environment are compared， and different sand concentrations are calculated and analyzed. The calculation results show that： in a dust environment， when the speed of the trainis constant， the aerodynamic resistance of the entire train， the aerodynamic resistance of the first car， and the aerodynamic resistance of the tail car gradually increase with the increase of the sand concentration. It is approximately linearly related to the concentration of sand particles; for aerodynamic lift， when the vehicle speed is constant， the absolute value of the aerodynamic lift of the head car increases with the increase of the concentration of sand particles， and the aerodynamic lift of the tail car decreases with the increase of the speed of the car. The aerodynamic characteristics of high-speed trains deteriorate. The research results provide a reference for the safe operation of high-speed trains in dusty environments.

Key words： high-speed train; sand particle concentration; aerodynamic characteristics; aerodynamic resistance; aerodynamic lift