旋转列车气流的大涡模拟概论
可逆转轮单通道流动大涡模拟
可逆转轮单通道流动大涡模拟
徐岚;崔桂香;许春晓;张兆顺;陈乃祥
【期刊名称】《水力发电学报》
【年(卷),期】2007(26)4
【摘要】为准确预测不可压复杂边界的湍流流动,本文应用高精度有限体积法对曲面边界湍流进行了大涡模拟。
空间离散采用有限体四阶紧致格式,时间推进采用四阶Runge-Kutta法,压力-速度耦合应用四阶紧致格式的动量插值,亚格子应力模式采用动态Smagorinsky模式,复杂边界的处理则应用了浸没边界法,成功地实现了弯槽湍流、NACA0012标准翼型绕流流动和可逆式水泵水轮机转轮内单流道流动的大涡模拟计算,所得结果与已有结果或实验结果吻合良好,表明该方法对于湍流大涡模拟方法在流体机械工程领域中的应用和发展具有重要意义。
【总页数】6页(P124-129)
【关键词】水力机械;可逆转轮;大涡模拟;有限体积法;浸没边界
【作者】徐岚;崔桂香;许春晓;张兆顺;陈乃祥
【作者单位】清华大学热能工程系;清华大学工程力学系
【正文语种】中文
【中图分类】TK72;TV131.2
【相关文献】
1.基于大涡模拟尾水管涡带模拟及流动特性研究 [J], 钟林涛
2.基于大涡模拟的抑制孔腔涡旋流动与脉动压力的流动控制方法研究 [J], DENG
Yu-qing;ZHANG Nan
3.混流式转轮中流场的大涡模拟 [J], 张昌兵;杨永全;鞠小明;桂林
4.大涡模拟研究两相燃烧的进展(Ⅱ)复杂气固流动和煤燃烧的大涡模拟(英文) [J], 周力行;胡瓅元
5.可逆转轮三维流动的涡动力学诊断研究 [J], 樊红刚;陈乃祥;杨琳
因版权原因,仅展示原文概要,查看原文内容请购买。
高速列车内气流组织的大涡模拟的开题报告
高速列车内气流组织的大涡模拟的开题报告
目的:
本文旨在通过大涡模拟(LES)模拟高速列车内气流组织,研究列车内空气流动
特性,为列车内空气质量的评估及改善提供理论基础。
背景:
近年来,随着高速列车的发展,越来越多的人选择高速列车作为出行方式。
然而,在高速列车内,由于车内人数众多,呼吸排放、食物烹饪、香烟等因素会使得车内空
气质量下降。
因此,对高速列车内气流组织进行研究,对于理解车内空气质量变化规律,为改善车内空气质量提供依据和方案。
方法:
采用计算流体力学(CFD)软件中的大涡模拟(LES)方法,模拟高速列车内气
流组织。
首先,选取高速列车内部为模拟区域,通过计算高速列车内的风压、温度、
湿度等参数,确定模拟边界条件。
然后,使用CFD软件建立高速列车内部模型,在模拟过程中考虑人员排放、空调系统、车窗开关等因素,模拟高速列车内部气体流动情况。
最后,对模拟结果进行分析并与实际观测结果进行对比。
预期结果:
通过大涡模拟,本研究预计获得高速列车内部不同区域的气流组织情况,揭示列车内流场特征,探讨不同因素对空气质量的影响。
同时,结果分析将为改善高速列车
内部空气质量提供参考,为制定车内通风方案、调整温度湿度等因素提供理论基础。
结论:
通过大涡模拟,本研究将为理解高速列车内气流组织、探究车内空气流动特性、评估空气质量等提供较为准确的数值模拟结果,并为改善列车内部空气质量提供理论
基础和实践指导。
旋转列车气流的大涡模拟
等 的 和加 密 的 3种 计 算 网格 。加 密 网 格 的 计 算 结 果 与 试 验 数 据 吻 合 较 好 。运 用 大 涡 模 拟 获 得 了 不 同 的 流 动 区域 。研 究 表 明 , 以柱 面形 式 支 撑 的风 挡 和车 下 复 杂 结 构 对 气 流 的速 度 有 很 大 的影 响 。在 合 适 的雷 诺 数 范 围
H a s n He da, t a . ( . ) s a mi e 1 U K.
Ab ta t sr c :Th l s r a o i h s e d t an wa n e t a e sn a g — d y s u a i n ( es i te m fa h g p e r i s i v si t d u ig l r e e d i l t p g m o LES ). e Th
m e h a e i f il g e m e ih t x rm e a t . Dif r ntfow e on e e o ane i g t s r n ary a r e ntw t he e pe i nt lda a fe e l r gi s w r bt i d usn he LES.
w a a e a w o Re no dsnu be s o 7 00 nd 94 0 0 b e he h i ft r n an t p e s m d tt y l m r f7 0 a 0 as d on t eghto hetai d iss e d. Thr e e
内 , 流 流 速 与列 车速 度 近 似 地 呈 线 性 关 系 。 气
关 键 词 : 车 ;旋 转 ; 流 ;大 涡模 拟 列 气
中 图 分 类 号 :U2 0 1 1 7 . 文 献 标 识 码 :B
挡风墙背风侧高速列车周围流场非定常特性大涡模拟研究
挡风墙背风侧高速列车周围流场非定常特性大涡模拟研究挡风墙背风侧高速列车周围流场非定常特性大涡模拟研究近年来,高速列车在交通运输领域起着越来越重要的作用,然而其在高速行驶时所面临的背风问题对行车安全和旅客舒适度带来了一系列挑战。
为了解决这一问题,科学家们开始采用计算流体力学(CFD)方法来研究列车周围的流场特性,进而改善列车的设计和操作。
本文旨在使用大涡模拟(LES)方法,对挡风墙背风侧高速列车周围的流场非定常特性进行研究。
通过模拟和分析,我们将深入探讨不同参数对流场的影响,从而为改善列车背风问题提供重要的理论依据。
首先,我们需要对列车及其周围流场进行建模。
为了简化问题,我们假设列车为一个理想化的长方体,并采用三维本构方程来描述流体运动。
列车周围的空气被假设为不可压缩、粘性和理想气体。
在数学方面,我们将采用标准的控制方程(Navier-Stokes方程),并应用LES方法来模拟非定常流体运动。
接下来,我们将根据实际情况设定流场的边界条件和初始条件。
边界条件包括列车表面的速度和压力分布,以及周围空气的入口速度和出口静压。
初始条件将决定流场在初始状态下的速度和压力分布。
在模拟过程中,我们将使用合适的数值计算方法来求解控制方程。
其中,我们将采用有限体积法来离散方程,并结合显式时间推进方法来求解非定常问题。
通过合理选择网格大小和时间步长,可以提高计算的准确性和效率。
当流场的计算结果得到后,我们将进行详细的结果分析和讨论。
首先,我们将关注流场中的速度和压力变化,了解列车背风侧产生的大涡结构以及其对列车尾部气流的影响。
然后,我们将通过比较不同参数条件下的计算结果,分析其对流场的影响。
最后,我们将根据模拟结果提出一些建议和改进方案,以减小列车背风问题带来的负面影响。
这些建议可能包括调整挡风墙的位置和形状,优化列车表面的流线型设计,以及改善列车尾部的气动性能等。
总之,本研究采用大涡模拟方法,对挡风墙背风侧高速列车周围的流场非定常特性进行了深入研究。
大型客机复杂可压缩流的大涡模拟主要研究方法
大型客机复杂可压缩流的大涡模拟主要研究方法一、大涡模拟基础1. 大涡模拟简介大涡模拟是一种将流场分解成小尺度湍流和大尺度湍流的方法。
在LES中,大尺度结构通过直接数值模拟来求解,而小尺度结构则通过子网格模型(sub-grid model)进行建模。
由于小尺度结构不再需要直接求解,因此可以使用更粗的网格来进行计算,从而减少计算量。
同时,LES还能够提供更加真实的湍流统计数据,如湍流强度、湍流长度等。
2. LES的优点和局限性与其他流体力学方法相比,LES有以下几个优点:(1)能够考虑湍流中的时间和空间尺度差异,提供更加真实的湍流信息;(2)计算结果对于网格的依赖性相对较小,使得计算可以在较粗的网格上进行;(3)LES能够模拟复杂流场,如湍流燃烧、多相流等。
虽然LES具有很多优点,但它也有一些局限性:(1)计算量较大,需要使用高性能计算机进行计算;(2)由于需要建立子网格模型,LES的结果可能受到模型误差的影响;(3)由于直接数值模拟只考虑了大尺度结构,因此对于小尺度结构的预测可能存在误差。
二、大涡模拟在大型客机流场研究中的应用1. 大涡模拟在飞行器气动力学研究中的应用大型客机的外形复杂,流场也非常复杂。
对于这样的流场,传统的计算流体力学方法可能无法准确地预测气动力学行为。
因此,大涡模拟成为研究大型客机流场的一种重要方法。
在大涡模拟中,通过将流场分解成大尺度结构和小尺度结构,可以更加准确地模拟大型客机流场中的湍流现象。
大涡模拟还能够提供更加真实的气动力学数据,如升阻比、气动力矩等。
这些数据对于飞机设计和优化非常重要。
2. 大涡模拟在飞行器噪声研究中的应用随着人们对噪声污染的关注度不断提高,飞机噪声研究也越来越受到关注。
大型客机飞行时产生的噪声主要来自于引擎和机翼表面的湍流。
由于湍流现象非常复杂,传统的计算流体力学方法无法准确地预测噪声的产生和传播。
因此,大涡模拟成为研究飞机噪声的一种重要方法。
通过大涡模拟,可以更加准确地模拟湍流现象,从而预测噪声的产生和传播方式。
超声速横侧射流混合特性的大涡模拟
超声速横侧射流混合特性的大涡模拟曹长敏;赵马杰;周涛涛;叶桃红【期刊名称】《中国科学技术大学学报》【年(卷),期】2015(000)008【摘要】采用大涡模拟方法研究了 Gamba 超声速燃烧室内的横侧射流流场中大尺度涡旋结构以及当地混合特性.超声速来流受到音速射流流体的阻碍,形成了复杂的激波和涡旋结构.由计算结果中的平均马赫数分布图可以清楚地看到激波结构,包括弧形激波、λ激波、桶状激波以及马赫盘;采用Q 准则表征三维涡旋结构,可以看到稳态的反向旋转涡对(CVP)、尾迹反向旋转涡对(TCVP)、马蹄涡以及非稳态的射流剪切层涡等结构;此外,由平均流线图可以看到,TCVP 结构与CVP 结构的旋转方向相反,不对称的 CVP 结构导致燃料质量分数展向分布不均匀.引入概率密度函数方法分析当地混合特性,结果表明射流近场混合主要发生在射流出口上游回流区以及桶状激波下侧和射流剪切层下侧的射流尾迹区内;射流远场混合分数的概率密度分布从β分布逐渐过渡为高斯分布.研究射流浓度衰减特征,结果表明氢气质量分数沿射流浓度最大迹线呈指数-0.7衰减.%Large eddy simulation was performed to study the large scale vertical structures and the mixing characteristic of supersonic transverse injection flow in the Gamba combustor . The supersonic flow is obstructed by the sonicjet ,resulting in very complex three dimensional shock waves and vertical structures .All of the shock waves ,including bow shock ,λ shock ,expansion shock ,barrel shock and Mach disk ,can be seen in the mean Mach number contour .The three‐dimensional vertical structures ,such as steady counter‐rotating vortex pairs (CVP) ,trailing counter‐rotating vortex pairs (TCVP) ,horseshoe vortices and unsteady jet shear layer vortices ,can be visualized by the iso‐surface of Q‐criterion . The TCVP structure rotates in the opposite direction of the CVP ,which can be characterized by the average streamline .Moreover , the asymmetric CVP structure leads to non‐uniform distribution of fuel mass fraction in span‐wise direction . The probability density functions (PDFs) were introduced to study the mixing characteristic in the near field ,far field and the recirculation regions in the upstream of the jet .The results show that significant mixing occurs in the recirculation regions ,the bottom of the barrel shock and the wake below the jet shear lay er .The PDFs of mixture fraction obey βdistribution in the near jet flow‐field and develop into Gauss distribution in the downstream .The study of the concentration decay shows that the H2 mass fraction obeys exponential decay along the maximum centerline with the index being about - 0.7 .【总页数】9页(P665-673)【作者】曹长敏;赵马杰;周涛涛;叶桃红【作者单位】中国科学技术大学热科学和能源工程系,安徽合肥 230027;中国科学技术大学热科学和能源工程系,安徽合肥 230027;中国科学技术大学热科学和能源工程系,安徽合肥 230027;中国科学技术大学热科学和能源工程系,安徽合肥 230027【正文语种】中文【中图分类】TK472+.6【相关文献】1.超声速高度欠膨胀冲击射流的大涡模拟 [J], 刘海;姚朝晖2.RBCC燃烧室超声速反应混合层特性的大涡模拟 [J], 魏祥庚;曹东刚;秦飞;吴继平3.超声速气流中液体横向脉冲射流一次破碎的大涡模拟 [J], 林森; 沈赤兵; 肖锋; 朱元昊4.激励射流在超声速流场中的混合特性研究 [J], 顾声龙;陈立红;顾洪斌;李飞;卫喆;张新宇5.三维超声速混合层内标量混合的大涡模拟 [J], 孙明波;范晓樯;梁剑寒;王振国因版权原因,仅展示原文概要,查看原文内容请购买。
大涡模拟简单介绍
《粘性流体力学》小论文题目:浅谈大涡模拟学生姓名:***学生学号:*********完成时间:2010/12/16浅谈大涡模拟丁普贤(中南大学,能源科学与工程学院,湖南省长沙市,410083)摘要:湍流流动是一种非常复杂的流动,数值模拟是研究湍流的主要手段,现有的湍流数值模拟的方法有三种:直接数值模拟、大涡模拟和雷诺平均模型。
本文主要是介绍大涡模拟,大涡模拟的思路是:直接数值模拟大尺度紊流运动,而利用亚格子模型模拟小尺度紊流运动对大尺度紊流运动的影响。
大涡模拟在计算时间和计算费用方面是优于直接数值模拟的,在信息完整性方面优于雷诺平均模型。
本文还介绍了对N-S方程过滤的过滤函数和一些广泛使用的亚格子模型,最后简单对一些大涡模拟的应用进行了阐述。
关键词:计算流体力学;湍流;大涡模拟;亚格子模型A simple study of Large Eddy SimulationDING Puxian(Central South University, School of Energy Science and Power Engineering, Changsha, Hunan,410083)Abstract:Turbulent flow is a very complex flow, and numerical simulation is the main means to study it. There are three numerical simulation methods: direct numerical simulation, large eddy simulation,Reynolds averaged Navier-Stokes method. Large eddy simulation (LES) is mainly introduced in this paper. The main idea of LES is that large eddies are resolved directly and the effect of the small eddies on the large eddies is modeled by subgrid scale model. Large eddy simulation calculation in computing time and cost is superior to direct numerical simulation, and obtain more information than Reynolds averaged Navier-Stokes method. The Navier-Stokes equations filtering filter function and some extensive use of the subgrid scale model are simply discussed in this paper. Finally, some simple applications of large eddy simulation are told.Key words:computational fluid dynamics; turbulence; large eddy simulation; subgrid scale model0 引言无论是在自然界还是在工程中,流体的流动很多都是湍流流动,例如,山中的流水,飞流直下的瀑布,飞机机翼旁边的气体流动,喷嘴的射流,炉内的气体流动等等。
列车进入隧道时产生涡流的数值模拟
文章编号:167320291(2007)0120063206列车进入隧道时产生涡流的数值模拟骆建军1,王梦恕1,高 波2(1.北京交通大学土木建筑工程学院,北京100044;2.西南交通大学土木学院,四川成都610031)摘 要:给出了列车穿越隧道时在隧道入口形成涡流场的数值模拟过程.控制方程为三维黏性、可压缩、等熵、非定常流的N-S (Navier-Stokes )方程,空间离散采用了中心有限体积法格式,时间采用预处理二阶精度多步后差分格式进行离散,对列车与隧道之间的相对移动采用移动网格技术处理.计算结果与国外的试验结果基本一致.研究结果表明,在整个涡流的形成过程中,可以将这个过程分为连续的4个阶段,即涡的形成、发展、传播及破坏阶段.研究结果还表明,涡受到列车头部形状的影响要比受到列车速度的影响要大;喷射流速度与列车速度有关系.关键词:列车;空气动力学;隧道;涡;数值模拟中图分类号:U45113 文献标志码:ANumerical Simulation of the V ortex Inducedby a T rain E ntering into a TunnelL UO Jian-j un 1,W A N G Meng-shu 1,GA O Bo2(1.School of Civil Engineering and Architecture ,Beijing Jiaotong University ,Beijing 100044,China ;2.S outhwest Jiaotong University Civil engineering Department ,ChengDu Si-chuang 610031,China )Abstract :A cell-centered finite volume method for numerical integration of N-S equations has been pre 2sented here for three-dimensional unsteady viscous compressible isentropic flow to simulate the forma 2tion of vortex produced by the train entering a tunnel.The numerical procedure is based on precondi 2tioned second-order accurate backward difference in time ,the relative motion between the train and tunnel treated with sliding mesh technique.The results of numerical simulation agree with that of for 2eign.The results show that the formation of the vortex constitutes four phases ,the pre-vortex phase ,the v ortex development phase ,the vortex convection phase and the v ortex breakdown phase.The duration of each of these phases has its own features ,and further more ,the effect of the train nose geometry is greater than the train speed ’s on the vortex ev olution.and the jet velocity is related to the train speed.K ey w ords :train ;aerodynamics ;tunnel ;vortex ;numerical simulation 列车进入隧道入口时伴随着两个重要的空气动力学现象,一是与列车运动方向相反的朝隧道入口运动的出口流,其次是朝列车前方运动的被压缩的空气流并形成压缩波.出口流与隧道内的压缩空气流有着密切的关系.对隧道入口形成的出口流和隧道内压缩波的研究是对由于隧道-列车这一体系所引起的空气动力学现象的较全面的了解.对隧道出口流的研究,国外是通过使用烟流发生装置和高速摄影仪进行可视化研究,没有进行数值模拟.国内对隧道空气动力学的研究,主要还集中在压缩波的研究,对隧道出口流的研究还没有涉及[1-4].本文作者在大型流体软件Fluent 的基础上进行了二次开发,收稿日期:2005209220基金项目:中国博士后科学基金资助项目(中博基字[2005]193)(CC06001);国家自然科学基金资助项目(50178060)作者简介:骆建军(1971-),男,湖南汨罗人,博士.em ail :jj -luo @第31卷第1期2007年2月 北 京 交 通 大 学 学 报JOURNAL OF BEI J IN G J IAO TON G UN IV ERSIT Y Vol.31No.1Feb.2007对列车进入隧道产生出口流进行三维数值模拟计算[5],对出口流现象进行了比较详细的研究,为后续隧道空气动力学的深入研究打下基础.1 数值模拟计算111 控制方程在笛卡尔坐标系下,对于空间任意控制体V,积分形式的三维可压缩性N-S方程为55tµW d V+κ[F-G]・d S=0(1)并且有W=ρρvxρvyρvzρe,F=ρνρνvx+p^iρνvy+p^jρνvz+p^kρνe+pν,G=τxiτyiτziτijv j+q(2)式中,ρ为气体的密度,v为气体的速度矢量,p为气体压力,e为单位体积的总能量,τ为黏性应力张量,q为热流通密度,ν为流体的运动粘滞系数,S 为面积,^i、^j、^k为单位方向向量.式(1)中共有8个未知数,5个方程,为了封闭上述方程组,引进方程p=ρR TH=e+p/ρp=(γ-1)[e-ρ(v2x+v2y+v2z)/2]h=C p T(3)式中,H为焓,R为气态常数,γ为比热比,T为气体温度,C p为比压.热传导系数κ由分子导热和湍流导热两部分组成,即κ=κc +κt=γRγ-1μcPr c+μtPr t(4)式中,κc为分子热传导系数,κt为代数湍流热传导系数;Pr c为分子Prantl常数,Pr t为湍流Prantl常数(在壁面处Pr t=0185);μc为分子黏性系数,由Sutherland式计算得到,μt为代数湍流黏性系数,采用Baldwin-Lomax提出的代数湍流黏性模型计算. 112 空间和时间离散对计算区域流场采用有限体积法进行离散.由于列车进入隧道时所诱发的气体速度比声速c小得多,由式(1)所表示的N-S方程在低马赫数下变得越来越“硬化”,在这种情况下会导致方程的解不能很快地收敛.采用时间导数预处理方法[6-7]对式(1)进行处理,这样就能按比例对c重新放大或缩小,以减小“硬化”,得到很好的收敛效果.为了能够给经过预处理后的方程提供一个有效的时间精度解,采用对偶时间步、多阶段格式并在式(1)中引进预处理拟时间微分项.在求解过程中,物理时间步长Δt根据所需要的非定常流的临时精度确定,而拟时间步长Δt′则根据稳定性条件确定.计算过程中,计算步数为374578.113 网格的划分由于列车与隧道之间是一种相对运动,当列车向前运动时,计算区域不断地发生变化,对于这种情况可以采用移动网格来处理.整个区域划分成不同的子区域并各自单独进行网格划分.在实际的流场中,由于列车上有列车车轮、电动机;在隧道底板上有道床、枕木和轨道,因此流场是很复杂的.列车底部以及隧道的地面都简化为平板.最后对整个计算区域进行网格划分后得到107420个单元.114 边界条件的确定及流场的初始化首先,列车的侧壁表面及隧道表面的法向速度为零,其次列车表面上压力梯度的法向分量为零.压力远场边界条件是基于黎曼不变量的无反射边界条件,计算区域内产生的压力波在该处不是反射波而是透射波.在边界面处,流体流出的一面的切向速度分量和熵值都通过插值函数从流场内部得到;而流体流入的一面,这些值则指定为自由流场的值.流场的初始化,首先在隧道入口处,列车进入隧道前可对列车周围的流场进行稳定流的计算,计算收敛时的流场即可作为非定场流场的初始条件.时间步长Δt=010005s,这个时间步长对于非定常流场来说足够小,能够满足计算的需要[8].115 方程求解对流体分析软件Fluent进行二次开发,并进行三维模拟.对列车进入隧道的全过程采用滑移网格法进行模拟.湍流方程采用标准k-ε模型.对三维可压缩性N-S方程,可采用耦合和分离两种方法进行求解.耦合解法即同时求解连续性方程、动量方程、能量方程;而分离求解法是分别求解连续性方程、动量方程和能量方程.本文作者采用耦合求解法,其一般过程如图1所示.116 计算模型的验证为了与英国Ueberlinger西部隧道现场量测实验结果进行比较,在本文的计算过程中,完全按照隧道及列车的实际尺寸进行模拟.隧道的横截面积为28m2,隧道长度为948m,列车横截面积为1013 m2,列车长度为122m,列车速度V=2213m/s,阻46北 京 交 通 大 学 学 报 第31卷图1 方程求解过程示意图Fig.1 Sketch of solution of equations塞比为01368.当列车进入隧道时,由于列车头部排挤其前方的气体,使其产生扰动向前传播.图2给出了隧道内距离洞口117m 处相对压力随着时间变化的值,可以看出实测值与数值模拟吻合得较好,从而可以通过该模型进行隧道出口流的计算研究.图2 隧道内压力的历时变化Fig.2 Curve of internal pressure in tunnel2 计算结果及分析图3为列车以20m/s 的速度进入隧道时,在隧道入口10m 处所形成的压力波示意图.图4表示涡流的形成过程,涡的发展顺序按图中标注的数字顺序进行.整个涡流的形成过程可分为4个阶段.图3 测点的压力变化Fig.3 Pressure curve of measure point(1)涡的形成阶段此阶段是压缩波的形成阶段,并且压缩波达到最大值,此时列车头部完全进入隧道,列车头部前方的气体一部分在压缩波压力作用下向隧道前方运动,另一部分气体则通过列车与隧道之间的环状空间向列车后方的隧道出口方向溢出隧道.随着列车头部进一步进入隧道,无论是列车前方的气体还是溢出隧道的气体,都得到加速,并存在一个最大的喷射流速度.当气体从隧道中溢出时,由于边界层的影响,在隧道内壁形成涡层,当涡层向后运动离开隧道口时,它形成涡线并立即卷曲形成涡,涡的直径逐渐增大并向周围扩展.(2)涡的发展阶段图5中,在隧道入口处距离隧道入口端墙d =50mm 的同一截面上布置4个不同的测点,设列车头部进入到隧道内z =z c 处时,其相对应的时间为t =t c (t =0表示列车头部到达隧道入口平面处).列车进入隧道的过程中,各个测点的速度如图6所示.速度的变化可以分成两个阶段.①t ≈3s 阶段,各测点处的速度逐渐增大,这个时间大于整个列车头部进入隧道所需要的时间.在这一阶段,在-019m <x <0范围内(列车与隧道之间的环状空间为118m ),喷射流速度大于其它点的速度.②t ≥3s 阶段.对于-0145m ≤x ≤0(测点4),它表示隧道壁面的边界层.在列车头部进入隧道后,出口处的喷射流速为常数,这个常速度用V jet 表示.图6中的喷射速度近似为V jet =1712m/s.在x =-1135m 和x =-1171m (即测点2和测点1)的速度-时间曲线上出现不连续的起伏,这些起伏表示在隧道外紊流边界层中紊流的特性.当t ≥3s 时,喷射流速度基本上保持为一个常数,随着时间的推移,喷射流速度稍微有点上升,这是由于列车头部断面积的不断增大和隧道边界层厚度的增加的缘故.在计算过程中,喷射流速度的增长不会超过列车的速度,即V jet <V .(3)涡的传播阶段这个阶段的特点是,在隧道入口平面处,压力的变化基本上为零.同时无论列车的头部形状和列车速度如何,涡的厚度变化很微弱.同时,涡以常速度向前传播.(4)涡的破坏阶段这个阶段的特点是,涡的厚度发生极大的变化,涡发生破坏并最终衍变成紊流状态.56第1期 骆建军等:列车进入隧道时产生涡流的数值模拟图4 列车进入隧道时形成紊流涡示意图Fig.4 Formation of the vortex when a train entering tunnel66北 京 交 通 大 学 学 报 第31卷图5 测点平面布置示意图Fig.5 Sketch of measure points inplane图6 喷射流速度历时变化曲线Fig.6Duration curve of jet velocity3 列车速度及头部形状对涡的影响为了研究列车速度以及列车头部形状对涡的影响,有必要对涡中心离开隧道端墙的距离z n ,涡中心半径r 以及涡的厚度δ进行研究.3个参数的定义见图7.图8表示列车速度对δ、r 、z n 的影响.图8(a )表示涡的厚度随着列车不断进入而发生的变化,从图8(a )中可以看出,当0<t <5s ,涡的厚度随着列车不断进入而逐渐增大,但基本上不受到列车速度的影响,此时涡处于发展阶段.在图8(b )和图8(c )中,当t <2s 时,r 、z n 的变化与列车速度无关.图7 对形成的涡各参数进行定义Fig.7 Parameter definition of the vortex图9表示列车头部形状(见图10)对δ、r 、z n 的影响.在图9(a )中,对于给定速度的列车,涡的厚度也不受到列车头部形状的影响.在图9(b )和图9(c )中,随着列车头部长细比a/b 的增大,在列车进入隧道的过程中,涡形成的时间延长,并且长细比越小的列车头部所形成的涡的厚度δ和涡中心半径r 越大. (a )涡厚度与列车速度的关系 (b )涡中心半径的变化 (c )涡中心距隧道入口的距离图8 列车速度对涡的影响Fig.8 Influence of the train speed on the vortex图11表示列车头部形状与喷射流速V jet 之间的关系.从图11(a )中可以看出,当列车速度V =20m/s 时,在3种列车头部形状下,喷射流速并不相同,而在列车速度V =45m/s 时,得到的喷射流速度相同,如图11(b )所示.显然,喷射流速度与列车速度有关系.在列车低速的情况下,这种差异是由于不同的列车头部形状所引起的.在高(较高)速以及高雷诺数下,边界层处于完全紊流状态.76第1期 骆建军等:列车进入隧道时产生涡流的数值模拟图9 列车头部形状对涡的影响Fig.9 Influence of the shape of train on thevortex 图10 列车头部形状示意图图11 喷射流速度与列车头部形状的关系Fig.10 Shapes of trainFig.11Relationship between V jet and the shape of train4 结论(1)列车进入隧道时所诱发的流动是一种低速可压缩性气体流动.在低马赫数下,当流场的流速的大小与声速比起来相差很大时,控制方程中的对流项变得“硬化”,并且时间前进方法收敛很慢;同时由于扩散项的存在,使得收敛更加困难.因此,通过引进预处理矩阵和人工压缩性方法消除了上述“硬化”现象.(2)列车进入隧道时在隧道入口附近形成涡流.在整个涡流的形成过程中,可以将这个过程分为连续的4个阶段,各个阶段具有自己的特点.(3)涡的厚度随着列车不断进入而逐渐增大,但基本上不受到列车速度的影响,而受到列车头部形状的影响比较大.(4)喷射流速度与列车速度有关系.在列车低速的情况下,这种差异是由于不同的列车头部形状所引起的.在高(较高)速以及高雷诺数下,边界层处于完全紊流状态.参考文献:[1]Auvity B ,Bellenoue M ,K ageyama T.Experimental Studyof the Unsteady Aerodynamic Field Outside a Tunnel Dur 2ing a Train Entry[J ].Exp.Fluids ,2001,30(2):221-228.[2]G lezer A.The Formation of Vortex Rings[J ].Phys Fluids1988,A31:3532-3542.[3]Gharib M ,Rambod E ,Shariff K.A Universal Time Scalefor Vortex Ring Formation [J ].Fluid Mech ,1998,360:121-140.[4]Nitsche M ,Krasny R.A Numerical Study of Vortex RingFormation at the Edge of a Circular Tube[J ].Fluid Mech ,1994,276:139-161.[5]骆建军.高速列车进入隧道产生压缩波的数值模拟及试验研究[D ].成都:西南交通大学,2003.Luo jian 2jun.The Experiment and Numerical Simulation on the Pressure Produced By a High-Speed Train Entering into a Tunnel[D ].Chengdu :S outhwest Jiaotong Universi 2ty ,2003.(in Chinese )[6]Turkel E.Preconditioned Methods for S olving the Incom 2pressible and Low Speed Compressible Equations[J ].Jour 2nal of Computational Physics ,1987,72:277-298.[7]Weiss J M ,Smith W A.Preconditioning Applied to Vari 2able and Constant Density Flows[J ].AIAA Journal ,1995,33(11):2050-2057.[8]Woods W A ,Pope C W.Secondary Aerodynamic E ffectsin Rail Tunnels During Vehicle Entry[C]∥In Proc 2th Int Symposium on the Aerodynamic and Ventilation of Vehicle Tunnels (Cambridge ,U K ,March ,1976,23-25)Cran 2field ,U K ,BHRA ,Fluid Engineering ,1976,C5:71-86.86北 京 交 通 大 学 学 报 第31卷。
浅谈大涡模拟
浅谈大涡模拟摘要:湍流流动是一种非常复杂的流动,数值模拟是研究湍流的主要手段,现有的湍流数值模拟的方法有三种:直接数值模拟、大涡模拟和雷诺平均模型。
本文主要是介绍大涡模拟,大涡模拟的思路是:直接数值模拟大尺度紊流运动,而利用亚格子模型模拟小尺度紊流运动对大尺度紊流运动的影响.大涡模拟在计算时间和计算费用方面是优于直接数值模拟的,在信息完整性方面优于雷诺平均模型.本文还介绍了对N—S方程过滤的过滤函数和一些广泛使用的亚格子模型,最后简单对一些大涡模拟的应用进行了阐述。
关键词:计算流体力学;湍流;大涡模拟;亚格子模型ﻬA simple study of Large EddySimulationDING Puxian(Central South University, School of Energy ScienceandPowerEngineering,Changsha, Hunan, 410083)Abstract:Turbulentflow isa very complex flow,andnumerical simulation is the m ain means to study it.There are three numerical simulation methods:direct numer ical simulation,largeeddy simulation,Reynolds averagedNavier—Stokes method. La rge eddysimulation (LES)is mainly introduced in this paper。
The main idea of LESis that large eddiesareresolved directlyandtheeffect of thesmalleddies on thelarge eddies is modeledby subgrid scalemodel。
高速列车脉动压力的大涡模拟及小波分解
高速列车脉动压力的大涡模拟及小波分解
王亚南;陈春俊;何洪阳
【期刊名称】《机械设计与制造》
【年(卷),期】2015(000)008
【摘要】高速列车表面脉动压力是引起气动噪声的主要根源,研究车体表面脉动压力对噪声控制等方面有重要意义.采用大涡模拟(LES)仿真计算高速列车运行时头车和尾车外流场的脉动压力,利用二进正交db小波将脉动压力分解为能量互不重叠的正交频带,并分析脉动压力在各频带上的能量分布规律.数值仿真结果表明:列车表面脉动压力由平均压力和在平均压力附近上下波动的脉动部分组成,脉动压力在全频带均有分布,且主要集中在低频区域;随着列车运行速度的提高,车体表面脉动压力幅值迅速增大,主要能量向高频区域移动;头车、尾车脉动压力变化趋势相似,且头车脉动压力大于尾车脉动压力.
【总页数】4页(P86-88,93)
【作者】王亚南;陈春俊;何洪阳
【作者单位】西南交通大学机械工程学院,四川成都610031;西南交通大学机械工程学院,四川成都610031;西南交通大学机械工程学院,四川成都610031
【正文语种】中文
【中图分类】TH16
【相关文献】
1.高速列车车头曲面脉动压力的大涡模拟 [J], 肖友刚;张晓缝;康志成
2.基于小波变换的高速列车表面脉动压力提取 [J], 陈春俊;何洪阳;邵云龙
3.边界层人工转捩的大涡模拟及壁面脉动压力研究 [J], 徐嘉启;梅志远;刘志华
4.应用EEMD和小波包分解的压力脉动信号时域特征提取方法 [J], 李瑞;谷立臣;赵鹏军;郭西惠
5.跌坎突扩型消力池脉动压力大涡模拟研究 [J], 卢洋亮;尹进步;张曙光;杨钊因版权原因,仅展示原文概要,查看原文内容请购买。
高速列车舱内气流分布的数值模拟
要:高速列车的高密闭性要求空调座舱满足驾驶员和旅客的热舒适性要求, 而传统的舱内通风设计依赖自
由射流的经验公式, 确定座舱内温度场和速度场的方法, 因无法考虑送排风气流和室内障碍物的影响, 通风设计 很大程度上依赖于模型实验。本文对旅客车厢采用传统的有限差分法, 将送排风气流与车厢形状及座椅作为一 体来考虑, 通过计算机来模拟车厢内温度场和速度场, 以及确定出送风速度和送风温差; 对驾驶室, 根据高速列车 机车形状, 采用贴体网格生成法和计算机仿真技术模拟了三维驾驶室形状及室内气流分布。 关键词:高速列车;通风设计;气流分布;数值计算 中图分类号: 文献标识码:
, ,
( , , , ; ! " # $ " % " % $ $ & ’ ) * & * $ & % & * $ & ’ , ( + ( , , , ) * * * " * $ $ & , & % & $ $ &, & $ % $ $ &’ , + + . (
[ ] 张量形式的通用输运方程 , 然后进行如下的坐标变
式中, , , , , , 和 和 对 , 和 的一阶导数。
, , 分别表示 ,
三维计算网格的生成
建立了合适的控制方程后需将微分方程离散到差 分网格上求解。差分网格的设计, 在很大程度上确定 了数值计算的精度及敛散性, 尤其是对复杂的边界形 状引起的复杂的流动现象。 笛卡尔直角坐标下的差分网格为矩形网格, 网格 划分的粗密程度视计算精度和计算成本而定。图 为
低雷诺数翼型层流分离现象大涡模拟方法
低雷诺数翼型层流分离现象大涡模拟方法朱志斌; 刘强; 白鹏【期刊名称】《《空气动力学学报》》【年(卷),期】2019(037)006【总页数】9页(P915-923)【关键词】低雷诺数; 翼型; 层流分离; 转捩; 再附; 大涡模拟【作者】朱志斌; 刘强; 白鹏【作者单位】中国航天空气动力技术研究院北京 100074【正文语种】中文【中图分类】V211.30 引言随着近年来临近空间低速飞行器和高性能微小型飞行器的快速发展,翼型在低雷诺数范围内的气动问题日益受到关注。
低雷诺数条件下(1×104~1×105)翼型绕流流动会出现层流分离现象,其形成和演化会严重损害翼型气动性能,并直接影响飞行器续航能力以及操纵性和稳定性[1]。
准确预测翼型层流分离流动现象,并分析认识其气动特性及流动机理,对于低速临近空间飞行器和微型飞行器的设计研制具有重要意义。
对翼型层流分离流动的研究可追溯至20世纪60年代,Gaster和Horton提出了经典层流分离泡模型[2-3],认为层流附面层在逆压梯度作用下分离,在空间发生转捩后,外部高能量气流的引入导致再附,形成层流分离泡。
大量的风洞试验[4-5]和数值模拟[6-7]研究发现,翼型层流分离流动存在显著的非定常特征和非线性效应。
层流分离流场中存在流动分离、转捩、再附等非定常流动结构,并不存在稳定的分离泡,而是一系列的大尺度旋涡结构时均结果。
基于雷诺平均假设的RANS(Reynolds Averaged Navier-Stokes)方法难以准确预测层流分离流动中的转捩现象。
为此,Yuan[8]、Radespiel[9], Lian[10]分别发展了基于稳定性分析的eN方法与RANS联合求解的方法,用于预测工程问题中的低雷诺数转捩过程;此外,Choudhry[11]等采用γ-Reθ和k-kL-ω转捩模型模拟了NACA0021翼型低雷诺数流动,对翼型层流分离泡结构、特性及随迎角、雷诺数、湍流度的变化规律展开了细致的研究。
铁道机车车辆通用旋转式涡流制动_RWB_
4/200646电力牵引足,近年来试制的RWB制动系统已在日常中使用。
2.1RWB与LWB的对比ICE3采用LWB制动就是利用了涡流制动原理的优点。
对于自300 km/h速度起的制动,尤其是在超过40‰下坡区段上,为使列车达到期望的制动功率(控制列车下坡运行时的较大动能)采用LWB制动是绝对必要的。
ICE3的驱动功率(也包括发电制动功率)每200 m单位达到8 MW。
ICE3的LWB的最大制动功率为13MW。
这两种制动的电制动功率为21 MW。
如果在同一地段连续使用LWB制动,那么在LWB制动时就会使钢轨发热,钢轨焊接处就会产生附加应力。
这种附加应力要求有相应抗扭曲的线路(重型碎石道床或整体道床的线路),且不得超过规定的温度值(在整体道床轨道上约为80℃)。
此外,线路中的电气开关装置必须耐受涡流电流。
与LWB相比,RWB有如下一系列优点:(1)任意类型的铁道机车车辆都可补装RWB系统;(2)可部分替代摩擦制动;(3)在某一速度范围内可用RWB作为常用制动;(4)根据至今所掌握的技术,在所有铁路线上可不受限制使用RWB。
2.2RWB与摩擦制动的比较使用常规的闸瓦制动或盘形制动,会出现以下现象:(1) 车轮有很大的热负荷,特别是在长大下坡道上制动时(每一车轮的平均热功率达50 kW);(2) 声音扩散值很高,尤其是货物列车时,在速度v =100 km/h时声音扩散值达91 dB(A)。
这里没有考虑车辆装备吸音器。
例如改变摩擦副可减小盘形制动噪声和磨耗。
此外,高速列车制动的摩擦副使用的是钢制动盘和烧结合金制动盘,这些制动盘的价格比常用材料的要贵1倍,烧结合金闸片的价格比常用材料的要贵9倍。
铁道运输中使用RWB的目标是:在保持或改进性能的前提下,减少或克服机械闸瓦制动或盘形制动系统固有的缺点,从而进一步获得安全性、环境友好性和经济性的效果。
RWB可以取代摩擦制动或补充摩擦制动。
在载重汽车上RWB与摩擦制动的比例为1∶4 ̄1∶10。
高速列车车头曲面脉动压力的大涡模拟(1)
摘
要:为了得出车头曲面脉动压力的分布规律, 利用参数变换技术, 生成了高速列车车头曲面的四边形贴体网
格, 进而生成了高速列车头部流场的六面体网格 。对车头流线型曲面的脉动压力进行了大涡模拟, 得出了脉动压 力的分布特性:在流场中某点处产生的压力主要影响该处本地, 而对其余点的影响作用则随着距离的增加而迅速 无明显的主频率。各点脉动压力的频谱在低频时幅值较大, 随着频率升高, 幅值以负 减弱;脉动压力的频带很宽, 随着列车运行速度提 指数规律持续下降。脉动压力 1 /3 倍频程频谱的主要能量集中在 20 ~ 500 Hz 频率范围内, 高, 频谱的主要能量范围有向高频移动的趋势 。 关键词: 高速列车;车头曲面;脉动压力;大涡模拟 中图分类号:U238 文献标识码:A
228
四川大学学报( 工程科学版)
第 42 卷
[7 ] 法估算了受电弓脉动压力的分布 ; 者基 [8 ] 于参数变换技术 , 生成了高速列车车头曲面的四 边形贴体网格, 并建立了计算车头曲面脉动压力的
( 1 ,v1 , ( u i, 的边 界 点 坐 标 是 ( 0 ,v0, j ), j ), 0 ,0 ) , 1) , ( u i, 利用这些边界点进行线性插值, 可得相对 1, 弧长参数空间的表面节点坐标 ( u i, j ,v i, j ) 的计算公 式: u i, j = ( 1 - v i, j ) u i, 0 + v i, j u i, 1 v i, j = ( 1 - u i, j ) v0 , j + u i, j v1 , j (2) 对上述方程组进行迭代求解, 就可以得到相对弧长 v i, 参数空间的表面节点坐标( u i, j, j)。 1. 3 参数网格向物理空间的映射 得到相对弧长参数空间内的参数网格 u i, j ,v i, j 后, 还需把该网格向物理空间映射, 以得到物理空间 y i, 内的表面网格( x i, j, j ,z i, j ) 。 该过程如下 : 相对弧 0) , (0, 1) , (1, 长参数空间中四个角点的坐标是( 0 , 0) , (1, 1) , v i, 对点 ( u i, 对应的 4 个边界点坐标 j, j), 0) , ( u i, 1) , ( 0 ,v0,, ( 1 ,v1 , 是 ( u i, 将这 0, 1, j), , j), 些角点和边界点在物理空间中对应的坐标向量记为 R0 , R0 , R1 , R1 , Ri , Ri , R1,j , , 可以 0, 1, 0, 1, 0, 1 ,R0 , j, 得到 ( u i, j ,v i, j) 点映射到物理空间中的坐标向量 R i ,j 。 R i, j = ( 1 - u i, j ) R0, j + u i, j R1, j + ( 1 - v i, j ) R i, 0 + v i, j R i, 1 - ( 1 - u i, j ) ( 1 - v i, j ) R0, 0 - u i, j ( 1 - v i, j ) R1, 0 - ( 1 - u i, j ) v i, j R0, 1 - u i, j v i, j R1, 1 (3) 采用上述方法, 笔者开发了车头曲面流场计算的网 格生成软件, 不但解决了自由型曲面网格的生成问 题, 而且实现了车头曲面网格生成的程序化。 使用 可自动控制网格疏 者根据程序提示输入相关参数, 密, 并自动生成车头曲面的四边形贴体网格 , 进一步 可自动生成计算车头曲面外部流场的六面体贴体网 格。图 1 是车头曲面的初始局部网格, 图 2 是插值 并映射后车头曲面的局部网格, 图 3 是该软件输出 的车头曲面的四边形贴体网格。图 4 是该软件输出 的车头曲面流域的六面体网格。
基于大涡模拟方法的三维旋转爆轰流场结构研究
基于大涡模拟方法的三维旋转爆轰流场结构研究
雷特;武郁文;徐高;邱彦铭;康朝辉;翁春生
【期刊名称】《兵工学报》
【年(卷),期】2024(45)1
【摘要】为研究环形燃烧室中边界层、黏性、湍流模拟方法对旋转爆轰流场结构的影响,采用开源计算流体动力学软件OpenFOAM,以氢气为燃料、空气为氧化剂,基于大涡模拟(Large Eddy Simulation,LES)方法、RANS方法、Euler方法,分别结合滑移和无滑移边界,对三维旋转爆轰发动机模型进行数值模拟,分析对比不同计算方法下旋转爆轰流场结构。
着重讨论以LES方法得到的流场结构。
研究结果表明:当采用滑移边界时内、中、外截面的流场温度无太大差异,当采用无滑移边界时内外壁面温度高于中间截面,边界层会影响近壁区域气体的流动速度,导致内外壁面爆轰波高度低于中间截面,还会影响燃烧产物的流动状态轴向截面上爆轰波波头产生变形;不同湍流计算方法得到的旋转爆轰流场结构存在相似性,黏性是影响旋转爆轰流场结构的主要原因。
研究结果对于揭示边界层和黏性对旋转爆轰的影响机制具有一定的科学意义。
【总页数】12页(P85-96)
【作者】雷特;武郁文;徐高;邱彦铭;康朝辉;翁春生
【作者单位】南京理工大学瞬态物理国家重点实验室
【正文语种】中文
【中图分类】V231.22
【相关文献】
1.气体与液体两相连续旋转爆轰发动机爆轰波传播特性三维数值模拟研究
2.连续旋转爆轰三维流场的数值模拟
3.连续旋转爆轰发动机隔离段流场数值模拟研究
4.基于OpenFOAM的三维H2/Air连续旋转爆轰流场数值模拟
5.基于煤油的旋转爆轰发动机流场数值模拟
因版权原因,仅展示原文概要,查看原文内容请购买。
基于大涡模拟的高速列车横风运行安全性研究
基于大涡模拟的高速列车横风运行安全性研究第33卷第4期2011年4月铁道学报JOU R NA L OF T H E CH IN A RA ILW A Y SO CI ET Y V ol 133 N o 14A pr il 2011收稿日期:2009-10-27;修回日期:2010-08-20基金项目:国家重点基础研究发展计划(2007CB714701);/十一五0国家科技支撑计划(2009BAG12A01-C09)作者简介:刘加利(1985)),男,四川成都人,博士研究生。
E -mail:liujl_123@1631com文章编号:1001-8360(2011)04-0013-09基于大涡模拟的高速列车横风运行安全性研究刘加利, 于梦阁, 张继业, 张卫华(西南交通大学牵引动力国家重点实验室,四川成都 610031)摘要:结合高速列车空气动力学和多体系统动力学,研究横风对高速列车运行安全性的影响。
首先采用大涡模拟计算方法,研究了不同横风风速下高速列车非定常气动载荷的时域及频域特性,列车周围流场结构及相应的非定常流场特性。
然后建立高速列车多体系统动力学模型,将得到的气动力作为外加载荷作用于列车上,研究了不同横风风速下定常气动力和非定常气动力对直线上高速列车运行安全性的影响特性,计算结果表明,与定常气动力相比,作用于车身上的非定常气动力使列车的振动加剧。
最后参照高速列车的安全运行标准,对高速列车的安全运行进行分析,为横风下高速列车的安全运行提供参考。
关键词:高速列车;横风;大涡模拟;非定常气动力;运行安全性中图分类号:U 270111 文献标志码:A doi:1013969/j 1issn 11001-8360120111041003Study on Running Safety of High -speed Train under Crosswindby Large Eddy SimulationLIU Jia -li, YU M eng -g e, ZH ANG J-i ye, ZH ANG We-i hua(State Key Laboratory of Traction Pow er,S ou th w est Jiaotong University,Chengdu 610031,China)Abstract:The aerody namics and mult-i body system dy namics of the hig h -speed train are combined in this paper to study the effects o f crossw inds on the running safety of the high -speed tr ain 1Fir stly,the larg e eddy simula -tion (LES)is used to study the time and frequency dom ain characteristics o f the unsteady aero dynamic forces on the hig h -speed train,the structure of the flow field around the train and the char acteristics of the unsteady flow field under different crossw ind speeds 1Then the mult-i body system dynam ics m odel of the hig h -speed train is established,in w hich aerody nam ic forces are dealt w ith as the external loads acting on the train 1With this model,the influence of steady and unsteady aerodynam ic fo rces o n the running safety of the high -speed train un -der cr ossw inds on straight lines is studied 1The r esults show that the vibration of the tr ain under unsteady aero -dynamic forces is more intense than that under steady aer ody namic for ces 1In the end,w ith reference to the safe operation criteria for high -speed trains,the running safety o f the high -speed train is studied,w hich is helpful to safe operation of high -speed trains under cro ssw inds 1Key words:high -speed train;cr oss w ind;large eddy sim ulation(LES);unsteady aerody namic force;running safety随着列车运行速度的不断提高,许多在低速时被合理忽略的问题都逐渐凸显,并且在很大程度上影响着列车的提速[1]。
高速空调列车内气流组织的大涡模拟
高速空调列车内气流组织的大涡模拟匡骁;齐朝晖;李超;韦婷婷【期刊名称】《铁道学报》【年(卷),期】2009(031)003【摘要】高速空调列车车厢内气流的温度场和速度场研究是空调列车内气流组织设计的重要基础,也是空调车内舒适环境评价的依据.考虑到高速列车车厢内复杂的几何结构、数值模拟的边界条件以及空气流场的湍流特性,本文采用湍流大涡模拟方法对高速列车车厢内空气流动与热质传递过程建立了数学模型,采用有限容积法进行区域离散,应用均匀六面体网格划分车厢,并在同位网格的基础上采用SIMPLEX 算法,考虑车内座椅、行李架等障碍物以及车体各壁面辐射、车窗热流、乘客散热和高速列车运行特性等因素的影响,对高速列车车厢内空气流场和温度场进行数值模拟.其结果对高速列车的空调效果及车内舒适环境的优化提供了依据.【总页数】6页(P94-99)【作者】匡骁;齐朝晖;李超;韦婷婷【作者单位】中南大学,制冷空调研究所,湖南,长沙,410075;青岛四方车辆研究所有限公司,山东,青岛,266031;中南大学,制冷空调研究所,湖南,长沙,410075;中南大学,制冷空调研究所,湖南,长沙,410075;中南大学,土木建筑工程学院,湖南,长沙,410075【正文语种】中文【中图分类】U271;O357【相关文献】1.关于空调列车车厢内空气品质的现状分析与对策 [J], 胡春莲;张登春2.空调列车中的气流组织与人体热舒适性 [J], 陈焕新;刘瑞;杨培志3.空调列车车厢内确定新风量的指标分析 [J], 杨培志;陈焕新4.基于乘客热舒适性的空调列车软卧包厢送风气流组织优化研究 [J], 刘志永;尹海国;孙翼翔;李安桂5.空调列车室内气流组织的三维数值模拟 [J], 史自强;卢纪富;靳谊勇;胡世兴因版权原因,仅展示原文概要,查看原文内容请购买。
运行速度对列车火灾温度分布影响的大涡模拟与分析
运行速度对列车火灾温度分布影响的大涡模拟与分析张丽荣;马纪军;夏冬宁【期刊名称】《铁道车辆》【年(卷),期】2012(050)002【摘要】以某列车车厢实际尺寸为依据,采用美国NIST开发的大涡模拟软件FDS 建立列车火灾模型,以列车运行速度为变量,进行列车火灾数值模拟与分析,重点分析了典型低速和典型高速时的温度场分布情况.结果表明,列车运行速度在120 km/h 以下时,车厢内的平均温度在240℃以下,最高温在车厢外表面;运行速度在120 km/h以上时,列车速度越高,车厢内的平均温度越高.%With the actual dimensions of a certain train compartment as the basis, with the use of the large eddy simulation software FDS developed by the American NIST, the model of train on fire is set up, with the train operation speed as the variable, the digital simulation and analysis of train on fire are carried out, the temperature field distributions at typical low speed and typical high speed are mainly analyzed. The result shows that, when the train running speed is below 120 km/h and the average temperature inside compartment is below 240 , the maximum temperature is at the outside surface of the compart-ment; when the running speed is above 120 km/h, the higher the train speed is, the higher the average temperature inside compartment will be.【总页数】4页(P8-11)【作者】张丽荣;马纪军;夏冬宁【作者单位】唐山轨道客车有限责任公司,河北唐山063035;唐山轨道客车有限责任公司,河北唐山063035;唐山轨道客车有限责任公司,河北唐山063035【正文语种】中文【中图分类】U271【相关文献】1.隧道列车火灾的顶棚射流平均温度分布2.基于大涡模拟的地铁区间火灾横向温度分布研究3.高速铁路列车运行速度提高对车体加速度的影响分析4.地铁列车运行速度对客室内火灾烟气扩散的影响5.运行速度对高速列车火灾热释放速率的影响因版权原因,仅展示原文概要,查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
旋转列车气流的大涡模拟Hassan Heniida.Naliia Gil,Chris Baker摘要利用大涡模拟(LES)方法研究高速列车的气流问题,釆用标准的Smagorinsky模型模拟亚格子应力。
列车模型是由4辆车组成的1/25比例的ICE 2型列车。
该模型被放置于直径为3.61m的旋转试验台上.基于列车的高度和速度,分别对雷诺数77 000和94 000 的大涡进行了模拟。
模拟中运用了粗糙的、中等的和加密的3种讣算网格。
这三种网格分别由6X106, 10X 106 ,和15X106个节点组成。
加密网格的计算结果与试验数据吻合较好。
运用大涡模拟获得了不同的流动区域:上流区、鼻端区、边界层区、风挡区、尾流区和远尾流区。
在靠近列车鼻端区域从气流的最大速度幅值中可以出现局部的速度峰值。
面压力的最大值和最小值分别出现在黑近鼻尖区域的顶面和底面。
所有的湍流结构都产生于列车的顶部,并被列车外侧的径向速度分量所掠过。
在列车的外侧,主要是大结构的大湍流占据主要地位。
研究表明,以柱面形式支撑的风挡和车下复杂结构对气流的速度有很大的影响。
在合适的雷诺数范圉内,气流流速与列车速度近似地呈线性关系。
1.引言列车在空气中运行时,会导致列车两侧以及尾部的气流产生重要的气流流速。
这种现象会对乘客和铁路沿线工作人员的安全造成威胁,同时也会给婴儿车以及手推车带来很多问题。
鉴于对外部环境所造成的影响,铁路安全与标准学会(RSSB)⑴近期已将其确定为亟待解决的课题,各种研究工作也需要开展。
RSSB最新的一项研究显示,在英国,与其他危险相比,所有与列车气流相关的危险所占比例较小。
然而,如得不到有效的组织管理,列车气流会对站台乘客以及铁路沿线工作人员的安全造成很大的威胁。
自1972年以,英国大陆地区已经报道了24起事件,这些事件不但涉及到气流产生的作用力对静止站台上的婴儿车、手推车所造成的伤害,而且也有对乘客及其物品的伤害[29。
作用于站台或沿线物体的气动力取决于物体的外形,并与气流流速的平方成正比。
气流流速本身通常与列车的速度成一定的比例。
自新一代列车通过提高速度减少运行时间来适应铁路运营客流高峰以来,气流造成的危险性变得更大,急需探究高速列车气流的产生机理。
当前有两种途径来研究列车的气流问题,全尺寸模型或缩尺模型。
所有这些方法在文献[4]〜文献[5]中都有所阐述。
文献[5]运用热膜的方法来测试列车周圉的速度,该热膜覆盖于放置在移动设备上的山4辆成的1/25比例模型的表面,在他们的工作中也进行了一些粒子图像速度测量方法(PIV)的试验。
文献[4]开展了全尺寸货车和客车的模拟工作。
然而,山于列车的气流是大湍流的瞬态现象,为获得足够的结果以得出时间的平均和标准偏差,需要进行大量的试验。
这样就会使这两种试验方法耗费大量的时间和财力。
基于此,伯明翰大学的铁路研究中心运用直径为3.16m的旋转台来测量气流的流速。
这种新的试验方法使得测量沿列车长度的某些点的平均速度变为可能。
然而,这种方法仍旧不能得到流动域和过渡部分气流以及压力场的全貌,超出了试验所能探讨的范圉⑹。
过去,单纯依靠数值仿真预测列车周圉随时间变化的气流流动,即使雷诺数相对比较低的流动也都是不可行的。
最近儿年,讣算机配置的提升使得这些模拟(适宜的雷诺数)变成了可能。
对比高昂的、无法实现或难以实现的试验研究,计算流体动力学(CFD)的应用更具优势。
通常,试验需要大量的时间和财力。
曲于受到雷诺数的制约,通过缩尺的风洞试验来弥补这些不足是不可能的。
同样对于非同寻常雷诺数的流动算例,在风洞中去除边界和/或条件在操作上都是无法实现的。
况且,空旷空气的试验主要受到大气条件的制约,试验重现的模拟显得极为困难。
伴随着列车的提速以及新一代高速列车产品的问世,很多在低速完全可以被忽略的空气动力学问题,如气动噪声、列车交会产生的作用力、隧道出口压力波以及车内乘坐不舒适性等,都将显现。
这些是限制列车系统提速的主要因素⑺。
因此,在过去20年中, 高速列车空气动力学得到了很多研究者的重视,并一直试图改善其气动性能Sa]。
文献[13]中,利用CFD模拟列车气动模型(ATM)周圉的流动,采用分涡模拟(DES)来分析列车尾流的流动状况,并将CFD的讣算结果与试验数据进行了对比,运用分涡模拟方法获得了儿乎所有工况下的讣算结果。
同时,尾部和边界层的速度量峰值采用DES方法进行了评价。
最近,采用标准的SmagoHnsky模型的大涡模拟方法(LES), 已经被成功地应用于硏究列车周围侧风的流动[1446]o本文针对旋转模型周围气流的瞬态行为,釆用标准的Smagorinsky模型的LES方法进行研究。
研究中所用到的77 000和94 000这2个雷诺(Re),是基于列车的高度和速度而得到的。
利用计算软件CFX”】的滑移网格技术来模拟模型的旋转。
2. 试验方法旋转轨道台由一个直径为3.16m的铁路轨道框架组成,最高转速可以达到118r/niin,相当于轨道运行速度为22m/s (图1)。
起初设想将标准轨道段置于框架上,并借助于新的激光技术〔习使研究丄作成为可移动的有机合成体,用以研究售票员鞋子静止结冰问题,售票员鞋子保持在旋转轨道上温度低于零摄氏度以下的条件。
LI前的研究中,高度H=70mm、长度Z=0.5m的1/25的ICE 2型列车的模型,通过工装和螺钉(每节车2套)安放在旋转轨道上,保证两者的运动是一体的(见图1 (b))。
该模型由3个风挡连接组成。
典型风挡的深度和宽度 炉相等,同为列车高度的%。
为了减少轮子 旋转所引起的湍流,同时也为了模拟列车运动的地面层,比轨道稍大的带有圆形狭槽的 木质平台被放置于装配台上。
列车和地平面之间是紧密的刷子。
根据列车置于轨道上的 方法可以使车辆很方便地编挂或者解编。
这样,任何编组形式的列车都可以进行流场的 测试。
该试验台最为重要的方面当然是多种列车通过方案可以很快地得以实现(每次只 需旋转一次试验台),同理,与现有的全尺和缩尺技术相比,该试验方法更利于快速操 作。
(b )支撑轨道和列车的装置示意图3. 列车模型列车模型山4辆弯曲的车辆组成。
平直的车辆所带来的问题是列车模型长约50cm, 相对于曲率半径R=180.5cm 的试验台来说尺寸过大,这样就会导致在列车方向每节车厢 的后部都存在明显的不连续性。
并且,探针与列车表面之间的距离也不是恒定量,因此, 探针的测量位置在列车的顶部和侧面会对测试结果产生非常明显的影响。
为了解决这些 问题,构造一列弯曲的列车,使其在靠近车辆的弯曲线上进行测量成为可能(见图1(a ))。
此外,这种新型的列车被设计成简单的1/25比例的ICE 2列车,以便于所得到的测试结 果可以与TRAIR 试验台获得的结果进行对比。
列车模型111 1/25比例的4辆ICE 2列车 组成。
很显然,这些结果将会受到列车曲率的影响(试验车辆越多影响越大),凸面侧 要比凹面侧边界层增加得多(见文献[18]〜文献[19])。
在该试验中,气流和尾流速度釆 用静止的眼镜蛇探针(湍流所使用的P/L 仪器)进行测试,四孔压力传感器用以测试速 度的3个分量,能够测试的最高速度为lOOnVso 探针的详细信息见文献[6]。
4. 数值方法运动列车周用的流动是非常复杂的,山大范围的湍流所组成。
据证实,大湍流是影 响列车空气动力学的主要因素。
然而,尾流和边界层受大湍流结构所主导。
因此,一种(a) (b)(a )列车模型和测试探针兰州交通大学毕业设计(论文)(5) 解决大结构的模拟方法例如LES 方法,能够更好地模拟列车周圉的流动。
尽管LES 方 法比平均雷诺方程Navier-Stokes (RANS)法计算更为昂贵,但是它可以提供更为准确 的时间均值结果,同时也可以得到瞬态流动的信息,这些都是RANS 方法无法媲美的。
LES 法已经被证实,在预测简化列车和非流线体周圉气流流动方面是一种可黑的技术 [14-16]通常,LES 法将流动的结构分解成大的和小的等级。
将小等级对大等级流动的影响 建成模型后,大活动级的流动就可以直接被模拟。
因此,LES 法是一种介于RANS 法和 直接数值模拟法之间的一种方法。
对于不可压缩流的简化的连续性方程和动量方程如下:式(2)中:乞讣和诸分别是滤波后的速度、滤波后的压力和残余应力。
残余应力为未知量,需要模拟。
线性涡粘性模型中亚格子应力与应变率张量爲的关式(3)中:比例系数卩的是亚格子运动的涡粘性系数。
V SGS =仏)*2$泸》 式(4)中:A 是当量宽度;Cs •是Smagoriiisky 系数;/是近墙的阻尼函数。
SGS 模型系数值•此处取值为O.lo III 于壁面的应变率张量不为零,因此,模型系数值应当适当考虑壁面的影响。
近墙的阻尼函数/•用于控制靠近壁面的湍流长度。
该函数可表 达为:式(5)中:厂是标准的壁面距离。
ANSYS-CFX 中关于LES 的更多内容参见文献[17]。
系式如下: (4)5. 计算区域和边界条件数值条件的设置尽量与试验条件相一致。
计•算域由2部分组成:静止部分和旋转部 分(图2 (a) ) o 静止部分正方形外形的边长为6m,高度为6H°静止区域的地面采用 无滑移的边界条件,侧面和顶部采用开放的边界条件。
旋转部分距离列车两侧面的距 离为2R,高度为4H (见图2 (b) ) o 列车模型、轨道以及支撑柱面都位于可旋转的部 分。
旋转部分顺时针方向旋转的角速度等同于列车的运行速度。
旋转部分的底面也采用 无滑移的边界条件。
为了同试验具有相似的边界条件,轨道与列车旋转的速度其大小和方向均一致,平 台旋转方向与旋转框架的旋转方向相反。
两部分间的瞬态网格交界面在每个时间步长内 都进行信息交换。
模型壁面采用无滑移的边界条件。
静止和旋转部分模拟的初始速度均 为零。
图2 (b)显示的是笛卡尔坐标系。
区域的圆心位于平台的水平面上。
x-y 坐标平 面对应于水平面(平台平面),z 轴方向垂直于x-y 平面向上。
计算结果本节介绍基于列车高度和速度获得的77 000和94 000两个雷诺数下的LES 模拟计 算结果。
气流的速度为列车的运行速度,列车的高度为边界层的厚度。
大部分是对低雷 诺数Re 时的气流分析,在本节的最后将会介绍雷诺数Re 对气流结构的影响。
6.1数值模拟精度估算LES 方程采用三维有限体积法进行求解。
为了研究网格质量对汁算结果和数值精度 的影响,进行了粗糙的、中等的和加密的三种网格质量工况下的计算。
运用商业网格划 分软件ICEM-CFD 进行实体建模和网格划分。
利用ICEM-CFD 在模型周围生成六面体 网格。
列车周围厚度0.1H 处釆用O 型网格,最终在各个方向生成光滑网格(图3 (b))。