基于Lilley方程射流噪声研究.pdfVIP

中固]==程热物理学会 热机气动热力学 2007年学术会议论文 编号：072024 i II 基于L ey方程射流噪声的研究术 甘加业吴克启 华中科技大学能源与动力工程学院，武汉430074 Tel：02783 63．com Email：g酬iaye0775@1 摘要：本文基于声源模型理论，研究了射流气动噪声频谱。此模型理论基于L订ley方程来构造伴随格 林函数(adjointgreenfunction)，且模型关联，声源和远场声压。声场的边界通过施加完全匹配层 避免声波产生数值反射。采用高阶预测一校正的计算方法，数值汁算控制方程。计算结果与前人的实 验相一致，有效预测了射流中小尺度湍流噪声的特性。 关键词：Lilley方程；格林函数；完全匹配层边界条件； 湍流 引言 众所周知，来自高速射流的湍流混合噪声由两部分组成，分别是大尺度湍流结构形 成的噪声，以马赫波辐射；及指向性弱的小尺度湍流噪声。马赫波辐射控制着一卜-游的声 压频谱，而小尺度湍流噪声则控制着射流两侧和上游方向。预测射流噪声的主要困难在 丁，人们对湍流形成机理仍然缺乏足够的认识，当前射流噪声理论主要还是基于半经验 模型【l】o 学类比法得到不同改进怛一j。其中，“lleyl4J将声的输运、散射项从源项中抽出而移至波 动算子中，获得了更趋实际的控制方程。此方程为三阶波动方程，通过减小源项的舍入 误差，提高了噪声预测精度，适用于预测高频声波的散射和反射现像，研究声对流动的 反作用过程。Goldstein【’J回顾了Lilley方程的理论计算方法。 求解Lilley方程通常是引入格林函数，可以考虑到平均流动中声的反射效应，而此 效应的计算独立于声源的分布，同时把经验系数的误差局限在声学类比法中的源项求解 上【61。Tam【7】发展了一种关联声源和远场声压的格林函数，在此称为伴随格林函数。此 方法的优点之一是：在射流场中所有声源点位置的格林函数可以只通过一次计算得到， 火火。弦省了计算时间，目前此方法在射流噪声研究领域备受关注pj。 本文基于前人的实验数据，主要研究射流中小尺度湍流噪卢。采用伴随格林函数法， 同时引入完全匹配层边界条件，构建声传播方程。控制方程的求解采用高精度的预测一 校正计算方法，即采用时间6阶、空间10阶。 2． 控制方程 ey方程 2一L．1 考虑如’F图1所示的平行剪切均匀流动。 p=po=c9即s，■=谚1u(x2，x3)，丁=瓦(x2，屯) +本文为国家自然科学基金资助项日(No、高校博J‘学科点专项科研綦金资助项目 (No．20060487036)及华中科技大学创新基会资助项H的前期7{：J|究成果 744 对丁一此流动，Lilley的三阶波动方程将线性化【5】 皿兰旦旦一2型堕型 Dt瓠l瓠j瓠i瓠每j 式中，x三(x1，屯，屯)为笛卡尔坐标，五 为流动方向，p为压力，p=p’+夕o， 如果脉动压力夕7相对于平均风是小量， 则自变量n三p’／(y风)，c2=y尺瓦 f 为平均声速，D}Dt三8}8t+U8ax、 为对流导数算子。’ 2．2伴随格林函数 3．数值方法 3．1计算网格 计算网格的生成采用高斯一梯形混合积分法则【0】，这种方法有利于提高计算效率及 精度。网格梯度参数为7，相应控制点参数计算结果为： a(1)=O．099193378414510，《2)=0．507659266964552，a(3)=1．1 1869310， a(4)=2．047493467134072，《5)=3 a(7)=5．000007879022339，a(8)=6