  Traumatic brain injury studies on the living human brain are experimentally
infeasible. We present a simulation approach that models ultrasound propagation
in the human brain while it is moving due to the complex shear shock wave
deformation from a traumatic impact. Finite difference simulations can model
ultrasound propagation in complex media such as human tissue. Recently, we have
shown that a finite difference approach can be used to represent displacements
that are much smaller than the grid size, such as the motion encountered in
shear wave propagation from ultrasound elastography. However, this
subresolution displacement model, called impedance flow, was only implemented
and validated for acoustical media composed of randomly distributed scatterers.
Here we propose a generalization of the impedance flow method that describes
the continuous subresolution motion of structured acoustical maps, and in
particular of acoustical maps of the human brain. It is shown that the average
error in the subresolution displacement method is small compared to the
wavelength. The method is applied to acoustical maps of a moving human brain,
imposed by the propagation of a shear shock wave. Then the Fullwave simulation
tool is used to model transmit-receive imaging sequences based on an L7-4
imaging transducer. The simulated radiofrequency data is beamformed using a
conventional delay-and-sum method and a normalized cross-correlation method
designed for shock wave tracking is used to determine the tissue motion. It is
shown that the proposed generalized impedance flow method accurately captures
the shear shock wave motion. This method can lead to improvements in image
sequence design that takes into account the aberration and multiple reflections
from the brain and in the design of tracking algorithms that can more
accurately capture the complex brain motion that occurs during a traumatic
impact.
