二维振荡机翼含激波跨声速非定常绕流广义变分原理的普遍形式及其派生族

本文建立了二维机翼作任意周期振荡时的非定常绕流的广义变分原理的更普遍形式及其派生族。文中充分发挥了“变域变分”与“自然界面条件”的有力作用,把振荡激波和尾涡面上的间断条件(例如Rankine-Hugoniot条件)也转化为自然界面条件,并兼顾了翼面的喷(吸)气作用。本文主要是为引进有限元等变分解法提供一个更普遍、更完善的理论基础,并可推广到三维机翼、二维及三维旋转叶栅上去。     

绕任意刚性翼型,机翼非定常流动的跨声速欧拉方程解

本文在有限体积方法的基础上引入变系数隐式残值光顺及固接网格法求解了绕二维翼型及三维机翼非定常流动的跨声速欧拉方程,计算结果与同类结果及实验结果吻合。     

机翼跨声速抖振数值模拟及模态分析

绕机翼的跨声速抖振流动是典型的复杂不稳定流动，对其非定常特性及失稳机制的研究具有重要的工程和学术价值。通过非定常雷诺平均Navier-Stokes（URANS）仿真方法和动模态分解（DMD）分析手段，研究了CRM（Common Research Model）等典型机翼的跨声速抖振流动特性及其主要失稳模态。数值仿真结果表明机翼的跨声速抖振表现为多失稳模式下的宽频特性。除了激波的弦向失稳，还会伴随发生激波的展向失稳，它们都表现为低频特性。翼梢处的高频响应可能是由激波诱导的低频失稳与翼尖涡相互耦合形成。DMD分析结果显示机翼展长和后掠因素诱导了激波展向失稳模态。本研究对抖振流动的物理建模、控制及理解相关的气动弹性现象具有指导意义。     

跨声速机翼非定常气动力的全位势粘位迭代计算

采用Ｃ－Ｈ型网格,守恒型非定常全位势方程的时间精确近似因式分解差分地计算二维,三维的跨声速非定常位势流,用准定常,准二维方法计算边界层位移厚度,通过粘位迭代得到的跨声速翼型,机翼的非定常气动力,所得结果与实验数据吻合很好。     

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

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

用粘性／非粘性相互作用法求解机翼跨声速分离流

提供了一种计算机翼的跨声速绕流的粘性／非粘性相互作用的计算方法,包括无粘流场计算,混合边界层计算及两者之间的相互作用,其中三维混合边界的计算包括了层流边界层,转捩区,湍流边界层和分离流的积分方法计算了,特别是在靠近分离的区域采用边界层反方法计算,无粘流场由全速势方程计算得到,通过粘流无粘流耦合迭代求得了Ｍ６机翼跨声速绕流的收敛解,与实验结果比较,吻合得较好,本方法能够计算出激波诱导分离泡和后缘分离     

跨声速翼型上激波／边界层干扰的自适应控制计算

在激波区使用自适应壁对跨声速翼型上的激波/边界层干扰进行控制,可改变机翼的气动性能。这种被动控制可通过在翼型的激波区开一凹腔,其上覆盖一弹性橡胶膜柔壁来实现。本文给出用N-S方程数值模拟这一自适应控制翼型的跨声速粘性绕流,提出了一个适用于本特殊情况(物面边界局部地区在求解过程中有变化)的处理办法。并探讨了自适应柔壁对当代跨声速翼型绕流的影响。     

机翼和叶栅非定常流的Kutta条件

本文在现有的实验研究和数值计算的基础上，对经典Ｋｕｔｔａ条件在非定常流动中的适用性作了进一步探讨，并对文献〔１，２〕中提出的一种新的Ｋｕｔｔａ条件－Ｇｉｅｓｉｎｇ－Ｍａｓｋｅｌｌ（Ｇ－Ｍ）模型进行了分析与比较，本文给出了二维振荡机翼和叶栅绕流的变域变分有限元解法中Ｋｕｔｔａ条件的处理方法。本文旨在为机翼和叶栅非定常绕流的数值计算提供更合理可靠的尾缘边界条件。     

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

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

吸气式高超声速推进助推段内流振荡及其抑制

为了掌握高超声速飞行器盲腔流场的激波振荡机理,采用"双时间步"方法对轴对称吸气式高超声速飞行器内流道盲腔流场进行了非定常数值模拟。结果表明,飞行器内流道出现了非定常激波振荡现象,内流道壁面压力出现了低频高幅的周期性变化。为了消除高超声速飞行器内流道的激波振荡现象,设计了级间段泄流方案,对飞行器进行了内外流一体化定常和非定常的数值模拟。结果表明,级间段泄流方案可以有效地消除内流道的激波振荡现象。     

计算机翼跨音速绕流的加罚Galerkin有限元法

1.引言 全速位Galerkin有限元法是跨音速有限元法中重要方法之一。处理超临界流时,为了自动捕捉到激波,通常采用上风技术,但因此导致了质量守恒条件得不到完全满足。文献[1]从Bateman变分原理出发,应用加罚函数方法,找出了强制满足质量守恒条件在Galerkin有限元方程中所应加入的附加项。本文将文献[1]的方法应用到超临界跨音速绕机翼流动,克服了一般Galerkin有限元法中由于采用上风技术而出现的激波不够陡和位     

单级跨音速压气机三维非定常粘性流动计算

本文以数值求解雷诺平均非定常N-S方程方法完成了单级跨音速压气机动静叶相干三维非定常粘性流动计算,给出了高分辨率NND格式的时间推进有限差分解。对某跨音速单级压气机进行了定常和非定常流的详细数值计算,设计转速下定常粘性流场计算所得到的性能参数与实验结果吻合较好。在此基础上,对最高效率状态进行了动/静叶相互干扰的非定常流动计算,并对结果进行了分析。      

机翼跨声速绕流的Galerkin有限元解法

本文给出了翼型和机翼跨声速绕流的Galerkin带权余数有限元解法。利用一系列适当形状的元素将计算区域离散,在元素内采用线性插值函数。对于超临界流动,文中采用了适合于有限元法的超声速上风技术,具有捕捉激波的能力。将这种技术应用于线松弛迭代解法中,在计算钝头翼型和机翼绕流时,都获得了成功。     

TIME ACCURATE 2D NAVIER-STOKES CALCULATIONS WITH A DUAL-TIME STEPPING METHOD

ＴＩＭＥＡＣＣＵＲＡＴＥ２ＤＮＡＶＩＥＲＳＴＯＫＥＳＣＡＬＣＵＬＡＴＩＯＮＳＷＩＴＨＡＤＵＡＬＴＩＭＥＳＴＥＰＰＩＮＧＭＥＴＨＯＤＹａｎｇＹｏｎｇ，ＱｉａｏＺｈｉｄｅ（Ｆａｃｕｌｔｙ５０３，ＮｏｒｔｈｗｅｓｔｅｒｎＰｏｌｙｔｅｃｈｎｉｃａｌＵｎｉ...     

翼型动失速问题研究

本文对俯仰振荡和快速上仰翼型的非定常动失速问题进行了计算研究。对绕流流场采用了高效率的分区精确数值模拟的研究方法,保证了计算结果的精度,大大节省了对非定常粘性流动至关重要的计算机时。文中单独划出前缘驻点区,对形成前缘动失速旋涡的前缘附近区域的流场进行精确计算。算例显示了动失速流场的形成、变化过程及其对气动特性的影响,计算结果和实验的比较是满意的。     

解跨音速升力翼型的有限元法

 在应用解全速位方程的最小压强积分有限元法求解绕升力翼型的跨音速流动时,将不可压流中求解绕升力翼型的耦合单位环量流动和无环量流动的解法推广到可压流中。为了确定环量,本文所用Kutta条件是:在后缘处,气流流向平行于后缘角二等分线。因有限元法对网格无正交性要求,因而可在椭圆变换前后进行剪切和延伸变换。这种网格生成法易于构成适用于复杂形状的有限元网格。通过计算并将其结果与文献中的数据比较,表明这种方法应用方便且有较快的计算速度和较高的计算精度。     

计算绕薄翼型跨音速非定常流的积分方程法

 本文从小扰动方程出发推导出绕薄翼型的跨音速非定常流动的积分方程,为使此方程适用于具有激波的流场引入人工粘性,并对其数值求解。通过算例讨论了人工粘性和计算域大小对计算结果的影响,与实验及其它数值计算结果比较,表明本方法的准确度令人满意。计算量小,收敛性好是此法的特点。     

欧拉方程计算非定常流动的一种数值方法

提出了利用欧拉方程计算非定常流动的一种数值方法。在利用中心差分的有限体积法和双时间推进法求解非定常欧拉方程过程中，通过在固定物面边界上满足完全的非定常边界条件，计算出非定常流场。算例结果表明，本文方法具有高效、简便的特点，而且能描述激波的非定常运动和变化情况。     

旋翼翼型非定常动态失速响应的计算

基于旋翼非定常翼型气动模型,给出了计算分离流和深度失速状态下的翼型非定常升力、俯仰力矩的数值计算方法。该方法采用半经验指数响应公式,利用数值离散方法来求解翼型的非定常法向力和俯仰力矩。分别计算了NACA0012和SC-1095翼型上的非定常气动载荷,并与可得到的试验结果进行了对比,验证了方法的有效性。文中还讨论了缩减频率和马赫数对动态失速响应的影响;然后,这个模型被改进以适用于后掠流下的翼型动态失速响应计算,分析了后掠角对翼型动态失速响应的影响。最后,得出了一些结论。