Directivity calculation with experimental verification for a conformal array of underwater acoustic transducers

The radiation directivity of a complicated conformal array of underwater acoustic transducers is presented based on the boundary ＂element method. It includes the element directivity of each transducer, the natural beam pattern and the controlled beam pattern of the transducer array. At first, the boundary element model of the conformal array is built up, and then the boundary condition is exerted on the model according to the design and environment in which the transducer array is used, and the radiation directivity of the conformal array is calculated. An experiment has been done to measure the directivity in an anechoic water tank. The calculated and the experimental results are compared and analyzed. They are consistent to each other. It is shown that the boundary element method together with the detailed calculations is successful to simulate and predict the radiation directivity of an underwater acoustic transducer array.     

Response of line focus acoustic microscope to specimen containing a subsurface crack

A numerical model is presented to calculate V(z) and V(x, z) curves for a line focus acoustic microscope and a specimen containing a subsurface crack. In this model, a Gaussian beam which is tracked through the lens into the coupling fluid and into the specimen, interacts with the crack. The numerical approach is based on the solution of singular integral equations by the boundary element method. The system of singular integral equations follows from the conditions at the interface of the coupling fluid and the specimen and on the faces of the crack. An electromechanical reciprocity relation is used to express the voltage at the terminals of the microscope's transducer in terms of the calculated incident and back-scattered fields. V(z) and V(x, z) curves are presented for various locations and orientations of the crack. The characteristic features of the V(z) and particularly the V(x, z) curves, as they relate to crack configuration, are discussed in some detail.     

Implementation of a constrained Ritz series modeling technique for acoustic cavity-structural systems

A prior study [Ginsberg, J. H. (2010b). J. Acoust. Soc. Am. 127, 2749-2758] used Ritz series in conjunction with Hamilton's principle to derive general equations describing the time domain response of an acoustic cavity bounded by an elastic structure. The equations of motion are supplemented by constraint equations that explicitly enforce velocity continuity at the cavity's surface. These constraints are imposed by the surface traction, which is represented by unknown coefficients of Ritz-type series. The resulting set of equations are differential-algebraic type. Three methods are presented to convert the governing equations to forms that are familiar to structural acoustics, including one that transforms them from differential-algebraic type to the standard ordinary differential equations associated with linear multi-degree-of-freedom vibratory systems. In cases where only the structure is excited, the formulation offers options as to how displacement/velocity boundary conditions on the nonstructural boundary are enforced, as well as whether zero pressure boundary conditions are enforced at all. An example of a one-dimensional waveguide that is closed at one end by an oscillator is used to explore the quality of solutions obtained from each of these options. Results for natural frequencies and mode functions are examined for accuracy and convergence.     

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 signature of a submarine hull under harmonic excitation

The structural and acoustic responses of a submarine under harmonic force excitation are presented. The submarine hull is modelled as a cylindrical shell with internal bulkheads and ring stiffeners. The cylindrical shell is closed by truncated conical shells, which in turn are closed at each end using circular plates. The entire structure is submerged in a heavy fluid medium. The structural responses of the submerged vessel are calculated by solving the cylindrical shell equations of motion using a wave approach and the conical shell equations with a power series solution. The far-field radiated sound pressure is then calculated by means of the Helmholtz integral. The contribution of the conical end closures on the radiated sound pressure for the lowest circumferential mode numbers is clearly observed. Results from the analytical model are compared with computational results from a fully coupled finite element/boundary element model.     

Influence of viscosity on the diffraction of sound by a periodic array of screens

The paper contains an analysis of the transmission of a pressure wave through a periodic grating including the influence of the air viscosity. The system of equations in this case consists of the compressible Navier-Stokes equations associated with no-slip boundary conditions on solid surfaces. The problem is reduced to two hypersingular integral equations for determining the velocity components along the slits. These equations are solved by using Galerkin's method with some special trial functions. The results can be applied in designing protective screens for miniature microphones realized in the technology of micro-electro-mechanical systems (MEMS). In this case, the physical dimensions of the device are on the order of the viscous boundary layer so that the viscosity cannot be neglected. The microfluidic model of the screen consists of a periodic array of slits in a substrate. The analysis indicates that the openings in the screen should be on the order of 10 microm in order to avoid excessive attenuation of the signal.     

Optoacoustic imaging with synthetic aperture focusing and coherence weighting

Optoacoustic imaging takes advantage of high optical contrast and low acoustic scattering and has found several biomedical applications. In the common backward mode a laser beam illuminates the image object, and an acoustic transducer located on the same side as the laser beam detects the optoacoustic signal produced by thermoelastic effects. A cross-sectional image is formed by laterally scanning the laser beam and the transducer. Although the laser beam width is generally narrow to provide good lateral resolution, strong optical scattering in tissue broadens the optical illumination pattern and thus degrades the lateral resolution. To solve this problem, a combination of the synthetic aperture focusing technique with coherence weighting is proposed. This method synthesizes a large aperture by summing properly delayed signals received at different positions. The focusing quality is further improved by using the signal coherence as an image quality index. A phantom comprising hair threads in a 1% milk solution was imaged with an optoacoustic imaging system. The results show that the proposed technique improved lateral resolution by 400-800% and the signal-to-noise ratio by 7-23 dB over conventional techniques.     

Acoustic diffraction by a half-plane in a viscous fluid medium

We consider the diffraction of a time-harmonic acoustic plane wave by a rigid half-plane in a viscous fluid medium. The linearized equations of viscous fluid flow and the no-slip condition on the half-plane are used to derive a pair of disjoint Wiener-Hopf equations for the fluid stresses and velocities. The Wiener-Hopf equations are solved in conjunction with a requirement that the stresses are integrable near the edge of the half-plane. Specific wave components of the scattered velocity field are given analytically. A Padé approximation to the Wiener-Hopf kernel function is used to derive numerical results that show the effect of viscosity on the velocity field in the immediate vicinity of the edge of the half-plane.     

The Green's matrix and the boundary integral equations for analysis of time-harmonic dynamics of elastic helical springs

Helical springs serve as vibration isolators in virtually any suspension system. Various exact and approximate methods may be employed to determine the eigenfrequencies of vibrations of these structural elements and their dynamic transfer functions. The method of boundary integral equations is a meaningful alternative to obtain exact solutions of problems of the time-harmonic dynamics of elastic springs in the framework of Bernoulli-Euler beam theory. In this paper, the derivations of the Green's matrix, of the Somigliana's identities, and of the boundary integral equations are presented. The vibrational power transmission in an infinitely long spring is analyzed by means of the Green's matrix. The eigenfrequencies and the dynamic transfer functions are found by solving the boundary integral equations. In the course of analysis, the essential features and advantages of the method of boundary integral equations are highlighted. The reported analytical results may be used to study the time-harmonic motion in any wave guide governed by a system of linear differential equations in a single spatial coordinate along its axis.     

Implementation of nonreflecting boundary conditions for the nonlinear Euler equations

Computationally efficient nonreflecting boundary conditions are derived for the Euler equations with acoustic, entropic and vortical inflow disturbances. The formulation linearizes the Euler equations near the inlet/outlet boundaries and expands the solution in terms of Fourier–Bessel modes. This leads to an ‘exact’ nonreflecting boundary condition, local in space but nonlocal in time, for each Fourier–Bessel mode of the perturbation pressure. The perturbation velocity and density are then calculated using acoustic, entropic and vortical mode splitting. Extension of the boundary conditions to nonuniform swirling flows is presented for the narrow annulus limit which is relevant to many aeroacoustic problems. The boundary conditions are implemented for the nonlinear Euler equations which are solved in space using the finite volume approximation and integrated in time using a MacCormack scheme. Two test problems are carried out: propagation of acoustic waves in an annular duct and the scattering of a vortical wave by a cascade. Comparison between the present exact conditions and commonly used approximate local boundary conditions is made. Results show that, unlike the local boundary conditions whose accuracy depends on the group velocity of the scattered waves, the present conditions give accurate solutions for a range of problems that have a wide array of group velocities. Results also show that this approach leads to a significant savings in computational time and memory by obviating the need to store the pressure field and calculate the nonlocal convolution integral at each point in the inlet and exit boundaries.     

On solution of the vector 3-D locally irregular waveguide problem

The vector 3-D problem of a point-source field in a plane waveguide with a large-scale local inhomogeneity on one of its walls is considered. The field components on the boundary surfaces comply with the Leontovich conditions, which are used as a basis for obtaining expressions for the derivatives of the field vectors normal to the boundaries; these expressions reflect the 3-D nature of the inhomogeneity. The problem is reduced to a system of 2-D integral equations allowing for overexcitation and depolarization of the field scattered by the irregularity. The system of 2-D integral equations is asymptotically transformed over the inhomogeneity region on the surface of the walls bounding the waveguide space into a system of linear integral equations, for which the integration contour is represented by the line between the source and observation point, as well as by the linear geometric contour of the irregularity.State University, St. Petersburg. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 38, No. 8, pp. 785–803, August, 1995.     

A NUMERICAL METHOD FOR SCATTERING FROM ACOUSTICALLY SOFT AND HARD THIN BODIES IN TWO DIMENSIONS

This paper presents a numerical method for predicting the acoustic scattering from two-dimensional (2-D) thin bodies. Both the Dirichlet and Neumann problems are considered. Applying the thin-body formulation leads to the boundary integral equations involving weakly singular and hypersingular kernels. Completely regularizing these kinds of singular kernels is thus the main concern of this paper. The basic subtraction-addition technique is adopted. The purpose of incorporating a parametric representation of the boundary surface with the integral equations is two-fold. The first is to facilitate the numerical implementation for arbitrarily shaped bodies. The second one is to facilitate the expansion of the unknown function into a series of Chebyshev polynomials. Some of the resultant integrals are evaluated by using the Gauss-Chebyshev integration rules after moving the series coefficients to the outside of the integral sign; others are evaluated exactly, including the modified hypersingular integral. The numerical implementation basically includes only two parts, one for evaluating the ordinary integrals and the other for solving a system of algebraic equations. Thus, the current method is highly efficient and accurate because these two solution procedures are easy and straightforward. Numerical calculations consist of the acoustic scattering by flat and curved plates. Comparisons with analytical solutions for flat plates are made.     

声测井中的声压-速度有限差分方法

为了解决薄互储层的声测井问题,提出了声压-速度有限差分方法:用声压和速度矢量做为场变量,分别描述井内流体和井外弹性固体或双相介质。这样选择场变量的优点是:处理脉冲点源(或线源)与套网格技术相比简单得多;在内边界上的差分格式稳定,精度得到了改进;人为边界上的吸收效果较好。用柱坐标分别给出了井壁上流体与弹性固体、流体与双相介质的声压-速度边界条件,并用守恒积分方法处理了井壁上的边界条件。通过用声压-速度有限差分方法模拟弹性固体和双相介质地层的声场,证明了声压-速度有限差分方法的有效性。     

Receiving sensitivity and transmitting voltage response of a fluid loaded spherical piezoelectric transducer with an elastic coating

A method is presented to determine the response of a spherical acoustic transducer that consists of a fluid-filled piezoelectric sphere with an elastic coating embedded in infinite fluid to electrical and plane-wave acoustic excitations. The exact spherically symmetric, linear, differential, governing equations are used for the interior and exterior fluids, and elastic and piezoelectric materials. Under acoustic excitation and open circuit boundary condition, the equation governing the piezoelectric sphere is homogeneous and the solution is expressed in terms of Bessel functions. Under electrical excitation, the equation governing the piezoelectric sphere is inhomogeneous and the complementary solution is expressed in terms of Bessel functions and the particular integral is expressed in terms of a power series. Numerical results are presented to illustrate the effect of dimensions of the piezoelectric sphere, fluid loading, elastic coating and internal material losses on the open-circuit receiving sensitivity and transmitting voltage response of the transducer.     

进一步改进边界元方法以克服振动声辐射计算中解的非唯一性

章利用基于三次B样条插值的边界元方法，对振动体外部声辐射问题进行了研究，对CHIEF法及其改进方法作了进一步的改进，提出在加权余量意义下，通过把内部Helmholtz积分方程与其对内点坐标取导后的方程式作线性叠加，在域外构作的一个小体积块上进行积分以形成补充方程，经与表面Helmholtz积分方程相结合，来求解任意频率下的声辐射问题，并以脉动球和摆动球作为算例，说明本提出的方法能够有效地克服在特殊频率处解的非唯一性问题。     

Prediction capabilities of classical and shear deformable beam models excited by a moving mass

In this paper, a comprehensive assessment of design parameters for various beam theories subjected to a moving mass is investigated under different boundary conditions. The design parameters are adopted as the maximum dynamic deflection and bending moment of the beam. To this end, discrete equations of motion for classical Euler-Bernoulli, Timoshenko and higher-order beams under a moving mass are derived based on Hamilton's principle. The reproducing kernel particle method (RKPM) and extended Newmark-β method are utilized for spatial and time discretization of the problem, correspondingly. The design parameter spectra in terms of the beam slenderness, mass weight and velocity of the moving mass are introduced for the mentioned beam theories as well as various boundary conditions. The results indicate the existence of a critical beam slenderness mostly as a function of beam boundary condition, in which, for slenderness lower than this so-called critical one, the application of Euler-Bernoulli or even Timoshenko beam theories would underestimate the real dynamic response of the system. Moreover, there would be a roughly linear relation between the weight of the moving mass and the design parameters for a certain value of the moving mass velocity in most cases of boundary conditions.     

Asymptotic study of three-dimensional flows of a viscous gas near blunt bodies with a permeable surface

Hypersonic flow of a viscous gas near blunt bodies with a permeable surface is investigated for large Reynolds numbers with intense blow-in (suction) of gas from the surface, including the case when the blow-in (suction) velocity vector is oriented in a different direction from the outer normal to the surface of the body, by asymptotic methods with the help of models of three-dimensional viscous shock and three-dimensional laminar boundary layers. Analytical solutions are presented for the velocity and temperature profiles across the layer and the coefficients of friction and heat transfer on the surface, and the manner in which the basic determining parameters of the problem influence the structure of the perturbed flow and the basic behavior of the integral flow characteristics, which are important for practical applications, and the region of existence of the solution of the initial equations are analyzed.Moscow Physicotechnical Institute. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 4, pp. 4–14, April, 1993.     

Nonlinear standing waves in 2-D acoustic resonators

This paper deals with 2-D simulation of finite-amplitude standing waves behavior in rectangular acoustic resonators. Set of three partial differential equations in third approximation formulated in conservative form is derived from fundamental equations of gas dynamics. These equations form a closed set for two components of acoustic velocity vector and density, the equations account for external driving force, gas dynamic nonlinearities and thermoviscous dissipation. Pressure is obtained from solution of the set by means of an analytical formula. The equations are formulated in the Cartesian coordinate system. The model equations set is solved numerically in time domain using a central semi-discrete difference scheme developed for integration of sets of convection-diffusion equations with two or more spatial coordinates. Numerical results show various patterns of acoustic field in resonators driven using vibrating piston with spatial distribution of velocity. Excitation of lateral shock-wave mode is observed when resonant conditions are fulfilled for longitudinal as well as for transversal direction along the resonator cavity.