插板式畸变发生器后非均匀流场结构数值模拟

随着计算机技术和计算方法的迅速发展,三维非定常方法已经被直接用来求解非均匀进口条件下的压气机流场.但是考虑到网格数量和计算周期等因素,计算往往以压气机前作为进口,对进口畸变形式进行简化.然而这种简化与实际压气机进口非均匀特性存在一定出入,降低了数值模拟的准确性.为了尽可能真实地反映进口畸变形式,同时又不增加计算域的网格数,针对插板式畸变发生器制造的非均匀流场结构开展数值研究,分析了插板后不同位置的流场结构和畸变度,并研究了压气机前总压分布随插板深度和来流总压的变化规律,进而给出了压气机前总压分布形式的规律性描述.     

激光吸收光谱技术测量非均匀燃烧流场研究进展

可调谐半导体激光吸收光谱技术(TDLAS)具有非侵入性、灵敏度高和时间响应快等优点,将激光吸收光谱技术与最小二乘法、计算机断层扫描重建技术(CT)相结合,可以实现对非均匀燃烧流场的分布测量。首先简要介绍了激光吸收光谱技术的发展历程及测量基本原理,然后分别对激光吸收光谱技术测量非均匀流场一维分布、二维分布的国内外研究现状及关键技术进行了综述,比较分析了二维非均匀流场诊断实验中旋转和固定两种安装模式的优缺点及相对应的光线布局,总结了用于流场二维重建的相关重建算法,最后讨论了激光吸收光谱技术测量非均匀流场研究工作的发展趋势和有待解决的相关问题。     

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

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

二维非均匀超音来流下最有利压缩面型面的随机选取法数值研究

采用近年发展起来的随机选取法，研究非均匀二维超音来流流过６种不同形状的楔型压缩面时的波后非均匀超音流场。结果表明来流的非均匀性对下游流场具有重大影响，对６种不同形状楔型压缩面流场的计算表明，单纯用改变压缩面型面的办法难以减少来流的非均匀性，但是可以部分改善下游流场的非均匀程度，在研究的６种型面中，凹圆弧压缩面及等熵压缩面在本文限定的非均匀来流下，其性能优于其它４种型面的压缩面。     

基于Mohring声类比的开式转子发动机噪声分析

为了系统地研究非均匀流场中开式转子的噪声特性,运用计算流体动力学(CFD)和计算气动声学(CAA)无限元法相结合的方法,对开式转子的频域噪声做了数值预测。通过采用基于可压缩修正的Spalart-Allmaras湍流模型对非均匀流场中开式转子气动特性进行分析,运用声学无限元方法,引入平均流场并将其插值到声学网格作为背景流场,对转子噪声和噪声传播进行数值模拟,并对特征点做频域分析。由噪声分布云图和远场特征点频谱曲线分析得出:非均匀流场中开式转子发动机的峰值辐射噪声主要集中于低频范围,噪声最大值为142dB,峰值噪声频率为115Hz,随着频率增加,交替出现多个局部峰值声压,噪声衰减速度减小。为下一代开式转子发动机噪声预测提供一种方法。      

亚声速风洞非均匀流场的修正计算

本文给出了一种修正风洞非均匀方向场和速度场的方法。对于任意的机翼-尾翼布局和风洞剖面的计算结果和实验数据作了比较,两者吻合得相当好。计算证明在非均匀流场中俯仰力矩曲线可能发生畸变。本方法为修正试验结果提供了一个手段。     

非均匀栅距对轴流压气机流动影响的数值研究

本文进行了非均匀栅距叶片排的轴流式压气机流场数值研究,进口导叶采用了不同的非均匀栅距结构,探讨了不同工作状态下对压气机性能及动叶非定常气动力的影响。数值研究中采用LU-SGS隐式解法与双时间方法相结合,求解圆柱面二维流动的雷诺平均非定常N-S方程。数值研究的结果表明,进口导叶非均匀栅距结构对压气机气动性能影响不大,但对动叶表面受到的非定常气动力有一定影响。采用不同的非均匀栅距进口导流叶片排结构型式,动叶表面的非定常气动力的分布是不同的。不同工作状态下的计算结果显示,即使采用相同的非均匀栅距结构,非定常气动力的分布特性也是不相同的。     

有限差分法的坐标变换诱导误差

有限差分法应用于具有复杂外形的网格时需要进行坐标变换,在此过程中经常会引入坐标变换诱导误差。在柱坐标系下使用均匀网格进行均匀流场计算,计算结果表明,即使物理坐标对计算坐标的变换函数连续可导、计算过程中坐标变换系数直接采用准确的解析式、网格完全正交并且充分光滑,也无法避免坐标变换诱导误差。理论分析表明,产生坐标变换诱导误差的机理是笛卡尔坐标系下的守恒型欧拉方程变换至贴体坐标系下后增加了源项。针对该问题,目前国内外学者通常采用几何守恒律,构建与差分格式相匹配的坐标变换系数计算方法来消除源项。本文介绍了从包含源项的离散等价方程基础上直接进行离散的新算法,在此基础上针对非等距网格条件下MUSCL类格式重构过程进行误差分析,理论推导表明重构中需要考虑非等距插值公式的影响系数,将变量转换至计算空间内进行MUSCL重构才能保证该过程具有均匀网格下的插值精度。通过理论分析及数值实验证明新算法对于均匀流场完全不会引入坐标变换误差。     

吸气式高超声速飞行器非均匀尾喷流试验

在中国空气动力研究与发展中心0.5m高超声速风洞中,开展了非均匀喷流条件下的吸气式高超声速飞行器后体尾喷流/外流干扰测压试验研究。采用非均匀内喷管,模拟飞行器尾喷管非均匀入流,测量了飞行器后体膨胀面及水平翼表面压力,采用高清纹影观测了喷流干扰区域的流场结构,获得了不同工况下非均匀入流对尾部及水平翼表面压力分布的影响规律。试验结果显示尾喷管非均匀入流对飞行器尾部壁面压力分布及流场结构有明显影响,喷管入流的非均匀特征在吸气式高超声速飞行器喷流模拟中不可忽视。非均匀喷流核心区压力分布明显高于均匀喷流时的结果;核心区域外,非均匀喷流的作用面积略小于均匀喷流,且非均匀喷流同外流交叉干扰区域的面积和强度要略小于均匀喷流;均匀喷流在喷管出口区域存在明显的膨胀波系,交叉干扰激波及剪切层的扩张角也大于非均匀入口条件时的结果。     

二维非均匀射流遇障碍时波系的变化

本文研究了非均匀射流在二维管道内的钝头体绕流问题。所采用的Godunov格式适合于不规则区域上双曲守恒系统的流场计算,数值计算的结果给出该流场的主要特征和激波结构,对于认识这一物理现象和了解波的运动规律是有帮助的。     

轴流压气机机匣变形对多排转子流场特性的影响

在实际轴流压气机加工、装配、使用过程中,机匣的圆度控制、叶片的高度控制、转子和机匣之间的同轴度控制等都普遍存在误差,意味着叶尖间隙的不均匀性必然地存在于每台压气机中。随着对轴流压气机流动研究的不断深入,叶尖周向非均匀间隙成为了需要考虑的问题。针对由机匣变形造成的非均匀叶尖间隙,研究了其对多排转子流场特性的影响。首先介绍了表征机匣变形程度的新参数,进而引进了非轴对称压气机模型的建模方法。针对3种不同的间隙周向布局,采用定常和非定常方法进行了数值模拟,对比了均匀与非均匀间隙下多排转子的气动性能,分析了机匣变形对气动损失分布和传播的影响,并研究了非均匀间隙下非定常压力和气动力的脉动特征。结果表明：非均匀间隙会降低转子性能,转子叶排叶尖区域的流场在非均匀间隙下呈现明显的周向非对称性,在上游叶排“尾迹”的干扰下,后排转子叶尖流场的周向非对称性要大于前排转子。非均匀间隙布局与各叶片气动力分布呈现明显的对应关系,小间隙区域叶片的气动力高,大间隙区域叶片气动力相对较低,从而增加了非均匀间隙下转子叶片的气动力幅值。     

入口非均匀流对亚燃燃烧室性能影响数值研究

为研究人口非均匀流对亚燃燃烧室性能的影响,基于进气道出口试验测量数据,构造三种典型非均匀流型,对燃烧室进行三维数值计算,研究非均匀流对马赫数分布、温度场、燃油分布及燃烧效率的影响。结果表明：基于均匀流供油律,非均匀流影响燃烧室性能的机理是受流速主导作用,低速区域当量比偏高,高速区域当量比偏低。导致扩压段边区燃烧,烧蚀外环阻塞结构,而内圈高速区域燃烧组织效果欠佳。匹配供油律后,非均匀流燃烧性能有所改善,但难以达到均匀流水平。所以,提出了改进入口流均匀性及减弱非均匀流对燃烧不良影响的措施。     

考虑进口非均匀的喷管设计及冷流条件下的性能研究

超燃冲压发动机尾喷管与燃烧室直接相连,由于没有几何喉道和收缩段的整流作用,实际工作过程中尾喷管的进口气流是非均匀的,因而有必要研究非均匀进口条件下的尾喷管设计方法。以超燃冲压发动机尾喷管非均匀进口马赫数分布为条件,采用有旋特征线设计了超燃冲压发动机非对称尾喷管的等熵膨胀型线。利用数值模拟和试验相结合的方法研究了冷流条件下喷管的气动性能。结果表明,在相同的进口条件下,相对于假定进口马赫数均匀分布设计得到的喷管,考虑进口马赫数非均匀分布设计的喷管推力增加0.6%~2.0%,负升力降低可达82.0%,俯仰力矩增加8.6%~13.0%,这说明在喷管的设计过程中考虑进口参数分布的非均匀性是有必要的。     

非均匀来流下三维激波反问题的微元密切轴对称解法

传统的密切轴对称理论被广泛应用于均匀来流下的三维密切曲面激波反设计，为解决非均匀来流条件下的三维曲面激波反问题，提出了一种微元密切轴对称流场（MOA）求解方法。该方法沿激波面的周向和流向构建一系列微元密切面，在每个微元面内进行三维向二维流动的等效转换，从而突破了传统密切方法中不能有横向波后流动的限制。利用该方法编写设计程序，分别基于带攻角来流条件和外锥型流来流条件重构了标准内锥曲面激波，并与数值仿真结果进行了比较。结果表明，非均匀来流下激波曲面的三维形状均与预设形状完全一致，实现了非均匀来流下曲面激波形状可控。MOA方法在吸气式高超声速推进领域中前体/进气道一体化设计方面有重要应用前景。     

Stationary flow fields prediction of variable physical domain based on proper orthogonal decomposition and kriging surrogate model

In this paper a new flow field prediction method which is independent of the governing equations, is developed to predict stationary flow fields of variable physical domain. Predicted flow fields come from linear superposition of selected basis modes generated by proper orthogonal decomposition(POD). Instead of traditional projection methods, kriging surrogate model is used to calculate the superposition coefficients through building approximate function relationships between profile geometry parameters of physical domain and these coefficients. In this context,the problem which troubles the traditional POD-projection method due to viscosity and compressibility has been avoided in the whole process. Moreover, there are no constraints for the inner product form, so two forms of simple ones are applied to improving computational efficiency and cope with variable physical domain problem. An iterative algorithm is developed to determine how many basis modes ranking front should be used in the prediction. Testing results prove the feasibility of this new method for subsonic flow field, but also prove that it is not proper for transonic flow field because of the poor predicted shock waves.     

Adjoint boundary sensitivity method to assess the effect of nonuniform boundary conditions

The performance of turbomachinery is largely affected by the nonuniform boundary conditions caused by the coupling between neighboring parts, such as the inlet distortion and hot streak. Existing works study this problem by comparing the flow fields with uniform and nonuniform boundary conditions, which is cost extensive. In this work a new and efficient method is developed by computing the sensitivities of arbitrary performance metric relative to the boundary conditions and this method quantita...     

基于基尔霍夫面的四维声学频域公式

流动诱发噪声问题是实际工程领域极为普遍的问题之一,经典的声比拟模型仅以声学压力为参考来评估声场特征分布还远远不够。从声压和声学速度矢量为变量的四维线性波动方程出发,选择包围非线性声源的基尔霍夫面为积分面,并结合对流格林函数,给出均匀运动介质的四维声学频域积分公式;针对静止、旋转单极子源和偶极子源开展数值预测研究。结果表明：本文获得的声压和声学速度分布与解析解吻合,均匀来流情形下静止点源的声场分布表现出典型的对流效应;受均匀来流、点源的自激频率、谐波阶次和旋转频率的共同影响,旋转点源的声场分布则表现出明显的多普勒效应和对流效应。本文对均匀流矢量气动噪声开展的精细化研究能够为声能量的评估和传输路径预测提供技术支撑,进一步为降噪研究提供理论依据。     

基于物理知识约束的数据驱动式湍流模型修正及槽道湍流计算验证

对湍流摩擦阻力的精准预测是学术界和工业界普遍关心的重要问题，而数据驱动式的湍流模型修正方法对此显示出较大的潜力和前景。提出了一种基于物理知识约束的数据驱动式湍流模型修正方法，根据湍流摩擦阻力分解获得先验物理知识，在S-A湍流模型的生成项中引入非均匀分布的修正因子，以修正因子为设计变量，设定包含物理知识约束的目标函数，利用离散伴随方法求解目标函数与设计变量之间的梯度关系，通过高效率的迭代求解获得修正因子的分布。以槽道湍流为例，验证了包含物理知识约束的数据驱动式建模方法的优势，并分析了物理知识约束对湍流摩擦阻力预测精度的影响，结果表明引入物理知识约束可进一步提高湍流摩擦阻力的预测精度。     

转动坐标系中三维跨声速欧拉流的有限体积 TVD格式

在非惯性转动坐标系中,本文采用贴体网格、有限体积法离散和修正数值通量技术,将Harten的一维TVD格式推广到三维。由于转动使方程出现源项,文中通过对源项的巧妙处理,使修改后的格式能用于非齐次双曲守恒律方程组高分辨率的数值计算。为了加速解的收敛,提高显式时间推进的CFL数,本文采用隐式残值光顺技术。三维跨声速带非齐次源项欧拉方程的典型算例表明:捕捉的激波分辨率较高;激波前后没有发现大的数值波动和伪振荡现象;所得的跨声速流场解与实验较接近。