Microcavitation: the key to modeling blast traumatic brain injury?

Estimating Brain Injury Risk from Shipborne Underwater Blasts Using a High-fidelity Finite Element Head Model

INTRODUCTION

Assessing the survivability of, and potential injury to, a ship's crew from underwater blast is crucial to understanding the operating capability of a military vessel following blast exposure. One form of injury that can occur and affect a crew member's ability to perform tasks is traumatic brain injury (TBI). To evaluate the risk of TBI from underwater blasts, injury metrics based on linear head acceleration have traditionally been used. Although these metrics are popular given their ease of use, they do not provide a direct measure of the tissue-level biomechanical responses that have been shown to cause neuronal injury. Tissue-based metrics of injury, on the other hand, may provide more insight into the potential risk of brain injury. Therefore, in this study, we assess the risk of TBI from underwater blasts using tissue-based measures of injury, such as tissue strain, strain rate, and intracranial pressure, in addition to the more commonly used head acceleration-based injury metrics.

MATERIALS AND METHODS

A series of computational simulations were performed using a detailed finite element (FE) head model to study how inertial loading of the head from underwater blast events translates to potential injury in the brain. The head kinematics loading conditions for the simulations were obtained directly from Floating Shock Platform (FSP) tests where 3 Anthropomorphic Test Devices (ATDs) were positioned at 3 shipboard locations (desk, bulkhead, and bench), and the head acceleration was directly measured. The effect of the position and orientation of the ATDs and the distance of the underwater blast from the FSP (20-50 ft) on the risk of brain injury were assessed from the FE analysis.

RESULTS

The head accelerations and estimated TBI risk from the underwater blasts highly depend on the positioning of the ATDs on the FSP and decrease in severity as the charge standoff distance is increased. The ATD that was seated at a desk had the largest peak linear head acceleration (77.5 g) and negative intracranial pressure (-51.8 kPa). In contrast, the ATD that was standing at a bulkhead had the largest computed 95th percentile maximum principal strain (19%) and strain rate (25 s-1) in the brain. For all tested conditions, none of the ATDs exceeded the Head Injury Criterion (HIC-15) threshold of 700 for serious or fatal brain injury; however, the predicted tissue strains of the bulkhead ATD at the 20-ft charge standoff distance were within the range of proposed strain thresholds for a 50% risk of concussive injury, which illustrates the added value of considering tissue-level measures in addition to head acceleration when evaluating brain injury risk.

CONCLUSIONS

In this work, we assessed the risk of brain injury from underwater blasts using an anatomically detailed subject-specific FE head model. Accurate assessment of the risk of TBI from underwater explosions is important to evaluate the potential injury risk to crew members from underwater blast events, and to guide the development of future injury mitigation strategies to maintain the safety of crew members on military ships.

Converging-diverging shock-driven instabilities along soft hydrogel surfaces

Intense surface eruptions are observed along the curved surface of a confined cylindrical film of hydrogel subject to laser-induced converging-diverging shock loading. Detailed numerical simulations are used to identify the dominant mechanisms causing mechanical instability. The mechanisms that produce surface instability are found to be fundamentally different from both acoustic parametric instability and shock-driven Richtmyer-Meshkov instability. The time scale of observed and simulated eruption formation is much larger than that of a single shock reflection, in stark contrast to previously studied shock-driven instabilities. Moreover, surface undulations are only found along external, as opposed to internal, soft solid boundaries. Specifically, classic bubble surface instability mechanisms do not occur in our experiments and here we comment only on the new surface undulations found along the outer boundary of solid hydrogel cylinders. Our findings indicate a new class of impulsively excited surface instability that is driven by cycles of internal shock reflections.

Non-Lethal Blasts can Generate Cavitation in Cerebrospinal Fluid While Severe Helmeted Impacts Cannot: A Novel Mechanism for Blast Brain Injury

Cerebrospinal fluid (CSF) cavitation is a likely physical mechanism for producing traumatic brain injury (TBI) under mechanical loading. In this study, we investigated CSF cavitation under blasts and helmeted impacts which represented loadings in battlefield and road traffic/sports collisions. We first predicted the human head response under the blasts and impacts using computational modelling and found that the blasts can produce much lower negative pressure at the contrecoup CSF region than the impacts. Further analysis showed that the pressure waves transmitting through the skull and soft tissue are responsible for producing the negative pressure at the contrecoup region. Based on this mechanism, we hypothesised that blast, and not impact, can produce CSF cavitation. To test this hypothesis, we developed a one-dimensional simplified surrogate model of the head and exposed it to both blasts and impacts. The test results confirmed the hypothesis and computational modelling of the tests validated the proposed mechanism. These findings have important implications for prevention and diagnosis of blast TBI.

Shock-Induced Damage Mechanism of Perineuronal Nets

The perineuronal net (PNN) region of the brain’s extracellular matrix (ECM) surrounds the neural networks within the brain tissue. The PNN is a protective net-like structure regulating neuronal activity such as neurotransmission, charge balance, and action potential generation. Shock-induced damage of this essential component may lead to neuronal cell death and neurodegenerations. The shock generated during a vehicle accident, fall, or improvised device explosion may produce sufficient energy to damage the structure of the PNN. The goal is to investigate the mechanics of the PNN in reaction to shock loading and to understand the mechanical properties of different PNN components such as glycan, GAG, and protein. In this study, we evaluated the mechanical strength of PNN molecules and the interfacial strength between the PNN components. Afterward, we assessed the PNN molecules’ damage efficiency under various conditions such as shock speed, preexisting bubble, and boundary conditions. The secondary structure altercation of the protein molecules of the PNN was analyzed to evaluate damage intensity under varying shock speeds. At a higher shock speed, damage intensity is more elevated, and hyaluronan (glycan molecule) is most likely to break at the rigid junction. The primary structure of the protein molecules is least likely to fail. Instead, the molecules’ secondary bonds will be altered. Our study suggests that the number of hydrogen bonds during the shock wave propagation is reduced, which leads to the change in protein conformations and damage within the PNN structure. As such, we found a direct connection between shock wave intensity and PNN damage.

Computational modeling investigation of pulsed high peak power microwaves and the potential for traumatic brain injury

When considering safety standards for human exposure to radiofrequency (RF) and microwave energy, the dominant concerns pertain to a thermal effect. However, in the case of high-power pulsed RF/microwave energy, a rapid thermal expansion can lead to stress waves within the body. In this study, a computational model is used to estimate the temperature profile in the human brain resulting from exposure to various RF/microwave incident field parameters. The temperatures are subsequently used to simulate the resulting mechanical response of the brain. Our simulations show that, for certain extremely high-power microwave exposures (permissible by current safety standards), very high stresses may occur within the brain that may have implications for neuropathological effects. Although the required power densities are orders of magnitude larger than most real-world exposure conditions, they can be achieved with devices meant to emit high-power electromagnetic pulses in military and research applications.

Evolution of the Laser-Induced Spallation Technique in Film Adhesion Measurement

Laser-induced spallation is a process in which a stress wave generated from a rapid, high-energy laser pulse initiates the ejection of surface material opposite the surface of laser impingement. Through knowledge of the stress-wave amplitude that causes film separation, the adhesion and interfacial properties of a film-on-substrate system are determined. Some advantages of the laser spallation technique are the noncontact loading, development of large stresses (on the order of GPa), and high strain rates, up to 10⁸/s. The applicability to both relatively thick films, tens of microns, and thin films, tens of nm, make it a unique technique for a wide range of materials and applications. This review combines the available knowledge and experience in laser spallation, as a state-of-the-art measurement tool, in a comprehensive pedagogical publication for the first time. An historical review of adhesion measurement by the laser-induced spallation technique, from its inception in the 1970s through the present day, is provided. An overview of the technique together with the physics governing the laser-induced spallation process, including functions of the absorbing and confining materials, are also discussed. Special attention is given to applications of laser spallation as an adhesion quantification technique in metals, polymers, composites, ceramics, and biological films. A compendium of available experimental parameters is provided that summarizes key laser spallation experiments across these thin-film materials. This review concludes with a future outlook for the laser spallation technique, which approaches its semicentennial anniversary.

Shape stability of a gas bubble in a soft solid

Predicting the onset of non-spherical oscillations of bubbles in soft matter is a fundamental cavitation problem with implications to sonoprocessing, polymeric materials synthesis, and biomedical ultrasound applications. The shape stability of a bubble in a Kelvin-Voigt viscoelastic medium with nonlinear elasticity, the simplest constitutive model for soft solids, is analytically investigated and compared to experiments. Using perturbation methods, we develop a model reducing the equations of motion to two sets of evolution equations: a Rayleigh-Plesset-type equation for the mean (volume-equivalent) bubble radius and an equation for the non-spherical mode amplitudes. Parametric instability is predicted by examining the natural frequency and the Mathieu equation for the non-spherical modes, which are obtained from our model. Our theoretical results show good agreement with published experiments of the shape oscillations of a bubble in a gelatin gel. We further examine the impact of viscoelasticity on the time evolution of non-spherical mode amplitudes. In particular, we find that viscosity increases the damping rate, thus suppressing the shape instability, while shear modulus increases the natural frequency, which changes the unstable mode. We also explain the contributions of rotational and irrotational fields to the viscoelastic stresses in the surroundings and at the bubble surface, as these contributions affect the damping rate and the unstable mode. Our analysis on the role of viscoelasticity is potentially useful to measure viscoelastic properties of soft materials by experimentally observing the shape oscillations of a bubble.

Mechanisms of cell damage due to mechanical impact: an in vitro investigation

The dynamic response of cells when subjected to mechanical impact has become increasingly relevant for accurate assessment of potential blunt injuries and elucidating underlying injury mechanisms. When exposed to mechanical impact, a biological system such as the human skin, brain, or liver is rapidly accelerated, which could result in blunt injuries. For this reason, an acceleration of greater than > 150 g is the most commonly used criteria for head injury. To understand the main mechanism(s) of blunt injury under such extreme dynamic threats, we have developed an innovative experimental method that applies a well-characterized and -controlled mechanical impact to live cells cultured in a custom-built in vitro setup compatible with live cell microscopy. Our studies using fibroblast cells as a model indicate that input acceleration (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${a}_{in}$$\end{document}ain) alone, even when it is much greater than the typical injury criteria, e.g., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${a}_{in}>1{,}000$$\end{document}ain>1,000 g, does not result in cell damage. On the contrary, we have observed a material-dependent critical pressure value above which a sudden decrease in cell population and cell membrane damage have been observed. We have unambiguously shown that (1) this critical pressure is associated with the onset of cavitation bubbles in a cell culture chamber and (2) the dynamics of cavitation bubbles in the chamber induces localized compressive/tensile pressure cycles, with an amplitude that is considerably greater than the acceleration-induced pressure, to cells. More importantly, the rate of pressure change with time for cavitation-induced pressure is significantly faster (more than ten times) than acceleration-induced pressure. Our in vitro study on the dynamic response of biological systems due to mechanical impact is a crucial step towards understanding potential mechanism(s) of blunt injury and implementing novel therapeutic strategies post-trauma.

Characterization of Astrocytic Response after Experiencing Cavitation In Vitro

When a traumatic brain injury (TBI) occurs, low-pressure regions inside the skull can cause vapor contents in the cerebral spinal fluid (CSF) to expand and collapse, a phenomenon known as cavitation. When these microbubbles (MBs) collapse, shock waves are radiated outward and are known to damage surrounding materials in other applications, like the steel foundation of boat propellers, so it is alarming to realize the damage that cavitation inflicts on vulnerable brain tissue. Using cell-laden microfibers, the longitudinal morphological response that mouse astrocytes have to surrounding cavitation in vitro is visually analyzed. Astrocytic damage is evident immediately after cavitation when compared to a control sample, as their processes retract. Forty-eight hours later, the astrocytes appeared to spread across the fibers, as normal. This study also analyzes the gene expression changes that occur post-cavitation via quantitative polymerase chain reaction (qPCR) methods. After cavitation a number of pro-inflammatory genes are upregulated, including TNFα, IL-1β, C1q, Serping1, NOS1, IL-6, and JMJD3. Taken together, these results confirm that surrounding cavitation is detrimental to astrocytic function, and yield opportunities to further the understanding of how protective headgear can minimize or eliminate the occurrence of cavitation.

Modeling tissue-selective cavitation damage

The destructive growth and collapse of cavitation bubbles are used for therapeutic purposes in focused ultrasound procedures and can contribute to tissue damage in traumatic injuries. Histotripsy is a focused ultrasound procedure that relies on controlled cavitation to homogenize soft tissue. Experimental studies of histotripsy cavitation have shown that the extent of ablation in different tissues depends on tissue mechanical properties and waveform parameters. Variable tissue susceptibility to the large stresses, strains, and strain rates developed by cavitation bubbles has been suggested as a basis for localized liver tumor treatments that spare large vessels and bile ducts. However, field quantities developed within microns of cavitation bubbles are too localized and transient to measure in experiments. Previous numerical studies have attempted to circumvent this challenge but made limited use of realistic tissue property data. In this study, numerical simulations are used to calculate stress, strain, and strain rate fields produced by bubble oscillation under histotripsy forcing in a variety of tissues with literature-sourced viscoelastic and acoustic properties. Strain field calculations are then used to predict a theoretical damage radius using tissue ultimate strain data. Simulation results support the hypothesis that differential tissue responses could be used to design tissue-selective treatments. Results agree with studies correlating tissue ultimate fractional strain with resistance to histotripsy ablation and are also consistent with experiments demonstrating smaller lesion size under exposure to higher frequency waveforms. Methods presented in this study provide an approach for modeling tissue-selective cavitation damage in general.

Multi-Scale Modeling of Head Kinematics and Brain Tissue Response to Blast Exposure

Injuries resulting from blast exposure have been increasingly prevalent in recent conflicts, with a particular focus on the risk of head injury. In the current study, a multibody model (GEBOD) was used to investigate the gross kinematics resulting from blast exposure, including longer duration events such as the fall and ground impact. Additionally, detailed planar head models, in the sagittal and transverse planes, were used to model the primary blast wave interaction with the head, and resulting tissue response. For severe blast load cases (scaled distance less than 2), the translational head accelerations during primary blast were found to increase as the height-of-burst (HOB) was lowered, while the HOB was found to have no effect for cases with scaled distance greater than 2. The HOB was found to affect both the magnitude and direction of rotational accelerations, with increasing magnitudes as the HOB deviated from the height of the head. The choice of ground contact stiffness was found to greatly affect the predicted head accelerations during ground impact. For a medium soil ground material, the kinematics during ground impact were greater for scaled distances exceeding 1.5, below which the primary blast produced greater kinematic head response.

Damage and Failure of Axonal Microtubule under Extreme High Strain Rate: An In-Silico Molecular Dynamics Study

As a major cytoskeleton element of the axon, the breaking of microtubules (MTs) has been considered as a major cause of the axon degeneration. High strain rate loading is considered as one of the key factors in microtubule breaking. Due to the small size of microtubule, the real-time behavior of microtubule breaking is hard to capture. This study employs fully-atomistic molecular dynamics (MD) simulation to determine the failure modes of microtubule under different loadings conditions such as, unidirectional stretching, bending and hydrostatic expansion. For each loading conditions, MT is subjected to extreme high strain rate (108-109 s-1) loading. We argue that such level of high strain rate may be realized during cavitation bubble implosion. For each loading type, we have determined the critical energy for MT rupture. The associated rupture mechanisms are also discussed. We observed that the stretching has the lowest energy barrier to break the MT at the nanosecond time scale. Moreover, the breakage between the dimers starts at ~16% of total strain when stretched, which is much smaller compared to the reported strain-at-failure (50%) for lower strain rate loading. It suggests that MT fails at a significantly smaller strain states when loaded at higher strain rates.