Simulation of blast-induced early-time intracranial wave physics leading to traumatic brain injury

The objective of this modeling and simulation study was to establish the role of stress wave interactions in the genesis of traumatic brain injury (TBI) from exposure to explosive blast. A high resolution (1 mm3 voxels) five material model of the human head was created by segmentation of color cryosections from the Visible Human Female data set. Tissue material properties were assigned from literature values. The model was inserted into the shock physics wave code, CTH, and subjected to a simulated blast wave of 1.3 MPa (13 bars) peak pressure from anterior, posterior, and lateral directions. Three-dimensional plots of maximum pressure, volumetric tension, and deviatoric (shear) stress demonstrated significant differences related to the incident blast geometry. In particular, the calculations revealed focal brain regions of elevated pressure and deviatoric stress within the first 2 ms of blast exposure. Calculated maximum levels of 15 KPa deviatoric, 3.3 MPa pressure, and 0.8 MPa volumetric tension were observed before the onset of significant head accelerations. Over a 2 ms time course, the head model moved only 1 mm in response to the blast loading. Doubling the blast strength changed the resulting intracranial stress magnitudes but not their distribution. We conclude that stress localization, due to early-time wave interactions, may contribute to the development of multifocal axonal injury underlying TBI. We propose that a contribution to traumatic brain injury from blast exposure, and most likely blunt impact, can occur on a time scale shorter than previous model predictions and before the onset of linear or rotational accelerations traditionally associated with the development of TBI.