剪切流的气动声学数值模拟中边界处理地研究.pdfVIP

2003空气动力学前沿研究论文集 剪切流的气动声学数值模拟中 边界处理的研究 邢超罗纪生 (天津大学力学系，天津300072) 摘要基于特征分析，利用完整的流体力学方程给出计算扰动的新的边界处理方法， 在静止流场中模拟扰动的演化验证了所提边界的处理方法的适用性。采用所提边界 的处理方法，模拟了不同来流速度的均匀流场中扰动及均匀剪切流中扰动的演化，计 算结果表明本文给出的边界条件的处理方法是合适的。 关键词 均匀剪切流 气动声学 扰动 数值模拟 边界条件 1 引言 航空和工业应用中经常会碰到各种形式的由于流动中的脉动引起的噪声问题，典型的例 子有各种形式叶轮机械的气动噪声、燃气噪声、喷气噪声等。气动噪声不仅会污染环境，在军 事上还会影响飞机的防御性能，在某些情况下甚至会造成结构的疲劳和破坏。所以，理解、预 测并最终控制气动噪声的工作越来越重要，也引起学术界和工业应用部门的日益重视…1。 在气动声学的研究中，Lighthill理论被许多学者所采用。此理论将整个声场分为近场和 远场两个部分，近场是声源区域，远场为声音的传播区域。Lighthill理论认为远场的气动声场 后，Lighthill理论可以对远场作出精确的预测，但如何给出近场的声源分布，该理论的结果就 相对比较粗糙HJ。近年来有学者把三层结构理论用于气动声学领域，该理论指出渐进展开的 前三项可以分别作为近场四极子，偶极子和单极子声源，并进一步指出三种声源对于远场的声 波都有贡献[“。 为了研究剪切流中声波的传播机理，本文用数值方法研究了强迫压力脉动产生的扰动的 传播特征。由于计算区域总是有限的，边界条件的处理应使得整个计算在流场中不产生非物 理的数值反射。由于本文的计算是有剪切的流动，因此，没有完全采用Lde等提出的一维无 黏假设的边界条件_6』，而是从完整的流体力学方程出发，结合特征分析给出了新的边界条件 的处理方法，计算结果是比较理想的。 本文以静止流场中加压力强迫脉动产生扰动为例，数值模拟了扰动的演化来验证程序的 可信性。对均匀流动，均匀剪切流动中扰动的演化进行了计算，得到了比较好的结果。本文还 分析了各种流动情况下扰动演化的特征及与流动的关系。 2数值计算方法 研究中，控制方程采用直角坐标系下二维守恒形式的N-S方程，连续性方程和能量方程 剪切流的气动声学数值模拟中边界处理的研究487 掣+差+芒：孥+孥 …7(2川 at az。aYaz。aY 其中，u=(P，∥，p，陋。)’为守恒型通量；adE。 a时F，，十流项及压力项；鲁和考为黏性项 及热传导项。具体形式参考文献L2]。计算中对流项采用特征分裂，对正、负通量在空间上采 用5阶弱迎风紧致格式，对黏性项采用6阶中心型紧致格式进行计算，具体形式如下： 4轭 一 3E2 +r +r) fz訾+s等=忐(一蹈+12Ej+。墙Ej — 2 2 一 一 “E 3E + 一： 2 l lz等+s等：一志(一E--2+12E-一。+36