超临界压力下航空煤油圆管流动和传热的数值研究

对超临界压力下大庆RP-3航空煤油在小管道内的流动、传热过程进行了数值研究。湍流模拟采用了RNG k-ε两方程模型和Wolfstein一方程模型结合的两层模型;同时,采用煤油的10组分替代模型以及NIST Super-trapp程序库对大庆3号航空煤油的热物理和输运特性进行了确定。圆管传热的计算条件为:入口压力4 MPa,入口温度300 K,质量流量范围:0.06～0.12 kg/s,壁面热流密度范围:300～700 kW/m2。计算结果显示,煤油的流动和传热特性比水、二氧化碳等简单化合物复杂得多。在超临界压力下,煤油的吸热升温导致其热物理特性以及流动特性均发生剧烈变化,其中,雷诺数沿管道方向上升了至少一个量级,而普朗特数下降了一个量级。在加热开始段,煤油的对流传热系数迅速上升;当壁面温度超过其拟临界温度后,对流传热系数略有所回落;随着煤油温度的进一步上升,传热系数又得到明显增强。计算表明,煤油对流换热特性的变化与煤油复杂的高温热物理特性以及湍流流动在近壁区的增强和抑止有关。     

点阵单元类型和排布角度对传热特性的影响

基于瞬态液晶测试技术和数值仿真,在相同孔隙率条件下,研究了桁架型点阵结构的单元类型和排布角度对矩形通道壁面对流传热特性的影响.结果表明:点阵单元类型对壁面传热性能有决定性的影响,其中交叉排布的体心立方(body centered cubic-0,BCC-0)和Kagome-0点阵归因于出色的传热强化效果,雷诺数在5 7...     

超临界压力低温甲烷波纹管内强化换热数值研究

以发动机主动再生冷却系统为研究对象,建立了碳氢燃料热物性高精度计算方法,在此基础上对超临界压力下低温甲烷在水平波纹管内的流动换热现象展开数值研究,初步分析了波纹管强化换热机理。进一步系统研究了波纹管节高比、管壁材料导热系数、壁面热流密度、入口压力以及雷诺数对强化换热和阻力特性的影响,并采用综合换热性能评价准则对各种因素的影响进行评价。研究表明:在超临界压力下合理选择波纹管可以显著提升换热能力,消除传热恶化现象,并且不会带来明显的压降损失;存在一个最优波高和最佳雷诺数,使波纹管具有最优的综合换热性能;增大管壁材料导热系数和甲烷入口压力可提高换热能力。     

具有凹陷涡发生器的网格冷却结构内的流动与传热特性

对一种具有U 形子通道以及凹陷涡发生器的新型网格冷却结构内的流动和传热特性进行了实验和数值计算研究,通过数值计算可以较准确地获得网格冷却结构的传热性能和流阻性能.在数值计算模型得到验证的基础上,还分别对具有U 形子通道和凹陷涡发生器的网格冷却结构,以及具有U 形和矩形子通道的网格冷却结构内的流动和传热性能进行了对比研究.对比研究表明:与矩形子通道相比,U 形子通道结构能够改善子通道底部的流动状况,降低流动阻力,但同时也略微降低了总体传热性能;在U 形子通道壁面上布置凹陷涡发生器能够显著增强底部流动扰动,大幅提高传热性能,但同时也增加了流动阻力.      

飞机发动机冷气道与隔热层的耦合传热分析

数值研究了某飞机发动机外侧冷气道与隔热层的耦合传热过程。采用低雷诺数k-ε模型与SIMPLEC算法计算通道内可压缩变物性气流的湍流对流换热,采用蒙特卡罗法求解通道壁面间的辐射换热。通道内气流湍流对流换热、壁面间辐射换热与隔热层内导热耦合求解。通过模拟计算,分析了通道与隔热层的耦合传热机制,考察了相关参数的影响。结果表明,在所考虑的通道结构与空气流条件下,冷气道外环壁面的温度高于气流温度,气流对内外环壁面均起冷却作用;在隔热层参数不变条件下,壁面间的辐射换热与气流的对流冷却是该传热过程的控制机制,增大冷气流量、降低壁面发射率均可显著降低隔热层的外壁面温度。     

超临界压力下正癸烷在水平矩形冷却通道内的流动传热数值模拟

为了研究浮升力效应对超燃冲压发动机燃烧室不同位置冷却通道传热特性的影响,对超临界压力下正癸烷在水平矩形冷却通道内的流动传热-固体导热耦合过程展开了详细的数值模拟研究。重点考察了燃烧室不同位置的冷却通道中浮升力对温度以及热流分配的影响及其机理。结果表明：由浮升力引起的二次流动使燃烧室不同位置的冷却通道温度分布和热流分配呈现显著的差异;浮升力提升了燃烧室不同位置冷却通道的换热效果,其中受热方向和重力作用方向相同的冷却通道换热性能提升得最多;修正的Jackson&Hall经验公式不能预测浮升力对冷却通道壁面对流换热的影响,需要寻找其他的经验公式或使用CFD手段进行计算分析来解决这一问题。     

滑油射流冲击旋转壁面的流动换热数值模拟研究

为研究射流流速与发动机转速对轴承腔旋转壁面冷却效果的影响,采用数值模拟方法研究了滑油射流冲击旋转壁面的流动与换热特性,重点关注和分析了中高转速下射流雷诺数Rej和旋转雷诺数Reω对于滑油流动和壁面整体及局部传热的影响。为体现不同工况对壁面换热的影响,分别采用了壁面平均传热量、壁面努赛尔数、局部平均努塞尔数进行了综合换热分析。研究结果表明,随着射流流速增加,壁面附着油量增多,壁面平均传热量与局部换热分别增加约32%和40%;而随着壁面转速增加,壁面附着油量下降,壁面平均传热量减少约13%,局部换热增加约62%。     

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

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

肋开孔高度对大宽高比矩形通道流动传热的影响

为了获得大宽高比矩形通道内开孔肋的流动与传热特性,采用数值模拟的方法在开孔肋的肋高e/Dh=0.188、开孔率β=0.131时研究了通道的摩擦系数、努塞尔数等参数随孔排高度和雷诺数的变化规律。结果表明:与实体肋相比,开孔肋能够减小通道的摩擦系数并提高其壁面温度分布的均匀性,但是传热增强因子有所减小;随着孔排高度的提升,通道的摩擦系数单调减小,壁面温度分布的不均匀度增大,而传热增强因子则先增大后减小,因此存在一最优孔排高度(h/e)opt=0.65使开孔肋的强化传热综合指标达到最大值;随着雷诺数的增大,开孔肋通道的摩擦系数缓慢减小、换热逐渐增强,这与实体肋的变化规律是类似的。     

分形树状通道换热器内的流动换热特性

建立了分形树状通道换热器中层流流动与传热的三维稳态模型,采用流固耦合计算方法对入口水力直径为4mm的矩形截面树状通道内流动换热进行了数值模拟,重点研究了分叉效应对传热的强化机理和换热器受热面的温度分布。研究结果表明:分叉处形成的二次流能有效地强化换热;与传统的蛇形通道相比,分形树状通道换热器具有温度均匀性好、压降小的明显优势。在相同入口雷诺数时,分形树状通道换热器受热面的最大温差远小于蛇形通道换热器,另外,分形树状通道的层流流动压降较之蛇形流道可减小50%以上。同时,加工了分形树状通道换热器及蛇形通道换热器各一套,对数值模拟结果进行了实验验证。实验值与模拟值能较好地吻合,证明了所建流动换热三维数值模型正确可信。     

梯形和矩形通道内短扰流柱排流动与换热计算

为了得到逐渐收缩的梯形通道内扰流柱排的流动换热的规律,对梯形通道内扰流柱排的端壁换热和压力损失进行了数值计算,并与矩形通道进行比较。计算结果表明:(1)梯形通道与矩形通道的端壁总平均换热系数相差不大,但是梯形通道内每排扰流柱的Nu数相差较大。(2)相同来流Re数条件下,梯形通道的压力损失系数远大于矩形通道。所以,在实际计算逐渐收缩的梯形通道内扰流柱排的平均换热时,可近似采用矩形通道内扰流柱的实验关联式,并且将每排扰流柱分别计算。在计算压力损失时,不能将梯形通道近似成矩形通道。     

涡轮叶片尾缘出流对带叉排扰流柱阵列的内冷通道传热特性的影响

本文用实验的方法模拟研究了涡轮叶片尾缘出流对带叉排扰流柱阵列的内冷通道的对流换热特性的影响。实验模型通道为矩形通道,扰流柱的排列参数S_n/d为6.5,S_p/S_n为1.2,其中d为扰流柱的直径,S_n为扰流柱横向间距,S_p为2倍的纵向间距。矩形模型通道的长边与转轴的夹角为20°。研究了出流孔总面积A₀与通道最大截面积A的比值,雷诺数Re,转动数R₀等参数对传热特性的影响。并比较了尾缘出流与叶尖出流的换热特性。      

轴承腔内壁与油膜换热的数值模拟与试验

航空发动机后轴承腔内壁与滑油的换热分析是轴承腔热防护结构设计的基础。对航空发动机轴承腔内壁换热模拟试验件开展试验与数值模拟研究,得到了滑油油膜对轴承腔内壁的换热影响。通过测量试验件外壁面、内壁面以及滑油油膜的温度得到了试验件内壁面换热热流密度与换热系数的分布;结合CLSVOF(Coupled Level Set and Volume Of Fluid)油/气两相流以及热-流-固耦合计算方法对试验件进行了换热分析,并将内壁对滑油的局部热流量的计算结果与试验结果进行了对比,结果显示两者在各个工况下均吻合较好。通过将局部换热系数计算值与当地的滑油流动雷诺数Rel进行对比分析,结果显示内壁局部努赛尔数Nuw与Rel的0.7次方成正比关系。另外,对转速对换热的影响进行分析得到Nuw与旋转雷诺数Re_(rot)的0.345次方呈正比关系。     

涡轮叶片内部冷却通道传热和压力 分布特性的实验

针对带45°肋片圆形弯头U型通道进行了传热和压力分布特性的实验研究.实验中进口雷诺数变化范围在30000~55000之间.研究表明,肋片引起的横向二次流在弯头前区域是强化传热的主要因素.不同雷诺数下传热分布和压力分布的趋势完全一致,但弯头区后半段上传热分布有所差异.通道中的最大压降出现在弯头区,进口雷诺数越大,弯头损失系数和总损失系数都有所减小.流体在圆形弯头中,加速效果没有矩形弯头强烈,弯头区和弯头后的传热较矩形弯头都有较大幅度地减小.      

带前缘冲击复杂叶片内冷通道换热特性实验

采用瞬态热色液晶方法对不同雷诺数下静止状态涡轮工作叶片带前缘冲击复杂蛇形通道的换热特性进行了考察,得到了不同位置吸、压力面的表面传热系数分布.实验结果表明:相同雷诺数下,叶片通道内表面带肋区域两肋间的局部换热强于光滑区域,相邻两肋间的传热系数沿流动方向呈单峰分布.随着雷诺数增加,内表面传热系数提高,在弯头下游的光滑区域离心作用加强,带来该处局部换热分布不均匀,而弯头下游肋的出现可以削弱这种效果.此外,通道的收缩扩张导致主流速度的变化也给局部换热的数值和分布带来了影响.      

带射流的收缩型通道内部换热特性液晶瞬态实验

将航空发动机进气道支板冲击腔简化和放大为一带射流的收缩型通道,并采用最新窄带热色液晶全表面瞬态测温技术对其内表面进行冲击换热实验,具体研究了射流雷诺数、孔径及孔间距变化对努赛尔数分布及大小的影响。实验结果表明:射流雷诺数的增大、孔径的增大及孔间距的减少均使通道内部换热得以加强,但努赛尔数分布的变化及平均努赛尔数的增幅不尽相同;并且侧壁换热受孔间距影响最大,前缘换热则受射流雷诺数的影响最大。     

用高效矢通量分裂格式求解二维叶栅的跨音流动

提出了一种高效隐式矢通量分裂格式。该格式通过在显格式上添加合适的高阶小量, 避免了近似因子分解和系数矩阵运算, 大大地减少了计算量。同时该格式在收敛性及激波捕捉精度等方面均保持了矢通量分裂格式所具有的优势。将该格式用到了二维叶栅的无粘跨音流场的数值模拟中去, 结果令人满意, 证实了无矩阵运算、无近似因子分解的隐格式构造有一定的推广意义。      

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

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

锥柱形药柱环形槽内传热现象的数值研究

采用标准k-ε湍流模型和离散传递热辐射模型对不同来流雷诺数下锥柱形药柱前向环形槽内的流场与传热进行了数值研究。对槽内各表面与热气流之间的换热量进行了比较,并对努塞尔数沿环形槽各表面的变化进行了分析。结果表明来流雷诺数的大小影响环形槽内的旋涡数和旋涡区域的分布,进而影响槽内温度的分布,最终影响槽内装药表面与点火热气流之间的换热量。整个计算结果可为装药设计和分析推进剂的热解燃烧提供参考。     

有预旋进气转静盘腔中的流动和换热特性数值研究

运用RNG湍流模型对有静盘外缘预旋进气和轴向中心进气、转盘外缘轴向出气和外围屏径向出气的转-静盘腔中流动和传热过程进行了数值研究。在盘腔结构中,静止内隔片和旋转内隔片将盘腔分为内转-静盘腔室和外预旋腔室。由于盘腔流动结构具有周向周期性,取盘腔的1/30作为计算域。研究发现:旋转雷诺数Re_θ、预旋喷嘴进气无量纲流动速度C_w,_p、静盘中心轴向进气的无量纲流动速度C_w,_d和预旋比β_p都不同程度地影响盘腔内流动和换热。预旋和隔片的协同作用可以很好地起到阻止热燃气入侵转-静盘腔和提供低静温冷却气体流入涡轮叶片冷却通道。