一种快速计算螺旋桨气动声学特性的数值方法-空气动力学学报.pdf

第24卷第2期 空气动力学学报 V01．24．No．2 2006年06月 ACTAAERODYNAMICAs】【NICA Jun．．碱 文章编号：0258．1825(2006)02．0187．07 一种快速计算螺旋桨气动声学特性的数值方法 高永卫，乔志德 (西北工业大学翼型研究中心，陕西西安710072) 摘要：在低速螺旋桨性能计算方面，涡格法是快速有效的方法。但是在计算表面压力分布时，传统的涡格法遇到 了无法克服的困难。本文推导出在叶片表面布置分布涡强的方法，既保持了涡格法计算快的优点，又较好地计算 出叶片表丽压力分布。以此为基础，本文采用阿一H方程忽略四极子效应的时域解法，使用与气动计算相同的表 面网格，快速解出了螺旋桨的远场声学特性，计算结果与实验结果符合程度较好。说明该一体化方法可以用于螺 旋桨噪声机理分析和低噪声螺旋桨设计，从而为分析螺旋桨噪声机理和设计低噪声螺旋桨提供简便工具。 关键词：螺旋桨；气动声学；数值计算 中圈分类号：V212．4 文献标识码：A O 引 言 1螺旋桨表面压力分布的数值计算方 法 螺旋桨发出的噪声可以分为旋转噪声(离散噪 声)和宽带噪声两部分。旋转噪声是叶片与来流周期 螺旋桨总的气动特性(如推力、扭矩、效率等)可 性相互作用引起的，而宽带噪声是叶片与周围流场的 以使用片条理论和升力线等简单方法计算，但要准确 随机脉动相互作用产生的。大量研究表明，旋转噪声 计算桨叶表面的压力分布就必须使用升力面，全速势 是螺旋桨的主要噪声，所以对螺旋桨噪声的研究工作 方程，Euler方程甚至N．S方程等更精细的方法。本 主要是围绕确定旋转噪声而展开的。旋转噪声又可 论文主要研究速度较低的螺旋桨，桨尖Ma数一般小 分为：厚度噪声、负荷噪声和非线性噪声。厚度噪声 于0．7。文献[1]指出在这样的范围里完全可以使用 是由叶片的有限厚度引起的旋转单极子噪声。负荷 桨叶表面布涡的升力面方法来模拟流动的基本特性。 噪声是叶片升力和阻力引起的旋转偶极子噪声。非 1．1传统的涡格法 线性噪声又称为四极子噪声，包括非线性噪声源和非 线性传播效应两项，也是以叶片通过频率为基频的单 传统的升力面方法通常采用涡格法，具体做法是 音噪声，但四极子噪声只有在叶尖工作在跨音速和超 将桨叶沿展向分成若干段，每一段沿弦线方向在桨叶 音速条件下才显示出重要性。这样，对于亚音速螺旋 上下表面均匀布置四边形涡环，由最后一些面元拖出 桨噪声来讲，厚度噪声和负荷噪声是主要的噪声，其 马蹄涡在尾缘处相遇，伸向无穷远。展向附着涡位于 中负荷噪声是最重要的噪声。 每个面元的前缘，两个面上的第一根附着涡在桨叶前 要计算负荷噪声必须知道叶片上力的分布情况， 缘重合。 即需要知道叶片表面的压力分布情况。用Euler方程 边界条件为桨叶表面法向分速为零，即： 或N．S方程计算螺旋桨表面压力分布需要的计算机 V·露=0 (1) 资源和计算时间都是巨大的。本文试图利用升力面 若桨叶表面的方程为： 方法发展一种既简便又具有一定精度的计算螺旋桨 彳=，(菇，Y) (2) 表面压力分布的方法，从而达到快速模拟螺旋桨声学 则表面法线方向由下式确定： 特性的目的。