基于控制理论和N—S方程的二维叶栅气动优化算法.pdfVIP

第 3O卷第 12期 工 程 热 物 理 学 报 Vo1．3O．No．12 2009年 12月 JOURNAL OF ENGINEERING THERM OPHYSICS Dec．，2009 基于控制理论和 N—S方程的二维叶栅气动优化算法 厉海涛 宋立明 李颖晨 丰镇平 (西安交通大学叶轮机械研究所，陕西 西安 710049) 摘 要 基于控制理论的气动优化方法的计算量与设计变量无关，可快速精确地完成控制变量的灵敏度分析 。本文 以出 口熵增最小为 目标函数，详细推导 了二维 N—s方程伴随系统的偏微分方程组及其相应的边界条件和敏感性导数的表达式。 以拟牛顿算法为优化求解器，利用 CFD方法求解流动变量，采用时间推进方法求解伴随方程，建立了基于黏性伴随方法 和 N—s方程的二维叶栅气动优化设计算法 。 关键词 气动优化；控制理论；黏性伴随方程；熵增 中图分类号：TK124 文献标识码：A 文章编号 2D AER0DYNAM IC oPTIM IZATIoN ALG0RITHM FoR TURBINE BLADESBASED 0N C0NTRoL THE0RY AND N—SEQUAT10N LIHal—-Tao SONG Li--M ing LIYing·C·hen FENG Zhen--Ping (Instituteofff-hrbomachinery,Xi’ai1JiaotongUniversity,Xi’ 710049，China) Abstract A computationalcomplexity ofdesign optimization ofaerodynam ic shape based on controltheoryisindependentofthenumberofdesignvariables，anditcanrealizethequickandexact sensitivityanalysisofcontrolvariables．Beingconsidered asthecostfunction，them inimizationof outletentropygenerationisappliedtoadirectdesign，andtheadjointpartialdifferentialequation (P．D．E)ofthe2D N—Sequation，aswellasitsboundaryconditions，andthefinialexpressionof sensitivitygradientarededucedindetail．Combiningquasi—Newtonmethodasanoptim izationsolver， usingCFDmethodtocalculateflowvariables，introducingtime—marchingmethodtosolveanadjoint equation，wedevelopdesignoptimizationofallaerodynamicshapealgorithmbasedonviscousadjoint methodandN—Sequationfor2D turbineblades． Keywords aerodynamicoptimization；controltheory；viscousadjointequation；entropygeneration 0 引 言 目前常用的气动优化设计方法大致可分为全局 伴随方法最先是 由Pironneau[5J引入到流体动 优化方法和梯度优化方法。全局优化方法包括进化 力学形状反问题求解 中，Jameson[。J最早将其应用 算法 (EAs)[、实验设计方法 (DOE)_2J以及响应面 在航空领域 ，其后得到了发展。近年来国内也展开 方法 (RSM)[。】等，这类气动优化方法有着优秀的鲁 了基于控制理论的 2D、3D翼型气动反设