亚、跨音速机翼携带外挂物气动力干扰的混合差分计算

本文应用跨音速、定常、小扰动势流的混合差分方法,计算了机翼—挂架—外挂物气动力干扰。利用线超松弛改进迭代方法,在物理空间网格点上满足x向大扰动速势方程,在物面上满足精确的边界条件,在涡面上满足库塔条件,远场处取速势方程的线性解。计算网格点上的速势值ψ、下洗速度ψ_y、侧洗速度ψ_z:的分布,以及所有部件、组合体的压强分布、气动力系数,及机翼、外挂物各自所受的气动干扰量。 本程序用BCY语言在上海华东计算技术研究所655机上进行计算。文内三个算例均得到收敛或接近收敛的结果。与可以找到的风洞实验比较尚一致。     

三角形机翼前缘扭曲问题的理论和实验研究

本文研究了对较小展弦比的三角形机翼采取适当的扭转和弯曲,推迟气流前缘分离,使之在设计状态下实现最小诱导阻力的问题。文中理论部分,应用亚音速锥型流理论,提出了一种设计最小诱阻三角形扭曲机翼的解析方法,并对两个不同设计状态的算例进行了数值计算。实验部分,对带有三角形前缘扭曲机翼（设计状态:M_∞=0,C_L=0.20）和三角形平面机翼的两种翼身组合体,进行了纵向三分力试验,作了理论和实验结果的比较分析。结果表明:按本文方法设计的前缘扭曲机翼,其翼身组合体在设计状态下,实现了理论预期的最小诱导阻力值,组合体的升阻特性有了显著的改善。     

旋转体跨音速零升波阻和压力的差分计算

应用Murman-Cole的有限差分法,求解具有纵向大扰动而横向小扰动的跨音速轴对称速势方程,由此计算旋转体跨音速零升力时的压力分布和波阻力,以及激波位置。物面边界条件被转移到物体轴上。远场边界条件由无穷远处的条件近似代替。计算物面压力系数时,用细长体理论进行物面速势插值。 速势的差分方程用沿半径方向线超松弛改进迭代求解。网格取62×16,迭代初场取零,达到收敛的迭代次数对M_∞<1,M_∞>1以及M_∞≈1分别大约为150,40和300次。松弛因子取为:M_∞<1时,0.9≤ω_b≤1.7,0.9≤ω_p<1.0;M_∞≥1时,0.8≤ω_b≤0.9,0.8≤ω_p≤0.9,这里ω_b,ω_p分别为局部亚音速点和超音速点的松弛因子。 算例为七种不同外形的细长体,计算结果与实验符合尚好。 文中对网格、初场、迭代方法、松弛因子等有关收敛性、收敛速度问题进行了探讨。在局部线化条件下,对定常小扰动轴对称势流的差分方程,进行了线超松弛迭代的稳定性和收敛性分析。数值计算经验与理论分析所得结论相符。     

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

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

机翼跨音速无粘绕流数值计算

本文介绍机翼跨音速定常无粘绕流的一种数值计算方法。选用精确速势方程作为问题的数学模型。在直角坐标系中,经过适当的坐标变换,先把后掠翼变成矩形翼,再把无限的物理空间变成有限的计算空间。在计算空间中,用有限差分格式(亚音速区用中心差分,超音速区用旋转差分)对精确速势方程离散化。差分方程形成的代数方程组用线松弛迭代求解。     

翼型跨音速绕流的数值计算

本文选用精确速势方程作为翼型跨音速无粘绕流的数学模型。在变换过的直角座标中,用有限差分(亚音速区用中心差分,超音速区用旋转差分)对精确速势方程离散化。差分方程形成的代数方程组用列松弛迭代法求解。 为了计入粘性效应,利用精确速势方程和附面层动量积分方程进行联合迭代求解。这一点对超临界翼型的计算显得特别重要。 算例与其它数值解及实验作了比较,对于一般翼型和超临界翼型,本文结果是良好的。     

简化自适应机翼的气动外形优化设计

自适应机翼技术的研究对于超声速飞机设计具有重要意义。本文在二维翼型自适应的研究基础上,用Powell法优化计算了对前缘后掠角为35°,展弦比λ=3.9,梢根比η=0.17,机翼剖面为NACA65006翼型的梯形翼的前后缘舵面偏转角,从而获得了在亚跨声速时升阻比大而在超声速时阻力系数小的简化自适应机翼的最优气动外形。本文采用了并行遗传算法,计算了要求亚跨声速升阻比大同时超声速阻力小的气动双目标优化机翼的外形。讨论了优化机翼相对于原始机翼的气动增益。与二维一样,三维数值算例也证明了自适应机翼可获得明显的气动增益。     

跨音速组合体纵向空气动力的差分计算

 本文建立机翼-机身-平尾-立尾组合体跨音速纵向空气动力的混合差分线松弛迭代计算方案。在物理空间采用直角座标网格,列出组合体的远场边界条件,在流场内部和物面、涡面网格点上,满足速势方程,利用物面边界条件,将速势延拓到机身内部。做了两个算例,得到了收敛或接近收敛的结果,与实验尚符合。对收敛过程和超松弛加快收敛的方法作了一些研究,发现超音速点的松弛参数应取大于1。     

绕翼型跨音速大扰动对称势流的插值混合差分法

 本文将跨音速定常小扰动势流混合差分法推广到跨音速大扰动定常势流,提出了在局部速度坐标系中求解跨音速精确势流方程的插值混合差分法。作为算例,计算了双圆弧翼型和NACA0015翼型对称问题压力分布,并与已知实验值和双圆弧翼型小扰动混合差分法计算值进行比较,结果接近。试算表明,本文提出的插值混合差分格式是稳定和收敛的。本文解决了M_∞趋近于1的计算难点。     

亚、跨音速导弹气动力的有限差分计算

本文介绍亚、跨音速导弹纵向气动力特性计算的有限差分方法。采用园柱座标系,用变步长的混合差分格式进行线超松弛改进迭代计算。速势方程考虑了x向大扰动,并把精确边界条件嵌入到头部的速势方程中,以适用于计算钝头外形的要求。对尾部收缩的弹体,翼涡沿其平面伸展到收缩部分,以模拟涡面对尾翼和尾段的影响。 所建立的程序适用于计算任意形状旋成体弹身,翼身组合体及翼身尾组合体的表面压力分布和纵向气动力系数。 为鉴定方法的有效性和程序的正确性,对三个不同的气动外形进行了计算,与风洞实验结果比较,基本符合。     

后掠机翼跨音速绕流的二级近似方法

本文采用TSDH方程计算三维后掠机翼的跨音速绕流,考虑了适用于机翼钝前缘的前缘边界条件和前缘速势方程。采用Jameson格式在不等步长格网中的推广形式,把TSDH方程离散化为差分方程组。然后,在整个计算空间内布置稀网,在机翼附近再布置密网,进行稀密网的交替迭代,以加速收敛和提高计算精度。对ONERA M6机翼的超临界无激波和有激波情况的计算表明,TSDH解与FVP解和风洞试验符合良好。     

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

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

亚超音速机翼最佳弯扭综合设计计算方法

本文介绍一种亚超音速机翼最佳弯扭综合设计的计算方法,它应用了有限基本解方法。分别在亚超音速各选取一个设计点(M数和C_L),进行机翼弯扭设计,其目的是减小与升力相关的阻力。在此基础上,顾及亚音速和超音速这两个设计点的气动力特性,还要兼顾到飞机其它性能和结构上实现的可能性,进行机翼的综合设计。本文分别给出了亚音速最佳弯扭设计,超音速最佳弯扭设计和综合设计的计算结果。经过分析表明,计算结果是合理的。     

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

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

求解平面翼型亚,跨声速绕流的一个新方法

本文提出了求解平面翼型亚、跨声速绕流的一个新方法。引入流函数和Ｖｏｎ Ｍｉｓｅｓ变换后，亚、跨声速绕机翼无旋流动的基本方程组被化为以流线纵坐标ｙ为未知量的单个二阶偏微分方程－流线控制方程。并通过变换将物理平面上的无限域变为计算平面上有限的矩形域，而后在计算平面采用有限差分线松弛迭代法求解。作为算例，计算了对称翼型ＮＡＣＡ００１２－３４和非对称翼型ＮＡＣＡ４４１２的亚、跨声速有攻角绕流，所得数值结果     

绕三维物体亚音速振荡位流及其在飞行器动导数计算上的应用

 本文介绍处理不同外形三维物体亚音速振荡绕流一种统一的方法。本方法的主要特点是采用有限元素法直接解由格林定理导出的物面速势积分-微分方程以求得物面的速势分布,然后再用有限差分法对速势进行微分求物面的压力分布。 由于本方法理论上比较严格,适用于复杂外形物体绕流的计算,所得结果又比较准确,因此近几年来在国外得到越来越广泛的应用。在本文中,采用了与有关文献相同的基本方程,但在气动影响系数的计算上略有不同,本文并将这一方法应用于飞行器动导数的计算,所得的结果与实验结果符合。     

三维跨音速流有限差分计算中的某些问题

一、基本方法简介 1.三维跨音速小扰动方程 柱坐标下的三维跨音速小扰动方程为     

纵向大扰动方程的变参数型变分裂法

 本文用变参数型变分裂法,解定常跨音速纵向大扰动方程。计算翼型气动力。对变参数型变分裂法以及它的参数选取方法作了必要的介绍。对NACA0012翼型,M_∞=0.8,α_A=0°;M_∞=0.95,α_A=0°;M_∞=0.75,α_A=1°,进行了计算。所得结果与其它结果进行了比较,表明本方法是快速收敛的、可靠的。     

伸缩机翼变体飞机气动特性研究

伸缩机翼变体飞机通过机翼伸缩调整机翼展长,从而改变机翼面积和展弦比,改变飞机的气动布局和机翼的气动特性,满足多任务点的设计要求。简要介绍伸缩机翼变体飞机的发展历史,重点研究一种采用伸缩机翼设计的超音速飞机的气动特性变化。研究结果表明：亚音速时机翼展长伸长,展弦比增大,飞机诱导阻力降低,升阻比提高,可以明显提高飞机的航程;超音速时机翼展长缩短,展弦比减小,飞机的波阻降低,升阻比增大,提高了超音速飞行性能。伸缩机翼概念用于超音速飞机设计时能很好地兼顾亚音速巡航和超音速冲刺。