Abstract
The traditional Born series (TBS) and convergent Born series (CBS) methods to numerically solve the time-independent inhomogeneous photoacoustic (PA) wave equation are discussed. The performance of these algorithms is examined for a circular PA source (a disk of radius, $a = 5\,\,\unicode{x00B5}{\rm m}$) in two dimensions. The speed of sound within the source region was gradually decreased from ${v_s} = 1950$ to 1200 m/s, but the same quantity for the ambient medium was fixed to ${v_f} = 1500 \;{\rm m/s}$. The PA fields were calculated over a large frequency band from $f = 7.3$ to 2000 MHz. Accordingly, the wave number (${k_f} = 2\pi f/{v_f}$) varied from ${k_f} = 0.03$ to ${8.38}\;\unicode{x00B5}{{\rm m}^{- 1}}$. The TBS method does not offer converging solutions when ${k_f}a \ge 25$ for ${v_s} = 1950\; {\rm m/s}$ and ${k_f}a \ge 9$ for ${v_s} = 1200 \;{\rm m/s}$. These have been observed in both the near and far fields. However, the solutions for the CBS technique converge in all cases. Both methods facilitate accurate solutions if the computational domain contains a collection of monodisperse/polydisperse disks considered in this study. Our numerical results suggest that the CBS protocol can provide accurate solutions under various test conditions.
© 2020 Optical Society of America
Full Article | PDF ArticleMore Like This
Benjamin Krüger, Thomas Brenner, and Alwin Kienle
Opt. Express 25(21) 25165-25182 (2017)
Anuj Kaushik, Avijit Paul, and Ratan K. Saha
J. Opt. Soc. Am. A 37(12) 1895-1904 (2020)
R. E. Kleinman, G. F. Roach, and P. M. van den Berg
J. Opt. Soc. Am. A 7(5) 890-897 (1990)