Dear K-wave toolbox users and admins,

At the moment I am working on a project in which I want to reproduce the study by Cox and Beard (2009) on Photoacoustic tomography with a single detector in a reverberant cavity. I am particularly interested in the Modal Reconstruction method.

The coding can be split in two parts, first the simulation part and second the reconstruction part. As a first check I want to compare my simulated photoacoustic time series with the first part of a photo acoustic time series from a single source as presented in the study (fig 5). To do this I have simulated a rectangular reverberant cavity by defining a 3 cm thick metal box as the walls of my cavity. What I find is that acoustic pressure of my time series does decrease over time, while the time series presented in the paper does not. This is explained by the fact that they build a simulation assuming a reflection coefficient of 1, while in my case some of the acoustic pressure will propagate though the metal box and end up in the perfect matched layer.

Does anyone have an idea how to to model a reflection coefficient of 1 using the K-wave toolbox? Or will this inherently be impossible because of how the toolbox is build? And if so, what direction would you suggest me to investigate?

Thank you very much for your help!

Inge