应用改进的低雷诺数湍流模型预测转捩流动

基于Launder-Sharma(LS)低雷诺数k-ε两方程湍流模型发展了一种改进的具有转捩敏感性的低雷诺数湍流模型.针对LS模型的可实现性(realizability)问题和前缘滞止点湍动能预测过大的不足,模型进行了改进.改进模型只使用当地物理量,不需要求解壁面距离、y+和边界层积分参数.改进模型能够适用于广泛的流动,且容易应用到通用的计算流体动力学(CFD)程序中.对具有详细数据的零压力梯度平板转捩边界层T3A实验的模拟结果显示,改进模型能够预测转捩流动,并能对自由湍流变化给出合理的响应.      

气膜冷却恢复区的湍流模型

在 L am- Bremhorst低雷诺数 k- ε模型的基础上 ,发展了气膜冷却恢复区的湍流模型。原模型中的常数 cμ 和 c1代之以沿 x、y方向变化的参数 ,低雷诺数函数 fμ、f1和 f2 沿 x方向有所变化。修正后的模型可以在 0 .4～ 1.0的吹风比范围内 ,较好地描写恢复区内时均速度、湍动能和湍应力分布 ,并能正确模拟向不吹风情况过渡的恢复过程。     

用非线性k-ε-k_p两相湍流模型模拟旋流两相流动

建立了两相湍流的代数应力模型 ,并由此出发 ,导出非线性k ε kp 两相湍流模型 ,目的是合理地模拟各向异性较强的旋流两相流动 ,保持二阶矩模型的优点 ,同时比二阶矩模型简单 .文中得到了气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k ,kp,以及两相脉动关联湍动能kpg的方程联立 ,就构成非线性k ε kp 模型 .将该模型用于模拟旋流两相流动 ,给出两相时均速度场及雷诺应力各分量 ,并且将模拟结果和实验数据以及二阶矩模型的模拟结果对照 .研究结果表明 ,该模型预报旋流两相流动的能力和二阶矩模型的能力相差不多 ,但计算量比二阶矩模型的小     

扩缩管充分发展层流的强化传热场协同分析

采用二维轴对称不可压缩层流N-S方程,利用有限体积方法和SIMPLE算法,在雷诺数100～1300范围内,对轴对称周期性扩缩管的充分发展层流流动和传热特性进行了数值模拟,重点研究了不同雷诺数、几何外形与换热强度的关系,并利用场协同理论对计算的结果进行了分析,获取了层流态下扩缩管对流换热强度随雷诺数、几何尺寸变化的规律,对进一步研究扩缩管换热问题、探讨强化传热理论和发展换热技术均有实际意义.     

用非线性k-ε-kp两相湍流模型[1*9/9]模拟旋流两相流动

建立了两相湍流的代数应力模型,并由此出发, 导出非线性k-ε-kp两相湍流模型, 目的是合理地模拟各向异性较强的旋流两相流动, 保持二阶矩模型的优点, 同时比二阶矩模型简单. 文中得到了气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式.这些代数式和两相各自的湍动能k, kp, 以及两相脉动关联湍动能kpg的方程联立, 就构成非线性k-ε-kp模型. 将该模型用于模拟旋流两相流动, 给出两相时均速度场及雷诺应力各分量,并且将模拟结果和实验数据以及二阶矩模型的模拟结果对照. 研究结果表明,该模型预报旋流两相流动的能力和二阶矩模型的能力相差不多,但计算量比二阶矩模型的小.     

矩形横截面螺旋管湍流流动与传热特性的数值研究

 采用了剪切应力输运(SST)k-ω两方程湍流模型并考虑近壁低雷诺数的修正对矩形横截面螺旋管内冷却水流动和传热特性进行了数值研究。数值分析了在不同入口雷诺数、曲率半径以及扭距条件下,螺旋管内的温度、速度场以及流线的变化,讨论了螺旋管内、外壁面对流传热系数的差异及产生机理,同时与直通管道传热性能进行了比较。研究发现由于离心力的作用,螺旋管内存在显著的二次流动,管内、外侧壁面对流传热存在差异。旋转一周后,螺旋管即进入了流动稳定状态,入口雷诺数可以显著提升螺旋管整体的对流换热效率,扭矩和曲率对内外壁面传热效果的影响不大,而窄高型的横截面构型可以显著改善螺旋管的传热效果。研究结果对应用矩形横截面螺旋管的冷却设计提供参考。     

矩形与收敛通道内高实度扰流柱换热特性研究

采用SST k-ω湍流模型,对矩形与收敛两种通道内高实度(45%)圆形扰流柱的传热与流动特性进行数值模拟。通过比较两种通道内的流场结构,分析其内部冷却气体的流动机理,进而探求扰流柱阵列的传热性能与压降变化。结果表明:通道类型对内部流场结构有很大影响。收敛通道内流体的速度相比矩形通道的呈现沿程增加的趋势,其对自身的流动损失及换热效果影响加大;两种通道内扰流柱的平均换热水平随进口雷诺数的增大均呈指数上升趋势,相比较收敛通道的更高,但差距逐渐缩小;两种通道的整体压力损失系数均呈指数下降趋势。同一进口雷诺数下,矩形通道在中游出现换热峰值,收敛通道的换热效果沿流向持续提升,两者虽在相同排列产生换热波动点,但单排换热差距逐渐增大。矩形通道内单排扰流柱压力损失系数沿流向先降低后升高,收敛通道内则持续上升。     

应用一个新的具有转捩敏感性两方程低雷诺数湍流模型预测转捩流动

本文基于Launder-Sharma(LS)低雷诺数的 两方程湍流模型发展了一种新的具有转捩敏感性的低雷诺数湍流模型。针对LS模型的可实现性(realizability)问题和前缘滞止点湍动能预测过大的不足,新模型进行了改进。新模型只使用当地物理量,不需要求解壁面距离、Y+和边界层积分参数。新模型能够适用于广泛的流动,且容易应用到通用的CFD程序中。对具有详细数据的零压力梯度平板转捩边界层T3A实验的模拟结果显示,新模型能够预测转捩流动,并能对自由湍流变化给出合理的响应。 关键词：航空、航天推进系统;转捩模型;转捩流动;数值模拟     

部分平均Navier-Stokes方法模拟超声速斜面空腔流动

采用部分平均Navier-Stokes (PANS)方法对超声速斜面空腔流动进行了数值模拟研究,评估其对于超声速非定常湍流模拟的性能,并与雷诺平均N-S(RANS),脱体涡模拟(DES)的计算结果及实验数据进行对比.研究表明:模化湍动能比例全场可变的PANS方法预测的速度、壁面压力和摩擦力系数分布与DES的结果非常接近,都与实验值吻合得较好,优于RANS的计算结果;非定常雷诺平均N-S(URANS)计算的自由剪切层近乎为二维稳态的,而DES和PANS方法可以求解出更为丰富的流动结构;模化湍动能比例全场统一的PANS方法虽然可以求解出相比RANS更多尺度的流动结构,但在近壁区不能回归到RANS模型,预测的湍流边界层的速度型偏离对数律,后续的流动计算也偏离实验数据.      

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

高位垂直进气径向出流旋转盘腔换热的实验研究

用实验的方法对高位垂直进气转静系旋转盘附近冷气换热特性的基础数据进行补充和拓展研究,建立了高速旋转换热实验台,并运用遥测技术采集温度信号,得到了转盘表面温度、局部努赛尔特数的分布及平均努赛尔特数,主要结论为:(1)转盘表面温度呈外高内低分布.(2)转盘表面的局部努赛尔特数在r＜0.4b的区域基本不变,r＞0.4b的区域里,局部努赛尔特数随半径的增加而迅速增大.(3)其它条件不变时,随流量系数的提高,平均努赛尔特数大幅度上升.     

基于内外传热耦合的热气防冰系统仿真计算

基于内外传热耦合原理,建立了热气防冰系统的AMESim仿真模型,研究热气防冰系统在不同引气状态下,管路流量、压力及蒙皮温度的变化规律.采用所建立的仿真方法,计算飞机机翼热气防冰系统的内部流动特性和内外耦合传热特性,将计算结果与热气防冰系统流量分配试验结果进行对比.结果表明:管路流量、压力和蒙皮温度的仿真计算结果与试验结果最大误差为9%,验证了该仿真模型的正确性.在此基础上,对飞机短舱热气防冰系统进行了仿真,分析了飞行包线下系统内外的换热、温度变化、加热效率等关键参数的瞬态特性.仿真结果为热气防冰系统的设计、分析与优化提供依据.      

表面凸起对机翼热气防冰腔内换热强化的影响

采用数值模拟对比研究了光滑表面和具有表面凸起结构热气防冰腔内湍流流动的换热特性。机翼防冰腔内笛形管具有三排射流孔,射流孔角度有0°±45°组合以及0°±30°组合。为了强化射流冲击光滑表面的流动换热,在防冰腔内表面正对射流孔的射流冲击区,设计了表面凸起结构,用来强化射流对壁面的冲击换热效果并起到引流作用。通过改变射流孔射流角度研究了射流角度对传热特性的影响。计算结果表明:与光滑防冰腔内表面射流冲击换热相比,表面凸起结构可以将均匀发散的壁面射流集中为高速壁面射流,提高壁面射流区的对流换热系数,从而增强射流冲击换热效果,机翼前缘的强化换热效果尤为明显。     

气冷涡轮叶片三维换热问题计算

建立了三维气冷涡轮叶片换热计算的热-流耦合计算模型, 将辐射热量作为源项加入到方程中, 采用热-流耦合方法完成了某气冷涡轮叶片带辐射的三维换热计算.辐射热量的计算采用了P-1辐射模型, 比较了计算辐射和不计辐射时叶片表面的温度值及温度分布的差异.计算表明对于1600 K这样的高温燃气, 辐射热量不可忽略, 考虑辐射热量后, 叶片温度明显升高, 温度分布也产生了显著差异.      

基于流-固耦合传热的热气防冰系统干空气飞行蒙皮温度场计算研究

以机翼热气防冰系统为研究对象,建立了包含热气防冰系统防冰腔内外流场对流换热和固体结构导热的三维稳态流-固耦合传热物理模型,对整个计算区域生成混合网格,边界条件为第三类边界条件,采用计算流体力学方法以 FLUENT 软件为工具,对干空气飞行状态下流-固耦合传热模型进行了求解,获得防冰腔蒙皮内外表面对流换热系数分布和温度场结果,并对计算结果进行了分析。结果表明：防冰腔铝合金蒙皮沿展向和厚度方向导热显著,温度分布较均匀,防冰引气温度为200℃时,防冰腔蒙皮内外表面上最高温度为101℃,最低温度为21℃,3 mm厚的蒙皮同一点处内外表面最大温差仅为4℃,防冰腔排气口处气体的平均温度为63℃。热气防冰系统蒙皮温度场计算方法和计算结果,能够为热气防冰系统干空气飞行试验设计和测试中温度传感器的选型与布置提供依据。     

燃气热力参数计算方法及在涡轮叶栅对流换热中的应用

为了研究燃气热力参数对对流换热的影响,发展了一种计算航空发动机碳氢燃料与空气燃烧的燃气热力参数计算方法。应用该方法给出的一部分燃气热力参数作为输入条件,对不同入口温度下的平板和涡轮平面叶栅的对流换热进行了数值模拟,并与空气作为工质的对流换热数值模拟结果进行了对比。结果表明:在其他定解条件相同时,对平板和涡轮平面叶栅的对流换热,工质为多组分燃气时的壁面当地努塞尔数均大于入口工质为空气的努塞尔数,且随着入口温度的升高,空气与燃气的当地努塞尔数的差别变大;在给定的入口温度范围内,与工质为空气相比,工质为燃气的平板当地努塞尔数从28%增大到30%,涡轮叶栅则从725%增大到969%。     

振荡管流换热的解析解及强化换热特性

经求解随时间正弦变化的非定常压力梯度下,管内流体振荡时的运动方程和能量方程,得到了振荡管流换热的速度场、温度场以及当量导温系数的解析解.通过分析可知,影响管内振荡流轴向换热的无量纲参数主要有:无量纲振荡频率,无量纲振幅和流体普朗特数.其中,无量纲振荡频率直接影响流动的速度分布,且随着无量纲振荡频率增加,强化换热效果增强.随着无量纲振幅的增大,管内流动有效导热面积增加,等温面法向温度梯度也增加,使得强化换热效果增强.强化换热效果同样也随流体普朗特数的增加而增强.      

Effects of flow parameters on thermal performance of an inner-liner anti-icing system with jets impingement heat transfer

The multiple jets impingement heat transfer is widely applied in the wing anti-icing system. It is challenging to apply the similarity criterion to carry out the anti-icing experiments due to the complex flow and heat transfer behavior. In the present study, the full-scale slat model is used to carry out anti-icing experimental researches in a 2 m × 3 m icing wind tunnel of China Aerodynamics Research and Development Center. The effects of icing parameters Liquid Water Content (LWC) and Median Volume Diameter (MVD) and hot air parameters (mass flow rate and temperature) on the thermal performance of an inner-liner anti-icing system with jets impingement heat transfer are studied. The effects of the experimental parameters are analyzed in detail by combining impingement and evaporation heat transfer mechanisms. The impingement hot air mass flow rate dramatically affects the heat transfer performance of the impingement stagnation region within the range of the experimental parameters. The temperature of impingement hot air and that of wing skin are approximately linear correlated. The experimental results show the effects of LWC and MVD on water film formation and runback ice accretion. The formation of water film is analyzed by an analytical method based on the wing skin temperature difference of dry and wet air conditions.     

Energy-efficient intensifiers of laminar heat transfer

In this paper, calculations of some thermo-hydraulic parameters of different intensified channels have been carried out and a grooved channel that is optimal in mass and energy expenditure is determined. A new variant of the intensified channel is proposed and recommendations for engineers are presented.     

Helmholzt振荡腔对光滑圆管换热的强化实验

气流通过Helmholzt振荡腔产生脉动气流.实验固定振荡腔进出口直径和振荡腔直径, 在改变振荡腔腔长的情况下, 在圆管进口Re为22000到82000范围内对等热流加热条件下的圆管换热实验进行了研究.实验将脉动气流所产生的圆管换热效果与定常气流产生的换热效果进行对比.结果表明:气流流经Helmholzt振荡腔后对圆管换热有强化作用, 其强化作用主要集中在沿圆管轴向前10倍直径的长度上, 强化比约为1.1-1.9.随着距离增大换热强化比迅速减小.振荡腔腔长对换热效果有较强的影响.