一种新型的直角网格生成方法及应用研究

提出了一种基于背景网格的直角网格新型生成方法。背景网格是用来提供网格加密的信息．通过网格尺度与所得计算参数的比较，定量判断确定网格是否需要加密；而在求解Euler方程计算中，采用的是Jameson的有限体积法，它可以适合任意形状的网格单元，以及四步Runge—Kutta时间推进。本文对NACA0012翼型和复杂多段翼型等问题进行了数值实验。结果表明，这种网格生成方法易于推广到三维情况中去，具有网格生成时间短．收敛速度快，易于处理复杂边界的优点。     

直升机旋翼计算流体力学的研究进展

依据直升机旋翼计算流体力学(CFD)的发展，分别介绍了小扰动位势方程、全位势方程、Euler方程和Navier—Stokes方程在旋翼流场计算中的研究现状和发展趋势，着重从方程离散、数值算法、网格生成、计算时间和计算精度等方面分析了不同旋翼CFD方法的特点，指出了旋翼的尾迹在旋翼流场计算中的重要性，并针对位势方程和Euler／N—S方程分别讨论了求解的边界条件。最后，对旋翼CFD的发展提出了几点展望。     

混合网格在求解三维N—S方程中的应用

用三棱柱／四面体组成的混合网格求解三维Navier-Stokes方程，在物面附近采用三陵柱网格，而其他地方则采用四面体网格，与单一类型网格相比，混合网格易于调节粘性层内网点数目，并能减少存储量，提高解的效率，Navier-Stokes方程的求解采用有限体积法，对人工耗散项进行了改进，使得求解的结果更接近真实的物理解，本文对上述方法开展了数值实验研究。结果表明，用三棱柱／四面体组成的混合网格求解N－S方程是非常有效的。     

围绕旋翼飞行器的三维结构化运动嵌套网格生成方法

考虑到旋翼飞行器的旋翼各种运动和机身复杂外形特点,建立了一种鲁棒性高的三维结构化嵌套网格生成方法。采用三维Poisson方程对初始网格进行迭代,获得了高质量的旋翼桨叶三维结构化网格。以此方法为基础,以复杂外形几何体表面网格为初始网格,建立了一种抛物型法向外推方法。推进过程中,使用经抛物化处理的新Poisson方程进行光顺迭代,源项由Hilgenstock法确定,有效克服了椭圆型方程生成网格时边界点较难调整的困难。通过网格最小正交度的检测,表明了该抛物型方法的先进性。进一步将高效的挖洞法("Top map"法)和贡献单元搜寻法("Inverse map"法)引入上述三维结构化网格生成中,建立了围绕旋翼飞行器的运动嵌套网格生成方法,采用RANS方程分别对Caradonna-Tung旋翼和V-22倾转旋翼机流场进行计算流体力学(Computational fluid dynamics,CFD)计算,表明该网格生成方法能够满足CFD计算的需求。     

模型燃烧室中两相化学反应流动的大涡模拟

对带V型火焰稳定器的模型燃烧室中气液两相化学反应流动进行大涡模拟(LES)，研究中采用了欧拉一拉格朗日方法。利用大涡模拟来研究气相流场，k方程亚网格尺度模型来模拟亚网格涡粘性。化学反应速率由亚网格EBU燃烧模型来确定，两步亚网格紊流燃烧模型被用来确定CO浓度，NO亚网格污染物生成模型被用来评估NO生成速率，辐射热量传递采用了热通量辐射模型。并在交错网格体系下气相采用SIMPLE算法和混合差分格式求解，液相在拉格朗日坐标系下采用PSIC算法对其进行求解。通过大涡模拟，能得到气相和液相之间的质量、动量和能量交换。通过计算值和实验值的比较可知，出口温度和污染物成分浓度实验结果与计算结果吻合较好。因此采用大涡模拟方法来模拟两相化学反应流动和污染物生成过程是可行的。     

粘性非结构网格生成方法及绕流计算

运用非结构网格进行计算时，需要在物体附面层区域内生成高展弦比的三角形网格。本文给出了一种基于物理学中弹簧原理的粘性非结构网格生成新方法－－挤压法，它能有效地实现粘性区域内位伸网格到附面层外一般非拉伸网格的光滑过渡。然后在非结构粘性网格上对Ｎ－Ｓ方程运用有限体积法进行了数值计算，ＮＡＣＡ－００１２翼型以及平板层流的数值试验结果较为满意。     

计入桨叶运动的旋翼CFD网格设计技术

提出了一个考虑桨叶实际的挥舞和变距运动的旋翼CFD网格生成方法。该方法首先采用改进的几何法在桨叶表面生成初始网格，通过求解Poisson方程得出贴体正交网格，然后结合嵌套网格和动态网格的思想，给出了一种运动嵌套网格来计入桨叶的挥舞和变距运动．并针对嵌套网格中关键的贡献单元搜寻问题，提出了一种新的“伪贡献单元搜寻法(PSDE)”。该搜寻方法不仅可以适用于在结构网格中搜索贡献单元，而且不作修改即可应用于非结构网格，还可克服“InverseMap”法在结构网格中需要单独处理的奇异线(面)问题。     

三维Euler方程的自适应八叉树结构直角网格算法

给出了三维自适应八叉树结构直角网格的生成方法，采用了以几何外形为基础的网格加细方法，并解决了网格生成过程中的若干难点，使网格生成更简便迅速，产生的网格更适合于Ｅｕｌｅｒ方程计算。在计算中采用以中心差分为基础的Ｊａｍｅｓｏｎ的有限体积法和以面为基础的通量计算方法，减少了工作量。在物面边界处采用通量分解的方法实现边界条件，简单易行。本文对前机身、单个及多个外挂等问题进行了数值实验。结果表明，计算与实验符合较好。     

双流混合边界的正交贴体网格生成

通过数值实验表明，用泊松方程变换形式生成的贴体坐标网格，其边界点的β值与源项呈连续单调的函数关系。根据这种关系，可以用方程迭代求根方法中的增值寻根法和对分法，使网格边界正交化。由薄壁结构相隔的两种流体混合边界正交网格生成的实例表明，这种方法对于生成适应于红外抑制器全流场计算的网格是很有效的。     

任意回转面叶栅C网格的数值生成

首先尝试采用任意曲面椭圆方程，数值生成任意回转面上叶栅Ｃ网格，经过对边界条件的仔细处理，生成的网格可保证在周期外边界和重合内边界上网格线导数连续，在叶片表面网格线相互垂直，提出了结合初始网格线分布，控制环绕叶片网格线分布的方法。     

Aerodynamic and Acoustic Optimization for Multi-element Airfoils

The paper is to integrate aerodynamic and aero-acoustic optimizatiom design of high lift devices,especially for two-element airfoils with slat.Aerodynamic analysis on flow field utilizes a high-order,high-resolution spatial differential method for large eddy simulation(LES),which can guarantee accuracy and efficiency.The aeroacoustic analysis for noise level is calculated with Ffowcs Williams-Hawkings(FW-H)integration formula.Fidelity of calculation is verified by standard models.Method of streamline-based Euler simulation(MSES)is used to obtain the aerodynamic characters.Based on the confirmation of numerical methods,detailed research has been conducted for the leading edge slat on multi-element airfoils.Various slot parameter influences on noise are analyzed.The results of the slot optimization parameters can be used in multi-element airfoil design.     

Euler方程无网格算法及布点技术

研究了Bitina的显式无网格算法。在该算法的基础上,给出了Euler方程无网格离散形式,运用RungeKutta显式时间推进格式推进求解。文中的耗散项采用了与非结构网格上类似的守恒型耗散算子,边界条件的处理借鉴了结构化网格处理技术。此外,本文还描述了一种区域离散布点方法,研究了点云生成的选点准则,并成功地数值模拟了二维翼型的典型绕流。     

基于结构网格点的Euler方程无网格自适应求解算法

提出了一种自适应无网格Euler方程求解算法。构造了基于无网格点云结构的布点加密技术，实现了借助压力梯度变化准则的流场局部加密；借助流场结构化剖分技术，给出初始无网格节点分布，并采用四步显式Runge—Kutta时间格式推进计算，求解了绕翼型的典型流动问题。数值算例表明，本文提出的方法能有效地提高流场的分辨率，如捕捉激波等流动特征清晰明了。     

尖前缘任意翼型的超音速和高超音速俯仰安定性导数

本文推广了文[1]的欧拉方程摄动法,讨论了尖前缘任意翼型俯仰安定性导数的计算方法,并应用该方法计算了双圆弧翼型的俯仰安定性导数,其数值结果与文[2]的超音速线化位流方程有限差分法的数值结果比较,更接近于实验值。在激波附体的条件下,本文的方法可用于计算大迎角尖前缘任意翼型的超音速和高超音速俯仰安定性导数。     

自适应计算网格的生成和优化方法

本文提出一种简单通用的自适应计算网格生成方法。配合文[1]的网格优化方法,能在生成自适应网格的同时,保证网格具有较好的光滑性和正交性。本文还阐述了单元加权函数W的取值和边界点的确定对自适应网格性质的影响。作为应用实例,本文将自适应网格用于二维无粘跨音速叶栅流场计算。结果表明,采用自适应网格能明显改善流场的计算精度。     

Highly-Efficient Aerodynamic Optimal Design of Rotor Airfoil Using Viscous Adjoint Method

In order to overcome the efficiency problem of the conventional gradient-based optimal design method,a highly-efficient viscous adjoint-based RANS equations method is applied to the aerodynamic optimal design of hovering rotor airfoil.The C-shaped body-fitted mesh is firstly automatically generated around the airfoil by solving the Poisson equations,and the Navier-Stokes(N-S)equations combined with Spalart-Allmaras(S-A)one-equation turbulence model are used as the governing equations to acquire the reliable flowfield variables.Then,according to multi-constrained characteristics of the optimization of high lift/drag ratio for hovering rotor airfoil,its corresponding adjoint equations,boundary conditions and gradient expressions are newly derived.On these bases,two representative rotor airfoils,NACA0012 airfoil and SC1095 airfoil,are selected as numerical examples to optimize their synthesized aerodynamic characteristics about lift/drag ratio in hover,and better aerodynamic performance of optimal airfoils are obtained compared with the baseline.Furthermore,the new designed rotor with the optimized rotor airfoil has better hover aerodynamic characteristics compared with the baseline rotor.In contrast to the baseline airfoils optimized by the finite difference method,it is demonstrated that the adjoint optimal algorithm itself is practical and highly-efficient for the aerodynamic optimization of hover rotor airfoil.     

一种基于块结构的非结构网格快速生成方法

以阵面推进法为基础，提出了一种用预制模块的方法实现在平面域内生成非结构网格的新方法，即模块方法。模块是由一组三角形构成，文中以三角形模块为例，介绍了网格模块的预制方法。网格的生成是根据空间步长的背景信息，从预存的网格模块中选取当地合适的模块，直接复制或变形到计算区域实现的。由于生成过程中每一步引入的是预制的网格模块而不是传统的单个三角形单元，因此可期望新方法要比传统方法快得多。算例显示了新方法的快速网格生成过程，该方法具有推广应用前景。     

多段翼型绕流压缩性修正的实用方法

本文对前人提出的各种二维流动压缩性修正方法进行了大量分析、计算和比较,提出了多段翼型复杂混合流场的压缩性修正方法,把翼型绕流划分为外部流场和内部流场(翼段间缝隙)。外部流场采用Nrstrud方法进行修正,内部流场采用Lieblein方法修正,并对这两种修正方法进行了局部改进使其更加应用方便。本文翼型绕流压缩性修正方法的应用,使现有的多段翼型失速特性计算程序有了更为广泛的应用范围。     

以风轮气动性能为目标的风力机翼型优化设计

提出了一种优化风力机翼型气动特性的设计方法。在现有风力机翼型的基础上,通过改变其最大厚度位置至后缘部分的外形进行优化计算。采用N-S方程计算翼型的流场及气动特性,以提高风力机切向力同时减小其法向力为目的,用二维线性插值的方法进行多目标优化。以DU 93W 210翼型为例,针对风力机翼型气动性能本身的矛盾性,根据实际需求提出几种方案加以优化。结果表明,在不同的性能需求下,本文所提出的设计方法能取得较好的效果。     

几何尺度对GAW-1两段翼型升力特性的影响

在大型飞机和高速飞机的起降过程中增升装置空气动力性能的好坏对飞机有着举足轻重的影响.飞机增升装置通常采用多段翼型形式.在飞机设计阶段,一般是采用风洞实验的方法对多段翼型的气动性能进行评估.由于风洞的尺寸和动力所限,实验的模型和实验雷诺数都小于实际飞行情况.在数据的使用时,一般认为实验数据进行雷诺数修正即可.通过实验研究发现,即使雷诺数相同,几何尺度对两段翼型升力特性的影响也是非常显著的.建议在风洞实验中重视几何尺度对气动性能的影响.该研究结果和研究思路可供增升装置的设计和实验研究人员参考使用.