主动冷却通道内正庚烷热裂解与换热过程的数值模拟 *

为了深入研究碳氢燃料在主动冷却通道中的传热特性,采用详细基元反应机理对微小通道内正庚烷热解反应耦合流动换热过程进行了数值模拟。根据燃料的热裂解反应和壁面传热特性可将通道分为三个换热区域。通过对不同的入口温度、运行压力、壁面温度和通道直径等典型工况进行模拟,探索了正庚烷的热裂解反应和流动换热过程耦合的基本规律。研究结果表明,入口温度对通道整体的换热性能影响不大;保持其它参数不变,通道直径从2mm减小到0.5mm,其整体换热量降低约10%,单位面积的换热量提升约4倍。     

流场不均匀对液体火箭发动机纵向振荡燃烧影响及控制

针对常温推进剂发动机推力室再生冷却和撞击式喷注器结构,分析了推力室身部与喷注器对接部位的流场特性,对流场均匀性进行了实验测量。结果表明：推力室身部再生冷却通道出口压力存在约0.15 MPa周向不均匀。身部出口节流显著提高局部流速,使喷注器面氧化剂湍流度和不均匀性增加,进而改变燃烧特性。通过撞击喷注单元雾化试验,获得了18 m/s的推进剂入口边界流速。基于喷注器流场均匀性,提出控制推进剂流速,降低不均匀性,进而抑制纵向高频燃烧不稳定性的控制方法。发动机热试结果表明,采用（15±1） m/s的推进剂入口流速,控制方法抑制了纵向高频燃烧不稳定性。     

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

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

钝头双锥喷流致冷流场结构及密度脉动特性

采用大涡模拟方法对钝头双锥喷流致冷流场开展了数值模拟,研究了超声速喷流混合流场结构特征及密度脉动特性。大涡模拟方法基于隐式亚格子模型,空间离散采用高精度通量限制型紧致格式,时间推进采用显式Runger-Kutta方法。数值模拟清晰地捕捉到了流场波系结构,精细地预测了流动发生失稳、转捩以及发展为充分发展湍流的物理过程,直接获得了流场密度脉动特性。通过有、无喷流状态对称面流场的对比,发现超声速喷流能够有效冷却光学窗口;喷流与主流形成的混合层不稳定,很快发生失稳和转捩,形成大尺度湍流结构,进而引起强烈的密度脉动。此外,获得了钝头双锥整体模型喷流致冷流场的空间发展形态特征。      

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

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

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

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

轴对称内锥流动中马赫盘的形成与演化

针对轴对称内收缩直锥流场中的马赫反射问题，采用无黏数值模拟结合理论分析，在来流马赫数6及直锥入口半径R恒定的条件下，以模型壁面前缘角度θ_w和壁面长度w/R为几何参数，研究了内锥形激波汇聚增强过程以及马赫盘下游流动对马赫盘位置的影响。结果表明，壁面前缘角度θw通过改变轴对称内锥形入射激波汇聚增强过程中达到von Neumann强度的位置，限定了马赫盘可能存在的位置范围；马赫盘位置还受其下游流动的影响，在不同的θ_w和w/R条件下，马赫盘下游流管内声速喉道的形成机制可以划分为依赖或不依赖壁面尾缘膨胀波两种类型；对于不依赖于壁面尾缘膨胀波的情况，声速喉道形成之前流管内的压力仅与反射激波的波后非均匀压力相匹配，马赫盘的位置不随w/R而变化；对于依赖于壁面尾缘膨胀波的情况，声速喉道形成之前流管还需要匹配壁面尾缘膨胀波引起的压力变化；根据反射激波下游与马赫盘下游流管的压力匹配关系，可以求解马赫盘位置。     

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

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

圆盘结构下旋转爆震波的二维数值研究

为揭示圆盘形燃烧室内旋转爆震波的传播特性,以2H_2+O_2+3.76N_2为反应混合物,开展该结构下的二维数值研究。实现了旋转爆震波的成功起爆并得到了两种稳定的流场结构;详细分析了两种流场的特点,以及不同直径处爆震波参数及出口马赫数的变化规律。结果表明,旋转爆震波在壁面形成的反射激波对流场中的温度、压力分布以及流线都产生影响;随着直径的不断减小,旋转爆震波压力、温度以及传播速度都不断减小;流场稳定后,出口马赫数随时间的变化趋势与压力趋势相符,呈现"双波峰"特征。喷注总压为0.4MPa时,喷注入口的阻塞比为21%,超声速流动占整个出口流动的52.17%,出口流动的平均马赫数为0.98;当喷注总压降为0.2MPa时,入口阻塞比为20%,超声速流动占21.99%,出口流动的平均马赫数为0.95。     

超临界二氧化碳超声速流动中的激波问题研究

对于理想气体而言，激波前后的流动参数变化受马赫数和气体比热比的影响，但对于超临界二氧化碳(S-CO₂)，其物性规律在临界点附近变化显著，使得激波后的流动参数难以预估。本文通过推导激波关系式，并结合S-CO₂物性数据库，建立激波计算的数值迭代程序，对S-CO₂超声速流动中的正激波和斜激波进行研究，并将理论分析结果和CFD仿真结果进行对比。结果表明：激波压缩不会造成气态CO₂向液态CO₂转变，但会使得气态或液态CO₂向超临界态转变；CO₂正激波前后压比随马赫数变化主要受来流相态影响，温度比变化则主要受来流初始参数影响；S-CO₂正激波波后参数在不同来流条件下的变化趋势具有相似性，当来流压力在临界压力附近时，波后流动参数变化剧烈；马赫数越高，激波出现最小熵增时对应的来流温度-压力曲线越靠近准临界线；在准临界线两侧，S-CO₂的斜激波脱体角度随来流压力增大分别呈现出“类气体”和“类液体”的变化特性。     

方形通道内超临界碳氢燃料传热恶化数值研究

为了深入理解火箭发动机再生冷却过程中碳氢燃料的流动传热规律,采用RNG k-ε湍流模型结合增强壁面处理方法对非对称受热(上壁面外侧加热)方形冷却通道内超临界压力甲烷的对流换热进行了数值研究。重点考察了加热壁面内侧的传热恶化以及由传热恶化和固壁热传导共同作用引起的热流异常传递现象,拟合得到了传热恶化的临界热流密度和起始内壁温关系式、加热侧内壁面和侧壁面内侧平均热流密度的预测关系式。结果表明:当外壁热流密度和质量流速比值高于0.288 k J/kg时,近壁流体流动加速诱发了加热侧内壁面的传热恶化;同时,固体区域温度畸变导致加热侧内壁面热流密度减小,热流更多地向侧壁面内侧传导。运行压力越低,该现象越显著。     

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

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

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

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

高压涡轮导叶内冷通道流动特性计算分析

涡轮叶片内部冷却特性直接影响到叶片的壁面温度分布,必须对其展开详细的研究。本文利用网络法计算研究了某高压涡轮导叶内部冷却通道中的流动和换热特性,并进行了相应的流量特性试验研究。结果表明,改进后的计算程序可以较为准确地模拟出叶片内部冷却通道中的流量分配和压力、温度分布,以及冷气沿径向和弦向的压力损失。流量特性计算结果与试验结果吻合较好。      

引气结构对二元塞式喷管冷却和红外辐射特性的影响

为降低加力状态下二元塞锥表面温度和喷管红外辐射强度，对塞锥进行冷却结构设计。采用数值模拟的方法对比分析了引气结构、冷却通道高度和冷气入口总压比对塞锥冷却和喷管红外辐射特性的影响。结果表明：塞锥冷却后其表面温度和喷管红外辐射强度显著降低；引气腔内无冲击板时，引气角度的改变引起射流核心区位置的变化，造成塞锥头部和前缘展向温度分布差异明显，引气角度为90°时塞锥表面最高温度要比30°和60°的模型高50K；加装冲击板后，冷却通道内的流量分配和塞锥前缘的展向温度分布得到有效改善、塞锥头部的换热得以增强，但同时会引起较大的总压损失，因此相同入口总压比下，加装冲击板后冷却流量降低、塞锥外表面温度升高；随着冷却通道高度增大，冷气流量增加、流速降低，故存在一个最佳通道高度使得塞锥冷却效果最好；以塞锥无冷却为基准，入口总压比为1.0～1.8时，塞锥外表面最高温度降低了470～590K,0°探测角上红外辐射强度降低了25%～33%。     

等离子体炬电弧放电热膨胀不稳定性分析

针对层流电弧中的不稳定性以及电弧脉动的起因,利用简化的一维放电模型,分析了稳态放电非转移弧层流等离子体射流的放电热膨胀不稳定性表现为射流周期性震荡的宏观性质,提出电弧放电室进气流初始定向流速增加能够抑制等离子体炬射流放电不稳定性的发展。放电热膨胀不稳定性发生的主要机制是:由于放电气体密度的随机涨落与电场加热气体功率的反相关变化,导致电弧放电通道内温度高的区域随着加热功率的增加而正反馈增加,由此导致电弧放电热膨胀不稳定;放电气体的定向流速度为射流增加的一项动压与密度变化成正比,导致放电热膨胀不稳定性的发展具有阈值限制;增加定向流速度能够抑制放电热膨胀不稳定性。     

热辐射对气膜冷却叶片冷却性能影响

针对一种带有气膜冷却结构的涡轮一级导向叶片进行气-固-热耦合数值模拟,通过比较考虑/不考虑热辐射的温比和综合冷却效率,分析了多种辐射因素对叶片表面温度和冷却性能的影响。结果表明:入口黑体辐射温度在1200~1900K之间,叶片表面发射率在0.3~0.7之间时,考虑热辐射作用均会使叶片表面温度明显上升。入口黑体辐射温度1600K,叶片表面发射率为0.5时,叶片压力面温度整体上升约100K,叶片表面最高温度点（1350K）温度上升约50K;气体辐射对叶片吸力面和尾缘区域造成5%左右的温升;考虑辐射作用使得叶片综合冷却效率下降,叶片前缘和压力面尽管布置密集的气膜孔仍然难以满足冷却需求,综合冷却效率下降至0.3以下。      

局部抽吸对边界层流动稳定性的影响和控制研究

在边界层壁面上,设计局部抽吸结构,采用直接数值模拟的方法,获得稳定的三维基本流.在此基础上,研究稳定及最不稳定的二维扰动T-S波的时、空演化机制;进一步探讨了局部抽吸的形式、强度大小及分布结构对二维T-S波的非线性演化影响及其对增长率的贡献大小.结果表明,局部抽吸结构诱导产生的三维基本流是扰动波得以快速增长的一个关键性因素,这是由于平均流剖面的改变及展向速度的出现,增强了流体运动中的不稳定性、扩大了中性曲线的不稳定区域范围.在最不稳定的二维扰动T-S波的非线性演化过程中,由于非线性作用的不断增强,逐渐激发产生出三维扰动波及高次谐波,其三维扰动波的流向波数和频率与二维扰动波的流向波数和频率相同;同时展向速度的大小对二维扰动波的增长、流动的失稳、流向涡的形成等方面都起着激励的作用.随着时、空的不断发展和非线性作用的迅速加强,正、负相间的流向涡逐渐形成,强度逐渐增大,流向涡的影响区域也在不断扩大,涡的形状逐渐拉伸变长,并出现强的剪切层,流动开始失稳等其它机制;这些结论与文献[4、5]的结果相吻合.     

压气机旋转不稳定性的周向模态特性及其分解方法

基于实验和数值模拟的方法,研究压气机叶顶区域不稳定压力波的周向模态特性。使用基于空间傅里叶变换的周向模态分解方法对实验和数值计算结果进行了分析,给出了实验和数值计算结果的频谱模态分解结果,得到了压气机连续周向模态特征频谱峰值。对压气机叶顶区域各模态幅值和相位的时变特性做了详细分析,得到了叶顶旋转不稳定性周向模态传播速度和传播方向。结果表明:旋转不稳定性的主模态具有稳定的周向传播结构,而其他模态波的幅值和相位存在较大波动。利用多组有限测点组合方式消除了模态混叠现象影响,将可以分析的模态上限提高数倍,取决于测点数的选择;基于该方法可以用有限测点捕捉高阶模态,大大减少了实验所需的传感器数目。      

竖直圆管内浮升力对超临界碳氢燃料裂解传热传质特性的影响

为了探究再生冷却过程中,浮升力对竖直圆管内超临界碳氢燃料裂解传热传质特性的影响,基于详细裂解反应动力学模型,建立了同时考虑碳氢燃料流动传热和裂解吸热的耦合算法,在此基础上对竖直管道内,浮升力对超临界RP-3的流动、传热和裂解反应的影响展开了数值研究。计算结果表明：与不考虑浮升力的情况相比,在浮升力影响显著的条件下,浮升力增强了向下流动的碳氢燃料壁面处与中心流区域的传热传质过程,燃料温度和裂解率的径向分布更加均匀,燃料吸热能力增强,换热系数上升,同时可以有效地抑制管道壁面上结焦的生成;而对于向上流动的流体,浮升力不利于壁面处与中心流区域的传热传质,导致冷却通道内碳氢燃料温度和裂解率径向分布的不均匀性增强,燃料吸热能力降低,换热系数下降,同时增加了管道壁面上的结焦量;同时,为了更好地理解浮升力的影响,本文还对不同壁面热流密度下向上和向下冷却通道内超临界碳氢燃料的裂解传热特性进行了分析;判别式Bo^*<6.0×10^-7不能准确地预测竖直管道内浮升力对超临界碳氢燃料裂解换热的影响。