http://www.k-wave.org/forum/topic/how-to-simulate-a-reverberant-cavity-with-reflection-coefficient-1 Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 00:19:29 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/how-to-simulate-a-reverberant-cavity-with-reflection-coefficient-1#post-8433 Fri, 21 Jan 2022 17:06:50 +0000 Inge 8433@http://www.k-wave.org/forum/ <p>Dear Ben,</p> <p>Thank you very much for your response and helpful suggestions! I really appreciate it. I will implement the reflectivity with your tricks.</p> <p>Regards,</p> <p>Inge </p> http://www.k-wave.org/forum/topic/how-to-simulate-a-reverberant-cavity-with-reflection-coefficient-1#post-8409 Thu, 06 Jan 2022 16:45:12 +0000 bencox 8409@http://www.k-wave.org/forum/ <p>Hi Inge, </p> <p>The trick has two parts. First take the initial pressure distribution on the grid that you wish to simulate. Expand it by flipping copies of it about the right-hand side and the top side (say). You will then have a domain that is twice as large in each direction, and has four copies of the original initial pressure distribution. Second, use periodic boundary conditions on the boundaries of that larger domain. The combination of these two will simulate all the reflections in the original (smaller) domain. </p> <p>Probably the easiest way to do this (and what I think we did in that paper) is to use the function <code>kspaceSecondOrder</code>, which has periodic boundaries by default.</p> <p>Hope that helps.</p> <p>Best wishes<br /> Ben </p> http://www.k-wave.org/forum/topic/how-to-simulate-a-reverberant-cavity-with-reflection-coefficient-1#post-8354 Tue, 07 Dec 2021 17:18:34 +0000 Inge 8354@http://www.k-wave.org/forum/ <p>Dear K-wave toolbox users and admins,</p> <p>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. </p> <p>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. </p> <p>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? </p> <p>Thank you very much for your help!</p> <p>Inge </p>