https://www.k-wave.org/forum/topic/simulate-thin-layered-phantom-using-shear-modulus-input-directly Support for the k-Wave MATLAB toolbox en-US Thu, 18 Jun 2026 06:12:44 +0000 http://bbpress.org/?v=1.0.2 q https://www.k-wave.org/forum/search.php https://www.k-wave.org/forum/topic/simulate-thin-layered-phantom-using-shear-modulus-input-directly#post-9180 Wed, 12 Feb 2025 03:25:52 +0000 TWC 9180@https://www.k-wave.org/forum/ <p>Thank you for providing such a convenient and free simulation tool.</p> <p>Currently, I am using <code>pstdElastic3D</code> to simulate shear wave propagation in thin-layered phantoms and to calculate shear wave speed. However, it is known that when the thickness of the layer is smaller than the wavelength of the shear wave, guided waves can occur. In such cases, the shear modulus cannot be derived using the simplified formula: shear modulus = density*(shear wave speed)^2.</p> <p>In <code>pstdElastic3D.m</code>, we input <code>medium.sound_speed_shear</code> to calculate lame parameter using<br /> <code>mu = medium.sound_speed_shear.^2 .* medium.density;</code>. </p> <p>I would like to ask if it is possible to:</p> <p>1. Modify the input to directly provide the shear modulus (<code>mu</code>) instead of <code>medium.sound_speed_shear</code> while maintaining simulation accuracy.<br /> 2. Set <code>flags.kelvin_voigt_model</code> to <code>false</code> to simulate an anisotropic, incompressible, and nonlinear thin layer.</p> <p>Best ragards,<br /> Ting-Wei </p>