Fast calculation of pulsed photoacoustic fields in fluids usingk-space methods

Abstract

Two related numerical models that calculate the time-dependent pressure field radiated by an arbitrary photoacoustic source in a fluid, such as that generated by the absorption of a short laser pulse, are presented. Frequency-wavenumber (k-space) implementations have been used to produce fast and accurate predictions. Model I calculates the field everywhere at any instant of time, and is useful for visualizing the three-dimensional evolution of the wave field. Model II calculates pressure time series for points on a straight line or plane and is therefore useful for simulating array measurements. By mapping the vertical wavenumber spectrum directly to frequency, this model can calculate time series up to 50 times faster than current numerical models of photoacoustic propagation. As the propagating and evanescent parts of the field are calculated separately, model II can be used to calculate far- and near-field radiation patterns. Also, it can readily be adapted to calculate the velocity potential and thus particle velocity and acoustic intensity vectors. Both models exploit the efficiency of the fast Fourier transform, and can include the frequency-dependent directional response of an acoustic detector straightforwardly. The models were verified by comparison with a known analytic solution and a slower, but well-understood, numerical model.

An optimized ultrasound detector for photoacoustic breast tomography

“…To study the resolution and visibility of objects with a photoacoustic tomography (PAT) system employing the optimized detector, and comparing performances using detectors described in the literature, numerical simulations are performed using k-wave Matlab Tool-280 box [36]. Three detectors are simulated, from:…”

A new design of light illumination scheme for deep tissue photoacoustic imaging

“…Several numerical approaches, such as the finite-difference time domain (FDTD) method, the finite element (FE) method, and the pseudospectral (PS) method had also been used to simulate PA in different applications. In efforts to provide a fast and accurate prediction of the PA pressure field, a k-space method was proposed [11]. The thermal-pressure equation [12] was also employed to show the close relationship between the electromagnetic absorption and acoustic wave, but it did not fully describe the cause-end effect among several physical phenomena of PA.…”

Photoacoustic clutter reduction by inversion of a linear scatter model using plane wave ultrasound measurements

“…It is assumed that the wave is only scattered once (Born approximation). In the spatial frequency domain, the temporally evolving photoacoustic wave for t ≥ 0 can be described by [16,17]:…”

“…It has been shown in several publications that, in a homogeneous medium, the photoacoustic source p 0 relates to the photoacoustic measurement p h as [15][16][17][18]:…”