二维流热耦合数值模拟程序及方法研究

本文发展了一套二维流热耦合数值模拟程序,流体域中采用结构化网格,固体域中采用非结构化网格,流体域和固体域耦合面上的网格为非匹配连接。为了使该程序能够同时模拟低速不可压流动,程序中还加入了预处理方法。此外,本文还发展了一种能较好保持通量守恒的界面插值技术,来处理流体域和固体域耦合面上的非匹配网格。数值计算结果与实验值或商业软件结果的对比表明:本程序中的预处理方法有效地改善了低速不可压流动下的计算收敛性和收敛精度,使本程序能够同时模拟高速和低速不可压流动;扩压器和叶栅的流热耦合计算表明本文程序在流热耦合计算方面具有较好的模拟能力,同时也表明本文发展的界面插值技术的可靠性。     

基于自适应采样的高超声速飞行器气动热全局快速预示

高超声速飞行器热防护系统设计中高精度气动热分析模型使得设计计算成本不断增加，基于数据驱动的气动热环境预示方法受到广泛关注。针对有限高精度模型计算成本下提升全局预示精度的问题，提出一种基于模糊聚类的批量自适应采样方法。根据预示误差分布特征通过聚类采用超球分割构建采样影响域，兼顾误差较大的重点采样域与全局探索；通过当地误差评分系数加权构建采样拒绝域，减小新增样本冗余；结合maxmin准则在综合确定的重点采样空间中新增样本，提升采样质量，进而实现预示模型全局精度快速提升。数值测试算例表明，所提方法与One-Shot、APSFC、CV-Voronoi方法相比能有效降低所需采样规模，加速提升预示精度。通过类HTV-2飞行器气动热快速预示实例，验证了方法的有效性与工程实用性。     

回流燃烧室流动特性试验

为了揭示有/无燃烧状态下燃烧室热态和冷态流场的特征和流动特性,针对某型回流燃烧室单头部试验件,使用粒子图像测速仪（PIV）,测量燃烧室燃烧状态下不同截面处的热态流场,以及没有燃烧状态下不同截面处的冷态流场,探讨不同总压损失系数对回流燃烧室热态/冷态流场特征及流动特性的影响。研究表明:随着总压损失系数的增大,冷态条件下各截面流场结构基本保持不变,如射流孔穿透深度、射流角度、回流区位置及大小、流线等基本保持一致,但是各位置点速度大小逐渐增大。热态条件下各截面流场随着总压损失系数增大,流场结构也基本保持不变;相同总压损失系数时,热态流场与冷态流场存在差异,燃油喷射与气流的相对运动将会对燃烧室头部的流场结构造成影响,速度较冷态流动时略微增大。      

热气泡微机械加速度计

本文论述了一种以热气泡为主的微机械加速度计,技术上不同点在-t-]s~速度计中唯一可移动的元件是一个小的热气泡,在微型舱体中用高磁通量的加热器蒸发液体产生了该热气泡。分析和讨论了该传感器的基本物理特征,包括热量传递和液体流动。用数字仿真验证上述加速度计特性,并用特殊的测试设备演示该性能,其结果是响应良好,灵敏度高。     

球面收敛二元扩张矢量喷管热射流特征的实验研究

为了获得球面收敛二元扩张矢量喷管热射流特征，在小型热射流实验台架上运用红外热像显示和流场测量，在总温673～823K、喷管落压比(NPR)1．46～2．25状态下，对不同喉道宽高比、俯仰角和偏航角的热射流特征进行了实验研究。在本研究范围内，NPR对矢量喷管上壁面附面层流动影响不大；在扩张上壁附面层分离的情况下．壁面附近的流线被上抬，形成“马鞍”形异型总压分布；喷121出121下游的激波呈现一连串的“葫芦”状，随喷管矢量角的增大，流场的不对称性越明显；在收扩喷管喉道宽高比（ARi）为2时，热射流被分成两股。     

碳氢燃料热声不稳定流动特性和边界实验分析

为保障空-油换热器安全稳定运行,对竖直上升圆管内超临界压力航空燃料的热声不稳定流动进行了实验研究。阐述了管壁温度的振荡特征和热声不稳定流动的诱发原因。考察了运行参数对流动稳定工况和流动不稳定工况界限的影响。取相邻流动稳定工况和流动不稳定工况的热流密度平均值作为边界,分析了该边界随进口拟过冷度的变化情况。通过临界拟相变数和拟过冷度数的比值表征,建立了热声不稳定流动边界与进口雷诺数和相对压力关联的预测准则。结果表明:拟沸腾效应导致了热声不稳定流动现象。进口温度和运行压力越低,热声不稳定流动越容易出现。提出的边界预测准则具有较高的精度,与实验结果的相对偏差在±15%以内。     

离心力场下细微封闭通道内R12的热驱动换热研究

本文以试验的方法,研究了一种利用彻体力场下封闭通道内流体的热驱动循环流动来强化换热的新型冷却技术。主要讨论了封闭循环通道宽度d为1mm,氟利昂R12为工质的热驱动流动和换热特性,重点分析了随转速(即离心加速度)的变化,R12的热驱动流动和换热变化规律,以及热端壁面的温度分布及变化规律。研究结果表明:离心力场下流体热驱动流动和换热特点同重力场下基本相似;随着转速的提高,R12的热驱动强度不断增加,热端壁面温度逐步下降,同时热端壁面的温度梯度也逐渐降低;在转速为1140r/min时,总平均传热系数达到了1550W/m2·K。      

乘波体特性研究现状及展望

乘波体特性研究是开展乘波飞行器设计的基础。本文对乘波体相关特性的研究现状进行了分析，重点总结了乘波体几何特性、气动特性、气动热特性等静态特性和横航向稳定性、纵向稳定性以及弹性效应对稳定性的影响等动态特性的研究进展；在此基础上，对乘波体特性研究中还存在的困难和问题进行了分析，并对未来进一步的研究重点进行了展望。研究表明，乘波体的几何特性决定了其适用于高超声速飞行器的前体设计；乘波飞行器适用于在设计点和小迎角工况下飞行，宽速域、大空域范围内气动性能还难以满足实际使用需求；在进行高超声速乘波飞行器气动外形设计时，应综合考虑气动力和气动热之间的平衡折中，采用合适的气动热解决方法来满足工程应用的需要；乘波体需具有上反角和“下凸”的外形特征，以满足稳定性要求。可为乘波飞行器的研究发展提供参考。     

阵列射流冲击冷却局部对流换热特性的数值计算与实验研究

运用数值计算的方法对具有初始横流的阵列射流在不同的排列方式、不同的冲击间距和不同的横流/射流质量流量比下的流动换热进行了三维数值研究,并采用热色液晶测试技术对阵列射流冲击的冷却表面温度分布进行了试验研究,获得了每一股射流的冲击冷却局部对流换热系数分布的特征。研究结果表明本文的计算结果与实验特征是基本吻合的。     

高温化学非平衡湍流边界层脉动量象限分析

高超声速飞行器在较低空域以极高马赫数飞行时，表面会同时存在湍流与化学非平衡流动，但目前针对此类高温化学非平衡湍流边界层流动特性的研究工作还比较有限，对不同湍流特征的主导流动机制的认识还有待于进一步深入。选取高超声速楔形体头部斜激波后的流动状态，设置3种不同的壁面温度，通过直接数值模拟对比了不同壁温条件下的边界层参数分布特性，并采用象限分解技术分析了边界层不同象限流动事件对雷诺剪切应力、湍流热流、湍流质量扩散的贡献。结果显示：在整个边界层中上抛和下扫运动对雷诺剪切应力的贡献占优；冷壁效应会使得流向和法向湍流热流通量的主导流动事件在温度峰值两侧发生改变。O原子组分流向湍流组分扩散主要受到高质量分数流体慢速运动事件和低质量分数流体快速运动事件的影响，而法向湍流组分扩散则主要受到高质量分数流体向上运动事件和低质量分数流体向下运动事件的影响。     

微小交错流道散热块传热性能

实验研究微小交错流道散热块的传热性能以及热阻系数.随着Re数增加,交错流道与直通流道内流动换热性能均会升高,前者的变化更加显著,对流热阻系数与Re数接近指数关系;相同Re数条件下,交错流道对流热阻系数约比直通流道下降18%.散热块底面温度分布明显不一致,流体轴向温升对温度分布影响显著,温度在空间与时间上均有波动;预测交错流道在Re数为700左右开始发生流动转捩.      

潜热型功能热流体在扁管内的对流换热实验

通过质量分数为5％～2o％的潜热型功能热流体在扁管内的对流换热特性实验,研究了相变材料微胶囊质量分数、质量流量等因素对流体换热性能的影响.实验发现,流量较小的情况下,功能流体传热性能优于水且随着质量分数的升高而增强,而当流量较大的情况下,接近入口处的功能流体的传热性能则弱于水.以改进的Nusselt数评价,入口段内,质...     

Application of Coolant Diffusion Equation in Numerical Simulation of Flow in Air-cooled Turbines

This article deals with time-averaged coolant diffusion equation in arbitrary, body-fitted coordinates. Schmidt number and Lewis number are used to obtain the mass diffusion coefficient. Other than the existing N-S solver, a new solver for the diffusion equation and the influences of cooling air on the thermal parameters, such as gas constant, constant-pressure specific heat, is established. Capable of simulating the flow of multi component fluid, this solver is used to simulate the flow in air-cooled two-stage turbines. The comparison of test results from the modified solver to those from the traditional N-S solver has evidenced that a marked difference present is in thermal parameters and flow field structure.     

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

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

考虑发动机冷却通道固壁内耦合导热影响的低温甲烷超临界压力传热研究

 在考虑发动机冷却通道固壁内耦合导热影响的情况下,开展了低温甲烷在矩形冷却通道中的超临界压力湍流换热数值模拟研究;仔细分析了热流密度及管道几何形状对低温甲烷超临界压力下的流动和传热的影响;得到了流体速度、壁面温度、壁面热流密度等参数的详细变化情况以及Nusselt数的变化规律。计算结果表明:在考虑流固耦合作用时,上壁面施加的热流有一部分会通过固体壁面内的热传导,经由侧壁面传入超临界压力流体,并且随着热流密度的增加,经侧壁面传导的热流所占的比例也会随之增大;减小冷却通道内截面的高宽比,可以提高超临界压力下的换热效果,但流动压降会大大增加,因此冷却通道高宽比的选择需综合考虑传热与压力损失的影响,可以引入热性能参数作为参考;修正的Jackson&;Hall对流换热关系式基本可以适用于本文中的各种工况。     

涡产生器倾角对换热器内流动与换热特性的影响

为优化涡产生器强化传热的性能,通过改变三角翼式涡产生器的倾角,数值分析了倾角变化对板式换热器通道内流动与换热特性的影响,倾角变化范围为-30°～30°。获得了不同倾角下通道内的纵向涡、涡量分布,换热及阻力特性。结果表明涡产生器倾斜方向对流动与换热的影响不同,存在最佳的涡产生器倾角-10°,使得换热通道内涡产生器强化换热的性能最优;改变倾角可进一步提高涡产生器的强化传热性能,不同倾角对应努塞尔数和热性能评价因子间的最大差值分别达到8.9%和5.7%。     

重力场下细微循环通道内流体热驱动现象研究

本文以试验结合数值模拟的方法,研究了一种利用彻体力场下封闭通道内流体的热驱动循环流动来强化换热的新型冷却方式。由显示实验清楚地拍摄到封闭通道内水的热驱动循环流动,同时通过实验和数值分析研究了水的运动和换热特性,实验和计算结果基本一致。研究表明在本文研究的参数范围内,只有当加热量超过30000W/m2时,流体的热驱动力才能克服阻力影响来驱使流体运动,同时加热端的平均换热系数达到了2500W/m2·K。此换热结构在高密度封装电子器件的冷却设计中将会有良好的运用前景。      

微通道热沉对流传热理论模型及实验

用理论及实验相结合的方法研究了微通道热沉流动与传热特性.首先,总结并提出了微通道热沉对流传热的理论模型;然后,实验测量并计算了微通道热沉的压降及努塞尔数,其理论值与实验值吻合较好,平均误差在10%左右;最后,分析了不同雷诺数及通道宽高比时的导热热阻、对流热阻及电容热阻占总热阻份额的大小.结果表明:对流热阻是影响微通道热沉传热性能的重要因素,当雷诺数为985,通道宽高比为1时,对流热阻占总热阻90%左右;而在雷诺数较小时,导热热阻占总热阻的份额小于10%,可以忽略不计;电容热阻占总热阻的份额随着雷诺数及通道宽高比的增大而降低.      

旋转封闭循环小通道内液态金属热驱动换热特性研究

高速旋转条件下,封闭循环小通道内的流体(热驱动介质)在浮升力的作用下会产生热驱动运动,将热量不断的从高温端携带至冷却端释放。本文通过数值模拟的方法,研究了离心力场下封闭循环小通道中,不同性质流体的流动和换热特性。重点分析了以液态金属锂、钠和钠钾合金为热驱动介质时,热驱动运动的换热特性Nu随Bu和Ro两个参数的变化规律。研究中发现在一定的Bu数变化范围内,以液态金属锂、钠以及钠钾合金为热驱动介质所得到的Nu均随Ro数的增大而减小;当Ro数相同时,Nu随Bu数的增加而变大。研究结果表明高转速,低加热量时,采用液态钠钾合金为热驱动介质可以得到最好的热驱动换热效果;在高转速,高加热量时以液态锂为热驱动介质时,热驱动换热能力最强。      

不同进气方式对机载电子设备气冷冷板性能影响

电子设备的热负荷占到了航空器全机热负荷的60%以上,高温失效是电子设备失效的主要形式,采取合理的散热手段及时带走电子设备产生的热量对飞行安全和电子设备工作的稳定性有着至关重要的影响。文章利用计算流体力学方法(CFD)对某机载电子设备用气冷冷板进行了传热模拟,通过改变入口速度以及冷空气的流动方式,得到不同工况下气冷冷板各位置的温度分布和工作性能。研究结果显示,入口速度增加会使冷却效果得到提升,但同时增加了流动阻力使得设备工作的稳定性下降。另外,流动方式一(中间进两侧出)的冷却效果优于流动方式二(两侧进中间出),流动方式二的工作噪音和流动阻力优于流动方式一,可根据实际需要进行选择。