http://www.k-wave.org/forum/topic/about-the-wave-equation Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:07:30 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/about-the-wave-equation#post-8419 Sun, 09 Jan 2022 17:37:32 +0000 Amita 8419@http://www.k-wave.org/forum/ <p>ok<br /> Thank you </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-8418 Sun, 09 Jan 2022 15:23:52 +0000 Bradley Treeby 8418@http://www.k-wave.org/forum/ <p>See <a href="https://discovery.ucl.ac.uk/id/eprint/10113895/1/2014_Treeby_IEEE_kWaveElastic.pdf">here</a>. </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-8417 Sun, 09 Jan 2022 15:16:37 +0000 Amita 8417@http://www.k-wave.org/forum/ <p>Hi</p> <p>It is written there in user manual that the tool is solving governing equations for fluid medium, then how k-wave is applicable for solid or viscoelastic medium when the wave equation itself is not for these medium?</p> <p>or</p> <p>which governing equations are being solved behind the command pstdElastic2D and pstdElastic3D?</p> <p>please help. </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-4988 Thu, 19 Feb 2015 19:39:20 +0000 Bradley Treeby 4988@http://www.k-wave.org/forum/ <p>Hi Amin,</p> <p>There are a few reasons why we solve the first-order equations rather than the second-order wave equation in the main simulation functions. First, it allows both mass and force sources to be easily included into the discrete equations. Second, it allows the PML to be easily defined. Third, calculation of the particle velocity allows quantities such as the acoustic intensity to be calculated. This is useful, for example, when modelling how ultrasound heats biological tissue due to acoustic absorption. Finally, there is some evidence that solving the first-order equations gives improved accuracy for heterogeneous media.</p> <p>If you're interested in homogeneous media, the function <code>kspaceSecondOrder</code> solves the second-order wave equation.</p> <p>Hope that helps,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-4985 Wed, 11 Feb 2015 13:52:26 +0000 aminzadeh 4985@http://www.k-wave.org/forum/ <p>Dear Dr. Treeby</p> <p>Hello and Thanks for your Toolbox.</p> <p>As you said k-Wave solve a system of coupled 1st-order equations. why you don't solve the single 2nd-order equation?<br /> I asked this question, because before I know about your Toolbox, I solved the single 2nd order wave eq. and now I see that results are in a good agreement with the k-Wave. </p> <p>Best regards,<br /> Amin </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-321 Wed, 22 Feb 2012 02:39:18 +0000 xiong 321@http://www.k-wave.org/forum/ <p>Hi Brad,</p> <p>This is exactly what I want.</p> <p>Thank you!</p> <p>Xiong </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-314 Fri, 17 Feb 2012 05:34:24 +0000 Bradley Treeby 314@http://www.k-wave.org/forum/ <p>Hi Xiong,</p> <p>Thanks for your question. You are correct about the units. The alpha(r) in Chao's equation actually represents a relaxation time. This is normally denoted tau and is often called Stokes' relaxation time. In the limit of w*tau &lt;&lt; 1, the absorption will be proportional to frequency squared, in other words:</p> <p> absorption = alpha_0 w^2.</p> <p>In this case, the relationship between tau and alpha_0 is given by: </p> <p> alpha_0 = tau/(2*c) </p> <p>Here c is the sound speed in m/s and w is the frequency in rad/s. Note, this is not the wave equation solved in k-Wave.</p> <p>I hope that helps!</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-305 Fri, 10 Feb 2012 00:45:10 +0000 xiong 305@http://www.k-wave.org/forum/ <p>Dear Chao and Ben,</p> <p>I also saw the equation Chao posted. I have a question about it. The unit of alpha(r) in it should be the same as that of time (second) to ensure the unit consistency of the equation. But the acoustic attenuation coefficient (use beta(r) to denote it) usually has the unit of Np/cm or Np/cm/MHz. I don't know the relationship between alpha(r) (ultrasonic absorption distribution) and beta(r) (acoustic attenuation coefficient). If possible, please provide the corresponding equation with beta(r) represents the loss of acoustic wave.</p> <p>Thank you.</p> <p>Xiong </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-101 Wed, 16 Feb 2011 11:48:20 +0000 bencox 101@http://www.k-wave.org/forum/ <p>Dear Chao,</p> <p>I think deconvolving your frequency response before you do the image reconstruction is probably the safest way, as you can then deal with the problems that the deconvolution presents (noise amplification essentially) before it gets tangled up in an image reconstruction. Deconvolution of a frequency response from a time series will be covered in many standard digital signal processing books.</p> <p>Kind regards,</p> <p>Ben </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-100 Wed, 16 Feb 2011 09:29:12 +0000 huangchao 100@http://www.k-wave.org/forum/ <p>Hi Dr. Treeby,</p> <p>Thanks for the clarification. </p> <p>One more question: If the signal is the convolution of the simulated pressure data and the sensor frequency response, and I want to use this signal for reconstruction by time reversal method, other than direct deconvolution before reconstruction (FFT and IFFT), is there any other way to remove the frequency response effect during reconstruction process?</p> <p>Thanks,<br /> Chao </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-99 Wed, 16 Feb 2011 05:39:26 +0000 Bradley Treeby 99@http://www.k-wave.org/forum/ <p>Hi Chao,</p> <p>You are more than welcome. At this stage, the toolbox only supports scalar inputs for the absorption variables. We are planning to modify this to accommodate matrix inputs in the near future.</p> <p>The sensor frequency response simply takes the FFT of the recorded time domain signals after the simulation is complete, multiplies the frequency domain signals by a zero phase Gaussian window, and then computes the inverse FFT (see kspaceFirstOrder_filterSensorData.m within the private folder).</p> <p>To implement another type of frequency response, you could manually filter the recorded time domain signals in a similar manner (or in any other way you like) after they are returned from the simulation.</p> <p>I hope that is of some help,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-98 Wed, 16 Feb 2011 05:23:23 +0000 huangchao 98@http://www.k-wave.org/forum/ <p>Hi, Dr. Treeby,</p> <p>Thank you for your explanation. It's very helpful.</p> <p>And I have two other questions:</p> <p>1. It seems that the absorption parameters, like alpha_power, alpha_coeff, can not be given in the matrix form. For example, when I used matrix to define the absorption parameters in the "2D Time Reversal With Compensation For Acoustic Absorption Example", it will prompt error. </p> <p>2. Can I define the sensor frequency response in other forms instead of Gaussian shape?</p> <p>If not, hope the toolbox will support these two functions in the future.</p> <p>Thanks,<br /> Chao </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-97 Tue, 15 Feb 2011 22:22:59 +0000 Bradley Treeby 97@http://www.k-wave.org/forum/ <p>Hi Chao,</p> <p>Thanks for your question. The first order simulation functions in k-Wave (kspaceFirstOrder1D, 2D, and 3D) are based on three coupled first equations as you say. Without the absorption term, these are given by:</p> <p><img src="http://www.k-wave.org/images/kwave_first_order_equations.gif" /></p> <p>Here p and rho are acoustic quantities, u and d are the particle velocity and displacement, and rho_0 and c_0 are the background parameters.</p> <p>These equations are valid for both homogeneous or heterogeneous media (in rho_0 and c_0). If you combine these in the normal fashion including a photoacoustic source, you will get the equation you mention without the absorption term (the absorption model used in k-Wave is based on the fractional Laplacian, see the references for more details).</p> <p>Note, within the actual k-Wave code, the grad(rho_0) terms are not computed because they exactly cancel. (More subtly, this means the calculated acoustic density is not quite the acoustic density, but that doesn't matter in this case because this parameter is not used for input or output.)</p> <p>For more information - see the Modelling Wave Propagation section under Getting Started in the k-Wave help files.</p> <p>I hope that helps. If not, please let us know!</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/about-the-wave-equation#post-96 Tue, 15 Feb 2011 21:49:11 +0000 huangchao 96@http://www.k-wave.org/forum/ <p>Hi,</p> <p>Thanks for this wonderful toolbox!</p> <p>The k-Wave solves 3 first order equations (or second order wave equation for homogeneous media), which, at least in my opinion, is valid when the ambient density rho_0 is constant.</p> <p>Therefore, I am wondering if you can solve the wave equation that models the sound speed, density and absorption heterogeneities. To my knowledge, it should be </p> <p><img src="http://img225.imageshack.us/img225/6199/86132079.jpg" /></p> <p>where α(r) is the ultrasonic absorption distribution.</p> <p>Thanks,<br /> Chao</p> <p>Edit: I don't know why the image does not show properly, but you can find it here:<br /> <a href="http://img225.imageshack.us/f/86132079.jpg/" rel="nofollow">http://img225.imageshack.us/f/86132079.jpg/</a> (ADMIN: Fixed - you need the direct image link) </p>