流-固耦合计算的应用研究

为了解决用常规算法求出的流体与固体对流换热系数的不真实性问题,采用了流固耦合计算方法计算固体壁面的温度.其中采用模型模拟湍流流动,用壁面函数法修正处于流场内部的固壁,通过固体和流体双向耦合换热计算,得出了整个流场、温度场包括(固体部分和流体部分)的分布.     

涡轮机匣壁面换热特性的数值计算

为了研究涡轮机匣的换热特性,采用分段计算技术,各段连接为耦合连接,对某型涡轮机匣模型进行数值模拟,减小了计算复杂结构的难度,详细分析了机匣流动区域的流动特性和壁面的换热特性。结果表明:孔出流对壁面有很强的冲击作用,导致冲击区壁面换热系数很大;腔内气流流动情况复杂,具有复杂的漩涡结构;气流涡作用于壁面,增强壁面局部换热,换热系数增大;计算结果与实验结果吻合很好,定性地反映了机匣的流动和换热特性。     

半封闭肋化通道射流冲击换热特性的数值计算与实验验证

 运用数值计算的方法对半封闭肋化通道单排孔射流冲击冷却进行了三维模拟,并采用红外热像测试技术对冲击靶面的温度分布进行了实验研究,揭示出各种流动参数和几何参数对冲击射流传热特性的影响规律。研究结果表明：数值计算与实验特征基本吻合;冲击间距对肋化部分的换热影响较大,在本文研究范围内,Z-n/d=1时壁面射流区的换热效果最好;由于产生了两个方向的二次流,V型肋化通道在壁面射流区获得了最佳换热效果。     

燃烧室辐射换热一维区域法模型研究

目前用于工程计算燃烧室辐射热流的Lefebvre零维辐射模型,不考虑燃气区域对附近壁面的辐射影响,过于简化,不能反映实际燃烧室火焰筒内的辐射换热过程。本文发展了一种机理严密,精度高的用于燃烧室壁温计算的区域法辐射换热模型。文中采用Lefebvre模型和区域法模型计算了一个圆筒型不冷却火焰筒结构的换热;对比分析了两种方法计算的壁面辐射热流和壁面温度分布情况。结果表明,使用Lefebvre模型计算时,由于将燃气按一维划分为定向辐射,在燃烧室前后壁面辐射热流较小,存在一定不足;而使用区域法模型计算时,辐射热流分布更趋合理。     

运用区域法模型计算燃烧室内一维辐射换热

采用一种机理严密,精度高的区域法辐射换热模型用于燃烧室壁温计算.采用Lefebvre模型和区域法模型计算了一个圆筒型不冷却火焰筒结构的换热;对比分析了两种方法计算的壁面辐射热流密度和壁面温度分布情况.结果表明,使用Lefebvre模型计算时,由于将燃气按一维划分为定向辐射,在燃烧室前后壁面辐射热流密度较小,存在一定不足;而使用区域法模型计算时,辐射热流密度分布更趋合理.      

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

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

冲击射流流动换热超大涡模拟研究

为了准确预测发动机热端部件中广泛采用的冲击射流冷却复杂的流动和换热特性,发展了基于BSL k-ω模型的超大涡模拟(VLES)高精度模拟方法,并对高雷诺数Re=4×104,两种不同射流距离2和6的单孔冲击射流及三孔冲击射流这一经典的流动传热问题进行三维非稳态高精度数值计算。同时,将分离涡方法 (DDES)和k-ωSST,RNG,Transition SST三种RANS方法的数值模拟和开发的超大涡模拟(VLES)方法进行对比。研究表明,VLES方法均能够准确捕捉冲击射流流场的复杂非稳态流动及传热特征,包括自由射流区、壁面射流区小尺度涡系和大尺度湍流结构的演化和破碎,同时冲击壁面的换热系数计算结果与实验值吻合较好。DDES方法未能准确捕捉流场复杂的小尺度湍流结构,壁面换热计算结果与实验值差异较大。RANS方法计算的换热结果与实验数据差异最大,基本未能预测到壁面换热特性。在相同的计算网格和计算方法下,VLES方法计算结果优于DDES方法,DDES方法一般好于RANS方法。这表明新开发的VLES方法能够准确地计算冲击射流相关的流动及换热问题。     

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

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

三维数值模拟再生冷却喷管的换热

为了解液体火箭发动机喷管再生冷却的换热特点，采用数值模拟的方法，对内喷管燃气、壁面和冷却液建立不同的三维控制方程，进行流动和传热的耦合计算。在计算中，假定喷管流动为冻结流动，考虑燃气向壁面的对流换热和辐射换热；采用二阶迎风格式离散控制方程，采用DO模型离散求解辐射换热方程，水蒸气的吸收系数根据Leckner公式计算。计算模型采用缩比热试车发动机，数值计算结果与实验结果吻合较好，较准确地模拟出了喷管的壁面热流密度，得到了喷管燃气和冷却液的流场和温度场，对高压再生冷却喷管的设计具有指导意义。     

表面凸起对机翼热气防冰腔内换热强化的影响

采用数值模拟对比研究了光滑表面和具有表面凸起结构热气防冰腔内湍流流动的换热特性。机翼防冰腔内笛形管具有三排射流孔,射流孔角度有0°±45°组合以及0°±30°组合。为了强化射流冲击光滑表面的流动换热,在防冰腔内表面正对射流孔的射流冲击区,设计了表面凸起结构,用来强化射流对壁面的冲击换热效果并起到引流作用。通过改变射流孔射流角度研究了射流角度对传热特性的影响。计算结果表明:与光滑防冰腔内表面射流冲击换热相比,表面凸起结构可以将均匀发散的壁面射流集中为高速壁面射流,提高壁面射流区的对流换热系数,从而增强射流冲击换热效果,机翼前缘的强化换热效果尤为明显。     

非对称交叉激波和湍流边界层相互作用的数值研究

针对7°×11°双尖鳍外形的非对称交叉激波与湍流边界层相互作用,采用Navier-Stokes方程和5种湍流模型进行了计算.主要考察对壁面压强、热传导、绝热壁面温度和壁面摩擦力线分布的计算精度.计算结果表明:两道斜激波相交后的区域的壁面压强和热传导都比较高;计算的压强和壁面摩擦力线与试验吻合很好,绝热壁面温度次之,热传导最差,峰值高达试验的3倍左右.湍流模型对壁面压强和壁面摩擦力线影响很小,对绝热壁面温度和热传导影响很大.在5种湍流模型中,TNT(turbulent/non-turbulent)和SST(shear-stress transport)模型表现较好.      

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

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

A new simpler rotation/curvature correction method for Spalart-Allmaras turbulence model

 A new and much simpler rotation and curvature effects factor, which takes the form of Richardson number suggested by Hellsten originally for SST k-x model, is presented for Spalart and Shur's rotation and curvature correction in the context of Spalart-Allmaras (SA) turbulence model. The new factor excludes the Lagrangian derivative of the strain rate tensor that exists in the SARC model, resulting in a simple, efficient and easy-to-implement approach for SA turbulence model (denoted as SARCM) to account for the effects of system rotation and curvature, techniquely. And then the SARCM is tested through turbulent curved wall flows: one is the flow over a zeropressure- gradient curved wall and the other is the channel flow in a duct with a U-turn. Predictions of the SARCM model are compared with experimental data and with the results obtained using original SA and SARC models. The numerical results show that SARCM can predict the rotation-curvature effects as accurately as SARC, but considerably more efficiently. Additionally, the accuracy of SARCM might strongly depend on the rotation-curvature model constants. Suggesting values for those constants are given, after some trials and errors.     

Numerical simulation of compression corner flows at Mach number 9

A hypersonic vehicle encounters a wide range of conditions during its complete flight regime. These flight conditions may vary from low to high Mach numbers with varying angles of attack. The near-wall viscous dissipation associated with flows at combined high Mach and Reynolds numbers leads to significant wall heat transfer rates and shear stresses. The shock wave/boundary-layer interaction results in a flow separation region, which commonly augments total pressure losses in the flow and lowers the efficiency of aerodynamic control surfaces such as fins installed on a vehicle. The standard turbulence models, when used to resolve such flows, result in incorrect separation bubble size for large separated flows. Therefore, it results in an inaccurate aerodynamic load, such as the wall pressures, skin friction distribution, and heat transfer rate. In previous studies, the application of the shock-unsteadiness correction to the standard two-equation k-ω turbulence model improved the separation bubble size leading to an accurate pressure prediction and shock definition with the assumption of constant Prandtl number. In the present work, the new shock-unsteadiness modification to the k-ω turbulence model is applied to the hypersonic compression corner flows. This new model with variable Prandtl number is based on the model parameter, which depends upon the local density ratio. The computed wall pressures, heat flux and flow field are compared to the experimental data. A parametric study is carried out by varying compression deflection angles, free stream Reynolds number and wall temperatures to compute the flow field and wall data accurately, particularly in the shock boundary layer interaction region. The new shock-unsteadiness modified k-ω model with variable Prandtl number shows an accurate prediction of initial pressure rise location, pressure distribution in the plateau region and heat flux in comparison to the standard k-ω model.     

扩展型面排孔气膜冷却的实验

设计了三种型面的气膜孔扩展结构,研究了气膜冷却效率、壁面换热系数及流量系数随吹风比的变化规律.研究结果表明,与典型的圆柱形气膜孔相比,有后向扩展型面的气膜孔具有相对较高的冷却效率和换热系数,流量系数略有降低.扩展型面的流向扩展角对冷却效率和换热系数的影响大于侧向扩展角的影响.侧向扩展角只在特定吹风比下才对提高气膜的侧向覆盖率、进而提高气膜冷却效率起作用,它使得位于二次流出流下游的低压区域范围增大,流量系数减小.      

水平喷流对叶顶气膜冷却效率影响的数值研究

为进一步提高燃气轮机叶顶区域的气膜冷却效率,在叶顶模型结构上对4种叶顶冷却结构的流动和传热进行了数值模拟,结果显示带有盖板的水平喷流结构能够有效地提高叶顶区域的气膜冷却效率,对于圆孔,可以将气膜冷却效率提高至原来的2倍以上.通过流场的分析发现:水平喷流结构可以调整冷却工质的流动,使得冷却工质的贴壁性更好、分布更加均匀,进而提高气膜冷却效率.此外还研究了间隙宽度、盖板尾缘厚度、盖度和整流栅对水平喷流结构气膜冷却效率的影响.结果表明:增加盖度和减小间隙宽度可以提高气膜冷却效率,但是盖度和间隙宽度的选取受到了强度、工艺和冷却工质入口压力等因素的限制.水平盖板结构的冷却效果比倾斜盖板结构的好.整流栅足够长时可以调整流动、提高气膜冷却效率.      

不同缝槽进气孔间距的纯气膜绝热温比试验研究

通过传热试验法对不同缝槽进气孔间距的纯气膜冷却效率进行研究。研究所用气膜缝槽结构无收缩,进气孔与舌片垂直。试验中主流湍流强度在10%以上,测量气膜绝热壁温以确定其重要冷却效率描述参数——绝热温比。试验结果表明:缝槽进气孔间距对气膜冷却效率影响很大,减小进气孔间距,可以减小缝槽几何特征参数MIX_N,明显提高气膜的绝热温比。      

主流湍流度对涡轮导向叶片气膜冷却特性影响的实验

采用基于窄带热色液晶的瞬态全表面传热测量技术,研究了主流湍流度对涡轮导向叶片吸力面圆柱形孔排气膜冷却特性的影响规律.结果表明:在实验工况范围内,主流湍流度从0.59%提高至6.85%,可以在气膜出流的上游区域促进气膜贴向壁面并扩大展向覆盖面积,从而改善气膜覆盖效果,但是在主流湍流度较大的工况下,气膜覆盖效果迅速变差;在气膜出流的下游区域,主流湍流度的提高使得气膜冷却效率逐渐降低;主流湍流度的增大,增强了无气膜冷却光滑叶片表面的对流换热;在气膜冷却条件下,气膜出流对叶片表面对流换热的增强效果随着主流湍流度的增大呈现出明显的区域性特点:表面传热系数比在上游区域是先增强后减弱;中游区域是逐渐减弱;下游区域则是逐渐增强.     

高压燃气涡轮径向内冷叶片气热耦合的数值分析

采用气热耦合方法对采用径向内冷方式的MarkⅡ型跨声速高压燃气涡轮金属导叶进行数值模拟,通过分析叶片通道内的传热和流动过程发现叶片表面附面层内流动非常复杂,包含层流流动、转捩和湍流流动状态,所以只有使用转捩模型计算的叶片附面层内流动与实际情况相符,叶片壁面温度和换热系数分布与实验结果吻合的较好,使用其他湍流模型由于不能准确描述附面层内流动而使得计算结果误差相对较大,但是所有的湍流模型都能很好的模拟附面层以外流动.      

涡轮叶栅流动和传热耦合计算

通过求解三维雷诺平均N-S(Navier-Stokes)方程和带转捩模型的二方程SST(Shear stresstransport)湍流模型,完成了带简单冷却的MARKⅡ涡轮导向叶栅流动和传热耦合计算.计算结果与试验数据的对比表明,发展的带转捩模型的数值方法明显地提高了叶栅温度和外换热系数的计算精度.同时研究发现,湍流模型对叶栅表面压力分布影响很小,而对叶栅表面温度和外换热系数影响很大;在其它边界条件相同的情况下进口湍流度对壁面压力几乎无影响,但对叶栅温度和外换热系数影响较大.