液体火箭发动机推力室冷却通道温度分层数值研究

为了研究冷却剂温度分层的形成机理及其对流动和换热的影响, 应用雷诺应力模型(RSM)对液体火箭发动机推力室再生冷却通道的流动与传热进行了三维数值模拟, 冷却剂为气氢, 考虑其物性随温度和压力的变化.所得结果表明:冷却剂在非流动方向会出现温度分层现象, 随着冷却剂的不断受热, 温度分层现象越明显, 由于喉部二次流加强了冷却剂间的混合, 在喉部区域温度分层被减弱, 温度分层对冷却剂温升及压降影响较小, 严重影响气壁温度的估算.      

通道参数对再生冷却通道流动换热的影响

应用RSM模型对冷却通道的流动与换热进行了三维数值模拟,冷却剂为气氢,考虑其物性随温度和压力的变化.所得结果表明:增加壁面粗糙度使冷却剂换热强化,但会增加流阻损失;在突扩突缩区域会出现旋涡,旋涡使局部流阻损失加大且使湍流加强,壁温在旋涡出现处降低;冷却通道内的流动发展不受入口湍流强度的影响;冷却剂离心力引发径向平面内的二次流动,二次流引起的冷却剂质量重新分布使传热在凹曲率段强化,凸曲率段恶化.      

考虑煤油裂解效应的超声速燃烧室再生冷却过程分析

为了分析宽马赫数飞行条件下超声速燃烧室再生冷却性能以及考虑燃料高温裂解效应对冷却的影响,发展了具有一定通用性的超声速燃烧室再生冷却系统气-固-液传热分析模型,对燃烧室内流、冷却剂流动以及冷却结构进行了气-固-液传热耦合计算.燃烧室内流计算模型无需实验测量的静压数据以及总温/释热分布假设,通过直接求解质量、动量、能量守恒微分方程并结合燃料混合及燃烧模型来获得内流参数分布.同时对燃烧室壁面传热进行了计算,将冷却结构内冷却剂的流动、换热与燃烧室内流耦合,并且着重考虑了煤油作为冷却剂,其物态随温度、压力变化以及高温时出现的热/催化裂解吸热化学反应.基于实验数据发展了煤油热/催化裂解总包反应模型,对煤油热裂解和催化裂解两种过程的化学吸热性能进行了对比,研究了热/催化裂解效应对再生冷却的影响.     

层板推力室再生冷却通道的传热特性分析

提出液体火箭发动机层板推力室再生冷却通道传热过程的数理模型。采用通用形式控制方程处理冷却剂紊流换热和通道材料导热的共轭传热问题，计算采用LVEL紊流模型，并考虑冷却剂（氢）的热物性参数随温度和压力的变化及层板材料热物性随温度的变化。结果表明，采用大高宽比、小气壁厚度的通道设计，可显著提高再生冷却能力，降低室壁温度和温差。采用对流换热系数和热物性为常数的简化处理会引起很大误差。     

再生冷却推力室耦合传热数值模拟

对大推力液体火箭发动机再生冷却推力室内部燃气、室壁和再生冷却剂进行了耦合传热数值计算.采用二维轴对称N-S方程描述推力室内部燃气的湍流流动与传热,对冷却剂流动采用简化的一维模型,通过室壁的偶合传热采用一维传热模型.N-S方程的求解采用贴体坐标系下的有限容积法,速度和压强的耦合采用可压缩的SIMPLE算法.湍流的模拟采用可压缩的标准κ-ε模型,辐射传热采用热流模型计算.研究表明,本文方法可较好地模拟燃气的二维流动,同时能快速计算壁面热流密度、壁温和冷却套温升,计算结果对大推力发动机推力室的设计具有一定的指导意义.     

燃气非平衡流再生冷却流动传热数值模拟

为准确预测液体火箭发动机推力室身部再生冷却换热状况,采用数值模拟方法,对燃气、推力室壁和超临界气氢进行三维流动和换热耦合计算。采用6组分9步反应动力学模型计算燃气的非平衡化学反应,采用DO模型计算燃气辐射换热,考虑超临界气氢物性随温度和压力的变化。获得了室壁温度场、燃气及冷却剂流场。结果表明,Redlich-Kwong方程、Peng-Robinson方程、Lucas法、TRAPP法能分别准确计算超临界氢的密度、定压比热容、粘度、导热系数,采用燃气非平衡流计算所得流场值更符合实际情况。     

通道深宽比对液体火箭发动机推力室再生冷却的影响

应用湍流模型对液体推进剂火箭发动机再生冷却推力室通道的流动与传热进行了三维数值模拟, 冷却工质为氢气, 其密度、导热系数、动力粘度随着温度和压力而变化, 冷却剂比热容及金属固体物性随着温度而变化.计算采用标准k-ε两方程湍流模型及气-固耦合算法.保持再生冷却通道个数及冷却工质进口流量不变, 通过改变通道肋壁厚度来改变冷却通道深宽比, 研究不同深宽比对推力室壁面再生冷却效果的影响规律.计算结果表明:增加通道深宽比对推力室壁面能够起到强化传热的作用, 但同时也增加了冷却通道的进出口压差.这是由于冷却工质流速的增高, 从而提高了推力室传热系数.随着深宽比不断增加, 推力室再生冷却效果趋于饱和, 而冷却工质进出口压降则不断上升.      

液体火箭发动机推力室发汗冷却传热过程的数值模拟(Ⅱ)数值方法与计算结果

对液体火箭发动机推力室发汗冷却传热过程的二维局部非热平衡模型进行了数值计算。计算中采用了正交曲线坐标系（贴体坐标）,并计及了冷却剂（氢）的热物性参数随温度和压力的剧烈变化及固体壁沿轴向的导热。结果表明：推力室多孔壁面中靠近燃烧室的部分温度梯度很大;固体骨架与冷却剂的温度差异在推力室内壁面上最大;推力室多孔壁面材料导热系数的提高有利于降低壁面温度及温度梯度;随着冷却剂流量的增大,推力室壁中的最高温度明显下降;若设计合理,发汗冷却所需要的冷却剂的量只占总流量的２％左右。     

液体火箭发动机推力室复合冷却流动与传热研究

为了预测液体火箭发动机推力室的复合冷却性能,建立了推力室再生冷却通道和超临界氢的三维仿真模型以及推力室内燃气和超临界氢膜的轴对称二维仿真模型。通过边界耦合发展了液体火箭发动机推力室复合冷却流动与传热的数值仿真方法。对航天飞机主发动机推力室内部燃气、超临界冷却膜、室壁和再生冷却剂进行了流动与传热耦合计算仿真研究。研究表明,仿真方法可较好地预测推力室燃气及再生冷却剂的流动和传热,计算得到航天飞机主发动机的燃气侧壁面最高热流密度为129MW/m2,最高壁温为885K,冷却剂温升为192K,压降为8.8MPa,结果与已有数据吻合较好。模型和仿真方法可用于液体火箭发动机推力室冷却系统传热计算和冷却结构的优化设计。     

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

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

液体推进剂火箭发动机推力室再生冷却通道三维流动与传热数值计算

应用湍流模型对液体推进剂火箭发动机再生冷却推力室通道的流动与传热进行了三维数值模拟,冷却工质为氢气,其密度、导热系数、动力粘度随着温度和压力而变化,冷却剂比热容及金属固体物性随着温度而变化。计算采用标准k-ε双方程湍流模型及气-固耦合算法。结果表明:推力室燃气侧壁面的温度和热流密度的最高点均发生在喉部附近,喉部横截面固体区域最大温度梯度靠近燃气,喉部附近氢气在垂直主流方向的截面上产生了二次流。气固耦合面最大热流密度及最大对流换热系数同样位于推力室喉部附近。      

液氧甲烷发动机台阶型冷却通道的耦合传热特性

为了研究液氧甲烷发动机再生冷却通道中跨临界甲烷的流动和传热特性,以及冷却通道较大幅度的突扩突缩对冷却效果的影响,采用整场直接耦合的方法对推力室三维耦合传热进行了数值模拟,考虑了燃气的非平衡流动.通过计算得到了推力室三维温度场和流场.计算结果表明：由于喉部截面附近存在较强的二次流,燃气侧壁面温度的最大值出现在喉部上游.由于突扩突缩处存在较强的旋涡运动,冷却剂的湍流强度增强,冷却剂侧表面传热系数显著提高,燃气侧壁面温度出现局部极小值,同时也产生了较大的局部损失.由于铜内衬热阻比镍外套热阻小得多,从燃烧室进入的大部分热量在冷却通道底面和侧面被冷却剂吸收.冷却通道底面的温度和热流密度沿程变化比顶面更剧烈.     

火箭发动机燃烧室液膜-再生复合冷却数值仿真

对液体火箭发动机燃烧室液膜-再生复合冷却进行了数值计算,针对液膜-燃气流场区多组分、轴对称Navier-Stokes(N-S)方程和再生冷却区单组分N-S方程进行求解,并使用k-ε方程求解湍流流动.对文献中的某液氧/煤油火箭发动机燃烧室进行了数值模拟,该模型的计算结果能够与文献中的计算结果较好地吻合.计算结果表明:①液膜-再生复合冷却能有效地减少壁面热流密度和降低壁面温度,且其形成的冷气边区覆盖了整个燃烧室及喷管壁面;②再生冷却液入口质量流量越大,复合冷却作用越明显,壁面温度越低;③随再生冷却液质量流量的不同其温升在450~600K之间,且质量流量越大,再生冷却液的温升越小.④壁面煤油的质量分数不断下降,在喷管出口壁面处达到最低值,但含有煤油的区域不断变大.      

内冷通道横流条件下气膜冷却特性

为了研究内冷通道横流条件下气膜冷却的流动和换热特性,采用窄带瞬态液晶测量技术获得了内冷通道横流条件下吹风比分别为0.5,1,2时气膜孔下游冷却效率和表面传热系数云图,并通过数值模拟得到了气膜孔内及下游区域流场的详细信息.结果表明:内冷通道横流对气膜孔下游冷却效率和表面传热系数分布有重要的影响.横流增强了气膜孔射流的展向分布能力,增强了高吹风比时气膜冷却效果.另外,气膜孔下游涡的分布出现明显的不对称性,涡的结构更加复杂.      

冷却剂不同流动方式对膨胀循环推力室再生冷却换热的影响

为了解液体火箭发动机膨胀循环推力室再生冷却换热特性,对某一参考发动机推力室和另外两种面积比的膨胀循环推力室建立三维计算模型,采用数值模拟的方法,考察冷却剂的温升、冷却通道压降以及推力室内壁面温度和热流密度的分布情况.重点比较了不同燃烧室圆柱段长度、冷却剂不同流动方式以及不同面积比对以上结果的影响.计算过程中采用二阶迎风格式离散控制方程.计算结果表明:采用逆流冷却时,通过加长推力室圆柱段长度使推力室受热面积增加70%后,冷却剂温升提高了一倍左右;对膨胀循环推力室进行再生冷却时,采用顺流冷却要比逆流冷却的冷却通道压降低,但同时冷却剂温升也较低,并且对喉部壁面的冷却效果较差.     

摆盘流动与换热性能数值仿真

针对受燃烧室出口高温燃气冲刷的摆盘装置的冷却水回路,建立二维、三维物理仿真模型,模拟了冷却水在内部流道的 流动与换热过程,考察了冷却水进口压力、冷却水流量对流阻和换热性能的影响.结果表明:①摆盘冷却水进口压力由2.3×105Pa提高到8.3×105Pa,摆盘壁温变化微小,结构1壁温升高5K,结构2壁温升高7K;②冷却水进口速度由0.5m/s提高到5.3m/s,结构1壁温降低约120K,结构2壁温降低约100K,冷却效果明显;③结构2通过缩小流道的流通面积,能在更小的冷却水流量的工况下得到更好的换热效果.