基于比拟理论的翼型扰流声场数值模拟

从二维模型方程的全离散形式出发,重点分析了差分格式的色散特性和各向异性效应,证实迎风紧致格式比对称格式有更好的色散和各向同性特性,故有利于声场的数值模拟,并采用三阶迎风紧致格式（ＵＣＤ３）和四阶对称紧致格式（ＳＣＤ４）计算了绕ＮＡＣＡ００１２翼型的可压缩非定常流场,并将此流场作为近场声源,运用声学比拟理论对气动声进行模拟。     

考虑地面效应的高速列车远场气动噪声计算方法研究

为研究高速列车远场气动噪声的计算方法,根据高速列车近地面运行的实际情况,利用半自由空间的Green函数求解FW-H方程;建立考虑地面效应时的远场声学积分公式,并研究地面效应对高速列车远场气动噪声的影响.研究表明,由于存在地面效应,原来的自由声场变成了相当于真实列车声场与镜像列车声场的叠加,并且作用在镜像列车上的力源和法向运动速度与真实列车上的相同.当列车运动速度为350 km/h时,不考虑地面效应时,远场测点的等效连续A计权声压级的最大值为90.76 dB;考虑地面效应之后,远场测点的等效连续A计权声压级的最大值为94.72 dB.     

Use of reciprocity considerations for the two-dimensional BEM analysis of wave propagation in an elastic half-space with applications to acoustic emission

A simple numerical treatment of the infinite boundary in the BEM analysis of two-dimensional wave propagation problems in elastic half-spaces is proposed to avoid the spurious reflections of non-decaying Rayleigh waves introduced by the truncation of the boundary. The proposed method exploits the knowledge of the far-field asymptotic behavior of the solution to adequately correct the BEM displacement system matrix for the truncated problem to account for the contribution of the omitted part of the boundary. The reciprocal theorem of elastodynamics is used for a convenient computation of this contribution exclusively in terms of the boundary integrals of the original BEM system. The method is applied to the study of the acoustic emission from nucleating and propagating surface-breaking and buried cracks in a two-dimensional elastic half-space. It is shown to be particularly advantageous since it allows for an accurate calculation of the generated signal even when the observation point is located far from the acoustic emission source.     

Harmonic radiation from a liquid-filled spherical acoustic lens with an internal eccentric spherical source

Radiation of sound from a modally vibrating shell-encapsulated (eccentric) spherical source is analyzed in an exact manner using the classical method of separation of variables. The proposed model is a realistic idealization of a spherical acoustic lens with focal point inside the lens when used as a sound projector. The analytical results are illustrated with a numerical example in which the modal acoustic radiation impedance load on the source and the radiated far-field pressure are evaluated for representative values of the parameters characterizing the system. Numerical results clearly illustrate that in addition to frequency, surface velocity distribution and eccentricity of the source, the dynamic interaction of the encapsulating shell can be of great consequence in sound radiation.     

Synthetic Turbulence Generators for RANS-LES Interfaces in Zonal Simulations of Aerodynamic and Aeroacoustic Problems

The paper presents the detailed formulation and validation results of simple and robust procedures for the generation of synthetic turbulence aimed at providing artificial turbulent content at the RANS-to-LES interface within a zonal Wall Modelled LES of attached and mildly separated wall-bounded flows. There are two versions of the procedure. The aerodynamic version amounts to a minor modification of a synthetic turbulence generator developed by the authors previously, but the acoustically adapted version is new and includes an internal damping layer, where the pressure field is computed by “weighting” of the instantaneous pressure fields from LES and RANS. This is motivated by the need to avoid creating spurious noise as part of the turbulence generation. In terms of pure aerodynamics, the validation includes canonical shear flows (developed channel flow, zero pressure gradient boundary layer, and plane mixing layer), as well as a more complex flow over the wall-mounted hump with non-fixed separation and reattachment, with emphasis on a rapid conversion from modeled to resolved Reynolds stresses. The aeroacoustic applications include the flow past a trailing edge and over a two-element airfoil configuration. In all cases the methodology ensures a very acceptable accuracy for the mean flow, turbulent statistics and, also, the near- and far-field noise.     

Stability/dispersion analysis of the discontinuous Galerkin linearized shallow‐water system

The frequency or dispersion relation for the discontinuous Galerkin mixed formulation of the 1‐D linearized shallow‐water equations is analysed, using several basic DG mixed schemes. The dispersion properties are compared analytically and graphically with those of the mixed continuous Galerkin formulation for piecewise‐linear bases on co‐located grids. Unlike the Galerkin case, the DG scheme does not exhibit spurious stationary pressure modes. However, spurious propagating modes have been identified in all the present discontinuous Galerkin formulations. Numerical solutions of a test problem to simulate fast gravity modes illustrate the theoretical results and confirm the presence of spurious propagating modes in the DG schemes. Copyright © 2005 John Wiley & Sons, Ltd.     

带旋转自由度的精化非协调平面四边形等参元

对于C~1类不协调元文[1]提出了一种精比直接刚度法。本文进一步将其应用到C~0类问题,建立一种带旋转自由度的不协调平面四边形单元,其协调部分是用人Allman插值法建立的单元函数,不协调部分用了四个内部自由度。该单元能保证通过分片试验,保持了单变量有限元列式简单、性能可靠(无多余零能模式及坐标不变性)等长处,同时,还具备多变量有限元(杂交/拟协调元)高精度的优点。算例表明,本文提出的单元收敛、可靠、高精度且高效率。     

PIV速度场坏矢量的本征正交分解处理技术

介绍了一种针对粒子图像测速(PIV)基于本征正交分解(POD)的速度场后处理技术.该技术改变了现在后处理技术将速度场坏矢量识别和修正分开实现的局面,通过迭代方法有效地实现了速度场坏点统一的识别和修复算法.算法利用POD分解的低阶模态信息重构出可以用于坏矢量识别的参考速度场,利用该参考速度场对全流场进行坏点识别并完成修正.通过对一套光滑的PIV速度场数据引入高斯分布的随机误差,测试验证了该POD方法的优越性.在坏矢量识别方面新方法较归一化中值检验有更高的正确性,能识别大面积出现的坏矢量区域.在坏矢量修补的插值算法中,新方法的计算效率又高于传统Gappy POD方法,且计算精度优于常见的矢量场内插数学方法.特别是在数据缺失的大连通区域,该方法对物理流场有很好的预测效果.     

A new formulation of the effective elastic–plastic response of two-phase particulate composite materials

In this paper the double-inclusion model, originally developed to determine effective linear elastic properties of composite materials, is reformulated and extended to predict the effective nonlinear elastic–plastic response of two-phase particulate composites reinforced with spherical particles. The resulting problem of elastic–plastic deformation of a double-inclusion embedded in an infinite reference medium subjected to an incrementally applied far-field strain is solved by the finite element method. The proposed double-inclusion model is evaluated by comparison of the model predictions to the available exact results obtained by the direct approach using representative volume elements containing many particles. It is found that the double-inclusion formulation is capable of providing accurate prediction of the effective elastic–plastic response of two-phase particulate composites at moderate particle volume fractions.     

Vibroacoustic comparisons of composite laminated cylindrical shells according to three shear deformation shell theories

Whether the first-order and Reddy third-order shear deformation shell theories are able to evaluate the vibroacoustic responses of laminated cylindrical shells with normal deformation in the high frequency range or not is examined by comparison with a 3D higher-order shear deformation shell theory. The implicit governing equations of arbitrary angle-ply laminated cylindrical shells are derived from the 3D higher-order and Reddy third-order shell theories, and solved on the basis of the Fourier transform. The Reddy third-order shell theory can be obtained as a special case from the 3D higher-order shell theory. The first-order and Reddy third-order shell theories almost give rise to the same vibrational and acoustic results. These two simple shear deformation shell theories can be used to study far-field acoustic radiation from laminated cylindrical shells from the low to high frequency range, but they show some differences from the 3D higher-order shell theory in high frequency vibration of shells. Nevertheless, the differences of vibrational responses seem not to be distinct. The helical wave spectra of the higher-order radial displacements are nearly separate from those of the low-order radial displacement and play a minor role in far-field acoustic radiation, which makes the two simple shell theories applicable in prediction of acoustic power of the cylindrical shells in the much higher frequency range. Moreover, it also results in the fact that far-field sound is least sensitive in comparison with near-field sound and vibration of shells.     

Acoustic radiation induced by bubble motion in compressible fluid

Based on the theory of compressible fluid, a three-dimension boundary element method is utilized to research the motion of bubble. The far-field noise radiation during the growth and contraction is calculated by the Kirchhoff formula and the Ffowcs Williams-Hawkings （FW-H） formula with a fixed radiation surface being arranged at the near-field of bubble as a new acoustic source. The results show that the amplitude of the sound pressure induced by non-spherical bubble is lower than that of spherical bubble in the contraction phase. The retardance effect is more obvious when the observer is farther away from the bubble. In the anaphase of contraction, the observer with the maximum amplitude of sound pressure moves up with the obvious jet. Larger buoyance parameters will generate lower sound pressure amplitudes in the anaphase, while larger intensive parameters will cause higher sound pressure amplitudes in the whole procedure of bubble motion.     

Three-dimensional modeling of elastic guided waves excited by arbitrary sources in viscoelastic multilayered plates

This paper provides a modal solution for the three-dimensional modeling of Lamb and SH waves excited by sources of arbitrary shape. This solution is applicable to elastic and viscoelastic plates, in the far-field as well as in the near-field regions, under the assumption of transverse isotropy about the thickness direction. The theoretical developments are conducted based on a semi-analytical finite element formulation. This formulation yields a one-dimensional modal problem, fast from a computational point of view, and allows to readily handle heterogeneous materials having depth-varying properties (multilayered, piecewise or continuously varying, functionally graded). The modal solution is shown to be expressed in terms of Hankel functions of multiple order thanks to a proper application of inverse transforms and Cauchy residue calculus. The link between the proposed formulation and a fully analytical approach is discussed. The solution of this paper is then successfully compared to literature results and degenerates to the point source case. Formulas are presented to calculate point source excitabilities from lines sources. These formulas remain valid for non-propagating modes, viscoelastic materials and account for the near-field contribution. Finally, the example of a viscoelastic bilayer waveguide excited by a rectangular source is considered in order to check the theoretical results.     

一种求解瑞利波散射问题的修正边界元方法

边界元法的一大优势是用于求解半空间等无限域问题,然而对于弹性波的传播问题,传统边界元法在采用全平面或全空间格林函数时,在截断边界处仍会产生虚假的反射回波,直接影响到散射场的求解准确性。因此,本文在传统边界元法基础上提出一种修正边界元法,用于计算无限大半平面中的弹性波场问题。该方法以瑞利波形式的远端散射场代替原本因截断而舍去的部分,通过互易定理建立单位瑞利波和全平面格林函数的积分方程,求得修正系数,并代入修正边界元矩阵,计算出瑞利波的散射场。为验证本文所提方法,文中将多个算例的结果与解析解对比,并用该方法计算了不同缺陷的散射场。这些对比结果表明,本文所提修正边界元法可准确求解瑞利波散射场,为基于表面波的缺陷反演问题研究提供了有效的正演途径。     

基于Navier-Stokes／Kirchhoff方程的悬停旋翼跨声速气动声学计算

研究了Kirchhoff积分面是否有盖有底,以及是否计及旋翼网格上的流场值,这两个因素对噪声预测结果的影响.发展了一种基于重叠网格的计算悬停旋翼远场噪声的数值方法.数值计算过程分为流场模拟和声场模拟两部分.悬停旋翼流场的数值模拟是在两个相互重叠的网格上进行的:在高质量的旋翼网格上求解Navier-Stokes方程,用于模拟旋翼附近的粘性流动和近场尾涡的捕捉;在远离粘性区域处布置符合悬停流场物理特征的圆柱形背景网格,控制方程为Euler方程,用于远场尾涡的捕捉.计算得到的流场信息插值到用于声场计算的Kirchhoff积分面上.观测点处的噪声可以认为是由这个完全包含桨叶的Kirchhoff积分面上的面元(声源)发声得到.远场声波的传播由Kirchhoff积分公式描述.计算结果表明:采用有盖有底的Kirchhoff积分面并且同时计及旋翼网格流场值时,计算得到的HSI噪声与实验值吻合最好.     

Coherent structures and wavepackets in subsonic transitional turbulent jets

A large eddy simulation (LES) is performed for two subsonic jets with a Reynolds number of Re = 105, which have different core temperatures, i.e., the cold and hot jet. The far-field overall sound pressure levels (OASPL) and noise spectra are well validated against previous exper-imental results. It is found that the OASPL is raised by heating at shallow angles. The most energetic coherent struc-tures are extracted with specified frequencies using the filter based on the frequency domain variant of the snapshot method of proper orthogonal decomposition (POD). The m = 0, 1 modes have high coherence of near-field pres-sure for both jets, while the coherence of m = 0 modes is enhanced greatly by heating. Based on the coherent struc-tures, spatial wavepackets are educed and the characteristics of growth, saturation and decay are analyzed and compared between the two jets in detail. The results show that heat-ing would enhance the linear growth rate for high frequency components, and nonlinear growth rates for low frequency components in general, which are responsible for higher OASPL in the hot jet. The far-field sound generated by wavepackets is computed using the Kirchhoff extrapolation, which matches well with that of LES at shallow angles. This indicates that the wavepackets associated with coherent structures are dominant sound sources in forced transitional turbulent jets. Additionally, the present POD method is proven to be a robust tool to extract the salient features of the wavepackets in turbulent flows.     

Consistent Behavior of Eulerian Monte Carlo fields at Low Reynolds Numbers

A slightly different version of the Eulerian Monte Carlo method (EMC) (Valiño, Flow Turbul. Combust. 60, 157–172, (1998)) is presented in this paper. The EMC method is an effective stochastic numerical approach to solve the Probability Density Function (PDF) of reacting species in turbulent flows. In contrast with the original formulation, the spurious Wiener term associated with the molecular diffusion is removed, by splitting the micro-mixing into mean gradient and fluctuating contributions. The evolution of the EMC fields representing the PDF in the proposed formulation is then consistent in the laminar limit: the EMC fields follow the same standard convection-diffusion equation, without any stochastic terms.     

A stream function–vorticity formulation‐based immersed boundary method and its applications

A new stream function–vorticity formulation‐based immersed boundary method is presented in this paper. Different from the conventional immersed boundary method, the main feature of the present model is to accurately satisfy both governing equations and boundary conditions through velocity correction and vorticity correction procedures. The velocity correction process is performed implicitly based on the requirement that velocity at the immersed boundary interpolated from the corrected velocity field accurately satisfies the nonslip boundary condition. The vorticity correction is made through the stream function formulation rather than the vorticity transport equation. It is evaluated from the firstorder derivatives of velocity correction. Two simple and efficient ways are presented for approximation of velocity‐correction derivatives. One is based on finite difference approximation, while the other is based on derivative expressions of Dirac delta function and velocity correction. It was found that both ways can work very well. The main advantage of the proposed method lies in its simple concept, easy implementation, and robustness in stability. Numerical experiments for both stationary and moving boundary problems were conducted to validate the capability and efficiency of the present method. Good agreements with available data in the literature were achieved. Copyright © 2011 John Wiley & Sons, Ltd.     

高速列车车体长度对气动噪声影响的数值研究

高速列车具有细长形状, 数值评估气动噪声往往需要巨大的计算量.目前对高速列车气动噪声的数值模拟大多基于对简化短编组列车的评估,而实际列车通常具有较长的8$\sim$16节编组.如何基于现有条件合理评价真实长度列车的气动噪声,是一个急需探讨的问题. 本文应用非线性声学求解器(NLAS)和FW--H声学比拟法的混合算法, 先求解噪声积分面上的声场脉动,再进行远场积分, 引入多噪声面积分技术,通过对三种不同长度(3节、4节、6节)列车模型的气动性能和噪声数值模拟,分析了车体长度对列车气动噪声的影响. 结果表明,同一列车模型的各节车厢具有相似的沿线噪声分布,其噪声曲线在量值上十分接近,只是主峰位置会随着车厢空间位置的不同而相应地发生偏移;不同长度编组列车对应部位之间的远场噪声特性具有较强的关联性,它们的远场噪声具有接近的总声压级和噪声频谱.通过利用短编组计算数据进行分解、平移和叠加,成功重构了4编组和6编组列车远场噪声特性,与直接计算结果相比误差在可接受范围内.由此发展了基于短编组列车噪声的数值结果,重构长编组列车沿线噪声的近似评估方法.      

Application of the 2-D constant strain assumption to FEM elements consisting of an arbitrary number of nodes

The formulation for the constant strain element is revisited to develop multi-noded elements that can be used as transition elements in a finite element mesh. Although the constant strain approach is computationally attractive, spurious force-free displacement modes arise for elements consisting of more than three nodes. These unstable mode shapes are typified by “hourglassing” which develops in quadrilateral elements when constant strain is assumed within the element. The means for stabilizing spurious mode shapes for quadrilateral elements is well documented in the literature, however, in this paper, a general formulation for stabilizing forces is presented for elements having an arbitrary number of nodes and therefore is not restricted to quadrilateral elements. This paper examines the use of meshes consisting of constant strain elements created from polygons having differing numbers of element nodes. The effectiveness of the stabilization procedure is illustrated along with “patch test” examples to assess the consistency of the approximation. The elements are shown to be surprisingly robust, yielding reasonable results even when poorly designed mesh transitions are used.     

The effects of initial perturbation to mixing-layer noise

The far-field noise radiated from mixing layers is determined by the near-field flow dynamics which is sensitive to the initial perturbation of instability introduced physically or numerically. This study focuses on the effects of the phase delay in two initial perturbations, one at the fundamental wave number and the other at its subharmonic both calculated from linear instability analysis, on the sound generation in mixing layers. When different phase delays φ₁ changing from zero to 2π is applied on the fundamental mode, we observe different vortex merging processes (e.g. vortex pairing or tearing). The strong nonlinear interaction in the merging process generates most of the noise from mixing layers. There shows a pattern in a period of 2π for the response of far-field sound to the change of φ₁. Similar effects on the dynamics and acoustics can be achieved by adding different phase delays φ₂ to the subharmonic mode instead, however, the response repeats in a period of only π for φ₂. The effects of the combination of different phase delays to other parameters, including the amplitude and wave number for each perturbations, are also investigated. All the results indicate a critical role of nonlinearity in the sound generation mechanism of mixing layers.