Traumatic Brain Injury (TBI) studies on the living human brain are experimentally infeasible due to ethical reasons and the elastic properties of the brain degrade rapidly post-mortem. 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 the Fullwave finite difference approach can also 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 sub-resolution 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 sub-resolution motion of structured acoustical maps, and in particular of acoustical maps of the human brain. It is shown that the average error in simulating subresolution displacements using impedance flow is small compared to the acoustical wavelength (¿/1702). The method is then applied to acoustical maps of the human brain with a motion that is imposed by the propagation of a shear shock wave. This motion is determined numerically with a custom piecewise parabolic method that is calibrated to ex vivo observations of shear shocks in the porcine brain. 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. This overall process is an in silico reproduction of the experiments that were previously performed to observe shear shock waves in fresh porcine brain. It is shown that the proposed generalized impedance flow method accurately captures the shear wave motion in terms of the wave profile, shock front characteristics, odd harmonic spectrum generation, and acceleration at the shear shock front. We expect that this approach will 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. These methods of modeling ultrasound propagation in moving media can also be applied to other displacements, such as those generated by shear wave elastography or blood flow.