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

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

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

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

翼型—扰流片的分离气动特性计算

本文介绍了用涡面元法模拟带扰流片的翼型低速无粘分离绕流。在翼型和扰流片的面元上分布线性变化涡。在扰流片后的上下分离流线的面元上分布等强度的涡。上分离流线始自扰流片的梢部，下分离流线自翼型的后缘引出。分离所泡由两离散涡结尾。气泡内总压为常值，它与涡强大小一同求解。分离气泡的形状在迭代求解过程中确定。压强分布和升力系数的计算值与现存文献的数值结果和实验数据是一致的。     

多段翼型增升小翼纳维-斯托克斯方程分析

装有扁平板增升小翼的多段翼型流场数值分析已经完成。研究的小翼装在接近主翼段后缘，高度范围是0．25％～1．25％参考翼型弦长。二维数值模拟在嵌入式结构拓扑网格上采用不可压缩纳维-斯托克斯（N-S）方程求解分析。新的网格细化提高了网格边界重叠部位的计算精度。在固定雷诺数条件下，对装有窄缝襟翼的二段翼型研究了各种小翼的影响。计算和实验结果都表明，相对未装小翼的基本翼型，主翼段凹陷处装有小翼提高了最大升力和升阻比。计算得出的流线分布揭示，正是小翼的存在引起了襟翼流场分离的减弱。也研究了三段翼型在一个雷诺数范围内的流场情况。对优化的襟翼来说，计算和实验所得的雷诺数影响结果一致。当襟翼从最有利的位置移动时，数值分析表明，小翼有助于翼型重新优化，大概能起到所优化襟翼效益的1％量级的作用。     

多段翼型失速分离特性的计算

本文用高阶奇点分布板块法(Panel Method》及边界层理论,对多段翼型的失速性能进行计算。采用分离尾迹的模型,按尾迹是流线用迭代法求出分离尾迹边界(此为内迭代)。在尾迹上配置奇点,然后把尾迹考虑在内进行粘/位流迭代(此为外迭代)。收敛后可得多段翼型空气动力特性,且可计算到失速以后的情况。算例有单段翼型GA(W)-1,NACA 4412,和F111飞机的三段翼剖面,所得结果与实验值能很好地符合。     

柔性自适应壁风洞的翼型实验技术

给出一种新的二维跨音速柔性自适应壁风洞实验迭代方案,计算方法和实验验证结果。根据实验时实验段上下壁和模型上下表面实测压强分布对风洞内、外流场进行非线性数值模拟;计算流线化壁面的形状,进行自适应实验。用该迭代方案,在堵塞比ε=8％,实验段高与翼型弦长比值H/c=1.5情况下,对NACA-0012翼型进行了高亚音速验证性实验。实验结果与国外大风洞无干扰实验结果吻合很好。实验时迭代次数仅需1～2次。实验结果展示了自适应壁风洞实验技术用于翼型跨音速实验的前景。     

跨声速流场计算的时间分裂有限体积法

本文用时间分裂有限体积法计算了翼型和平面叶栅的跨声速流场,计算与实验值的对比表明,数值结果是比较准确的。计算应用了滞止焓沿流线不变的关系来代替能量方程后,给出了差分方程的稳定条件。分析了跨声速流计算中,差分格式引进人工耗散的必要性。此外,本文还对加速计算收敛的方法及其效果作了探讨。     

翼型失速的非线性动力学特性及其控制

以NACA0012翼型为研究对象,应用CBS( Characteristic-Based Split)有限元方法对翼型绕流流场进行数值模拟,着重对翼型失速现象的动力学机理和特征进行研究.首先给出了静止翼型在不同迎角下流场的瞬时流线图以及静态翼型的迎角-升力系数关系曲线,并且分析出此曲线具有滞后、突跳等典型的非线性动力学...     

扑动翼型的低雷诺数气动特性分析

通过求解引入拟压缩项的不可压Navier-Stokes方程,数值模拟了绕扑动翼型的低雷诺数非定常流动。针对厚度在4%-12%之间的NACA对称翼型,分析了翼型厚度等参数对扑动翼型气动特性的影响。在低雷诺数条件下,对于纯俯仰运动,随着翼型厚度的减小,平均阻力系数也变小。而对于纯沉浮运动,发现翼型厚度对气动特性的影响和俯仰运动有很大的差别,平均阻力系数随着翼型厚度的减小而变大。通过对沉浮运动一个周期流线图的分析,认为这是翼型前缘涡的影响造成的。由于前缘涡的影响,翼型厚度增加,平均压差阻力系数变小,甚至会出现负值。雷诺数的影响研究表明,随着雷诺数的增加,扑动翼型的阻力系数减小的趋势越缓慢。     

变形体超临界翼型动态气动特性

为实现飞行器外形根据来流条件变化做出相应调整,需深入了解变形体非定常气动现象和机理。以典型超临界翼型为研究对象,分别制定了翼型厚度及弯度变形方案,实现了翼型柔性变形。利用数值仿真方法在高雷诺数条件下开展了翼型厚度、弯度连续变形带来的非定常效应。结果表明,翼型厚度、弯度连续变形均会形成明显的升、阻力系数迟滞环,且频率越高、幅度越大,变形导致的非定常效应也越强;相对而言,翼型弯度变形产生的非定常效应较厚度变形产生的非定常效应更明显。最后通过分析流线图及压力系数分布图发现,流动结构随几何变形的迟滞性会导致非定常效应的产生。     

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

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

Gurney襟翼对翼型流场影响的数值模拟及研究

本文采用计算流体力学的方法计算了NACA23012翼型以及安装四个不同高度（1％C、2％C、3％C、4％C）Gumey襟翼翼型的流场,并比较了不同来流马赫数、气流攻角条件下的气动性能．包括升阻比、翼型表面压力系数分布和马赫数分布。计算结果表明,安装Gumey襟翼后翼型的升阻比得到提高,升力分布趋于平均,在所计算的气流条件下安装3％C高度Gumey襟翼的翼型获得了最高的升阻比。     

导体与非导体介质翼型的缩比模型RCS数值计算分析

应用二维电磁散射有限元算法理论考虑了几组导体和非导体介质翼型的缩比计算模型,计算出缩比模型的RCS并与1:1模型进行对比.结果显示,导体翼型的缩比模型RCS计算结果与推及的二维导体外形缩比模型测量定律完全符合;非导体翼型的缩比模型在等效阻抗增大时,按缩比测量定律计算的RCS与1:1模型间的相对误差越来越大,但在等效阻抗小于0.5时缩比定律是有效的.     

低雷诺数下吹吸气射流对NACA0012翼型气动性能的影响

基于Favre过滤的大涡模拟方法,对雷诺数Re=104,迎角α=6°下的NACA0012翼型上表面吹吸气射流进行了数值模拟,从翼型周围流场流线图、速度场云图、上下表面压力系数曲线以及上表面边界层位移厚度等多角度地分析了射流位置以及速度变化对翼型气动性能的影响。结果表明:射流位置对翼型气动性能影响较大,且吸气射流要明显优于吹气射流。对于吸气射流,前缘吸气要明显优于中后缘吸气,可有效增升减阻,并减小翼型尾部流动分离,抑制翼型气动参数扰动,其最佳吸气位置随着速度的增大逐渐向下游移动;而吹气射流对翼型气动系数的作用效果较差,但中后缘的吹气射流可减小飞行过程中的气动扰动量,且吹气越大,效果越明显。     

考虑转捩影响的翼型动态失速数值模拟

用数值求解雷诺平均NS方程的方法来对考虑转捩影响的中低雷诺数下振荡翼型动态失速进行数值模拟,计算中采用了k-ω SST两方程湍流模式,并加入Chen-Thyson转捩模型来模拟流动中的转捩效应.采用该方法分别对雷诺数Re=1.35×105和Re=7.7×104情况下NACA0012翼型的动态失速进行了数值模拟.计算结果显示:计算出的翼型动态失速气动力系数迟滞曲线与实测结果符合较好;对于Re=7.7×104的工况,实测的升力系数迟滞曲线中出现了高频振荡,计算结果有效地捕捉到了这一流动现象,并通过分析瞬时流线的计算结果揭示出尾缘涡的涡脱落是引起该高频振荡的主要原因.此外,通过数值计算分析了中低雷诺数下减缩频率对升力系数迟滞曲线的影响,从结果中看到,随着减缩频率的增加,翼型的失速攻角值会增加,升力系数峰值会增加;当减缩频率减小时,升力系数曲线中的高频振荡频率会增加.通过进一步计算分析知,升力曲线中高频振荡的产生不仅取决于减缩频率,还取决于流动的雷诺数,只有在中低雷诺数、较小的减缩频率下翼型动态失速的升力迟滞曲线中才有可能出现高频振荡.     

不可压缩势流绕多段翼型流动的快速迭代解法

利用保角变换计算不可压缩势流绕翼型的流动,具有计算速度快和要求计算机内存小的优点,但主要用于单段翼型问题,而对于多段翼型,尤其是三段和三段以上的翼型则很难应用。本文在单段翼型的保角变换的基础上,提出了一种利用迭代算法计算任意多段翼型的不可压缩势流的绕流问题的方法。本文方法和目前常用的边界元法相比,具有计算速度快的优点。本方法还可进一步推广到考虑无粘流/边界层的迭代计算中。     

冰形表面粗糙度对翼型的失速特性影响分析

结冰翼型气动影响分析是飞机结冰问题重要的研究内容之一,可帮助进行飞机结冰安全分析和容冰设计.为填补现阶段关于冰形表面粗糙度对翼型气动影响的研究空缺,以冰形的粗糙度参数(冰形表面粗糙高度和冰形表面粗糙分散度)为依据,提出了一种冰形粗糙度叠加方法;并通过有限体积法求解空间离散的RANS方程,结合S-A湍流模型,分别针对流线...     

基于Kriging模型的翼型气动性能优化设计

将基于Kriging模型的近似技术引入气动优化设计。通过计算两个测试函数的全局极值以及翼型几何形状的重构,比较了基于Kriging模型的优化方法与传统的梯度类优化方法、标准遗传算法的特点。结果表明基于近似模型的优化方法不仅全局性良好,而且有效节约了计算资源。用此方法进行了翼型的气动优化设计,目标函数为来流马赫数为0.7、攻角为3°时升阻比最大。设计结果表明,与参考翼型相比,优化翼型升阻比提高了81%,而翼型截面积减少不到1%,优化效果非常明显。     

可压缩性粘流绕翼型的摩阻计算

本文研究了一种可压缩性粘流绕翼型的摩擦阻力计算方法。这种方法包括用修正了的Thwaites~[1]法计算头部层流边界层以及按Green~[4]法计算湍流边界层。按照上述方法计算了在不同雷诺数和不同马赫数下两个翼型剖面的摩阻系数。经与文献[5]比较,吻合较好。     

神经网络在高空长航时无人机翼型多点优化中的应用

NSGAII算法在翼型多点设计中有着广泛的应用价值,然而其巨大的计算资源和计算成本限制了它的使用。为了解决这个问题,本文引入具有较强非线性映射能力的神经网络代理模型,采用实验设计结合BP法训练神经网络响应面来代替N-S方程求解翼型的性能。在实验点的数值模拟中,为了进一步节省计算资源,提高计算效率,采用网格的变形代替网格的重新划分,使得计算网格的更新速度提高了约50%。在翼型的参数化过程中,采用改进的PARSEC方法,用较少的参数实现了翼型的精确控制。为了增强神经网络的泛化能力,采用12-7-4-3-1的隐层结构。对NLF1015翼型的多点优化算例表明,此方法不仅显著降低了整个优化过程的计算量,而且对翼型的气动性能预测也具有较高的可信度,在高空长航时无人机的翼型设计中具有一定的潜力。