方腔中的涡斑运动发声.pdfVIP

方腔中的涡斑运动发声 刘晓宙刘应中李谊乐缪国平 (上海交通大学海洋工程国家重点实验室，上海，200030) 摘要 本文应用涡声原理研究了育腔中的涡斑运动发声，对其声压随阿间距离的女化，频谱垃 指向性逆行了探讨。结果表明：涡斑的声篷剁与涡斑的运动有关．在不考虑j=晏条件的影阿时，涡斑的 发声呈四极子的特点．不同的涡斑结构的舂辐利有再自趵特点，涡讲的声辐牙j{j现丁涡疆之间的非线性 相互作用。 关键词涡斑，涡运动，声辐射 动以来_遇左垄适得到很大的发展。但能清楚说明声与涡运动关系机理的论文并不多见。 et ’Fang al【3】研究了二维涡系的声辐射，但考惠的是无粘性的情况。我们已对圆柱外的 点涡系的低Mach数运动发声作了系统的研究hj。但二维点涡系统模型距粘性流体中真实 的有核结构的涡结构相距甚远，因此考虑具体确定满嚣分布的集中i呙(涡斑)系统是非 常必要的，它更接近真实的流动情况。本文研究方腔中的二维涡斑运动发声，考虑粘性 的影响，对其发声规律进行探讨。 2运动方程及计算方法 本文首先考虑不可压缩粘性流体四周为嗣壁边界的二维方腔的满斑运动，(如图1) 其中r1．r：分别为两涡斑的环量。 l y b _(- 图1 ·53· 不可j{：缩粘性流体的NaYier-Stokes方程及连续性方程为 o“ 。=一-LP。+¨j， 。 a， -·f ，P {j|月{ “J√=0 其中沩运动粘性系数。 引进满量函数破流函数y — d“2 0Ⅳl ∞=‘o一一2 苏。 舐， V2∥=一厶) (4) 蜥、屹分别为坐标方向的速度分量，方腔的四周为滑移条件． 初始条件为：在静止的方腔 内放置两个等量，旋转方向相同的涡斑。 由初始涡量分布解方程(4)可得流函数妒 再由公式《，5)： d∥ —二-=“． 却 d∥ —o=一“， 以 得到初始速度分“。、u2 由初始速度和无滑移边界条件解Navier—Stokes方程和连续性方程，可得各个时刻U，n 由(3)可计算各时刻方腔中各点的涡量。 为减小方腔腔壁对涡斑运动发声的影响，方腔的长宽取为1．涡斑的半径为0．0I，并且 将方腔四周边缘的涡量不计入发声之内。 3声场计算 按涡声理论．整个流场为一非齐次的波动方程： 孑t矿a2p—V2p=岛V·(。xⅣ) (6) 其中：。为声速，p为声压，肭流体平均密度，劝涡量分布，Ⅳ为流体质点速度．研究 运动发声，将整个流场分为内场和外场，即有旋流场和辐射声场，并利用匹配渐进展开 法进行分析。 ·54· 3．1内部有旋流场 在不可压缩的假定下，(6)式退化为Poi sson方程 V2P，=一风V．(coד) 只为内场中的压力脒动。 引入矢量函数占(丑-I)可得： p，=p。昙p(，!，)·G(x，y)咖 3．2外部无旋声场 由于外部声场无旋且微弱可压缩，(6)式退化为齐次波动方程 2po一甲2风=。 i1可a 其中Pu代表声场中的压力脉动即声压。 对上式进行反变换和多极展开。一!： p。=曩。^¨．a卅，(当)