Design of multi-frequency acoustic kinoforms

Complex diffraction limited acoustic fields can be generated from a single element transducer using inexpensive 3-D printable acoustic kinoforms. This is extremely promising for a number of applications. However, the lack of ability to vary the field limits the potential use of this technology. In this work, this limitation is circumvented using multi-frequency acoustic kinoforms for which different acoustic fields are encoded onto different driving frequencies. An optimisation approach based on random downhill binary search is introduced for the design of the multi-frequency kinoforms. This is applied to two test cases to demonstrate the technique: a kinoform designed to generate the numerals “1,” “2,” and “3” in the same plane but at different driving frequencies, and a kinoform designed to generate 3 sets of eight foci lying on a circle with a driving-frequency-dependent radius. Field measurements from these samples confirmed that multi-frequency acoustic kinoforms can be designed that switch between different arbitrary, pre-designed, acoustic field patterns in the target plane by changing the driving frequenc

Generating arbitrary ultrasound fields with tailored optoacoustic surface profiles

Acoustic fields with multiple foci have many applications in physical acoustics ranging from particle manipulation to neural modulation. However, the generation of multiple foci at arbitrary locations in three-dimensional is challenging using conventional transducer technology. In this work, the optical generation of acoustic fields focused at multiple points using a single optical pulse is demonstrated. This is achieved using optically absorbing surface profiles designed to generate specific, user-defined, wavefields. An optimisation approach for the design of these tailored surface profiles is developed. This searches for a smoothly varying surface that will generate a high peak pressure at a set of target focal points. The designed surface profiles are then realised via a combination of additive manufacturing and absorber deposition techniques. Acoustic field measurements from a sample designed to generate the numeral “7” are used to demonstrate the design method

Control of broadband optically generated ultrasound pulses using binary amplitude holograms

In this work, the use of binary amplitude holography is investigated as a mechanism to focus broadband acoustic pulses generated by high peak-power pulsed lasers. Two algorithms are described for the calculation of the binary holograms; one using ray-tracing, and one using an optimization based on direct binary search. It is shown using numerical simulations that when a binary amplitude hologram is excited by a train of laser pulses at its design frequency, the acoustic field can be focused at a pre-determined distribution of points, including single and multiple focal points, and line and square foci. The numerical results are validated by acoustic field measurements from binary amplitude holograms, excited by a high peak-power laser

Single Pulse Illumination of Multi-Layer Photoacoustic Holograms for Patterned Ultrasound Field Generation

A new method for the creation of patterned, focused, optically generated acoustic fields using a single optical pulse is introduced. This utilises multi-layer `holograms' composed of several spatially separate absorbing layers. Each layer is individually patterned so as to focus at a set of targeted points. To create the patterns, a ray-tracing model was implemented to calculate the impulse response of pixels within each absorbing layer to a set of targeted points. An optimisation approach was then used to find the optimal pattern for each layer to create a field evenly focused at each of the target points. The method was validated using both numerical simulations and acoustic field measurements. It was demonstrated that a 3×3 array of acoustic foci could be generated from a 3-layer hologram using a single laser pulse

Rapid Spatial Mapping of Focused Ultrasound Fields Using a Planar Fabry-Pérot Sensor

Measurement of high acoustic pressures is necessary in order to fully characterise clinical high-intensity focused ultrasound (HIFU) fields, and for accurate validation of computational models of ultrasound propagation. However, many existing measurement devices are unable to withstand the extreme pressures generated in these fields, and those that can often exhibit low sensitivity. Here, a planar Fabry-Pérot interferometer with hard dielectric mirrors and spacer was designed, fabricated, and characterised and its suitability for measurement of nonlinear focused ultrasound fields was investigated. The noise equivalent pressure of the scanning system scaled with the adjustable pressure detection range between 49 kPa for pressures up to 8 MPa and 152 kPa for measurements up to 25 MPa, over a 125 MHz measurement bandwidth. Measurements of the frequency response of the sensor showed that it varied by less than 3 dB in the range 1 - 62 MHz. The effective element size of the sensor was 65 μm and waveforms were acquired at a rate of 200 Hz. The device was used to measure the acoustic pressure in the field of a 1.1 MHz single element spherically focused bowl transducer. Measurements of the acoustic field at low pressures compared well with measurements made using a PVDF needle hydrophone. At high pressures, the measured peak focal pressures agreed well with the focal pressure modelled using the Khokhlov-Zabolotskaya-Kuznetsov equation. Maximum peak positive pressures of 25 MPa, and peak negative pressures of 12 MPa were measured, and planar field scans were acquired in scan times on the order of 1 minute. The properties of the sensor and scanning system are well suited to measurement of nonlinear focused ultrasound fields, in both the focal region and the low pressure peripheral regions. The fast acquisition speed of the system and its low noise equivalent pressure are advantageous, and with further development of the sensor, it has potential in application to HIFU metrology

Accurate simulation of transcranial ultrasound propagation for ultrasonic neuromodulation and stimulation

Non-invasive, focal neurostimulation with ultrasound is a potentially powerful neuroscientific tool that requires effective transcranial focusing of ultrasound to develop. Time-reversal (TR) focusing using numerical simulations of transcranial ultrasound propagation can correct for the effect of the skull, but relies on accurate simulations. Here, focusing requirements for ultrasonic neurostimulation are established through a review of previously employed ultrasonic parameters, and consideration of deep brain targets. The specific limitations of finite-difference time domain (FDTD) and k-space corrected pseudospectral time domain (PSTD) schemes are tested numerically to establish the spatial points per wavelength and temporal points per period needed to achieve the desired accuracy while minimizing the computational burden. These criteria are confirmed through convergence testing of a fully simulated TR protocol using a virtual skull. The k-space PSTD scheme performed as well as, or better than, the widely used FDTD scheme across all individual error tests and in the convergence of large scale models, recommending it for use in simulated TR. Staircasing was shown to be the most serious source of error. Convergence testing indicated that higher sampling is required to achieve fine control of the pressure amplitude at the target than is needed for accurate spatial targeting

Time domain reconstruction of sound speed and attenuation in ultrasound computed tomography using full wave inversion

Ultrasound computed tomography (USCT) is a non-invasive imaging technique that provides information about the acoustic properties of soft tissues in the body, such as the speed of sound (SS) and acoustic attenuation (AA). Knowledge of these properties can improve the discrimination between benign and malignant masses, especially in breast cancer studies. Full wave inversion (FWI) methods for image reconstruction in USCT provide the best image quality compared to more approximate methods. Using FWI, the SS is usually recovered in the time domain, and the AA is usually recovered in the frequency domain. Nevertheless, as both properties can be obtained from the same data, it is desirable to have a common framework to reconstruct both distributions. In this work, an algorithm is proposed to reconstruct both the SS and AA distributions using a time domain FWI methodology based on the fractional Laplacian wave equation, an adjoint field formulation, and a gradient-descent method. The optimization code employs a Compute Unified Device Architecture version of the software k-Wave, which provides high computational efficiency. The performance of the method was evaluated using simulated noisy data from numerical breast phantoms. Errors were less than 0.5% in the recovered SS and 10% in the AA. V

ElasticMatrix: A MATLAB toolbox for anisotropic elastic wave propagation in layered media

Simulating the propagation of elastic waves in multi-layered media has many applications. A common approach is to use matrix methods where the elastic wave-field within each material layer is represented by a sum of partial-waves along with boundary conditions imposed at each interface. While these methods are well-known, coding the required matrix formation, inversion, and analysis for general multi-layered systems is non-trivial and time-consuming. Here, a new open-source toolbox called ElasticMatrix is described which solves the problem of acoustic and elastic wave propagation in multi-layered media for isotropic and transverse-isotropic materials where the wave propagation occurs in a material plane of symmetry. The toolbox is implemented in MATLAB using an object oriented programming framework and is designed to be easy to use and extend. Methods are provided for calculating and plotting dispersion curves, displacement and stress fields, reflection and transmission coefficients, and slowness profiles

Regularization of image reconstruction in ultrasound computed tomography

We propose two regularizations techniques for a bent-ray (BR) tracing algorithm to reconstruct the speed of sound maps of breast tissues in an Ultrasound Computed Tomography (USCT) system. When high frequencies are employed, the use of BR is a good approximation to describe the propagation of the front of the pressure wave. The quantitative accuracy of the images reconstructed with the BR algorithm was evaluated without any kind of regularization, and with two regularization methods. The regularizations were based on some available a priori information, namely the known higher and lower values of the speed of sound expected in the breast tissues, and the maps of the internal structures obtained from the standard reflection ultrasound (US) imaging. The use of the proposed regularizations in the implemented algorithm improves the convergence and quality of the resulting images, although further improvements are still possible. These methods will help obtaining quantitative US images in a reasonable amount of time, expanding the possibilities and applications of this technique