论文部分内容阅读
随着社会经济的快速发展,航空工业也将快速的发展,而飞机噪声是制约航空发展的一个至关重要的因素。现代涡扇发动机的风扇噪声在飞机起飞和降落过程中占主要成分,被动声衬处理目前仍是控制风扇噪声的主要技术手段。为了设计和优化飞机发动机进气道内的噪声传播问题,急需发展能够考虑声衬阻抗边界的计算气动声学计算框架。本文将变阶次高精度度保色散差分(DRP)方法和时域阻抗边界条件相结合,提出一套计及阻抗边界的时域多重网格计算气动声学框架,并运用于简单缓变截面管道以及复杂的涡扇发动机的进气道气动声传播过程的计算。利用Fourier-Laplace变换研究了显式保色散高阶高分辨率有限差分格式的性能,对DRP优化格式的积分区间选择的重要性做了清晰的阐述,进一步简化了积分区间自动选择的算法。详细地推导了大模板变截断精度显式差分格式的构造公式,提出了通用的格式简化优化算法,并比较了优化算法中的优化区间,优化目标函数等对优化格式的影响,最终得到一系列变精度DRP格式。随后采用1D对流谐波问题,2D高斯脉冲标准算例,以及一个相对比较复杂的多圆柱声散射问题验证了提出格式的精度和性能。这些格式将被用于后续的管道应用算例中。值得指出的是这些公式均借助现代计算代数系统(CAS)计算完成,在附录3中作者给出了相应的代码。这些代码对分析标准以及复杂的差分格式具有相当的通用性,经小的修改将可以用于众多差分格式的分析,可以大大减小格式分析的工作量和错误。基于滤波器理论对?zy?rük时域阻抗模型与实验数据拟合的初始值选择问题进行了详细的讨论,针对NASA CT73实验声衬的公开数据验证了该初始值选择的合理性,从而降低?zy?rük模型在直接拟合和使用中的难度,提高拟合的准确度,增进对模型参数对应滤波特性的物理理解。研究了扩展的Helmholtz共振器阻抗模型(EHRM)与物理模型参数之间的对应关系,论述了EHRM模型参数初选的基本原则,并采用非线性约束方法对模型进行直接拟合。基于余切函数的连续分式级数展开形式,找出EHRM模型的五个参数和?zy?rük模型多个参数之间的对应关系,从而建立了EHRM模型和?zy?rük间的对应关系,这为EHRM模型转换为?zy?rük提供了途径。由于EHRM模型转换为?zy?rük模型为中低频近似,在转换之后需要再进行一次非线性拟合,通过两次非线性约束优化的拟合方法,可以将EHRM模型转换为?zy?rük模型,并保持与实验数据的高精度拟合。最后,本文从参数拟合、实现难度以及计算效率方面对EHRM阻抗模型和?zy?rük进行了详细的对比。其中EHRM到?zy?rük格式的转换采用计算代数系统推导完成。将变精度显式差分格式和EHRM转换的?zy?rük模型用于Rienstra提出的缓变管道模型,与相应的解析解以及有限元解进行了详细的对比,验证了本文提出的DRP差分格式以及时域阻抗模型的精度和有效性,并简要的讨论了缓变截面管道声传播的特点。将论文提出的差分格式和时域阻抗模型用于实际的(有轴对称简化)高涵道比发动机进气道的噪声传播问题的模拟。探讨了管道中不同流速的时均流,不同频率的入射声波,不同的径向模态数以及有无阻抗等工况下气动噪声的传播特点,比较了气动声传播的远场特点。研究了开式的Helmholtz共振器共振频率以及频响的参数数值计算方法,以便运用于复杂及通用的Helmholtz共振器频率和响应的精确估计。针对考虑阻抗边界的管道气动声传播问题,本文在显式格式算法,时域阻抗边界的建模、拟合、实现、校验以及应用方面做了多方面的研究,进一步完善和推进基于DRP算法的CAA(Computational Aeroacoustics)在实际工程中的应用。