基于响应面法低速翼型气动优化设计.pdf

第28卷第4期 空气动力学学报 Vol-28。No．4 ACTA SINICA 2010年8月 AERODYNAMICA Aug．，2010 文章编号：0258—1825(2010)04—0430-06 基于响应面法的低速翼型气动优化设计 邓 磊1，乔志德1，杨旭东1，熊俊涛2 (1．西北工业大学翼型、叶栅空气动力学国家级重点实验室，陕西西安710072} 2．机械与航空与太空工程系，美国加州大学) 摘要：响应面方法相较于其它直接优化方法有其高效、实用的优势，此前的研究更多地将响应面方法用于超音速 和跨音速翼型的减阻优化设计中。本文将此方法应用于低速翼型优化设计中，进行了基于RANS(Reynolds—Aver— agedNavier-Stokes)方程和自由转捩预测耦合求解的低速翼型气动优化设计。通过计算附面层方程得到附面层参 数并用e”方法计算转捩位置，并考虑了T_S波和层流分离造成的转捩。RANS方程计算中，使用了转捩过渡区模 型，以保证附面层外边界压力分布的精度。RANS方程和转捩预测迭代进行至转捩位置收敛。在响应面模型计算 中，使用不含二阶交叉项的二阶多项式模型．减少了构造模型所需的计算量；合理的选择设计空间保证了构造的响 应面模型具有较高的精度。使用三个设计点的多目标优化设计，保证了设计的合理性。通过对NACA64(1)一112 翼型优化计算结果表明，本文的方法可以有效地进行低速翼型的气动优化，各设计点上转捩位置也得到了改善，有 较好的工程实用前景。 关键词：响应面方法；RANS方程；转捩预测；低速翼型；优化设计 中图分类号：V211．41+2文献标识码：A 分布作为附面层方程求解的输入，求解附面层方程得 O 引 言 到附面层参数和速度型。当有层流分离时，使用交互 基于CFD的气动优化设计已成为当前计算流体 式附面层求解方法，避免了附面层的解在分离点的奇 力学研究中的一个重要领域。近年来，响应面方法由 异性，又满足了分离点之后回流的存在[1]。得到附面 于其高效、实用的特点，被研究者应用于各种不同领 层参数和速度型之后，使用e“方法求解出翼型上的 域。响应面方法可以通过较少的试验获得设计变量 转捩位置，带入RANS方程的求解中。如此耦合迭 与性能间足够精确的相互关系，并用简单的代数关系 代求解，直至转捩位置收敛。在RANS方程求解中， 式展现；同时，在进行气动优化设计中能消除高频噪 由于点转捩模型引入了局部扰动，使压力分布在转捩 声，因此可以期望得到全局的近似最优解，并在优化 位置不连续，这时需要使用转捩过渡区模型，消除局 设计过程中针对不同目标函数和约束条件，不需要增 部扰动对压力分布的影响，保证下一次迭代中转捩位 加额外的计算量，给设计带来很大的方便。 置计算的精度。 在低速翼型设计中，翼型阻力相当大的部分是粘 在响应面模型构造过程中，如果使用完全二阶多 性阻力，而粘性阻力与转捩位置有重要关系。在 项式，则构造模型所需进行的流场计算次数和设计变 RANS方程求解程序中，一般使用全层流或者全湍量个数的平方成正比，限制了响应面方法进行多变量 流或者固定转捩计算，不能满足低速翼型优化设计中 和多目标优化。如果去掉二阶交叉项，则流场计算次 对阻力计算的精度要求，因此必须使用一定精度的自 数和设计变量个数成线性关系，使进行多设计变量和 由转捩预测。基于小扰动理论的线性稳定性分析的 多目标优化设计所需的计算量可以接受。 e“方法由于其计算花费少、适应范围广等优点，广泛