http://www.k-wave.org/forum/topic/simulation-time Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:27:55 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/simulation-time#post-7084 Fri, 18 Oct 2019 19:09:52 +0000 Bradley Treeby 7084@http://www.k-wave.org/forum/ <p>Hi Michael,</p> <p>If the screen goes black, it's likely an instability. Unfortunately, the stability criterion for the elastic model means that with high absorption values, you need to use a small CFL. You can try minimising the simulation time by switching to single precision, reducing the size of your domain if possible, using different hardware, etc.</p> <p>Hope that helps,</p> <p>Brad </p> http://www.k-wave.org/forum/topic/simulation-time#post-7026 Wed, 28 Aug 2019 13:00:09 +0000 genny_p 7026@http://www.k-wave.org/forum/ <p>Hi Michael,</p> <p>I had a similar problem with a black screen forming (mine was due to some instability) when I increased the attenuation coefficients. I wonder if we are facing the same problem. In my case, I found the whether the simulation goes black is highly sensitive to the grid spacing. Increasing dx seemed to make the simulation more stable (as with increasing the cfl number), but then I couldn't confirm convergenge. I tried putting the parameters you supplied in the Shear Waves And Critical Angle Reflection Example and everything seems to run fine, even with the higher attenuation coefficients, so perhaps something else in your code is different (grid spacing/size probably considering the longer simulation time). Can you post your complete code (or a big enough subset of it to reproduce the problem) so I can see if your problem seems similar to mine? </p> <p>I posted my code at <a href="http://www.k-wave.org/forum/topic/instability-in-shear-wave-absorption-simulation" rel="nofollow">http://www.k-wave.org/forum/topic/instability-in-shear-wave-absorption-simulation</a> a few months ago, and am still searching for a solution. Perhaps you can try increasing the grid spacing in your simulation and see if a coarser grid helps your simulation like mine (though this would only confirm we are facing the same problem, and still need a solution to confirm for convergence).</p> <p>Best,<br /> Genny </p> http://www.k-wave.org/forum/topic/simulation-time#post-7025 Tue, 27 Aug 2019 18:23:23 +0000 elliotP 7025@http://www.k-wave.org/forum/ <p>Hi,</p> <p>I have a question about the relationship between elastic simulation time and coefficients.<br /> For instance, when I use below constants:<br /> -----------------------------------------<br /> % water<br /> cp1 = 1540; % compressional wave speed [m/s]<br /> cs1 = 0; % shear wave speed [m/s]<br /> rho1 = 1000; % density [kg/m^3]<br /> alpha0_p1 = 0; % compressional wave absorption [dB/(MHz^2 cm)]<br /> alpha0_s1 = 0; % shear wave absorption [dB/(MHz^2 cm)]<br /> % skull<br /> cp2 = 3000; % compressional wave speed [m/s]<br /> cs2 = 1500; % shear wave speed [m/s]<br /> rho2 = 1850; % density [kg/m^3]<br /> alpha0_p2 = 3; % compressional wave absorption [dB/(MHz^2 cm)]<br /> alpha0_s2 = 6; % shear wave absorption [dB/(MHz^2 cm)]<br /> % create the time array<br /> cfl = 0.13;<br /> t_end = 40e-6;<br /> -----------------------------------------<br /> estimated simulation time 1hours 15min 41.4135s...<br /> -----------------------------------------</p> <p>When I change the skull's attenuation values to alpha0_p2 = 8.83, alpha0_s2 = 19.5 (belongs to your paper(Modeling power law absorption and dispersion in viscoelastic solids using a split-field<br /> and the fractional Laplacian)). Unfortunately, simulation gives a black screen. I think I have to change the cfl number. But, when I decrease the cfl coefficient, simulation time takes much longer.<br /> How can I minimize the simulation time under these circumstances?</p> <p>Best,<br /> Michael </p>