压气机二维叶栅涡脱落的数值模拟

对某压气机二维叶栅的非定常分离流场进行了数值模拟,通过对多种工况的计算,进行了频谱分析,对叶栅非定常流动的流场结构和流动机理做了初步的探讨。在来流均匀,定常边界条件下,叶栅内流动仍然表现出强烈的非定常性;在大攻角下有类卡门涡的脱落;旋涡脱落的频率,随着攻角和马赫数的变化而变化;分离点的位置随着攻角和马赫数的变化而变化。     

三维扩压叶栅非定常流动机理研究的频谱分析

计算了三维直叶栅在不同攻角、不同马赫数下的流动情况,得到流场的非定常解,并进行了频谱分析,对叶栅非定常流动的流场结构和流动机理做了初步的探讨。分析计算结果表明:在来流均匀,定常边界条件下,叶栅内流动仍然表现出强烈的非定常性。分离区和尾迹中的流动,以旋涡的有规律周期性脱落为主要的运动形式。旋涡脱落的频率,随着攻角和马赫数的变化而变化:同马赫数下,攻角越大,频率越低;同攻角下,马赫数越高,频率越高。同时,在同一工况下,旋涡频率沿叶高呈非均匀分布,叶中区域频率相对低,靠近端壁区频率相对高。      

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

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

基于当地流活塞理论的翼-身组合体飞行器大攻角超声速颤振分析

基于当地流活塞理论推导了旋成体机身和任意外形三维机翼非定常气动力计算公式,结合CFD数值仿真,将当地流活塞理论推广用于复杂外形飞行器全机大攻角超声速颤振计算。以翼-身组合体面对称外形飞行器为例,进行带攻角的超声速颤振计算,并与非定常气动力采用CFD全数值仿真求解得到的颤振计算结果比较,验证了本文方法的正确性和效率。     

计算超音速非定常气动力的一种新方法及其应用

本文提出了一种计算超音速非定常气动力的新方法,它将超音速非定常的活塞理论和定常的锥形流理论结合起来,使得被活塞理论所忽略的翼面上各点间的相互影响用锥形流理论来近似计入。文中运用该方法对十个不同平面形状的机翼进行了二十九种状态的颤振计算,并和风洞颤振实验结果作了比较。比较的结果是令人满意的。特别是和用活塞理论所做的颤振分析相比,计算精度有了明显提高。这表明本方法用于超音速颤振分析是确实可行的,它具有方法简洁,计算精度高,计算时间少,程序编制容易的优点。     

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

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

带操纵面机翼的非定常跨声速流数值解法

本文介绍一种带操纵面机翼的非定常跨声速流的有限差分计算方法。采用的方程是三元非定常跨声速修正小扰动位势方程,使用时间积分法,求解格式是近似因式分解交替方向隐式(ADI)格式。使用这种方法计算了F-5机翼在来流马赫数为0.9和0.925时的定常气动力和操纵面振荡的非定常气动力。计算结果与国外的NLR试验结果和XTRAN3S方法的计算结果进行了比较,表明计算是成功的。     

超音速、高超音速飞行器动导数的高效计算方法

将超音速、高超音速气动力工程计算方法综合起来,推广到非定常气动力的计算,并用于求解任意外形飞行器作强迫运动时的非定常气动力.在此基础上进一步发展了一种适合于超音速、高超音速大迎角状态的飞行器动导数计算方法和程序.计算结果与国际标准实验数据作了比较,证明了该程序的正确性和精度.研究了减缩频率和基本迎角对动导数的影响,结果表明基本迎角对动稳定性导数影响较大,而减缩频率在小迎角范围内对动导数影响很小.     

可变形翼型超声速非定常气动力的研究

由小扰动线化方程出发,利用基于格林公式的边界元方法求解二维超声速流动速度场,进而对变形翼型气动特性进行理论分析和数值计算。计算结果显示虚拟质量力导致的非定常附加升力对厚度变化远比弦长或弯度变化敏感;弦长或弯度变化诱导的非定常效应影响很小,此时非定常气动升力和准定常气动升力近似相等,和不可压缩或亚声速时相类似,变形翼型的非定常气动升力近似等于准定常计算结果叠加上非定常附加升力,与具体的变形历史过程无关。同时还分析了超声速时翼型往复变形时的气动特性,和亚声速情况完全不同,随着马赫数升高,完成往复变形所需要外界输入的功逐渐减小,趋近一个常值。另外翼型完成往复变形所需做的功将不依赖于来流攻角,且与翼型参数的变化加速度a存在简单的数学关系:A ∝a～(1/2)。     

新型高超音速静压探针的设计和实验

参考了近40年来国外有关超音速和高超音速静压探针的研究成果,设计了多种并发展了二种高超音速静压探针。它们可精确地用于持续工作流场仅为30ms的流场静压测量。实测静压与自由流静压比p_(s,m)/p_(s,∞)=1.009;马赫数比M_(∞m)/M_(∞,)=0.993。实验在加拿大多伦多大学宇航研究院的高超音速炮风洞里完成。炮风洞喷管出口气流特性是:(空气)自由流总压比_(t,∞)=25.5MPa,总温T_(t,∞)=1000K;马赫数M_∞=8.30;自由流静压p_(s,∞)=2.13kPa;雷诺数R_e=3.2×10~7。以探针直径D为基准的R_(e,D)=34,000～52,800。这种探针己被用于“管内高超音速流场发展实验研究”中了。     

超音速、高超音速有攻角的三维薄翼非定常二次理论及其应用

 动导数对于飞行器的动态品质分析和控制系统的设计是必不可少的参数,它直接影响飞行器的动态特性。因此,国外自七十年代以来,无论在风洞实验、飞行试验,还是分析、计算等方面都作了很多工作,取得了较大的进展。 动导数计算说到底主要是非定常气动力计算问题,已有的线化理论不能满足需要,对考虑非线性效应的动导数计算方法的需求则日显迫切。二次理论是计算非线性效应的一种方法,它是由第二次叠代或物面偏角参数的二次展开法构成,主要非线性效应由二次方程算出。     

跨音速压气机动静叶排相干的三维非定常流动数值分析

采用时空二阶精度ＮＮＤ显式差分格式求解三维Ｅｕｌｅｒ方程，对跨音速压气机级中动静叶排相干形成的非定常流动进行了数值分析。除动静叶排流动的控制方程体系、数值方法外，着重讨论了动静叶排计算域、网格以及边界条件的处理。对某压气机第１级动静叶排相干的非定常流场的计算，给出了叶片表面的马赫数分布、压力波动幅值，以及叶片槽道中瞬态密度等值线等，并分析了结果的合理性。     

低超声速来流中机翼跨声速非定常气动力数值解法

本文给出一种低超声速来流中机翼跨声速非定常流的计算方法。该方法使用跨声速小扰动理论,正确地计入了沿流向平面流动特征线的压力波的运动边界条件。流动方程是使用时间精确交替方向隐式(ADI)有限差分解法求解。试算了M=1.1的矩形翼和M=1.1、1.34的F5机翼的定常和非定常气动力,并与国外的风洞试验和计算结果作了比较,符合很好。表明本方法是有效的。     

二维非定常跨音速时间线化积分方程的推导

1.引言 求解非定常跨音速流动的主要困难是非线性问题。对于微幅翼型振动问题(同时引起激波的微幅振动)可做时间线化简化处理,得到时间线化微分方程定解问题。时间线化积分方程是由时间线化微分方程推导出来的。     

一种基于动模态分解的翼型流动转捩预测新方法

考虑自由转捩的定常/非定常流动Navier-Stokes方程数值求解,对于翼型流动细节的精确模拟和气动力的精确预测均具有十分重要的意义。采用动模态分解(DMD)方法进行流动稳定性分析,再结合e N方法,提出了一套适用于翼型绕流的转捩预测新方法,称为DMD/e N方法。相比于传统的线性稳定性分析方法,DMD方法不需要求解附面层方程和线性稳定性方程,也没有引入平行流假设,具有更好的理论适用性和算法鲁棒性。开展了NLF0416、S809和SD7003等翼型的转捩预测数值验证研究,通过与实验结果以及与传统的基于线性稳定性分析的e N方法的比较,验证了本文所发展的转捩预测新方法在预测翼型的定常流动和非定常流动转捩方面的正确性,也表明了该方法具有解决含层流分离泡的翼型绕流转捩预测的能力。     

二维机翼二自由度跨音速颤振分析

 <正> On the assumption that airfoils take small amplitude ha-rmonic oscillating and middle reduced frequency, the two-dimensional, inviscid, small-disturbance nonlinear unsteady transonic flow equation can be split into two parts. One is the nonlinear steady small-disturbance equation, the other is the time-linearized unsteady small-disturbance equation. The steady equation is solved by Carlson code. The unsteady equation is solved by integral method. The standard v-g method is used to solve the flutter eigenvalue equations. A transonic flutter analysis is performed for a NACA 64A006 airfoil with pitching and plunging degrees of freedom.The Mach numbers considered are 0.7, 0.8, and 0.85. The aerodynamic coefficients are obtained by the above method. For each Mach number, the flutter speed and the corresponding flutter reduced frequencies are achieved by varying the airfoil-airmass ratio, plunge-pitch frequency ratio. The results are compared with other computational method and good agreement has been observe     

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

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

预处理法求解定常/非定常混合网格的全速流场

运用时间导数预处理法在混合网格上求解三维定常/非定常N-S方程。预处理后的N-S方程,用二阶迎风格式的有限体积法离散,隐式、LU-SGS迭代求解,双时间步推进求解非定常流。本文以一维预处理欧拉方程为基础进行分析,提出了一种特征变量的边界条件以适应预处理后的特征系统。从极低速到跨声速范围的众多算例表明,本文的边界条件处理稳定有效,预处理法在各种速度下都具备同样快速的收敛性,在较宽的使用范围内均能得到理想的计算结果。     

任意回转面亚、跨音流场计算方法

本文所设计的任意回转面亚、跨音流场的计算软件 ,能根据叶型和进出口条件分析流动特性 ,并自动选择精确和节省机时的合适的数值解法。用户在给出叶型坐标和定解条件后 ,可以得到叶面马赫数分布或压力分布 ,全流场的等马赫数或等压图