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

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

液体火箭发动机推力室冷却通道流动与传热数值研究

采用气固耦合算法对液体火箭发动机推力室再生冷却通道的流动与传热过程进行了三维湍流流动与传热数值模拟，冷却工质为氢气，其密度、导热系数、动力粘度随着温度和压力而变化。应用大涡模拟及标准k-ε双方程模型两种湍流模型分别进行数值模拟，详细揭示了再生冷却通道固体区和流体区内的速度场和温度场，并在不同的计算网格数目下对两种湍流模型的计算结果进行了对比。结果表明，在相同的网格条件下，标准k-ε双方程模型与实验数据的吻合精度比大涡模拟模型更好，且满足工程计算精度。随着网格数的增加，大涡模拟的计算精度逐渐得到改善。     

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

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

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

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

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

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

液体火箭发动机推力室冷却通道传热优化计算

采用标准K-ε两方程湍流模型对液体火箭发动机推力室再生冷却通道三维湍流流动与传热过程进行了数值预测，冷却工质为氢气，其密度、导热系数、动力粘度随着温度和压力而变化，通过两种优化方案来改变推力室冷却通道的深宽比。方案一为保持冷却通道的深度及肋宽不变，通过改变推力室壁面通道个数来改变通道的深宽比，方案二为保持通道数目不变，通过增加或降低通道高度来改变通道的深宽比。以此计算在不同通道深宽比下推力室壁面的传热特性，并进行了优化分析。计算结果表明：存在着一个最佳冷却通道个数，使得推力室壁面再生冷却效果达到最佳；在相同质量流量下，降低通道高度能够强化推力室传热，但同时增加了进出口压差。     

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

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

基于L-M算法的液体火箭推力室壁面热流密度测量方法

为获取液体火箭推力室热防护方案构建的基础数据,结合单点瞬态法,提出了基于L-M算法测量壁面热流密度的方法。通过模拟热流密度分析算法中关键参数对求解精度的影响,给出关键参数的选取建议。在此基础上,对过氧化氢煤油推力室典型试验工况进行热流密度测量。结果表明：圆筒段末端壁面热流密度为3.35MW/m²,喉部位置单组元工况下壁面热流密度为2.21MW/m²,双组元工况下壁面热流密度为10.59MW/m²,热流密度变化过程与室压变化过程完全对应。变工况过程中近壁处燃气温度存在恢复过程,燃气温度稳定后,壁面热流密度也达到稳定状态。本文提出的方法可迅速、准确地测量推力室壁面热流密度。     

液体火箭发动机推力室发汗冷却传热过程的数值模拟：（II）数值方法与计算结

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

推力室喉部层板发汗冷却段的结构设计分析

根据层板单元的热分析结果和层板发汗冷却推力室固有的结构特点,提出这种先进发动机冷却方案的设计原理和结构参数的计算公式,结合一台50kN发动机推力室喉部改再生冷却为发汗冷却的改形设计,分析其层板发汗冷却段喉部的设计方法和主要结构尺寸的计算结果。还比较了全再生冷却和发汗冷却两种冷却方式下发动机推力室的温度、热流密度和重量。对先进层板发汗冷却推力室的结构设计提供了参考。     

气壁镀镍和冷却剂入口对再生冷却的影响

考虑冷却剂入口模型及气壁镀镍,对液体火箭发动机推力室再生冷却通道和冷却剂进行三维流动与传热耦合计算.采用经验公式计算燃气侧对流及辐射换热,冷却剂为甲烷,考虑其物性随温度和压力的变化.所得结果表明:冷却剂入口二次流及突扩流场的叠加,使局部压力损失变大,影响进口下游较大区域流动传热状况;气壁镀镍能提高气壁的耐高温性,降低气壁锆铜的温度.      

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

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

双向旋转变喉面流量调节阀热应力分析

以ANSYS软件为工具,时均化N-S方程、RNGkε-湍流模型和傅里叶热传导模型为基础,建立三维面对称纯气相流动和固体热传导模型;采用流固耦合的方法,设置流场和固体区域交界面为热耦合面,分析了喷管内流场特性和喷管体与轴阀固体区域的温度、应力分布。结果表明在喷管喉部圆柱段的温度最高,调节阀的等效喉部段和等效扩张段出现了应力集中的现象,两者都随轴阀转角的增加、工作时间的延长而增大。     

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

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

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

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

高超声速飞行器热环境与结构传热的多场耦合数值研究

为了准确预测高超声速飞行器面临的严峻气动热/力环境以及结构的热力响应,发展了高超声速流动与结构传热耦合框架。采用分区求解方法,通过耦合界面的实时数据传递,实现了基于Navier-Stokes方程的高超声速化学非平衡计算流体力学（CFD）求解器与结构的热力全耦合有限元法（FEM）求解器的多场耦合计算,建立了高超声速飞行器的多场耦合数值分析方法。首先对经典高超声速圆柱绕流实验进行了耦合计算,结果与实验值吻合良好。然后针对典型的超高温陶瓷（UHTC）材料的耦合传热问题进行了数值研究,考虑热传导效应对气动热环境和结构热响应预测的影响,结果表明对于复杂外形且热导率相对较高的UHTC材料,结构内部热传导对热环境和表面温度分布的影响不可忽略。最后针对UHTC材料热物性（比热和热导率）非线性对高超声速流动传热过程的影响进行了研究,结果表明当比热和热导率处于合理的误差范围内时,材料表面温度响应对其变化并不敏感。     

双脉冲固体发动机喷管传热烧蚀特性

为了研究双脉冲固体发动机喷管的传热烧蚀特性,由燃烧室压强及发动机推力试验曲线得到了喷管喉径的瞬变值,由FLUENT流体计算软件进行流固耦合传热烧蚀计算,得到了喷管瞬态温度分布、绝热材料热解炭化情况及碳/碳(C/C)喉衬瞬态烧蚀率,分析了脉冲工作过程及脉冲间隔时间对喷管传热烧蚀的影响.计算结果表明,脉冲工作过程中,绝热材料热解线、炭化线向材料内部扩展,喉衬烧蚀率不断增大;脉冲间隔时间内,喷管材料内部的导热使各处温差减小,温度趋于一致;第一脉冲的传热烧蚀与脉冲间隔的材料导热使第二脉冲工作时喉衬整体热沉小、内壁初始温度高、表面粗糙度大,从而导致较高烧蚀率.