轴向通流旋转盘腔流动换热稳态实验与数值模拟

通过开展稳态实验及数值模拟探究了轴向通流旋转盘腔的流动结构与换热特性。通过改变流量系数、旋转雷诺数等参数,探究不同工况下旋转盘两侧及盘罩内侧壁面温度和努塞尔数的径向分布规律。结果表明:在轮缘加热的状态下,旋转盘两侧温度径向分布均呈凹函数形态,且旋转盘迎风面换热强度普遍高于背风面;后轴颈盘罩向两端旋转盘导热,其壁面温度径向分布呈"中间高、两侧低"的状态;随着轴向流量系数的增大,盘腔内部气体对流加剧,径向臂及涡对结构更加明显,旋转盘及轴颈表面换热效果增强;旋转盘腔内的流动换热特性受强迫对流和类Rayleigh-Benard对流2种机理的共同影响。      

轴向通流旋转盘腔内流动不稳定性研究

为了研究旋转腔内流动不稳定性问题,采用旋转坐标系稳态方程,用数值模拟的方法分析了轴向通流旋转盘腔内的流动,得到了旋转系下哥氏力和离心浮升力分别对流动不稳定的作用,以及这两个力的非线性综合作用对流动不稳定的影响.结果表明:哥氏力是恢复力,不影响流动稳定性;离心浮升力是造成流动不稳定的主要因素,盘腔内的流动是离心力场下的Rayleigh-Benard对流和强迫对流的混合流,随着浮升力的增强流动由稳态发展为非稳态;哥氏力与离心浮升力的综合作用加剧了流动不稳定性,盘腔内r-θ面出现了明显的旋向相反的对涡.     

旋转光滑U形通道内流动和换热的数值模拟

数值模拟了旋转状态下涡轮叶片U形内冷通道湍流流场和温度场的分布,分析了流阻和换热的变化规律.结果表明,旋转状态下哥氏力、离心力和浮升力的共同作用使得流场发生了复杂的变化.旋转强化了换热,减小了流阻.但旋转使得换热在各个面换热能力分布不均,增加了温度梯度.      

涡轮叶片尾缘复合冷却通道换热的数值模拟

用数值模拟的方法对一种典型的涡轮叶片尾缘冷却结构进行了数值计算,该冷却结构由横肋通道和交错扰流柱排构成.通过对不同进气雷诺数的计算结果分析可以看出:对所研究的冷却结构,扰流柱通道顶部及底部这两处冷气流速较高的区域的换热系数也比较高,温度较低;尾缘冷气出流孔处的换热系数也比较高;但扰流柱区的中部则由于冷气量明显不足导致明显的换热系数降低.      

航空发动机一维数值仿真

分析了航空发动机一维气动热力数值仿真模型,此模型能够用于航空发动机整机数值仿真.进行了整机性能数值仿真,求解过程考虑了粘性效应以及气体流动损失对于整机性能的影响;采用了具有TVD性质的高阶精度Godunov格式,此格式具有良好的数值稳定性和激波捕捉能力;应用了时间推进法,交替使用显/隐格式.计算了某型带加力的混排涡扇发动机的速度特性以及节流特性并且与设计数据进行了比较分析,数值仿真结果与设计数据比较吻合.分析了典型工作点的计算过程,模型收敛速度较快.     

带肋变截面回转通道内流动与换热的数值模拟

开发了三维流动换热的通用计算程序,数值研究了带肋变截面回转通道内流动与换热的特性.湍流模型采用低雷诺数k-ε模型.通道肋间距为25mm,肋高分别为1mm,1.5mm,2mm,冷气进口雷诺数Re分别为7500,12500,18500,25000.计算结果表明:①通道的平均努赛尔数均随进口雷诺数的增大而增大;②对于Re=7500和12500,肋高越高,换热越强;对于Re=18500和25000,肋高为1.5mm的通道换热最强;③局部雷诺数的不同和离心力的影响导致通道内各区域的局部换热随肋高的变化趋势并不一致;在进口段,肋高越高,换热越强;在出口段,当Re=7500和12500时,肋高越高,换热越强,而当Re=18500和25000时,存在最佳肋高1.5mm.     

旋转附加力对方通道内流动与换热的影响机理

在雷诺数为25000,旋转数为0~0.24,温度比为0~0.22的范围内,数值模拟了旋转光滑径向出流通道的内流动与换热分布,分析了哥氏力与离心浮升力对旋转管流的耦合作用机理.计算结果表明,切向哥氏力推动了通道截面内的双涡二次流,径向哥氏力则使得近侧壁流体加速和中心流体减速.离心浮升力对流动与换热的作用效果与哥氏力场的分布密切相关.换热计算结果从定性趋势上吻合公开文献中的实验现象,反映了旋转附加力的基本影响规律.     

水下弹体空化流动的数值模拟

采用均质流模型对水下弹体的空化流动进行了数值模拟,控制方程的求解采用有限体积法,对流项采用迎风通量差分格式,时间步的推进则采用了LU-SGS隐式方法.使用质量传输模型模拟汽-液之间的质量转换.为了更好的计算低速流动,在计算中对对流项进行了预处理,并对质量传输源项进行了点隐式处理.对绕圆柱钝头体空化流动进行了模拟,雷诺数Re＝136 000,计算得到了圆柱形钝头体在空化数σ＝0.1,0.2,0.3时的壁面压力系数,并与试验结果进行了对比,符合良好.在此基础上研究了弹体带尾迹的超空化现象,空化数σ＝0.2,0.1,0.06,在σ＝0.06下得到了水下弹体的超空化尾迹.      

旋转状态下涡轮叶片前缘的流动与换热

用数值模拟的方法对旋转状态下涡轮叶片前缘冷却结构进行了数值研究,该结构由进气腔、叶片尾缘块和前缘块构成,对此结构不同的旋转速度情况进行了计算,根据计算结果分析了旋转对涡轮叶片前缘流动与换热的影响.计算结果表明,旋转状态下带气膜出流的冲击流动中,前尾缘冲击面的换热随着转速的增加而减小,且尾缘冲击面的换热比前缘冲击面的换热要好;同时前尾缘冲击面换热的差别随着转速的增加将越来越小.     

上游双扰动对圆柱气动特性影响的数值模拟

通过用商用计算流体力学软件FLUENT 6.0中的SIMPLE方法求解二维非定常层流模型,对雷诺数Re=200时上游双扰动体对下游圆柱气动特性影响进行了数值模拟研究.计算结果表明当双扰动体相隔较近(H/d=1.1,1.33)时,其后形成单一的Karman涡街,它们对下游圆柱的干扰作用类似于单扰动体时的情况,但是流动结构发生变化的临界间距Lcri较之有效宽度相同的单扰动体的要大,因此可以推断出在该种布置下双扰动体和圆柱绕流会出现2种流动模式.在空穴流动模式下面,下游圆柱的阻力和升力脉动最多减小94%和79%,系统阻力(包括圆柱和扰动体)也能减小33%.而当双扰动体尾流为双Karman涡街时(双扰动体相隔较远H/d=3.33,5.0),圆柱迎风面部分区域仍受到来流的冲击作用,气动力减小有限.     

矩形通道中亚尺度肋片的流动换热数值分析

对装有不同结构亚尺度肋片矩形通道的流动和换热进行了数值模拟,获得了通道流场、温度场分布以及平均努塞尔数的基本特征,对各种肋片通道的换热特性进行了对比分析.计算结果表明:亚尺度肋片距主肋片根部越近越有利于散热;在远离热面区域,亚尺度肋片的比表面积越大传热效果越好.当扩展表面积相同时,亚尺度肋片的长宽比越大换热效果越好;计算结果及分析揭示了控制肋片最优几何形状的统一原则——广义温度梯度均匀化原则,并依照此原则确定了可以强化换热的较优的肋片结构.计算分析表明广义温度梯度均匀化原则可以控制对流换热过程的传热强化.      

轴向通流多级旋转盘腔内的稳态换热特性

针对航空发动机多级压气机盘腔的冷却问题,对具有轴向通流的三腔模型进行了系统的研究,通过对局部换热特性与平均换热特性的分析,可以表征旋转盘的热应力水平和相对冷却效果,结果表明表征旋转盘热应力水平的局部换热特性在旋转盘的半径上呈多峰分布,表征相对冷却效果的平均换热特性与旋转雷诺数Re_w、进气雷诺数Re以及格拉晓夫数Gr有关,当Re及Gr增加时,平均努塞尔数 u随之增加,而Re_w的增加使之减小.      

高位垂直进气旋转盘非稳态换热的实验研究

用实验的方法对高位垂直进气的转静系旋转盘在非稳态情况下的流动与换热特性进行了研究,主要研究了冷气流量系数和旋转雷诺数的变化对盘面温度和平均努塞尔数的影响.结果发现:冷气流量的改变对盘面的换热影响非常明显,盘面各点温度随时间的变化率是相同的;旋转雷诺数增加使盘面平均努塞尔数增大,旋转雷诺数变化的方式导致盘缘区域温度随时间的变化率与中心区域的不同.      

涡旋微槽流动与传热特性数值模拟

涡旋微槽散热器具有传输高热流密度的潜力,在解决航空航天高功率密度器件热控制方面具有广泛应用前景.在实验研究的基础上,采用有限体积法对不同体积流量和槽道结构的涡旋微槽中的流动与传热特性进行了数值模拟研究.对涡旋微槽流动的稳定性进行了分析,给出了摩擦因子和Nu数沿流动方向的变化曲线,并采用场协同原理对涡旋微槽强化传热的机理进行了探讨.计算得到的微槽平均传热系数和摩擦阻力系数与实验数据进行了对比.结果表明:涡旋微槽中二次流的出现是涡旋微槽强化传热的机理所在.     

Y形构形微通道流动换热特性的数值分析

针对矩形芯片提出了Y形构形微通道的芯片内部冷却结构,通过三维数值计算比较了Y形构形微通道与传统直微通道中的流体层流流动和换热特性,并对Y形构形微通道的重要结构参数进行了优化.研究表明:在相同对流换热面积和水利直径的条件下,Y形构形微通道具有比直微通道更高的传热效率和更低的进出口压降,可以满足高热负荷需求.在研究范围内通道级数为3级、分叉角度为60°、分支个数为2时,散热器的综合换热性能最好.以上结论为Y形构形微通道的结构设计提供了重要依据.     

旋转方位角对斜肋通道内部流动与换热的影响

数值模拟了静止及旋转状态下径向出流60°斜肋通道的内流动与换热分布,研究了3种旋转方位角:0°,-45°及+45°.分析了哥氏力与肋片对旋转通道内流及换热的耦合作用机理.计算结果表明:旋转方位角对通道的内流与换热影响显著.在0°旋转方位角下,通道前后缘的换热差异最大;在-45°方位角下,通道整体换热性能最差;在+45°方位角下,通道具有最佳的整体换热性能.