不可压缩翼型绕流数值研究

使用标准k-ω模型及与色散模型相耦合的k-ω模型分别计算了NACA0012翼型和NACA 4412翼型的低速绕流问题.NACA 0012翼型计算了其来流雷诺数为2.88×106,攻角从0°到15°范围内的流动结构、翼型表面压力分布和升力、阻力特性;NACA 4412翼型计算了临界雷诺数为1.52×106,攻角为13.87°时的流动分离和翼型表面压力系数 ,并与实验数据进行对比.结果表明:在同等条件下,使用与色散模型相耦合的k-ω模型计算得到的NACA 0012翼型的升力和阻力系数比标准k-ω模型提高精度约5%,NACA 4412翼型的表面压力系数精度提高了约3%,进一步验证了其可信性,可将其进一步应用到低速飞行器的气动计算中.      

跨声速翼型优化设计的复变量方法应用研究

应用复变量方法计算目标函数关于设计变量的灵敏度,以此进行优化设计.该方法尽管计算量较大,但求灵敏度简单易行,适合目标函数比较复杂的情形,而且精度很高,几乎不受扰动步长影响,数值优化过程特别健壮,所以非常适用于设计变量不多或气动力计算很快的情况.     

低雷诺数下微型飞行器主动变形机翼非定常气动特性数值分析

以微型飞行器的气动外形设计为研究背景,通过数值求解N-S方程分析了主动变形柔性机翼的气动特性。空间离散格式采用中心有限体积法,时间推进采用双时间推进方法,其中子迭代过程由隐式LU-SSOR方法完成。在动网格技术的基础上,模拟分析了薄翼面上作行波运动的流场。计算结果发现,该行波翼型的气动特性相比于刚性翼型有明显改善,有效地抑制了大攻角下大尺度流动分离,同时升阻比提高了约35%,起到了显著的增升减阻的作用。本文的工作对于进一步的柔性机翼优化具有良好的理论研究与数值模拟的参考价值。     

基于N-S方程的翼型气动优化设计

本文简要回顾了气动优化设计的最新发展，并采用连续伴随方法对粘性条件下的翼型气动外形进行了优化设计。设计过程中，采用B样条曲线的控制顶点作为设计变量，通过求解伴随方程得到目标函数关系设计变量的梯度，最后用最速下降法和BFGS变尺度法对二维翼型进行了外形优化，验证了本方法直接、快速的设计能力。     

非结构混合网格鲁棒自适应技术

数值格式、湍流模型和计算网格是影响CFD数值模拟精度的3个主要因素。结合流场信息的网格自适应技术具备动态优化计算网格的能力,被NASA列为未来CFD发展的一项关键技术。本文针对非结构混合网格,发展了网格单元分布优化、表面网格几何投影和空间网格协调匹配3项关键技术,建立了高鲁棒性几何保真的网格自适应系统。首先,为了提高自适应方法的鲁棒性和通用性,发展了基于标准面网格的多面体网格单元分布优化方法。其次,发展了仅依赖表面网格信息的局部曲面重构技术,采用参数点映射方法实现了新增表面网格点的几何投影,消除了自适应系统对几何CAD系统的依赖。再次,采用改进的距离函数方法实现了空间网格与投影后表面网格的快速匹配。最后,结合基于流场特征的自适应探测器,采用二阶格式的有限体积方法,开展了30P30N三段翼绕流和三角翼大迎角绕流的网格自适应数值模拟。结果表明,通过网格自适应对网格单元的分布进行优化后,流场求解的收敛性和模拟精度都得到了显著提高。     

基于DFFD技术的翼型气动优化设计

开展了直接操作自由变形（DFFD）技术在翼型参数化及翼型气动外形优化设计中的应用研究,应用该方法可以对翼型形状进行直接操纵和精细的局部修型,从而在一定程度上克服了自由变形（FFD）技术无法直接指定几何外形变形量的局限性。通过最小二乘模式根据翼型表面直接操作点的位移求解各个FFD控制点相应的位移,将翼型设计参数从FFD控制点转化为翼型表面的直接操作点,从而有效地减少了高阶FFD控制体进行翼型参数化时的设计参数个数。算例表明,相比于FFD方法,DFFD方法不仅具备直接操纵翼型几何外形的能力,更具物理直观性,并且比FFD方法具有更好的局部变形特性。运用该技术结合遗传算法对RAE2822翼型进行了气动减阻设计,显著减小了设计状态下翼型的阻力,并且可以有效施加如前后梁位置翼型厚度等工程实用的几何约束,证明了该方法的有效性。     

基于遗传算法多段翼型外形优化设计

为提高多段翼型增升效能,开展襟翼外形和缝道参数同时优化设计研究。优化算法采用遗传算法,以求解RANS方程为气动特性分析方法,通过椭圆方程控制生成多段翼型外形,同时优化缝宽、搭接量、襟翼偏角等位置参数和襟翼外形控制参数,实现多段翼型优化设计。设计实践表明,与只优化位置参数相比,同时优化襟翼外形和位置参数得到的多段翼型有更大的升力系数,方法是可行的,具有一定的工程应用前景。     

基于数值模拟的NSDBD等离子体激励器防冰特性

飞行器表面在一定气象条件下会产生积冰,积冰会使飞行器气动性能下降,是危害飞行安全的重要因素之一。常见的气热及电热防冰系统已经广泛运用于现有飞行器上。近些年,在纳秒脉冲阻挡介质放电（NSDBD）等离子体激励器的相关研究中发现NSDBD等离子体激励器可对周围流场进行快速加热,考虑到这种热效应可能作为飞机防冰的一种新方式。本文用数值方法对NSDBD等离子体激励器防冰特性开展了研究。首先,建立了基于Messinger模型的积冰模型,对典型积冰条件进行了验证计算;其次,耦合唯象学等离子体模型与非定常雷诺平均Navier-Stokes方程,计算等离子体对空气流场的影响;最后,将NSDBD等离子体激励器布置在NACA0012翼型前缘防冰区,结合积冰模型与唯象学等离子模型,对其防冰特性进行了研究。计算结果表明等离子体加热的热气流会覆盖在翼型表面防冰区。在相同的霜冰条件下,开启等离子体激励器时机翼前缘没有出现积冰,说明等离体子激励器应用于机翼防冰是有效的。针对不同的激励器参数对防冰特性的影响规律进行了研究,总体上防冰效果与峰值电压、激励器频率有关,从防冰效果和能耗方面考量,在给定计算条件下,存在最优电压值和最优激励器频率值。激励器分布方式对防冰特性的影响与其具体流场有关,需要具体分析。     

合成双射流控制NACA0015翼型大攻角流动分离试验研究

设计了一种卧式合成双射流激励器(DSJA),并对其在翼展中段控制NACA0015翼型大攻角流动完全分离进行试验研究,分析了合成双射流激励器两射流出口位置及射流能量对控制机翼流动分离的影响规律。结果表明:合成双射流激励器对机翼大攻角流动分离具有很强的控制能力,可显著提高机翼流动分离攻角;合成双射流激励器两射流出口相对分离点的位置是影响控制效果的重要参数;合成双射流激励器两出口任一出口位于分离点之前,且越靠近分离点,其对边界层分离的控制效果越好,并且当分离点位于合成双射流激励器两出口之间,且离第一出口位置较近时,合成双射流"接力"控制机翼分离的效果更加明显;与合成射流"单射流"相比,合成双射流"两射流"对分离点位置的有效控制区域明显增大。此外,提高合成双射流激励器的射流能量,其控制机翼流动分离的能力提高。     

翼型动态特性数值模拟及其影响因素分析

采用非结构动网格技术模拟了NACA 0012翼型绕1/4弦长位置俯仰的强迫振荡过程,验证了翼型的周期性动态气动特性。分析了不同的湍流模型、单周期内步数和内迭代次数等因素对计算结果的影响;研究了马赫数和减缩频率对动态气动迟滞特性的影响。结果表明, Transtion SST模型得到的结果更准确;单周期内计算步数大于200时,计算结果收敛性更好;随着马赫数变大,迟滞效应先增大后减小;在亚音速区,迟滞环沿逆时针方向,动导数小于零,模型动稳定;而在跨音速区,迟滞环沿顺时针方向,动导数大于零,模型动不稳定;减缩频率越大,翼型动态迟滞效应越明显。