Accelerating numerical methods for nonlinear acoustics using nested meshes
share
The numerical simulation of nonlinear ultrasound is important in the treatment planning for high-intensity focused ultrasound (HIFU) therapies in the abdomen. However, the large domain sizes and generation of higher harmonics at the focus make these problems extremely computationally demanding. Numerical methods typically employ a uniform mesh fine enough to resolve the highest harmonic present in the problem, leading to a very large number of degrees of freedom. This paper proposes a more efficient strategy in which each harmonic is approximated on a separate mesh, the size of which is proportional to wavelength of the harmonic. The increase in resolution required to resolve a smaller wavelength is balanced by a reduction in the domain size. This nested meshing is feasible owing to the increasingly localised nature of higher harmonics near the focus. Numerical experiments are performed for HIFU transducers in homogeneous media in order to determine the size of the separate meshes required to accurately represent the harmonics. In particular, a fast volume potential (VP) approach is proposed and employed to perform convergence experiments as the computation domain size is modified. The VP approach allows each harmonic to be computed via the evaluation of an integral over the domain. Discretising this integral using the midpoint rule allows the computations to be performed rapidly with the FFT. It is shown that at least an order of magnitude reduction in memory consumption and computation time can be achieved with nested meshing.
READ FULL TEXT