一种跨音速翼型和机翼设计方法的新进展

 讨论一种基于N-S方程、欧拉和速势方程加混合附面层修正的跨音速流动准确控制方程逆解法，及其在跨音速翼型和机翼设计中的应用。分析了两种紊流模型对分离的适应性，探讨了N-S和欧拉方法的阻力估算精度问题，列举了若干超临界和层流翼型及机翼的设计实例以及软件验证。     

跨音速翼型气动优化设计方法

 提出了一种有效的跨音速翼型气动优化设计方法。翼型的流场解由欧拉方程的数值解提供。带约束条件的优化计算分别采用了直接法和间接法两种优化算法。算例结果表明,本方法提供了一种跨音速翼型改型设计及新翼型设计的有效工具。     

基于粒子群算法的翼型优化设计

采用粒子群算法（PSO）对层流翼型进行了以提高升阻比为目标的优化设计。翼型的设计达到了设计要求,优化设计后的翼型气动特性也有显著地改善,这表明了粒子群算法应用于翼型气动优化设计的可行性。在优化设计的过程中,粒子采用递减惯性权重,以加强粒子初期的全局搜索能力与后期的局部搜索能力。翼型由解析函数线性叠加法表示,目标函数和粒子的适应度由基于二维欧拉方程的流场数值解来提供。     

基于气动数值模拟的翼型反设计方法

将计算流体动力学（CFD）与反设计技术相结合,通过数值求解欧拉方程,对翼型绕流流场作出数值模拟,再用几何和流动控制方程,反复迭代求得满足给定流场的翼型。以NACA0012为初始翼型,RAE-2822为目标翼型,选取两种工况,都取得了满意的结果。     

Roe格式在水滴欧拉方程数值求解中的应用

目前普遍采用欧拉法求解水滴运动控制方程以模拟结冰表面水滴撞击特性。水滴欧拉方程近似没有压力项的空气欧拉方程,但已应用于空气欧拉方程的空间离散格式多数不能直接应用于水滴欧拉方程,有研究者将经典Roe格式应用于水滴欧拉方程时出现了奇点问题。通过对水滴欧拉方程进行形式改造,基于中心格式验证改造方法的有效性,并应用Roe格式结合线性重构与限制器函数实现对水滴欧拉方程二阶精度的求解。将NACA23012翼型计算结果与FENSAP-ICE软件进行对比,两者水滴收集系数分布基本一致。     

基于遗传算法的翼型气动优化设计

采用遗传算法进行跨声速翼型的反设计与阻力和升阻比的优化设计。翼型的反设计达到了设计要求,优化设计后的翼型其气动特性也有显著的改善,这表明了遗传算法应用于翼型气动优化设计的可行性。在优化设计的过程中,翼型由解析函数线形叠加法表示,目标函数和个体的适应值由二维欧拉方程的流场解来提供。     

基于欧拉方程的机翼设计方法及应用

以Takanashi提出的三维机翼设计理论为基础，研究与发展了一个基于欧方程和“正－反迭代、余量修正原理”的机翼设计方法。用改进的无限插值方法生成绕机翼的O－O贴体网格，采用三维欧拉方程作为充动分析计算的基本方程，该设计方法已用于某无人机机翼和一个超临界机翼设计，设计结果达到了预期的目标。     

基于动态边界控制的气动优化设计方法研究

本文根据控制变形物体边界变化的思想,建立了一种基于动态边界控制的气动优化设计方法,此方法通过在求解非定常流动方程过程中计算设计目标函数随设计变量变化的敏感性导数,这样在优化过程中不再需要求解定常流动方程,或共轭方程,从而使得优化过程得到简化并且便于在不同条件下应用.文中以二维翼型为例,用非定常欧拉方程为基本流动控制方程,针对两种优化设计问题,开展了翼型的优化设计,取得了良好的优化设计结果.     

基于N-S方程的直升机旋翼翼型反设计方法

为了在实际旋翼翼型设计中考虑粘性影响,采用N-S方程为主控方程,建立了旋翼翼型流场求解方法。使用以Poisson方程为控制方程的网格生成程序生成围绕旋翼翼型的N-S粘性贴体网格,该网格生成程序能在反设计迭代过程中自动更新过渡翼型的粘性网格。在翼型流场计算及网格生成基础上,采用了MGM方程作为翼型反设计方程,建立了一套能够考虑粘性影响的直升机旋翼翼型的反设计方法。应用该方法,分别对二维NACA系列翼型、OA系列翼型、超临界翼型进行了反设计分析,获得了满足要求的二维旋翼翼型,并与目标值(压力分布、升力系数、阻力系数、力矩系数)吻合良好。     

欧拉方程多段翼型绕流计算的嵌套网格方法

本文采用嵌合体埋入技术，将主段翼型和襟翼分成两子域，对多段翼型的复杂构型进行了分区迭代求解，在有限体体积离散的基础上，求解欧拉方程计算各子域流场的解，通过各流场的耦合迭代到多段翼型绕流的解，典型的带３０％襟翼ＧＡ（Ｗ）－１翼型算例表明，该计算方法稳定，能较地模拟多段翼型中主段翼型与襟翼交接处的复杂流动现象计算结果与实验值有较好的一致性。     

旋成面叶栅气动反命题的新解法(Ⅱ)

本文利用映象平面ξ-η的特点, 将旋成面叶栅气动正命题的微分-积分解法与叶型修正的“喷气模型”结合, 建立了一种适用于旋成面叶栅气动反命题的新解法。叶片表面给定的气动参数分布作为反命题的边界条件, 在求解过程中直接得到满足。通过对一例实际透平叶型的计算验证了本文解法的可靠性和实用性。      

机翼跨声速绕流的Galerkin有限元解法

本文给出了翼型和机翼跨声速绕流的Galerkin带权余数有限元解法。利用一系列适当形状的元素将计算区域离散,在元素内采用线性插值函数。对于超临界流动,文中采用了适合于有限元法的超声速上风技术,具有捕捉激波的能力。将这种技术应用于线松弛迭代解法中,在计算钝头翼型和机翼绕流时,都获得了成功。     

后缘连续变弯度对跨声速翼型气动特性的影响

针对后缘连续变弯度对跨声速翼型气动特性的影响进行了研究。首先不考虑翼型后缘连续变弯度,基于搭建的优化设计系统对跨声速翼型进行气动减阻优化设计,通过添加不同的约束优化得到两种跨声速翼型:无激波翼型和超临界翼型。然后在这两种翼型的基础上,以后缘偏转角度为设计变量、以阻力系数最小为目标,针对不同的升力系数分别进行优化设计,并根据优化结果深入分析后缘连续变弯度对这两种翼型极曲线特性的影响机理。优化结果表明:无激波翼型与超临界翼型相比,其设计点处的气动特性较好,但鲁棒性较差;升力系数小于设计升力系数时,应用后缘连续变弯度后,无激波翼型的极曲线特性明显提高,减阻最高达到3.9%,而超临界翼型的极曲线特性提高不明显;升力系数大于设计升力系数时,应用后缘连续变弯度后,无激波翼型和超临界翼型的极曲线特性都明显提高,减阻分别达到2.4%~18.1%和1.7%~13.2%。     

现代自然层流翼型的设计方法

本文介绍一种设计跨声速自然层流翼型的计算流体力学(CFD)方法。本方法采用“正反迭代、余量修正”设计原理,通过将跨声速翼型设计软件NPU-TD2D中的反设计程序进行改进,并与含有层、湍流混合边界层修正的跨声速层流翼型计算程序DLRBGKWALZ耦合,实现了在跨声速粘性流动条件下直接设计层流翼型。亚、超临界的设计实例和风洞验证表明,本方法可以在几个设计迭代内设计出压力分布、转捩位置及气动参数均准确收敛于设计目标的新翼型,是一种设计现代自然层流翼型的有效而实用的CFD方法。     

Flow characteristics around airfoils near transonic buffet onset conditions

In transonic flow, buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation. The phenomenon is common in transonic flow, and it has serious impact on the structural strength and fatigue life of aircraft. In this paper, three typical airfoils: the supercritical OAT15A, the high-speed symmetrical NACA64A010, and the thin, transonic/supersonic NACA64A204 are selected as the research objects. The flow fields of these airfoils under pre-buffet and buffet onset conditions are simulated by Unsteady Reynolds Averaged Navier-Stokes (URANS) method, and the mode analysis of numerical results is carried out by Dynamic Mode Decomposition (DMD). Qualitative and quantitative analysis of the shock wave motion, shock wave intensity, shock foot bubble and trailing edge separation, and pressure coefficient fluctuation were performed to attain deep insight of transonic buffet flow features of different airfoils near buffet onset conditions. The results of DMD analysis show that the energy proportion of the steady mode of these airfoils decreases dramatically when approaching the buffet onset angle of attack, while the growth rate of the primary mode increases inversely. It was found that at the onset of buffet, there exist different degrees of merging behavior between shock foot bubble and trailing edge separation during one buffet cycle, and the instability of shock wave and separation induced shear layer are closely related to the merging behavior.     

Non-intrusive reduced-order model for predicting transonic flow with varying geometries

A Non-Intrusive Reduced-Order Model (NIROM) based on Proper Orthogonal Decomposition (POD) has been proposed for predicting the flow fields of transonic airfoils with geometry parameters. To provide a better reduced-order subspace to approximate the real flow field, a domain decomposition method has been used to separate the hard-to-predict regions from the full field and POD has been adopted in the regions individually. An Artificial Neural Network (ANN) has replaced the Radial Basis Function (RBF) to interpolate the coefficients of the POD modes, aiming at improving the approximation accuracy of the NIROM for non-samples. When predicting the flow fields of transonic airfoils, the proposed NIROM has demonstrated a high performance.     

跨声速翼型多目标优化设计方法

本文通过采用多目标分级优化方法并基于求解全速势方程的跨声速粘流翼型计算方法，研究发展了一种多参数、多约束和多目标的跨声速翼型数值优化设计方法，应用该方法可以从普通低速翼型和超临界民办型出发通过多目标优化后得到在跨声速区的多个马赫数下阻力系数最小化的翼型几何外形，设计实践表明，该方法具有收敛快，调用目标函数次数少等优点。     

一种求解核函数方程的新方法—辅助核法

本文提出了求解核函数方程的辅助核法。此法克服了偶极子栅格法对网格分布的限制,具有更大的适用范围。数值计算表明,此法对求解绕翼型的跨声速激波流动是很有利的。     

解跨音速升力翼型的有限元法

 在应用解全速位方程的最小压强积分有限元法求解绕升力翼型的跨音速流动时,将不可压流中求解绕升力翼型的耦合单位环量流动和无环量流动的解法推广到可压流中。为了确定环量,本文所用Kutta条件是:在后缘处,气流流向平行于后缘角二等分线。因有限元法对网格无正交性要求,因而可在椭圆变换前后进行剪切和延伸变换。这种网格生成法易于构成适用于复杂形状的有限元网格。通过计算并将其结果与文献中的数据比较,表明这种方法应用方便且有较快的计算速度和较高的计算精度。     

计算机翼跨音速绕流的加罚Galerkin有限元法

1.引言 全速位Galerkin有限元法是跨音速有限元法中重要方法之一。处理超临界流时,为了自动捕捉到激波,通常采用上风技术,但因此导致了质量守恒条件得不到完全满足。文献[1]从Bateman变分原理出发,应用加罚函数方法,找出了强制满足质量守恒条件在Galerkin有限元方程中所应加入的附加项。本文将文献[1]的方法应用到超临界跨音速绕机翼流动,克服了一般Galerkin有限元法中由于采用上风技术而出现的激波不够陡和位