跨音速翼型反设计的同伦方法研究.pdfVIP

跨音速翼型反设计的同伦方法 崔凯杨国伟 (中国科学院力学研究所高温气体动力学重点实验室，北京，100080) 摘要：跨音速翼型的反设计问题可以表达为复杂非线性方程组的零点求解问题，为增大求解的收敛范围，本 文依据同伦方法的思想，通过构造不动点同伦，将原问题的求解转化为其同伦函数的求解，并依据拟Sigmoid 函数调整同伦参数，进而构造出一种具有较高计算效率的大范围收敛反设计方法。数值算例以RAE2822翼型 的表面压力分布为拟合目标，分别采用B样条、PARSEC、正交形函数三种方法对翼型进行参数化，并分别 法，且具有较好的计算效率，从多个不同的初始翼型出发，经较少次数迭代后，均能与目标翼型很好的拟合， 是一种高效的大范围收敛方法。 1．前言 近年来，随着高性能计算机、数值计算方法及计算机可视化技术的飞速发展，在计算流 体动力学(CFD)和数值优化方法研究日趋成熟的基础上，航空航天飞行器气动外形设计方法的 研究已越来越引起人们的重视。目前气动外形设计的方法可主要分为两大类，即直接优化设 计方法和反设计方法。 直接优化设计方法是在满足一定约束的条件下，尽可能最大化或最小化某些设计目标(如 最大化升阻比、最小化阻力系数等)，从而将外形设计问题转化为一个约束优化问题求解。依 据所使用的求解方法的不同，直接优化设计方法可以概括分为全局方法和局部方法两类。全 局方法如遗传算法(GA)¨，2|、模拟退火算法(SA)等【3u，其主要优点是可能寻求到全局最优解， 适合于求解存在多个局部极小点的设计问题，而且这类方法也不需计算目标对设计变量的灵 敏度。但本质上而言，这类方法属于随机搜索方法，需要多次进行流场计算，因此计算效率 相对较低。局部方法即梯度型方法【4卅，是采用目标函数的梯度作为搜索方向的一种迭代方法， 虽然这种方法只能保证得到局部最优解，但所需的流场计算次数较少，因而具有较好的计算 效率，同时一般也可以得到一个较好的解，因此对于许多设计问题，特别是对于流场计算是 一个相当耗时的大型设计问题而言，梯度型方法仍是一种广泛使用的方法。 反设计方法[7qO】是以预先设计或要求的飞行器表面压力分布为目标，通过计算不断修正飞 行器的外形，使其表面压力分布的计算结果与目标值达到最佳拟合，从而得到飞行器外形的 一种设计方法。反设计方法的本质是将气动外形设计转化为一个反问题求解，一般也采用全 局或梯度型优化方法求解。当采用全局方法求解时，同样需要较大的计算消耗，而当采用一 般梯度型方法求解时，因反问题本身所具有的不适定性，特别是求解过程的不稳定性使求解 过程强烈地依赖于迭代初值的选取，若选取的初值与解相距较远，很可能造成迭代的发散。 针对反设计中使用梯度型方法出现的问题，为增大求解的收敛范围，本文依据同伦方法【ll】 的思想，通过同伦映射，将原反问题转化为其不动点同伦方程的求解，同时在迭代中依据拟 Sigmoid函数修正同伦参数，进而提出一种具有较高效率的大范围收敛反设计方法，并应用于 跨音速翼型设计。为验证该方法的大范围收敛性及对不同翼型参数化方法的适应性，我们分 别采用B样条方法【12l、PARSEC方法【131及正交形函数方法【14】等三种不同的翼型参数化方法， 翼型进行计算，计算结果显示了该方法的大范围收敛性和高效性。 2．流场计算的数值方法和网格变形技术 在反设计过程中，需要对不断修正的翼型进行流场计算以获得压力分布。本文中，流场 397 第十一届全国激波与激波管学术会议 三维程序用于二维流动模拟，需要在翼型展向拓展二层网格。为保证流动的二维特性，在此 方向的流场边界设置为对称边界条件。 所使用的三维薄层NS程序采用了包含子迭代的LU．SGS隐式时间推进方法和改进的 HLLEW空间离散格式。子迭代可以消除由于方程线性化，滞迟边界条件等带来的误差，使计 算更加稳定。改进的HLLEW格式在等熵流动区域，格式自动恢复到迎风差分Roe通量分裂 格式，在熵跳跃大的流动区域，格式变为标准的HLLEW格式。这样既可克服Roe格式模拟 激波流动产生的非物理振荡现象，也可克服HLLEW格式在连续流动区域表现出的较强的格 式粘性。提高流场整体计算精度。此外，程序中粘性项用中心差分离散，湍流模拟采用涡粘 性B-L代数摸型，计算网格为C型网格。初始网格通过代数方法和求解包含控制网格边界正 交的二维椭圆形方程光滑