H型叶栅通道流计算网格的改进

提出一种壁面正交、尾迹区网格节点分布逐渐均匀化叶栅通道网格生成方法。为了研究壁面处网格正交性对数值计算结果的影响，推导出壁面附近流动参数一阶导数离散误差与网格线交角的数学关系。采用N-S方程分别对平板紊流附面层流和叶栅通道流进行计算分析。结果表明：壁面处网格正交性对附面层内流动计算预测精度影响较大，而对势流区流动影响很小；尾迹区切向网格节点逐渐均匀化分布能更准确计算尾迹区流动。     

亚音速扩压器性能计算

 一、组合法 用Bower的方法计算具有附着和分离紊流附面层的轴对称和二元亚音速扩压器的有关性能参数时,如分离区较大,则与试验数据相比误差较大。而按联合积分法（UIM法）计算具有分离区的扩压器的有关性能参数,则有较大的优点。其基本点和Bower的方法相同,即将附面层方程和无粘核心流方程联立求解。但对分离区的处理比Bower的方法较为合理,它直接采用修正的Coles壁面-尾流定律的速度分布关系式,将附面层动量积分方程中的壁面摩擦系数与速度分布关系式中的摩擦速度联系起来，而不采用对分离区不适用的壁面摩擦系数经验公式，并考虑了附面层外部区域的非平衡流情况。但该方法没有考虑气流压缩性，计算高亚音速扩压器时误差较大。     

火箭发动机推力室附面层的研究

文本采用由壁面的切线和法线构成的正交曲线坐标推导出了纵向大曲率壁面的附面层方程.它适于轴对称和平面的,内和外流的可压缩湍流和层流附面层.作者从附面层方程出发,发展了内区涡动粘性系数和内区涡动热传导系数的表达式.此表达式本身已考虑了曲率的影响,不需要对它另作曲率修正.本方法针对火箭发动机模型,取不同的喷管收敛段曲率半径进行计算.此无因次曲率半径等于1.7时的喉部区峰值热流要比无因次曲率半径等于10时峰值热流高13%.此外在同样的壁型面和同样的压力分布下,考虑曲率和不考虑曲率的计算结果差别明显,喉部热流差别达9.8%.     

非均匀来流下等压力梯度曲面压缩系统的试验研究

根据激波和膨胀波理论,设计了Ma=5.3等压力梯度曲面压缩系统。在Ma=3.85和Ma=5.3的小高超风洞完成均匀来流试验,并利用风洞壁面自然发展的附面层进行非均匀来流等压力梯度曲面压缩系统性能试验。试验发现:均匀来流下,试验模型壁面压力均按照等压力梯度规律分布,但来流马赫数Ma=5.3时出口流场均匀性较差;非均匀来流下,大量的附面层吸入不仅影响壁面压力不再呈等压力梯度分布,而且对出口流场均匀性影响也极大。      

复杂外形气动力计算的阻力修正方法

基于气动力计算的位流线化理论,发展了一种亚音速附面层阻力修正方法。采用联合流场求解方法获得复杂外形的表面压力分布,在此基础上,将绕复杂外形物体沿流向面元划分为片条,在片条上,把三维附面层转换为二维附面层,对二维附面层进行数值求解,得到附面层粘性阻力修正。算例表明,该方法计算量小,与实验和N-S数值解结果吻合较好。为复杂气动外形飞行器初期设计、方案选型及评估提供了较为理想的工程设计手段。     

三维并联式TBCC发动机排气系统设计与实验

为实现三维并联式TBCC排气系统的设计,利用基于二维最大推力理论的密切法设计了圆转矩变截面非对称冲压发动机喷管,并在冲压发动机喷管上壁面为三维曲面的条件下,利用转动加滑动共用面的调节方式实现了涡喷发动机喷管的内型面设计及其喉道面积的调节,从而完成了排气系统设计。随后对设计结果进行了冷流实验验证。结果表明,这种排气系统设计方法具有可行性。在冲压单独工作状态,壁面压力分布的数值计算与实验结果吻合较好,相对方均根偏差低于10.8%;在涡喷单独工作状态及共同工作状态,上壁面压力在激波附面层干扰区域内偏差较大,最大相对偏差可达30.6%。但壁面压力分布的数值计算与实验结果的整体变化趋势一致。      

尾迹/势流干扰对跨声速风扇静叶附面层流动影响的数值研究

为揭示尾迹/势流干扰下压气机静叶附面层的非定常流动机理,采用数值方法对单级高负荷跨声速风扇中径处的非定常流场进行模拟,针对尾迹/势流干扰下的静叶附面层非定常流动特征进行研究。基于尾迹/势流与叶片附面层干扰模型,通过叶片壁面摩擦力、近壁面附面层湍动能和壁面静压脉动,详细分析了尾迹和势流干扰下尾迹对高负荷静叶附面层流动状态的影响。研究发现:在非定常条件下,尾迹干扰能够"刺穿"静叶附面层,使得静叶附面层发生跨越转捩现象。尾迹和受其诱导的势流碰撞静叶前缘产生的压力扰动波在压力面附面层内以声速传播,影响压力面整个弦长的静压和摩擦力分布,而在吸力面上,扰动波的传播仅局限在前缘区域部分弦长上,吸力面和压力面传播现象不同与逆压梯度和曲面凸、凹形状相关。     

用湍流附面层方程计算喷管内的对流换热

在流场计算的基础上,运用Patankar-Spalding方法积分求解了湍流附面层控制方程组,从而得到了固体火箭发动机喷管内燃气与壁面之间的对流换热。其中在附面层靠近壁面的底层区使用了Gouette流分析,在壁面边界和自由流边界处引入了滑移值的概念,因此保证了计算结果的相对准确性。最后将本文的计算结果与巴兹公式的计算结果和实验结果进行了比较。     

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

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

关于两方程湍流模型的考虑

通过对两方程湍流模型的分析,提出了在主流区或远壁区用两方程模型,近壁区用代数湍流模型的思路,即所谓的“双层”模型概念。通过对平板湍流附面层的计算表明,这种分层计算措施获得的速度分布、壁面剪切力较单纯应用低雷诺数两方程湍流模型时的结果准确。对非均匀来流条件下90°弯管内高湍流度流场的计算表明,双层模型可获得更加准确的速度分布。从而表明本文提出的分层计算措施是有效的,可满足工程计算要求。      

马赫数2.5超/亚声流流向涡掺混超燃室冷态实验

结合对应的数值计算,对一种带亚声速预燃室和流向涡掺混器的超声速燃烧模型燃烧室实验台,在其进口马赫数为2.5的来流条件下,进行了冷态流场的实验研究.实验测得其壁面静压分布和激波系结构与流场的CFD计算结果基本一致.实验结果表明,模型燃烧室全流场超声速,达到设计状态.马赫数2.5下的冷态实验数据和CFD计算数据为进行点火实验提供了依据.      

A study on turbulence transportation and modification of Spalart–Allmaras model for shock-wave/turbulent boundary layer interaction flow

It is of great significance to improve the accuracy of turbulence models in shock-wave/ boundary layer interaction flow. The relationship between the pressure gradient, as well as the shear layer, and the development of turbulent kinetic energy in impinging shock-wave/turbulent boundary layer interaction flow at Mach 2.25 is analyzed based on the data of direct numerical simulation(DNS). It is found that the turbulent kinetic energy is amplified by strong shear in the separation zone and the adverse pressure gradient near the separation point. The pressure gradient was non-dimensionalised with local density, velocity, and viscosity. Spalart–Allmaras(S–A) model is modified by introducing the non-dimensional pressure gradient into the production term of the eddy viscosity transportation equation. Simulation results show that the production and dissipation of eddy viscosity are strongly enhanced by the modification of S–A model. Compared with DNS and experimental data, the wall pressure and the wall skin friction coefficient as well as the velocity profile of the modified S–A model are obviously improved. Thus it can be concluded that the modification of S–A model with the pressure gradient can improve the predictive accuracy for simulating the shock-wave/turbulent boundary layer interaction.     

高超飞行器内流道激波振荡问题的数值研究及试验验证

采用数值模拟和风洞试验的方法,对吸气式高超声速飞行器盲腔状态的流动特性进行了研究。采用＂双时间步＂方法进行了内外流一体化的非定常数值模拟,利用彩色纹影系统对高超声速飞行器前体流场进行显示,并采用动态压力传感器测得了飞行器内流道的壁面压力随时间的变化。结果表明：盲腔构型在高超声速飞行中会出现周期性的激波振荡现象。数值模拟所得流场变化特征、内流道壁面压力振荡周期和壁面压力变化趋势与试验结果吻合良好。     

高超声速进气道-隔离段反压引起不起动计算

1引言隔离段是高超声速进气道不可缺少的组成部分,它是进气道与燃烧室之间的气动热力缓冲段。隔离段的主要作用是承受下游燃烧室的反压变化而不影响上游进气道的流态,避免进气道出现不起动。隔离段应在承受反压变化的同时不破坏上游进气道流场,组织好隔离段内所形成的激波串,高     

CFD/CSD-based flutter prediction method for experimental models in a transonic wind tunnel with porous wall

To predict the flutter dynamic pressure of a wind tunnel model before flutter test, an accurate Computational Fluid Dynamics/Computational Structural Dynamics (CFD/CSD)-based flutter prediction method is proposed under the conditions of a 2.4 m × 2.4 m transonic wind tunnel with porous wall. From the CFD simulations of the flows through an inclined hole of this wind tunnel, the Nambu’s linear porous wall model between the flow rate and the differential pressure is extended to the porous wall with inclined holes, so that the porous wall can be conveniently modeled as a boundary condition. According to the flutter testing approach for the current wind tunnel, the steady CFD calculation is conducted to achieve the required inlet Mach number. A time-domain CFD/CSD method is then employed to evaluate the structural response of the experimental model, and the critical flutter point is obtained by increasing the dynamic pressure step by step at a fixed Mach number. The present method is applied to the flutter calculations for a vertical tail model and an aircraft model tested in the current transonic wind tunnel. For both models, the computed flutter characteristics agree well with the experimental results.     

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

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

开槽叶片对大转角扩压叶栅性能的影响

采用从压力面向吸力面开槽的局部流动控制方法,设计了一种收敛转折型的槽道结构.实验对不同冲角下开槽叶栅的进、出口流场进行了测量,利用实验结果对数值模拟结果进行了校核,通过数值计算进一步得到了详细的叶栅通道内流场情况,并进行了结构静力分析.结果表明:在4°进气攻角下,开槽后叶栅尾迹区宽度减小了16.7%,总压损失系数峰值减小了6.07%;在6°进气攻角下,总压损失系数峰值减小了14.7%.叶片开槽从压力面吸入的气流可有效加速吸力面附面层流动,抑制吸力面分离,从而降低总压损失,增大静压比,扩大稳定工作范围.槽道前壁面的转折处存在应力集中,需要进行改进.      

大扩张角涡轮过渡段性能试验和数值研究

为研究某型大扩张角涡轮过渡段气动性能,对过渡段内部流场进行了详细的试验测量,同时采用CFD数值模拟对过渡段内部流场进行仿真,并与试验结果进行对比分析.结果表明:过渡段机匣表面流动受强逆压梯度影响,容易发生流动分离;轮毂表面流场受支板前缘冲击绕流的影响,呈现周向不均匀性.来流气流角使得过渡段内部流场向支板一侧偏斜,随着气流角的增大,过渡段总压损失增大.CFD模拟结果与试验测量结果吻合较好,均能很好地捕捉流场的细节特征;过渡段进、出口总压恢复系数随着来流气流角的增大而减小,CFD模拟和试验测量值的偏差约为0.2％.     

S型收缩流道对涡轮通流能力影响的研究

为进一步提高航空发动机涡轮通流能力，以小型跨声速涡轮为原型，研究了静子S型收缩流道对涡轮通流能力的影响。将S型流道直线段长度和收缩段内外圆半径比作为流道造型参数，分别针对原型涡轮静子内外端壁进行调整造型，得到一系列不同参数组合的S型内外端壁涡轮算例。保持膨胀比为设计值不变，利用CFD软件对原型和S型流道涡轮进行设计点模拟分析。结果表明，S型流道涡轮流量提升的原理在于增大的静子喉道面积。在相同造型参数下，S型外端壁涡轮的静子叶根损失被有效降低，流量提升明显，且流量高于S型内端壁涡轮1%左右，但由于最大外径增大使其质量通量提升效果减弱；与之相反，S型内端壁涡轮的质量通量提升明显，且高于S型外端壁涡轮3%左右。从提升涡轮质量通量并保证效率不低于原型的角度看，S型外端壁造型参数选取范围更广。     

基于OpenFOAM的三维H2/Air连续旋转爆轰流场数值模拟

为进一步研究旋转爆轰流场特征,基于开源计算流体动力学软件OpenFOAM,采用9组分19步的基元化学反应模型,对H2/Air连续旋转爆轰流场进行了三维数值模拟,得到了旋转爆轰波稳定传播时燃烧室内部流场的详细结构,研究了燃烧室头部激波的传播特性,分析了旋转爆轰燃烧室的压力增益性能。结果表明:旋转爆轰波后的第一道反射激波在由燃烧室外壁面向内壁面传播过程中反射激波的高度增加并在靠近内壁面附近与滑移线交汇形成局部高温高压区域;旋转爆轰波在外壁面位置处相位约落后于内壁面0.003rad～0.15rad,其相位差随燃烧室曲率差的增加而增大;燃烧室头部反射激波数目受到曲率差和进气总压的影响,燃烧室曲率差增大,反射激波数目减少,进气总压增大,反射激波数目增多;燃烧室压力增益保持在0.3以上,在进气总压一定的条件下,压力增益随着燃烧室曲率差的增大有增加的趋势。研究结果揭示了三维旋转爆轰流场的精细结构和燃烧室头部激波的传播规律。