快速随机粒子网格法的气动噪声预测方法

耦合随机湍流速度生成模型与线化欧拉方程技术,形成了一套具备模拟噪声在非均匀流场中传播能力的气动噪声混合预测方法。该混合方法的随机湍流速度生成模型采用了快速随机粒子网格法,为声传播模拟提供了可靠的源项。而噪声的传播计算选用线化欧拉方程,其空间离散采用9点5阶的色散保持关系格式,时间推进选用了高精度大时间步长的6级4阶龙格库塔格式,远场边界应用了无分裂形式的理想匹配层边界条件。首先,选用高斯脉冲传播算例对线化欧拉方程的时空离散格式、远场无反射边界条件进行了验证分析。然后,计算分析各向同性湍流的空间相关性验证湍流速度生成模型的可靠性。最后,基于已搭建的气动噪声混合预测方法进行了30P30N三段翼缝翼噪声的计算分析。计算分析可知:监测点处功率谱密度曲线、噪声指向性等计算结果与参考文献结果取得了较好的一致性。数值计算结果表明所建立的气动噪声混合预测方法能有效预测二维复杂构型的气动噪声问题。      

NACA-65-0012机翼的非定常反问题解

利用有限差分法解非定常欧拉方程是振荡机翼非定常气动反问题的一种求解方法。本文重点讨论了非定常气动反问题求解的计算方案与反问题固壁边界条件的处理,编写了相应的计算程序。本文以绕前缘点做小幅俯仰周期振动的NACA-65-0012翼型为初始机翼,在振幅分别为0°、2．51°的条件下进行了气动反问题计算,经过多次迭代后计算翼型的压力分布基本与目标压力吻合,这表明本文的反问题求解方案、反问题边界处理及计算程序是可行的。     

数值仿真生成的汽车后空调气动噪声预测及评价

为解决需要项目开发后期样车完成后才能进行整车空调气动噪声性能测试及评价的滞后性问题,针对某SUV整车后空调高档范围工作产生的气动噪声,提出了一种基于空调整车计算流体动力学和FW-H声类比方程的气动噪声仿真计算分析和神经网络的主观评价预测方法。首先采用计算流体动力学和FW-H声类比相结合的方法仿真计算和验证后空调高档运行时车内产生的气动噪声特性;然后将仿真得到的时域气动噪声样本转化成声频格式,并开发GUI程序界面进行噪声样本主观评价和客观参数计算;最后建立基于遗传算法优化的主客观映射神经网络预测模型以实现车内后空调气动噪声性能的预测评价。仿真计算及预测评价结果表明：该方法计算的气动噪声仿真误差在10%以内,主观预测误差在0.5分以内,可有效指导汽车空调气动噪声性能的前瞻性预测开发。     

应用控制理论的叶栅气动反设计

基于控制理论的气动设计方法作为一种基于梯度的优化方法,通过引入伴随系统计算目标函数的敏感性导数,大大降低设计成本.本文将基于控制理论的气动设计方法应用到透平叶栅的气动反问题中,应用Euler方程研究了二维叶栅的压力反设计问题,并讨论了该方法具体实施中的关键问题,包括采用非均匀B样条进行二维叶栅造型;应用Thompson时间相关边界条件理论进行伴随方程特征分析;研究伴随方程的数值求解方法,构造伴随方程的耗散通量.通过算例证明了该气动设计方法适用性好,速度快,可以大大节约计算成本.     

Fourier拟谱配点方法对第一类边界积分方程的应用

本文用Fourier拟谱配点方法求解有广泛应用的以对数核为主部的第一类边界积分方程,文中通过对积分算子的象征作拟谱插值来建立近似方程,利用快速Fourier变换将计算切换到频率空间进行。本文计算结果表明,用上述拟谱配点方法计算的数值精度较Galerkin配点法更为满意。     

对流占优扩散问题的高精度直线法

基于常微分方程边值问题的高精度求解器SEVORD对偏微分方程作半离散,提出了求解一维对流扩散方程的高精度直线法,并采用局部一维化方法(LOD)给出了求解二维对流扩散问题的高精度交替方向直线法。     

基于气动反问题的风力机翼型优化设计

翼型气动优化设计的核心在于快速、准确的流动性能分析与快速、可靠的寻优算法。本文提出通过控制翼型表面预期流动分布,应用气动反问题方法,求解性能优化的翼型气动设计新方法。本文的翼型流动分析是基于位势流动与边界层积分方程的迭代解法。气动参数寻优采用了基于多变量搜索的加速POWELL算法,在确定的参数空间内,遍历搜索最佳性能点,可以保证最终优化解的全局性。气动反问题求解采用了壁面喷吸气模型。计算案例表明,本文方法的计算时间少,在流动不分离时具有与RANS同等计算精度,为快速开发低速风力机翼型提供了一个有效的设计方法。     

重心Lagrange插值配点法求解二维双曲电报方程

提出一种求解二维双曲电报方程的高精度重心Lagrange插值配点法.采用重心Lagrange插值构造包含时间和空间变量的近似函数.在给定Chebyshev-Gauss-Lobatto节点上,将多变量重心Lagrange插值近似函数代入双曲电报方程及其定解条件,得到离散代数方程组.包含狄里克雷和诺依曼边界条件的数值算例表明,本文方法程序实现方便并具有高精度,可应用于求解高维问题.     

非线性Schrdinger方程的直接间断Galerkin方法

讨论一维和二维非线性Schrdinger(NLS)方程的数值求解.基于扩散广义黎曼问题的数值流量,构造一种直接间断Galerkin方法(DDG)求解非线性Schrdinger方程.证明该方法L2稳定性,并说明DDG格式是一种守恒的数值格式.对一维NLS方程的计算表明,DDG格式能够模拟各种孤立子形态,而且可以保持长时间的高精度.二维NLS方程的数值结果显示该方法的高精度和捕捉大梯度的能力.     

A Direct Discontinuous Galerkin Method for Nonlinear Schr(o)dinger Equation

讨论一维和二维非线性Schr(o)dinger (NLS)方程的数值求解.基于扩散广义黎曼问题的数值流量,构造一种直接间断Galerkin方法(DDG)求解非线性Schr(o)dinger方程.证明该方法L2稳定性,并说明DDG格式是一种守恒的数值格式.对一维NLS方程的计算表明,DDG格式能够模拟各种孤立子形态,而且可以保持长时间的高精度.二维NLS方程的数值结果显示该方法的高精度和捕捉大梯度的能力.     

轴对称体声振耦合的边界子波谱与有限元耦合方法

探讨了子波在Helmholtz积分方程及声振耦合中的应用,在建立了求解轴对称Helmholtz积分方程的子波谱方法的基础上,构造了轴对称子波谱与轴对称有限元的耦合方法,该方法可以处理轴对称问题的任意边界条件.进行了声振耦合问题的模态分析.     

Fast prediction of aerodynamic noise induced by the flow around a cylinder based on deep neural network

Accurate and fast prediction of aerodynamic noise has always been a research hotspot in fluid mechanics and aeroacoustics. The conventional prediction methods based on numerical simulation often demand huge computational resources, which are difficult to balance between accuracy and efficiency. Here, we present a data-driven deep neural network (DNN) method to realize fast aerodynamic noise prediction while maintaining accuracy. The proposed deep learning method can predict the spatial distributions of aerodynamic noise information under different working conditions. Based on the large eddy simulation turbulence model and the Ffowcs Williams-Hawkings acoustic analogy theory, a dataset composed of 1216 samples is established. With reference to the deep learning method, a DNN framework is proposed to map the relationship between spatial coordinates, inlet velocity and overall sound pressure level. The root-mean-square-errors of prediction are below 0.82 dB in the test dataset, and the directivity of aerodynamic noise predicted by the DNN framework are basically consistent with the numerical simulation. This work paves a novel way for fast prediction of aerodynamic noise with high accuracy and has application potential in acoustic field prediction.     

Helicopter noise in hover: Computational modelling and experimental validation

The aeroacoustic characteristics of a helicopter rotor are calculated by a new method, to assess its applicability in assessing rotor performance in hovering. Direct solution of the Euler equations in a noninertial coordinate system is used to calculate the near-field flow around the spinning rotor. The far-field noise field is calculated by the Ffowcs Williams–Hawkings (FW–H) method using permeable control surfaces that include the blade. For a multiblade rotor, the signal obtained is duplicated and shifted in phase for each successive blade. By that means, the spectral characteristics of the far-field noise may be obtained. To determine the integral aerodynamic characteristics of the rotor, software is written to calculate the thrust and torque characteristics from the near-field flow solution. The results of numerical simulation are compared with experimental acoustic and aerodynamic data for a large-scale model of a helicopter main rotor in an open test facility. Two- and four-blade configurations of the rotor are considered, in different hover conditions. The proposed method satisfactorily predicts the aerodynamic characteristics of the blades in such conditions and gives good estimates for the first harmonics of the noise. That permits the practical use of the proposed method, not only for hovering but also for forward flight.     

Prediction of flow-generated noise produced by acoustic and aerodynamic interactions of multiple in-duct elements

Flow-generated noise problem caused by in-duct elements is due to the complicated acoustic and turbulent interactions of multiple in-duct flow noise sources. The approach of partially coherent sound fields used previously by Mak and Yang [C.M. Mak, J. Yang, Flow-generated noise radiated by the interaction of two strip spoilers in a low speed flow ducts, Acta Acust united with Acustica 88 (2002) 861–868] and Mak [C.M. Mak, A prediction method for aerodynamic sound produced by multiple elements in air ducts, J Sound Vib 287 (2005) 395–403] is adopted to formulate the sound powers produced by interactions of multiple elements at frequencies below and above the cut-on frequency of the lowest transverse duct mode. The study indicates that the level and spectral distribution of the additional acoustic energy produced by the interactions of multiple elements can be predicted based on the measured data with respect to the interactions. The proposed method can form a basis of a generalized prediction method for flow-generated noise produced by multiple elements. The application of the proposed method is supported by two engineering examples.     

Automotive panel noise contribution modeling based on finite element and measured structural-acoustic spectra

The radiated noise contributions of automotive body panels to the interior sound pressure levels are modeled using an approximate spectral formulation that applies the theoretical interior acoustic sensitivity terms derived from a finite element model and measured spatial-averaged structural-acoustic spectra. The finite element calculation is validated by comparison to a set of experimental acoustic transfer functions. A measurement set-up for the sound intensity and structural-acoustic response is applied to acquire the cross and auto power spectra needed to predict the relative mean-squared velocity term of each control plane near the panel surface, and to obtain the individual panel contribution function. The proposed approach also computes the noise spectra in 1/12 octave band form at selected positions in the passenger compartment, which matches well with the overall experimental results. Through an actual passenger car application, the approximate computational scheme is proven to be generally quite robust and effective for analyzing higher frequency interior noise problems.     

多网格谱元法及其在高性能计算中的应用*

针对传统谱元法在每个单元内只能存在单一均匀介质,应用在复杂非均匀介质的波传播模拟中可能造成极大计算规模的问题,发展了多网格谱元法。该方法在谱元法单元内引入独立的辅助网格,用于精细描述单元内的介质和外力分布变化,在较稀疏的主网格上进行波场的求解。基于声波和弹性波方程推导了多网格谱元法公式,并对几种典型模型进行了波场的数值模拟。与传统谱元法的对比结果表明,此算法在复杂非均匀介质的弹性波传播模拟中可以利用较少的网格点数达到不低于传统算法的精度。此外,实现了并行化的多网格谱元法,获得了较好的并行效率。     

Computation of acoustic radiation from vibrating structures in motion

The noises from a vibrating structure in motion are often encountered in engineering practice, for example, tire noise and pass-by noise of moving vehicles. Consequently, research on the radiation characteristics of moving acoustic source is of significance. In this paper, a new computational method based on the wave superposition approach is developed for the acoustic field from a vibrating structure in motion. It inherits the advantages of the wave superposition approach in the acoustic computation, and in which a method of moving simple sources is used to eliminate the influence of the Doppler effect. By the proposed method, the acoustic radiation from the moving vibrating structure can be calculated easily with the same implementation process as the conventional wave superposition approach performed in the stationary acoustic field. Finally, the validity and feasibility of the proposed method are verified by the numerical results.     

Total annoyance from an industrial noise source with a main spectral component combined with a background noise

When living close to an industrial plant, people are exposed to a combination of industrial noise sources and a background noise composed of all the other noise sources in the environment. As a first step, noise annoyance indicators in laboratory conditions are proposed for a single exposure to an industrial noise source. The second step detailed in this paper involves determining total annoyance indicators in laboratory conditions for ambient noises composed of an industrial noise source and a background noise. Two types of steady and permanent industrial noise sources are studied: low frequency noises with a main spectral component at 100 Hz, and noises with a main spectral component in middle frequencies. Five background noises are assessed so as to take into account different sound environments which can usually be heard by people living around an industrial plant. One main conclusion of this study is that two different analyses are necessary to determine total annoyance indicators for this type of ambient noise, depending on the industrial noise source composing it. Therefore, two total annoyance indicators adapted to the ambient noises studied are proposed.     

A three-dimensional hybrid finite element/spectral analysis of noise radiation from turbofan inlets

This paper describes a new three-dimensional (3D) analysis of tonal noise radiated from non-axisymmetric turbofan inlets. The novelty of the method is in combining a standard finite element discretisation of the acoustic field in the axial and radial coordinates with a Fourier spectral representation in the circumferential direction. The boundary conditions at the farfield, fan face and acoustic liners are treated using the same spectral representation. The resulting set of discrete acoustic equations are solved employing the well-established BICGSTAB or QMR iterative algorithms and a very effective specialised preconditioner based on the axisymmetric mean geometry and flow field. Numerical examples demonstrate the suitability of the new method to engine configurations with realistic 3D features, such as relatively large degrees of asymmetry and spliced acoustic liners. The examples also illustrate the two advantages of the new method over a traditional 3D finite element approach. The new method requires a significantly smaller number of unknowns as relatively few circumferential Fourier modes in the spectral solution ensure an accurate field representation. Also, due to the effective preconditioner, the spectral linear solver benefits from stable iterations at a high rate of convergence.     

Extraction of the acoustic component of a turbulent flow exciting a plate by inverting the vibration problem

An improvement of the Force Analysis Technique (FAT), an inverse method of vibration, is proposed to identify the low wavenumbers including the acoustic component of a turbulent flow that excites a plate. This method is a significant progress since the usual techniques of measurements with flush-mounted sensors are not able to separate the acoustic and the aerodynamic energies of the excitation because the aerodynamic component is too high. Moreover, the main cause of vibration or acoustic radiation of the structure might be due to the acoustic part by a phenomenon of spatial coincidence between the acoustic wavelengths and those of the plate. This underlines the need to extract the acoustic part. In this work, numerical experiments are performed to solve both the direct and inverse problems of vibration. The excitation is a turbulent boundary layer and combines the pressure field of the Corcos model and a diffuse acoustic field. These pressures are obtained by a synthesis method based on the Cholesky decomposition of the cross-spectra matrices and are used to excite a plate. Thus, the application of the inverse problem FAT that requires only the vibration data shows that the method is able to identify and to isolate the acoustic part of the excitation. Indeed, the discretization of the inverse operator (motion equation of the plate) acts as a low-pass wavenumber filter. In addition, this method is simple to implement because it can be applied locally (no need to know the boundary conditions), and measurements can be carried out on the opposite side of the plate without affecting the flow. Finally, an improvement of FAT is proposed. It regularizes optimally and automatically the inverse problem by analyzing the mean quadratic pressure of the reconstructed force distribution. This optimized FAT, in the case of the turbulent flow, has the advantage of measuring the acoustic component up to higher frequencies even in the presence of noise.