弯曲叶片造型对涡轮叶栅作用力影响的非定常数值研究

为了研究弯曲叶片对涡轮叶片非定常气动力的控制作用,采用非定常数值模拟方法对高压涡轮两级的内部流场进行了计算,研究了第一级静叶采用正弯叶片时,汽轮机涡轮叶片非定常气动力和非定常压力场的变化.结果表明,采用弯曲静叶的设计方法,通过改变静叶出口处径向压力梯度和涡量分布,能够有效抑制动静叶片表面气动力的非定常脉动作用.但弯叶片给非定常压力场带来的改善是局部的,其对叶片中间截面改善最为显著.从整体上来看,正弯曲的静叶设计对于减小非定常流场中的压力脉动是有益的.     

基于动态边界的跨声速压气机过渡态三维模拟方法

针对航空涡轮发动机部件的过渡态,提出了一种基于压气机部件三维模型的过渡态性能计算方法,采用动态边界条件,以转速和出口静压为过渡过程的控制参数,二阶向后欧拉法求解三维非定常雷诺平均Navier-Stokes方程组,通过三维非定常模拟获得压气机的过渡态性能和内部非定常流场特征。以NASA跨声速压气机Rotor 67转子为模型,采用该方法模拟了压气机从60%设计转速加速至100%设计转速的过渡过程,获得了压气机过渡态性能及转子内部详细的激波结构与演变过程。对比通用特性结果,在整个过渡曲线上,总压比相对误差最大值小于6%,绝热效率相对误差最大值小于2%,验证了该方法的可靠性。结论表明:基于动态边界的三维非定常模拟方法能够准确模拟压气机的过渡态性能,并反映出过渡态的非定常流场详细信息。     

涡轮级轮缘封严内非定常流场的准三维LDV测量

研究涡轮级封严内流动和主流的相互作用、相互影响,进一步认识盘腔内燃气入侵等非定常现象,急需流场实验数据验证。针对轮缘封严内空间狭窄、光路安排困难等问题,本文发展了一种用于测量涡轮级轮缘封严内非定常流场的准三维LDV技术,即测量分两步,首先在机匣上安置二维测速系统用符合方式测出径向、切向速度及湍流度,再在轮毂内安置一维探头重复测量,从两次测量数据中计算出轴向速度;该技术借助锁相采样等技术能测出不同动静叶相对位置情况下封严及盘腔内的周期性非定常流场细节,本文给出了部分典型的动态流场测量结果。      

非定常环境中叶栅总压恢复系数的数值研究

利用以Ni格式和隐式双时间步法求解二维非定常N-S方程为基础的数值模拟方法研究了流动的非定常性对涡轮叶栅总压恢复系数的影响,并与文中提出的定常模型和准定常模型进行了对比。结果发现,非定常流动能够提高叶栅的总压恢复系数,但这是建立在与准定常流相比较的层面上,而和定常流比较时,非定常流只会使总压恢复系数下降。      

非线性谐波法在叶轮机械非定常计算中的应用

非线性谐波法(NLH)进行叶轮机械非定常计算,具有占用资源少、计算时间短等优点,本文探讨了该方法在叶轮机械三维瞬态非定常流场计算中的应用。采用非线性谐波法和传统双时间步法对跨风速风扇和液力涡轮分别进行非定常计算分析。分析表明,在相同计算机硬件条件下,两种方法达到相同的计算精度,非线性谐波法可节省40%的计算时间。探讨了关键参数谐波阶次对计算精度的影响,结果表明,对于高雷诺数进口条件下的跨风速风扇采用三阶谐波逼近即可达到与传统双时间步法相当的计算精度;对于低雷诺数进口条件下的不可压缩液力涡轮,采用二阶谐波逼近即可达到与双时间步法相当的计算精度。     

一种快速模拟振荡叶栅非定常流的数值方法

为了提高叶轮机内周期性非定常流数值计算效率,发展了一种计算振荡叶栅非线性非定常流的谐波平衡方法,将周期性非定常流动参数采用傅里叶级数表示,从而将传统的定常流动数值计算方法应用于非定常流问题求解中,大大提高了非定常流计算效率。计算结果表明,发展的谐波平衡方法与线性方法、非线性时域法计算结果吻合较好,且其计算效率比时域法高1个数量级。另外,系统研究了谐波阶次对计算效率与计算精度的影响,认为在气流未出现严重分离时,应采用低阶谐波平衡方法,以减少计算耗时。此外,从压气机振荡叶栅非定常流场计算结果可以看出,叶片吸力面诱导涡的生成频率与振荡频率并不一致,涡生成频率要高于振荡频率。     

类升力体外形俯仰阻尼特性数值研究

采用有限差分方法求解薄层近似的非定常Navier-Stokes方程，定常流场采用交替方向隐式分解的NND格式，非定常流场采用四步Runge-Kutta方法，在保证时、空二阶精度前提下引入变系数残值光顺技术提高非定常流场计算效率，复杂带翼外形的空间网格通过求解抛物化的椭圆型方程生成，最后在Etkin理论下给出球锥及类升力体外形的俯仰阻尼导数计算结果。     

涡轮级内湍流度影响的非定常数值模拟

通过求解二维雷诺平均Navier-Stokes方程，数值模拟了不同进口湍流度下的涡轮级中的非定常流场及温度场。结果表明，增大湍流度对叶栅通道内时均压力分布影响很小，但是对涡轮级中的温度分布会产生很大影响，当涡轮进口有热斑存在时，增大湍流度将会使热斑中的热流体温度扩散更快，从而叶栅通道中的温度场分布趋向均匀，同时动叶压力面侧与吸力面侧的时均温度差别减小，沿动叶统长方向上的温度分布也更平缓。     

叶轮机内部转/静干涉作用的数值模拟

本文对叶轮机械内部转/静干涉问题的非定常数值模拟进行了研究.采用双重时间步方法和Dornain-Scaling方法求解二维N-S方程组,数值模拟了某涡轮叶栅内部的非定常流场.通过计算展示了在转/静干涉作用下的叶型尾迹,转子尾迹的传输和尾迹的切割等非定常现象.     

基于机器学习的非定常流场网格自适应

现有针对非定常流场数值模拟的网格自适应方法,通常每隔一段时间步就进行一次网格调整,增加了计算复杂度和精度损失。针对这一问题,本文基于间断伽辽金（discontinuity Galerkin, DG）有限元法提出了结合BPNN(backpropagation neural network)和MMPDE(moving mesh partial differential equation)的非定常流场网格自适应方法。该方法首先采用DG有限元法对Navier-Stokes方程进行非定常计算,得到统计意义上的网格间断量;然后以初始网格和间断量训练BPNN回归模型,用于预测任意位置节点的间断量;接着使用MMPDE变分法移动网格节点,使其符合统计意义的间断量分布;最终通过Laplacian网格平滑法保证网格单元质量。圆柱绕流非定常流场算例的验证结果表明,该方法能够在不改变网格拓扑结构和不增加节点数的情况下完成一次性网格自适应,显著提高了非定常流场数值计算的精度和效率。     

COMPUTATION OF COMPRESSIBLE FLOW PAST SLENDER WING－BODY USING FULL－POTENTIAL EQUATION

ＣＯＭＰＵＴＡＴＩＯＮＯＦＣＯＭＰＲＥＳＳＩＢＬＥＦＬＯＷＰＡＳＴＳＬＥＮＤＥＲＷＩＮＧ－ＢＯＤＹＵＳＩＮＧＦＵＬＬ－ＰＯＴＥＮＴＩＡＬＥＱＵＡＴＩＯＮＨｕａｎｇＭｉｎｇｋｅ（ＤｅｐａｒｔｍｅｎｔｏｆＡｅｒｏｄｙｎａｍｉｃｓ，ＮａｎｊｉｎｇＵｎｉｖｅ...     

用全位势方程计算细长翼-身组合体的可压缩绕流

用Ｈ－Ｏ网格、守恒型全位势方程、差分和隐式近似因式分解迭代算法计算绕大后掠细长翼－身组合体的可压缩流动。对ＡＧＡＲＤ－Ｂ翼－身组合体的计算表明,计算的自由流Ｍａｃｈ数可从亚音速直到低超音速。算出的升力与力矩特性与实验数据符合较好     

不同轴长比椭球体的绕流计算与分离形态分析

用拟压缩性方法求解不可压Ｎ－Ｓ方程,对４∶１旋成体椭球和其他轴长比椭球的层流绕流进行数值模拟。分别用二阶中心差分格式和三阶迎风紧致差分格式来逼近无粘项。从计算所得的表面摩擦应力线、空间流线和横截面流场特性,分析了１０°≤α≤７０°迎角范围内的分离和脱体涡系形态,算出了相应的气动特性。计算结果表明,三阶迎风紧致差分比二阶中心差分可更好地模拟分离涡结构。迎角直到４０°时,算出的流态仍然对称,而迎角为７０°时,算出的流态为非对称非定常,与实验流态基本相符。椭球体横截面形状越扁平,开式分离就越提前产生。椭球体横截面形状横向越窄,分离区就越小。由于改变椭球体横截面形状具有改变分离和脱体涡系的效果,本文的研究结果有助于改善三维物体的气动力特性,为改进机体设计提供理论依据     

翼型和叶栅前后缘局部流场的新精化解法

本文提出一种利用前后驻点自动定位来求前后缘局部流场精化解的方法 ,即在全流场解的基础上 ,把局部流场解分为两部分的叠加 ,其中之一是驻点邻区的近似解析解 ,它吸收了原解中变化急剧的部分 ,使剩余部分成为变化比较平缓、具有更好的解析性态的分布 ,因而易于求准确。算例表明 ,本文方法为提高亚声绕统计算精度 ,特别是带攻角绕流工况 ,提供了有效的手段 ,本文思想也可同变分差分法结合应用     

隐式无网格算法及其应用研究

本文的主要目的在于研究求解Euler方程的隐式无网格算法，并应用于复杂的流场计算。采用无网格算法，计算区域用点云离散代替通常的网格划分；计算点上的空间导数，用当地点云上引入的二次极小曲面逼近。求解的Euler方程用隐式时间后差离散，结合用Roe的近似Riemann解确定通量，并用LU－SGS算法分步计算，数值求解了单翼型或双翼型模拟的复杂绕流。     

后缘喷流对三角翼绕流影响的N-S方程数值分析

本文用拟压缩性方法求解不可压流雷诺平均拟压缩Ｎ－Ｓ方程组,对带有后缘喷流的三角翼粘性绕流进行了数值模拟,求解中采用了Ｂｅａｍ－Ｗａｒｍｉｎｇ隐式近似因子分解格式以及ＭＭＬ代数湍流模型。计算结果说明,后缘喷流使涡核压强降低,使涡核速度增大,从而对三角翼前缘分离涡有稳定作用,并能增大上翼面的负压值和下翼面的正压值,从而可以增加部分升力。计算结果还说明,喷口面积或喷流下偏会使上述作用增强。     

二维埋入式进气道的数值分析

采用并改进了适用于分离流动的 Ｊｏｈｎｓｏｎ Ｋｉｎｇ 紊流模型,结合 Ｌ ＵＳ Ｇ Ｓ近似隐式格式和基于 Ｍ Ｕ Ｓ Ｃ Ｌ非线性插值的 Ｒｏｅ 格式编写了雷诺平均的二维 Ｎ-Ｓ方程的求解程序。然后借助于该流场计算程序完成了二维埋入式进气道的数值研究,分析了几何设计参数对进气道性能的影响。     

紧致格式数值模拟超音速粘性绕流问题

通过数值求解可压Navier-Stokes方程的方法计算了球锥粘性绕流问题。方程中粘性项按通常办法逼近,为改善收敛速度利用了文[1,2]中之算子附加修正方法。方程中之流向导数利用了四阶精度的紧致差分。这是由于该方向上流动参量光滑和边界条件易于处理的缘故。在法向方向上利用了中心差分。对差分方程之隐式部分利用了近似因式分解法。由于选用了文[3]中之特殊Jacobiall系数矩阵分裂法而使计算工作量大为减少。用这种方法计算了超音速粘性绕流问题。计算结果与实验结果相符很好。     

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。     

基于近似熵谱的碳纤维复合材料层压板拉伸损伤声发射分析

采用近似熵理论分析了碳纤维环氧树脂层压板在拉伸过程中产生的三种典型声发射信号。研究 表明基体开裂、脱胶以及纤维断裂这三种信号的近似熵谱特征不同,结合经验模式分解方法解释了三种损伤破 坏机制。此外,在经验模式分解的基础上考察了三种声发射信号本征模函数的频谱分布,进一步说明了近似熵 在本质上是衡量信号发生新振动模式产生的概率以及在复合材料声发射检测与损伤识别中使用近似熵谱分析 的有效性。