跨音速机翼数值优化设计方法

优化设计方法由机翼外形和计算网格生成方法、跨音速机翼全位势方程数值解法和数值优化算法等方法组成。采用独特的机翼外形曲线拟合函数摄动方法来改变机翼。数值实验表明,只须在半翼展中部剖面上选择两个结点处的翼面坐标进行摄动,即可获得气动上满意的优化解答,并可节省大量机时。     

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

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

计算跨音速三维机翼非定常载荷的核函数法

 本文给出了在跨音速流动中,计算三维机翼非定常载荷分布的核函数法和计算结果。此方法是基于线化升力面理论,利用在亚、超音速情况下,联系载荷分布与下洗的积分方程的核函数,采用局部线化法来处理跨音速的非定常绕流问题。在计算中引进了线分布的冲波偶极子,并嵌入适当的正冲波边界条件。 我们从统一的超、亚音速核函数表达式出发,对亚音速区数值积分作了新的简便的处理,对于亚、超音速区相互的诱导下洗给出了自己的数值积分过程。我们算了三个例题,和文献上的数值符合较好。     

机翼跨声速抖振研究进展

跨声速飞行中,激波附面层干扰会引起激波周期性自激振荡,这种现象称之为跨声速抖振。跨声速抖振引起的脉动载荷有可能造成结构疲劳甚至引发飞行事故,该问题一直是航空领域的研究难点。本文从风洞试验和数值模拟两个方面总结了跨声速抖振研究的主要方法,重点综述了近几年在跨声速抖振机理研究方面以及结构运动对抖振响应特性的影响方面取得的新进展,简要总结了跨声速抖振的被动控制及主动控制方法。在此基础上,提出了跨声速抖振后续研究的重点方向:抖振机理的进一步探究,结构弹性特征对抖振特性的影响以及闭环主动控制等。     

关于“跨音速绕流的边界积分法”

文献[1]是一篇参加国际会议的文章,由第二作者在会上宣读。当时,我们只准备在会议上透露一信息:我们已可建立可压缩非线性位流方程的边界元解法。对方法的中心思想无意细说,拟在此项研究工作结束后再正式发表专文系统叙述。所以文献[1]写得比较简单显得晦涩。现仅就文献[2]和文献[3]对文献[1]所提的问题作如下补充说明,我们今后发表的正式文章中对此将作更全面的论述。     

绕钝前缘机翼跨音速势流的松弛计算及其稳定性与收敛性讨论

本文在机翼钝前缘处用精确速势方程和精确的边界条件,其他地方用纵向大扰动而横向小扰动的速势方程和相应的边界条件,联立求解。数值算例1为矩形机翼,展弦比λ=12,翼剖面为NACA0012,自由流的马赫数M_∞=0.63,迎角α=2°,翼根剖面压力分布的计算结果与二元亚音速精确数值解（Sells,1968）接近。算例2为NACA RM A51G31实验的机翼,垂直于1/4弦线的翼剖面为NACA64A010,其后退角χ_1/4=45°,λ=3,根梢比η=2,M_∞=0.4,0.8,0.9,α=2°。计算与实验接近。 本文建立跨音速定常小扰动速势差分方程的线松弛改进迭代在局部线化假设下的稳定性条件和松弛解收敛到原来的微分方程解的条件。这些条件大多数与数值实验相符。     

跨音速机翼设计——非设计状态的效率和虚拟气体参数的影响

 改善下一代民航机空气动力效率的技术已成为航空科学中重要课题之一。高气动力效率要求在超临界Mach数下飞行,并要避免激波的出现或使激波变得足够弱。跨音速设计方法可以实现无激波或接近无激波的流动。解决跨音速设计问题往往是和发展可靠的流动分析方法相联系的。“椭圆连续”或“虚拟气体”方法是一种有效的设计方法。此方法及其应用于翼型和机翼的设计已概述于文献[1,2]。这种概念首先由H.Sobieczky给以说明:用一种数值计算方法来求解绕基本翼流动的一组虚拟气体的方程组,这些方     

关于“跨音速流的边界积分法”一文的再商榷

<正> 航空学报第10卷第3期上刊载了三篇文章,展开了对杨岞生、吴镇远的“跨音速流的边界积分法”一文的讨论。文献[4]对文献[2,3]提出的问题进行了答辩。在答辩的文章中主要谈了文献[1]只是以冲量场出发建立边界上可压流与不可压流ρν之间的关系。我们认为,对于任何流体力学问题,制约流场的基本方程和边界条件是紧密相联的整体。原作者在文献[4]中承认:“对于可压位流,实际定常流场的ρν与冲量流场的ρν不同”,但仍利用冲量流场来建立边界上可压流与不可压流ρν值的关系。因此,边界     

关于“跨音速流的边界积分法”一文的商榷

 近年来,杨岞生、吴镇远提出了一种新的解决非线性可压缩流的边界元法,他们在引进压力冲量之后将一个非线性问题转变为线性问题,然后利用边界积分表示式给出了二维、三维机翼跨音速绕流的计算方法。若干数值算例说明该方法的结果与其它数值方法和实验结果相比符合得较好。由于该方法能够将一个非线性可压缩位势流转换为对压力冲量的线性方程,因此理所当然地引起了国内外学者的关注。     

关于“跨音速流的边界积分法”一文的某些问题

 最近杨岞生教授和J.C.Wu教授提出了一种跨音速流的边界积分法。该方法以脉冲压强为基础建立非线性可压缩位流的边界积分表示式,利用边界元法给出了二维、三维机翼跨音速位流的计算方法。我们认为文献[1]在推导该方法的基本方程（即主管脉冲压强的方程）和求此方程的基本解这两个方面存在一些问题。     

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

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

机翼跨声速绕流的压强极小积分有限元算法及应用

本文从全速位方程出发。利用压强极小积分,给出了机翼跨声速绕流的有限元算法,并编制出相应的FORTRAN语言程序。文中给出了适于有限元法的网格生成技术,采用了线松弛迭代法和逐步加密法进行数值求解,提供了M6机翼和一个三角翼的计算结果。     

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

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

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

跨音速湍流边界层测量

 热线探头在中国空气动力研究与发展中心(CARDC)FL-23跨超音速风洞中进行了校准。结果表明;热线过热比a_ω≥0.6时,密度灵敏度系数S_ρ基本不随a_ω变化;当a_ω=0.8时,在Ma=0.6～1.16范围内S_ρ基本不随Mach数变化。用校准后的热线探头测量了FL-23风洞试验段侧壁边界层的平均气流和脉动气流特性。平均气流特性与皮托管测得的结果一致性很好;速度脉动u′/和密度脉动p′/随距洞壁距离y的增加不单调减小。     

跨音速叶栅流场的无粘-边界层迭代计算

二维跨音速叶栅流场的计算是在无粘时间相关法计算程序的基础上采用无粘-边界层迭代而得。马赫数从亚音速到超音速范围的几种数值试验均与测试结果吻合一致。在边界层无强烈脱流时可获得很好的结果。     

大尖削机翼非定常跨音速绕流的数值解法(英文)

本文给出一种计算飞机机翼上定常和非定常跨音速气动力的数值解法,使用了一种特殊设计的坐标变换,是用时间精确交替方向隐式(ADI)有限差分算法来求解非定常跨音速修正三元小扰动位势方程。给出了F5战斗机机翼的数值结果并与XTRAN3S,ATRAN3S及试验结果进行了比较,表明本方法是有效的和经济的。     

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

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

关于跨声速流动的敏感性

本文从大量实验事实出发,分析了细长锥-柱体跨声速流场结构随来流马赫数M_∞的变化。指出这种结构外形圆柱段上的流动随来流M_∞数的变化呈现一种特殊的敏感特征,由此提出了跨声速上、下敏感马赫数的概念。研究表明,圆柱段上的流动并非在整个跨声速范围都十分敏感,其真正敏感区域是介于上、下两个敏感马赫数M_(us)、M_(ls)之间的区域。这一特征对跨声速流动中的若干结果有一定概括性。对跨声速领域理论分析、数值计算、开展实验及实际飞行均有一定指导和参考意义。     

考虑粘性相互作用的跨音速机翼分析

本文将一些已有的三维全位势方程计算方法和一种三维积分边界层方法作了弱相互作用的耦合计算。计算了三种机翼并与已有的试验数据作了比较。