用全位势方程计算细长翼-身组合体的可压缩绕流

用Ｈ－Ｏ网格、守恒型全位势方程、差分和隐式近似因式分解迭代算法计算绕大后掠细长翼－身组合体的可压缩流动。对ＡＧＡＲＤ－Ｂ翼－身组合体的计算表明,计算的自由流Ｍａｃｈ数可从亚音速直到低超音速。算出的升力与力矩特性与实验数据符合较好     

攻角直到30°的超音速翼身组合体法向力和俯仰力矩的计算

本文基于“冲击流比拟”和“吸力比拟”提出了一种计算超音速翼身组合体的非线性法向力和俯仰力矩的方法。大量的计算与实验结果比较表明,法向力偏差约为±5％,压力中心偏差约为±1％弹长。本方法使用简便,可适用于α_∞≤30°的锐缘小展弦比、极小展弦比无滚转的十字形翼身组合体。     

翼—身—尾组合体绕流的Euler方程计算

本文采用重迭网格技术和Ｅｕｌｅｒ方程计算翼－身－尾组合体绕流。对翼－身与尾－身部分采用各自的Ｈ－Ｏ型网格。Ｅｕｌｅｒ方程求解采用Ｊａｍｅｓｏｎ有有限体积法，即中心差分挖和显式Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进。采用前后区交替迭代使前后两区通过重迭区交换信息。本文用ＮＡＣＡＴＮ４０４１翼－身－尾模型为例，计算的空气动力特性与实验符合较好。     

用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。     

DLR-F6/FX2B翼身组合体构型高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6和DLR-F6_FX2B 2种翼身组合体构型的高阶精度数值模拟,计算外形来自AIAA第三届阻力预测研讨会。主要目的是确认WCNS模拟跨声速典型运输机构型和预测局部构型变化引起的气动特性变化量的能力。在固定升力系数条件下,采用粗、中、细3套网格开展了网格收敛性研究,从气动力系数、压力系数分布、表面流态等方面研究了网格规模对DLR-F6和DLR-F6_FX2B翼身组合体数值模拟结果的影响;采用中等网格开展了来流迎角对2种翼身组合体气动特性的影响研究。通过与National Transonic Facility(NTF)的试验结果和CFL3D的计算结果对比,表明采用高阶精度计算方法得到了网格收敛的数值模拟结果,较好地模拟了DLR-F6翼身组合体局部修型引起的微小气动特性变化和翼身结合部流动特性的差异。     

翼-身组合体绕流的Euler方程数值模拟

用代数方法生成翼-身组合体H-O型网格,并用有限体积法研制出翼-身组合体绕流三维Euler方程计算程序。该方法的特点是改进了机身与机翼表面网格点分布,机翼后缘有后掠时也能保证后缘与网格线一致。程序除能按常规提供横流截面展向压强分布外,还能提供弦向压强分布。对NASA TND-712的翼-身组合体模型的计算结果与实验符合很好。     

跨音速三维机翼及翼身组合体非守恒型全位势流AF3算法

 本文对跨音速非守恒型全位势流,用Baker的快速有限差分法,改变其网络拓扑,研制出分析机翼及翼身组合体跨音速绕流的实用程序。其中包括:对“C-H”型贴体坐标网格提出AF3算法的具体实施方案取得了快速收敛效果;采用张量形任意曲线坐标上的非守恒全位势方程,使得程序中跨音速求解部分既能用于分析单独机翼又能用于翼身组合体跨音速绕流计算。     

大迎角体涡,翼涡干扰的Euler方程数值模拟

本文用Ｅｕｌｅｒ方程的简化区算法计算翼－身组合体绕流。当大迎角时，计算能自动捕捉到大后掠前缘的分离涡，机射的分离涡由强加Ｋｕｔｔａ条件生成，形成体涡与翼涡的相互干拢。本文对加长前机身的ＡＧＡＲＤ－Ｂ翼－射组合体模型，给出机射有分离和无分离时，横流截面的速度矢量分布、法向力和压力中心，分析体涡的影响。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法．在背景网格建立和三维初始推进面生成的前提下，利用推进阵面法生成了翼身组合体的三维空间网格，并应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

翼－身－尾组合体绕流的Euler方程计算

本文采用重迭网格技术和ＥＪｅｒ方程计算翼－身－尾组合体绕流。对算身与尾。身部分采用各自的Ｈ。Ｏ型网格，Ｅ＊ｅ／方程求解采用１ａｍｅｓｏｎ的有限体积法，即中心差分近似和显式Ｒｕｎｇｅ·Ｋｕｔｔａ时间推进。采用前后区交替迭代使前后两区通过重迭区交换信息。本文用ＮＡＣＡＴＮ４０４１翼－身，尾模型为例，计算的空气动力特性与实验符合较好。     

有翼飞行器大攻角分离涡模拟

基于势流方程,采用涡格法和三维离散涡相结合,模拟具有升力面分离涡的大攻角弹冀、弹体、尾翼组合体绕流。该方法适用于高亚音速、小攻角或大攻角有翼飞行器的压力分布、法向力和压力中心等气动特性的数值计算,并具有显示三维涡流场的能力。     

鸭式布局飞行器的翼体摇滚特性风洞试验

为研究鸭式布局飞行器摇滚特性,设计了一种包括鸭翼、脊型前体、边条翼、主翼和垂尾的模型,进行了自由滚转、扰动滚转、静动态测力和烟线流场显示多种技术手段相结合的风洞试验。通过自由滚转和扰动滚转试验得到了该模型翼体摇滚的时间历程,静态测力和动导数测定验证了非极限环运动形式摇滚的发生。结果表明该鸭式布局模型摇滚不仅同侧存在多个摇滚平衡点,而且在临界俯仰角,摇滚过程中可能出现从一摇滚平衡点跳动至同侧另一摇滚平衡点的突变。通过流场显示技术得到该鸭式布局模型复杂流场的基本形态分布,并对滚转角为0°时的全机涡系干扰和摇滚形成机理进行了简要分析。     

文氏效应在叶片根部马蹄涡控制中的应用

为研究文氏管在叶片根部马蹄涡控制作用的应用,针对简化了的叶片-平板结构在叶片上游平板上的某位置开狭缝,在平板的另一侧构造一文氏管(文氏管的喉道位于狭缝处),引射叶片根部的流体,进而控制马蹄涡.数值模拟了不同狭缝位置的湍流流动;结果表明了改变狭缝相对于叶片的位置能很好地控制马蹄涡,甚至能完全消除叶片前马蹄涡.提出的方法可以利用外势流的能量而无需额外消耗能量,为一种被动的控制手段,有很大的实际应用前景.利用文氏效应的控制方法简单实用,因此工程中采用该法是不错的一种选择.      

旋涡发生器对叶片根部马蹄涡的影响

针对简化了的叶片-平板结构,在叶片上游平板上的某位置布置旋涡发生器影响马蹄涡.数值模拟了不同高度旋涡发生器作用下的角区湍流流动.结果表明:任何一种旋涡发生器都能不同程度地削弱马蹄涡系,高度为2mm的旋涡发生器能使这种削弱作用达到最大;另外,随着高度的增加,叶片两侧的涡腿有向叶片聚拢的趋势.提出的方法无需额外消耗能量,为一种被动的控制手段,有很大的实际应用前景.布置旋涡发生器的控制方法简单实用,工程中采用该法是不错的一种选择.      

COMPUTATION OF COMPRESSIBLE FLOW PAST SLENDER WING－BODY USING FULL－POTENTIAL EQUATION

ＣＯＭＰＵＴＡＴＩＯＮＯＦＣＯＭＰＲＥＳＳＩＢＬＥＦＬＯＷＰＡＳＴＳＬＥＮＤＥＲＷＩＮＧ－ＢＯＤＹＵＳＩＮＧＦＵＬＬ－ＰＯＴＥＮＴＩＡＬＥＱＵＡＴＩＯＮＨｕａｎｇＭｉｎｇｋｅ（ＤｅｐａｒｔｍｅｎｔｏｆＡｅｒｏｄｙｎａｍｉｃｓ，ＮａｎｊｉｎｇＵｎｉｖｅ...     

前体涡诱导双极限环摇滚流动特性的实验研究

通过摇滚/PIV/压力同步测量实验,对翼身组合体前体涡诱导的双极限环摇滚过程中流动特性及演化规律进行了系统的研究,并分析了前体涡诱导翼-身组合体双极限环摇滚的流动机理。实验结果表明,前体涡与机翼翼面流动的相互作用使模型在正负滚转相位处分别出现极限环摇滚运动;正负滚转相位过渡是模型运动惯性与气动力共同作用的结果。     

旋翼桨尖涡模型及其在自由尾迹分析中的影响

 建立了一个旋翼桨尖涡初始涡核半径计算模型,用于消除旋翼自由尾迹分析中桨尖涡涡核半径依靠经验设定的现状。该方法利用涡量和涡矩守恒关系,考虑旋转效应的影响,并修正Vatistas涡量分布函数,得到了一个关于桨尖涡初始涡核半径与桨叶附着环量分布之间的非线性关系。将新建立的初始涡核模型同先前建立的涡核扩散与拉伸模型以微分方程初值问题的形式组成一个完整的桨尖涡模型,并将其应用于自由尾迹分析方法中,分别对两副模型旋翼的悬停尾迹进行分析。计算得到的桨尖涡初始涡核半径和展向位置与实验测量值的对比表明本文模型具有较高的准确性。此外,讨论了涡核湍流动量扩散效应在自由尾迹分析中的影响,分析表明该效应对自由尾迹分析的收敛性有重要影响。     

Transonic wing stall of a blended flying wing common research model based on DDES method

Numerical simulation of wing stall of a blended flying wing configuration at transonic speed was conducted using both delayed detached eddy simulation(DDES) and unsteady Reynolds-averaged Navier-Stokes(URANS) equations methods based on the shear stress transport(SST) turbulence model for a free-stream Mach number 0.9 and a Reynolds number 9.6 × 10~6. A joint time step/grid density study is performed based on power spectrum density(PSD) analysis of the frequency content of forces or moments, and medium mesh and the normalized time scale0.010 were suggested for this simulation. The simulation results show that the DDES methods perform more precisely than the URANS method and the aerodynamic coefficient results from DDES method compare very well with the experiment data. The angle of attack of nonlinear vortex lift and abrupt wing stall of DDES results compare well with the experimental data. The flow structure of the DDES computation shows that the wing stall is caused mainly by the leeward vortex breakdown which occurred at x/x_(cr)= 0.6 at angle of attack of 14°. The DDES methods show advantage in the simulation problem with separation flow. The computed result shows that a shock/vortex interaction is responsible for the wing stall caused by the vortex breakdown. The balance of the vortex strength and axial flow, and the shock strength, is examined to provide an explanation of the sensitivity of the breakdown location. Wing body thickness has a great influence on shock and shock/vortex interactions, which can make a significant difference to the vortex breakdown behavior and stall characteristic of the blended flying wing configuration.