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

针对两种双尖鳍外形的对称交叉激波与湍流边界层相互作用,采用N-S方程和两种湍流模型进行了计算。研究了网格收敛性、鳍的角度和湍流模型对壁面压强、Stanton数和壁面摩擦力线的影响。弱相互作用的计算结果较好,强相互作用的壁面压强和摩擦力线的计算结果与试验吻合较好,而Stanton数的计算结果较差,峰值高达试验的2.5倍左右。随着鳍的角度的增加,壁面压强和Stanton数的分布从单调分布发展为M型分布,两者的峰值不在相同的位置。湍流模型对壁面压强和壁面摩擦力线影响很小,对Stanton数计算的影响很大,SST模型比BSL模型表现好一些。     

非线性EASM在激波与湍流边界层相互作用中的数值应用

采用基于k-ω湍流模型的非线性显式代数应力模型(EASM)对超燃冲压发动机常用的超声速凹槽、压缩拐角和侧壁压缩进气道简化模型的激波与湍流边界层的相互作用进行了计算,主要研究了EASM 模型对壁面压强、摩阻、Stanton数和壁面摩擦力线的计算精度,计算结果与SST(shear stress transport)模型进行...     

三级压缩楔前缘半径对壁面静压及热流分布影响研究

为了研究三级压缩楔前缘半径对壁面静压和热流分布的影响,在0.6m激波风洞上开展了试验测量,模型长约0.6m前缘半径为0mm3mm试验名义马赫数为5.98。研究结果表明,试验得到的壁面静压和热流数据重复性很好。采用Fluent软件进行了二维和三维流场参数辅助模拟分析,三种不同湍流模型获得的壁面静压分布差别不大,均与试验结果吻合较好;不同湍流模型获得的壁面热流分布差异较大,采用标准k-ε模型得到的结果与试验吻合较好。试验和数值模拟结果均表明,在第二和第三压缩面上,经过激波后,壁面静压逐渐上升到一个压强平台;壁面热流逐渐上升到一个局部极大值,然后在同一压缩面持续下降。随前缘半径增加,壁面静压和壁面热流整体减小,压强平台值和热流局部极大值也减小,而达到压强平台和热流局部极大值需要的长度增加,显示激波边界层干扰影响区域增大。     

进气道中激波-边界层干扰特性研究

针对高超进气道中激波/边界层干扰等现象,选取了24°压缩拐角及二维高超进气道DLR-GK01等算例,比较研究了混合网格并行数值模拟软件中SA、SST、TNT三种湍流模型对激波/边界层干扰等现象的模拟能力,主要考察了不同湍流模型对流场结构、壁面压强以及壁面摩阻的影响。结果表明:三种湍流模型在模拟大分离流场时存在不同的流动分离特性,TNT湍流模型较SST和SA湍流模型模拟的分离起始位置、分离区大小及边界层厚度与实验吻合度最好。     

Hyshot超燃冲压发动机的CFD数值模拟

以Hyshot超燃冲压发动机试验为研究对象,分别对冷流和燃烧工况进行了CFD数值模拟.冷流模拟计算得到的壁面压力与试验结果吻合良好,且数值计算压力分布对采用的湍流流动模型不敏感;但是采用不同的湍流流动模型计算的湍流参数(湍动能和耗散率)差别很大,会对氢气燃料与气流的掺混产生重要影响,进而影响燃烧模拟结果.采用SST(shear stress transport)k-ω湍流流动模型和EDM(eddy dissipation model)湍流燃烧模型得到的模拟结果与试验结果基本吻合,但是计算压力峰值略靠后.对8个不同氢气燃料喷射角度的工况进行数值模拟结果表明喷射角度为99°和114°时燃料和空气混合最好、燃烧效率可达到最高.     

压缩拐角激波/边界层干扰的可压缩湍流模型研究

为研究压缩拐角激波/边界层干扰问题,抓住可压缩流动中的密度变化特性,利用构造的可压缩Von Karman尺度,基于KDO(Kinetic Dependent Only)湍流模型,发展出可压缩湍流模型CKDO(Compressible Kinetic Dependent Only)。通过对8°,16°,20°和24°压缩拐角算例的数值模拟,测试了其对可压缩、激波/边界层干扰这一湍流难题的求解能力。计算结果表明,总体上CKDO模型对壁面压力和壁面摩擦阻力系数的捕捉能力优于其它模型,并且随着压缩拐角角度的增大,其描述该流动的能力更加突出。CKDO模型在24°压缩拐角处计算的分离区大小仅比实验大10%左右,明显比其它模型结果好。这表明CKDO模型在模拟激波/边界层干扰这一类流动中有较好的适用性。     

三维非对称斜劈高速流动计算方法研究

交叉激波是高超声速进气道发展的关键气动问题之一。为实现对此类高速流动的准确数值模拟,基于RANS方程的自研多块结构网格计算软件,对经典三维非对称斜劈构型的斜激波交叉干扰流动特征开展了计算研究。通过对比分析,系统评估了计算网格、湍流模型及高精度格式对空间流场、底板摩擦力流线、压强分布、绝热壁温分布等的影响,获得了适用于三维非对称高速流动问题的网格生成准则和计算策略,同时验证了自研软件的模拟能力。     

壁温对氢燃料超燃燃烧室流动特性的影响研究

为了研究超燃燃烧室壁面换热与超声速燃烧之间的关系,运用FLUENT软件,采用RNG k-ε湍流模型、有限速率/涡耗散燃烧模型、密度基AUSM+隐式算法对中国空气动力研究与发展中心的双模态燃烧室模型开展三维冷态和热态流场计算,模拟条件来流马赫数为2.05,总温T_t为1870.9K,分别模拟了壁面温度为500K,700K,900K,1300K以及绝热条件下的超燃燃烧室的燃烧场。燃烧室壁面压力数值模拟结果与实验结果吻合较好,壁面温度为500K,700K,900K,1300K和绝热时,平均误差分别为8.89%,5.78%,14.41%,13.97%,16.53%。通过对比分析发现:随着壁面温度的降低,壁面压力趋势大致不变,但壁面压力值降低,同时壁面压力的压升起始点大幅后移;燃烧所产生的激波串逐渐向燃烧室下游移动,激波串结构发生改变,但激波串前端均为X形激波;燃烧室内马赫数有所升高;燃烧场高温区域面积减小;燃烧室燃烧模态由亚燃模态逐渐向超燃模态转换。     

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

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

TSTO并联分离激波/边界层干扰流动特性分析

针对两级入轨飞行器的缩比模型，通过试验与数值模拟相结合的方式，在马赫数6条件下开展典型级间距状态的激波/边界层干扰流场研究，详细分析干扰区壁面及空间的流动结构与特性。结果表明：试验中模型壁面边界层在激波入射之前为层流状态，在强激波干扰后迅速转捩为湍流状态，因此试验结果在第一道激波作用结束之前与层流计算结果吻合，而在第一道激波作用结束之后与湍流计算结果一致；激波/边界层干扰呈现复杂的三维流动特征和明显的开放结构，强激波在壁面形成的高压区呈弧状向下游展开，轨道级头部产生的入射激波在级间来回反射，强度依次递减；同时，干扰区内存在展向弯曲的主分离线与再附线、沿流向排列的二次分离线与再附线、流动剪切形成的旋涡结构以及包括鞍点、结点、焦点在内的临界点；层流边界层受到激波作用形成的分离区明显大于湍流边界层，同时开放特征更为显著。     

两种湍流模式对绕椭球体三维分离流计算的比较

将新近的Shih Lumley(S L)模式应用于雷诺数为7 2×106的6∶1椭球绕流的全N S方程计算,作为参考也计算了Jonse Launder(J L)模式。在没有发生分离时,两种湍流模式计算的平均流场与实验结果符合较好。在分离发生后,J L模式计算的分离结构与实验差别很大,没有计算出二次分离线;S L模式的计算结果较好,计算得到的表面摩擦力线图谱中主要分离结构与实验结果较接近,分离区的表面压力系数分布也与实验值定性相符;在S L模式计算结果中,壁面附近出现小尺度的分离再附结构.     

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.     

流量突变条件下煤油流动与对流传热特性实验

对流量突变条件下RP-3航空煤油的圆管流动与传热过程进行了实验研究。燃料流动与传热实验是在超临界压力为3MPa、燃料温度为300～650K、雷诺数为3000～60000范围内进行的。实验同步测量了圆管煤油温度、壁温、油压的沿程分布,并应用基于动量与能量守恒的非定常控制体分析方法,获得了航空煤油壁面摩阻系数与传热努塞尔数的分布。实验结果表明:当煤油温度未超过临界值（约650K）、处于过压液态时,该实验条件下的入口流量的突增或者突减过程对摩阻系数与努塞尔数的影响很小。      

气体模型对高超声速再入钝体气动参数计算影响的研究

本文采用完全气体模型、振动激发气体模型、平衡气体模型、一温度非平衡气体模型、两温度非平衡气体模型和三温度非平衡气体模型进行了钝体高超声速绕流流场的数值计算。分析了各种气体模型对激波脱体距离、壁面热流、温度分布和密度分布等的影响。结果表明：采用两温度（或三温度）非平衡气体模型计算的激波脱体距离更接近实验情况；采用平衡气体模型计算的壁面热流量高，一温度非平衡气体模型完全非催化的壁面热流最低，而完全气体模型的壁面热流在两者之间。最后给出了各种气体模型的密度等值线比较图。     

竖直圆管内超临界碳氢燃料换热恶化的直径效应

利用Fluent对超临界压力下直径对碳氢燃料换热恶化的影响进行数值研究,湍流模化采用Launder-Sharma（LS）低雷诺数模型,物性采用广义对应态法则对RP-3替代燃料计算。计算条件：系统压力为3 MPa,进口温度为573 K,热流密度为500 kW/m²,质量流量为0.001 5和0.003 0 kg/s,直径范围为1~10 mm。正常换热条件下的计算壁温和实验结果基本吻合,证明了计算方法的准确度。结果表明：强制对流下小质量流量时直径越大,换热恶化程度更高且更提前发生,换热恶化是由定压比热容处于极大值后的急剧下降区导致的,大质量流量时直径与壁温成正比,无换热恶化发生;浮升力效应仅在小质量流量下起作用,随着直径增大而加强,给出RP-3流动换热时浮升力起作用的判据和不同直径下换热恶化的边界。     

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

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

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

短周期风洞中导叶表面压力和换热测量

在发动机典型雷诺数和压比状态下对一种放大导叶叶型进行了表面静压和换热测量.雷诺数对表面压力系数的影响较小,压比增大使压力系数减小,并且吸力面压力系数最低点后移.雷诺数增大时叶片表面传热系数增加,并且吸力面上边界层转捩位置提前.压比主要影响吸力面传热系数,小雷诺数时压比增大会推迟吸力面上边界层转捩点位置,大雷诺数且吸力面后半段为超声速流动时,增大压比使该区域传热系数降低.保持主流总温不变,叶片表面绝热壁温随叶栅压比增大而降低,相同压比下,叶片表面处于层流状态时绝热壁温比处于湍流状态时低.