基于有限元/边界元耦合方法的管道进口声传播及声辐射模型

发展了一种基于有限元和边界元耦合方法的管道进口声传播及声辐射计算模型.该模型将整个声场分为内部有限域和外部无界域,分别用有限元和边界元方法求解控制方程,在两者之间的界面上使用具有物理意义的声阻抗参数进行匹配,并通过一种快速迭代方法实现全声场求解.这种迭代方法可以保证有限元刚度矩阵等带宽以及对称的特性不被破坏,有助于提高计算效率.该模型先得到了Levine-Schwinger标准解的检验,进而在无流动情况下对于简化的航发短舱进口管道模型进行了噪声辐射现象的数值模拟,最后基于计算结果分析了声衬对远场声辐射的影响.      

快速声散射方法在变截面管道中的应用

 快速声散射方法(FSM)是一种基于无流动Helmholtz方程边值问题的声散射预测工具,具有快速、灵活的特点.以航空发动机消声短舱的声学设计为背景,用该方法对变截面管道声传播特性进行数值模拟研究,用直接边界元方法(DBEM)进行数值求解,避免了求解管口反射系数,有效地提高了计算速度.数值研究了刚性壁面和不同位置壁面声衬组合对变截面圆环管道形状声传播的影响,并分析了管道厚度对散射声场的影响,数值结果与声类比方法结果进行了相互验证.最后,还对一种真实转子声源进行了管道声散射的数值研究,结果表明本方法在航空发动机声学设计中具有工程应用价值.     

一种可用于声衬优化设计的管道声传播有限元模型

以航空发动机消声短舱声学设计为目标, 发展了一种管道声传播有限元模型, 它包含了壁面声阻抗、管道截面积和形状变化、以及管内非均匀流动等典型影响因素.应用该模型针对一种真实的高速转子声源进行了管道消声的数值模拟, 不仅得到了细致的管内声场分布情况, 更利用其快速的特点完成了声衬声阻抗的优化计算.结果表明该模型具有工程应用价值.      

一种新的有限长等截面航空发动机短舱声学模型及数值结果

基于航空发动机短舱声学设计的需要,发展一种与以往工作不同的有限长管道声学模型声场的计算方法——直接边界元和传递单元法相结合。本方法对管道内部几何条件和声衬的敷设情况有更少的限制,可以很好地用于分段声衬的优化设计,并且提供求解有限长变截面管道声学模型的基础。     

变截面风扇管道快速声散射计算

针对航空发动机消声短舱声学设计的需要,发展了一种结合快速声散射和直接边界元计算发动机远场声散射的方法,用于计算包含了壁面声阻抗影响的复杂截面管道声散射问题.应用该方法针对某风扇实验台进行了管道消声的数值模拟,得到了远场声散射结果,进行了声衬声阻抗影响的计算,并通过与文献结果的对比表明该方法具有工程应用价值.      

一种基于谱方法的三维热声不稳定模型

研究流体的流动参数,包括热源热释放的周向不均匀分布对热声不稳定现象的影响,将欧拉方程进行线化处理之后,对轴向和径向方向进行基于切比雪夫多项式的谱离散,周向方向则进行基于傅里叶多项式的谱离散,建立了一种三维热声不稳定数值模型.通过对圆柱管道共振频率以及声场的求解,进行了一系列算例的验证,证实了该方法的可靠性,并且通过热释...     

考虑变截面影响的航空发动机短舱声学模型及数值结果

针对航空发动机短舱声学设计的需要,发展一种快速计算方法——传递单元法(TEM),并与直接边界元方法相结合,计算发动机远场声辐射。TEM以等价分布源方法为基础,该方法的优势在于易于计算变截面管道内的声传播,同时,对管道内部声衬位置及长度没有任何限制。给出更接近于真实发动机形状的声辐射算例。     

非紧致气动噪声半空间传播基本解的边界积分方法

为了获得低马赫数流动诱发的非紧致气动噪声在半空间内传播的基本解,结合复等效源方法和边界元方法建立了半空间精确格林函数的边界积分方程,当半空间边界为刚度型阻抗边界时可避免奇异性积分。同时基于等效源方法提出一种半空间二维非紧致圆柱声散射模型,推导了静止介质中声散射基本解的理论表达式。对静止介质中的二维圆柱声散射,数值解在研究的频率与观察点处与理论解一致。采用数值方法计算二维圆柱绕流诱发的半空间声场基本解,结果显示半空间边界强化了声源的声辐射能力,小于马赫数0.2的介质运动对声传播的影响可以忽略。     

基于KZK方程的声爆频域预测法

建立了一种快速预测声爆传播特性的频域方法,基于传统Khohklov-Zabolotskaya-Kuznetsov (KZK)方程,描述声爆沿激波波阵面法线方向的传播。为了验证模型正确性,以NASA TD N-161中试件C为对象,首先基于不同网格,利用计算流体力学(CFD)方法得到超音速流场;将CFD结果作傅立叶变换后代入模型方程,快速求解声爆传播特性。预测结果与NASA实验结果符合很好,研究表明：预测方法能够捕捉声爆的非线性传播,声衍射项使得声场在远场趋于轴对称分布,远距离传播后声能量集中于低频分量。     

基于Khokhlov-Zabolotskaya-Kuznetsov方程的声爆频域预测法

建立了一种快速预测声爆传播特性的频域方法,基于传统Khokhlov-Zabolotskaya-Kuznetsov(KZK)方程,描述声爆沿激波波阵面法线方向的传播.为了验证模型正确性,以NASA TD N-161中试件C为对象,首先基于不同网格,利用计算流体动力学(CFD)方法得到超声速流场;将CFD结果作傅里叶变换后代入模型方程,快速求解声爆传播特性.预测结果与NASA实验结果符合很好,研究表明:预测方法能够捕捉声爆的非线性传播,声衍射项使得声场在远场趋于轴对称分布,远距离传播后声能量集中于低频分量.      

Transfer element method with application to acoustic design of aeroengine nacelle

In the present survey, various methods for the acoustic design of aeroengine nacelle are first briefly introduced along with the comments on their advantages and disadvantages for practical application, and then detailed analysis and discussion focus on a kind of new method which is called ‘‘transfer element method'(TEM) with emphasis on its application in the following three problems: turbomachinery noise generations, sound transmission in ducts and radiation from the inlet and outlet of ducts, as well as the interaction between them. In the theoretical frame of the TEM, the solution of acoustic field in an infinite duct with stator sound source or liner is extended to that in a finite domain with all knows and unknowns on the interface plane, and the relevant acoustic field is solved by setting up matching equation. In addition, based on combining the TEM with the boundary element method(BEM) by establishing the pressure and its derivative continuum conditions on the inlet and outlet surface, the sound radiation from the inlet and outlet of ducts can also be investigated. Finally, the effects of various interactions between the sound source and acoustic treatment have been discussed in this survey. The numerical examples indicate that it is quite important to consider the effect of such interactions on sound attenuation during the acoustic design of aeroengine nacelle.     

一种预测平板尾迹噪声的时域无网格方法

基于离散涡方法和涡声理论建立了一种预测二维平板尾迹发声的时域无网格方法。该方法应用解耦方式完成声场计算,首先使用离散涡方法计算了均匀来流作用下的平板尾迹流场,得到了流场中点涡的涡量、位置和速度等关键参数,然后基于涡声理论建立了自由空间中点涡发声模型,并引入了时域边界元方法来模拟平板表面对声场的散射作用,计算得到了平板尾迹涡发声的偶极子声场分布和指向性等关键特征。通过对上下表面涡以及平板散射对声场贡献的深入分析表明,从平板尾缘上下角点脱落并卷起的涡团均为偶极子源,平板的散射作用使得声场在一定程度上得到加强,并且使声场具有极大值方向垂直于平板表面的偶极子指向性特征。所建立的无网格方法计算快速,能同时获得流场和声场分布的关键特征,可提升对气动噪声产生机理的基本认识,同时还为尾迹噪声的理论研究提供了一种具有工程应用价值的可靠计算方法。     

燃气轮机叶轮非线性瞬态温度场的边界元分析

首先采用克希荷夫变换将包含变物性及混合和辐射边界条件的热扩散方程转换成线性方程，进一步由加权余量法和离散化将方程转变成边界元方程组。用所开发的二维线性—非线性瞬态温度场边界元程序分析计算了燃气轮机叶轮起动工况下的非稳态温度场，与有限元数值计算结果比较表明，边界元法应用于这类非线性工程问题的分析是可靠、初始精度高，数值稳定性好及计算输入数据简单等优点。     

压气机轮盘孔边应力的边界元分析

 本文根据压气机轮盘的轴对称问题,取出轮盘包含孔的一个扇形部分,以有限元法计算的结果作为扇形部分的边界条件,采用边界元法计算轮盘孔边的应力。这与直接用有限元法计算结果对比表明,边界元法具有数据准备工作量小、方便,计算点选择灵活,计算精度高等优点,用于计算孔边的应力集中问题是有效的。文中还推导丁常单元的二阶差分逼近式,用以求解边界上的应力。