跨音速定常小扰动压力方程的混合差分法

本文提出直接求解跨音速定常小扰动压力方程的数值方法。对于研究某些洞壁干扰问题,与传统的速势方程相比,用压力方程作为求解跨音速流场的主管方程,边界条件为Dirichlet形式,易于处理,且待求变量为压力,可减少积累误差,提高计算精度。 本文采用混合差分法求解压力方程,通过数值试验,确定合适的差分格式及迭代线化方法。其收敛解与相应的以速势方程为主管方程求到的解相比吻合得比较好,从而证实了本文方法的可行性。 最后给出应用本文方法计算鉴定跨音速翼型风洞壁干扰以及由给定的压力分布计算翼型外形的典型算例。     

中心锥式进气道跨音速内外流场的混合差分计算

一、引言 文献[1]曾用跨音速定常势流的混合差分法成功地计算了皮托管式轴对称进气道内外跨音速流场,本文进一步用上述方法计算了中心锥式进气道的跨音速流场。中心锥式进气道的通道形状大多为收缩-扩张型,内通道壁面扰动较强。同时,由于中心体的存在,使传递进气道出口反压扰动的横向松弛线被隔断,因此,在皮托管式进气道计算中行之     

绕钝前缘翼型定常跨音速小扰动势流的一种计算方法

提高钝前缘翼型的跨音速压力分布计算方法的精度与效率对翼型设计十分重要。国内对跨音速小扰动势流方法进行了研究,但对M_∞>0.8的超临界情况尚未计算。国外文献[3]指出当M_∞>0.9时完全速势方程的一种方法失败了。 本文试图改进经典的小扰动势流方程,探索比较稳定的迭代方法和超松弛方法,以克服超临界流计算常常不易收敛的困难,使小扰动势流计算的应用范围扩大到更高的M     

关于跨音速零升力翼型绕流的松弛和两种自修正风洞的收敛性

本文采用横向小扰动而纵向大扰动速势方程,计算了跨音速零升力翼型的绕流。在线松弛的数值实验中,φ_γ的差分式用简单迭代和φ_(xx)的差分式用改进迭代时,稳定性较好。此结论与文献的线化理论分析相符。 本文用混合差分法数值模拟,证明了基于两个控制面上的静压和基于一个控制面及翼面上的静压的跨音速零升力翼型自修正风洞的收敛性。对前一种方案,NACA0012翼型,M∞=0.9,RAE104翼型,M∞=0.8,对后一种方案,NACA0012翼型,M∞=0.72,0.8,在迎角为零和风洞高度与翼弦之比为3时,均能收敛到无洞壁干扰的自由流。     

包括粘性影响的跨音速纵向大扰动非定常翼型绕流计算

本文提出一种计算效率高、并改进小扰动理论的二维跨音速定常和非定常流的计算方法——非定常纵向大扰动流速势方程和边界条件的数值解。本方法还考虑了包括边界层位移厚度以及激波-边界层干扰的粘性影响。文中给出了NACA 0012翼型和NLR 7301超临界翼型绕流的算例,计算结果与实验作了比较。     

正交流线坐标系中计算二维跨音速流的流线迭代法

本文提出了流线迭代法,用以实现文献[7]所提出的正交流线坐标系中的二维定常无粘气体动力学方程组的数值求解。 此方法的基本做法是:从物面起到流场远边界逐条流线地求出流动速度值,同时求出下一条流线的形状;按此方式对全流场进行迭代。迭代的收敛性很好。 用本文方法在零攻角和跨音速绕流的条件下计算了双圆弧翼型、NACA 0012翼型和R.A.E.101翼型。双圆弧翼型结果与实验数据符合较好。钝头翼型结果与其它方法结果相比是令人满意的。     

LTRAN2的一种改进方案

本文提出了LTRAN2的一种改进方案BTRAN2用来分析翼型的各种频率的非定常运动。用E-O调转换格式和ADI方法求解了完全的二维非定常跨音速小扰动位流方程,并用单调转换的AF2格式计算了定常跨音速小扰动方程,以此作为非定常计算的初场。本文给出了带后缘正弦振荡襟翼的NACA 64 A006翼型的绕流和做正弦俯仰振荡的NACA 64 A010翼型的绕流计算结果,它们与Euler方程解或实验数据很吻合。     

介绍一种适合于计算跨音速流动的混合差分格式

本文介绍MacCormack格式和迎风格式的一种混合差分格式。它适用于跨音速流计算。本文用它计算了绕圆弧翼型的跨音速流场,并详述了处理方法和结果。     

在正交流线坐标系中数值计算二维升力翼型跨音速绕流

本文将文献[1]提出的流线迭代法进一步应用于计算正交流线坐标系中有攻角时的二维跨音速绕流。用该算法计算了α=1°情况下的6％双圆弧翼型和NACA0012翼型,计算结果与实验数据符合较好。     

二元定常跨音速流有限差分计算中的初场问题

为研究初场对计算结果的影响,本文从二元定常跨音速小扰动方程出发,用混合有限差分松弛迭代法对NACA 0012翼型作了计算。计算结果表明,在来流为超临界的情况下,用非守恒差允格式捕捉激波时,以低熵值的初场(零值、小于计算攻角或马赫数的计算结果)进行迭代计算,就能较好地获得正确的收敛解,否则就可能得出错误的收敛解。在亚临界情况下,流场不出现激波,所以计算不受初场选取的影响,即收敛解是唯一的。     

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

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

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

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

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

 <正> 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     

激波控制鼓包提高翼型跨声速抖振边界

翼型抖振边界是仅次于升阻比的一项重要气动指标.采用定常雷诺平均Navier-Stokes方程,以升力线斜率平缓及激波位置振荡作为基本判据确定了RAE2822翼型在指定跨声速来流条件下的抖振边界.通过大量计算流体力学(CFD)验证,针对RAE2822翼型设计了一种特定外形的激波控制鼓包并确定了其具体安装位置.该激波控制鼓...     

非结构化多重网格流场数值计算

本文对二维非结构网格无粘流场计算的多重网格方法作了探讨。研究了非结构化多重网格的生成方法,插值算子的构造以及流场计算．以绕流ＮＡＣＡ００１２翼型无粘流场为算例,研究了多重网格法在非结构网格上的实现方法及其加速效果表明本文的方法有显著的加速收敛效果。     

全速位积分方程跨音速绕流的数值解

从全速位方程出发,利用Green公式将其化为激波捕获积分方程和激波装配积分方程,然后离散进行数值解。流场出现激波时,对激波捕获积分方程应用上风技术捕捉到激波,然后应用激波装配技术计算,得到了满意的结果。经算例考核,该方法具有计算区域小,收敛快和CPU时间较少等优点。     

基于特征正交分解的跨声速流场重构和翼型反设计方法研究

在二维流场的重构问题中应用特征正交分解(Proper Orthogonal Decomposition,POD)数据处理方法。利用CFD技术计算得到的流场快照对气动力模型进行降阶,然后利用基于POD的降阶模型(Reduced Order Model,ROM)对所需的流场参数进行重构,在快照范围内可以得到高精度的结果,且具有一定的外插能力。在翼型反设计问题中该方法仍然是成功的,通过修正快照向量,利用基于POD降阶模型的数据重构方法,由已知的翼型表面压力分布通过反设计就能够高效精确地得到对应的最优翼型形状,这极大地简化了翼型反设计问题。本文分别在跨声速范围对RAE 2822翼型的流场重构和Korn翼型及NACA 63212翼型的反设计进行了验证,证明了基于POD的流场重构和翼型反设计方法的有效性和高效性。     

AF-2格式的加速收敛参数研究

一、引言 Ballhaus等人提出的求解跨音速定常流动的隐式近似因式分解法（AF-2）具有收敛快速,占计算机存储少等优点。文献[1,2]分别对于完全小扰动方程和全位势方程的计算给出了很好的例证。文献[3]表明上述结论对于横向小扰动位势方程也是成立的。AF-2格式的快速收敛主要是通过加速收敛参数的循环取值来实现的。因此,合理正确地选取加速收敛参数序列对AF-2格式至关重要。     

超燃冲压发动机一维稳态跨声速流动奇异初值问题的一个解法

超燃冲压发动机燃烧室一维稳态跨声速流动方程在临界声速点存在奇异初值问题，现有的基于L’Hospital法则的求解方法在原理上存在较大的初值误差，影响一维稳态跨声速计算的精度。为此，本文提出了一种基于变量代换的改进算法，通过定义新的流动变量W=-Ma^2＋2Ma，构造出了非奇异的一维稳态跨声速流动方程。消除了微分方程的奇异性，有效的解决了一维稳态跨声速流动计算过程中存在的奇异初值问题。