二维涡格法的超收敛性

本文提出并研究了二维涡格法的超收敛性。在离散化误差分析中,与薄翼理论的二维平板精确解进行比较,根据Chebychev多项式理论,导出了涡格法中具有超收敛性的离散化格式。并证明了这种离散化格式对二维抛物弧板和三次曲线弧板等问题也具有超收敛性。     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

Multhopp网格在弯曲涡层计算中的超收敛性问题

本文从理论上分析了Ｍｕｌｔｈｏｐｐ网桥弯曲涡面积分方程求解中的数值方法和超收敛性问题;并通过算例验证,证实Ｍｕｌｔｈｏｐｐ网格应用于弯曲涡层时,的确同样具有远大于其它任何网格计算的超收敛速度。     

二维可压无粘流的自适应流量修正有限元解

从非定常Ｅｕｌｅｒ方程组出发,采用非结构三角形网格,研究二维无粘流动问题的自适应流量修正有限元解。将ＦＥＭ－ＦＣＴ与网格自适应相结合计算了ＮＡＣＡ００１２翼型的三个算例。计算结果与实验及其他计算结果比较,亚音速和超音速绕流符合很好;跨音速绕流情况存在一些差别。     

柱后旋涡运动离散涡方法数值模拟中的几个问题

本文简短评述用离散涡方法对圆柱绕流非定常旋涡运动进行数值模拟中的四个问题,即离散方法及涡面稳定性;边界层分离及新生涡位置;二次涡问题;尾迹涡中的环量减少问题等。文中还提出了尚需深入研究的课题。     

超声速空气介质自由扩散混合流涡的演化

在不同谐波扰动促发方式下,观察二维超声速空气/空气自由扩散混合流中涡的形成及演化.考虑含空气组分扩散效应的可压缩Navier-Stokes方程,对流项采用3阶迎风紧致格式离散,输运项采用6阶对称紧致格式离散,非定常时间推进采用3阶紧致存储显式Runge-Kutta方法.在无相差情形下,获得了大尺度基频涡结构的饱和、一次对并、二次对并、三涡对并等现象.在Ma_c=0.3低对流Mach数下,基频涡较饱满,但流向尺度较小.受空气真实气体特性影响,在Ma_c=0.8高对流Mach数下未发现涡/小激波结构.      

浸入边界法模拟小圆柱对主圆柱涡脱落的抑制

使用浸入边界法研究了小圆柱对主圆柱涡脱落的抑制.方法中使用非贴体笛儿尔网格,易于处理包含复杂边界的流动问题.采用离散附加力直接加入边界条件方法对虚单元进行重构,使边界条件在浸入边界上精确满足.使用隐式分步法解二维非定常不可压Navier-Stokes方程,通过速度和压力解耦提高计算效率.数值模拟单圆柱绕流及不同位置小圆柱和主圆柱的流动干扰,通过分析流场涡结构和升、阻力系数,得到小圆柱对主圆柱涡脱落的延迟和抑制作用.计算结果与已有实验结论和数值结果对比,计算误差不超过5%,说明浸入边界法可以简单有效地处理圆柱涡脱落抑制这类流动干扰问题.     

圆锥脉冲等离子体诱导流场特性分析

对等离子体诱导流场特性进行研究,有利于解决双稳态非对称分离涡带来的连续比例控制困难的问题。在封闭光学玻璃箱体内,应用介质阻挡放电等离子体对20°顶角圆锥附近静止大气进行了定常和脉冲循环控制,对等离子体诱导的圆锥截面绕流速度场进行了二维PIV测量,对定常控制和脉冲循环控制下最大绕流速度及最大轴向涡量进行了比较分析。实验结果表明：相对于定常控制模式,脉冲循环控制下沿垂直于圆锥截面对称面径线分布的时间平均切向速度和轴向涡量范围较广;在脉冲循环控制下,动量传递的主要表现在离散涡的形成而不是气流的加速。     

五种湍流涡粘模型在二维方柱绕流数值模拟中的对比研究

为研究不同雷诺时均Navier-Stokes(RANS)模型在求解钝体绕流场的差异,采用五种两方程湍流涡粘模型计算了二维方柱的绕流场(雷诺数Re=22 000),得到了不同湍流模型在不同网格离散方案下的方柱气动力与流场特性。结合以往文献数据,通过对比不同湍流模型计算结果的差异,揭示了不同湍流模型的特点。研究结果表明:网格离散方案对方柱绕流数值模拟结果有重要影响;RNGk-ε湍流模型的计算效率最高,SSTk-ω湍流模型的计算效率最低;Standard k-ε湍流模型的计算结果准确程度整体弱于其余湍流模型;RNGk-ε湍流模型、Realizable k-ε湍流模型与Standard k-ω湍流模型的计算结果大致相当,较接近大涡模拟结果;SSTk-ω湍流模型的模拟结果优于其余湍流模型,其尾流区速度场与试验结果吻合较好;两方程湍流模型计算的二维方柱绕流结果与试验结果以及大涡模拟结果相比,存在不可忽略的差异。     

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

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

一种高负荷叶型优化设计方法

采用3段3次贝塞尔曲线拟合叶型中弧线弯角分布,两段3次贝塞尔曲线拟合厚度分布,进行流线面造型,并结合优化算法进行寻优.S1流场评价工具为MISES软件,寻优算法为遗传算法（GA）.将传统遗传算法进行改进,引入并行模块,较大程度地提高了计算效率.实现了高负荷叶型的自动优化,探讨了高负荷叶型的设计准则,并将准则应用于超声叶型设计.结果表明:相比于原始叶型,优化后的叶型Ⅰ在总压损失系数略降低和出口气流角不变的条件下,可用攻角范围拓宽约2°,叶型Ⅱ的总压损失系数明显降低,可用攻角范围拓宽约0.6°,落后角降低约3.9°.      

外边界位置对数值计算翼型绕流的影响

用数值方法计算翼型绕流时,外边界位置一般只取在距翼型有限的距离。本文用不可压位流理论找出该问题的精确解,从而得到外边界位置对计算的升力、阻力和表面压力分布的影响。     

Unsteady oscillation of a deformable airfoil section in incompressible flow

The unsteady aerodynamic loads (pressure increment and generalized forces) acting on a deformable thin airfoil section are determined for its arbitrary motion in the incompressible flow with transverse gusts. In this case, the exact Küssner solution in series is used for harmonic oscillation of an airfoil, in which the Theodorsen function is approximated by a sum of fractional functions with poles. In the time domain, these functions are replaced by the unknown functions that satisfy the one-type ordinary first order differential equations. The equations obtained are combined with the differential equations of the airfoil motion in the generalized coordinates and are used for analyzing the aeroelastic system stability and calculating its dynamic response to gust loads.     

Aerodynamic shape optimization for alleviating dynamic stall characteristics of helicopter rotor airfoil

In order to alleviate the dynamic stall effects in helicopter rotor, the sequential quadratic programming(SQP) method is employed to optimize the characteristics of airfoil under dynamic stall conditions based on the SC1095 airfoil. The geometry of airfoil is parameterized by the class-shape-transformation(CST) method, and the C-topology body-fitted mesh is then automatically generated around the airfoil by solving the Poisson equations. Based on the grid generation technology, the unsteady Reynolds-averaged Navier-Stokes(RANS) equations are chosen as the governing equations for predicting airfoil flow field and the highly-efficient implicit scheme of lower–upper symmetric Gauss–Seidel(LU-SGS) is adopted for temporal discretization. To capture the dynamic stall phenomenon of the rotor more accurately, the Spalart–Allmaras turbulence model is employed to close the RANS equations. The optimized airfoil with a larger leading edge radius and camber is obtained. The leading edge vortex and trailing edge separation of the optimized airfoil under unsteady conditions are obviously weakened, and the dynamic stall characteristics of optimized airfoil at different Mach numbers, reduced frequencies and angles of attack are also obviously improved compared with the baseline SC1095 airfoil. It is demonstrated that the optimized method is effective and the optimized airfoil is suitable as the helicopter rotor airfoil.     

一种计算螺旋桨性能的高效升力线法

本文提供了一种能有效计算螺旋桨性能的方法。方法的特点在于:(1)边界条件中利用了中弧面上的局部斜率以代替平面边界条件,使本方法适用于具有较大弯度翼型的螺旋桨;(2)用近似解析式与数值积分结合起来计算螺旋形尾涡产生的诱导速度,在保证精度的条件下,使计算工作量降低到原有的 10 %;(3)利用翼型的升阻特性考虑翼型阻力及压缩性影响。利用本方法对某型飞机的螺旋桨进行了计算,并与实验数据作了比较,在进距比不太小的情况下,计算结果与实验数据符合很好。      

轴流压气机静子声辐射场的理论分析

本文利用随机脉动调制理论，建立了轴流压气机静子离和噪声和宽频噪声通用的计算模型，对静子声源及其辐射场进行了理论研究，依据现有粘性尾迹经验公式，基于薄翼阵风响应理论，给考虑翼型弯度及攻角影响时静叶叶片和上动叶片尾迹相互作用的脉动升力响应计算公式，分析了压气机几何参数和气动参数对声辐射的影响，并在声功率级频谱图上展示其变化趋势。本文对压气机气动噪声机理的探讨结果将为低噪声轴流压气机的设计提供理论依据。     

可连续光滑偏转后缘的变弯度翼型气动特性分析

以变弯度翼型为研究对象,计算了其六种外形的绕流流场,分析了不同的连续光滑变形翼型与传统偏转翼型的气动特性,讨论了变形参数对气动特性的影响,研究了气动特性的产生机理;与此同时,以形状记忆聚合物柔性蒙皮和机械结构实现了可连续光滑偏转后缘的变弯度翼型,并在风洞实验中测试了其气动特性。计算和实验结果表明:可连续光滑偏转后缘的变弯度翼型能改进传统主翼-简单襟翼翼型的气动特性和流场分离特性;可变形段范围、转轴位置、后缘偏转角度、后缘高度等变形参数对气动特性具有显著影响。     

重型燃气轮机压气机高雷诺数前转捩叶型设计

针对重型燃气轮机压气机雷诺数高而导致的转捩位置前移,开发了一种比可控扩散叶型(CDA)损失更小、工作范围更宽的前转捩叶型.采用正问题优化设计方法,将叶型几何参数化、叶片到叶片流场分析与遗传算法相结合,实现了叶型的自动优化.优化目标综合权衡了叶型损失和攻角范围,为减少优化变量的数目,应用了一种特别的叶型几何模型,将厚度分布与中弧线之间进行了一定的关联.优化得到的前转捩叶型的主要特征是吸力面速度峰值的位置前移至距前缘约10%弦长处,叶型中后部的速度变化更为平缓.最后根据优化结果总结了前转捩叶型的设计规律.      

混合遗传算法及其在叶片自动优化设计中的应用

 在遗传算法中 ,采用单纯形法寻优取代变异运算构成混合遗传算法 ,以提高遗传算法局部搜索能力。算例表明混合遗传算法可有效提高搜索效率和对最优解的逼近程度。采用三次多项式和多圆弧方法生成叶型中弧线 ,三次多项式分布叶型厚度 ,对叶型进行参数化。将N S方程正问题流场数值计算与混合遗传算法相结合 ,构成叶轮机械叶型自动优化设计。尝试由气流转角、总压损失和叶型型面面积构成目标函数 ,对压气机叶型进行自动优化设计。