反旋进气盘腔内流动与换热的数值模拟

应用RNGk-ε湍流模型对围屏上带反旋喷嘴的径向内流旋转涡轮盘腔内的流动与换热进行了数值模拟,揭示了盘腔内的压力损失及冷气流量、旋转雷诺数、湍流参数等因素对盘腔内流动换热的影响.计算结果表明:盘腔内的流动结构主要由湍流参数决定;在某一旋转雷诺数下盘腔内压力损失随冷气流量的增大而呈现S型变化;反旋喷嘴的应用能有效地降低盘腔内的压力损失;转盘附近的努赛尔数随冷气流量及旋转雷诺数的增大而增大.      

超临界压力下正庚烷的湍流传热数值研究

 应用修正的对应态方法、基本的热力学关系式和Soave-Redlich-Kwong(SRK)状态方程计算了超临界压力下流体的传输物性和热力学性质,并将这些普适性的物性计算方法与计算流体力学（CFD）数值模型相结合,针对火箭推进系统、高超声速飞行器相关的超临界传热和发动机冷却现象,在准确考虑热力学性质和传输物性变化情况下,系统地研究了一种典型的碳氢燃料——正庚烷在超临界压力下的湍流传热现象,详细分析了超临界压力对流动和传热的影响,揭示了对流换热努塞尔数的变化规律,并将数值模拟计算结果与已有经验公式进行了比较。计算结果显示在超临界压力下,努塞尔数随着压力的降低而降低,并且在临界区域附近会出现湍流传热的恶化和波动现象。     

平行平板内连续布置肋块流场与传热的数值模拟

为了研究平行平板内连续布置8个肋块的流场结构和传热现象, 利用计算流体力学软件FLU-ENT6.2内置的标准k-ε湍流模型、RNGk-ε湍流模型和Realizablek-ε湍流模型对此进行了数值模拟并对结果作了比较.流场结果显示在相同来流条件下, 由Realizablek-ε模型预测的肋块间逆时针方向回流区最大, 标准k-ε模型最小.数值结果表明从第2个肋块开始加热段表面压力系数Cp和努塞尔数Nu开始呈周期性变化.同时将第7个肋表面的努塞尔数数值解与实验值进行了比较, 结果为三种湍流模型都能准确模拟出努塞尔数的变化趋势.计算结果可为工程设计和使用提供参考.      

航空刹车机轮自然对流的数值模拟

采用商业软件Fluent对简化的航空机轮内部腔体的自然对流进行数值模拟,得到了中间环形狭缝瑞利数在106～1011范围下腔体内部的流动换热规律。对腔体内部的复杂的流动和换热进行分析,获得了机轮内部腔体关键部位的流动特征,并拟合出了刹车盘面与轮毂间对流换热的关联式。当环缝内瑞利数达到108,流动进入湍流状态。同时将数值模拟的结果与有关文献做了对比,证明了本文计算方法的正确性。     

稀薄流高超声速钝化前缘逆向射流流场特征研究

为准确预测稀薄过渡区逆向射流对稀薄大气的干扰流动特征，本文采用直接蒙特卡洛算法（DSMC），在自由流马赫数为7的情况下，对稀薄流中逆向射流干扰下的高超声速流动进行了数值模拟。计算中考虑了60～90 km高度的典型大气环境，研究了压比和自由流克努森数对压力、剪切力和热流密度的影响规律。结果表明，逆向射流改变了流场结构，且具有明显的热防护性能，但减阻效果不明显。逆向射流对高超声速稀薄来流的影响随压比的增加明显增强，形成的马赫盘大小也随之增加。高度的增加伴随稀薄效应不断增强，自由流克努森数增大，导致弓形激波厚度随高度增加而增加，壁面处的回流区范围不断减小至消失。这项研究初步揭示了稀薄流中逆向射流/高超声速流动相互作用的机制。     

低位预旋进气转静系盘腔流动特性研究

采用实验和数值模拟的方法对低位预旋进气转静系盘腔的流动特性进行研究,考察了流量与动盘转速变化对腔室内流动的影响,结果表明腔内径向静压差随流量增大而增大,随转速升高而增大;引入腔内有效角动量的概念,用角动量系数表征其与盘腔入口角动量之比,通过分析来流角动量与盘面转矩的相对大小,推导得到影响角动量系数的两个主要无量纲参数为修正的湍流流动参数与预旋来流雷诺数;对于给定结构形式的预旋系统,当改变转速时,角动量系数随修正湍流流动参数增大而减小,当修正的湍流流动参数大于0.0646时,角动量系数变化缓慢,腔内有效角动量约为入口角动量的70%～90%;当改变来流速度时,角动量系数随修正湍流流动参数增大呈现先减小后增大的变化趋势。     

转子叶片径向受限的“冲击-气膜出流”冷却结构流动与换热

为揭示转子叶片径向受限的"冲击-气膜出流"冷却结构流动换热规律,以某型双层壁叶片肋化分割形成的冷却单元为研究对象,通过数值模拟的方式,对冲击雷诺数Rej,旋转数Ro,无因次温比(Tw-Tf)/Tw等参数变化下流场和换热特性变化规律展开研究。结果表明:在哥氏力和离心力作用下,受限空间内存在射流偏转、径向二次流动以及二次冲击等现象;流动的径向受限可抑制射流偏转,强化冲击换热;相同的旋转数Ro下,逆转向冲击(叶背区)换热努赛尔数Nu比顺转向冲击(叶盆区)高8%。在研究的参数范围内,数值模拟和试验结果说明径向受限周向出流结构能有效的抑制旋转对换热的削弱。     

旋转弯管内耦合对流换热特性研究

由于离心力，科氏力，和离心型浮力（Centrifugal-type Buoyancy force）同时出现，使得旋转弯管内的流动和换热特性更为复杂。本文考虑了温度与密度的耦合关系，对旋转弯管的对流换热问题进行了数值研究。分析了不同F数（科氏力和离心力之比），Ran数（旋转Rayleigh数）,以及Ro数（Rossby数）下的二次流动、轴向主流、温度场、摩擦系数比、努塞特数比变化规律。为工程上解决旋转弯管道内耦合对流换热问题提供了可靠的理论依据。     

涡轮叶片内部沿周边分布的竹节孔换热特性研究

根据发动机涡轮叶片内部沿周边分布的竹节孔冷却通道结构建立了简化的传热分析模型。首先通过无量纲分析得到了竹节孔的平均努塞尔数不仅与雷诺数和普朗特数有关,还与“竹节”形状、通道外壁形状以及流体和固体的导热系数比等参数有关。然后,对雷诺数和冷却通道几何参数对竹节孔平均努塞尔数的影响进行了数值模拟研究,计算并对比了上述参数对平均努塞尔数的影响规律,结果表明：（1）“竹节”结构使平均努塞尔数显著增大,而通道外壁形状主要影响了局部热流密度和局部努塞尔数的分布情况;（2）随着通道内的雷诺数和“竹节”的节高/孔径比的增加,平均努塞尔数单调增大。但以增大节高/孔径比的途径来提高平均努塞尔数时,阻力系数迅速增大;（3）可通过增大雷诺数和优化节距/节高比来提高竹节孔的平均努塞尔数,计算表明最佳节距/节高比约为10。     

涡轮叶片内部沿周边分布的竹节孔换热特性

根据发动机涡轮叶片内部沿周边分布的竹节孔冷却通道结构建立了简化的传热分析模型.首先通过无量纲分析得到了竹节孔的平均努塞尔数不仅与雷诺数和普朗特数有关,还与"竹节"形状、通道外壁形状以及流体和固体的导热系数比等参数有关.然后,对雷诺数和冷却通道几何参数对竹节孔平均努塞尔数的影响进行了数值模拟研究,计算并对比了上述参数对平均努塞尔数的影响规律,结果表明:①"竹节"结构使平均努塞尔数显著增大,而通道外壁形状主要影响了局部热流密度和局部努塞尔数的分布情况;②随着"竹节"的节高/孔径比的增加,平均努塞尔数单调增大,阻力系数也随之增大;③可通过增大雷诺数和优化节距/节高比来提高竹节孔的平均努塞尔数,计算表明最佳节距/节高比约为10.      

液体火箭发动机推力室冷却通道流动与传热数值研究

采用气固耦合算法对液体火箭发动机推力室再生冷却通道的流动与传热过程进行了三维湍流流动与传热数值模拟，冷却工质为氢气，其密度、导热系数、动力粘度随着温度和压力而变化。应用大涡模拟及标准k-ε双方程模型两种湍流模型分别进行数值模拟，详细揭示了再生冷却通道固体区和流体区内的速度场和温度场，并在不同的计算网格数目下对两种湍流模型的计算结果进行了对比。结果表明，在相同的网格条件下，标准k-ε双方程模型与实验数据的吻合精度比大涡模拟模型更好，且满足工程计算精度。随着网格数的增加，大涡模拟的计算精度逐渐得到改善。     

绕任意轴旋转的管道流动和换热特性研究

当管道绕任意轴旋转时，管道内的轴向速度分布发生改变，截面上产生二次流动，管道的对流换热特性将引起变化。本文对绕任意轴旋转的管道流动与换热进行了数值模拟，研究结果表明：绕任意轴旋转管道内的对流换热特性主要由Re（雷诺数），Ro(Rossby数），Pr（普朗特数）数，和α（旋转轴与管道轴线的夹角）四个参数决定。文中分析了在不同情况下流场结构与摩擦系数比以及温度场与管道Nu数比率随Re数和Ro数的变化规律，为工程实际应用提供了理论依据。     

超临界压力下低温甲烷的湍流传热数值研究

通过系统的数值模拟计算,在准确确定甲烷的热力学和传输物性变化的情况下,详细分析了两种热流密度下超临界压力对低温甲烷的湍流传热过程的影响,揭示了对流换热Nusselt数的变化规律。计算结果表明:在超临界压力下,热力学和传输物性对湍流传热现象会造成很大的影响,尤其在甲烷的临界区域附近,由于物性的剧烈变化会导致传热过程的恶化现象;在高热流密度情况下(如7MW/m2),增大管内压力有利于提高对流换热强度;现有的常用变物性湍流传热公式不能适用于超临界压力下低温甲烷的对流换热计算。     

空腔流动特性数值模拟研究

本文采用数值方法模拟研究了空腔的长深比（L／D）对流动特性的影响。控制方程采用Navier-Stokes方程，计算方法采用有限体积法，湍流模型采用Baldwin-Lomax模型。通过分析三种典型空腔的流线图谱、压力云图、底部压力分布曲线，得到空腔流动特性与其几何参数和来流参数间的相应关系。从数值模拟的结果来看开式空腔流动的空腔内压力分布均匀，流场结构相对稳定，有利于内埋式武器安全平稳地分离。     

地效飞行器大迎角近地面飞行分离流动数值研究

 本文利用Fluent软件求解定常可压缩流动的质量加权平均N-S方程和标准 湍流模型，数值模拟地效飞行器大迎角近地面飞行的全机流场，研究地面效应引起的分离流动。研究结果给出了大迎角近地面飞行机翼上翼面分离流动随飞行高度的变化规律，分析了飞行高度降低引发流动分离的原因，阐述了流动分离对气动性能的影响规律。     

瑞利-伯纳德热湍流中的多湍流态现象

热湍流现象是自然界以及生产生活中很普遍的现象，它源于温差导致的密度差引起的自然对流。瑞利-伯纳德（Rayleigh-Bénard,RB）热湍流系统是从众多自然现象、工程实际中抽象出来的研究热湍流的经典流体力学模型系统。在过去的几十年间，瑞利-伯纳德热湍流系统被广泛用来研究湍流热输运、流动结构演化等特性。最近的研究发现即使在相同的控制参数下，RB热湍流系统中的大尺度流动存在不同的湍流态（结构），且不同的湍流态对应着不同的热输运效率。因此，理解不同的湍流态（结构）的产生机制有助于调控系统的热输运效率。本文就RB热湍流系统中的多湍流态现象,着重评述近年来的若干研究新进展,并对今后的研究工作进行展望。     

无轴向通流对转轴间的流动与传热特性

为研究航空发动机对转轴间的流动与换热,建立了无轴向通流的对转环腔湍流边界层的修正卡门三层结构速度分布模型。基于该速度分布模型并依据动量传递与热量传递的比拟理论,推导了外轴内壁面换热努塞尔数的表达式及准则关系式。与CFD数值计算的速度分布、外轴内壁面换热努塞尔数进行对比验证,结果表明,对转环腔湍流边界层的修正卡门三层结构模型给出了一致的速度分布和换热规律。理论与CFD计算得到的外轴内壁面换热努赛尔数在低旋转雷诺数下符合较好,高旋转雷诺数下偏离不超过20%。     

湍流模型在单/双缝引射矢量喷管中的验证分析研究

利用数值模拟方法,分别求解以SA、SSTk-ω、EASMk-ω或k-ζ湍流模型封闭的RANS方程,针对单/双缝引射的激波控制射流推力矢量喷管展开研究,在多个主喷流压比NPR(4.6,7.0,8.78,10.0)和次主流压比SPR(0.7,1.0)下,系统考察了四种不同湍流模型对射流推力矢量喷管内外流场的流场结构,性能参数及主喷管内壁面压力分布的预测能力,探讨了由激波控制产生的推力矢量随不同参数的变化规律。数值计算结果表明,四种湍流模型都能比较准确地预测出射流推力矢量喷管的性能参数和喷管流场结构,双缝引射的推力矢量效率并不一定优于单缝引射。     

变物性与热辐射对发汗冷却过程的影响规律研究

采用两层k-ε湍流模型对高温条件下有发汗冷却时的矩形槽道内湍流流动和换热进行了数值模拟。计算结果表明:随着冷却气体流量的增加,在发汗冷却壁面附近主流气体的温度梯度大大减小,壁面温度、局部对流换热系数都迅速下降;变物性条件下同时考虑热辐射时数值计算得到的平均St相对值St/S_t0随修正注入率F/St0的变化比已知的常物性条件下的值小;随着高温气体温度的上升,变物性和热辐射对与高温流体接触的表面温度有较大影响,热辐射对壁面温度的影响比变物性更大,而对流换热系数所受的影响不大;当主流为高温气体(如:2500K)、而壁面温度较低时(如:800K)热辐射的影响不大,发汗冷却效果仍然很明显。      

小流量下转静系盘腔传热特性

采用热敏液晶(TLC)瞬态传热实验技术对小流量下的中心进气转静系盘腔传热特性进行研究,介绍了同步触发拍摄系统,考察了流量与转速变化对动盘面传热的影响.结果表明:在所研究工况范围内,流量变化对动盘面中低半径处传热影响较大,对高半径处影响较小;转速升高使动盘面传热效果整体增强;当湍流流动参数较小时,动盘面中高半径区域努塞尔数随半径增大而增大;采用湍流流动参数对动盘面平均努塞尔数进行归纳,实验结果与雷诺相似分析结果的偏差小于15%,一定程度上验证了雷诺相似对转静盘腔传热问题的适用性.