Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strain rate”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Further Development of the Tamped Richtmyer-Meshkov Instability Method and Application to Molybdenum Dynamic Strength Calibration and Tabulation

The high pressure and high strain rate dynamic strength of Mo is experimentally and computationally investigated in the 3–20 GPa stress, 50–600 C temperature, and 10 5 –10 6 /s strain rate regimes using a modified tamped Richtmyer-Meshkov instability (RMI) method. Modifications to the tamped RMI method include a method to determine loading states during strain, a new strength calibration function based on interface shape, and a robust uncertainty quantification method. These modifications improve fidelity of the tamped RMI method, allowing evaluation of the compensating effects of pressure hardening, strain rate hardening, strain hardening, and thermal softening. The new calibration function based on interface shape is not limited to sinusoidal corrugations and could be applied to additional interface shapes. Plate impact experiments are performed at Argonne National Laboratory’s Advanced Photon Source’s Dynamic Compression Sector operated by Washington State University (DCS), driving a planar shock front through a corrugated Mo-D 2 O or Mo-C 8 F 18 interface, forcing the corrugation to significantly deform. The extent of interfacial deformation, RMI growth, is experimentally observed using X-ray phase contrast imaging at the DCS. RMI jet lengths and jet shapes are extracted from the experimental radiographs, then used to calibrate numerical simulations performed with the Sandia National Laboratories (SNL) hydrocode CTH. Mo yield strength, Y, as a function of shock pressure, P, strain rate, $\dot{\varepsilon }$, accumulated strain, ϵ, relative volumetric compression, RD , and temperature, T , is determined for each impact experiment and presented. The calibrated Mo yield strength values range 1.2–1.8 GPa, with strength generally decreasing as the impact stress increases. This trend is likely caused by thermal softening or strain localization. The tabular yield strength versus loading condition data presented in this paper can be used to fit complex strength models.

Voorhees, T. J. [Sandia National Lab. (SNL-CA), Li↗
DOE OSTI2022

Highly stretchable ionically crosslinked acrylate elastomers inspired by polyelectrolyte complexes

Dynamic bonds are a powerful approach to tailor the mechanical properties of elastomers and introduce shape-memory, self-healing, and recyclability. Among the library of dynamic crosslinks, electrostatic interactions among oppositely charged ions have been shown to enable tough and resilient elastomers and hydrogels. In this work, we investigate the mechanical properties of ionically crosslinked ethyl acrylate-based elastomers assembled from oppositely charged copolymers. Using both infrared and Raman spectroscopy, we confirm that ionic interactions are established among polymer chains. We find that the glass transition temperature of the complex is in between the two individual copolymers, while the complex demonstrates higher stiffness and more recovery, indicating that ionic bonds can strengthen and enhance recovery of these elastomers. We compare cycles to increasing strain levels at different strain rates, and hypothesize that at fast strain rates ionic bonds dynamically break and reform while entanglements do not have time to slip, and at slow strain rates ionic interactions are disrupted and these entanglements slip significantly. Further, we show that a higher ionic to neutral monomer ratio can increase the stiffness, but its effect on recovery is minimal. Lastly, taking advantage of the versatility of acrylates, ethyl acrylate is replaced with the more hydrophilic 2-hydroxyethyl acrylate, and the latter is shown to exhibit better recovery and self-healing at a cost of stiffness and strength. The design principles uncovered for these easy-to-manufacture polyelectrolyte complex-inspired bulk materials can be broadly applied to tailor elastomer stiffness, strength, inelastic recovery, and self-healing for various applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the flow-stress model of Preston, Tonks, and Wallace - Observations, parameter constraints, and implementations (Rev. 1)

A detailed discussion about the PTW flow-stress model has been presented, including theoretical considerations, model calibrations to experimental data, and model implementations in LANL ASC hydro-codes. The main findings are: In order to reproduce the intended behavior of the PTW model at high strain rates, the following conditions on the model parameters need to be valid: y 1 ≥ s 0 and y 2 ≥ β ; When the conditions above are not satisfied, PTW and PTW_Mod1 have slightly different behavior in the phonon-drag regime, and should be treated as different models. PTW exhibits a discontinuity above the critical strain rate. PTW_Mod1 does not show such a discontinuity but implies a different strain-rate dependence at high strain rates. Examples of this behavior include older calibrations for stainless steel and depleted uranium; In the thermal-activation regime, some flexibility was introduced in the original PTW model to allow for different rate-dependent behaviors at intermediate strain rates. Example: copper. However, many material calibrations do not require this extra flexibility, and the complexity and ambiguity of the intended behavior can be partially blamed for the inconsistent model implementations across ASC codes; The limit for the material constant p going to zero is well defined and leads to a simplified flow-stress equation, discontinuity free; Modern calibrations of the PTW parameters are done using Impala. Constraints on the PTW parameters are enforced; Implementation of the PTW model in LANL ASC hydro-codes suffered of several shortcom ings. At the time of writing, FLAG and xRAGE contain consistent and correct implemen tations of the PTW and PTW_Mod1 models. Corresponding changes to Pagosa have been suggested.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. II. Kinematics

This second part of our research explores the kinematic aspect of fine-scale bubble-induced turbulence (BIT) to (i) present the effect of bubble breakup and coalescence and (ii) compare it against the universal kinematic fine-scale turbulence characteristics reported in the literature. To this end, we simulate a dilute bubbly system of 0.5% void fraction using two distinct numerical simulations. In the volume-of-fluid (VoF) simulation, bubbles undergo breakup and coalescence. In the immersed boundary method (IBM) simulation, however, they act as rigid spheres. We also perform a simulation of classical homogeneous isotropic turbulence (HIT). The first important outcome of this study is that BIT is radically different from HIT in terms of its kinematic fine-scale characteristics. In the vorticity-dominating regions, BIT exhibits a weak vortex stretching. This weak vortex stretching is due to (a) the intermediate strain-rate eigenvalues skewed weakly to positive and (b) the extensive strain-rate eigenvector aligning perpendicular to the vorticity vector. The BIT has, on average, not only a weak enstrophy production but also a weak strain production in strain-dominating regions. The weak strain production is due to (a) the presence of vortex stretching in highly strained fluid elements and (b) the absolute magnitude of compressive strain-rate eigenvalue being as close to the extensive strain-rate eigenvalue. Thus, none of the kinematic fine-scale HIT characteristics is noted for BIT. The second important conclusion is that bubble breakup and coalescence play little to no influence on the kinematics of fine-scale BIT as VoF and IBM simulations produce similar results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Dislocation density-based plasticity model from massive discrete dislocation dynamics database

In this work, we present a dislocation density-based strain hardening model for single crystal copper through a systematic coarse-graining analysis of more than 200 discrete dislocation dynamics (DDD) simulations of plastic deformation under uniaxial tension. The proposed constitutive model has two components: a generalized Taylor relation connecting resolved shear stresses to dislocation densities on individual slip systems, and a generalized Kocks-Mecking model for dislocation multiplication. The DDD data strongly suggests a logarithmic dependence of flow stress on the plastic shear strain rate on each slip system, and, equivalently, an exponential dependence of the plastic shear strain rate on the resolved shear stress. Hence the proposed generalized Taylor relation subsumes the Orowan relation for plastic flow. The DDD data also calls for a correction to the Kocks-Mecking model of dislocation multiplication to account for the increase of dislocation density on slip systems with negligible plastic shear strain rate. This is accomplished by allowing the multiplication rate on each slip system to include contributions from the plastic strain rates of the two coplanar slip systems. The resulting constitutive model successfully captures the strain hardening rate dependence on the loading orientation as predicted by the DDD simulations, which is also consistent with existing experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Pressure effects on reactivity and extinction of n -dodecane diffusion cool flame

Low-temperature combustion at high pressure has gained increasing interest owing to the growing demand for clean and efficient propulsion and energy systems. In particular, cool flame structure, extinction limits, and fuel reactivities at high pressure are critical in affecting the engine performance at near-limit conditions. Here, this paper investigates the effects of pressure on the cool flame extinction limit, structure, radical index, reactivity, and oxygen concentration dependence by experiments, analysis, and modeling. A large n-alkane, namely n-dodecane (n C 12 H 26 ), is selected to study its diffusion cool flame dynamics and reactivity in a high-pressure counterflow burner up to 10 atm. The experimental results show that higher pressure increases the cool flame extinction strain rates, and that the pressure-weighted extinction strain rate (aP) is proportional to the square of pressure, P 2 . A scaling analysis explains the relationship between the dependence of flame structure, heat release rate, and pressure-weighted cool flame strain rate on pressure. Furthermore, radical indexes at different pressures are measured by isolating the thermal and transport effects from the chemical contribution to diffusion cool flame extinction. The radical index clearly shows that the low-temperature reactivity increases with pressure. In addition, due to the critical role of multiple oxygen addition reactions in low-temperature chemistry, the relationship between the cool flame extinction limit and the oxygen concentration is explored. It is found that the cool flame extinction limits are proportional to the n th power of the oxygen concentration, [O 2 ] n , and increasing pressure leads to stronger extinction limit dependence (larger n) on the oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions, while the scaling analysis explains well the experimental results.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Towards the Ultimate Strength of Medium‐Entropy Alloys Through Pulsed Lasers

The tensile strength of metals at extreme strain rates is a key predictor of their performance in ballistic and structural impact applications. An important experimental method to reach these extreme strain rates is the use of high-amplitude, short-duration pulsed lasers. The Jupiter Laser Facility at the Lawrence Livermore National Laboratory enabled probing for the first time the mechanical response of several promising High Entropy Alloys at times on the order of nanoseconds (strain rates of ∼10 7 and ∼10 9 s −1 ). The measured strength is in the range of 6 to 10 GPa, ten times the quasistatic value. The mechanisms of plastic deformation and failure were identified and quantified through analysis and molecular dynamics simulation. The reflected wave amplitudes, obtained by VISAR, were used to determine the tensile (spall) stress. The high tensile strength obtained is due to two factors: the strain-rate dependence of plastic flow and the kinetics of void nucleation, growth, and coalescence. The experimental results are compared with an analytical prediction considering both grain-interior and grain-boundary void initiation. Molecular dynamics simulations, conducted at strain rates of 10 8 and 10 9 s −1 , rationalize the experimental results. In conclusion, they provide valuable information about the process of failure evolution, and reveal that grain boundary separation plays a pivotal role in spalling.

Materials science↗
DOE OSTI2024

Crystal plasticity modeling and analysis for the transition from intergranular to transgranular failure in nickel-based alloy Inconel 740H at elevated temperature

The precipitation-strengthened Nickel alloy Inconel® 740H® (IN740H) exhibits increased ductility at higher applied strain rates during quasi-static tensile tests at an elevated temperature of 760°C. The examination of fracture surfaces in this context reveals a noteworthy transition of underlying fracture mechanisms from transgranular to intergranular fracture as the applied strain rate decreases from 1×10 -3 / s to 0.83×10 -4 / s . To thoroughly understand the mechanical response of IN740H under these conditions, this study develops a crystal plasticity finite element (CPFE) model. Further, this model incorporates various deformation mechanisms including dislocation slips, climb, and grain boundary sliding, which are relevant to the test conditions. The model is calibrated using data from both tensile tests at different strain rates and creep tests across a broad stress range at 760°C, enabling the accurate determination of model parameters for each mechanism. Simulation results well captured the experimental observations of different failure modes. At higher strain rates, the model shows a dominance of dislocation slip leading to heterogeneous plastic deformation and formation of transgranular shear bands causing the failure, while at lower strain rates, an increased activity of grain boundary sliding causes grain boundaries crack leading to intergranular failure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Mechanochemically responsive polymer enables shockwave visualization

Abstract Understanding the physical and chemical response of materials to impulsive deformation is crucial for applications ranging from soft robotic locomotion to space exploration to seismology. However, investigating material properties at extreme strain rates remains challenging due to temporal and spatial resolution limitations. Combining high-strain-rate testing with mechanochemistry encodes the molecular-level deformation within the material itself, thus enabling the direct quantification of the material response. Here, we demonstrate a mechanophore-functionalized block copolymer that self-reports energy dissipation mechanisms, such as bond rupture and acoustic wave dissipation, in response to high-strain-rate impacts. A microprojectile accelerated towards the polymer permanently deforms the material at a shallow depth. At intersonic velocities, the polymer reports significant subsurface energy absorption due to shockwave attenuation, a mechanism traditionally considered negligible compared to plasticity and not well explored in polymers. The acoustic wave velocity of the material is directly recovered from the mechanochemically-activated subsurface volume recorded in the material, which is validated by simulations, theory, and acoustic measurements. This integration of mechanochemistry with microballistic testing enables characterization of high-strain-rate mechanical properties and elucidates important insights applicable to nanomaterials, particle-reinforced composites, and biocompatible polymers.

Science & Technology - Other Topics↗
DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Kinetics and extinction of non-premixed cool and warm flames of dimethyl ether at elevated pressure

Here, the growing demand of clean and efficient propulsion and energy systems has sparked an interest in understanding low-temperature combustion at high pressure. Cool flame transition and extinction limits as well as oxygen concentration dependence at elevated pressures provide insights of the low-temperature and high-pressure fuel reactivity. A new experimental high-pressure counterflow burner platform was designed and developed to achieve the studies of high-pressure cool flames. Dimethyl ether (DME) was chosen to study its non-premixed cool flame in high-pressure counterflow burner at pressure up to 5 atm, perhaps for the first time. This paper investigates the effects of pressure on cool flame structure, extinction and transition limits, and oxygen concentration dependence as well as ozone assisted warm flames of DME in experiments and numerical simulations. The results show that the reignition transition from cool flame to hot flame occurs either with the decrease of the strain rate at a given fuel concentration and pressure or with the increase of fuel mole fraction or pressure at a given strain rate. Furthermore, it is shown that the higher pressure shifts the cool flame to higher strain rates and results in higher cool flame extinction strain rates. However, the existing kinetic model of DME fails in predicting the cool flame extinction limit at elevated pressures. Besides, the cool flame extinction limits are proportional to nth power of the oxygen concentration, [O 2 ] n , and the increase of pressure leads to stronger extinction limit dependence (larger n) on oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions. In addition, stable warm flame was first experimentally observed by using DME at elevated pressure with ozone sensitization.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

What Can Plastic Flow Fields Tell Us About Heat Sources in Deformation Processing?

We characterize primary (shape-change) and secondary (friction) deformation, and associated temperature fields, in metal cutting and forming processes, using in situ imaging and simulation. The experimental configurations enable access to the deformation zones and die contact-interfaces, for measuring deformation, temperature, and frictional drag. Infra-red thermography reveals that the plastic strain-rate field is an excellent proxy for the deformation-induced heat sources. Both spatially confined and diffuse strain-rate fields occur, depending on the initial workpiece deformation state. When the strain rate is confined, as in pre-hardened material, the temperature modeling is much simplified, as the heat source is also now localized. However complex, microstructure-driven deformation modes, like sinuous flow in annealed metals, result in spatially diffuse strain-rate and body heat sources, more challenging to analyze. Furthermore, our unified measurements should be of value for accurately estimating the fraction of plastic dissipation that is converted into heat in large-strain deformation processes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Impala's horn applied to posterior samples of Ti-6Al-4V strength model parameters

We have generated strength model parameters for the Preston-Tonks-Wallace (PTW) model of plastic deformation based on a collection of Ti-6Al-4V alloy data and a Bayesian framework. This data collected from the literature includes different chemistries and processes. The experiments include quasistatic compression, Split Hopkinson Pressure Bar (SHPB) compression, and Taylor cylinders. Here we use posterior parameter distributions applied to a hypothetical impact referred to as the Impala's horn. The Impala's horn is like a Taylor cylinder, but modified to have a conical pro le that is truncated before coming to a point. The smaller diameter at the impact end results in higher calculated strain rates and strains than the same velocity impact with the traditional Taylor profile. The modified geometry exercises the behavior of the strength model outside the range of conditions found in the calibration data. Here we sort posterior strength models based on their flow stress at characteristic values for this hypothetical deformation. Then we compare the performance of these posterior models based on two simple observable final state properties of the Impala's horn.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Bead–Spring Simulation of Ionomer Melts—Studying the Effects of Chain-Length and Associating Group Fraction on Equilibrium Structure and Extensional Flow Behavior

Ionomers are associative polymers with diverse applications ranging from selective membranes and high-performance adhesives to abrasion- and chemical-resistant coatings, insulation layers, vacuum packaging, and foamed sheets. Within equilibrium melt, the ionic or associating groups are known to form thermally reversible, associative clusters whose presence can significantly affect the system’s mechanical, viscoelastic, and transport properties. It is, thus, of great interest to understand how to control such clusters’ size distribution, shape, and stability through the designed choice of polymer architecture and the ionic groups’ fraction, arrangement, and interaction strength. In this work, we represent linear associating polymers using a Kremer–Grest type bead–spring model and perform large-scale MD simulations to explore the effect of polymer chain-length (l) and fraction (f s ) of randomly placed associating groups on the size distribution and stability of formed clusters. We consider different chain-lengths (below and above entanglement), varying fractions of associating groups (represented by ‘sticky’ beads) between 5 and 20%, and a fixed sticky–sticky nonbond interaction strength of four times that between regular non-associating beads. For all melts containing associating groups the equilibrium structure factor S(q) displays a signature ionomer peak at low wave vector q whose intensity increases with increasing f s and l. The average cluster size Nc increases with f s . However, the effect of chain-length on Nc appears to be pronounced only at higher values of f s . Under extensional flows, the computed stress (and viscosity) is higher at higher f s and l regardless of strain rate. Beyond a critical strain rate, we observe fragmentation of the associative clusters, which has interesting effects on the stress/viscous response.

36 MATERIALS SCIENCE↗
DOE OSTI2022

In-Situ Observation of FCC→HCP Transformation-Induced Plasticity Behavior During Dynamic Deformation of CoCrNi Multi-principal Element Alloys

CoCrNi multi-principal element alloys (MPEAs) have potential as high toughness, blast resistant materials, but limited studies of their high strain rate behavior have been performed. Their high toughness has been attributed to low stacking fault energies and the activation of toughness enhancing deformation mechanisms, such as twinning- and transformation-induced plasticity (TWIP and TRIP, respectively), that may be affected by increased strain rate. This work examines a series of CoCrNi MPEAs tested with a miniaturized pressure bar, while simultaneous synchrotron x-ray diffraction and imaging are performed, and compares these results to those collected during quasi-static deformation at elevated temperatures. Here, the collected diffraction indicates that TRIP was active at elevated strain rates for the most Co-rich alloys containing 50 or 55 at. pct. Co, while TRIP behavior was suppressed at elevated strain rates in a 40 at. pct. Co alloy anticipated to experience TRIP. The suppression of TRIP in the lower Co alloy at high strain rates is attributed to adiabatic heating, consistent with suppression of TRIP as a result of elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental study of vorticity-strain interactions in turbulent premixed counterflow flames

Here, the effects of heat release on interactions between vorticity (ω) and strain rate (s) in turbulent premixed CH 4 /O 2 /N 2 counterflow flames are investigated using simultaneous OH laser-induced fluorescence (LIF) and tomographic particle image velocimetry (TPIV) measurements. A comparison between the flames and a corresponding turbulent non-reacting variable density N 2 -vs-products counterflow reveals the impact of heat release on vorticity-strain rate alignment statistics. Vorticity and strain rate statistics in the flames and non-reacting flow are conditioned on distance from the local flame front and gas mixing layer interface (GMLI) contours, respectively. The magnitude, alignment, and spatial distribution of the vorticity and principal strain rates (s 1 , s 2 , s 3 ) are rather different when heat release is present. Density variations without heat release enhance the ω-s 2 alignment while significantly reducing the ω-s 3 alignment and modestly reducing the ω-s 1 alignment. In contrast, heat release at the flame front further reduces the ω-s 1 alignment but increases the ω-s 3 alignment and suppresses the preferential ω-s 2 alignment. Furthermore, increasing turbulence diminishes the effect of heat release on this preferential alignment. In regions with the largest vorticities, both the reacting and non-reacting counterflows show an increase in the probability of ω-s 2 alignment. All counterflow cases have a net positive vortex-stretching contribution to the enstrophy production with a peak production rate at the flame front or GMLI, but the peak values depend on the density variation, heat release, and turbulence level. Elucidation of the complex interplay between these factors contributes to the understanding of the dynamics of turbulence-flame interactions.

42 ENGINEERING↗
DOE OSTI2021

Localizing Clinical Patterns of Blast Traumatic Brain Injury Through Computational Modeling and Simulation

Blast traumatic brain injury is ubiquitous in modern military conflict with significant morbidity and mortality. Yet the mechanism by which blast overpressure waves cause specific intracranial injury in humans remains unclear. Reviewing of both the clinical experience of neurointensivists and neurosurgeons who treated service members exposed to blast have revealed a pattern of injury to cerebral blood vessels, manifested as subarachnoid hemorrhage, pseudoaneurysm, and early diffuse cerebral edema. Additionally, a seminal neuropathologic case series of victims of blast traumatic brain injury (TBI) showed unique astroglial scarring patterns at the following tissue interfaces: subpial glial plate, perivascular, periventricular, and cerebral gray-white interface. The uniting feature of both the clinical and neuropathologic findings in blast TBI is the co-location of injury to material interfaces, be it solid-fluid or solid-solid interface. This motivates the hypothesis that blast TBI is an injury at the intracranial mechanical interfaces. In order to investigate the intracranial interface dynamics, we performed a novel set of computational simulations using a model human head simplified but containing models of gyri, sulci, cerebrospinal fluid (CSF), ventricles, and vasculature with high spatial resolution of the mechanical interfaces. Simulations were performed within a hybrid Eulerian—Lagrangian simulation suite (CTH coupled via Zapotec to Sierra Mechanics). Because of the large computational meshes, simulations required high performance computing resources. Twenty simulations were performed across multiple exposure scenarios—overpressures of 150, 250, and 500 kPa with 1 ms overpressure durations—for multiple blast exposures (front blast, side blast, and wall blast) across large variations in material model parameters (brain shear properties, skull elastic moduli). All simulations predict fluid cavitation within CSF (where intracerebral vasculature reside) with cavitation occurring deep and diffusely into cerebral sulci. These cavitation events are adjacent to high interface strain rates at the subpial glial plate. Larger overpressure simulations (250 and 500kPa) demonstrated intraventricular cavitation—also associated with adjacent high periventricular strain rates. Additionally, models of embedded intraparenchymal vascular structures—with diameters as small as 0.6 mm—predicted intravascular cavitation with adjacent high perivascular strain rates. The co-location of local maxima of strain rates near several of the regions that appear to be preferentially damaged in blast TBI (vascular structures, subpial glial plate, perivascular regions, and periventricular regions) suggest that intracranial interface dynamics may be important in understanding how blast overpressures leads to intracranial injury.

59 BASIC BIOLOGICAL SCIENCES↗