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

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

二维弹性扑翼沉浮运动流动特性的数值研究

采用非定常势流板块法和Euler-Bernoul1i梁振动微分方程,对后缘带弹性薄板的翼型在作沉浮运动时翼型绕流和薄板弹性振动之间的流-固耦合过程进行数值模拟,并计算流动特性、非定常气动载荷,尤其是沉浮运动的推力效应.用有限差分法求解薄板的弹性变形运动,气动弹性分析采用弱耦合迭代方式.通过对翼型非定常气动力(升力和推力)、推进效率、尾迹脱落涡计算结果的分析发现,通过调整弹性薄板的刚度,可以使得翼型产生最佳的非定常推力.     

用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

翼型纵向俯仰非定常气动力研究

飞行器飞行时受到的非定常气动力是飞行器设计以及机动飞行中须考虑的一个重要因素。采用CFD方法数值模拟了NACA0012翼型在马赫数从0.2~5.0范围变化时做强迫振荡的非定常流场。通过积分法计算了在不同马赫数时翼型的俯仰动导数。计算结果表明,迟滞效应随马赫数的增大先增大后减小,在跨声速区域迟滞效应达到最大;迟滞环的绕向在跨声速区域为顺时针方向;迟滞环顺时针绕向时,动导数大于零,模型动不稳定,逆时针绕向时,动导数小于零,模型动稳定;旋转角速度越大,翼型表面非定常气动效应越明显。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

高超声速验证飞行器助推分离段流场数值研究

为了模拟有相对运动的多体非定常绕流，发展了基于弹簧近似的非结构动网格方法及耦合动网格的Euler解算器，作为方法验证，模拟了俯爷振动的NACA0012翼型绕流，计算结果与实验及文献结果非常接近。采用二维近似外形，对高超声速验证习行器助推分离段流场进行了数值研究，得到了不同时刻由于多体相对运动形成的干扰流场结构以及分离过程的气动力参数。     

振荡翼型非定常流动数值模拟研究

进一步完善了非定常雷诺平均N-S方程求解程序.并借助运动网格生成方法,对周期性俯仰运动下NACA0012翼型的非失速绕流和动态失速问题进行了数值模拟研究,分析了振荡翼型的迟滞效应以及俯仰运动各个参数对迟滞效应的影响.数值模拟结果与实验结果对比表明,本文发展的数值模拟方法在模拟非失速绕流时是可靠的;结果还表明,随着振幅、折合频率的增加迟滞效应增强,平均攻角对迟滞效应无影响;保持无量纲参数--振幅和折合频率的乘积--恒定时,迟滞效应相同.在模拟动态失速问题时.本文方法在深度失速下结果与实验符合较好,而在轻度失速下符合较差.经分析表明这是不同条件下非定常流动对湍流模型和转捩的敏感程度不同所致.     

含运动边界二维多体干扰流场的数值模拟

本文采用二阶精度Godonov扩展方法求解非定常：Euler方程，以基于格心的自适应非结构网格的显式有限体积法为基础，空间和时间都是二阶精度，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，对作俯仰振动的NACA0012翼型绕流及激波管中运动、静止圆柱的流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步对含多体分离流场进行数值模拟打下了基础。     

多段翼型非定常地面效应数值模拟

发展了一种多段翼型着陆、起飞过程中非定常地面效应的数值模拟方法。针对着陆过程中翼型大位移运动问题,提出尽量保持翼型周围和尾迹区域网格随翼型做刚体运动的动网格策略,并采用一种简单的网格重构方法以保证翼型在整个下降过程中具有合理的网格分布。进一步建立起多段翼型非定常地面效应N-S方程计算方法,可应用于起降过程中多段翼型考虑非定常地面效应的实时数值模拟。NACA0012翼型和GAW-(1)两段翼型匀速下降过程的非定常升力计算结果均和相应准定常状态(下降速度折合为来流迎角)相吻合,说明了方法有效性。GAW-(1)两段翼型非定常和下降速度折合迎角加至翼型的准定常计算比较表明:多段翼型升力均随离地高度减小而降低,而非定常地面效应的影响随着离地高度的减小先小于准定常地面效应,而后有所增大。     

粘性流绕俯仰振动翼型流动的数值模拟

本文直接从N-S方程出发,利用LU-ADI格式和Baldwin-Lomax代数湍流模型研究了粘性流绕俯仰振动翼型的流动。数值实验表明,本文的数值结果同实验吻合较好。通过数值模拟手段研究了振动对流场中激波和分离这两大主要特征的影响,结果发现;(1)激波滞后于攻角的变化,如当NACA0012翼型在负攻角状态下上翼面存在激波而下翼面无激波。(2)由于翼型的振动,分离被减轻升力相对提高。(3)随翼型攻角和振幅的变化,翼型振动的升力回线走向可不同。     

NPU翼型的气动力分析和改进设计

 在飞行器设计中用计算方法设计超临面翼型已完全取代了选用现成翼型的设计方法。为考察已设计出的NPU翼型是否满足飞行器设计要求我们对其进行了全面气动分析,发现这些翼型尚有不足之处,有必要进行改进设计。     

多段翼型低速流场与气动特性数值计算研究

在结构化网格求解RANS方程的基础上,将嵌套网格技术应用在低速绕流气动特性计算中,并以NLR-7301两段翼型和某运输机5段翼型在起飞构型下对低速绕流流场进行数值模拟,将气动特性及流场特性的数值模拟数据与对接网格计算以及风洞实验数据进行了比较.结果表明,该气动特性计算方法有着较高的精度和可靠性.     

基于太阳能飞机应用的低雷诺数翼型研究

以太阳能飞机为背景,对低雷诺数翼型FX 63-137进行了折线型建模以拟合典型晶硅太阳能电池片对气动外形的影响,开展了气动数值模拟分析。首先引用"拟合优度"的定义描述本文折线型翼型轮廓与基准翼型(Baseline)的吻合度,并以此参数为变量建立5种折线型翼型模型;然后,采用计算流体力学(CFD)方法计算分析了不同雷诺数下各折线型翼型的气动特性,并着重研究了低雷诺数下折线型翼型的绕流机理;最后,基于工程应用实际的需求,提出了晶硅太阳能电池片的铺设方法也即折线型翼型设计思想准则,并进行算例验证。研究结果表明:低雷诺数条件下,折线型翼型升阻性能相比光滑翼型在一定程度上表现出了优势,但随着雷诺数的增加,升阻方面的优势逐渐消失;折线型翼型压力分布受各折线段长度影响,前缘吸力峰值、压力平台范围以及压力恢复区分布特征是决定折线型翼型气动性能的主要因素;通过设计的算例验证了本文提出的折线型翼型设计思想的可行性。     

油膜干涉测量翼型壁面摩阻低速风洞试验技术

油膜干涉试验可以定量测量得到模型壁面的摩擦应力,是开展摩阻特性研究的有力工具。通过对油膜干涉试验技术中的试验原理、试验方法、误差修正、数据结果分析等关键问题的研究,初步建立了试验系统,具备了开展模型壁面摩阻特性研究的试验能力。综合分析表面热膜试验结果与数值计算结果,在以模型弦长为特征长度,试验雷诺数 Re ＝1．3×106的条件下,翼型表面摩擦应力沿模型弦向的分布规律具有较好的一致性。     

应用嵌入式H型网格技术计算孤立翼型绕流场

提出一种嵌入式H型网格技术,用于生成孤立翼型和叶栅的计算网格。与传统的H型网格相比,采用嵌入式H型网格时,翼型头部及尾部附近的网格畸变很容易得到控制,有益于提高计算的准确性。      

重型燃气轮机压气机多截面叶型优化设计

出于对叶型沿叶高方向的不同以及不同叶高处气动边界的不同两方面考虑,开发了一套叶片多截面叶型优化的设计流程.以单级高负荷压气机的进口导叶和大折转角静叶为优化对象,以多截面关联的性能参数为优化目标,采用考虑附面层转捩的MISES计算程序对叶型气动性能进行分析,并通过遗传算法对叶片进行优化.结果显示:在有效工作范围内,各截面叶型经优化后，前部负荷提高,中部及后部负荷降低,性能得到一定提升.三维数值模拟结果显示,在无大尺度二次流动的条件下,多截面叶型优化方法对提升叶片性能切实可行.      

Calculation of the transonic dip of airfoils using viscous-inviscid aerodynamic interaction method

A numerical procedure for the calculation of the transonic dip of airfoils in the time domain is presented. A viscous-inviscid aerodynamic interaction method is taken to calculate the unsteady aerodynamic loads. In the present case the integral boundary layer equations are coupled with the Transonic Small Disturbance (TSD) Potential Equation. The coupling between structural motion and aerodynamic loads is carried out using State Space equation. It is solved by State Transition Matrix technique. Results are presented for NACA 64A010 and NLR 7301 airfoils with structural data from Isogai and DLR, respectively. Comparisons show good agreement with other numerical results. Certain deviations of experimental data taken from literature need more insight in the detailed test conditions.     

二维振荡机翼跨声速非定常绕流的变域变分有限元解

本文首先对文献（１－３）提出的二维振荡机翼含激波跨声速非定常绕流的变分原理进行了改造，使之能适用于时间推进法；构造了三维时空的可变节点有限元来捕获自由尾涡面，而跨声流中的激波采用人工密度法捕获。在远场边界上使用简化的无反射条件和新型非定常Ｋｕｔｔａ条件。用本文方法求解了作俯爷振荡的ＮＡＣＡ６４Ａ０１０翼型的跨声速绕流，计算结果令人满意。本文方法可以推广到三维机翼及二维和三维叶栅的同类气动力问题上去     

临近空间螺旋桨低雷诺数高效翼型数值分析

针对临近空间翼型边界层变厚更容易发生分离的特点,在修正了湍流模型后,通过分析比较13种低雷诺数翼型的升阻系数得到了比较适合用于临近空间螺旋桨叶素的高效翼型。结果显示,攻角在-6°～8°、8°～14°、14°～22°三种工况下,分别采用SST k-ω、RNG k-ε、Realizable k-ε湍流模型,可以得到与国外实验比较接近的合理结果;工况条件对翼型气动性能好坏的影响很大,翼型S-1223和FX63-137在所研究的工况内都具有较好的气动性能;因此可以选择这两种翼型作为临近空间螺旋桨用高效翼型。     

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

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