Finite-difference time-domain solution of second-order photoacoustic wave equation

A finite-difference time-domain numerical solution is presented for solving a single second-order photoacoustic equation, instead of solving three coupled first-order equations. In this way, we are able to insert the heating function to the simulation directly instead of initial pressure. Results are validated using k-Wave simulation and show a good agreement for future development. The perfectly matched layer boundary condition has been implemented for a second-order photoacoustic equation and results are compared to Dirichlet, Neumann and Mur boundary conditions.