Mathematical modelling of acoustic radiation force in transient shear wave elastography in the heart

¹ Inria, Université Paris-Saclay, Palaiseau, France
² LMS, Ecole Polytechnique, CNRS, Université Paris-Saclay, Palaiseau, France
³ Computational Cardiology Laboratory, Institute of Biophysics, Medical University of Graz, Graz, Austria

^* Corresponding author: federica.caforio@inria.fr

Received: 5 April 2019
Accepted: 17 March 2020

Abstract

The aim of this work is to provide a mathematical model and analysis of the excitation and the resulting shear wave propagation in acoustic radiation force (ARF)-based shear wave cardiac elastography. Our approach is based on asymptotic analysis; more precisely, it consists in considering a family of problems, parametrised by a small parameter inversely proportional to the excitation frequency of the probes, the viscosity and the velocity of pressure wave propagation. We derive a simplified model for the expression of the ARF by investigating the limit behaviour of the solution when the small parameter goes to zero. By formal asymptotic analysis – an asymptotic expansion of the solution is used – and energy analysis of the nonlinear elastodynamic problem, we show that the leading-order term of the expansion is solution of the underlying, incompressible, nonlinear cardiac mechanics. Subsequently, two corrector terms are derived. The first is a fast-oscillating pressure wave generated by the probes, solution of a Helmholtz equation at every time. The second corrector term consists in an elastic field with prescribed divergence, having a function of the first corrector as a source term. This field corresponds to the shear acoustic wave induced by the ARF. We also confirm that, in cardiac mechanics, the presence of viscosity in the model is essential to derive an expression of the shear wave propagation from the ARF, and that this phenomenon is related to the nonlinearity of the partial differential equation.