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

Nonlinear 3-D simulation of high-intensity focused ultrasound therapy in the Kidney

Kidney cancer is a severe disease which can be treated non-invasively using high-intensity focused ultrasound (HIFU) therapy. However, tissue in front of the transducer and the deep location of kidney can cause significant losses to the efficiency of the treatment. The effect of attenuation, refraction and reflection due to different tissue types on HIFU therapy of the kidney was studied using a nonlinear ultrasound simulation model. The geometry of the tissue was derived from a computed tomography (CT) dataset of a patient which had been segmented for water, bone, soft tissue, fat and kidney. The combined effect of inhomogeneous attenuation and soundspeed was found to result in an 11.0 dB drop in spatial peak-temporal average (SPTA) intensity in the kidney compared to pure water. The simulation without refraction effects showed a 6.3 dB decrease indicating that both attenuation and refraction contribute to the loss in focal intensity. The losses due to reflections at soft tissue interfaces were less than 0.1 dB. Focal point shifting due to refraction effects resulted in -1.3, 2.6 and 1.3 mm displacements in x-, y- and z-directions respectively. Furthermore, focal point splitting into several smaller subvolumes was observed. The total volume of the secondary focal points was approximately 46% of the largest primary focal point. This could potentially lead to undesired heating outside the target location and longer therapy times

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

Simulating Focused Ultrasound Transducers using Discrete Sources on Regular Cartesian Grids

Accurately representing the behaviour of acoustic sources is an important part of ultrasound simulation. This is particularly challenging in ultrasound therapy where multielement arrays are often used. Typically, sources are defined as a boundary condition over a 2D plane within the computational model. However, this approach can become difficult to apply to arrays with multiple elements distributed over a non-planar surface. In this work, a grid-based discrete source model for single and multi-element bowl-shaped transducers is developed to model the source geometry explicitly within a regular Cartesian grid. For each element, the source model is defined as a symmetric, simply-connected surface with a single grid point thickness. Simulations using the source model with the opensource k-Wave toolbox are validated using the Rayleigh integral, O'Neil's solution, and experimental measurements of a focused bowl transducer under both quasi continuous wave and pulsed excitation. Close agreement is shown between the discrete bowl model and the axial pressure predicted by O'Neil's solution for a uniform curved radiator, even at very low grid resolutions. Excellent agreement is also shown between the discrete bowl model and experimental measurements. To accurately reproduce the near-field pressure measured experimentally, it is necessary to derive the drive signal at each grid point of the bowl model directly using holography. However, good agreement is also obtained in the focal region using uniformly radiating monopole sources distributed over the bowl surface. This allows the response of multi-element transducers to be modelled, even where measurement of an input plane is not possible

On the Adjoint Operator in Photoacoustic Tomography

Photoacoustic Tomography (PAT) is an emerging biomedical "imaging from coupled physics" technique, in which the image contrast is due to optical absorption, but the information is carried to the surface of the tissue as ultrasound pulses. Many algorithms and formulae for PAT image reconstruction have been proposed for the case when a complete data set is available. In many practical imaging scenarios, however, it is not possible to obtain the full data, or the data may be sub-sampled for faster data acquisition. In such cases, image reconstruction algorithms that can incorporate prior knowledge to ameliorate the loss of data are required. Hence, recently there has been an increased interest in using variational image reconstruction. A crucial ingredient for the application of these techniques is the adjoint of the PAT forward operator, which is described in this article from physical, theoretical and numerical perspectives. First, a simple mathematical derivation of the adjoint of the PAT forward operator in the continuous framework is presented. Then, an efficient numerical implementation of the adjoint using a k-space time domain wave propagation model is described and illustrated in the context of variational PAT image reconstruction, on both 2D and 3D examples including inhomogeneous sound speed. The principal advantage of this analytical adjoint over an algebraic adjoint (obtained by taking the direct adjoint of the particular numerical forward scheme used) is that it can be implemented using currently available fast wave propagation solvers.Comment: submitted to "Inverse Problems

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

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