阵风/叶栅干涉噪声的一种前馈控制策略研究

 基于反问题的研究思路,探讨了阵风/ 叶栅干涉噪声的一种前馈式主动控制策略。具体研究思路是试图通过在静子叶栅表面引入附加的“二次级声源”,精确地消去所感兴趣区域的声音,二次级声源强度基于阵风/ 叶栅干涉气动声学反问题模型获得,并通过数值实验的方法探寻了二次级声源合理布放策略。研究结果不仅证实了在叶栅表面布放偶极子声源的可行性,而且提出了在尽量不影响气动性能前提下的前馈主动噪声控制方法。     

基于浸入式边界法的叶栅颤振数值模拟

基于浸入式边界法建立了求解振荡叶栅非定常流动的快速计算模型,并在模型中添加湍流模型,使之能够考虑更加接近实际流动的情况。具体的方法是求解目前为止最常用的雷诺时均Navier-Stokes方程并添加低雷诺数湍流模型(Lam-Bremhorst模型,k-ε模型的一种)进行计算,为了验证该方法的正确性,分别对层流边界层和湍流边界层进行了数值模拟,计算结果与布拉修斯解和壁面率吻合得很好,证明了湍流模型的可靠性。在此基础上,对高雷诺数条件下的振荡叶栅进行了数值模拟。结果表明,折合速度是影响叶栅振荡的重要因素,这与层流计算的结论类似。值得注意的是耦合过程没有生成贴体网格,减小了计算时间,可以准确快速地模拟真实的叶栅流动情况。     

基于浸入式边界法的振荡转子叶片数值模拟

对于流固耦合问题,传统的数值模拟方法由于需要不断的网格重构使得问题变得复杂,本文通过利用浸入式边界法建立了在单一坐标系和网格下求解振荡转子叶片的快速计算模型,从而避免了传统方法中由于需要不断地重构贴体网格造成的数值模拟的复杂性。为了验证该方法的正确性,分别对低KC(Keulegan-Carpenter)数下的振荡圆柱和两自由度振荡圆柱进行了数值模拟,计算结果与以往的试验结果和数值结果吻合得很好,证明了浸入式边界方法的可靠性。在此基础上,对振荡转子叶片进行了数值模拟。结果表明折合速度和叶栅稠度是影响振荡转子叶片的重要因素。并且,值得注意的是这种耦合过程并没有生成贴体网格,因此大大减小了计算时间,这样可以更加快速、准确地模拟真实的叶栅流动情况。     

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

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

错频对叶片的气动弹性稳定性影响

采用求解动网格下的非定常可压缩Navier-Stokes方程模拟了振荡叶片下的气动弹性问题,研究了叶片的频率错频对叶片气动弹性稳定性的影响.通过数值模拟平面叶栅的气动弹性第10标准算例,验证了气动弹性的数值模拟方法,计算了不同叶片间相位角和折合频率下的气动阻尼系数,研究了叶片振动频率改变对叶片气动弹性稳定性的影响.计算结果说明:频率错频是提高气动弹性稳定性的有效方法,其主要作用是减小叶片间振动的耦合效应和叶片间相位角的影响,并且随着错频量增大叶片稳定性增强;通过模拟三维环形叶栅的气动弹性第4标准算例,计算了气动阻尼系数随错频量变化的规律,验证了错频量和气动弹性稳定性增强的规律.      

组合压气机中尾迹与势流同频/异频干涉的叠加特性

为揭示多级压气机中上下游叶轮对中间叶片叠加气动影响特性,阐述不同叠加干涉情况下下游叶轮进气角度变化,采用数值方法模拟了一级轴流和一级离心组成的组合压气机非定常流场.详细讨论了上游动叶尾迹和下游动叶势流对中间导流叶栅段气流非定常流动的异频和同频叠加干涉特性,依据计算结果,直观地展示了静叶通道中两种干涉间相互激励和抑制作用的位置和时间,与数学公式的推导结果相互印证.研究结果表明:当上下游动叶对中间静叶段异频干涉时,干涉的激励、抑制区域的轴向位置随时间发生变化;当上下游动叶对中间静叶干涉频率相同时,干涉的相互激励、抑制区域的轴向位置不随时间发生变化,但干涉的激励、抑制区域的轴向位置受时序位置影响.另外,上游动叶尾迹与下游离心叶轮势流的不同叠加情况,决定着下游离心叶轮进口相对气流角的大小及波动幅值.     

尾缘冷气喷射对超声涡轮叶栅性能的影响

采用试验与数值模拟相结合的方法,研究了某超声速涡轮导向叶栅尾缘冷气喷射对叶栅流场结构的影响。数值模拟时,使用环形叶栅模型近似模拟平面叶栅内的流动。研究结果表明:数值模拟结果与试验结果吻合较好;尾缘冷气喷射可减少主气流在尾缘停滞区的能量耗损,削弱叶栅尾缘处的内边缘激波,叶栅气动效率随冷气量的增加先增大后减小;尾缘冷气喷射对叶栅出口附近气流角的周向分布有影响,但对质量平均的叶栅出口气流角基本无影响。     

翼型和机翼阵风响应的数值模拟

采用有限体积法来数值求解非定常Euler和NS方程,并在非定常计算中通过引入"网格速度"来计入阵风的影响,实现了对刚性翼型和机翼阵风响应的有效数值模拟。首先采用该方法来对翼型攻角突然增大的阵风、陡边型阵风这两种典型阵风情况下的翼型气动力响应进行了计算,结果表明:计算结果与现有的一些理论结果都符合较好,与文献中计算结果也吻合,NS方程的计算结果与Euler方程计算结果在变化趋势上是一致的,但具体的数值有差别。最后针对攻角突然增大的阵风,对不同马赫数下不同展弦比平直机翼的阵风响应进行了数值模拟,得出的计算结果与文献中的计算结果一致。     

阵风发生装置流场测量与分析

以小型阵风发生装置为研究对象,采用二维热线风速仪,测量了不同工况下的阵风装置流场参数,给出了二维热线仪测量方案和在线角度修正方法,以及阵风幅值和波形稳定性计算方法。结果表明：采用风洞在线角度修正,可以提高二维热线的y向测量精度和测试效率;采用波形相位分析方法,可以满足阵风发生器产生的周期性波形的幅值和流场脉动量分析;阵风幅值与4～15Hz的叶栅摆动频率、8°～30°的叶栅摆动角度、50～100mm的叶栅弦长成正比关系;阵风波形的不稳定性（等相位速度脉动量）与叶栅摆动角度、来流速度有一定的线性关系,在本试验范围内随叶栅摆动角度、来流速度的增加而增加,同时,适当的叶栅摆动频率可以降低叶栅的波形不稳定性。研究结论对阵风发生装置研制及其流场测量具有一定的参考意义。     

涡轮叶栅叶型损失的数值模拟方法

给出了一个计算亚、跨音涡轮叶栅叶型损失的数值计算方法。主流采用时间推进有限体积法求解积分型欧拉方程,并采用了局部网格修正技术;附面层采用全隐格式求解有限差分方程;在叶栅出口与远后方均匀流之间进行了叶片尾迹与主流的掺混损失计算。算例表明本文的数值方法可准确地预测涡轮叶栅的叶型损失。     

激波／边界层干扰数值模拟的格式效应

用不同的数值格式,包括Ｒｏｅ的ＦＤＳ格式和ｖａｎＬｅｅｒ的ＦＶＳ格式并配备不同限制器,求解二维Ｎ－Ｓ方程,模拟研究了激波／边界层干扰的流动问题,分析了模拟出的流动现象,仔细研究了其中的格式效应,得出同前人不尽一致的结论：ｖａｎＬｅｅｒ格式对限制器的依赖性比Ｒｏｅ格式更强;在本文情况下,ｍｉｎｍｏｄｌｉｍｉｔｅｒ的激波、膨胀波的分辨率比ａｌｂａｄａｌｉｍｉｔｅｒ还好     

电弧放电等离子体激励控制超声速压气机叶栅激波/边界层干扰仿真研究

为研究电弧放电等离子体激励对超声速压气机叶栅激波/边界层干扰的控制作用,建立了模拟等离子体激励作用效果的唯象学模型,进一步以ARL-SL19超声速叶栅为对象,通过数值仿真研究了电弧放电等离子体与叶栅通道内部流动的相互作用及其对叶栅流动损失的影响。结果表明:等离子体唯象学模型能够较好模拟电弧放电等离子体诱导产生冲击波的气动特性。电弧放电等离子体激励对叶栅通道内部流动主要具有三种作用效果:在放电区,注入的热量会产生阻塞效应,增加近壁面气流的流动损失;在激波/边界层相互作用区,能够改变激波系结构,减小激波损失;在尾迹区,冲击波会诱导产生脱落涡。     

不同壁温及差分格式下超燃冲压发动机的仿真

针对圆截面超燃冲压发动机相关实验进行数值仿真,采用了不同的壁面温度条件和差分格式,并探究其对仿真结果的影响规律。结果表明采用常用的绝热壁面边界条件会导致激波边界层交互位置处的回流区变大,进而导致流道收缩,壁面压力升高。采用高精度格式不仅不能更为精确地计算回流区的大小,反而会使回流区面积变得更大。通过改变壁面边界条件,采用300 K的等温壁面,可以使激波边界层交互位置处的回流区变得极小,使仿真结果与实验结果十分吻合。计算结果表明,针对当前构型与实验条件,盲目提高数值格式的精度并不一定会得到更好的结果,相反可能会使仿真与实验更加偏离,而适当修改壁面温度条件,即使数值格式精度较低,仍然可以得到很好的仿真结果。最后,针对与实验结果吻合的仿真结果,分析了圆截面超燃冲压发动机的流场特性,重点研究了凹腔剪切层及其质量交换特性。     

旋翼BVI噪声的理论模拟与分析

ＢＶＩ噪声是旋翼产生的一种非常强烈的脉冲式噪声,本文应用理论模拟的方法分析了ＢＶＩ噪声的规律。对叶片表面非定常力的计算采用了有限翼展叶片对倾斜正弦式阵风的响应函数,声场计算应用了基于Ｌｉｇｈｔｈｉｌｌ声类比理论的频域方法。应用本文方法计算的ＢＶＩ噪声声压信号与实验结果比较吻合较好。     

欠膨胀火箭多股喷流干扰流场的数值模拟

本文利用高精度的TVD有限体积差分格式,对欠膨胀火箭多股自由喷流的干扰流场进行了数值研究,给出了两股自由射流的干扰流场的复杂激波系结构,包括内、外层反射激波、相交激波等波系结构。数值结果与实验结果较为吻合。     

一种新的边界处理方法在笛卡尔网格中的应用

采用自适应笛卡尔网格方法求解Euler方程,固壁边界条件通过一种离散强迫浸入边界方法引入,边界切割网格与流场网格保持一致,从而有效解决了小切割网格单元的时间步长限制问题.针对自适应笛卡尔网格的特点提出了新的虚拟点反射方法和插值方法,解决了数值吸力峰值问题.对RAE2822翼型、双NACA0012翼型和Suddhoo三元翼型的流动状况进行了无黏数值模拟,并与现有的理论解、单域以及分区结构网格解进行了对比.结果表明:该离散强迫浸入边界方法有效地引进了固壁边界条件,结合自适应笛卡尔网格技术能够准确模拟复杂几何外形下的流动状况.      

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

跨声速压气机数值模拟格式效应研究

针对跨声速叶轮机械复杂内部粘性流动,开发了全三维高精度RANS数值求解器。分别采用FVS格式与AUSM+格式结合多种限制器对Rotor 67跨声速压气机转子进行数值模拟,分析和比较了各格式的计算效果并与试验结果对比。结果表明:与FVS格式相比,AUSM+格式的数值粘性更小,边界层的模拟精度更高;Hemker限制器的综合表现最优,Van Albada限制器的粘性分辨率略低于Hemker限制器,Minmod限制器对流动分离现象的捕捉能力较差,Van Leer限制器容易引入色散误差。     

Time-domain nonlinear aeroelastic analysis and wind tunnel test of a flexible wing using strain-based beam formulation

A theoretical formulation for time-domain nonlinear aeroelastic analysis of a flexible wing model is presented and validated by wind tunnel tests. A strain-based beam model for nonlinear structural analysis is combined with the Unsteady Vortex Lattice Method (UVLM) to form the complete framework for aeroelastic analysis. The nonlinear second-order differential equations are solved by an implicit time integration scheme that incorporates a Newton-Raphson sub-iteration technique. An advanced fiber optic sensing technique is firstly used in a wind tunnel for measuring large structural deformations. In the theoretical study, the nonlinear flutter boundary is determined by analyzing the transient response about the nonlinear static equilibrium with a series of flow velocities. The gust responses of the wing model at various gust frequencies are also studied. Comparisons of the theoretical and experimental results show that the proposed method is suitable for determining the nonlinear flutter boundary and simulating the gust response of flexible wings in the time domain.