牛顿流理论及其在谐振机翼上的应用

根据牛顿流理论,文中讨论了高马赫数、任意迎角、一般剖面在均匀来流中作沉浮和俯仰谐振的翼型和机翼。导得翼型和机翼表面的牛顿撞击压强和安定性导数公式。给出楔、平板三角翼、双抛物线剖面三角翼俯仰刚性导数(-C_(mθ)和阻尼导数(-C_(m))的数值例子。数值结果与现有的牛顿-布兹曼流理论、薄激波层理论、欧拉方程摄动理论的结果作了比较,算例表明本文的结果令人满意。     

超声速和高超声速机翼俯仰导数—当地流活塞理论解法

本文给出计算超声速和高超声速尖前缘三维机翼俯仰导数的一个近似解析方法。该方法是根据当地流活塞理论和片条理论导得的。在激波附着于机翼前缘情况下,本方法可用于计算任意迎角,颇为一般的平面形状、尖前缘剖面机翼的俯仰导数。中小平均迎角的数值结果与有限的试验数据和其它理论作了比较,示例表明本方法的结果令人满意。     

非定常超音速侧滑机翼的升力、力矩系数及动导数的数值计算

 本文利用分布在机翼平面上的非定常超音速点源的特征线网格法给出非定常超音速侧滑机翼的数值解。本文给出侧滑机翼作沉浮、俯仰、滚转谐振时的振型以及计算升力、力矩和动导数（包括交叉导数）的计算方法。按本文用非定常理论数值解对典型的无侧滑机翼计算出的动导数与常用的按准定常理论解析公式计算结果进行比较和鉴定,结果表明本文方法是令人满意的。     

超,高超声速机翼俯仰导数的欧拉方程摄动解法

应用片条理论，将文献（１）的尖前缘任意翼型的俯仰稳定性导数的计算方法作、推广、用于计算任意平面形状的三维机翼。在各种特殊的情况下，其计算结果与现有的理论都符合得很好。在激波附体的条件下，本文的计算方法可用于计算任意迎角。任意平面形状机翼的超、高速声速俯仰稳定性导数。     

计算高亚声速情况下三维机翼动稳定性导数的面元法

本文应用我们提出的跨-亚音速三维谐振机翼的局部线化面元法,同时采用了与之配偶的亚临界定常流绕三维厚度机翼的面元法迭代解,对几种不同平面形状的机翼,进行了俯仰和滚转阻尼导数的计算和分析,并和现有的实验数据以及线化理论的计算结果进行了比较,从而表明,本方法可以得到满意的动导数计算结果;还表明,当接近临界马赫数时,动导数计算必须计及跨音速非线性效应。     

近高超声速流中旋成体俯仰阻尼导数的计算

本文用改进的二阶激波-膨胀法导出了旋成体法向力静导数的二阶Eggers公式,并计算了俯仰尼阻导数。它与修正的Newton理论相比,由于在压力分布的计算中考虑了流线上前方气动参数的影响,理论上更加完善。计算结果与实验数据比较表明,这是一种比较有效的工程计算方法。     

计算机翼跨声速绕流的全速位积分方程数值解

本文利用格林定理将全速位方程转换为全速位积分方程，应用高斯定理对沿激波面的显式积分进行简化。将全速位积分方程的解表示成为面元项和场元项之和，然后离散进行数值求解。本文将激波捕捉方法和激波装配技术结合起来，对机翼跨声速绕流进行数值实验，计算结果令人满意。     

二维非定常欧拉方程的摄动解在超声速矩形翼上的应用

本文利用超声速线化理论给出的矩形翼载荷系数解析解作为权函数,对沉浮和俯仰谐振的二维平板翼压强系数进行加权修正,获得矩形翼谐振气动力的近似解,计算结果与线化理论数值结果符合得很好,本方法是超声速和高超声速矩形翼非定常气动力的一种统一和快速的计算方法。     

升力浮力复合型飞艇动导数分析

计算了球体和椭球体的平移、转动加速度导数,数值模拟了NACA0015翼型的俯仰振荡,验证了计算流体动力学(CFD)方法计算动导数的有效性.模拟升力浮力复合型飞艇的强迫振荡运动,得到了全艇以及机翼、艇身、平尾和垂尾的动导数,并与常规飞艇进行了比较分析.机翼对速度导数的贡献导致复合型飞艇升沉速度导数、滚转角速度导数较常规飞艇显著增大.常规飞艇加速度导数以惯性力成分为主,但黏性力成分以及机翼对艇身和尾翼的干扰对复合型飞艇升沉、滚转和俯仰3个方向的加速度导数有重要影响,不能将加速度导数简化为附加质量系数.      

弹性机翼的跨音速压强分布计算

本文采用迭代求解的方法计算弹性机翼的跨音速压强分布。用跨音速小扰动势流的二级近似方法求解气动力;结构变形则采用一维简单梁理论或二维平面结构的柔度影响系数矩阵进行计算。对M6弹性机翼所作的示范性计算表明,计算弹性翼的迭代次数与刚性翼计算相当,弹性变形对冀面气动载荷分布和激波强度均有明显影响。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

Experimental investigation of shock-buffet criteria on a pitching airfoil

An experimental investigation of the shock-buffet phenomenon subject to unsteady pitching supercritical airfoil around its quarter chord has been conducted in a transonic wind tunnel. The model was equipped with pressure taps connected to the fast response pressuretransducers. Measurements were conducted at different free-stream Mach number from 0.61 to0.76. The principle goal of this investigation was to experimentally discuss the shock-buffet criterion over a SC(2)-0410 supercritical pitching ...     

低雷诺数下翼面局部振动增升机理研究

采用计算流体力学(CFD)方法研究低雷诺数下翼面局部振动对翼型气动特性及其流动特征的影响规律。建立局部振动激励的力学模型,并采用任意拉格朗日-欧拉坐标系下的特征线有限元(ALE-CBS)方法对局部振动激励下翼型绕流问题进行模拟,分析局部振动对非定常流动演化的影响规律,揭示其增升机理。研究结果表明:翼面局部变形的增加会有效降低翼型上表面的前缘压力;非定常流动分离中旋涡之间的距离及其演化频率与振动频率的关系是影响翼型翼面局部振动增升效果的重要因素。当流场主频率与振动频率相同,次要频率为主频率的2倍,即发生锁频时,翼面振动产生的移动分离泡能够使分离区从主流获取更多的能量,使翼型上表面保持较低的压力,有效提高翼型升力。     

用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

带连续变弯度后缘操纵面机翼的动态失速减缓

主动偏转后缘操纵面可以减小翼型动态失速对气动特性产生的不利影响。研究连续变弯度后缘操纵面在减缓翼型动态失速方面的性能,利用CFD 结合动网格方法,计算NACA 0012 翼型在大幅度俯仰振荡时的非定常气动力;从减缓效果和能量需求两个方面,对比传统刚性操纵面和两种连续变弯度操纵面的动态失速减缓性能。结果表明：当后缘操纵面按正弦脉冲规律偏转时,可以推迟前缘涡的产生,加速后缘涡的发展,降低压力分布在后缘的峰值,进而减小动态失速时翼型的低头力矩极值;后缘操纵面的弯度构型会影响减缓效果,在相同的偏转策略下,弯度描述函数为2 阶多项式的连续变弯度操纵面的减缓效果最好,且能量需求最小。     

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

俯仰振荡翼型推阻力转变滞后机制数值研究

扑翼产生的反卡门涡街被认为是一种推力型尾迹,但已有研究指出,随着斯特劳哈尔数(St)增大,低雷诺数下俯仰振荡翼型的净推力产生明显滞后于反卡门涡街的出现。为探究该现象背后的物理机制,对NACA0012翼型在雷诺数1000条件下作简谐俯仰运动的流场进行了数值模拟。采用翼型表面积分方法和基于有限控制体的气动力估计方法分别研究了翼面分布力特性和尾迹流场特性变化对阻力-推力转变临界点的影响。翼面分布力积分结果表明,当翼型振荡参数进入对应反卡门涡街尾迹形态区域时,翼面压力分布产生的推力分量无法克服剪切层的黏滞阻力,是造成俯仰翼型推力产生滞后于反卡门涡街出现的主要原因。对尾迹流场及相应的推阻力特性变化的分析表明,尽管反卡门涡街会产生“喷流效应”,但在St较小时,其产生的动量诱导推力无法克服反卡门涡自身的涡致阻力和尾迹流场压力损失引入的附加阻力,因此即使存在反卡门涡街也不能产生净推力,进而从流场分析的角度进一步解释了这一滞后现象的发生机制。     

翼型舵面偏转非定常流动数值模拟

针对NACA0012翼型舵面偏转问题,数值模拟了不同参数对翼型气动特性的影响。基于非结构动网格技术,采用ALE有限体积描述下的二维可压缩非定常N-S方程,计算通量采用Vanleer格式、时空二阶格式,利用Venkatakrishnan限制器抑制数值振荡。非定常计算结果表明,NACA0012翼型绕1/4弦点作周期性俯仰振动的升力系数和俯仰力矩系数结果与实验数据吻合良好,验证了数值方法的准确性;在翼型舵面表面有分离区产生,升力系数和俯仰力矩系数形成滞回环,在亚声速情况下,滞回环幅值较小,进入超声速阶段以后,幅值增大,随着翼型间缝隙宽度逐渐增加,翼型升力系数和俯仰力矩系数与无缝翼型相比逐渐降低。     

多段翼型局部主动变形流动控制的非定常数值模拟

对30P30N三段翼型失速攻角附近的分离流动进行了数值模拟研究。为了抑制大攻角时背风区的流动分离,在主翼段上表面引入了行波壁变形模型和抛物型局部主动振动模型,利用作者以往发展的动态混合网格技术和相应的非定常计算方法,对变形过程中的非定常分离流动进行了数值模拟,分析了各种变形参数对流动分离的影响。计算结果表明,在适当的条件下,局部主动变形能够抑制翼型背风区的分离,由此可以起到增升减阻的作用,改善翼型的气动性能。