高速受限脉动射流发声数值模拟研究.pdfVIP

第十二届全国空气弹性学术交流会会议论文集 CARS-2011-2202 高速受限脉动射流发声的数值模拟研究 周鹤峰赵云曾新吾 (国防科学技术大学湖南长沙410073) 摘要作为调制气流声源研究中的核心问题，高速受限脉动射流的发声问题没有得到足够的认识。一般认 为调制气流声源和人发声器官的换能机理类似，通过FLUENT软件建立了人声道简化结构的仿真模型，利用动 网格技术实现了人发声过程的动态模拟，对不同入口总压条件下声道内非定常流场进行分析，揭示了脉动声源 的形成过程。模拟结果符合气动声学相关理论，可以为气流声源的改进设计提供有益的指导。 关键字调制气流声源、高速受限脉动射流、发声、数值模拟、动网格 受限脉动射流的发声问题在生理学和工程中都占有重要地位。在人发声过程中，声带振动产生 了通过声门(声带之间的狭长缝隙)的脉动射流【I】。声带运动如图l所示。开始时声门几乎闭合， 肺部产生的高压气流使声门分开并形成了起始射流。根据伯努利方程，声带周围的压力降低。接着 弹力和低压使得声门再次闭合并重复这一过程。一个发声周期内声门呈现出三种典型形状：收缩、 等截面和扩张。脉动射流产生的声波构成了发声的主要声源。调制气流声源(AMS)作为一种通过 调制高速气流发声的大功率高压流体声源[21，在强声波传播和效应研究等工作中处于基础性地位。 一般认为，A_MS和人发声器官的换能机理基本相同，两者的发声过程都可以归结为调制气流问题pJ。 八 八． 八、 八． (a)闭合 (b)收缩 (c)等截面 (d)扩张 图l声门运动示意图 值得注意的是，相对于自由射流和受限定常射流的研究工作，针对受限脉动射流的研究较为缺 乏，并且以人发声机理这类低速流动为主，流体可压缩性的影响可以忽略M1。由于AMS中气室压 力远高于上述流动，喷口流速增加至声速甚至更高，即使不考虑其与人发声问题在管道几何形状和 流固耦合作用等方面的差异，流体的可压缩性也会对受限脉动射流流场的演化和声能的转化产生重 要影响嘲。作为AMS基础研究中的核心问题，高速受限脉动射流的发声问题属于气动声学研究领域 的前沿，以往的数值模拟或实验测量并没有对其进行专门的研究，涉及的物理过程较为复杂，因此 我们需要在相关问题的研究基础上建立理论分析模型，系统研究不同参数条件下的流场、声场特性。 本文基于FLUENT软件，采用动网格技术，实现了人发声过程的动态数值模拟，对不同入口总 压条件下通过声门的高速受限脉动射流进行了研究，为下一步探索新的声源性能提升途径和预测声 源辐射声场奠定基础。 1控制方程 气室压力较高时，声道内部为跨声速流动，基本控制方程为时均化Reynolds方程，湍流模型采 用Spalart．Allmaras模型，控制方程如下： 连续方程 587 周鹅峰．赵i等：高速受限脉动射流技卢的数值楔撑日究 詈m(川=o ㈣ 动量方程 巫笋胁(以”)=div(ggradu．)一妾+s (2) 能量方程 掣州肼叫寺删r卜 ㈣ 其中，P，t、H、m f、c。分别表示密度、时间、速度矢量、动力粘度、温度、比热窖·p是流 体微元体上的压力．I=1，2，3分别表示侑卡尔坐标系的三个坐标分量。S是动量方程的广义源项， S．为流体的内热源及由f牯性作用流体机械能转换为热能的部分。 工作舟质为空气，多方指数为I4，密度按理想气体执态方程求解，粘度采用Sutherland三系数 方法计算，气体温度300K。采用有限体积法对控制方程进行离散．可保证物理量的积分守恒在任意 fl