高亚声速二维位势流动的气动力计算

本文提供了一个精确表达式,通过它可将亚声速位势流动仿射变换成不可压缩流动。同时给出了新的气动力计算公式。在流过Joukovski翼型的气动力计算中,它与用计算机提供的精确数值计算结果相比较,二者很符合。     

亚声速二维风洞收缩段的设计

本文给出亚声速二维风洲收缩段的一种设计方法。对两端为均匀流,风洞壁上无压力中断和逆压力梯度的条件下,设计出壁为最短的亚声速二维风洞的收缩段,借以减少边界层的影响。这个方法对低速风洞和高亚声速风洞收缩段都是适用的。对开发风力资源有一定的参考价值。     

旋成体轴对称跨声速全位势流的高效差分算法

本文用守恒型全位势方程,贴体坐标网格,对旋成体轴对称跨声速绕流的差分数值计算方法进行了研究;根据最佳收敛准则,提出了轴对称情形的AF2迭代算法,并将此算法应用于半球头柱体、弹体等各种外形的旋成体。与一般方法仅适用于亚声速自由流不同,本计算可从亚声速、跨声速自由流一直到低超声速自由流。计算结果表明,本文方法收敛快,与实验及其它方法的结果符合较好。     

大型亚声速风洞试验段设计的新概念

用一个亚声速引导性风洞的研究数据揭示了常规开口试验段设计中存在的问题及一些解决方法;介绍了国外亚声速风洞开口试验段设计的新概念,使风洞试验段允许的模型堵塞度从10%提高到20%,大幅度降低风洞的建造和运营成本,为国内未来大型亚声速风洞建设的创新设计提供参考.     

用全位势方程计算机翼的亚声速,跨声速和超声速绕流

对大后掠小展弦比细长机翼,本文对机翼纵轴垂直的每一横流截面生成O型网格,形成对机翼流场的H-O型网格,用守恒型全位势方程、差分和隐式近似因式分解迭代算法计算绕机翼的可压缩位流。自由流可从亚声速直到低超声速的全部跨声速范围。本算法要求机翼前缘有大后掠角,后缘可稍许后掠或前掠。本文算例表明,所研制的计算程序已可提供工程实用。     

跨音速非等熵全位势方法及粘流/无粘流相互作用计算

跨音速流动是现代航空军机和民机设计中必须计算的一种流动。目前采用求解N-S方程的方法尚不能提供很经济而又非常精确的解,也就不能成为工程设计使用的常规手段。在实际飞行中,高雷诺数的事实表明,粘流/无粘流相互作用计算的方法可以成为工程设计计算的有力工具。 欧拉方程是描述无粘流的准确方程,但与全位势方程相比其求解过程要复杂得多,计算费用也高。为取得最优外形工程设计,要作大量的流动模拟计算,故推广全位势方法并提高其准确度是很有意义及具吸引力的。全位势方法假定流动是等熵和无旋的,故只适用于无激波或仅有较弱激波的流场计算。为改进传统全位势方法,应对其作两方面的修正。分析表明,旋度修正比熵修正更为高阶,作为改     

快速CFD计算工具在民机概念优化设计中的应用

高亚声速民机概念设计和优化阶段需进行大量的反复迭代,其中气动力计算的准确性与时效性是一对很难处理的矛盾体。Navier-Stokes(N-S)方程可进行较为精确的计算,但需要付出的时间和资源代价很大;工程估算方法可以快速得到估算结果,但精度较差,尤其表现在跨声速阶段的波阻估算。本文采用基于全速势方程的快速CFD计算工具与工程估算相结合的气动计算方法,避免了N-S方程耗时长和工程估算精度差的问题,实现了在民机概念设计和优化程序中的集成和应用。     

跨声速风洞扩散段数值模拟

运用时间相关的全隐ＡＦ方法，计算二维及轴对称亚声速高速扩散段内流场。流动控制方程使用全Ｎ－Ｓ方程，湍流模选择了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数流模型，和两层ｋ－ε二方程模型。     

高空舱内大流量进气收缩型面气动特性数值研究

应用CFD方法对不同收缩型面设计方案进行了数值模拟,并对结果进行了对比分析。结果表明:①采用伯努利双扭线进气结构,收缩段后直段区域径向马赫数梯度较大,收缩段出口附近壁面静压变化剧烈,压力损失较大;②无伯努利双扭线进气结构,双三角函数收缩型面(DTC)壁面静压变化率较小;③不同工况下,伯努利双扭线和圆柱-四次曲线-圆锥-四次曲线组合收敛曲线(CQCQ)型线组合、伯努利双扭线和DTC型线组合、仅CQCQ型线、仅DTC型线四种结构,壁面附面层厚度在收缩段出口基本一致;④进行大流量发动机试验,高空台可参考选用仅CQCQ型线和仅DTC型线作为亚声速进气结构。     

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

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

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

本文在机翼钝前缘处用精确速势方程和精确的边界条件,其他地方用纵向大扰动而横向小扰动的速势方程和相应的边界条件,联立求解。数值算例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°。计算与实验接近。 本文建立跨音速定常小扰动速势差分方程的线松弛改进迭代在局部线化假设下的稳定性条件和松弛解收敛到原来的微分方程解的条件。这些条件大多数与数值实验相符。     

一种设计亚音速叶型的解析方法

本文在近似于势、流函数的正交面上采用解析方法进行亚音速叶型设计。将改写后的无旋条件与连续方程联立求解,可在计算面上将叶型设计问题用标准拉普拉斯方程来表示。这一方法的优点是在流场计算中用解析解取代数值迭代计算过程,求解简便,节省计算机时,并为从理论上分析各种不同叶型的设计特性提供了便利条件。文中提供的算例表明,方法是有效和可行的。      

二维定常亚音速无旋气流绕物体的流动

本文采用小扰动法给出了定常亚音速无旋气流绕物体流动的解析解,同时也给出了一个将单位圆域保角映射成一个任意封闭翼型的简单计算公式。此外还导出了压强系数的新的近似表达式。在流过鼓包的流动问题上,它与由Kaplan提出的精确解析结果相比,二者极为符合。     

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

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

旋翼三维非定常跨音速绕流的数值模拟

 利用小扰动位势方程，用有限差分方法求解了旋翼跨音速流场。导出了严格的三维非定常位势方程远场无反射边界条件，并成功地用于数值计算中，缩小了物理求解域的空间尺度。采用了贴体坐标和网格加密技术，以得到合理的非均匀网格分布。应用较少的网格点数成功地计算了无升力情况下旋翼的三维非定常亚、跨音速绕流流场，得到了与实验结果符合良好的数值解     

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

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

高压压气机叶栅的高速模型和低速模型的相似变换准则

本文以相似理论为基础,从线性小扰动方程、跨音小扰动方程和全扰动势方程出发,推导出一组在进口马赫数不同时,保持两组叶栅气动力相似的变换准则。变换涉及跨音与亚音,跨音和不可压,亚音与不可压及其全扰动的情形。应用二维和准三维无粘数值计算方法进行数值验证,证明了变换关系在无粘下的正确性。文中使用了局部线化概念和一些假设,其合理性一并由数值试验得到证实。本文在理论和数值试验上为压气机叶栅的高速模型和低速模型之间的相似变换提供了必要的基础。      

非线性波流场和波流载荷的数值模拟

本文用数值模拟的方法研究了大尺度物体上所受的非线性波与水流的共同作用力。非线性水波与水流相互作用场的计算采用有限差分法,引入坐标变换获得了固定计算域,用松弛迭代法求解差分方程。波流联合作用力的计算用时间步进法,每一时刻采用简单格林函数的边界元法求解。入射势采用以上波流场的计算结果,开边界选取离物体足够远,其上的速度势为入射势。