利用COMSOL软件对圆管声模态的频域模拟.docx

利用 COMSOL 软件对圆管声模态传播的模拟 刘济洲 (北京航空航天大学能源与动力工程学院，102206,北京) 摘要：本文主要介绍利用COMSOL 软件中的频域线化 Euler 方程求解器模拟圆形管道内声 模态传播时的相关理论和具体参数设置方法。在一个简化的无限长圆形管道内，通过在某一 截面加载声模态条件，模拟了不同声模态在无、有平均流动条件下的传播。其中，核心问题 包括圆形管道声模态面内波数的求解以及声模态边界条件的加载。本论文介绍的方法可为模 态声学相关研究提供借鉴和参考。 关键字：管道声模态， COMSOL 软件，数值模拟 0 引 言 管道声模态理论在航空发动机声学、叶轮机声学等领域均得到广泛的应用。该理论由 Tyler和 Sofrin于1962年总结成文[1],提出了管道内线性声波的一种分解方式，即声模态， 并给出了声模态传播规律的分析方法。根据这一理论，可以快速判断某一压气机或涡轮所产 生的噪声的是否可以传播，也可以求解出任意声模态的具体数学表达式。 COMSOL Multiphysics 是一个包含多种物理求解器的商用软件，其中的声学模块可以为 广泛的声学问题提供解决方案，如：机械振动噪声、气动噪声、水声、压电声学等等。目前， 利用COMSOL 软件求解管道声模态的教程[2]采用了声学模块中的频域势流求解器，该方法 可以快速准确地求解圆形或者环形管道截面处的声模态面内波数，并可以将得到的声模态作 为边界条件，进而模拟该模态的传播。但是，该解决方案只能模拟单一声模态的传播，无法 得到两个或者更多模态同时作用下声波的传播结果。在实际问题中，当需要给出多个模态共 同作用的结果时，则需要寻找其它的解决方案。其次，在给定模态强度时，频域线性势流方 法中用于表征强度的物理量为速度势，而在多数声学问题中，更广泛地应用声压、振动速度 等声学概念描述。因此，当该方法与其它方法衔接共同解决声模态相关问题时，在不了解实 际换算公式时，需要进行单位间的校准。 声学模块中的频域线化Euler方程求解器可以求解存在背景平均流动条件下线性声波的 传播，十分适用于管道声模态的传播问题。在该方法中，声模态可通过加载边界条件 (Prescribed acoustic field) 实现，未知数为波动密度、波动速度和波动压力，这与一般声学 问题中的未知数直接对应。其次，可以根据模态的叠加原理，可同时加载多个声模态，从而 模拟多个声模态同时传播问题。相比于频域线性势流的方法，该解决方案在应用时存在以下 两个关键点。第一，求解声模态的面内和面外波数；第二，正确加载声模态作为边界条件。 由于该解决方案中的声模态是通过公式给定的，因此对于确定的几何构型，需要能够求解其 面内、面外波数的求解方法。对于圆形管道，面内波数的求解涉及对贝塞尔函数根的求解， 需要相应的求解程序。而对于圆环形截面管道，其面内波数求解相对更加复杂。在声模态理 论中，被分解的物理量是压力，而声波振动速度和密度通常没有被考虑。但是，在以边界条 件给定声场时，需要在确定压力场的同时给出正确的声波速度和密度，这需要质量和动量方 程的辅助。 在本文中，详细介绍了采用频域线化Euler 方程求解器模拟圆形管道中不同声模态在有、 无平均流动条件下的传播。着重分析了声模态边界条件的正确加载方式，并给出了计算不同 截面形式面内波数的一个外部程序源。 1模型的建立 计算模型简图如图1所示为一个长度为L=1m, 半径为R=0.1m 的圆柱，管道两侧均 布置有PML 层作为无反射条件，给定声源面在z=0m 处定义。通过这样的模型布置，可以 模拟声波在无限长圆形管道内的传播。 计算时选择3维频域线化Euler方程求解器，在求解器设置时可以通过给定z 方向平均 速度值模拟平均流动对声波传播的影响。在z=0m 处，添加 Prescribed acoustic field 边界条 件作为声模态源截面。 图1计算模型简图 2声模态求解与边界条件设置 本小节介绍采用3维频域线化Euler 方程求解器模拟声模态时的两个关键问题。 2 . 1声模态求解 根据管道声模态理论，圆形管道内某一截面的声压可以分解为一系列声模态之和。 假设壁面为硬壁面，则声压在壁面处所需满足的条件为 (1) (2) 根据这一边界条件，可以求解出所有声模态的面内波数 μm。 但是具体求解时，由于涉及 Bessel 函数零点的求解，过程相对复杂。为了方便起见，可采用作者提供的分解程序[3]。 对于截面为圆环形的管道，其在声模态分解时的基函数中包含两类 Bessel 函数，并有 两个壁面边界条件，这将求解面内波数的问题转化为求解 Bessel 函数交叉乘积(Cross product Bessel function)问题的零