Curle's equation and acoustic scattering by a sphere

Recent papers have initiated interesting comparisons between aeroacoustic theory and the results of acoustic scattering problems. In this paper, we consider some aspects of these comparisons for acoustic scattering by a sphere. We give a derivation of Curle's equation for a specific class of linear acoustic scattering problems, and, in response to previous claims to the contrary, give an explicit confirmation of Curle's equation for plane wave scattering by a stationary rigid sphere of arbitrary size in an inviscid fluid. We construct the complete solution for scattering by a rigid sphere in a viscous fluid, and show that the neglect of viscous terms in Curie's equation yields an incomplete prediction of the far field dipole pressure. We also consider the null field solution of the sphere scattering problem, and give its extension to the vorticity modes associated with viscosity. Finally, we construct a solution for an elastic sphere in a viscous fluid, and show that the rigid sphere/null field solution is recovered from the limit of infinite longitudinal and shear wave speeds in the elastic solid.     

Acoustic scattering by a rigid elliptic cylinder in a slightly viscous medium

A complete solution is obtained for the two-dimensional diffraction of a time-harmonic acoustic plane wave by an impenetrable elliptic cylinder in a viscous fluid. Arbitrary size, ellipticity, and angle of incidence are considered. The linearized equations of viscous flow are used to write down expressions for the dilatation and vorticity in terms of products of radially and angular dependent Mathieu functions. The no-slip condition on the rigid boundary then determines the coefficients. The resulting computations are facilitated by recently developed library routines for complex input parameters. The solution for the circular cylinder serves as a guide and a differently constructed solution for the strip is also given. Typical results in the "resonant" range of dimensionless wave number, displaying the surface vorticity and the far-field scattering pattern are included, with the latter allowing comparison with the inviscid case.     

Modal vibrations of a cylindrical radiator over an impedance plane

The problem of acoustic radiation from an infinite cylinder undergoing harmonic modal surface vibrations near a locally reacting planar boundary is considered. The formulation utilizes the appropriate wave field expansions, the classical method of images, and the translational addition theorem for cylindrical wave functions, along with a simple local surface reaction model involving a complex amplitude wave reflection coefficient applied to simulate the relevant boundary conditions for the given configuration. The analytical results are illustrated with a numerical example in which the cylindrical surface is immersed near a layer of fibrous material set on an impervious rigid wall. The numerical results reveal the important effects of interface local surface reaction and source position on the computed modal impedance component values and the radiated on-axis far-field pressure. The benchmark solution presented can lead to a better understanding of acoustic radiation from near-interface two-dimensional sources, which are commonly encountered problems in outdoor acoustics and noise control engineering. Eventually, it could be used to validate those found by numerical approximation techniques.     

谐振管形状对热声发动机内非线性声场的影响

研究了谐振管一端受活塞声源激励,另一端刚性封闭条件下,管道形状对热声发动机谐振管内部非线性声场的影响。基于流体力学基本方程建立了渐变截面谐振管内一维非线性声场的模型,考虑了黏性耗散及非线性效应的影响。利用伽辽金法数值求解了该模型的速度势方程,分析了谐振管形状、活塞振动速度及激励频率对管内声场的影响。将双曲形、指数形、锥形、正弦形等四种变截面谐振管内的非线性声场与圆柱形直管的情况进行了比较。结果反映了谐振管内声场的压力波动受活塞振动速度及谐振管形状的影响;显示了当活塞振动幅度较大时,谐振管内出现的波形畸变、频率曲线偏移、共振频率滞后等非线性现象;揭示了变截面谐振管在抑制管内的高阶谐波及提高压比等方面的优越性。      

Study of the second harmonics in nonlinear scattering of interacting acoustic waves by a rigid cylinder

The field formed by the second harmonics in the case of the nonlinear scattering of interacting plane acoustic waves by a rigid cylinder is considered. The method of successive approximations is used to obtain the solutions to the inhomogeneous wave equation in the first and second approximations. Asymptotic expressions are derived for the components of the total acoustic pressure of the second harmonics, and the scattering diagrams for these components are presented.     

Influence of resonator shape on nonlinear acoustic field in a thermoacoustic engine

The influence of the resonator shape on nonlinear acoustic field in a thermoacoustic engine is studied.The resonator of themoacoustic engine is boundary driving by a piston at one end,and the other end of it is rigid closed.A one-dimensional wave equation that accounts for gas dynamic nonlinearities and viscous dissipation in the resonator is established based on the governing equations of viscous hydromechanics.The nonlinear wave equation is solved using approximate Galerkin method.The nonlinear acoustic field in four different types of shaped resonators including hyperbolical,exponential,conical and sinusoidal are obtained and compared with that of a cylindrical resonator.It is found that the amplitude and waveform of the pressure are strongly affected by the resonator shape,the driving amplitude and the oscillation frequency of the piston.Waveform distortion,resonance frequency shift and hysteresis are observed,when the piston oscillation amplitude is large enough.The advantages of shaped resonator for thermoacoustic engine lie in inhibition of higher order harmonics and improvement of pressure ratio,etc.     

Theoretical and numerical calculations for the time-averaged acoustic force and torque acting on a rigid cylinder of arbitrary size in a low viscosity fluid

In this paper, theoretical calculations as well as numerical simulations are performed for the time-averaged acoustic force and torque on a rigid cylinder of arbitrary size in a fluid with low viscosity, i.e., the acoustic boundary layer is thin compared to the cylinder radius. An exact analytical solution and its approximation are proposed in the form of an infinite series including Bessel functions. These solutions can be evaluated easily by a mathematical software package such as mathematica and matlab. Three types of incident waves, plane traveling wave, plane standing wave, and dual orthogonal standing waves, are investigated in detail. It is found that for a small particle, the viscous effects for an incident standing wave may be neglected but those for an incident traveling wave are notable. A nonzero viscous torque is experienced by the rigid cylinder when subjected to dual orthogonal standing waves with a phase shift even when the cylinder is located at equilibrium positions without imposed acoustic forces. Furthermore, numerical simulations are carried out based on the FVM algorithm to verify the proposed theoretical formulas. The theoretical results and the numerical ones agree with each other very well in all the cases considered.     

Mean force on a small sphere in a sound field in a viscous fluid

A mean force exerted on a small rigid sphere by a sound wave in a viscous fluid is calculated. The force is expressed as a sum of drag force coming from the external steady flow existing in the absence of the sphere and contributions that are cross products of velocity and velocity derivatives of the incident field. Because of the drag force and an acoustic streaming generated near the sphere, the mean force does not coincide with the acoustic radiation pressure, i.e., the mean momentum flux carried by the sound field through any surface enclosing the sphere. If the sphere radius R is considerably smaller than the viscous wave penetration depth delta, the drag force can give the leading-order contribution (in powers of delta/R) to the mean force and the latter can then be directed against the radiation pressure. In another limit, delta< or =R, the drag force and acoustic streaming play a minor role, and the mean force reduces to the radiation pressure, which can be expressed through source strengths of the scattered sound field. The effect of viscosity can then be significant only if the incident wave is locally plane traveling.     

The use of the virtual source technique in computing scattering from periodic ocean surfaces

In this paper the virtual source technique is used to compute scattering of a plane wave from a periodic ocean surface. The virtual source technique is a method of imposing boundary conditions using virtual sources, with initially unknown complex amplitudes. These amplitudes are then determined by applying the boundary conditions. The fields due to these virtual sources are given by the environment Green's function. In principle, satisfying boundary conditions on an infinite surface requires an infinite number of sources. In this paper, the periodic nature of the surface is employed to populate a single period of the surface with virtual sources and m surface periods are added to obtain scattering from the entire surface. The use of an accelerated sum formula makes it possible to obtain a convergent sum with relatively small number of terms (～40). The accuracy of the technique is verified by comparing its results with those obtained using the integral equation technique.     

Influence of viscosity on the diffraction of sound by a circular aperture in a plane screen

The linearized equations of viscous fluid flow are used to analyze the diffraction of a time-harmonic acoustic plane wave by a circular aperture in a rigid plane screen. Arbitrary aperture size and arbitrary angle of incidence are considered. Sets of dual integral equations are derived for the diffracted velocity and pressure fields, and are solved by analytic reduction to sets of linear algebraic equations. In the case of normal incidence, numerical results are presented for the fluid velocity in the aperture and the power absorption due to viscous dissipation. The theoretical results for power absorption are compared to previously obtained results from high amplitude acoustic experiments in air. The conditions under which the dissipation predicted by linear theory becomes significantare quantified in terms of the fluid viscosity and sound speed, the acoustic frequency, and the aperture radius.     

Acoustic scattering by multiple elliptical cylinders using collocation multipole method

This paper presents the collocation multipole method for the acoustic scattering induced by multiple elliptical cylinders subjected to an incident plane sound wave. To satisfy the Helmholtz equation in the elliptical coordinate system, the scattered acoustic field is formulated in terms of angular and radial Mathieu functions which also satisfy the radiation condition at infinity. The sound-soft or sound-hard boundary condition is satisfied by uniformly collocating points on the boundaries. For the sound-hard or Neumann conditions, the normal derivative of the acoustic pressure is determined by using the appropriate directional derivative without requiring the addition theorem of Mathieu functions. By truncating the multipole expansion, a finite linear algebraic system is derived and the scattered field can then be determined according to the given incident acoustic wave. Once the total field is calculated as the sum of the incident field and the scattered field, the near field acoustic pressure along the scatterers and the far field scattering pattern can be determined. For the acoustic scattering of one elliptical cylinder, the proposed results match well with the analytical solutions. The proposed scattered fields induced by two and three elliptical–cylindrical scatterers are critically compared with those provided by the boundary element method to validate the present method. Finally, the effects of the convexity of an elliptical scatterer, the separation between scatterers and the incident wave number and angle on the acoustic scattering are investigated.     

Nonlinear dynamics of a cavitation bubble pair near a rigid boundary in a standing ultrasonic wave field

The dynamics of a micrometer-sized bubble pair in water near a rigid boundary under standing ultrasonic wave excitation is investigated in this study. The viscous effect in the boundary layer at the air-water interface is considered following the viscous correction model. The evolution of the bubble surface at the collapsing stage of the bubble pair is presented for different parameter sets. The field pressure near the rigid boundary, which is induced by the oscillating bubble pair, and the high-speed water jet at the collapse stage, form the main focus of the analysis. This reveals that a horizontal configuration of the bubble pair retards the strength of the bubble jet towards the boundary, whilst a vertical configuration, especially with differently-sized bubbles, can enhance the bubble collapse. This study may help to understand the interaction of multiple bubbles in an acoustic field and its application to surface cleaning.     

Computational simulation of wave propagation problems in infinite domains

This paper deals with the computational simulation of both scalar wave and vector wave propagation problems in infinite domains. Due to its advantages in simulating complicated geometry and complex material properties, the finite element method is used to simulate the near field of a wave propagation problem involving an infinite domain. To avoid wave reflection and refraction at the common boundary between the near field and the far field of an infinite domain, we have to use some special treatments to this boundary. For a wave radiation problem, a wave absorbing boundary can be applied to the common boundary between the near field and the far field of an infinite domain, while for a wave scattering problem, the dynamic infinite element can be used to propagate the incident wave from the near field to the far field of the infinite domain. For the sake of illustrating how these two different approaches are used to simulate the effect of the far field, a mathematical expression for a wave absorbing boundary of high-order accuracy is derived from a two-dimensional scalar wave radiation problem in an infinite domain, while the detailed mathematical formulation of the dynamic infinite element is derived from a two-dimensional vector wave scattering problem in an infinite domain. Finally, the coupled method of finite elements and dynamic infinite elements is used to investigate the effects of topographical conditions on the free field motion along the surface of a canyon.     

A characterization of the scattered acoustic intensity field in the resonance region for simple spheres

The properties of the scattered acoustic vector fields generated by simple spheres illuminated by monotonic continuous wave (CW) plane waves are investigated. Analytical solutions are derived from general acoustic pressure scattering models and analyzed for wave numbers in the resonance region. Of particular interest is the understanding of the characteristics of the scattered acoustic vector field in the near-to-far-field transition region. The separable active and reactive components of the acoustic intensity are used to investigate the structural features of the scattered field components. Numerical results are presented for the near and transition regions for a rigid sphere. A method of mapping nulls in the scattered intensity field components is described. The analysis is then extended to include a simple fluid-filled boundary and finally the evacuated thin-walled shell. Near field acoustic intensity field structures are compared against mechanical material properties of vacuous shells. The ability to extract scattered field features is illustrated with measurements obtained from a recent in-air experiment using an anechoic chamber and acoustic vector sensor probes to measure the scattered acoustic vector field from rigid spheres.     

On the sound fields of infinitely long strips

Exact solutions are derived for sound radiation from four kinds of infinitely-long strips: namely a rigid strip in a baffle of finite width, a resilient strip in free space, and a resilient or rigid strip in an infinite baffle. In one limit, the strip in a finite baffle becomes a rigid strip in free space and in the other, a line source in a finite baffle. Here "rigid" means that the surface velocity is uniform, whereas "resilient" means that the surface pressure is uniform, and the strip is assumed to have zero mass or stiffness, as if a force were driving the acoustic medium directly. According to the Babinet-Bouwkamp principle, radiation from a resilient strip in an infinite baffle is equivalent to diffraction of a plane wave through a slit in the same. Plots are shown for the radiation impedances, far-field directivity patterns, and on-axis pressure responses of the four kinds of strip. A simple relationship between the radiation admittance of the rigid strip in an infinite baffle and the resilient strip in free space is presented. The two-dimensional rectangular wave functions developed in this paper can be applied to related problems.     

Factors influencing electromagnetic scattering from the dielectric periodic surface

The scattering characteristics of the periodic surface of infinite and finite media are investigated in detail.The Fourier expression of the scattering field of the periodic surface is obtained in terms of Huygens' s principle and Floquet's theorem.Using the extended boundary condition method(EBCM) and T-matrix method, the scattering amplitude factor is solved,and the correctness of the algorithm is verified by use of the law of conservation of energy.The scattering cross section of the periodic surface in the infinitely long region is derived by improving the scattering cross section of the finite period surface.Furthermore, the effects of the incident wave parameters and the geometric structure parameters on the scattering of the periodic surface are analyzed and discussed.By reasonable approximation, the scattering calculation methods of infinite and finite long surfaces are unified.Besides, numerical results show that the dielectric constant of the periodic dielectric surface has a significant effect on the scattering rate and transmittance.The period and amplitude of the surface determine the number of scattering intensity peaks, and, together with the incident angle, influence the scattering intensity distribution.     

Ultrasonic scattering cross sections of shell-encapsulated gas bubbles immersed in a viscoelastic liquid: first and second harmonics

The oscillations of gas bubbles, without shell, immersed in viscoelastic liquids and driven by an acoustic wave have been the subject of several investigations. They demonstrate that the viscosity coefficient and the spring constant of the liquid have significant influence on the scattering cross section of the gas bubble. For shell-encapsulated gas bubbles, the investigations have been concentrated to bubbles immersed in a pure viscous liquid. This present work computes the ultrasonic scattering cross section, first and second harmonics, of shell-encapsulated gas bubbles immersed in a viscoelastic liquid. The theoretical model of the bubble oscillation is based on the generalized Rayleigh-Plesset equation of motion of a spherical cavity immersed in a viscoelastic liquid represented by a three-parameter linear Oldroyd model. The scattering cross section is computed for Albunex type of bubble (shell thickness=15 nm, shell shear viscosity=1.77 Pas, shell modulus of rigidity=88.8 MPa) irradiated by a 3.5 MHz ultrasonic pressure wave with an amplitude of 30 kPa. The results demonstrate that encapsulated bubbles respond independently of the surrounding liquid being pure viscous or viscoelastic as long as the surrounding liquid shear viscosity is as low as 10(-3) Pas. Nevertheless, for higher shear viscosities, the bubble responds differently if the surrounding liquid is pure viscous or viscoelastic. In general, the scattering cross sections of first and second harmonics are larger for the viscoelastic liquid.     

Gauss声束对离轴椭圆柱的声辐射力矩

利用部分波展开法求解得到了Gauss声束入射下刚性和非刚性椭圆柱的声散射系数,推导了一般情况下的声辐射力矩表达式.在此基础上,通过一系列数值仿真详细分析了离轴距离、入射角度和束腰半径对声辐射力矩的影响.结果表明:正向与负向声辐射力矩均可以在一定条件下存在;低频情况下刚性椭圆柱比非刚性椭圆柱更容易产生较强的声辐射力矩;特定频率的入射声场可以激发出非刚性椭圆柱不同阶的共振散射模式,因而非刚性椭圆柱的声辐射力矩峰值与频率的关系更密切;增加束腰半径有利于扩大散射截面,进而增加椭圆柱的声辐射力矩.该研究结果预期可以为利用声辐射力矩实现粒子的可控旋转和流体黏度的反演提供一定的理论指导.     

Errors when assuming locally reacting boundary condition in the estimation of the surface acoustic impedance

Numerical simulations usually require boundary conditions in terms of surface acoustic impedance. The surface acoustic impedance depends on the porous material acoustic properties (e.g., characteristic impedance and wave number) and its thickness as well as the type of wave front impinging on its surface. The locally reactive behaviour hypothesis is often assumed to simplify the choice of proper boundary conditions assigning a constant acoustic impedance value on the porous material surface at a given frequency and for each angle of sound incidence. This hypothesis is also used in measurement procedures or for the estimation of the edge effects.     

Effective dynamic properties of 3D composite materials containing rigid penny-shaped inclusions

The propagation of time-harmonic plane elastic waves in infinite elastic composite materials consisting of linear elastic matrix and rigid penny-shaped inclusions is investigated in this paper. The inclusions are allowed to translate and rotate in the matrix. First, the three-dimensional (3D) wave scattering problem by a single inclusion is reduced to a system of boundary integral equations for the stress jumps across the inclusion surfaces. A boundary element method (BEM) is developed for solving the boundary integral equations numerically. Far-field scattering amplitudes and complex wavenumbers are computed by using the stress jumps. Then the solution of the single scattering problem is applied to estimate the effective dynamic parameters of the composite materials containing randomly distributed inclusions of dilute concentration. Numerical results for the attenuation coefficient and the effective velocity of longitudinal and transverse waves in infinite elastic composites containing parallel and randomly oriented rigid penny-shaped inclusions of equal size and equal mass are presented and discussed. The effects of the wave frequency, the inclusion mass, the inclusion density, and the inclusion orientation or the direction of the wave incidence on the attenuation coefficient and the effective wave velocities are analysed. The results presented in this paper are compared with the available analytical results in the low-frequency range.