An Optimization Approach for Creating Application-specific Ultrasound Speckle Tracking Algorithms.

OBJECTIVE: Ultrasound speckle tracking enables in vivo measurement of soft tissue deformation or strain, providing a non-invasive diagnostic tool to quantify tissue health. However, adoption into new fields is challenging since algorithms need to be tuned with gold-standard reference data that are expensive or impractical to acquire. Here, we present a novel optimization approach that only requires repeated measurements, which can be acquired for new applications where reference data might not be readily available or difficult to get hold of. METHODS: Soft tissue motion was captured using ultrasound for the medial collateral ligament (MCL) of three quasi-statically loaded porcine stifle joints, and medial ligamentous structures of a dynamically loaded human cadaveric knee joint. Using a training subset, custom speckle tracking algorithms were created for the porcine and human ligaments using surrogate optimization, which aimed to maximize repeatability by minimizing the normalized standard deviation of calculated strain maps for repeat measurements. An unseen test subset was then used to validate the tuned algorithms by comparing the ultrasound strains to digital image correlation (DIC) surface strains (porcine specimens) and length change values of the optically tracked ligament attachments (human specimens). RESULTS: After 1500 iterations, the optimization routine based on the porcine and human training data converged to similar values of normalized standard deviations of repeat strain maps (porcine: 0.19, human: 0.26). Ultrasound strains calculated for the independent test sets using the tuned algorithms closely matched the DIC measurements for the porcine quasi-static measurements (R > 0.99, RMSE < 0.59%) and the length change between the tracked ligament attachments for the dynamic human dataset (RMSE < 6.28%). Furthermore, strains in the medial ligamentous structures of the human specimen during flexion showed a strong correlation with anterior/posterior position on the ligaments (R > 0.91). CONCLUSION: Adjusting ultrasound speckle tracking algorithms using an optimization routine based on repeatability led to robust and reliable results with low RMSE for the medial ligamentous structures of the knee. This tool may be equally beneficial in other soft-tissue displacement or strain measurement applications and can assist in the development of novel ultrasonic diagnostic tools to assess soft tissue biomechanics.

Simultaneous transmission and reception on all elements of an array: binary code excitation.

Pulse-echo arrays are used in radar, sonar, seismic, medical and non-destructive evaluation. There is a trend to produce arrays with an ever-increasing number of elements. This trend presents two major challenges: (i) often the size of the elements is reduced resulting in a lower signal-to-noise ratio (SNR) and (ii) the time required to record all of the signals that correspond to every transmit-receive path increases. Coded sequences with good autocorrelation properties can increase the SNR while orthogonal sets can be used to simultaneously acquire all of the signals that correspond to every transmit-receive path. However, a central problem of conventional coded sequences is that they cannot achieve good autocorrelation and orthogonality properties simultaneously due to their length being limited by the location of the closest reflectors. In this paper, a solution to this problem is presented by using coded sequences that have receive intervals. The proposed approach can be more than one order of magnitude faster than conventional methods. In addition, binary excitation and quantization can be employed, which reduces the data throughput by roughly an order of magnitude and allows for higher sampling rates. While this concept is generally applicable to any field, a 16-element system was built to experimentally demonstrate this principle for the first time using a conventional medical ultrasound probe.

Scattering of near normal incidence SH waves by sinusoidal and rough surfaces in 3-D: comparison to the scalar wave approximation.

The challenge of accurately simulating how incident scalar waves interact with rough boundaries has made it an important area of research within many scientific disciplines. Conventional methods, which in the majority of cases focus only on scattering in two dimensions, often suffer from long simulation times or reduced accuracy, neglecting phenomena such as multiple scattering and surface self-shadowing. A simulation based on the scalar wave distributed point source method (DPSM) is presented as an alternative which is computationally more efficient than fully meshed numerical methods while obtaining greater accuracy than approximate analytical techniques. Comparison is made to simulated results obtained using the finite element method for a sinusoidally periodic surface where scattering only occurs in two dimensions, showing very good agreement (<0.2 dB). In addition to two-dimensional scattering, comparison to experimental results is also carried out for scattering in three dimensions when the surface has a Gaussian roughness distribution. Results indicate that for two-dimensional scattering and for rough surfaces with a correlation length equal to the incident wavelength (λ) and a root mean square height less than 0.2λ, the scalar wave approximation predicts reflected pulse shape change and envelope amplitudes generally to within 1 dB. Comparison between transducers within a three-element array also illustrate the sensitivity pulse amplitude can have to sensor position above a rough surface, differing by as much as 17 dB with a positional change of just 1.25λ.

Study on the bubble transport mechanism in an acoustic standing wave field.

The use of bubbles in applications such as surface chemistry, drug delivery, and ultrasonic cleaning etc. has been enormously popular in the past two decades. It has been recognized that acoustically-driven bubbles can be used to disturb the flow field near a boundary in order to accelerate physical or chemical reactions on the surface. The interactions between bubbles and a surface have been studied experimentally and analytically. However, most of the investigations focused on violently oscillating bubbles (also known as cavitation bubble), less attention has been given to understand the interactions between moderately oscillating bubbles and a boundary. Moreover, cavitation bubbles were normally generated in situ by a high intensity laser beam, little experimental work has been carried out to study the translational trajectory of a moderately oscillating bubble in an acoustic field and subsequent interactions with the surface. This paper describes the design of an ultrasonic test cell and explores the mechanism of bubble manipulation within the test cell. The test cell consists of a transducer, a liquid medium and a glass backing plate. The acoustic field within the multi-layered stack was designed in such a way that it was effectively one dimensional. This was then successfully simulated by a one dimensional network model. The model can accurately predict the impedance of the test cell as well as the mode shape (distribution of particle velocity and stress/pressure field) within the whole assembly. The mode shape of the stack was designed so that bubbles can be pushed from their injection point onto a backing glass plate. Bubble radial oscillation was simulated by a modified Keller-Miksis equation and bubble translational motion was derived from an equation obtained by applying Newton's second law to a bubble in a liquid medium. Results indicated that the bubble trajectory depends on the acoustic pressure amplitude and initial bubble size: an increase of pressure amplitude or a decrease of bubble size forces bubbles larger than their resonant size to arrive at the target plate at lower heights, while the trajectories of smaller bubbles are less influenced by these factors. The test cell is also suitable for testing the effects of drag force on the bubble motion and for studying the bubble behavior near a surface.

High-temperature (>500°c) wall thickness monitoring using dry-coupled ultrasonic waveguide transducers.

Conventional ultrasonic transducers cannot withstand high temperatures for two main reasons: the piezoelectric elements within them depolarize, and differential thermal expansion of the different materials within a transducer causes them to fail. In this paper, the design of a high-temperature ultrasonic thickness gauge that bypasses these problems is described. The system uses a waveguide to isolate the vulnerable transducer and piezoelectric elements from the high-temperature measurement zone. Use of thin and long waveguides of rectangular cross section allows large temperature gradients to be sustained over short distances without the need for additional cooling equipment. The main problems that had to be addressed were the transmission and reception of ultrasonic waves into and from the testpiece that the waveguides are coupled to, and optimization of the wave propagation along the waveguide itself. It was found that anti-plane shear loading performs best at transmitting and receiving from the surface of a component that is to be inspected. Therefore, a nondispersive guided wave mode in large-aspect-ratio rectangular strips was employed to transmit the anti-plane shear loading from the transducer to the measurement zone. Different joining methods to attach the waveguides to the component were investigated and experiments showed that clamping the waveguides to the component surface gave the best results. The thickness of different plate samples was consistently measured to within less than 0.1 mm. Performance at high temperatures was tested in a furnace at 730°C for 4 weeks without signal degradation. Thicknesses in the range of 3 to 25 mm could be monitored using Hanning windowed tonebursts with 2 MHz center frequency.