Generalized Fourier transforms for the acoustic pressure and velocity in compressible dissipative stratified media

Generalized Fourier transforms are derived for the acoustic pressure and velocity in compressible dissipative plane stratified media. The acoustic source terms are accounted for in the equations of continuity and force. The acoustic pressure and velocity are each expressed as sums of two infinite (branch cut) integrals and a discrete term. In the far field the infinite integrals correspond to the direct and specularly reflected waves and the lateral wave. The discrete term associated with the pole of the reflection coefficient is the surface wave. The transforms provide a suitable basis for the expansion of the acoustic pressure and the velocity when the height of the interface, the adiabatic bulk modulus, the equilibrium density, and the absorption of the medium vary. Both exact boundary conditions and the approximate impedance boundary condition are considered in this work.     

Diffraction by a metallic wedge covered with a dielectric material

In the literature, the usual formulation which assumes a constant impedance boundary condition does not permit one to express both the value of the geometrical optics field and the actual excited surface wave, since the impedance associated with the material for the reflection and for the surface wave are generally different. To obtain these elements in the field expression, we replace the notion of constant impedance by a differential operator. The general solution presented here preserves enough degrees of freedom to satisfy the continuity of fields at the internal junction of the two materials which cover each face of the wedge. The numerical calculations are based upon a function, related to the Maliuzhinets function, in a form which is easily computable.     

Acoustic radiation force on a rigid cylinder in an off-axis Gaussian beam near an impedance boundary

The exact equations of the axial and transverse acoustic radiation force functions of a Gaussian beam arbitrarily incident on an infinite rigid cylinder close to an impedance boundary and immersed in an ideal fluid are deduced by expressing the incident wave, the scattering wave and the boundary reflected wave in terms of the cylindrical wave function. The effects of the beam waist, the sound reflection coefficient, the cylinder position and the distance from the impedance boundary on the acoustic radiation force are studied using numerical simulations. The simulation results show that the amplitude of the acoustic radiation force function increases with beam width. Moreover, the values of the acoustic radiation force in both the axial and transverse directions reach those of a plane wave when the beam width is considerably larger than the wavelength of the Gaussian beam. The properties of the impedance boundary and the position of the cylinder in the Gaussian beam have a considerable effect on the magnitude and direction of the force. The simulation results, particularly in the case of a transverse force, indicate the presence of a negative acoustic radiation force that is related to the nondimensional frequency and position of the cylinder in the Gaussian beam.     

Modal acoustic impedance force on a spherical source near a rigid interface

The modal acoustic radiation load on a spherical surface undergoing angularly periodic axisymmetric harmonic vibrations while immersed in an acoustic halfspace with a rigid (infinite impedance) planar boundary is analyzed in an exact fashion using the classical technique of separation of variables. The formulation utilizes the appropriate wave field expansions, the classical method of images and the appropriate translational addition theorem to simulate the relevant boundary conditions for the given configuration. The associated acoustic field quantities such as the modal impedance matrix and the modal acoustic radiation force acting on the spherical surface are determined. The analytical results are illustrated with a numerical example in which the spherical surface, excited in vibrational modes of various orders, is immersed near an impervious rigid wall. The presented solution could eventually be used to validate those obtained by numerical approximation techniques.     

广义有限差分法在含阻抗边界空腔声学分析中的应用

声学分析在噪声控制、室内隔音等工程计算中有着重要的作用.由于现实生活中的声学模型往往伴随着吸声材料,因此分析含阻抗边界条件的声学问题显得十分必要.广义有限差分法是一种新型区域型无网格数值离散方法,该方法基于多元函数泰勒级数展开式和加权最小二乘拟合,将控制方程中未知参量的各阶偏导数表示为相邻节点函数值的线性组合.本文首次...     

Artificial boundary conditions for out-of-plane motion in penta-graphene

In this paper, artificial boundary conditions are designed for out-of-plane waves in penta-graphene, a newly proposed allotrope of carbon. By matching the dispersion relation for acoustic branch phonons in the long wave limit, we determine parameters in proposed linear constraints among displacements and velocities at the boundary and nearby atoms. Reflection analysis for normal incidences and a numerical test demonstrate the effectiveness of the artificial boundary conditions. These conditions may be used for studying mechanical behaviours of the novel complex lattice of penta-graphene.     

Inverse scattering problem from an impedance crack via a composite method

In this paper we consider an inverse scattering problem from an open arc with impedance boundary conditions on both sides of the crack. Our aim is to recover both the impedance function and the unknown crack simultaneously from the far-field pattern with only one incident wave. Making the most out of the direct problem, a straightforward method of iterative nature is developed for the inverse problem. The ill-posedness of this problem is considered by incorporating the Tikhonov regularization. Numerical examples are provided at the end of the paper to show the feasibility of our method.     

Wave propagation in a fluid wedge over a solid half-space – Mesh-free analysis with experimental verification

Interaction between a bounded ultrasonic beam and a liquid wedge over a solid half-space is studied. A semi-analytical technique called distributed point source method (DPSM) is adopted for modeling the ultrasonic field in a wedge shaped fluid structure over a solid half space. This study is important for analyzing the ultrasonic waves used for the non-destructive inspections of partially immersed structures. It is also useful for studying the effect of underwater ultrasonic or acoustic wave experiments on marine lives near the shore. The problem geometry considers a bounded acoustic beam striking a fluid–solid interface between a fluid wedge and a solid half-space at steady-state. Solution of this problem is beyond the scope of the currently available analytical methods when the beam is bounded. Only numerical method (boundary element method (BEM) or finite element method (FEM)) based packages (e.g. PZFlex) are in principle capable of modeling ultrasonic fields in such structures. At high frequencies FEM and BEM based packages require huge amount of computation memory and time for their executions that DPSM technique can avoid. Effect of the angle of strike and the fluid wedge angle variation on the wave propagation characteristics is studied. Theoretical predictions are compared with some experimental results.     

Axial acoustic radiation force on a rigid cylinder near an impedance boundary for on-axis Gaussian beam

This work presents a theoretical model to calculate the acoustic radiation force on a rigid cylindrical particle immersed in an ideal fluid near a boundary for an on-axis Gaussian beam. An exact solution of the axial acoustic radiation force function is derived for a cylindrical particle by applying the translation addition theorem of cylindrical Bessel function. We analyzed the effects of the impedance boundary on acoustic radiation force of a rigid cylinder immersed in water near an impedance boundary with particular emphasis on the radius of the rigid cylinder and the distance from the cylinder center to impedance boundary. Simulation results reveal that the existence of particle trapping behavior depends on the choice of nondimensional frequency as well as the offset distance from the impedance boundary. The value of the radiation force function varies when the cylinder lies at the different position of the on-axis Gaussian beam. For the particle with different radius, the acoustic radiation force functions vary significantly with frequency. This study provides a theoretical basis for acoustic manipulation, which may benefit to the improvement and development of the acoustic control technology.     

Acoustic radiation force on an elastic cylinder in a Gaussian beam near an impedance boundary

Acoustic radiation force (ARF) is studied by considering an infinite elastic cylinder near an impedance boundary when the cylinder is illuminated by a Gaussian beam. The surrounding fluid is an ideal fluid. Using the method of images and the translation-addition theorem for the cylindrical Bessel function, the resulting sound field including the incident wave, its reflection from the boundary, the scattered wave from the elastic cylinder, and its image are expressed in terms of the cylindrical wave function. Then, we deduce the exact equations of the axial and transverse ARFs. The solutions depend on the cylinder position, cylinder material, beam waist, reflection coefficient, distance from the impedance boundary, and absorption in the cylinder. To analyze the effects of the various factors intuitively, we simulate the radiation force for non-absorbing elastic cylinders made of stainless steel, gold, and beryllium as well as for an absorbing elastic cylinder made of polyethylene, which is a well-known biomedical polymer. The results show that the impedance boundary, cylinder material, absorption in the cylinder, and cylinder position in the Gaussian beam significantly affect the magnitude and direction of the force. Both stable and unstable equilibrium regions are found. Moreover, a larger beam waist broadens the beam domain, corresponding to non-zero axial and transverse ARFs. More importantly, negative ARFs are produced depending on the choice of the various factors. These results are particularly important for designing acoustic manipulation devices operating with Gaussian beams.     

Simulation of Transverse Combustion Instability in a Multi-Injector Combustor using the Time-Domain Impedance Boundary Conditions

The time-domain impedance boundary condition (TDIBC) is used as a reduced-order model (ROM) in large-eddy simulation (LES) to study self-sustained transverse oscillations in an experimentally studied high-pressure, shear coaxial multi-injector combustor. This work is an extension of the recent study using ROM-LES to simulate a single-element combustor that exhibited longitudinal instability. Here, we focus on transverse instability in a seven-injector combustor. The fuel and oxidizer inlets are truncated and the conventional inflow boundary conditions at the original inlet are replaced by an impedance describing function (IDF) in the form of a reflection coefficient that couples with LES through characteristic based boundary conditions at the truncated inlet. The impedance model is also generalized to include the effects of entropy fluctuations at the inflow. The hybrid ROM-LES simulations are compared with LES simulations with the full combustor geometry. Results show very good agreement and confirm that the use of TDIBC within LES is a viable tool to account for complex acoustic/boundary interaction in a physical way without explicitly solving the full geometry at LES level. Some simplifications and approximations have to be invoked and these constraints are also discussed.     

Nonreflecting boundary conditions based on nonlinear multidimensional characteristics

Nonlinear characteristic boundary conditions based on nonlinear multidimensional characteristics are proposed for 2‐ and 3‐D compressible Navier–Stokes equations with/without scalar transport equations. This approach is consistent with the flow physics and transport properties. Based on the theory of characteristics, which is a rigorous mathematical technique, multidimensional flows can be decomposed into acoustic, entropy, and vorticity waves. Nonreflecting boundary conditions are derived by setting corresponding characteristic variables of incoming waves to zero and by partially damping the source terms of the incoming acoustic waves. In order to obtain the resulting optimal damping coefficient, analysis is performed for problems of pure acoustic plane wave propagation and arbitrary flows. The proposed boundary conditions are tested on two benchmark problems: cylindrical acoustic wave propagation and the wake flow behind a cylinder with strong periodic vortex convected out of the computational domain. This new approach substantially minimizes the spurious wave reflections of pressure, density, temperature, and velocity as well as vorticity from the artificial boundaries, where strong multidimensional flow effects exist. The numerical simulations yield accurate results, confirm the optimal damping coefficient obtained from analysis, and verify that the method substantially improves the 1‐D characteristics‐based nonreflecting boundary conditions for complex multidimensional flows. Copyright © 2009 John Wiley & Sons, Ltd.     

Diffraction of a solitary wave by a thin wedge

The diffraction of a solitary wave by a thin wedge with vertical walls is studied when the incident solitary wave is directed along the wedge axis. The method of multiple scales is extended to this problem and reduces the task to that of solving the two-dimensional KdV equation with proper boundary and initial conditions. The finite-difference numerical procedure is carried out with the fractional step algorithm in which difference schemes are all implicit. Except the maximum run-up at the wall, the results in this paper are found to corroborate the Melville's experiments not only qualitatively but also quantitatively. The maximum run-up of our results agrees well with Funakoshi's numerical one but it is considerably larger than that in Melville's experiment. An important reason for this discrepancy is believed to be the effect of viscous boundary layer on the vertical side wall.     

Effect of boundary layer on Mach reflection over a wedge surface

Abstract. In this paper, we consider the phenomenon of unsteady Mach reflection generated by a plane shock wave advancing over a straight wedge surface, with particular attention to the deviation of the flow field from the self-similar nature. We examine the observed change in angle between incident and reflected shocks, which is in contrast to the fact that the angle should remain constant with time in a self-similar flow. The effect of the boundary layer behind the advancing shock wave over the surface of the wedge is considered to cause this, and boundary layer theory is utilized to estimate the thickness of the layer. It is found that the thickness increases as to the time t compared with t by the overall expansion in the self-similar flow. Assuming that the thicker boundary layer is effectively equivalent to a change in wedge angle, the effect of the boundary layer on the flow field should be less in later stages with larger t values in accordance with the observation above. Received 6 March 2000 / Accepted 23 April 2001     

Acoustic Field Effect on Laminar Turbulent Transition on a Swept Wing in the Favorable Pressure Gradient Region

The development of disturbances in a three-dimensional boundary layer on a swept wing model is studied both under natural conditions and for artificial excitation of traveling waves by an acoustic field. It is found that steady-state streamwise structures are formed in the three-dimensional boundary layer; under natural conditions a wave packet leading to turbulence is detected. When the flow is exposed to the action of an acoustic field at a frequency from the wave packet, disturbances whose velocity along the streamwise structures is equal to 0.55 of the oncoming flow velocity are formed, while the laminar-turbulent transition is displaced upstream.     

Normal incidence of a plane acoustic wave on a rigid wall covered with a thin compressible layer

The one-dimensional unsteady problem of the variation of the pressure on a rigid wall covered with a thin compressible layer upon which a plane acoustic wave impinges is investigated. The investigation is carried out from two standpoints: without allowance for wave processes in the layer (in this case the layer is modeled by means of a special boundary condition [1] and the pressure on the wall is a continuous function of time) and with allowance for the waves transporting the pressure perturbation from the outer edge of the layer to the wall and back (in this case the pressure on the wall is a piecewise-continuous function of time). A criterion of the proximity of the results of the two approaches is the smallness of the acoustic impedance ratio before interaction begins. This holds true even when the high intensities of the incident waves lead to considerable compression of the layer and an increase in its acoustic impedance.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 139–148, July–August, 1988.     

Effects of surface impedance on current density in a piezoelectric resonator for impedance distribution sensing

We study the relationship between the surface mechanical load represented by distributed acoustic impedance and the current density distribution in a shear mode piezoelectric plate acoustic wave resonator. A theoretical analysis based on the theory of piezoelectricity and trigonometric series is performed. In the specific and basic case when the surface load is due to a local mass layer, numerical results show that the current density concentrates under the mass layer and is sensitive to the phy...     

On the stability and dissipation of wall boundary conditions for compressible flows

Characteristic formulations for boundary conditions have demonstrated their effectiveness to handle inlets and outlets, especially to avoid acoustic wave reflections. At walls, however, most authors use simple Dirichlet or Neumann boundary conditions, where the normal velocity (or pressure gradient) is set to zero. This paper demonstrates that there are significant differences between characteristic and Dirichlet methods at a wall and that simulations are more stable when using walls modelled with a characteristic wave decomposition. The derivation of characteristic methods yields an additional boundary term in the continuity equation, which explains their increased stability. This term also allows to handle the two acoustic waves going towards and away from the wall in a consistent manner. Those observations are confirmed by stability matrix analysis and one‐ and two‐dimensional simulations of acoustic modes in cavities. Copyright © 2009 John Wiley & Sons, Ltd.