基于SST k-ω的超声速气膜冷却湍流模型温度修正

针对现有湍流模型无法准确预测非等温超声速气膜冷却行为的问题,在现有SST k-ω可压缩修正模型基础上,以总温梯度为变量,完成了湍流模型的温度修正,并首先通过非等温可压缩自由剪切流动实验数据初步验证了其修正效果,在此基础上对温度修正模型预测超声速气膜冷却传热的准确性进行了验证。结果表明,基于剪切层总温变化的湍流模型修正效果显著,可准确预测大温度梯度下的自由剪切流动轴向速度分布;修正模型计算得到超声速气膜冷却壁面热流分布与对应的实验结果吻合;当用于剪切层温度大梯度变化的超声速气膜冷却数值模拟时,温度修正后的SST k-ω模型与可压缩修正的k-ω模型、SST k-ω模型相比,具有显著优越性。     

先进旋涡燃烧室多场协同分析

为了对先进旋涡燃烧室(AVC)的流动和传热性能进行综合评价,对不同来流速度、来流温度、壁面温度下AVC燃烧性能进行了数值模拟,对速度场、温度场、压力场进行了多场协同分析。结果表明:随着燃烧室来流速度的增大,速度与速度梯度的协同角α、速度与温度梯度的协同角β,速度与压力梯度的协同角θ均减小,温度梯度与速度梯度的协同角γ及压力梯度与速度梯度的协同角φ增大。随着来流温度的提高,α,θ增大,β,γ,φ减小。随着壁面温度的提高,α,θ减小,β,γ,φ增大。对于AVC湍流流场的传热强化问题,增大来流速度和来流温度,降低壁面温度会强化流动换热;增大来流速度和壁面温度,降低来流温度可以减少传热功耗;增大来流速度和壁面温度,降低来流温度能够提高强化传热的综合性能。     

简化三方程转捩模型在高超声速流动中的应用

在雷诺平均方法的基础上,通过耦合k_(-ω)SST湍流模型和间歇因子转捩模型,引入湍流模型和转捩模型的可压缩修正方法,对高超声速平板、双楔、尖锥三类模型边界层转捩流动开展了数值模拟研究。与实验结果的对比分析表明基于压力梯度表征参数T_w=R_(TΩ/ω)的简化三方程转捩模型,能够准确捕捉高超声速平板边界层流动的转捩起始位置、转捩区域长度以及湍流区壁面热流。而对于双楔、尖锥模型,改进前的简化三方程转捩模型由于受到流动可压缩效应的影响,边界层转捩后湍流区的壁面热流模拟预测结果明显高于实验值。在添加模型可压缩修正方法后,转捩区域长度和湍流区壁面热流模拟结果得到有效改善,与实验值吻合较好。可见,简化三方程转捩模型在添加可压缩修正方法后具备准确模拟预测高超声速边界层流动转捩的潜力。     

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

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

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

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

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.      

用于可压缩自由剪切流动的湍流混合长度

抓住可压缩流动变密度特性,构造出基于有效涡量的三维von Karman混合长度。湍流模型采用仅依赖湍动能k的单方程KDO(Kinetic Dependent Only)模型,引入新构造的混合长度替换旧尺度得到CKDO模型。为了验证其描述可压缩自由剪切湍流的能力,选择无壁面束缚、密度梯度大和可压缩效应强的自由剪切混合层为算例,其对流马赫数Ma_c=0.8。计算结果表明,KDO模型对混合层的速度分布有着良好的控制和模拟,而经可压缩修正后的CKDO模型与原模型及其他可压缩修正模型相比,所计算的速度分布、主雷诺剪切力和混合层厚度与试验结果更加接近,说明了该混合长度对可压缩混合层这种自由剪切湍流有着良好的刻画能力。     

扩缩管充分发展层流的强化传热场协同分析

采用二维轴对称不可压缩层流N-S方程,利用有限体积方法和SIMPLE算法,在雷诺数100～1300范围内,对轴对称周期性扩缩管的充分发展层流流动和传热特性进行了数值模拟,重点研究了不同雷诺数、几何外形与换热强度的关系,并利用场协同理论对计算的结果进行了分析,获取了层流态下扩缩管对流换热强度随雷诺数、几何尺寸变化的规律,对进一步研究扩缩管换热问题、探讨强化传热理论和发展换热技术均有实际意义.     

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

高温下加热/冷却管内参与性介质的辐射与对流耦合换热研究

本文通过数值模拟方法研究了高温下加热/冷却管内层流充分发展流动的参与性介质热起始段的辐射与对流耦合换热。采用离散坐标法求解圆柱坐标系下的非线性积分微分辐射传递方程,采用有限差分法求解辐射与对流的耦合换热能量方程。通过计算,分析了不同壁温及入口温度下,介质导热系数、壁面发射率、介质散射反照率及散射相函数等因素对局部努谢尔数、壁面热流密度及有效加热/冷却长度的影响,并同时采用努谢尔数和壁面热流密度评价了耦合换热中辐射对对流换热的影响。      

Numerical simulation of axial liquid film cooling in rocket combustor

Numerical simulation has been done for liquid film cooling in liquid rocket combustor.Multiple species of axial Navier-Stokes equations have been solved for liquid-film/hot-gas flow field,and k-ε equations have been used for compressible turbulent flow.The results of the model agree well with the results of software FLUENT.The results show that:(1) Liquid film can decrease the wall heat flux and temperature effectively,and the cold border area formed by the film covers the whole combustor and nozzle wall.(2) The turbulent viscosity is higher than the physical viscosity, and its biggest value is in the border area of the convergent area in nozzle.The effect of turbulent flow on the whole simulation field can not be ignored.(3) The mass fraction of kerosene at the film inlet is 1,but it decreases along the nozzle wall and achieves its lowest value at the outlet.However,the mass fraction of kerosene near the wall is the biggest at any axial location.      

高超声速飞行器-进气道一体化热流数值计算

采用CFD(计算流体动力学)技术, 开展了飞行器前体/发动机一体化气动热环境分析.对层流区、转捩区和湍流区分别采用计算模型, 在湍流区利用压缩性修正的SSGZ-Jk-ε湍流模型, 在转捩区引入代数型转捩因子模型描述边界层由层流逐渐过渡为完全湍流的流动过程.计算了前体和内通道的表面热流, 并与实验结果进行了对比.结果表明所采用的计算方法可以较好地预测前体及发动机内通道热流率, 流动状态、几何结构及激波入射对热流值影响较大.      

喷注壁面楔状流层流边界层速度与温度相似解

建立了壁面有喷注的楔状流层流边界层冷却数学模型,求解经相似变换得到的描述无量纲流函数和无量纲温度的常微分方程,获得了楔状流层流边界层无量纲速度和温度的相似解,给出了考虑壁面喷注边界条件的形式简洁、物理意义明确的无量纲流函数拟合式;用Runge-Kutta法求解无量纲速度与温度的常微分方程,获得了壁面有喷注的楔状流层流边界层无量纲速度和温度随楔形角、吹风参数、冷却介质温度的变化规律.通过计算0°,18°和36°楔形角的楔状流层流边界冷却速度与温度分布得到以下结论:楔状流速度边界层和温度边界层厚度都随楔形角的增大而变小,随着吹风参数的增大而增大;当冷却介质温度越低、吹风参数越小、楔形角越大时,靠近壁面处楔状流层流边界层内温度梯度越大.      

高强化活塞振荡冷却的场协同分析

数值模拟了高强化活塞内冷油腔中机油的振荡冷却过程,定义了有效速度的概念,通过数据处理得到了内冷油腔壁面的各场参数在时间和空间上的分布规律,分析了壁面表面传热系数与场协同性能间的联系.结果表明:在面相位为60°~120°和240°~300°范围内,场协同匹配性能较好;有效速度的分布规律与协同角余弦值的分布规律一致;内冷油腔壁面表面传热系数分布规律与协同角余弦值和有效速度的分布规律一致,表明可以使用场协同机理阐述机油振荡冷却的对流换热规律.      

碳/酚醛复合材料烧蚀性能的实验研究

总结了电弧加热器湍流导管试验装置上对碳／酚醛复合材料的烧蚀试验结果，利用多元线性回归方法。拟合出烧蚀材料的有效烧蚀焓与冷壁热流密度、壁面压力的关联式。在一定热流、压力范围内可以很方便地计算出碳／酚醛复合材料的有效烧蚀焓和质量烧蚀率。     

超临界压力下航空煤油水平管内对流换热特性数值研究

采用RNG(renormalization group)k-ε湍流模型和近壁区的Wolfstein一方程模型对超临界压力下大庆RP-3航空煤油在水平圆管内的流动和换热特性进行了数值研究.超临界压力下,由于航空煤油在拟临界点附近热物性的剧烈变化,浮升力将引起显著的二次流动.二次流动使得水平圆管的下表面湍流强度和对流换热增强,而上表面的湍流强度和对流换热减弱.最后分析了两种水平管内对流换热受浮升力影响判别标准在超临界流体中的适用性.      

超临界压力下低温甲烷的湍流传热数值研究

通过系统的数值模拟计算,在准确确定甲烷的热力学和传输物性变化的情况下,详细分析了两种热流密度下超临界压力对低温甲烷的湍流传热过程的影响,揭示了对流换热Nusselt数的变化规律。计算结果表明:在超临界压力下,热力学和传输物性对湍流传热现象会造成很大的影响,尤其在甲烷的临界区域附近,由于物性的剧烈变化会导致传热过程的恶化现象;在高热流密度情况下(如7MW/m2),增大管内压力有利于提高对流换热强度;现有的常用变物性湍流传热公式不能适用于超临界压力下低温甲烷的对流换热计算。     

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

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

碳氢燃料裂解吸热反应及超临界传热现象数值模型的构建与验证

 基于正癸烷在超临界压力下的热裂解反应机理,建立了考虑正癸烷裂解吸热反应及超临界传热现象的数值计算模型。采用该模型对正癸烷在3.45~11.38 MPa压力条件下的湍流传热过程进行了数值模拟计算研究,得到了详细的温度、速度、裂解转化率、物性及壁面热流密度的变化和分布情况。通过与已有的实验和数值计算结果的对比,对数值模型和计算软件进行了充分的验证。结果表明,本文的数值模型是准确、可靠的。这就为研究碳氢燃料的裂解吸热反应及超临界传热现象提供了一个有效的模拟计算工具。