Software development
For theoretical study of nonlinear wave fields generated by powerful focusing ultrasound sources, various numerical models are developed in the LIMU Laboratory. Mainly, nonlinear equations of Khokhlov-Zabolotskaya-Kuznetsov (KZK) and Westervelt are used for modeling, which are applied to describe unidirectional wave propagation. A special feature of the developed algorithms is the ability to calculate fields in cases where high-amplitude shock fronts are formed in the focal area. Studying such modes is necessary when developing new shock-wave protocols for influencing biological tissues.
Based on the numerical algorithms for solving the above equations, corresponding software packages have been developed, one of which, “HIFU beam”, has been made available for free use by researchers around the world. The “HIFU beam” software allows for calculations of axially symmetric beams generated by classical single-element spherical sources and ring arrays, and can also be used for transducers of a more complex design using the equivalent source method. “HIFU beam” is equipped with a graphical user interface and can be run on a desktop computer using parallel computing on several processor cores. More resource-intensive 3D tasks are solved on high-performance servers.

Acceleration of numerical experiments using GPU
Among other things, LIMU is working on acceleration of existing algorithms using parallel computing on graphics processors (GPU). The developed software packages are actively used to solve practical problems of calibrating fields of powerful medical transducers.
Numerical algorithms for the problems of parametric sound radiation
The numerical algorithms developed in the LIMU Laboratory are used not only to solve problems of medical ultrasound. In particular, the interaction of two high-frequency pump waves is calculated using the KZK equation, allowing to predict generation of a difference frequency wave (the so-called parametric sound emitter). To perform three-dimensional calculations here, it was necessary to optimize the algorithms by taking into account only the most significant spectral components of the nonlinear wave.
Numerical experiments in aeroacoustics
Numerical studies are being conducted related to the propagation of nonlinear waves in randomly inhomogeneous media. In particular, the sonic boom wave from a supersonic aircraft after passing through the atmospheric turbulent layer near the ground acquires a random character and its characteristics must be described statistically. Modeling based on the KZK equation for large random realizations of the acoustic field makes it possible to predict the statistical characteristics of the noise level.
Equations being solved
- one-dimensional equations of nonlinear acoustics in media with frequency dependent absorption
- Khokhlov-Zabolotskaya-Kuznetsov (KZK) type evolution equations
- Westervelt equation
LIMU tasks
- Development of new algorithms
- Software development
- Acceleration of numerical experiments using graphics processors (GPU)
- Implementation of numerical experiments for specific problems of medical and industrial ultrasound:
- Manipulation of objects with ultrasound for kidney stone propulsion
- Ultrasound-assisted biofabrication
- Thermal non-invasive ultrasound surgery (HIFU)
- Ultrasonic methods in diagnostics and therapy of the brain
- Histotripsy – mechanical non-invasive ultrasound surgery
- Ultrasound in cosmetology
- Parametric (directed) sound emission
- Acoustic holography
- Design of high-power ultrasonic phased arrays
- Industrial airborne ultrasound
- Acoustic radiation force
Patents
- P.V. Yuldashev, V.A. Khokhlova, O.A. Sapozhnikov, P.B. Rosnitskiy, M.M. Karzova.
“A graphical shell program for calculating nonlinear ultrasonic beams generated by powerful single-element transducers and annular arrays”
Federal Service for Intellectual Property (Rospatent)
Patent No. 2017613572
22 March - V.A. Khokhlova, S.A. Ilyin, P.B. Rosnitskiy, P.V. Yuldashev, L.R. Gavrilov, O.A. Sapozhnikov.
“A program with a graphical shell designed for calculating and analyzing the quality of ultrasonic fields of multi-element phased arrays”
Federal Service for Intellectual Property (Rospatent)
Patent No. 2015618574
12 August
Contacts
Details
- in our webinar
- in our talks by:
- in the papers below
[1] “HIFU beam”: a simulator for predicting axially symmetric nonlinear acoustic fields generated by focused transducers in a layered medium / P. V. Yuldashev, M. M. Karzova, W. Kreider et al. // IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control. — 2021. — Vol. 68, no. 9. — P. 2837–2852. DOI: 10.1109/TUFFC.2021.3074611
[2] Wide-angle parabolic approximation for modeling high-intensity fields from strongly focused ultrasound transducers / P. V. Yuldashev, I. S. Mezdrokhin, V. A. Khokhlova // Acoustical Physics. — 2018. — Vol. 64, no. 3. — P. 309–319. DOI: 10.1134/S1063771018030168
[3] Simulation of three-dimensional nonlinear fields of ultrasound therapeutic arrays / P. V. Yuldashev, V. A. Khokhlova // Acoustical Physics. — 2011. — Vol. 57, no. 3. — P. 334–343.
[4] Setting boundary conditions to the Khokhlov–Zabolotskaya equation for modeling ultrasound fields generated by strongly focused transducers / P. B. Rosnitskiy, P. V. Yuldashev, B. A. Vysokanov, V. A. Khokhlova // Acoustical Physics. — 2016. — Vol. 62, no. 2. — P. 151–159. DOI: 10.1134/S1063771016020123
[5] Simulating and measuring the acoustic radiation force of a focused ultrasonic beam on elastic spheres in water / A. V. Nikolaeva, M. M. Karzova, S. A. Tsysar et al. // Bulletin of the Russian Academy of Sciences: Physics. — 2019. — Vol. 83, no. 1. — P. 77–81. DOI: 10.3103/S1062873819010192
[6] Simulation of nonlinear trans-skull focusing and formation of shocks in brain using a fully populated ultrasound array with aberration correction / P. B. Rosnitskiy, P. V. Yuldashev, O. A. Sapozhnikov et al. // Journal of the Acoustical Society of America. — 2019. — Vol. 146, no. 3. — P. 1786–1798. DOI: 10.1121/1.5126685
[7] Fully nonlinear three-dimensional modeling of parametric interactions in the field of a dual-frequency acoustic array / A. V. Kvashennikova, P. V. Yuldashev, V. A. Khokhlova, I. B. Esipov // Journal of the Acoustical Society of America. — 2024. — Vol. 155, no. 3. — P. 1682–1693. DOI: 10.1121/10.0025049
[8] Three-dimensional wide-angle parabolic equations with propagator separation based on finite fourier series / P. V. Yuldashev, E. O. Konnova, M. M. Karzova, V. A. Khokhlova // Acoustical Physics. — 2024. — Vol. 70, no. 5. — P. 783–796. DOI: 10.1134/S1063771024602206
[9] Demodulation of pulsed acoustic signals in strongly nonlinear propagation regimes / A. V. Kvashennikova, M. S. Sergeeva, P. V. Yuldashev et al. // Acoustical Physics. — 2024. — Vol. 70, no. 5. — P. 797–807. DOI: 10.1134/S1063771024602279
[10] Impact of the trajectory of treatment on the rate of thermal ablation and ablated volume of biological tissue irradiated by shockwave focused ultrasonic exposure / P. A. Pestova, P. V. Yuldashev, V. A. Khokhlova, M. M. Karzova // Bulletin of the Russian Academy of Sciences: Physics. — 2024. — Vol. 88, no. 1. — P. 108–112.
[11] Thermal ablation of biological tissue by sonicating discrete foci in a specified volume with a single wave burst with shocks / P. A. Pestova, P. V. Yuldashev, V. A. Khokhlova, M. M. Karzova // Acoustical Physics. — 2024. — Vol. 70, no. 3. — P. 434–443.
[12] The use of focused ultrasound beams with shocks to suppress diffusion effects in volumetric thermal ablation of biological tissue / P. A. Pestova, M. M. Karzova, P. V. Yuldashev, V. A. Khokhlova // Acoustical Physics. — 2023. — Vol. 69, no. 4. — P. 448–458.
[13] Comparative characterization of nonlinear ultrasound fields generated by Sonalleve V1 and V2 MR-HIFU systems / M. M. Karzova, W. Kreider, A. Partanen et al. // IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control. — 2023. — Vol. 70, no. 6. — P. 521–537. DOI: 10.1109/TUFFC.2023.3261420
[14] Compensation for aberrations when focusing ultrasound through the skull based on CT and MRI data / D. D. Chupova, P. B. Rosnitskiy, O. V. Solontsov et al. // Acoustical Physics. — 2024. — Vol. 70, no. 2. — P. 288–298.
[15] Compensation for aberrations of focused ultrasound beams in transcranial sonications of brain at different depths / D. D. Chupova, P. B. Rosnitskiy, L. R. Gavrilov, V. A. Khokhlova // Acoustical Physics. — 2022. — Vol. 68, no. 1. — P. 1–10.
[16] Nonlinear ultrasound fields generated by an annular array with electronic and geometric adjustment of its focusing angle / E. M. Ponomarchuk, P. V. Yuldashev, D. A. Nikolaev et al. // Acoustical Physics. — 2023. — Vol. 69, no. 4. — P. 459–470. DOI: 10.1134/S1063771023600560
[17] Dependence of inertial cavitation induced by high intensity focused ultrasound on transducer F-number and nonlinear waveform distortion / T. Khokhlova, P. Rosnitskiy, C. Hunter et al. // Journal of the Acoustical Society of America. — 2018. — Vol. 144, no. 3. — P. 1160–1169. DOI: 10.1121/1.5052260
[18] Performance and accuracy analysis of nonlinear k-Wave simulations using local domain decomposition with an 8-GPU server / B. Treeby, F. Vaverka, J. Jaros // Proc. Mtgs. Acoust. / 2018. — Vol. 34. — P. 022002. DOI: 10.1121/2.0000883
[19] Shock formation and nonlinear saturation effects in the ultrasound field of a diagnostic curvilinear probe / M. M. Karzova, P. V. Yuldashev, O. A. Sapozhnikov et al. // Journal of the Acoustical Society of America. — 2017. — Vol. 141, no. 4. — P. 2327–2337. DOI: 10.1121/1.4979261
[20] On the possibility of using multi-element phased arrays for shock-wave action on deep brain structures / P. Rosnitskiy, L. Gavrilov, P. Yuldashev et al. // Acoustical Physics. — 2017. — Vol. 63, no. 5. — P. 531–541. DOI: 10.1134/S1063771017050104
[21] Viscoelastic nonlinear resonator with gas-filled cavities / T. B. Krit, V. G. Andreev, I. Yu. Demin // Acta Acustica united with Acustica — 2015. — Vol. 101, no. 5. — P. 915–919. DOI: 10.3813/AAA.918886