http://www.k-wave.org/forum/topic/scaling-fator-when-impose-dirichlet-boundary-conditions Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 01:03:52 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/scaling-fator-when-impose-dirichlet-boundary-conditions#post-9186 Wed, 19 Feb 2025 22:26:34 +0000 HD Y 9186@http://www.k-wave.org/forum/ <p>Hi all, </p> <p>when imposing pressure dirichlet boundary conditions, according to Section 2.2 Acoustic Source Terms in k-wave manual, equation (2.19), there is a scaling factor applied to source.p. </p> <p> 1/(c_02N) 2 c_0/\Delta x</p> <p>There is a constant 2 in the factor, is there any reason choosing this constant? </p> <p>In the implementation, kspaceFirstOrder3D.m, from line 930 on:</p> <p> switch source.p_mode<br /> case 'dirichlet'</p> <p> % enforce source values as a dirichlet boundary condition<br /> rhox(p_source_pos_index) = source.p(p_source_sig_index, t_index);<br /> rhoy(p_source_pos_index) = source.p(p_source_sig_index, t_index);<br /> rhoz(p_source_pos_index) = source.p(p_source_sig_index, t_index);</p> <p> .....</p> <p>However, here the source.p is directly assigned to the density without any scaling factor!<br /> Is there any thing missing in this source code? </p> <p>This seems to me a bug in the code. Can any developer provide some details on this? </p> <p>Thank you in advance. </p>