Abstract: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 postmortem. 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 subresolution displacement model, called impedance flow, was only implemented and validated for acoustical media composed of randomly distributed scatterers. Herein, 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 simulating subresolution displacements using impedance flow is small when compared to the acoustical wavelength ( $lambda $ /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 radio frequency data are 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.

Full-waveform ultrasound modeling of soft tissue-bone interactions using conforming hexahedral meshes

Image-based Meshing of Patient-Specific Data:Converting Medical Scans into Highly Accurate Computational Models

An Investigation of Two Finite Element Modeling Solutions for Biomechanical Simulation Using a Case Study of a Mandibular Bone.

Design and simulation of a 2-D array flexible ultrasound transducer system for tissue characterization – A pilot study

Nonlinear Medical Ultrasound Tomography: 3D Modeling of Sound Wave Propagation in Human Tissues

Biomechanical modelling of probe to tissue interaction during ultrasound scanning

Assessment and Improvement of a Novel ultrasound-based 3D Reconstruction Method: Registered for Lumbar Spine

Frequency-domain full-waveform inversion-based musculoskeletal ultrasound computed tomography

A forward model incorporating elevation-focused transducer properties for 3D full-waveform inversion in ultrasound computed tomography

Hexahedral mesh of anatomical atlas for construction of computational human brain models: Applications to modeling biomechanics and bioelectric field propagation

A fast full-wave solver for trans-abdominal high-intensity focused ultrasound treatment planning

A comparison study assessing the feasibility of ultrasound-initialized deformable bone models

A Forward Model Incorporating Elevation-Focused Transducer Properties for 3-D Full-Waveform Inversion in Ultrasound Computed Tomography

Model-based Inversion of Ultrasonic Guided Waves for Cortical Bone Properties

Modeling of Personalized Anatomy using Plastic Strains

Ultrasound imaging of the human body with three dimensional full-wave nonlinear acoustics. Part 1: simulations methods

High-resolution bone microstructure imaging based on ultrasonic frequency-domain full-waveform inversion

Modeling Ultrasound Propagation in the Moving Brain: Applications to Shear Shock Waves and Traumatic Brain Injury

Biomechanics applied to computer-aided diagnosis: examples of orbital and maxillofacial surgeries

In-vivo measurement of the human soft tissues constitutive laws. Applications to Computer Aided Surgery

Full-Waveform Inversion Imaging of Cortical Bone Using Phased Array Tomography