[论文解读] Simple analytic equation for airfoil shape description

David Ziemkiewicz|arXiv (Cornell University)|Dec 30, 2016

Advanced Numerical Analysis Techniques参考文献 9被引用 8

一句话总结

本文提出一种简单的解析方程，仅使用六个直观参数——基准形状（B）、厚度（T）、尖削指数（P）、中弧度（C）、中弧度指数（E）和反射度（R）——对二维机翼剖面进行参数化描述。该方法可实现高效的手动设计与基于优化的设计，参数最少，生成的气动外形平滑连续，非常适合基于XFOIL的气动分析，经进化优化后，升阻比（L/D）相较Clark Y参考翼型提高了20%。

ABSTRACT

We show a simple, analytic equation describing a class of two-dimensional shapes well suited for representation of aircraft airfoil profiles. Our goal was to create a description characterized by a small number of parameters with easily understandable meaning, providing a tool to alter the shape with optimization procedures as well as manual tweaks by the designer. The generated shapes are well suited for numerical analysis with 2D flow solving software such as XFOIL.

研究动机与目标

• 开发一种参数最少、具有解析形式的方程，用于描述具有直观几何意义的二维翼型轮廓。
• 实现气动优化中高效的手动与自动形状修改。
• 确保与XFOIL等二维流场求解器兼容，以实现性能的精确预测。
• 在保持平滑性与连续性的同时，实现形状生成的高灵活性。
• 支持逆向设计与优化工作流，并实现参数影响的可预测性。

提出的方法

• 使用θ ∈ [0, 2π] 范围内的参数方程X(θ)与Y(θ)定义翼型几何形状，通过绝对值与幂函数控制形状。
• 利用基准形状系数B控制前缘曲率：当B=2时为椭圆形状，B→1时趋近于矩形。
• 将厚度分布建模为Th(θ) = T × |sinθ|^B / sinθ × (1 - X^P)，其中P控制后缘尖削程度。
• 将中弧线定义为Cm(θ) = C × sin(X^E × π) + R × sin(X × 2π)，其中E控制中弧线位置，R可实现后缘反射或襟翼状结构。
• 通过曲线拟合校准参数，使生成的翼型与参考翼型（如Clark Y）的误差最小化。
• 使用进化算法优化，以在多个迎角范围内最大化升阻比（L/D）。

实验结果

研究问题

• RQ1 是否可通过最少参数的解析方程，以直观的几何控制方式描述广泛的翼型形状？
• RQ2 各个参数（B、T、P、C、E、R）如何影响翼型形状与气动性能？
• RQ3 该参数化模型在可接受的精度范围内，能在多大程度上复现已知翼型特性（如Clark Y）？
• RQ4 使用该模型进行进化优化，是否能显著提升升阻比，相较基线翼型？
• RQ5 所生成翼型的平滑性与连续性在多大程度上影响其与XFOIL及流场分析的兼容性？

主要发现

• 参数化翼型模型成功复现了Clark Y翼型，几何偏差极小，可实现可靠的基线匹配。
• 优化后生成的翼型在Re=10^6条件下，最大升阻比相较Clark Y参考翼型提高了20%。
• 优化后的翼型具有比Clark Y更低的俯仰力矩系数（Cm），提升了纵向稳定性。
• 在高升力系数（Cl ≈ 0.9）时，阻力显著降低，有助于提升升阻比性能。
• 进化算法（10代，每代9个个体）成功在广泛迎角范围内提升了升阻比。
• 该方法支持可预测、直观的形状修改：增大P可增强层流特性，而R可实现后缘反射以提升稳定性。