复杂腔体四极子矩阵三维CFD计算方法.docVIP

复杂腔体四极子矩阵三维CFD计算方法 　　摘要： 应用基于三维计算流体动力学（CFD）的时域方法计算了某复杂腔体四极子矩阵的两种形式：散射矩阵和传递矩阵；根据这两种矩阵和管道末端声学边界分别计算了该腔体的消声量。与实验测量结果对比，验证CFD法的有效性，并分析了两种四极子矩阵的特点，发现传递矩阵计算值更加准确，而且其计算方法更加简单快捷。进一步计算了存在与声波传播方向相反的平均流时该腔体的传递矩阵，得到有流条件下的消声量，分析了平均流对腔体消声量的影响。结果显示有流条件下，腔体消声量有所增加，且其频谱曲线向低频方向移动。 　　关键词： 计算流体动力学； 消声量； 传递矩阵； 空滤器； 流场 　　中图分类号：O351.2； TB535文献标识码： A文章编号：2013　　引言 　　发动机的进、排气系统除了有提供清洁空气和排出废气的作用，还能分别降低进、排气噪声。进、排气系统由管道和消声元件组成，整个系统往往非常复杂。通常将系统分解成一些简单的单元，并用一个二维矩阵表示单元内的声波传播特性，该矩阵即为四极子矩阵（4pole matrix）。将各个单元的四极子矩阵综合，即可得到整个系统的声学特性[1～3]。无论消声元件内部介质的稳态流场和瞬态波动如何复杂，只要进、出口处声波可用一维平面波表示，该消声元件就可以用四极子矩阵表示。 　　四极子矩阵一般有两种形式：一种是散射矩阵，反映了消声元件声学上、下游端正向声波和反向声波的关系；另一种是传递矩阵，也叫阻抗矩阵，反映了消声元件上、下游端声压和介质声学质量速度的关系[4]。两种形式的四极子矩阵可以互相转换。 　　对于一些由简单的膨胀腔、旁支消声器和管道组成的消声器，可用一维线性平面波法和一维有限体积法分别从频域和时域的角度计算其声学特性。其中一维有限体积法将三维实体简化成一维管道模型，通过求解非稳态时域的一维流场，求得频域上的消声性能。该方法在低频区域比较准确，在高频区域不理想[5]。三维有限元法和三维边界元法相比，一维方法能在更宽的频率范围内有效地计算复杂腔体的声学性能[6，7]，但很难考虑介质的粘性及稳态流场的影响。 　　三维计算流体动力学法（Computational Fluid Dynamics，简写为CFD）相比一维有限体积法，可以直接计算三维复杂消声元件的质量、动量和能量方程，得到进、出口处压力波动的时域结果，进而得到频域上的消声元件的声学特性。Middleberg等人对几种几何形状简单的消声器进行了二维CFD建模[8]，得到了不同网格大小对应下的传递损失结果，与测量结果对比发现，在2 kHz以下，8 mm的单元网格即能达到精度要求，但文中并未对复杂腔体进行分析。文献[9，10]用三维CFD方法计算了消声器在无流和有流条件下的传递损失，分析了稳态流场对消声器声学特性的影响。A Broatch等运用FLUENT中的远场压力边界作为声学无反射边界取代出口的长管模型，计算了几种消声器无流和有流条件下的传递损失，与实验结果吻合良好[11]。无反射边界的使用极大地简化了模型，缩短了计算时间。 　　本文在文献[8～11]的基础上，通过合理设置三维CFD方法中的边界条件，以某内燃机进气空滤器的腔体（不带滤芯）为例，计算了无平均流时复杂腔体的四极子矩阵的两种形式，即散射矩阵和传递矩阵。结合管道声学理论，进一步计算了复杂腔体的消声量。将消声量计算结果和实验测量结果对比，间接验证了CFD方法的有效性；并发现根据传递矩阵得到的消声量计算值相比根据散射矩阵得到的计算值，与测量值更加吻合。进一步计算了有流条件下空滤器腔体的传递矩阵和消声量，分析了稳态流场对该复杂腔体声学特性的影响。文中运用三维CFD方法计算复杂腔体的四极子矩阵，考虑了稳态流场和介质粘性对声波的影响，为进一步用三维CFD方法进行内燃机进、排气系统中的消声器设计提供了理论依据。 　　第2期刘联?，等： 复杂腔体四极子矩阵的三维CFD计算方法振 动 工 程 学 报第26卷1CFD方法计算原理〖2〗11散射矩阵在消声器的入口和出口管道内，声波可视为平面波。入射声波传到消声元件后，产生了反射波和透射波，如图1所示。 　　用CFD方法计算消声元件散射矩阵的原理为：分别在消声器的上、下游施加压力脉冲，记录两种情况下的消声器上、下游管道内的压力随时间的变化历程，此过程须保证声波出口端无反射或者不受反射信号的干扰。为了分离声波入口端的入射波和反射波，使用和消声元件上游或下游截面一样的直管，在直管入口施加与计算消声器时一样的压力脉冲，且保证入口与压力测点的距离和计算消声器时一样，记录测点的入射波信号，结合之前消声器的声波入口端的信号，即可分离出入射波和反射波[11]。 　　12传递矩阵 　　对空滤器腔体建模，网格大