Search

Ultrasound attenuation in soft tissue follows a power law as a function of the ultrasound frequency, and in medical ultrasound, power law attenuation is often described by fractional calculus models that contain one or more time- or space-fractional derivatives. For certain time-fractional models, exact and approximate time-domain Green's functions are known, but similar expressions are not available for the space-fractional models that describe power law attenuation. To address this deficiency, a numerical approach for calculating time-domain Green's functions for the Chen–Holm space-fractional wave equation and Treeby–Cox space-fractional wave equation is introduced, where challenges associated with the numerical evaluation of a highly oscillatory improper integral are addressed with the Filon integration formula combined with the Pantis method. Numerical results are computed for both of these space-fractional wave equations at different distances in breast and liver with power law exponents of 1.5 and 1.139, respectively. The results show that these two space-fractional wave equations are causal and that away from the origin, the time-domain Green's function for the Treeby–Cox space-fractional wave equation is very similar to the time-domain Green's function for the time-fractional power law wave equation.Ultrasound attenuation in soft tissue follows a power law as a function of the ultrasound frequency, and in medical ultrasound, power law attenuation is often described by fractional calculus models that contain one or more time- or space-fractional derivatives. For certain time-fractional models, exact and approximate time-domain Green's functions are known, but similar expressions are not available for the space-fractional models that describe power law attenuation. To address this deficiency, a numerical approach for calculating time-domain Green's functions for the Chen–Holm space-fractional wave equation and Treeby–Cox space-fractional wave equation is introduced, where challenges associated with the numerical evaluation of a highly oscillatory improper integral are addressed with the Filon integration formula combined with the Pantis method. Numerical results are computed for both of these space-fractional wave equations at different distances in breast and liver with power law exponents of 1.5 and ...

Harbor porpoises (Phocoena phocoena) use narrow band echolocation signals for detecting and locating prey and for spatial orientation. In this study, acoustic impedance values of tissues in the porpoise's head were calculated from computer tomography (CT) scan and the corresponding Hounsfield Units. A two-dimensional finite element model of the acoustic impedance was constructed based on CT scan data to simulate the acoustic propagation through the animal's head. The far field transmission beam pattern in the vertical plane and the waveforms of the receiving points around the forehead were compared with prior measurement results, the simulation results were qualitatively consistent with the measurement results. The role of the main structures in the head such as the air sacs, melon and skull in the acoustic propagation was investigated. The results showed that air sacs and skull are the major components to form the vertical beam. Additionally, both beam patterns and sound pressure of the sound waves through four positions deep inside the melon were demonstrated to show the role of the melon in the biosonar sound propagation processes in the vertical plane.

A finite element method based numerical indentation technique was used to quantify the effect of the material stiffness variation and the subglottal convergence angle of the vocal fold on the vertical stiffness difference of the medial surface. It was found that the vertical stiffness difference increased with the increasing subglottal angle, and it tended to saturate beyond a subglottal angle of about 50°. The material stiffness variation could be as important as the subglottal angle depending on the actual material properties.

The Bayesian framework for parameter inference provides a basis from which subject-specific reduced-order vocal fold models can be generated. Previously, it has been shown that a particle filter technique is capable of producing estimates and associated credibility intervals of time-varying reduced-order vocal fold model parameters. However, the particle filter approach is difficult to implement and has a high computational cost, which can be barriers to clinical adoption. This work presents an alternative estimation strategy based upon Kalman filtering aimed at reducing the computational cost of subject-specific model development. The robustness of this approach to Gaussian and non-Gaussian noise is discussed. The extended Kalman filter (EKF) approach is found to perform very well in comparison with the particle filter technique at dramatically lower computational cost. Based upon the test cases explored, the EKF is comparable in terms of accuracy to the particle filter technique when greater than 6000 particles are employed; if less particles are employed, the EKF actually performs better. For comparable levels of accuracy, the solution time is reduced by 2 orders of magnitude when employing the EKF. By virtue of the approximations used in the EKF, however, the credibility intervals tend to be slightly underpredicted.

In certain cavitation-based ultrasound techniques, the relative importance of thermally vs mechanically induced damage is unclear. As a first step to investigate this matter, a numerical model for bubble dynamics in tissue-like, viscoelastic media is presented in which full thermal effects are included inside and outside the bubble, as well as interdiffusion of vapor and non-condensible gas inside the bubble. Soft tissue is assumed to behave according to a Kelvin-Voigt model in which viscous and elastic contributions are additive. A neo-Hookean formulation, appropriate for finite-strain elasticity, accounts for the large deformations produced by cavitation. Numerical solutions to problems of relevance to therapeutic ultrasound are examined, and linear analysis is used to explain the underlying mechanisms. The dependence between the surrounding medium's elasticity (shear modulus) and the extent to which the effects of heat and mass transfer influence bubble dynamics is quantified. In particular, the oscillation properties are related to the eigenvalues determined from linear theory. Regimes under which a polytropic relation describes the heat transfer to sufficient accuracy are identified, for which the complexity and computational expense associated with solving full partial differential equations can be avoided.

The evolution of reduced-order vocal fold models into clinically useful tools for subject-specific diagnosis and treatment hinges upon successfully and accurately representing an individual patient in the modeling framework. This, in turn, requires inference of model parameters from clinical measurements in order to tune a model to the given individual. Bayesian analysis is a powerful tool for estimating model parameter probabilities based upon a set of observed data. In this work, a Bayesian particle filter sampling technique capable of estimating time-varying model parameters, as occur in complex vocal gestures, is introduced. The technique is compared with time-invariant Bayesian estimation and least squares methods for determining both stationary and non-stationary parameters. The current technique accurately estimates the time-varying unknown model parameter and maintains tight credibility bounds. The credibility bounds are particularly relevant from a clinical perspective, as they provide insight into the confidence a clinician should have in the model predictions.

The attenuation of ultrasound propagating in human tissue follows a power law with respect to frequency that is modeled by several different causal and noncausal fractional partial differential equations. To demonstrate some of the similarities and differences that are observed in three related time-fractional partial differential equations, time-domain Green's functions are calculated numerically for the power law wave equation, the Szabo wave equation, and for the Caputo wave equation. These Green's functions are evaluated for water with a power law exponent of y = 2, breast with a power law exponent of y = 1.5, and liver with a power law exponent of y = 1.139. Simulation results show that the noncausal features of the numerically calculated time-domain response are only evident very close to the source and that these causal and noncausal time-domain Green's functions converge to the same result away from the source. When noncausal time-domain Green's functions are convolved with a short pulse, no evidence of noncausal behavior remains in the time-domain, which suggests that these causal and noncausal time-fractional models are equally effective for these numerical calculations.

Impedance maps (ZMs) have been proposed as a tool for modeling acoustic properties of tissue microstructure. Three-dimensional (3D) ZMs are constructed from a series of adjacent histological tissue slides that have been stained to emphasize acoustic impedance structures. The power spectrum of a 3DZM can be related to the ultrasound backscatter coefficient, which can be further reduced to a form factor. The goal of this study is to demonstrate the ability to estimate form factors using two-dimensional (2D) ZMs instead of 3DZMs, which have reduced computational and financial cost. The proposed method exploits the properties of isotropic media to estimate the correlation coefficient from slices before estimating the 3D volume power spectrum. Simulated sparse collections of objects were used to study the method by comparing the results obtained using 2DZMs to those predicted by theory. 2DZM analysis was conducted on normal rabbit liver histology and compared to 3DZM analysis of the same histology. The mean percent error between effective scatterer diameter estimates from 2DZMs and 3DZMs of rabbit liver histology was

An acoustic tonometer that measures shifts in resonance frequencies associated with intraocular pressure (IOP) could provide an opportunity for a type of tonometer that can be operated at home or worn by patients. However, there is insufficient theoretical background, especially with respect to the uncertainty in operating frequency ranges and the unknown relationships between IOPs and resonance frequencies. The purpose of this paper is to develop a frequency function for application in an acoustic tonometer. A linear wave theory is used to derive an explicit frequency function, consisting of an IOP and seven other physiological parameters. In addition, impulse response experiments are performed to measure the natural frequencies of porcine eyes to validate the provided function. From a real-time detection perspective, explicitly providing a frequency function can be the best way to set up an acoustic tonometer. The theory shows that the resonance oscillation of the eyeball is mainly dominated by liquid inside the eyeball. The experimental validation demonstrates the good prediction of IOPs and resonance frequencies. The proposed explicit frequency function supports further modal analysis not only of the dynamics of eyeballs, but also of the natural frequencies, for further development of the acoustic tonometer.

The relative role of the various structures in the head of the baiji (Lipotes vexillifer) is examined. A finite element approach was applied to numerically simulate the acoustic propagation through a dolphin's head to examine the relative role of the skull, air sacs, and melon in the formation of the biosonar beam in the vertical plane. The beam pattern obtained with the whole head in place is compared with the beam pattern when the air sac is removed and the other structures (skull and melon) are in place, with only the skull removed, and finally with only the melon removed. The beam pattern with the air sacs and skull intact and the melon removed closely resembled the beam pattern for the complete head, suggesting that the melon has a minor role in the formation of the beam. The beam pattern for the other two cases had very little resemblance to the beam pattern for the whole head. The air sacs seem to have a role of directing propagation of the signal toward the front and the skull prevents the sound propagating below the rostrum. The beam patterns along with a correlation analysis showed that the melon had only a slight influence on the shape and direction of the beam. The resultant beam exiting the head of the dolphin is the result of complex reflection processes within the head of the animal.

Absolute pressure measurements of acoustic emissions by single-element, focused passive cavitation detectors would be facilitated by improved wideband receive calibration techniques. Here, calibration methods were developed to characterize the absolute, frequency-dependent receive sensitivity of a spherically focused, single-element transducer using pulse-echo and pitch-catch techniques. Validation of these calibration methods on a focused receiver were made by generating a pulse from a small diameter source at the focus of the transducer and comparing the absolute pressure measured by a calibrated hydrophone to that of the focused transducer using the receive sensitivities determined here.

The problem of nonlinear acoustic to articulatory inversion mapping is investigated in the feature space using two models, the deep belief network (DBN) which is the state-of-the-art, and the general regression neural network (GRNN). The task is to estimate a set of articulatory features for improved speech recognition. Experiments with MOCHA-TIMIT and MNGU0 databases reveal that, for speech inversion, GRNN yields a lower root-mean-square error and a higher correlation than DBN. It is also shown that conjunction of acoustic and GRNN-estimated articulatory features yields state-of-the-art accuracy in broad class phonetic classification and phoneme recognition using less computational power.

An acoustical/mechanical model of normal adult human middle-ear function is described for forward and reverse transmission. The eardrum model included one component bound along the manubrium and another bound by the tympanic cleft. Eardrum components were coupled by a time-delayed impedance. The acoustics of the middle-ear cleft was represented by an acoustical transmission-line model for the tympanic cavity, aditus, antrum, and mastoid air cell system with variable amounts of excess viscothermal loss. Model parameters were fitted to published measurements of energy reflectance (0.25-13 kHz), equivalent input impedance at the eardrum (0.25-11 kHz), temporal-bone pressure in scala vestibuli and scala tympani (0.1-11 kHz), and reverse middle-ear impedance (0.25-8 kHz). Inner-ear fluid motion included cochlear and physiological third-window pathways. The two-component eardrum with time delay helped fit intracochlear pressure responses. A multi-modal representation of the eardrum and high-frequency modeling of the middle-ear cleft helped fit ear-canal responses. Input reactance at the eardrum was small at high frequencies due to multiple modal resonances. The model predicted the middle-ear efficiency between ear canal and cochlea, and the cochlear pressures at threshold.

The measurement of stress in a structure presents considerable interest in many fields of engineering. In this paper, the diagnostic potential of nonlinear elastic guided waves in a prestressed plate is investigated. To do so, an analytical model is formulated accounting for different aspects involved in the phenomenon. The fact that the initial strains can be finite is considered using the Green Lagrange strain tensor, and initial and final configurations are not merged, as it would be assumed in the infinitesimal strain theory. Moreover, an appropriate third-order expression of the strain energy of the hyperelastic body is adopted to account for the material nonlinearities. The model obtained enables to investigate both the linearized case, which gives the variation of phase and group velocity as a function of the initial stress, and the nonlinear case, involving second-harmonic generation as a function of the initial state of stress. The analysis is limited to Rayleigh-Lamb waves propagating in a plate. Three cases of initial prestress are considered, including prestress in the direction of the wave propagation, prestress orthogonal to the direction of wave propagation, and plane isotropic stress.

An alternative to the spectral overlap assessment metric (SOAM), first introduced by Wassink [(2006). J. Acoust. Soc. Am. 119(4), 2334-2350], is introduced. The SOAM quantifies the intra- and inter-language differences between long-short vowel pairs through a comparison of spectral (F1, F2) and temporal properties modeled with best fit ellipses (F1 × F2 space) and ellipsoids (F1 × F2 × duration). However, the SOAM ellipses and ellipsoids rely on a Gaussian distribution of vowel data and a dense dataset, neither of which can be assumed in endangered languages or languages with limited available data. The method presented in this paper, called the Vowel Overlap Assessment with Convex Hulls (VOACH) method, improves upon the earlier metric through the use of best-fit convex shapes. The VOACH method reduces the incorporation of "empty" data into calculations of vowel space. Both methods are applied to Numu (Oregon Northern Paiute), an endangered language of the western United States. Calculations from the VOACH method suggest that Numu is a primary quantity language, a result that is well aligned with impressionistic analyses of spectral and durational data from the language and with observations by field researchers.

A planewave, incident on a panel, produces an acoustic trace wavelength that propagates along the surface of the panel. The trace wavelength excites the panel into vibration, creating structural waves within the panel that propagate. These structural waves can be purposely Bragg scattered, creating replicas of the trace wavenumber. The replicas are shifted in wavenumber precisely by the inverse of the periodic separation distance l. Hence, in principle, it should be possible to resolve the acoustic trace wavelength from one of the shifted replicas of the panel's response. The incident angle can then be ascertained from the replicated trace wavelengths.

International audience; This paper presents two models of sound production in flute-like instruments that allow time-domain simulations. The models are based on different descriptions of the jet flow within the window of the instrument. The jet-drive model depicts the jet by its transverse perturbation that interacts with the labium to produce sound. The discrete-vortex model depicts the jet as two independent shear layers along which vortices are convected and interact with the acoustic field within the window. The limit of validity between both models is usually discussed according to the aspect ratio of the jet W/h, with W the window length and h the flue channel height. The present simulations, compared with experimental data gathered on a recorder, allow to extend the aspect ratio criterion to the notion of dynamic aspect ratio defined as λ/h where λ is the hydrodynamic wavelength that now accounts for geometrical properties, such as W/h, as well as for dynamic properties, such as the Strouhal number. The two models are found to be applicable over neighboring values of geometry and blowing pressure.

The baiji (Lipotes vexillifer) of the Yangtze River possesses a sophisticated biosonar system. In this study, a finite element approach was used to numerically investigate the propagation of acoustic waves through the head of the Yangtze River dolphin, which possesses an inhomogeneous and complex structure. The acoustic intensity distribution predicted from models with and without the melon and/or skull showed that the emitted sound beam was narrow and formed a highly directed acoustic beam, and the skull and melon significantly enhanced the directional characteristics of the emitted sound. Finally, for a short duration impulsive source, the emitted sound pressure distributions were also simulated at different propagation times. The results provide useful information for better understanding the operation of the biosonar system in this rare and perhaps extinct animal.

The biosonar beampatterns found across different bat species are highly diverse in terms of global and local shape properties such as overall beamwidth or the presence, location, and shape of multiple lobes. It may be hypothesized that some of this variability reflects evolutionary adaptation. To investigate this hypothesis, the present work has searched for patterns in the variability across a set of 283 numerical predictions of emission and reception beampatterns from 88 bat species belonging to four major families (Rhinolophidae, Hipposideridae, Phyllostomidae, Vespertilionidae). This was done using a lossy compression of the beampatterns that utilized real spherical harmonics as basis functions. The resulting vector representations showed differences between the families as well as between emission and reception. These differences existed in the means of the power spectra as well as in their distribution. The distributions were characterized in a low dimensional space found through principal component analysis. The distinctiveness of the beampatterns across the groups was corroborated by pairwise classification experiments that yielded correct classification rates between ~85 and ~98%. Beamwidth was a major factor but not the sole distinguishing feature in these classification experiments. These differences could be seen as an indication of adaptive trends at the beampattern level.

Head-related transfer functions (HRTFs) describe the filtering of the incoming sound by the torso, head, and pinna. As a consequence of the propagation path from the source to the ear, each HRTF contains a direction-dependent, broadband time-of-arrival (TOA). TOAs are usually estimated independently for each direction from HRTFs, a method prone to artifacts and limited by the spatial sampling. In this study, a continuous-direction TOA model combined with an outlier-removal algorithm is proposed. The model is based on a simplified geometric representation of the listener, and his/her arbitrary position within the HRTF measurement. The outlier-removal procedure uses the extreme studentized deviation test to remove implausible TOAs. The model was evaluated for numerically calculated HRTFs of sphere, torso, and pinna under various conditions. The accuracy of estimated parameters was within the resolution given by the sampling rate. Applied to acoustically measured HRTFs of 172 listeners, the estimated parameters were consistent with realistic listener geometry. The outlier removal further improved the goodness-of-fit, particularly for some problematic fits. The comparison with a simpler model that fixed the listener position to the center of the measurement geometry showed a clear advantage of listener position as an additional free model parameter.

An exact series method is presented to analyze classical Langevin transducers with arbitrary boundary conditions. The transducers consist of an axially polarized piezoelectric solid cylinder sandwiched between two elastic solid cylinders. All three cylinders are of the same diameter. The length to diameter ratio is arbitrary. Complex piezoelectric and elastic coefficients are used to model internal losses. Solutions to the exact linearized governing equations for each cylinder include four series. Each term in each series is an exact solution to the governing equations. Bessel and trigonometric functions that form complete and orthogonal sets in the radial and axial directions, respectively, are used in the series. Asymmetric transducers and boundary conditions are modeled by using axially symmetric and anti-symmetric sets of functions. All interface and boundary conditions are satisfied in a weighted-average sense. The computed input electrical admittance, displacement, and stress in transducers are presented in tables and figures, and are in very good agreement with those obtained using atila-a finite element package for the analysis of sonar transducers. For all the transducers considered in the analysis, the maximum difference between the first three resonance frequencies calculated using the present method and atila is less than 0.03%.

Although vocal folds are known to be anisotropic, the influence of material anisotropy on vocal fold vibration remains largely unknown. Using a linear stability analysis, phonation onset characteristics were investigated in a three-dimensional anisotropic vocal fold model. The results showed that isotropic models had a tendency to vibrate in a swing-like motion, with vibration primarily along the superior-inferior direction. Anterior-posterior (AP) out-of-phase motion was also observed and large vocal fold vibration was confined to the middle third region along the AP length. In contrast, increasing anisotropy or increasing AP-transverse stiffness ratio suppressed this swing-like motion and allowed the vocal fold to vibrate in a more wave-like motion with strong medial-lateral motion over the entire medial surface. Increasing anisotropy also suppressed the AP out-of-phase motion, allowing the vocal fold to vibrate in phase along the entire AP length. Results also showed that such improvement in vibration pattern was the most effective with large anisotropy in the cover layer alone. These numerical predictions were consistent with previous experimental observations using self-oscillating physical models. It was further hypothesized that these differences may facilitate complete glottal closure in finite-amplitude vibration of anisotropic models as observed in recent experiments.

A direct numerical simulation of flow-structure interaction is carried out in a subject-specific larynx model to study human phonation under physiological conditions. The simulation results compare well to the established human data. The resulting glottal flow and waveform are found to be within the normal physiological ranges. The effects of realistic geometry on the vocal fold dynamics and the glottal flow are extensively examined. It is found that the asymmetric anterior-posterior laryngeal configuration produces strong anterior-posterior asymmetries in both vocal fold vibration and glottal flow which has not been captured in the simplified models. It needs to be pointed out that the observations from the current numerical simulation are only valid for the flow conditions investigated. The limitations of the study are also discussed.

A linear three-dimensional (3D) elasto-acoustic finite element model was used to simulate the occlusion effect following mechanical vibration at the mastoid process. The ear canal and the surrounding soft and bony tissues were reconstructed using images of a female cadaver head (Visible Human Project(®)). The geometrical model was coupled to a 3D earplug model and imported into comsol Multiphysics (COMSOL(®), Sweden). The software was used to solve for the sound pressure at the eardrum. Finite element modeling of the human external ear and of the occlusion effect has several qualities that can complement existing measuring and modeling techniques. First, geometrically complex structures such as the external ear can be reconstructed. Second, various material behavioral laws and complex loading can be accounted for. Last, 3D analyses of external ear substructures are possible allowing for the computation of a broad range of acoustic indicators. The model simulates consistent occlusion effects (e.g., insertion depth variability). Comparison with an experimental dataset, kindly provided by Stenfelt and Reinfeldt [Int. J. Audiol. 46, 595-608 (2007)], further demonstrates the model's accuracy. Power balances were used to analyze occlusion effect differences obtained for a silicone earplug and to examine the increase in sound energy when the ear canal is occluded (e.g., high-pass filter removal).

A data-based matched-mode source localization method is proposed in this paper for a moving source, using mode wavenumbers and depth functions estimated directly from the data, without requiring any environmental acoustic information and assuming any propagation model. The method is in theory free of the environmental mismatch problem because the mode replicas are estimated from the same data used to localize the source. Besides the estimation error due to the approximations made in deriving the data-based algorithms, the method has some inherent drawbacks: (1) It uses a smaller number of modes than theoretically possible because some modes are not resolved in the measurements, and (2) the depth search is limited to the depth covered by the receivers. Using simulated data, it is found that the performance degradation due to the afore-mentioned approximation/limitation is marginal compared with the original matched-mode source localization method. The proposed method has a potential to estimate the source range and depth for real data and be free of the environmental mismatch problem, noting that certain aspects of the (estimation) algorithms have previously been tested against data. The key issues are discussed in this paper.

This article focuses on computations of low frequency sound propagation from an off-shore wind farm. Two different methods for sound propagation calculations are combined with meteorological data for every 3 hours in the year 2010 to examine the varying noise levels at a reception point at 13 km distance. It is shown that sound propagation conditions play a vital role in the noise impact from the off-shore wind farm and ordinary assessment methods can become inaccurate at longer propagation distances over water. Therefore, this paper suggests that methodologies to calculate noise immission with realistic sound speed profiles need to be combined with meteorological data over extended time periods to evaluate the impact of low frequency noise from modern off-shore wind farms.

An ultrasonic pulse propagating in cancellous bone can be separated into two waves depending on the condition of the specimen. These two waves, which are called the fast wave and the slow wave, provide important information for the diagnosis of osteoporosis. The present study proposes to utilize a signal processing method that extracts the instantaneous frequency (IF) of waveforms from multiple spectral channels. The instantaneous frequency was expected to be able to show detailed time-frequency properties of ultrasonic waves being transmitted through cancellous bone. The employed method, termed the multichannel instantaneous frequency (MCIF) method, showed robustness against background noise as compared to the IF that was directly derived from the original waveform. The extracted IF revealed that the frequency of the fast wave was affected by both the propagation distance within the specimen and the bone density, independently. On the other hand, the alternation of the center frequency of the originally transmitted wave did not produce proportional changes in the extracted IF values of the fast waves, suggesting that the fast wave IF mainly reflected the thickness of the specimens. These findings may provide the possibility of obtaining a more precise diagnosis of osteoporosis.

Two-dimensional (2D) numerical simulations of vocal tract acoustics may provide a good balance between the high quality of three-dimensional (3D) finite element approaches and the low computational cost of one-dimensional (1D) techniques. However, 2D models are usually generated by considering the 2D vocal tract as a midsagittal cut of a 3D version, i.e., using the same radius function, wall impedance, glottal flow, and radiation losses as in 3D, which leads to strong discrepancies in the resulting vocal tract transfer functions. In this work, a four step methodology is proposed to match the behavior of 2D simulations with that of 3D vocal tracts with circular cross-sections. First, the 2D vocal tract profile becomes modified to tune the formant locations. Second, the 2D wall impedance is adjusted to fit the formant bandwidths. Third, the 2D glottal flow gets scaled to recover 3D pressure levels. Fourth and last, the 2D radiation model is tuned to match the 3D model following an optimization process. The procedure is tested for vowels /a/, /i/, and /u/ and the obtained results are compared with those of a full 3D simulation, a conventional 2D approach, and a 1D chain matrix model.

Deep back cavities are usually required for micro-perforated panel (MPP) constructions to achieve good low frequency absorption. To overcome the problem, a close-box loudspeaker with a shunted circuit is proposed to substitute the back wall of the cavity of the MPP constructions to constitute a composite absorber. Based on the equivalent circuit model, the acoustic impedance of the shunted loudspeaker is formulated first, then a prediction model of the sound absorption of the MPP backed by shunted loudspeaker is developed by employing the mode solution of a finite size MPP coupled by an air cavity with an impendence back wall. The MPP absorbs mid to high frequency sound, and with properly adjusted electrical parameters of its shunted circuit, the shunted loudspeaker absorbs low frequency sound, so the composite absorber provides a compact solution to broadband sound control. Numerical simulations and experiments are carried out to validate the model.

This article presents an inverse method for estimating the electromechanical parameters of a moving-coil loudspeaker with or without the eddy current and suspension creep effects. With known voice-coil displacement, voice-coil current, and stimulus signal as inputs, four calculation procedures for the direct problem, adjoint problem, sensitivity problem, and conjugate gradient method are involved in inversely solving the unknown electromechanical parameters. The proposed method features high efficiency in solving the direct problem through a hybrid spline difference method. It requires a small number of iterations for the computational algorithm, while offering excellent accuracy in parameter estimations. Analysis results demonstrate small differences between the estimated and measured electromechanical parameters under a variety of stimulus signals, excitation times, and initial guesses. The results are also confirmed by experimental measurements. These results indicate that the proposed method has a strong potential for estimating the electromechanical parameters of moving-coil loudspeakers.

Two effective medium models are presented and used to predict complex reflection and transmission coefficients of finite periodic arrays of resonant elastic shells as well as their effective density and bulk modulus at low frequencies. Comparisons with full multiple scattering theory and measurements show that the self-consistent model fails to correctly predict the shape of the transmission/reflection curves when scatterer resonances are close to the first Bragg bandgap. The low frequency grating model, which neglects the evanescent modes and considers scattered wave propagation only in the same direction as the incident one, gives a much better agreement with both measurements and the full multiple scattering theory. Moreover, because it does not require the wavelength to strongly exceed the size of scatterers, the model gives reliable predictions even at frequencies around the first periodicity related bandgap. In contrast to the self-consistent model, the low frequency grating model is applicable when the resonant scatterers have more than two low frequency resonances.

This article presents a simulation study on two integrated tuned vibration absorbers (TVAs) designed to control the global flexural vibration of lightly damped thin structures subject to broad frequency band disturbances. The first one consists of a single axial switching TVA composed by a seismic mass mounted on variable axial spring and damper elements so that the characteristic damping and natural frequency of the absorber can be switched iteratively to control the resonant response of three flexural modes of the hosting structure. The second one consists of a single three-axes TVA composed by a seismic mass mounted on axial and rotational springs and dampers, which are arranged in such a way that the suspended mass is characterized by uncoupled heave and pitch-rolling vibrations. In this case the three damping and natural frequency parameters of the absorber are tuned separately to control three flexural modes of the hosting structure. The simulation study shows that the proposed single-unit absorbers produce, respectively, 5.3 and 8.7 dB reductions of the global flexural vibration of a rectangular plate between 20 and 120 Hz.

Fully-sampled two-dimensional (2D) arrays can have two-way focusing of the ultrasound beam in both lateral directions leading to high quality, real-time three-dimensional (3D) imaging. However, fully-sampled 2D arrays with very large element counts (>16 000) are difficult to manufacture due to interconnect density and large element electrical impedance. As an alternative, row-column or crossed electrode arrays have been proposed to simplify transducer fabrication and system integration. These types of arrays consist of two one-dimensional arrays oriented perpendicular to each other. Using conventional delay-and-sum beamforming, each array performs one-way focusing in perpendicular lateral directions which yield higher sidelobe and acoustic clutter levels compared to fully-sampled 2D arrays with two-way focusing. In this paper, the use of spatial matched filters to improve focusing of row-column arrays is investigated. On receive, data from each element are first spatial match filtered in the elevation direction. After summation, the data are filtered again in the azimuth direction. Beam widths comparable to one-way focusing are seen in azimuth and beam widths comparable to two-way focusing are achieved in elevation. 3D beam patterns from computer simulation results using a 7.5 MHz 128 × 128 row-column array are shown with comparison to a fully sampled 2D array.

11 pages; International audience; Theoretical and experimental results are presented into the sound absorption and transmission properties of multi-layer structures made up of thin micro-perforated panels (ML-MPPs). The objective is to improve both the absorption and insulation performances of ML-MPPs through impedance boundary optimization. A fully coupled modal formulation is introduced that predicts the effect of the structural resonances onto the normal incidence absorption coefficient and transmission loss of ML-MPPs. This model is assessed against standing wave tube measurements and simulations based on impedance translation method for two double-layer MPP configurations of relevance in building acoustics and aeronautics. Optimal impedance relationships are proposed that ensure simultaneous maximization of both the absorption and the transmission loss under normal incidence. Exhaustive optimization of the double-layer MPPs is performed to assess the absorption and/or transmission performances with respect to the impedance criterion. It is investigated how the panel volumetric resonances modify the excess dissipation that can be achieved from non-modal optimization of ML-MPPs.

The main purpose of the paper is to contribute at presenting an analytical and a numerical modeling which would be relevant for interpreting the couplings between a circular membrane, a peripheral cavity having the same external radius as the membrane, and a thin air gap (with a geometrical discontinuity between them), and then to characterize small scale electrostatic receivers and to propose procedures that could be suitable for fitting adjustable parameters to achieve optimal behavior in terms of sensitivity and bandwidth expected. Therefore, comparison between these theoretical methods and characterization of several shapes is dealt with, which show that the models would be appropriate to address the design of such transducers.

Microbubble clustering may occur when bubbles become bound to targeted surfaces or are grouped by acoustic radiation forces in medical diagnostic applications. The ability to identify the formation of such clusters from the ultrasound echoes may be of practical use. Nonlinear numerical simulations were performed on clusters of microbubbles modeled by the modified Keller-Miksis equations. Encapsulated bubbles were considered to mimic practical applications but the aim of the study was to examine the effects of inter-bubble spacing and bubble size on the dynamical behavior of the cluster and to see if chaotic or bifurcation characteristics could be helpful in diagnostics. It was found that as microbubbles were clustered closer together, their oscillation amplitude for a given applied ultrasound power was reduced, and for inter-bubble spacing smaller than about ten bubble radii nonlinear subharmonics and ultraharmonics were eliminated. For clustered microbubbles, as for isolated microbubbles, an increase in the applied acoustic power caused bifurcations and transition to chaos. The bifurcations preceding chaotic behavior were identified by Floquet analysis and confirmed to be of the period-doubling type. It was found that as the number of microbubbles in a cluster increased, regularization occurred at lower ultrasound power and more windows of order appeared.

In this article, interleaved finite difference time domain methods are developed for the purpose of simulating the dynamics of the acoustic bore, using, as a starting point, an impedance formulation of wave propagation in an acoustic tube; attention is focused here on modeling of viscothermal and radiation losses in the time domain. In particular, in contrast to other methods, the bore, including the mouth-piece and bell, is treated as a unit, and is not subdivided into smaller units such as cylindrical or conical segments. Numerical simulations of input impedances are then compared with measurement for a variety of brass instruments.

International audience; This study investigates the influence of blood flow on temperature distribution during high-intensity focused ultrasound (HIFU) ablation of liver tumors. A three-dimensional acoustic-thermal-hydrodynamic coupling model is developed to compute the temperature field in the hepatic cancerous region. The model is based on the nonlinear Westervelt equation, bioheat equations for the perfused tissue and blood flow domains. The nonlinear Navier-Stokes equations are employed to describe the flow in large blood vessels. The effect of acoustic streaming is also taken into account in the present HIFU simulation study. A simulation of the Westervelt equation requires a prohibitively large amount of computer resources. Therefore a sixth-order accurate acoustic scheme in three-point stencil was developed for effectively solving the nonlinear wave equation. Results show that focused ultrasound beam with the peak intensity 2470 W/cm2 can induce acoustic streaming velocities up to 75 cm/s in the vessel with a diameter of 3 mm. The predicted temperature difference for the cases considered with and without acoustic streaming effect is 13.5 °C or 81% on the blood vessel wall for the vein. Tumor necrosis was studied in a region close to major vessels. The theoretical feasibility to safely necrotize the tumors close to major hepatic arteries and veins was shown.

In the spring of 2009, broadband transmissions from a ship-suspended source with a 284-Hz center frequency were received on a moored and navigated vertical array of hydrophones over a range of 107 km in the Philippine Sea. During a 60-h period over 19,000 transmissions were carried out. The observed wavefront arrival structure reveals four distinct purely refracted acoustic paths: One with a single upper turning point near 80 m depth, two with a pair of upper turning points at a depth of roughly 300 m, and one with three upper turning points at 420 m. Individual path intensity, defined as the absolute square of the center frequency Fourier component for that arrival, was estimated over the 60-h duration and used to compute scintillation index and log-intensity variance. Monte Carlo parabolic equation simulations using internal-wave induced sound speed perturbations obeying the Garrett-Munk internal-wave energy spectrum were in agreement with measured data for the three deeper-turning paths but differed by as much as a factor of four for the near surface-interacting path.

International audience; The generation of ultra-low frequency acoustic noise (0.1 to 1 Hz) by the nonlinear interaction of ocean surface gravity waves is well established. More controversial are the quantitative theories that attempt to predict the recorded noise levels and their variability. Here a single theoretical framework is used to predict the noise level associated with propagating pseudo-Rayleigh modes and evanescent acoustic-gravity modes. The latter are dominant only within 200 m from the sea surface, in shallow or deep water. At depths larger than 500 m, the comparison of a numerical noise model with hydrophone records from two open-ocean sites near Hawaii and the Kerguelen islands reveal: (a) Deep ocean acoustic noise at frequencies 0.1 to 1 Hz is consistent with the Rayleigh wave theory, in which the presence of the ocean bottom amplifies the noise by 10 to 20 dB; (b) in agreement with previous results, the local maxima in the noise spectrum support the theoretical prediction for the vertical structure of acoustic modes; and (c) noise level and variability are well predicted for frequencies up to 0.4 Hz. Above 0.6 Hz, the model results are less accurate, probably due to the poor estimation of the directional properties of wind-waves with frequencies higher than 0.3 Hz.

Infrasound can propagate very long distances and remain at measurable levels. As a result infrasound sensing is used for remote monitoring in many applications. At local ranges, on the order of 10 km, the influence of the presence or absence of forests on the propagation of infrasonic signals is considered. Because the wavelengths of interest are much larger than the scale of individual components, the forest is modeled as a porous material. This approximation is developed starting with the relaxation model of porous materials. This representation is then incorporated into a Crank-Nicholson method parabolic equation solver to determine the relative impacts of the physical parameters of a forest (trunk size and basal area), the presence of gaps/trees in otherwise continuous forest/open terrain, and the effects of meteorology coupled with the porous layer. Finally, the simulations are compared to experimental data from a 10.9 kg blast propagated 14.5 km. Comparison to the experimental data shows that appropriate inclusion of a forest layer along the propagation path provides a closer fit to the data than solely changing the ground type across the frequency range from 1 to 30 Hz.

A pressure compensated metal diaphragm based fiber laser hydrophone configuration that can provide good sensitivity, large bandwidth, and sea state zero noise floor is proposed in this paper. A simplified theoretical model of the proposed sensor configuration is developed in which the acoustic elements of the sensor configuration are modeled using a four-pole acoustic transfer matrix and the structural elements are modeled as second order single degree of freedom elements. This model is then used to optimize the design parameters of the sensor system to achieve the performance objectives. An axisymmetric finite element analysis of the sensor configuration is also carried out to validate the results from the simplified theoretical model. Prototype sensors were fabricated and hydrostatic testing in a pressure vessel validated the static pressure compensation performance of the sensor. Frequency dependent sensitivity of the sensor system was measured through acoustic testing in a water tank. The prototype sensor gave a flat frequency response up to 5 kHz and experimental results compared well with theoretical predictions. The sensor has an acceleration rejection figure on the order of 0 dB ref 1 m/s(2) Pa and the pressure compensation approach worked reasonably well up to a hydrostatic pressures equivalent to a depth of 50 m.

The atmospheric wind and temperature can be estimated through the traveltimes of infrasound between pairs of receivers. The traveltimes can be obtained by infrasonic interferometry. In this study, the theory of infrasonic interferometry is verified and applied to modeled stratospherically refracted waves. Synthetic barograms are generated using a raytracing model and taking into account atmospheric attenuation, geometrical spreading, and phase shifts due to caustics. Two types of source wavelets are implemented for the experiments: blast waves and microbaroms. In both numerical experiments, the traveltimes between the receivers are accurately retrieved by applying interferometry to the synthetic barograms. It is shown that microbaroms can be used in practice to obtain the traveltimes of infrasound through the stratosphere, which forms the basis for retrieving the wind and temperature profiles.

Deep ocean acoustics, in the absence of shipping and wildlife, is driven by surface processes. Best understood is the signal generated by non-linear surface wave interactions, the Longuet-Higgins mechanism, which dominates from 0.1 to 10 Hz, and may be significant for another octave. For this source, the spectral matrix of pressure and vector velocity is derived for points near the bottom of a deep ocean resting on an elastic half-space. In the absence of a bottom, the ratios of matrix elements are universal constants. Bottom effects vitiate the usual "standing wave approximation," but a weaker form of the approximation is shown to hold, and this is used for numerical calculations. In the weak standing wave approximation, the ratios of matrix elements are independent of the surface wave spectrum, but depend on frequency and the propagation environment. Data from the Hawaii-2 Observatory are in excellent accord with the theory for frequencies between 0.1 and 1 Hz, less so at higher frequencies. Insensitivity of the spectral ratios to wind, and presumably waves, is indeed observed in the data., Comment: 14 pages, 10 figures

Wave source localization from a sensor array has long been the most active research topics in both theory and application. In this paper, an explicit and time-domain inversion method for the direction and distance of a monopole source from a circular array is proposed. The approach is based on a mathematical technique, the weighted integral method, for signal/source parameter estimation. It begins with an exact form of the source-constraint partial differential equation that describes the unilateral propagation of wide-band waves from a single source, and leads to exact algebraic equations that include circular Fourier coefficients (phase mode measurements) as their coefficients. From them, nearly closed-form, single-shot and multishot algorithms are obtained that is suitable for use with band-pass/differential filter banks. Numerical evaluation and several experimental results obtained using a 16-element circular microphone array are presented to verify the validity of the proposed method.

Clutter is related to false alarms for active sonar. It is demonstrated that, in shallow water, target-like clutter in reverberation signals can be caused by nonlinear internal waves. A nonlinear internal wave is modeled using measured stratification on the New Jersey shelf. Reverberation in the presence of the internal wave is modeled numerically. Calculations show that acoustic energy propagating near a sound speed minimum is deflected as a high intensity, higher angle beam into the bottom, where it is backscattered along the reciprocal path. The interaction of sound with the internal wave is isolated in space, hence resulting in a target-like clutter, which is found to be greater than 10 dB above the mean reverberation level.

Broadband underwater acoustic communication signals undergo either a compression or dilation in the presence of relative motion between a source and a receiver. Consequently, underwater acoustic communications with a mobile source/receiver require Doppler compensation through resampling. However, resampling may not be necessary when a channel-estimate-based time reversal approach is applied with frequent channel updates. Using experimental data (20-30 kHz), it is demonstrated that the performance of time reversal communication without resampling is similar to the case with resampling, along with the benefit of a modest computational saving.

Measurements of acoustic ray travel-times in the ocean provide synoptic integrals of the ocean state between source and receiver. It is known that the ray travel-time is sensitive to variations in the ocean at the transmission time, but the sensitivity of the travel-time to spatial variations in the ocean prior to the acoustic transmission have not been quantified. This study examines the sensitivity of ray travel-time to the temporally and spatially evolving ocean state in the Philippine Sea using the adjoint of a numerical model. A one year series of five day backward integrations of the adjoint model quantify the sensitivity of travel-times to varying dynamics that can alter the travel-time of a 611 km ray by 200 ms. The early evolution of the sensitivities reveals high-mode internal waves that dissipate quickly, leaving the lowest three modes, providing a connection to variations in the internal tide generation prior to the sample time. They are also strongly sensitive to advective effects that alter density along the ray path. These sensitivities reveal how travel-time measurements are affected by both nearby and distant waters. Temporal nonlinearity of the sensitivities suggests that prior knowledge of the ocean state is necessary to exploit the travel-time observations.

This letter investigates the use of multiple line arrays (LAs) in a Time-Reversal Mirror for localizing and characterizing multipole aeroacoustic sources in a uniform subsonic mean flow using a numerical Time-Reversal (TR) method. Regardless of the original source characteristics, accuracy of predicting the source location can be significantly improved using at least two LAs. Furthermore, it is impossible to determine the source characteristics using a single LA, rather a minimum of two are required to establish either the monopole or dipole source nature, while four LAs (fully surrounding the source) are required for characterizing a lateral quadrupole source.

Killer whales in British Columbia are at risk, and little is known about their winter distribution. Passive acoustic monitoring of their year-round habitat is a valuable supplemental method to traditional visual and photographic surveys. However, long-term acoustic studies of odontocetes have some limitations, including the generation of large amounts of data that require highly time-consuming processing. There is a need to develop tools and protocols to maximize the efficiency of such studies. Here, two types of analysis, real-time and long term spectral averages, were compared to assess their performance at detecting killer whale calls in long-term acoustic recordings. In addition, two different duty cycles, 1/3 and 2/3, were tested. Both the use of long term spectral averages and a lower duty cycle resulted in a decrease in call detection and positive pod identification, leading to underestimations of the amount of time the whales were present. The impact of these limitations should be considered in future killer whale acoustic surveys. A compromise between a lower resolution data processing method and a higher duty cycle is suggested for maximum methodological efficiency.