基于DPM模型的高速列车转向架区域风雪流数值模拟.pdf

第 20 卷 第 8 期2023 年 8 月铁道科学与工程学报Journal of Railway Science and EngineeringVolume 20 Number 8August 2023基于DPM模型的高速列车转向架区域风雪流数值模拟谢菲1，高广军2，何侃2，陈功3(1.湖南铁道职业技术学院，湖南 株洲 412001；2.中南大学 交通运输工程学院，湖南 长沙 410075；3.中车株洲机车电力机车有限公司，湖南 株洲 412001)摘要：为分析高速列车转向架区域积雪结冰分布规律，采用三维非定常雷诺平均(unsteady Reynolds-averaged Navier-Stokes,URANS)数值模拟方法及离散相模型(discrete phase model,DPM)研究高速列车转向架区域的压力分布、流场流线、雪粒子运动轨迹及积雪堆积情况，并通过对比数值计算与风洞试验结果，验证数值计算方法的准确性。数值计算结果表明：转向架区域由于受气流冲刷，轮对、电机、制动钳夹等突出部件在迎风侧出现明显正压，背风侧则为负压；空气流经转向架迎风侧时会向上扬起，进入转向架的轮对、电机、制动钳夹等关键部位，空气流经转向架背风侧时会冲刷转向架区域的设备舱后端板；空气中的雪粒子不仅具有随动性，还具有惯性，会随气流运动到转向架轮对、电机、制动钳夹等关键部位，并在惯性作用下黏附在其上；积雪堆积主要集中在轮对、电机、制动钳夹以及端板等部件底部迎风区域，而在转向架顶部和背风区域堆积较少。从而得出结论：风场的数值计算结果与风洞试验观测结果具有较高一致性，适用于转向架区域积雪结冰研究；在转向架迎风侧，雪粒子对于气流的跟随性较强，雪粒子的运动轨迹与流场流线分布相符，随气流冲刷至转向架，容易形成积雪结冰；在转向架背风侧以及转向架上方区域，雪粒子对于绕流的跟随程度减弱，不易形成积雪结冰。关键词：DPM模型；数值模拟；高速列车；转向架；积雪中图分类号：U238 文献标志码：A 开放科学(资源服务)标识码(OSID)文章编号：1672-7029（2023）08-2814-10Numerical simulation of the flow and snow particles in the bogie area of a high-speed train based on DPM modelXIE Fei1,GAO Guangjun2,HE Kan2,CHEN Gong3(1.Hunan Railway Professional Technology College,Zhuzhou 412001,China;2.School of Traffic&Transportation Engineering,Central South University,Changsha 410075,China;3.CRRC Zhuzhou Locomotive Co.,Ltd.,Zhuzhou 412001,China)Abstract:In order to analyze the distribution law of snow and ice in the bogie area of high-speed trains，the pressure distribution,the flow field streamline,the trajectory of snow particles and the snow accumulation in the 收稿日期：2022-08-30基金项目：2022年度湖南省教育厅科学研究项目(22C1111)；2022年度湖南铁道职业技术学院校级课题(KJ202203)通信作者：何侃(1991)，男，湖南长沙人，讲师，博士，从事地面交通空气动力学研究；E-mail：DOI:10.19713/ki.43-1423/u.T20221688第 8 期谢菲，等：基于DPM模型的高速列车转向架区域风雪流数值模拟bogie area of a high-speed train were studied by using a combination of three-dimensional unsteady Reynolds-averaged Navier-Stokes(URANS)numerical simulation method and discrete phase model(DPM).The numerical calculation method was verified by comparing the numerical calculation results with the wind tunnel test results.Numerical results show that the wheels,the motor,the braking caliper and other protruding parts have obvious positive pressure on the windward side and negative pressure on the leeward side due to the air flow scouring.In addition,the air rises and impinges the wheels,the motor,the braking caliper and other key parts of the bogie when the flow through the windward side of the bogie.Then,the flow impinges the equipment compartment rear end panel of the bogie when the flow through the leeward side.The snow particles in the air have both inertia and randomness,moving to the key parts of the bogie with the airflow and adhering to these parts,which results from the action of inertia.Snow accumulation is mainly concentrated in the windward area at the bottom of the wheels,the motor,the braking caliper,the panel,and so on,while it is less accumulated at the top and leeward area of the bogie.It is concluded that the numerical results are in good agreement with the wind tunnel test results,and are suitable for the study of snow accumulation and freezing in the bogie area.On the windward side of the bogie,the snow particles have a strong follow-up to the air flow.The movement trajectory of the snow particles is consistent with the streamline distribution of the flow field.It is easy to form snow and ice when the air flow scours the bogie,and in the leeward side of the bogie and the area above the bogie,the following degree of snow particles to the surrounding flow is weakened,so it is not easy to form accumulation.Key words:DPM model;numerical simulation;high-speed train;bogies;snow accumulation 高速列车在寒冷的冬季运行时，轨道上的积雪会受到列车风影响，在列车的底部和尾部飞扬1。此外，列车转向架上的电机、制动钳夹等设备在运行时会产生热能。当空气中运动的雪粒子接触这些发热设备时，将会融化为水，进一步在高寒低温条件下凝结成冰。因此，列车在长期运行过程中，积雪会在转向架区域持续堆积，最终形成厚重的冰雪混合物，如图1所示。这一现象将增加转向架配重，危害转向架各关键部件的工作性能，进而对列车行车安全产生极大隐患。针对列车积雪结冰问题，JEMT2采用对环境无害的加热丙二醇作为材料加入除冰设备中，来改善瑞典等北欧国家的轨道积雪问题。PAULUKUHN等3针对俄罗斯冬季低温、多雪环境，为保证列车能够经得起无法避免的冰雪侵害，在除冰后的列车底部喷射水和丙烯类混合液，从而减少结冰。BETTEZ4针对日本铁路线的列车运行对高效除冰的要求，采用喷洒器作为除冰装置，有效解决了因积雪结冰而导致列车延误的问题。综上所述，来自于冬季处于高寒环境国家的研究人员提出了一系列的列车防积雪应对策略。然而，针对转向架区域的积雪结冰研究鲜有文献提及。高广军等58采用化学发热材料或物理结构改进等方式来改善转向架区域的积雪情况。丁叁叁等911均对列车转向架区域积雪结冰问题开展了流场数值模拟研究。然而，这些研究仅从流场进风量的角度间接预测转向架区域的积雪情况，并未考虑雪粒子在空气中的分布。雪粒子虽对流场运动具有一定跟随性，但气流出现扰动或者流速降低时，雪粒子的运动惯性将使得其对流场的跟随性减弱。在转向架区域，由于复杂几何模型影响，气流较为紊乱，存在大量湍流结构，且流速相对车身周围流场较低。与真实情况相比，仅从流场的角度考虑转向架区域积雪分布有失偏颇。因此，本文试图从风雪两相流的角度探索转向架区域积雪分布情况。关于风雪两相流，FUJII等12针对日本铁路积雪情况进行研究，采用实验的方式分析了雪粒子运动的临界风速与雪密度及含水量之间的关系，提出雪粒子初始密度范围为36132 kg/m3。TOMINAGA 等13采用数值仿真对雪粒子的堆积、侵蚀等情况进行模拟，并通过实验验证得出雪粒子的密度为100 kg/m3时数值仿真与实验结果具有良好的一致2815铁 道 科 学 与 工 程 学 报2023 年 8月性。桑嘉宾14建立风雪流的数值模型，研究了混合粒径，得出雪粒子直径的影响因素，指出雪粒子直径一般分布在 00.6 mm 之间，且主要集中在0.2 mm附近。莫华美15针对屋面积雪分布建立了风雪流数值模型，利用该数值模型对屋面上的积雪分布规律进行数值模拟，指出雪粒子直径主要分布在 00.4 mm 之间。TOMINAGA 等16采用CFD数值模拟仿真的方法模拟雪粒子的堆积深度，并且采用风洞试验验证了模拟方法的可行性。BEYERS等17采用数值仿真模拟了雪粒子在建筑物周围逐步形成堆积的过程，有助于在建筑设计时对潜在的积雪情况进行处理。SERINE等18将两维空间上的雪粒子数值模拟方法和风洞试验对比，并将其推广到三维空间上，从而确定雪粒子数值仿真的可行性。张强等19基于欧拉两相流理论和扩展传热模型，提出一种模拟飞机机翼积雪结冰的方法，并在实验中得到验证。GORDON 等2021针对严寒天气的积雪现象，对雪粒子的形状、大小、速度及密度等进行研究，得出雪粒子的相关物理性质。SATO等22针对在强风情况下的雪粒子破碎问题，采用风洞试验对积雪的密度进行研究，并在实地观测结果中验证其合理性。ZHOU等23采用数值模拟技术对机场航站楼屋顶的雪粒子现象进行模拟，并用风洞试验对比验证数值模拟技术的可靠性，从而得出湍流模型、风速及风向等对积雪的影响。BEYERS等24采用数值方法模拟立方体周围的积雪情况，并通过实验进行对比验证。韩俊臣等25采用三维非定常雷诺时均模拟高速动车组转向架区域流场和雪粒子分布情况。娄振26采用准瞬态的方法模拟高速列车拖车转向架区域的积冰过程，研究转向架区域的积冰分布特征。上述文献表明，采用风雪两相流数值模拟仿方法可以再现积雪分布实际情况，具有较高的可行性。因此，本文采用数值模拟方法，复现列车转向架区域的流场流动情况，结合DPM模型预测了雪粒子在空气中的运动轨迹，揭示了雪粒子在转向架上的分布情况，并从雪粒子随流场运动的角度描述了转向架区域积雪形成的原因，为高速列车转向架区域积雪结冰研究提供参考。1 研究方法1.1物理模型、计算域及边界条件设置本文基于转向架区域流场，建立了如图2所示的带有单个转向架区域的高速列车车体模型。考虑到实际运行中，头车的结构和形状与其他车不一致，不具有普遍性。因此，选择中间车型作为研究对象，在实践中能更好地估计一般转向架上的积雪结冰情况。该简化模型忽略了列车头、尾车等流线型结构对真实列车宏观流场产生影响的因素，进而着眼于转向架区域。车体的高度为3.6 m，宽度为3.2 m，长度为17.6 m。转向架区域距离车体迎风侧端面距离为6.9 m，使得车体底部流场在进入转向架区域之前得到充分发展，符合实际运行的列车转向架区域流场情况。同时，转向架区域距离车体背风侧端面同样为6.9 m，使得转向架区域的尾流得到充分发展，避免尾流对转向架区域流场产生干扰。此外，为真实模拟高速列车运行时车体底部的流场结构，增加了离地高度为0.2 m的钢轨模型。图1转向架区域冰雪混合物堆积Fig.1Accumulation of the ice and snow mixture in the bogie area单位：m图2车体模型Fig.2Geometry of the train2816第 8 期谢菲，等：基于DPM模型的高速列车转向架区域风雪流数值模拟为模拟车体周围流场，建立了长 60 m，宽20 m，高20 m的计算区域，如图3所示。车体上游处(平面 ABCD)设置为速度入口，给定速度为55.56 m/s，使得雷诺数为 1.325107。车体下游处(平面 EFGH)设置为压力出口，给定压力为 0 Pa。车体两侧(平面 ADHE 和平面 BCGF)及顶面(平面ABFE)均设置为固定壁面，模拟风洞环境。地面(平面DCGH)设置为滑移地面，给定与车速相同的滑移速度。车体及转向架壁面设置为固定壁面。车体距离气流入口15 m，车体两侧与计算域侧边界距离等宽。其中，模型的阻塞比约为2.8%，少于5%，满足一般钝体车辆计算需求27。计算时间步长设置为0.001 s，使得库朗数平均值小于1。1.2网格设置本文数值计算中采用六面体主导的网格。计算域采用开源计算流体力学分析软件OpenFoam进行 SnappyHex 网格离散，网格单元总数 2 000 万。网格采用切割单元(cut-cell)方式，以级数递增方式从计算域边界向车体壁面逐层加密。车体壁面网格尺寸为0.02 m，转向架壁面网格尺寸为0.01 m。为准确模拟边界层流场，在车体和转向架壁面上布置了8层附面层网格，近壁面第1层网格高度为0.6 mm，使得壁面无量纲化尺寸y+处于30300之间。图4展示了转向架表面网格以及车身周围计算域剖面网格。为了明确本文数值模拟不受网格密度的影响，基于相同的壁面法向y+，在流向和展向使用了不同的网格尺度，并采用了中等密度网格和精细密度网格进行对比。其中，中等网格和精细网格数量分别是2 000万和4 000万。如图5所示，为从车体底面中心线的2个网格获得的时均压力系数Cp。在图中可以发现中等网格和精细网格的压力系数基本保持一致，说明本文采用的中等网格已经达到计算精度的收敛区间，具有较好的准确性。单位：m图3计算区域Fig.3Domain for simulation(a)转向架表面；(b)计算域纵截面；(c)计算域横截面图4数值计算网格Fig.4Mesh for simulation图5车体底面中心线压力系数对比Fig.5Time-averaged Cp along the centerline of the bottom surface of the train body2817铁 道 科 学 与 工 程 学 报2023 年 8月1.3流场数值模拟方法针对高速列车周围流场，HE等28指出采用大涡模拟(LES)能较好地捕捉湍流特性，但采用改进的延迟分离涡(IDDES)或者非稳态雷诺平均(URANS)方法同样能较为准确地模拟流场速度及压力分布。考虑到LES或IDDES对于网格分辨率要求较高，不适应于较高雷诺数条件下的数值仿真。此外，本文研究内容除了对流场运动方程求解之外，还需要额外求解雪粒子运动方程，计算量较大。因此，本文采用对网格尺度要求更低、在工程尺度得到广泛认可及应用的URANS方法对流场部分进行求解。Realizable k-湍流模型能较好地捕捉流动分离，保持雷诺应力与真实湍流一致，使计算结果更符合真实性。因此，为适配壁面网格尺寸，选用了Realizable k-湍流模型，其控制方程见参考文献29。1.4雪粒子数值模拟方法风吹雪是一种气固两相流，根据现代多相流计算理论，计算模型主要有欧拉两流体模型、混合模型以及拉格朗日粒子轨迹模型等30。拉格朗日粒子模型是另一种气固两相流计算方法，目前风吹雪数值模拟中应用较少。拉格朗日模型的优点是物理概念明确、结果直观，可计入粒子的轨道经历效应和历史效应，能得到粒子的详细运动信息31。因此，本文采用拉格朗日随机模型方法，对转向架区域的积雪分布进行数值模拟。在本文中，通过对粒子的力平衡方程进行积分，来预测离散相粒子的运动轨迹。其中，作用在粒子上的外力主要包括气动力、重力、曳力等32。除此之外，其余力对粒子运动影响相对较小，在分析中可忽略不计30，其运动方程见式(1)：dupdt=FD(u-up)+g(p-)p(1)其中：FD(u-up)为粒子单位质量受到的气动力，u为 流 体 相 速 度 矢 量，up为 粒 子 相 速 度，FD=18pd2pCDRe24，为动力学黏度，dp为粒子直径，p为粒子密度；为空气密度；g为惯性加速度矢量；Re为相对雷诺数，Re=dp|up-u；CD为阻力系数。对方程(1)积分就得到了颗粒轨道上每一个位置上的颗粒速度，见式(2)：dxdt=up(2)其中：x为粒子位移矢量，t为时间。粒子的轨迹是在离散的时间步长上逐步进行积分运算得到的，对式(2)进行积分，就可以得到粒子的轨迹。对于车体、轮对等不易积雪部件，设定雪粒子边界为反射条件；对于构架、电机、端板、制动钳夹等易积雪部件，设定雪粒子边界为捕捉条件，认为转向架区域雪粒子全部被捕捉，以此来考虑积雪最严重情况。其中，雪粒子密度为100 kg/m3，雪粒子直径采用均匀分布为 0.2 mm，雪 粒 子 材 料 为 water-liquid，列 车 运 行 速 度 为200 km/h。此外，时间步长t由库朗数Co最大值控制，定义为Co=vt/x。其中，x是网格尺寸，v是网格内的流体速度。此外，隐式算法中的库朗数范围为 0100。如图 6 所示，在本文数值模拟计算中，库朗数分布为011.85，且绝大部分集中在0和1之间，基本满足计算要求。当t=2.00 s时，气流多次通过该区域(风速55.56 m/s)，转向架中的风雪流场得到相对稳定和充分发展。并且，基于课题组多年计算转向架积雪的经验，2 s的时间足以使得风夹带的雪粒子数量在转向架区域达到相对稳定。因此，本文的计算时间设置为2 s，时间步长为0.001 s，共计2 000步。2 数值模拟方法验证关于本文中流场数值模拟方法的验证请参考图6库朗数等值线Fig.6Contours of the courant number2818第 8 期谢菲，等：基于DPM模型的高速列车转向架区域风雪流数值模拟文献33。本文仅针对数值模拟中雪粒子在转向架区域的分布情况进行验证分析，开展离散相风洞试验。该试验在中南大学高速列车研究中心的风洞实验室高速段进行。实验段的尺寸为 3.4 m1.0 m0.8 m(长宽高)。由于真实雪粒子存在难制备、难保存等问题，试验中采用了与雪粒子密度相似的木屑粒子(密度约为100 kg/m3)作为替代物。木屑均经过筛网严格过滤，挑选出与雪粒子直径(0.2 mm)大小类似的粒子进行试验。试验段主要安装的装置包括高速列车原始转向架模型、高速列车原始车体模型、固定装置，其中固定装置又包括支撑杆、支撑座、分流板以及安装架等构成，如图7所示。本次两相流风洞实验需要模拟雪粒子在转向架区域各关键部件的黏附情况，为此，在转向架和车体模型表面涂抹蜂蜜以黏附吹过的木屑。如图8所示，为涂刷蜂蜜后的转向架以及转向架各个细节部位的示意图。为验证数值模拟方法的可靠性，且保障计算与试验的一致性，依据风洞试验室和试验模型的大小，建立了和风洞试验模型完全相同的数值模拟模型，且计算入口风速与风洞试验自由流风速一致。网格离散原则及数值模拟方法与本文第1节描述一致，由于简化了车体模型，网格总数为1 400万。图9展示了风洞试验和仿真计算的转向架粒子堆积情况。从图9可以看出，风洞试验中木屑分布特点与数值仿真计算的积雪分布规律一致，粒子的堆积均产生在转向架各部件的迎风侧。因此，本文采用的数值仿真计算方法能较为准确地模拟雪粒子在转向架区域的分布情况。3 数值仿真结果分析本节采用第1节所述流场及雪粒子数值模拟方图7模型的安装与固定Fig.7Installation and fixation of the model图8涂刷蜂蜜后的转向架Fig.8Bogie painted with honey2819铁 道 科 学 与 工 程 学 报2023 年 8月法，针对转向架区域进行数值模拟仿真，得出转向架区域压力分布、流线分布、雪粒子分布及转向架表面积雪分布情况。3.1转向架区域压力分布针对转向架区域主要部件的积雪结冰现象进行分析，在转向架区域作纵截面空间切片，如图10 所示，切片 1 和 4 在转向架轮对、制动钳夹位置，切片2和3在转向架电机位置。如图11所示，展示了转向架区域时均压力云图，来流方向为从左至右。由于切片1和4较为对称，为简略表达仅展示切片1，2和3上的时均压力分布。从图11可以看出，轮对、制动钳夹、电机等突出部件在迎风侧由于受气流冲刷出现明显正压，而在这些部件的背风侧则为负压。转向架区域的压力分布体现了复杂几何结构对流场分布具有较大影响。因此，在进一步探讨雪粒子轨迹之前，需对流线分布进行分析。3.2转向架区域流线分布图12展示了转向架区域流线分布，来流方向为从左至右。从图12可以看出，车体底部气流在转向架区域发生了偏转，部分气流在流经转向架前端时向上扬起进入转向架区域，并在转向架区域各个部件处形成多处漩涡，说明转向架区域流场具有较高湍流度。此外，气流进入转向架区域之后受复杂几何结构影响，速度明显降低，尤其在转向架轮对、制动钳夹以及电机等部位的流速接近为0。当流速降低时，雪粒子跟随性降低，当流速低于雪粒子起动风速时，重力占雪粒子所受合力比重增大。因此，当雪粒子随气流上扬进入转向架区域时，更易在低速区域附着至转向架各部件上形成堆积。值得注意的是，在转向架后端，部分气流在设备舱端板受阻，再次回流入转向架区域。(a)风洞实验积雪情况；(b)数值仿真积雪情况图9转向架底部表面积雪结果对比Fig.9Comparison of snow results on the bottom surface of bogie图10转向架区域切片位置Fig.10Position of slices in the bogie area1轮对；2制动钳夹；3电机(a)切片1；(b)切片2；(c)切片3图11转向架区域压力分布Fig.11Pressure distribution in the bogie area2820第 8 期谢菲，等：基于DPM模型的高速列车转向架区域风雪流数值模拟3.3转向架区域雪粒子分布通过对转向架区域净风场流场情况的分析，已经初步判断出转向架区易形成雪粒子堆积区域。为进一步直观体现雪粒子在转向架区域的分布情况，图13展示了雪粒子运动的快照，体现了雪粒子的运动历程。从图 13 可以看出，在 t=0.12 s 时刻，雪粒子进入到转向架区域的前部；在t=0.16 s时刻，雪粒子进入到转向架区域的中部，且集中在高度较低区域，对制动钳夹、电机等发热元件的冲刷不明显；在t=0.20 s时刻，雪粒子经过整个转向架区域，更多的雪粒子进入到转向架的内部；当时间到达t=2.00 s时刻，经过更长物理时间后，大量雪粒子进入转向架区域。此后，雪粒子分布趋于稳定。整体而言，在列车运行过程中，大量雪粒子跟随气流从列车车体底部流入，尽管大多数雪粒子到达并且通过了转向架区域，仍然有许多雪粒子随着上扬的空气进入了转向架区域。在转向架区域前方，雪粒子跟随气流运动，具有一定的跟随性；在转向架区域后方，雪粒子在运动惯性作用下和气流发生偏离。进入转向架区域的雪粒子流速降低、跟随性减弱，因此被转向架区域各部件捕捉。随着时间推移，雪粒子将会持续进入转向架区域，从而在各部件上产生堆积。3.4转向架表面积雪分布如图14所示，为展示转向架区域主要积雪部位，本文显示了转向架区域积雪浓度大于 11012 kg/m3以上时的积雪覆盖情况，并采用速度渲染显示冲击到转向架区域时雪粒子的速度情况。图14中，来流方向为从左至右。可以看出，在转向架底部区域产生大量积雪，且积雪更易堆积在转向架各部件迎风侧，而在各部件背风区域积雪(a)t=0.12 s；(b)t=0.16 s；(c)t=0.20 s；(d)t=2.00 s图13雪粒子运动时间序列快照Fig.13Snap shots of the movement of snow particles(a)切片1；(b)切片2；(c)切片3图12转向架区域流线分布Fig.12Streamline in the bogie area(a)底视图；(b)俯视图图14转向架表面积雪分布Fig.14Snow accumulation on the surface of the bogie2821铁 道 科 学 与 工 程 学 报2023 年 8月相对较少。此外，积雪主要集中在电机、制动钳夹、轮对以及端板等部件底部迎风区域，而在转向架部件顶部和背风区域堆积相对较少。4 结论1)数值计算结果与风洞试验观测结果具有较高一致性，说明本文数值模拟方法准确可靠，适用于转向架区域积雪结冰研究。2)空气流经转向架迎风侧时会向上扬起，进入转向架的轮对、电机、制动钳夹等关键部位，空气流经转向架背风侧时会冲刷转向架区域的设备舱后端板。3)雪粒子运动轨迹和积雪堆积主要集中在轮对、电机、制动钳夹以及端板等部件底部迎风区域，而在转向架顶部和背风区域堆积较少。本文的数值计算方法可适用于模拟高速列车在风雪天气下运行时转向架区域积雪结冰情况，对于高寒动车组防积雪结冰研究具有一定参考价值。参考文献：1KLOOW L.High-speed train operation in winter climateR.KTH Railway Group and Transail,2011.2JEMT T.De-icing solutionJ.International Railway Journal,2009,49(1):2425.3PAULUKUHN L,WU Xinmin.The low temperatures technology concepts and operational experience in Russian high speed train velaro RUS J.Foreign Rolling Stock,2012,49(3):1619.4BETTEZ M.Winter technologies for high speed railD.Sweden:Norwegian University of Science and Technology,2011.5高广军,汪馗,赵世越,等.一种动车组车载转向架防冰涂层:CN111234587AP.20200605.GAO Guangjun,WANG Kui,ZHAO Shiyue,et al.Anti-icing coating for vehicle-mounted bogie of motor train unit:CN111234587AP.20200605.6余志壮.携热微砂喷射式高速列车转向架快速除积雪结冰方法:CN109131433AP.20200131.YU Zhizhuang.Heat-carrying micro-sand jet type high-speed train bogie rapid snow removing and deicing system and method:CN109131433AP.20200131.7王石,薛源,周业明,等.一种轨道车辆用车轮罩及轨道车辆:CN211642186UP.20201009.WANG Shi,XUE Yuan,ZHOU Yeming,et al.Wheel cover for railway vehicle and railway vehicle:CN211642186UP.20201009.8余以正,杨明智,滕万秀,等.车辆转向架前端加装弧形导流槽对转向架积雪结冰情况的影响J.城市轨道交通研究,2020,23(2):4954.YU Yizheng,YANG Mingzhi,TENG Wanxiu,et al.Influence of installing arc air deflector at the vehicle bogie front end on the situation of snow and ice accumulation in bogie zoneJ.Urban Mass Transit,2020,23(2):4954.9丁叁叁,田爱琴,董天韵,等.端面下斜导流板对高速列车转向架防积雪性能的影响J.中南大学学报(自然科学版),2016,47(4):14001405.DING Sansan,TIAN Aiqin,DONG Tianyun,et al.Influence of inclined guiding plate on anti-snow performance of high-speed train bogieJ.Journal of Central South University(Science and Technology),2016,47(4):14001405.10 苗秀娟,何侃.高速列车转向架区域积雪形成原因及防积雪研究J.中南大学学报(自然科学版),2018,49(3):756763.MIAO Xiujuan,HE Kan.Cause analysis of snow packing in high-speed trains bogie regions and anti-snow packing designJ.Journal of Central South University(Science and Technology),2018,49(3):756763.11 冯永华,黄照伟,张琰,等.高寒动车组转向架区域积雪结冰数值仿真研究J.铁道科学与工程学报,2017,14(3):437444.FENG Yonghua,HUANG Zhaowei,ZHANG Yan,et al.Research of numerical simulation of the snow icy phenomenon of the high-speed train bogie areaJ.Journal of Railway Science and Engineering,2017,14(3):437444.12 FUJII T,KAWASHIMA K,IIKURA S,et al.Preventive measures against snow for high-speed train operation in JapanC/Cold Regions Engineering.Anchorage,Alaska,USA.Reston,VA:American Society of Civil Engineers,2002:448459.13 TOMINAGA Y,MOCHIDA A,OKAZE T,et al.Development of a system for predicting snow distribution in built-up environments:combining a mesoscale meteorological model and a CFD modelJ.Journal of Wind Engineering and Industrial Aerodynamics,2011,99(4):460468.14 桑嘉宾.混合粒径风雪流运动的数值模拟D.兰州:兰2822第 8 期谢菲，等：基于DPM模型的高速列车转向架区域风雪流数值模拟州大学,2012.SANG Jiabin.Numerical simulations of blowing snow with mixed diameterD.Lanzhou:Lanzhou University,2012.15 莫华美.典型屋面积雪分布的数值模拟与实测研究D.哈尔滨:哈尔滨工业大学,2011.MO Huamei.Numerical simulation and experimental study of snow distribution on typical roofsD.Harbin:Harbin Institute of Technology,2011.16 TOMINAGA Y,OKAZE T,MOCHIDA A.CFD modeling of snowdrift around a building:an overview of models and evaluation of a new approachJ.Building and Environment,2011,46(4):899910.17 BEYERS M,WAECHTER B.Modeling transient snowdrift development around complex three-dimensional structures J.Journal of Wind Engineering and Industrial Aerodynamics,2008,96(10/11):1603 1615.18 SERINE A,SHIMURA M,MARUOKA A,et al.The numerical simulation of snowdrift around a buildingJ.International Journal of Computational Fluid Dynamics,1999,12(3/4):249255.19 张强,曹义华,钟国.飞机机翼表面霜冰的三维数值模拟J.航空动力学报,2010,25(6):13031309.ZHANG Qiang,CAO Yihua,ZHONG Guo.Three-dimensional numerical simulation of rime ice accretions on an aircraft wing