The preliminary evaluation of a 1 MHz ultrasound probe for measuring the elastic anisotropy of human cortical bone

Our objective is to evaluate an ultrasound probe for measurements of velocity and anisotropy in human cortical bone (tibia). The anisotropy of cortical bone is a known and mechanically relevant property in the context of osteoporotic fracture risk. Current in vivo quantitative ultrasound devices measuring the velocity of ultrasound in long bones can only be applied in the axial direction. For anisotropy measurements a second direction for velocity measurements preferably perpendicular to the axial direction is necessary. We developed a new ultrasound probe which permits axial transmission measurements with a simultaneous second perpendicular direction (tangential). Anisotropy measurements were performed on isotropic and anisotropic phantoms and two excised human female tibiae (age 63 and 82). Anisotropy ratios (AI; ratio of squared ultrasound velocities in the two directions) were for the isotropic phantom 1.06 ± 0.01 and for the anisotropic phantom 1.14 ± 0.03 (mean ± standard deviation). AI was 1.83 ± 0.29 in the tibia from the older donor and 1.37 ± 0.18 in the tibia from the younger donor. The AIs were in the expected range and differed significantly (p < 0.05, t-test) between the tibiae. Measured sound velocities were reproducible (mean standard deviation of short time precision of both channels for phantom measurements 31 m/s) and in agreement with reported velocities of the phantom material. Our results document the feasibility of anisotropy measurements at long bones using a single probe. Further improvements in the design of the probe and tests in vivo are warranted. If this approach can be evaluated in vivo an additional tool for assessing the bone status is available for clinical use.     

Experimental study of underwater transmission characteristics of high-frequency 30 MHz polyurea ultrasonic transducer

In this paper, we present the transmission characteristics of a polyurea ultrasonic transducer operating in water. In this study, we used a polyurea transducer with fundamental resonance at approximately 30 MHz. Firstly, acoustic pressure radiated from the transducer was measured using a hydrophone, which has a diameter of 0.2 mm. The transmission characteristics such as relative bandwidth, pulse width, and acoustic sensitivity were calculated from the experimental results. The results of the experiment showed a relative bandwidth of 50% and a pulse width of 0.061 μs. The acoustic sensitivity was 0.60 kPa/V with good linearity, where the correlation coefficient R in the fitting calculation was 0.996. A maximum pressure of 13.1 kPa was observed when the transducer was excited at a zero-to-peak voltage of 21 V. Moreover, we experimentally verified the results. The results of the pulse/echo experiment showed that the estimated diameters of the copper wires were 458 and 726 μm, where the differences between the actual and measured values were 15% and 4%, respectively. Acoustic streaming was also observed so that a particle velocity map was estimated by particle image velocimetry (PIV). The sound pressure calculated from the particle velocity obtained by PIV showed good agreement with the acoustic pressure measured using the hydrophone, where the differences between the calculated and measured values were 12–19%.     

Numerical study of the synergistic effect of cavitation and micro-abrasive particles

In ultrasonic-assisted machining, the synergistic effect of the cavitation effect and micro-abrasive particles plays a crucial role. Studies have focused on the investigation of the micro-abrasive particles, cavitation micro-jets, and cavitation shock waves either individually or in pairs. To investigate the synergy of shock waves and micro-jets generated by cavitation with micro-abrasive particles in ultrasonic-assisted machining, the continuous control equations of a cavitation bubble, shock wave, micro-jet, and micro-abrasive particle influenced by the dimensionless amount (R/R₀), a particle size-velocity–pressure model of the micro-abrasive particle was established. The effects of ultrasonic frequency, sound pressure amplitude, and changes in particle size on micro-abrasive particle velocity and pressure were numerically simulated. At an ultrasonic frequency of 20 kHz and ultrasonic sound pressure of 0.1125 MPa, a smooth spherical SiO₂ micro-abrasive particle (size = 5 µm) was obtained, with a maximum velocity of 190.3–209.4 m/s and pressure of 79.69–89.41 MPa. The results show that in the range of 5–50 μm, smaller particle sizes of the micro-abrasive particles led to greater velocity and pressure. The shock waves, micro-jets, and micro-abrasive particles were all positively affected by the dimensionless amount (R/R₀) of cavitation bubble collapse, the larger the dimensionless quantity, the faster their velocity and the higher their pressure.     

Estimation of speed of sound in dual-layered media using medical ultrasound image deconvolution

The speed of sound in soft tissues is usually assumed to be 1540 m/s in medical pulse-echo ultrasound imaging systems. When the true speed is different, the mismatch can lead to distortions in the acquired images, and so reduce their clinical value. Previously we reported a new method of sound-speed estimation in the context of image deconvolution. Unlike most other sound-speed estimation methods, this enables the use of unmodified ultrasound machines and a normal scanning pattern. Our approach was validated for largely homogeneous media with single sound speeds. In this article, we demonstrate that sound speeds of dual-layered media can also be estimated through image deconvolution. An ultrasound simulator has been developed for layered media assuming that, for moderate speed differences, the reflection at the interface may be neglected. We have applied our dual-layer algorithm to simulations and in vitro phantoms. The speed of the top layer is estimated by our aforesaid method for homogeneous media. Then, when the layer boundary position is known, a series of deconvolutions are carried out with dual-layered point-spread functions having different lower-layer speeds. The best restoration is selected using a correlation metric. The error level (e.g., a mean error of −9 m/s with a standard deviation of 16 m/s) for in vitro phantoms is found to be not as good as that of our single-speed algorithm, but is comparable to other local speed estimation methods where the data acquisition may not be as simple as in our proposed method.     

Propagation of sound waves through a spatially homogeneous but smoothly time-dependent medium

The propagation of sound through a spatially homogeneous but non-stationary medium is investigated within the framework of fluid dynamics. For a non-vortical fluid, especially, a generalized wave equation is derived for the (scalar) potential of the fluid velocity distribution in dependence of the equilibrium mass density of the fluid and the sound wave velocity. A solution of this equation for a finite   transition period τ$\tau$ is determined in terms of the hypergeometric function for a phenomenologically realistic, sigmoidal change of the mass density and sound wave velocity. Using this solution, it is shown that the energy flux of the sound wave is not conserved but increases always   for the propagation through a non-stationary medium, independent of whether the equilibrium mass density is increased or decreased. It is found, moreover, that this amplification of the transmitted wave arises from an energy exchange with the medium and that its flux is equal to the (total) flux of the incident and the reflected wave. An interpretation of the reflected wave as a propagation of sound backward in time is given in close analogy to Feynman and Stueckelberg for the propagation of anti-particles. The reflection and transmission coefficients of sound propagating through a non-stationary medium is analyzed in more detail for hypersonic waves with transition periods τ$\tau$ between 15 and 200 ps as well as the transformation of infrasound waves in non-stationary oceans.     

Standardised acoustic characterisation of sonic crystals noise barriers: Sound insulation and reflection properties

This work presents sound insulation and sound reflection measurements conducted over sonic crystal noise barriers according to the European standards EN 1793-2, EN 1793-5 and EN 1793-6. In most of the reference literature, sound insulation and reflection properties of sonic crystals are measured or a diffuse sound field or in a direct sound field including the top and side edge diffraction effects together with the transmitted (or reflected) components. The aim of this work is to perform free-field measurements over a real-sized sample in order to window out all diffraction components and to verify the points of strength and weakness of the application of standardised measurements to sonic crystals. Diffuse field measurements in laboratory are also done for comparison purposes. Since the target frequency range for traffic noise spectrum is centred at around 1000 Hz, a finite element based parametric investigation is performed to design unit cells capable of generating band gaps in the one-third octave bands ranging from 800 Hz to 1250 Hz. Then, 3 × 3 m sonic crystal noise barriers are installed in the Laboratory of the University of Bologna and sound insulation and sound reflection measurements are performed according to the mentioned active standards for normal incidence. Sound insulation is measured for diffuse incidence too. The two methods give different results. The method more directly comparable to calculations is the free-field one. However, if on the one hand the application of a time window allows to compute the transmitted or reflected component only, on the other hand the time window itself limits the maximum width of the sample for which all reflections of the n-th order having a significant spectral content are included, and thus results critical in the analysis of this kind of noise barriers. Nevertheless, the standardised measurements allow a direct comparison between the performance of sonic crystals and common noise barriers.     

Noise characteristics of tractors and health effect on farmers

An investigation was conducted to determine the noise propagation (under stationary condition) and noise at operator’s ear level of popular 18.7 and 26.1 kW tractors and 4.6 and 6.7 kW hand tractors during field operations with various implements. It was observed that both tractors produced the noise of 92 dB(A) L_eq in the working zone of operator. The sound pressure level (SPL) of the hand tractor was about 2 dB(A) L_eq higher than that of the tractor. The SPL during field operations at operator’s ear level increased with increase in engine speed and forward speed. Furthermore, the SPL was higher for field operations corresponding to the implement requiring higher draft. It was observed that the SPLs of the tractors and hand tractors were more than the exposure limit of noise for 8-h workday recommended by ISO and OSHA. This may cause health problems to the farmers in the long run.     

3-D Model of sound pressure field in a meridinal section plane of fruit

A theoretical model was suggested for qualitative evaluation of a sound pressure field in fruit tissue, as affected by ultrasonic probe dimensions and fruit properties. The classic directivity pattern of an ideal fluid model, expressed by Bessel function of the first kind, was extended to include energy dissipation of a real material. The directional characteristics of wave propagation, as influenced by transmitter frequency and diameter, and by fruit properties, were discussed. The model indicates how to select the parameters of the ultrasonic transducer (transducer diameter, frequency and excitation power) to control the magnitude and directivity of the ultrasonic waves in the fruit tissue. The suggested theoretical model represented fairly well the experimental sound wave distribution over the half-cut surface of potato and avocado (R² > 0.862 and 0.977, respectively); the same theoretical model could not represent the sound wave distribution over a half-cut melon. Results of the study were applied in a new probe design for ultrasonic testing of whole fruit.     

Experimental implementation of a low-frequency global sound equalization method based on free field propagation

An experimental implementation of a global sound equalization method in a rectangular room using active control is described in this paper. The main purpose of the work has been to provide experimental evidence that sound can be equalized in a continuous three-dimensional region, the listening zone, which occupies a considerable part of the complete volume of the room. The equalization method, based on the simulation of a progressive plane wave, was implemented in a room with inner dimensions of 2.70 m × 2.74 m × 2.40 m. With this method, the sound was reproduced by a matrix of 4 × 5 loudspeakers in one of the walls. After traveling through the room, the sound wave was absorbed on the opposite wall, which had a similar arrangement of loudspeakers, by means of active control. A set of 40 digital FIR filters was used to modify the original input signal before it was fed to the loudspeakers, one filter for each transducer. The optimal arrangement of the loudspeakers and the maximum frequency that can be equalized is analyzed theoretically in this paper. The presented experimental results show that sound equalization was possible from 10 Hz to approximately 425 Hz in the listening zone. A flat frequency response with deviations within ±5 decibels from the desired value was achieved. A higher demanding performance with deviations within ±1.5 decibels from a flat frequency response was attained in the interval between 20 Hz and 280 Hz. At the same time, the impulse response was quite well approximated to a delayed delta function in the listening zone. Examples of the spatial distribution of the sound field are also shown.     

Theoretical and phantom based investigation of the impact of sound speed and backscatter variations on attenuation slope estimation

Quantitative ultrasound features such as the attenuation slope, sound speed and scatterer size, have been utilized to evaluate pathological variations in soft tissues such as the liver and breast. However, the impact of variations in the sound speed and backscatter due to underlying fat content or fibrotic changes, on the attenuation slope has not been addressed. Both numerical and acoustically uniform tissue-mimicking experimental phantoms are used to demonstrate the impact of sound speed variations on attenuation slope using clinical real-time ultrasound scanners equipped with linear array transducers. Radiofrequency data at center frequencies of 4 and 5 MHz are acquired for the experimental and numerical phantoms respectively. Numerical phantom sound speeds between 1480 and 1600 m/s in increments of 20 m/s for attenuation coefficients of 0.3, 0.4, 0.5, 0.6, and 0.7 dB/cm/MHz are simulated. Variations in the attenuation slope when the backscatter intensity of the sample is equal, 3 dB higher, and 3 dB lower than the reference is also evaluated. The sound speed for the experimental tissue-mimicking phantoms were 1500, 1540, 1560 and 1580 m/s respectively, with an attenuation coefficient of 0.5 dB/cm/MHz. Radiofrequency data is processed using three different attenuation estimation algorithms, i.e. the reference phantom, centroid downshift, and a hybrid method. In both numerical and experimental phantoms our results indicate a bias in attenuation slope estimates when the reference phantom sound speed is higher (overestimation) or lower (underestimation) than that of the sample. This bias is introduced via a small spectral shift in the normalized power spectra of the reference and sample with different sound speeds. The hybrid method provides the best estimation performance, especially for sample attenuation coefficient values lower than that of the reference phantom. The performance of all the methods deteriorates when the attenuation coefficient of the reference phantom is lower than that of the sample. In addition, the hybrid method is the least sensitive to sample backscatter intensity variations.