Multidimensional optimization of finite difference schemes for Computational Aeroacoustics

Because of the long propagation distances, Computational Aeroacoustics schemes must propagate the waves at the correct wave speeds and lower the isotropy error as much as possible. The spatial differencing schemes are most frequently analyzed and optimized for one-dimensional test cases. Therefore, in multidimensional problems such optimized schemes may not have isotropic behavior. In this work, optimized finite difference schemes for multidimensional Computational Aeroacoustics are derived which are designed to have improved isotropy compared to existing schemes. The derivation is performed based on both Taylor series expansion and Fourier analysis. Various explicit centered finite difference schemes and the associated boundary stencils have been derived and analyzed. The isotropy corrector factor, a parameter of the schemes, can be determined by minimizing the integrated error between the phase or group velocities on different spatial directions. The order of accuracy of the optimized schemes is the same as that of the classical schemes, the advantage being in reducing the isotropy error. The present schemes are restricted to equally-spaced Cartesian grids, so the generalized curvilinear transformation method and Cartesian grid methods are good candidates. The optimized schemes are tested by solving various multidimensional problems of Aeroacoustics.     

Direct and inverse scattering of transient acoustic waves by a slab of rigid porous material

This paper provides a temporal model of the direct and inverse scattering problem for the propagation of transient ultrasonic waves in a homogeneous isotropic slab of porous material having a rigid frame. This new time domain model of wave propagation takes into account the viscous and thermal losses of the medium as described by the model of Johnson et al. [D. L. Johnson, J. Koplik, and R. Dashen, J. Fluid. Mech. 176, 379 (1987)] and Allard [J. F. Allard (Chapman and Hall, London, 1993)] modified by a fractional calculus based method applied in the time domain. This paper is devoted to the analytical calculus of acoustic field in a slab of porous material. The main result is the derivation of the expression of the scattering operators (reflection and transmission) which are the responses of the medium to an incident acoustic pulse. In this model the reflection operator is the sum of two contributions: the first interface and the bulk of the medium. Experimental and numerical results are given as a validation of our model.     

基于积分方法的二维平行剪切层声远场预测

本文采用二维Ffowcs Williams＆Hawkings(FW-H)方程对平行剪切层远声场辐射特性进行了研究。近流场时间精确数据通过计算气动声学(Computational Aeroacoustics,CAA)技术数值模拟获得,声远场信息则通过FW-H方程对近流场内的可穿透积分面进行积分获得。该方法首先采用具有解析解的涡/尾缘干涉问题进行了校核,进一步采用CAA/FW-H匹配技术对二维平行剪切层声辐射问题进行了预测,计算结果表明,积分解与计算域内的CAA数值解吻合较好。     

三维FW-H方程与CAA数值模拟匹配技术研究

本文研究了三维FW-H方程与CAA数值模拟匹配技术。首先验证了适于亚音任意运动声源的Farassat时域公式和适于亚音匀速直线运动声源的Lockard频域公式,采用均匀流中单极子源声辐射问题对两类公式进行了校核。进一步,采用FW-H／CAA匹配技术对风扇／压气机前传声进行了预测。近声场基于轴对称三维CAA方法获得,远声场则基于近声场数据采用三维可穿透FW-H方程时域公式进行预测,并校核算例着重分析了不同积分面对远场声指向性的影响。本文研究证实了FW-H／CAA数值模拟匹配技术的可行性和解决工程实际问题的潜力。     

Coupled hydrodynamic-acoustic modeling of sound generated by impacting cylindrical water jets

A coupled hydrodynamic-acoustic model describing acoustic propagation in a fluid containing multiple bubbles is proposed and applied to simulate noise generated by impacting water jets. The total pressure is decomposed into a "hydrodynamic" part and an "acoustic" part and computed using different schemes. The hydrodynamic pressure field is calculated independently using a generalized hydrodynamic model, and the pressure variations serve as sources in the wave equation for the acoustic pressure. A numerical algorithm developed to calculate wave propagation in an irregular region is used to account for the existence of the cavities. Noise generated by the impact of two cylindrical water jets is predicted. The computed near-field pressure is compared with the experimental data.     

A three-dimensional coupled-mode model for the acoustic field in a two-dimensional waveguide with perfectly reflecting boundaries

This paper presents a three-dimensional(3D) coupled-mode model using the direct-global-matrix technique as well as Fourier synthesis. This model is a full wave, two-way three-dimensional model, and is therefore capable of providing accurate acoustic field solutions. Because the problem of sound propagation excited by a point source in an ideal wedge with perfectly reflecting boundaries is one of a few three-dimensional problems with analytical solutions, the ideal wedge problem is chosen in this work to validate the presented three-dimensional model. Numerical results show that the field results by analytical solutions and those by the presented model are in excellent agreement, indicating that the presented model can serve as a benchmark model for three-dimensional sound propagation problems involving a planar two-dimensional geometry as well as a point source.     

Discrete-time modeling of woodwind instrument bores using wave variables

A method for simulation of acoustical bores, useful in the context of sound synthesis by physical modeling of woodwind instruments, is presented. As with previously developed methods, such as digital waveguide modeling (DWM) [Smith, Comput. Music J. 16, 74-91 (1992)] and the multi convolution algorithm (MCA) [Martinez et al., J. Acoust. Soc. Am. 84, 1620-1627 (1988)], the approach is based on a one-dimensional model of wave propagation in the bore. Both the DWM method and the MCA explicitly compute the transmission and reflection of wave variables that represent actual traveling pressure waves. The method presented in this report, the wave digital modeling (WDM) method, avoids the typical limitations associated with these methods by using a more general definition of the wave variables. An efficient and spatially modular discrete-time model is constructed from the digital representations of elemental bore units such as cylindrical sections, conical sections, and toneholes. Frequency-dependent phenomena, such as boundary losses, are approximated with digital filters. The stability of a simulation of a complete acoustic bore is investigated empirically. Results of the simulation of a full clarinet show that a very good concordance with classic transmission-line theory is obtained.     

流管实验装置中声传播计算的模态方法

流管实验装置是测量有流动情况下航空发动机消声短舱内声衬声阻抗的主要装置。本文发展了一种解析的模态匹配方法进行在平均流有声衬条件下矩形流管中声传播的计算。用同伦方法求解特征值问题,并与用环绕积分求解的结果进行比较。声场通过轴向阻抗间断面的声压和声质点速度积分相等计算。第一个算例是无流动、硬壁、有限长、考虑端口反射的情况,并与北航流管实验台测量数据进行了对比;第二个算例为有流动情况下有限长声衬管道不考虑端口反射的声场计算,它与文献中NASA流管实验结果和CAA计算结果符合得很好。     

Benchmark solutions for sound propagation in an ideal wedge

Sound propagation in a wedge-shaped waveguide with perfectly reflecting boundaries is one of the few range-dependent problems with an analytical solution, and hence provides an ideal benchmark for a full two-way solution to the wave equation. An analytical solution for the sound propagation in an ideal wedge with a pressure-release bottom was presented by Buckingham and Tolstoy [Buckingham and Tolstoy 1990 J. Acoust. Soc. Am. 87 1511]. The ideal wedge problem with a rigid bottom is also of great importance in underwater acoustics. We present an analytical solution to the ideal wedge problem with a perfectly reflecting bottom, either rigid or pressure-release, which may be used to provide a means for investigating the sound field in depth-varying channels, and to establish the accuracy of numerical propagation models. Closed-form expressions for coupling matrices are also provided for the ideal waveguides characterized by a homogeneous water column bounded by perfectly reflecting boundaries. A comparison between the analytical solution and the numerical solution recently proposed by Luo et al. [Luo W Y, Yang C M and Zhang R H 2012 Chin. Phys. Lett. 29 014302] is also presented, through which the accuracy of this numerical model is illustrated.     

Sound Attenuation in a Circular Duct of a Viscous Medium in the Absence of Mean Flow

Abstract

An analytical method to study the effect of viscosity of a medium and the wave number on sound propagation and sound attenuation numbers in circular ducts has been presented. The method is based on the variation of parameters of the solution corresponding to the case of inviscid acoustic waves in circular ducts and axisymmetric modes. A mathematical model is constructed to describe the physical problem in general. Three basic assumptions have been considered, namely, each flow quantity has been written as the sum of a steady mean flow and an unsteady acoustic flow quantity. The effect of thermal conductivity of the gas has been neglected as well as no mean flow. The results for a wide range of wave numbers and Reynolds numbers show that for a viscous medium, the propagation number is a weak function of the Reynolds number, and as the Reynolds number increases, the propagation number approaches its inviscid value. Also the propagation number is independent of the wave number. For the attenuation number, it decreases monotonically with the increase of the Reynolds number and it vanishes when Reynolds number exceeds 10⁴.     

Direct computation of the noise radiated by a subsonic cavity flow and application of integral methods

The goal of this paper is to investigate the acoustic field generated by the flow over a cavity using two different and complementary numerical methods. First, a Direct Numerical Simulation of the 2-D compressible Navier-Stokes equations is performed to obtain directly the radiated noise. The results of the acoustic and aerodynamic fields are compared to the experimental data in the literature. Second, this reference solution is compared to solutions provided by hybrid methods using the flowfield computed inside the cavity combined with an integral formulation to evaluate the far-field noise. Numerical issues of three integral methods are studied: the Ffowcs Williams and Hawkings analogy that extends Lighthill's theory to account for solid boundaries and two Wave Extrapolation Methods from a control surface, the Kirchhoff and porous Ffowcs Williams and Hawkings methods. All methods show a good agreement with the Direct Numerical Simulation, but the first one is more expensive owing to an additional volume integral. However, the analogy can help in the analysis of wave patterns, by separating the direct waves from the reflected ones. The wave extrapolation methods from a surface are more efficient and provide a complementary tool to extend Computational Aeroacoustics near field to the very far field.     

非线性对大气介质中阵列聚焦声场分布影响的研究

基于时域有限差分算法将大气中近似到二阶微小项的非线性声波波动方程进行离散化,得到了模拟采用的差分波动方程.在此基础上,数值模拟了初始声压强弱不同的5个点声源组成的线阵列垂直或斜向辐射的连续正弦波在大气中传播时二维声场的分布情况.将线性条件下的模拟结果与非线性条件下的模拟结果进行比较后发现:弱非线性会对声场的分布和阵列聚焦增益产生一定的影响,使声场分布波形比线性条件下的声场分布波形更加靠近阵列,聚焦效果变差;强非线性会使波形发生更严重畸变,这是由于产生了基频以外的其他频率声波引起的;非线性对斜向传播时声场分布的影响与垂直传播时的影响效果基本相同,但由于斜向辐射时的声波几何扩展造成的轴向声压衰减要大于垂直辐射时的轴向声压衰减,因此聚焦增益和强非线性的影响都将小于垂直辐射时的情况.     

Mathematical model for characterizing noise transmission into finite cylindrical structures

This work presents a theoretical study of the sound transmission into a finite cylinder under coupled structural and acoustic vibration. Particular attention of this study is focused on evaluating a dimensionless quantity, "noise reduction," for characterizing noise transmission into a small cylindrical enclosure. An analytical expression of the exterior sound pressure resulting from an oblique plane wave impinging upon the cylindrical shell is first presented, which is approximated from the exterior sound pressure for an infinite cylindrical structure. Next, the analytical solution of the interior sound pressure is computed using modal-interaction theory for the coupled structural acoustic system. These results are then used to derive the analytical formula for the noise reduction. Finally, the model is used to predict and characterize the sound transmission into a ChamberCore cylindrical structure, and the results are compared with experimental data. The effects of incidence angle and internal acoustic damping on the sound transmission into the cylinder are also parametrically studied.     

A new combined stable and dispersion relation preserving compact scheme for non-periodic problems

A new compact scheme is presented for computing wave propagation problems and Navier–Stokes equation. A combined compact difference scheme is developed for non-periodic problems (called NCCD henceforth) that simultaneously evaluates first and second derivatives, improving an existing combined compact difference (CCD) scheme. Following the methodologies in Sengupta et al. [T.K. Sengupta, S.K. Sircar, A. Dipankar, High accuracy schemes for DNS and acoustics, J. Sci. Comput. 26 (2) (2006) 151–193], stability and dispersion relation preservation (DRP) property analysis is performed here for general CCD schemes for the first time, emphasizing their utility in uni- and bi-directional wave propagation problems – that is relevant to acoustic wave propagation problems. We highlight: (a) specific points in parameter space those give rise to least phase and dispersion errors for non-periodic wave problems; (b) the solution error of CCD/NCCD schemes in solving Stommel Ocean model (an elliptic p.d.e.) and (c) the effectiveness of the NCCD scheme in solving Navier–Stokes equation for the benchmark lid-driven cavity problem at high Reynolds numbers, showing that the present method is capable of providing very accurate solution using far fewer points as compared to existing solutions in the literature.     

Transverse plane wave analysis of short elliptical chamber mufflers: An analytical approach

Short elliptical chamber mufflers are used often in the modern day automotive exhaust systems. The acoustic analysis of such short chamber mufflers is facilitated by considering a transverse plane wave propagation model along the major axis up to the low frequency limit. The one dimensional differential equation governing the transverse plane wave propagation in such short chambers is solved using the segmentation approaches which are inherently numerical schemes, wherein the transfer matrix relating the upstream state variables to the downstream variables is obtained. Analytical solution of the transverse plane wave model used to analyze such short chambers has not been reported in the literature so far. This present work is thus an attempt to fill up this lacuna, whereby Frobenius solution of the differential equation governing the transverse plane wave propagation is obtained. By taking a sufficient number of terms of the infinite series, an approximate analytical solution so obtained shows good convergence up to about 1300 Hz and also covers most of the range of muffler dimensions used in practice. The transmission loss (TL) performance of the muffler configurations computed by this analytical approach agrees excellently with that computed by the Matrizant approach used earlier by the authors, thereby offering a faster and more elegant alternate method to analyze short elliptical muffler configurations.     

一种水平变化可穿透波导中声传播问题的耦合简正波方法

通过利用标准简正波程序KRAKEN计算本地简正波解及耦合矩阵, 进一步发展了求解水平变化波导中声场的全局矩阵耦合简正波方法(Luo et al., "A numerically stable coupled-mode formulation for acoustic propagation in range-dependent waveguides," Sci. China-Phys. Mech. Astron. 55, 572 (2012)), 使得该方法可以处理具有可穿透海底及随深度变化声速剖面等实际问题, 并提供声场的完全双向解. 本文还给出了双层波导中耦合矩阵的解析表达式, 并利用其验证了本方法中耦合矩阵数值算法的精度. 最后, 利用改善后的全局矩阵耦合简正波模型(DGMCM)计算了美国声学学会(ASA)提出的可穿透楔形波导标准问题, 将所得数值解与参考解比较, 结果表明DGMCM方法可以精确处理水平变化波导中声传播实际问题. 关键词： 耦合简正波理论 全局矩阵方法 可穿透楔形波导     

Numerical simulation of the nonlinear ultrasonic pressure wave propagation in a cavitating bubbly liquid inside a sonochemical reactor

We investigate the acoustic wave propagation in bubbly liquid inside a pilot sonochemical reactor which aims to produce antibacterial medical textile fabrics by coating the textile with ZnO or CuO nanoparticles. Computational models on acoustic propagation are developed in order to aid the design procedures. The acoustic pressure wave propagation in the sonoreactor is simulated by solving the Helmholtz equation using a meshless numerical method. The paper implements both the state-of-the-art linear model and a nonlinear wave propagation model recently introduced by Louisnard (2012), and presents a novel iterative solution procedure for the nonlinear propagation model which can be implemented using any numerical method and/or programming tool. Comparative results regarding both the linear and the nonlinear wave propagation are shown. Effects of bubble size distribution and bubble volume fraction on the acoustic wave propagation are discussed in detail. The simulations demonstrate that the nonlinear model successfully captures the realistic spatial distribution of the cavitation zones and the associated acoustic pressure amplitudes.     

Dispersion of acoustic waves in the plasma of a self-sustained gas discharge

The dispersion effects appearing during the propagation of acoustic waves through the plasma of a weakly ionized gas are studied. The main theoretical results are based on the equation of propagation of sound in the medium with the so-called Rayleigh energy release mechanism, which has been obtained earlier. Unlike the previous investigations, the problem of propagation of a perturbation from a source and not the problem of propagation of the initial perturbation is solved. In particular, the sources of an N-shaped shock wave and a wave in the form of a symmetrical step are analyzed in detail. It is shown that depending on the direction of wave propagation (along or across the electric field in a plasma), it degenerates either into a wave packet with a wave frequency lower than a certain frequency characterizing heating, or into a wave packet with a frequency higher than this value. In addition, a quantitative criterion is obtained, which makes it possible to estimate the plasma parameters for which it will be possible to observe the dispersion of acoustic waves in the plasma.