应用随机模型方法预测汽车风噪声

采用随机模型方法,对简化汽车头部外形进行风噪声数值模拟.将计算区域分为声源区域和传播区域,在声源区域采用随机模型构造湍流脉动速度场,传播区域通过求解带源项的线化欧拉方程实现声波向外传播得到声场解.同直接模拟方法相比,该方法具有计算量小、计算所需内存少等优点.数值模拟结果与实验数据吻合较好,验证了该方法预测汽车风噪声的可行性,为研究实际汽车外形的风噪声问题打下基础.     

可压缩Kelvin-Helmholtz不稳定性低耗散锐界面方法直接模拟

采用低耗散WENO(weighted essential non-oscillatory)格式及锐界面方法模拟可压缩Kelvin-Helmholtz不稳定性问题.由于物质界面被描述成一种接触间断, 该方法可精确求解切向速度间断.基于优化模板对原始光滑指标进行正规化后, 得到一种低耗散WENO格式.修正后的方法显著降低了普通流动区域的过衰减问题, 保持了良好的激波捕捉性能, 并可获得与混合格式相当的求解精度.不同于以往求解单一流体或易混界面时, 通过初始设定有限宽度的剪切层或快速数值耗散以抑制高波数模态, 该方法允许高波数扰动的发展.计算结果表明, 高波数扰动展现出与以往理想Kelvin-Helmholtz不稳定性问题数值模拟或线化理论结果不同的特征, 但与有限厚度的剪切层结果相符.      

液态金属快堆主管道中超声导波传播特性研究

讨论了钠冷快堆(Sodium-cooled Fast Reactor,SFR)主管道的整体温度和内部液态金属钠流动速度的变化对管道导波传播特性的影响。推导了充液管道中导波频散方程的一般形式,并给出了管道内液态金属钠处于流动状态下的导波频散方程。采用数值计算方法获得了管内液态金属钠处于不同温度和不同流速时的导波纵向模式频散曲线和导波时域波形。结果表明,温度变化对基阶纵向模式的影响较小,但对高阶纵向模式的影响较大;液态钠流速增大会使导波频散曲线向高频轻微移动,但在实际检测中可以忽路管内液体流动速度的影响。通过对时域接收波形的模拟计算,进一步考察了液态金属钠的温度及流动速度变化对导波传播的影响,并通过对比不同模态的激发特点和不同频段的导波时域波形特点,结合导波频散曲线,给出了适用于SFR管道超声无损检测的导波模态和声源激发频段选择方案。      

运动声源简正波模型稳相点的射线多普勒近似算法

从航空声源水下声场建模出发,提出了运动声源稳相点的近似算法。由于Hawker给出的运动声源声场简正波模型计算方法不适用于高速运动声源,针对速度相对较大的低空运动声源辐射噪声激发的水下声场建模,采用虚源概念,基于射线声学理论得到近似多普勒频率,结合简正波理论求解稳相点,较好地解决了高速运动声源声场的解算问题。几种波导环境下两种方法数值计算比较结果表明该算法的有效性。     

高精度有限差分法模拟Kelvin-Helmholtz不稳定性

 WENO有限差分格式有较高的分辨精度,适合复杂流场的计算,在国际上被广泛采用。本文利用WENO有限差分格式求解2维守恒型欧拉方程,实现了对无粘流体中Kelvin-Helmholtz不稳定性的数值模拟。速度剪切方向采用周期边界条件;扰动增长方向采用嵌边出流边界条件,一个不稳定波长分布64个网格。数值模拟给出的扰动幅值线性增长率与线性稳定性分析给出的结果很好符合,显示了该格式的有效性和精度。数值模拟给出了清晰的密度等值线,表明该方法还具有较好的界面变形捕捉能力。     

二维随机粗糙面上导体目标复合瞬态散射的混合算法

提出用时域积分方程法(TDIE)与时域基尔霍夫近似法(TDKA)的混合算法来求解二维导体随机粗糙面及其上方二维导体目标的复合瞬态散射,推导出了在TM波入射情形下显式及隐式格式的时间步进方程.将粗糙面与目标分别进行TDKA和TDIE计算,并考虑目标与粗糙面之间的耦合,对TDKA和TDIE进行混合迭代,既大大降低了粗糙面求解的复杂度,又保证了计算精度.数值算例中,考虑了角反射器(开放体)和圆柱(封闭体)两种目标,分别计算了目标表面电流响应和电场远场响应.计算结果表明,和单纯TDIE法相比,本文混合方法计算效率 关键词： 随机粗糙面 复合瞬态散射 时域积分方程法和时域基尔霍夫近似法 混合算法     

输油管道油流静电带电数值计算

对输油管道内油品流动带电问题的数值计算进行了研究.紊流条件下的电荷输运方程是一个对流占优的对流扩散反应方程,采用算子分裂法,将该方程分解为纯对流方程、纯扩散方程和纯反应方程,分别采用特征线法和差分法求解.算例证明,该方法能准确描述管道内电荷分布,因而提供了一种获取冲流电流的可靠方法.     

迎风紧致格式求解Hamilton-Jacobi方程

基于Hamilton-Jacobi(H-J)方程和双曲型守恒律之间的关系,将三阶和五阶迎风紧致格式推广应用于求解H-J方程,建立了高精度的H-J方程求解方法.给出了一维和二维典型数值算例的计算结果,其中包括一个平面激波作用下的Richtmyer Meshkov界面不稳定性问题.数值试验表明,在解的光滑区域该方法具有高精度,而在导数不连续的不光滑区域也获得了比较好的分辨效果.相比于同阶精度的WENO格式,本方法具有更小的数值耗散,从而有利于多尺度复杂流动的模拟中H-J方程的求解.     

自适应间断有限元方法求解三维欧拉方程

将龙格库塔间断有限元方法(RDDG)与自适应方法相结合,求解三维欧拉方程.区域剖分采用非结构四面体网格,依据数值解的变化采用自适应技术对网格进行局部加密或粗化,减少总体网格数目,提高计算效率.给出四种自适应策略并分析不同自适应策略的优缺点.数值算例表明方法的有效性.     

基于广义交替数值通量的局部间断Galerkin方法求解二维波动方程

波的传播往往在复杂的地质结构中进行,如何有效地求解非均匀介质中的波动方程一直是研究的热点.本文将局部间断Galekin(local discontinuous Galerkin, LDG)方法引入到数值求解波动方程中.首先引入辅助变量,将二阶波动方程写成一阶偏微分方程组,然后对相应的线性化波动方程和伴随方程构造间断Galerkin格式;为了保证离散格式满足能量守恒,在单元边界上选取广义交替数值通量,理论证明该方法满足能量守恒性.在时间离散上,采用指数积分因子方法,为了提高计算效率,应用Krylov子空间方法近似指数矩阵与向量的乘积.数值实验中给出了带有精确解的算例,验证了LDG方法的数值精度和能量守恒性;此外,也考虑了非均匀介质和复杂计算区域的计算,结果表明LDG方法适合模拟具有复杂结构和多尺度结构介质中的传播.     

Sound generation and transmission in flow ducts with axial temperature gradients

This article describes a one-dimensional, linearized, analysis of fundamental mode sound generation and propagation in rigid-walled flow ducts with axial temperature variation. An acoustic wave equation, including damping effects and volume sources, is derived and its solution (in the absence of sources) by a numerical technique and an approximate analytical method is discussed. The “forced” wave equation is then solved (the existence of an oscillating solution to the “unforced” equation being assumed) for sound generation by a side-branch volume source in an infinite duct, and the results are applied to a duct of finite length. Reasonably good agreement is obtained between measurements and predictions of the sound pressure field in a flow duct, away from the source region.     

Application of frequency-domain linearized Euler solutions to the prediction of aft fan tones and comparison with experimental measurements on model scale turbofan exhaust nozzles

Although it is widely accepted that aircraft noise needs to be further reduced, there is an equally important, on-going requirement to accurately predict the strengths of all the different aircraft noise sources, not only to ensure that a new aircraft is certifiable and can meet the ever more stringent local airport noise rules but also to prioritize and apply appropriate noise source reduction technologies at the design stage. As the bypass ratio of aircraft engines is increased - in order to reduce fuel consumption, emissions and jet mixing noise - the fan noise that radiates from the bypass exhaust nozzle is becoming one of the loudest engine sources, despite the large areas of acoustically absorptive treatment in the bypass duct. This paper addresses this ‘aft fan’ noise source, in particular the prediction of the propagation of fan noise through the bypass exhaust nozzle/jet exhaust flow and radiation out to the far-field observer. The proposed prediction method is equally applicable to fan tone and fan broadband noise (and also turbine and core noise) but here the method is validated with measured test data using simulated fan tones. The measured data had been previously acquired on two model scale turbofan engine exhausts with bypass and heated core flows typical of those found in a modern high bypass engine, but under static conditions (i.e. no flight simulation). The prediction method is based on frequency-domain solutions of the linearized Euler equations in conjunction with perfectly matched layer equations at the inlet and far-field boundaries using high-order finite differences. The discrete system of equations is inverted by the parallel sparse solver MUMPS. Far-field predictions are carried out by integrating Kirchhoff's formula in frequency domain. In addition to the acoustic modes excited and radiated, some non-acoustic waves within the cold stream-ambient shear layer are also captured by the computations at some flow and excitation frequencies. By extracting phase speed information from the near-field pressure solution, these non-acoustic waves are shown to be convective Kelvin-Helmholtz instability waves. Strouhal numbers computed along the shear layer, based on the local momentum thickness also confirm this in accordance with Michalke's instability criterion for incompressible round jets with a similar shear layer profile. Comparisons of the computed far-field results with the measured acoustic data reveal that, in general, the solver predicts the peak sound levels well when the farfield is dominated by the in-duct target mode (the target mode being the one specified to the in-duct mode generator). Calculations also show that the agreement can be considerably improved when the non-target modes are also included, despite their low in-duct levels. This is due to the fact that each duct mode has its own distinct directionality and a non-target low level mode may become dominant at angles where the higher-level target mode is directionally weak. The overall agreement between the computations and experiment strongly suggests that, at least for the range of mean flows and acoustic conditions considered, the physical aeroacoustic radiation processes are fully captured through the frequency-domain solutions to the linearized Euler equations and hence this could form the basis of a reliable aircraft noise prediction method.     

Unstable Surface Waves in Running Water

We consider the stability of periodic gravity free-surface water waves traveling downstream at a constant speed over a shear flow of finite depth. In case the free surface is flat, a sharp criterion of linear instability is established for a general class of shear flows with inflection points and the maximal unstable wave number is found. Comparison to the rigid-wall setting testifies that the free surface has a destabilizing effect. For a class of unstable shear flows, the bifurcation of nontrivial periodic traveling waves is demonstrated at all wave numbers. We show the linear instability of small nontrivial waves that appear after bifurcation at an unstable wave number of the background shear flow. The proof uses a new formulation of the linearized water-wave problem and a perturbation argument. An example of the background shear flow of unstable small-amplitude periodic traveling waves is constructed for an arbitrary vorticity strength and for an arbitrary depth, illustrating that vorticity has a subtle influence on the stability of free-surface water waves.     

ON TRANSMISSION OF SOUND IN A NON-UNIFORM DUCT CARRYING A SUBSONIC COMPRESSIBLE FLOW

The differential equations governing the transmission of one-dimensional sound waves in a non-uniform duct carrying a subsonic compressible mean flow have been the subject of a recent debate [1, 2]. Of the two formulations presented, one is considered to be non-acoustical and the other as neglecting the spatial variation of the speed of sound. The present paper shows that both formulations are acoustical and represent valid approximations to correct conditions for isentropic sound propagation in a subsonic low Mach number duct. Each formulation is associated with an “error wave”, which is essentially a hydrodynamic wave when the mean flow Mach number is small. Three-port modelling is required, however, to capture this wave when the Mach number of the mean flow is relatively large and a numerical matrizant approach is described which can be used for this purpose.     

The D'yakov-kontorovich instability of shock waves in real gases

In the 1950s, D'yakov and Kontorovich predicted that under certain conditions perturbed shock waves in nonideal gases can become unstable by emitting undamped sound and entropy-vortex waves. For the last 45 years, though, little progress has been made in the identification and numerical modeling of physical conditions for which this phenomenon might occur. Using a van der Waals equation of state, we present for the first time a dynamical simulation of a D'yakov-Kontorovich instability. The two-dimensional emission pattern of acoustic waves appearing in the simulation agrees with the prediction of a linearized theory.     

Secondary instabilities of interfacial waves due to coupling with a long wave mode in a two-layer Couette flow

We report experiments on the stability of interfacial waves in a two-layer Couette flow. As the shear rate is increased, the periodic wave train arising from the primary instability undergoes a secondary instability which results in wave coalescence or nucleation, after a long transient. This secondary instability crucially involves the coupling with a long wave mode, which corresponds to variations of the mean interface level. These observations are favourably compared to stability results on travelling wave solutions for a set of two coupled equations, one for the envelope of a weakly unstable wave packet, and the other for the marginal long wave mode with zero wave number. A physical mechanism for this instability is proposed, as well as an interpretation for the onset of chaos.     

Acoustical wave propagator for time-domain flexural waves in thin plates

In this paper, an explicit acoustical wave propagator technique is introduced to describe the time-domain evolution of acoustical waves in two-dimensional plates. A combined scheme with Chebyshev polynomial expansion and fast Fourier transformation is used to implement the operation of the acoustical wave propagator. Through this operation, the initial wave packet at t = 0 is mapped into the wave packet at any instant t > 0. By comparison of the results of the exact analytical solution and the Euler numerical method, we find that this new Chebyshev-Fourier scheme is highly accurate and computationally effective in predicting the acoustical wave propagation in thin plates. This method offers an opportunity for future study of dynamic stress concentration and time-domain energy flow in coupled structures.     

The wave equation with viscoelastic attenuation and its application in problems of shallow-sea acoustics

A suitable tool for the simulation of low frequency acoustic pulse signals propagating in a shallow sea is the numerical integration of the nonstationary wave equation. The main feature of such simulation problems is that in this case the sound waves propagate in the geoacoustic waveguide formed by the upper layers of the bottom and the water column. By this reason, the correct dependence of the attenuation of sound waves in the bottom on their frequency must be taken into account. In this paper we obtain an integro-differential equation for the sound waves in the viscoelastic fluid, which allows to simulate the arbitrary dependence of acoustic wave attenuation on frequency in the time domain computations. The procedure of numerical solution of this equation based on its approximation by a system of differential equations is then considered and the methods of artificial limitation of computational domain are described. We also construct a simple finite-difference scheme for the proposed equation suitable for the numerical solution of nonstationary problems arising in the shallow-sea acoustics.     

On the generation and maintenance of waves and turbulence in simulations of free-surface turbulence

One of the challenges in numerical simulation of wave–turbulence interaction is the precise setup and maintenance of wave and turbulence fields. In this paper, we investigate techniques for the generation and suppression of specific surface wave modes, the generation of turbulence in an inhomogeneous physical domain with a wavy boundary-fitted grid, and the generation and maintenance of waves and turbulence during the complex wave–turbulence interaction process. We apply surface pressure to generate and suppress waves. Based on the solution of linearized Cauchy–Poisson problem, we derive three pressure expressions, which lead to a δ-function method, a time-segment method, and a gradual method. Numerical experiments show that these methods generate waves as specified and eliminate spurious waves effectively. The nonlinear wave effect is accounted for with a time-relaxation method. For turbulence generation, we extend the linear forcing method to an inhomogeneous physical domain with a curvilinear computational grid. Effects of force distribution and computational grid distortion are examined. For wave–turbulence interaction, we develop an algorithm to instantaneously identify specific progressive and standing waves. To precisely control the wave amplitude in a complex turbulent flow field, we further develop an energy controlling method. Finally, a simulation example of wave–turbulence interaction is presented. Results show that turbulence has unique features in the presence of waves. Velocity fluctuations are found to be strongly dependent on the wave phase; variations of these fluctuations are explained by the pressure–strain correlation associated with the wave-induced strain field.     

含气泡液体中的非线性声传播

考虑到分布在液体中的气泡是声波在含气泡液体中传播时引起非线性的一个很重要的因素,本文研究了声波在含气泡液体中的非线性传播.将气体含量的影响引入到声波在液体中传播的方程中,从而得到声波在气液混合物中传播的数学模型.通过对该模型进行数值模拟发现,气体含量、驱动声场声压幅值及驱动声场作用时间均会影响到气液混合物中的声场分布及声压幅值大小.液体中的气泡会"阻滞"液体中声场的传播并将能量"聚集"在声源附近.对于连续大功率的驱动声场来说,液体中的气泡会"阻滞"气液混合物中声场及其能量的传播.