06基于NS方程的旋翼翼型反设计分析万俊明6.docVIP

第二十届（20）全国直升机年会论文南京航空航天大学摘要：方程为主控方程流场求解使用Poisson方程为控制方程网格采用MGM方程作为反设计方程，。应用方法，进行了反设计分析，满足要求的，并与目标吻合良好。1 引言 直升机是以旋翼为主要升力件的垂直起落飞行器，旋翼的气动特性的好坏决定了它的垂直起降、空中悬停等性能，而旋翼的气动特性又与旋翼桨叶的气动外形有着十分密切的关系。国外从上世纪七八十年代就专门针对直升机旋翼开展了翼型气动设计方法的研究，并且很多国家都已经有了自己的直升机专用翼型，如法国的OA系列，美国波音公司的VR系列，俄罗斯的TsAGI系列等发挥了重要的作用但目前翼型气动设计的方法主要有两种：优化设计方法和反设计方法。优化设计方法需要对流场求解器进行多次调用，使得该方法的计算量相当大，计算效率；反设计方法多余量修正法，其中流场求解与反设计迭代模块之间相对独立，因此反设计方法可与先进CFD（计算流体力学）流场相结合，在提高计算效率的同时，也能有效地提高翼型设计的。uler方程并采用余量修正反设计方法进行了旋翼翼型的反设计研究工作。但由于Euler方程没有考虑粘性，因此在计算受粘性影响较大的气动特性时精度不够，尤其是对翼型阻力的计算不够准确，影响翼型反设计的准确度。而N-S方程由于加入了粘性项，使得它能够更好的模拟流体的粘性运动，特别是翼型附面层的粘性运动，因此使用N-S方程能更准确地计算翼型的压力分布、升、阻力系数与力矩系数等气动特性参数，从而提高了翼型反设计的准确性和可靠性。 综上所述，本文采用Poisson方程为控制方程网格开展直升机旋翼翼型气动设计方法研究2 旋翼翼流场数值计算方法 2.1流场求解控制方程 本文反设计过程中的流场求解模块使用方程为主控方程 (1) 式中：,,,, ,，，，，其中、、、、、、分别为空气密度、方向的速度分量、、压、总焓,,分别为粘性系数、热传导系数和绝对温度，湍流模型采用Baldwin-Lomax模型。 2.2流场求解方法 对每个网格单元，方程（）可用通量表示为 ) 式中，代表守恒变量，代表无粘通量， 代表粘性通量，代表人工耗散项。 在时间上采用显式迭代格式进行时间离散，为减少计算量，在五步Runge-Kutta迭代过程中，人工耗散项只在第一步进行计算，其他步都取为相同值。为加快收敛速度和迭代过程的稳定性，当时间步长，并在五步Runge-Kutta方法的第一、三、五步进行隐式残值光顺等加速收敛措施。 .3旋翼翼型流场求解算例 为进一步验证流场求解程序的可靠性，本文选取，与验结果进行对比。，其中翼型上点数为189，法向上点数为53。 算例一：选用NACA0012翼型为验证翼型，来流马赫数为0.68，迎角为，雷诺数为。算例二：选用NACA0012翼型为验证翼型，来流马赫数为0.70，迎角为，雷诺数为。 图分别给出计算的压强系数分布与验值的对比，可以流场结果与验值，流场能够满足翼流场和旋翼翼型反设计的要求。 图压强系数对比 图压强系数对比3 网格生成方法 .1翼型网格生成方法 (3) 其中，，为源项，当p、q为0时，与分别为计算平面和物理平面的坐标。 图3 源项控制原理 由于N-S方程对附面层的要求，所生成的网格最内层网格间距要求在量级，因此生成的网格的网格单元数比Euler方程多，尤其是靠近翼型物面处网格相比Euler方程要更密集。图图图图图4 反设计方法 4.1翼型反设计Garabedian和McFadden发展了一种通用反设计方法：GM方法，但是在使用中发现，该方法设计出的新翼型激波附近不能很好与目标压力分布重合Malone[9]等人GM方法，得到MGM（Modified Garabedia -McFadden）方程并用于翼型的反设计。 (4) 式中，、和是与关的一组参数（文中称为反设计系数）为目标翼型的压强系数，为翼型表面的压强系数，为翼型坐标的量。4.2反设计方程的求解 在翼型上下表面分别对式（）进行离散，原反设计方程（）可表示为： (5) 对于下表面有 对于上表面有 4.3 旋翼翼型反设计本文旋翼反设计方法其迭代过程如图所示 图反设计计算过程示意图5 反设计算例分析及讨论 5.1一选用NACA00为初始翼型，NACA为目标翼型。来流马赫数0.7，迎角，。图9、图10和表1分别给出了初始翼型、目标翼型和计算翼型的外形几何形状、压强系数、升力系数(C