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

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

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

　　近几年，随着科学技术的不断进步，国家加快铁路网建设，驱动区域经济地协调发展。兰新高铁穿越我国西部高寒风沙区域，由于兰新高铁线路的特殊性，列车途经的百里风区、三十里风区是内陆大风天气频发的地区之一，时常会引发大风灾害性气象[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。

　　推荐阅读：工程科技专业论文好发表吗

期刊VIP网，您身边的高端学术顾问

文章名称： 沙尘环境下高速列车气动特性分析

文章地址： http://www.qikanvip.com/huanjingkexue/55041.html