一种基于共轭方程法求解黏性反问题的简化方法 An Approximate Method for Viscous Inverse Design Based on Adjoint Equations.pdf

2012Vol33No．4597—606 航空学报 Apt．25 ActaAeronauticaetAstronauticaSinica ISSN1000—6893CN 11．1929／V hkxb@buaa．edu．cn Mtp：#hkxb．buaa．(3du．cn 文章编号：1000—6893(2012)04—0597一10 一种基于共轭方程法求解黏性反问题的简化方法 杜磊，宁方飞’ 北京航空航天大学能源与动力工程学院，北京 100191 擒要：对于给定压力分布的黏性气动反问题，考虑到壁面微小扰动造成的压力变化主要由势流作用引起，因此可以简 化用以获得目标函数对设计变量敏感性导数的共轭方程。将黏性流场插值到粗网格中作为彻体力模型方程的解，则其 相应的共轭方程将以简单的源项取代原方程中复杂的黏性项。由于在粗网格中求解。网格数减少，同时收敛速度加快， 简化的共轭方程计算时间可以减少到黏性方程的十分之一。典型的算例结果表明，对于附着的边界层流动简化方法计 算得到的敏感性导数具有较高的精度，能够有效完成反问题设计且减少总的计算耗时。 美键词：黏性流动；彻体力模型；共轭方程；反问题设计；敏感性导数 ^ 中图分类号：V211．3 文献标识码：A 反问题设计作为一种重要的气动设计方法， 者根据连续的流动方程得到，然后离散求解，后者 可以解决给定壁面物理量分布求解其相应几何形 可直接由离散的流动方程推出。20世纪90年代 状的问题。给定的物理量可以是马赫数[1]、速 Jameson[8]最先提出了基于连续型共轭方程的气 度‘21、压力‘3。53或者负荷分布∽3等不同形式，因而动设计方法，随后将该方法应用到翼型的无黏和 存在多种不同的反问题求解方法，其中以偏微分 有黏反问题设计以及三维机体的无黏反问题设计 方程系统控制理论为基础的共轭方程法(AEM)中[9。“。同时离散型方法也得到了相应的发展， 获得了广泛的应用。该方法将给定的目标压力和 计算压力之间的差值作为目标函数，通过求解一 研究工作。Nadarajah等[1们对这两种技术途径进 次流场和一次与流场计算相当的共轭方程即可获 行了分析比较，结果表明只要加以很小限制两者 得所有设计变量的敏感性导数，再由梯度类优化 均可以获得精度相当的敏感性导数。早期该方法 算法求得目标函数最小值，实现反问题的求解。 主要被用来解决外流的气动设计问题，后来被拓 与传统的单独扰动每个设计变量然后再求解 展到压气机和涡轮等内流流动，包括带有掺混面 相应流场来获得敏感性导数的有限差分法 的多级情况[1¨引。国内也开展了共轭方程法用于 气动设计的相关研究[1纠”。 (FDM)相比，共轭方程法的最大优势在于获得 梯度所需的流场计算量与设计变量个数无关， Jameson推导黏性流动的连续型共轭方程时 特别适合解决多设计变量复杂几何外型的气 假定几何微小变化并不会引起黏性系数(分子黏 动反问题设计。 性和湍流黏性之和)的变化，并将其视为常数[2“， 共轭方程存在连续型和离散型两种形式，前 只考虑平均流方程，这样可以简化共轭方程的推 收稿日期：2011-06-21；退修日期：2011-09-23；录用日期：2011-10·14；网络出版时闻：2011．10．3110：58 cnkj 网络出版地址：vzr”／w net／kcms／detaiUll