SEARCH · Search NASA

Results for “elastic shear waves”

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.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS2000

Use of Guided Acoustic Waves to Assess the Effects of Thermal-Mechanical Cycling on Composite Stiffness

The introduction of new, advanced composite materials into aviation systems requires it thorough understanding of the long-term effects of combined thermal and mechanical loading. As part of a study to evaluate the effects of thermal-mechanical cycling, it guided acoustic (Lamb) wave measurement system was used to measure the bending and out-of-plane stiffness coefficients of composite laminates undergoing thermal-mechanical loading. The system uses a pulse/receive technique that excites an antisymmetric Lamb mode and measures the time-of-flight over a wide frequency range. Given the material density and plate thickness, the bending and out-of-plane shear stiffnesses are calculated from a reconstruction of the velocity dispersion curve. A series of 16 and 32-ply composite laminates were subjected to it thermal-mechanical loading profile in load frames equipped with special environmental chambers. The composite systems studied were it graphite fiber reinforced amorphous thermoplastic polyimide and it graphite fiber reinforced bismaleimide thermoset. The samples were exposed to both high and low temperature extremes its well as high and low strain profiles. The bending and out-of-plane stiffnesses for composite sample that have undergone over 6,000 cycles of thermal-mechanical loading are reported. The Lamb wave generated elastic stiffness results have shown decreases of up to 20% at 4,936 loading cycles for the graphite/thermoplastic samples and up to 64% at 4,706 loading cycles for the graphite/thermoset samples.

Seale, Michael D.↗

DOE OSTI2024

The role of annealing and grain boundary controls on the mechanical properties of limestones and marbles

Chemical and mechanical processes are coupled in many geological and geochemical environments. Reactive processes anneal defects and restructure grain boundaries, modifying their elastic properties, levels of internal friction, wave propagation rates and fracture behaviors. The nature of these changes is, however, contingent on the initial state of the rock. Here, in this study, impulse excitation was used to measure changes in mechanical properties as a function of dry and steam heating time at 300 °C for three carbonate rocks: Carrara marble, Carthage marble (Burlington Limestone), and Texas Cream limestone (Edwards limestone), with initial porosities of about 1 %, 3 %, and 27 %, respectively. Frequency-dependent phenomena along with mineral recrystallization were observed. This was coupled with small-angle X-ray and neutron scattering analysis to determine the relationship between changes in bulk mechanical and microstructural properties. An observed decrease in both the Young's and shear moduli in the experimental limestones as a function of heating time reflects and quantifies a reduction in the stiffness of the rock due to annealing. The internal friction of the samples first increases then decreases with reaction time, reflecting defect annealing, but this was balanced against grain boundary dissolution apparently driven by condensation of steam in the confined grain-boundary environment. This suggests an increase in the boiling point under confinement leading to dissolution, increased porosity and widening of the grain boundaries. In addition, comparison of small-angle X-ray and neutron scattering results suggests that it is inappropriate to assume that pores are empty for quantitative analysis of typical small-angle scattering samples.

58 GEOSCIENCES↗

NASA NTRS1999

Ultrasonic Nondestructive Characterization of Adhesive Bonds

Adhesives and adhesive joints are widely used in various industrial applications to reduce weight and costs, and to increase reliability. For example, advances in aerospace technology have been made possible, in part, through the use of lightweight materials and weight-saving structural designs. Joints, in particular, have been and continue to be areas in which weight can be trimmed from an airframe through the use of novel attachment techniques. In order to save weight over traditional riveted designs, to avoid the introduction of stress concentrations associated with rivet holes, and to take full advantage of advanced composite materials, engineers and designers have been specifying an ever-increasing number of adhesively bonded joints for use on airframes. Nondestructive characterization for quality control and remaining life prediction has been a key enabling technology for the effective use of adhesive joints. Conventional linear ultrasonic techniques generally can only detect flaws (delamination, cracks, voids, etc) in the joint assembly. However, more important to structural reliability is the bond strength. Although strength, in principle, cannot be measured nondestructively, a slight change in material nonlinearity may indicate the onset of failure. Furthermore, microstructural variations due to aging or under-curing may also cause changes in the third order elastic constants, which are related to the ultrasonic nonlinear parameter of the polymer adhesive. It is therefore reasonable to anticipate a correlation between changes in the ultrasonic nonlinear acoustic parameter and the remaining bond strength. It has been observed that higher harmonics of the fundamental frequency are generated when an ultrasonic wave passes through a nonlinear material. It seems that such nonlinearity can be effectively used to characterize bond strength. Several theories have been developed to model this nonlinear effect. Based on a microscopic description of the nonlinear interface binding force, a quantitative method was presented. Recently, a comparison between the experimental and simulated results based on a similar theoretical model was presented. A through-transmission setup for water immersion mode-converted shear waves was used to analyze the ultrasonic nonlinear parameter of an adhesive bond. In addition, ultrasonic guided waves have been used to analyze adhesive or diffusion bonded joints. In this paper, the ultrasonic nonlinear parameter is used to characterize the curing state of a polymer/aluminum adhesive joint. Ultrasonic through-transmission tests were conducted on samples cured under various conditions. The magnitude of the second order harmonic was measured and the corresponding ultrasonic nonlinear parameter was evaluated. A fairly good correlation between the curing condition and the nonlinear parameter is observed. The results show that the nonlinear parameter might be used as a good indicator of the cure state for adhesive joints.

Qu, Jianmin↗

DOE OSTI2021

Seismic Visibility of Eclogite in the Earth’s Upper Mantle—Implications From High Pressure‐Temperature Single‐Crystal Elastic Properties of Omphacite

Abstract Identifying and locating the geochemical and geophysical heterogeneities in the Earth’s interior is one of the most important and challenging tasks for the deep Earth scientists. Subducted oceanic crust metamorphizes into the dense eclogite in the upper mantle and is considered as a major cause of geochemical and geophysical heterogeneities in the deep Earth. In order to detect eclogitic materials inside the Earth, precise measurements of the high pressure‐temperature single‐crystal elasticity of major minerals in eclogite are thus exceedingly important. Omphacite, a Na,Al‐bearing clinopyroxene, constitutes up to 75 vol% of eclogite. In the present study, we performed the first high pressure‐temperature single‐crystal elasticity measurements of omphacite using Brillouin spectroscopy. Utilizing the finite‐strain approach, we obtained the following thermoelastic parameters for omphacite: K S0 ’ = 4.5(1), G 0 ’ = 1.53(5), ∂K S0 /∂ T = −0.029(5) GPa/K, ∂G 0 /∂ T = −0.013(5) GPa/K, with K S0 = 123(3) GPa, G 0 = 74(2) GPa, and ρ 0 = 3.34(1) g/cm 3 . We found that the seismic velocities of undeformed eclogite are similar to pyrolite at the depths of 200–300 and 410–500 km, thus eclogite is seismically invisible at these depths. Combined with the lattice‐preferred orientations of the omphacite in naturally deformed eclogites, we also modeled seismic anisotropy of eclogite at various pressure‐temperature conditions. A 10 km thick subducted eclogitic crust can result in ∼0.2 s shear wave splitting in the Earth’s upper mantle.

Hao, Ming↗

DOE OSTI2026

Atomically Resolved Acoustic Dynamics Coupled with Magnetic Order in a Van der Waals Antiferromagnet

Magnetoelastic coupling in van der Waals (vdW) magnetic materials enables a unique interplay between the spin and lattice degrees of freedom. Characterizing the elastic responses with atomic and femtosecond resolution across the magnetic transition is essential for guiding the design of magnetically tunable actuators and strain-mediated spintronic devices. Here, ultrafast X-ray diffraction employed at a free-electron laser reveals that the atomic displacements, wave vectors, and dispersion relations of acoustic phonon modes in a vdW antiferromagnet FePS3 are coupled with the magnetic order, by tracking both in-plane and out-of-plane Bragg peaks upon optical excitation across the Néel temperature (TN). One transverse mode shows that a quasi-out-of-plane atomic displacement undergoes a significant directional change across TN. Its quasi-in-plane wave vector is derived by comparing the measured sound velocity and the first-principles calculations. The other transverse mode is an interlayer shear acoustic mode whose amplitude is strongly enhanced in the antiferromagnetic phase, exhibiting eight times stronger amplitude than the longitudinal acoustic mode below TN. The atomically resolved characterization of acoustic phonon dynamics that couple with magnetic ordering opens opportunities for harnessing unique magnetoelastic coupling in vdW magnets on ultrafast timescales.

coherent acoustic phonons↗

NASA NTRS2004

Lamb Wave Stiffness Characterization of Composites Undergoing Thermal-Mechanical Aging

The introduction of new, advanced composite materials into aviation systems requires a thorough understanding of the long term effects of combined thermal and mechanical loading upon those materials. Analytical methods investigating the effects of intense thermal heating combined with mechanical loading have been investigated. The damage mechanisms and fatigue lives were dependent on test parameters as well as stress levels. Castelli, et al. identified matrix dominated failure modes for out-of-phase cycling and fiber dominated damage modes for in-phase cycling. In recent years, ultrasonic methods have been developed that can measure the mechanical stiffness of composites. To help evaluate the effect of aging, a suitably designed Lamb wave measurement system is being used to obtain bending and out-of-plane stiffness coefficients of composite laminates undergoing thermal-mechanical loading. The system works by exciting an antisymmetric Lamb wave and calculating the velocity at each frequency from the known transducer separation and the measured time-of-flight. The same peak in the waveforms received at various distances is used to measure the time difference between the signals. The velocity measurements are accurate and repeatable to within 1% resulting in reconstructed stiffness values repeatable to within 4%. Given the material density and plate thickness, the bending and out-of-plane shear stiffnesses are calculated from a reconstruction of the dispersion curve. A mechanical scanner is used to move the sensors over the surface to map the time-of-flight, velocity, or stiffnesses of the entire specimen. Access to only one side of the material is required and no immersion or couplants are required because the sensors are dry coupled to the surface of the plate. In this study, the elastic stiffnesses D(sub 11), D(sub 22), A(sub 44), and A(sub 55) as well as time-of-flight measurements for composite samples that have undergone combined thermal and mechanical aging for a duration of 10,000 hours are reported.

Seale, Michael D.↗

NASA NTRS1997

Ultrasonic wave velocity measurement in small polymeric and cortical bone specimens

A system was refined for the determination of the bulk ultrasonic wave propagation velocity in small cortical bone specimens. Longitudinal and shear wave propagations were measured using ceramic, piezoelectric 20 and 5 MHz transducers, respectively. Results of the pulse transmission technique were refined via the measurement of the system delay time. The precision and accuracy of the system were quantified using small specimens of polyoxymethylene, polystyrene-butadiene, and high-density polyethylene. These polymeric materials had known acoustic properties, similarity of propagation velocities to cortical bone, and minimal sample inhomogeneity. Dependence of longitudinal and transverse specimen dimensions upon propagation times was quantified. To confirm the consistency of longitudinal wave propagation in small cortical bone specimens (< 1.0 mm), cut-down specimens were prepared from a normal rat femur. Finally, cortical samples were prepared from each of ten normal rat femora, and Young's moduli (Eii), shear moduli (Gij), and Poisson ratios (Vij) were measured. For all specimens (bone, polyoxymethylene, polystyrene-butadiene, and high-density polyethylene), strong linear correlations (R2 > 0.997) were maintained between propagation time and distance throughout the size ranges down to less than 0.4 mm. Results for polyoxymethylene, polystyrene-butadiene, and high-density polyethylene were accurate to within 5 percent of reported literature values. Measurement repeatability (precision) improved with an increase in the wave transmission distance (propagating dimension). No statistically significant effect due to the transverse dimension was detected.

Non-NASA Center↗

DOE OSTI2024

OpenFAST Modeling of the T-Omega Wind Floating Offshore Wind Turbine System

To reduce the cost of offshore wind energy through a more efficient design of the floating support structure, T-Omega Wind developed a novel lightweight and shallow-draft platform concept that aims to achieve a wave-following behavior without resonance amplification. The rotor and generator are carried by four tower legs with each leg supported by a shallow-draft float at the base. In preparation for a more detailed analysis, a coupled aero-hydro-elastic model of the proposed concept is developed in OpenFAST. The detailed modeling approach is presented, including a novel application of the SubDyn substructure dynamics module of OpenFAST to approximate the loads on the axle tube at the tower top from the rotor hub bearings. Preliminary results obtained with the OpenFAST model by rigidifying the substructure and blades indicate that the original sizing of the design can lead to large hub acceleration in the axial/surge and pitch directions due to platform-pitch motion. The high axial acceleration also potentially leads to large bending moments in the four tower legs. To address these issues, an updated design is developed with, among other changes, increased float spacing at the tower base and improved shear-transmitting geometry of the tower legs. The new design suggests significantly reduced extreme hub accelerations under the same conditions even with structural flexibility and will be further analyzed in the future through full aero-hydro-servo-elastic simulations.

FOWT↗

DOE OSTI2024

OpenFAST Modeling of the T-Omega Wind Floating Offshore Wind Turbine System

ENGINEERING,WIND ENERGY↗

DOE OSTI2022

Structure of Shock Waves and Inelasticity in Shock-Compressed Cemented Tungsten Carbides

In this study, shock wave experiments are conducted on General Carbide cemented tungsten carbide with 3.7wt.% cobalt binder to determine its shock-induced compression behavior up to 100 GPa. The measured wave profiles indicate the cemented tungsten carbide to undergo elastic-plastic deformation during shock compression. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot Elastic Limit (HEL), an elastic-plastic ramp indicating substantial post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, the shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relation, and the longitudinal stress (σ x ) vs. specific volume (V) curve for the samples. The HEL of the material was determined to lie between 4.41 GPa and 4.58 GPa. The U s - u p relation was determined to be U s = 4.97 + 1. 457u p for particle velocities greater than 0.75 km/s. The measured plastic shock velocities for particle velocities less than 0.7 km/s were found to be larger than those predicted using the linear U s - u p Hugoniot relationship, indicating the cemented WC samples to preserve substantial shear strength in the post-yield deformation region. No phase transformation was observed up to 100 GPa.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Shock compression of cemented tungsten carbides to 100 GPa: Structure of shock waves, Hugoniot relations, and strength

Plate impact experiments are conducted on cemented tungsten carbides (WC) with a 3.7 and 6.0 wt. % cobalt binder to better understand their dynamic, high-pressure response to 100 GPa. The measured wave profiles show propagation of steady structured waves. Standard impedance matching procedures are used to determine the Hugoniot relations in the shock velocity–particle velocity (U s –v p ) and Hugoniot stress–specific volume (P–V/V o ) planes. The Hugoniot elastic limit of the samples is controlled by ductility of the Co binder and is determined to be 4.45 ± 0.29 GPa for cemented WC with 3.7 wt. % cobalt and 3.72 ± 0.24 GPa for a 6.0 wt. % cobalt binder. Both grades show a non-linear U s –v p relationship depending on whether the particle velocity is in the strength dominated or hydrodynamic regime. In the strength dominated regime, a non-linear decrease in U s is observed as v p increases from ambient to the material’s hydrodynamic limit. In the hydrodynamic regime, the U s –v p Hugoniot is linear and is determined to be U s = 4.97(±0.006)+1.446(±0.018)v p km/s for WC with 3.7 wt. % Co and U s = 4.93(±0.006)+1.463(±0.017)v p km/s for 6 wt. % Co. Both WC grades indicate shear-stress hardening with mean stress immediately after yield, followed by pressure softening, and then a sharp fall in stress carrying capacity as the mean stress is increased to ≈70 GPa (hydrodynamic limit) and beyond. This behavior is in contrast to pure WC ceramics, which show continued shear-stress hardening with mean stress to ≈80 GPa.

36 MATERIALS SCIENCE↗

NASA NTRS1998

Development of In-Fiber Reflective Bragg Gratings as Shear Stress Monitors in Aerodynamic Facilities

Bragg gratings centered at nominal wavelengths of 1290 nm and 1300 run were inscribed in a 9/125 microns germano-silicate optical fiber, using continuous wave frequency doubled Ar+ laser radiation at 244 nm. Such gratings have been used extensively as temperature and strain monitors in smart structures. They have, however, never been used for measuring aerodynamic shear stresses. As a test of their sensitivity as shear stress monitors, a Bragg fiber attached to a metal plate was subjected to laminar flows in a glass pipe. An easily measurable large flow-induced wavelength shift (Delta Lambda(sub B)) was observed in the Bragg reflected wavelength. Thereafter, the grating was calibrated by making one time, simultaneous measurements of Delta Lambda(sub B) and the coefficient of skin friction (C(sub f)) with a skin friction balance, as a function of flow rates in a subsonic wind tunnel. Onset of fan-induced transition in the tunnel flow provided a unique flow rate for correlating Delta Lambda(sub B) and (C(sub f) values needed for computing effective modulus of rigidity (N(sub eff)) of the fiber attached to the metal plate. This value Of N(sub eff) is expected to remain constant throughout the elastic stress range expected during the Bragg grating aerodynamic tests. It has been used for calculating the value of Cf at various tunnel speeds, on the basis of measured values of Bragg wavelength shifts at those speeds.

Parmar, Devendra S.↗

NASA NTRS1993

Assessment of damage in ceramics and ceramic matrix composites using ultrasonic techniques

This paper addresses the application of ultrasonic sensing to damage assessment in ceramics and ceramic matrix composites. It focuses on damage caused by thermal shock or oxidation at elevated temperatures, which often results in elastic anisotropy. This damaged-induced anisotropy is determined by measuring the velocities of ultrasonic waves in different propagation directions. Thermal shock damage is assessed in ceramic samples of reaction bonded silicon nitride (RBSN). Thermal shock treatment from different temperatures up to 1000 C is applied to produce the microcracks. Results indicate that most microcracks produced by thermal shock are located near sample surfaces. Ultrasonic measurements using the surface wave method are found to correlate well with measurements of degradation of mechanical properties obtained independently by other authors using destructive methods. Oxidation damage is assessed in silicon carbide fiber/reaction bonded silicon nitride matrix (SCS-6/RBSN) composites. The oxidation is done by exposing the samples in a flowing oxygen environment at elevated temperatures, up to 1400 C, for 100 hr. The Youngs' modulus in the fiber direction as obtained from ultrasonic measurements decreases significantly at 600 C but retains its original value at temperatures above 1200 C. This agrees well with the results of destructive tests by other authors. On the other hand, the transverse moduli obtained from ultrasonic measurements decrease continually until 1200 C. Measurements on the shear stiffnesses show behavior similar to the transverse moduli. The results of this work show that the damage-induced anisotropy in both ceramics and ceramic matrix composites can be determined successfully by ultrasonic methods. This suggests the possibility of assessing damage severity using ultrasonic techniques.

Chu, Y. C.↗

DOE OSTI2022

Dynamic Strength of Iron at High Pressures and Strain Rates

Accurate modeling of meteorite impacts, and deformation of planetary cores require characterization of the flow strength and in-elasticity of iron in its different phases. In this letter, we investigate the flow strength of both the ambient α-phase and high-pressure ϵ-phase of iron at strain-rates of 1 × 10 5 s –1 and pressures up to 42 GPa using high pressure - pressure shear plate impact experiments. We report the strength of the ϵ-iron to be significantly higher than α-phase but consequently one order smaller than the previously reported dynamic strength at high pressures. As a result, the complete stress-strain response of the ϵ-phase is reported for the first time.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Linear Seismic Source Equivalents in 3D Nonlinear Models: Effects of Embedded Small-Scale, Near-Source Structures

Gaining a proper understanding of how Earth structure and other near-source properties affect estimates of explosion yield is important to the nonproliferation mission. The yields of explosion sources are often based on seismic moment or waveform amplitudes. Quantifying how the seismic waveforms or estimates of the source characteristics derived from those waveforms are influenced by natural or man-made structures within the near-source region, where the wavefield behaves nonlinearly, is required to understand the full range of uncertainty in those yield estimates. We simulate tamped chemical explosions using a nonlinear, shock physics code and couple the ground motions beyond the elastic radius to a linear elastic, full waveform seismic simulation algorithm through 3D media. In order to isolate the effects of simple small-scale 3D structures on the seismic wavefield and linear seismic source estimates, we embed spheres and cylinders close to the fully- tamped source location within an otherwise homogenous half-space. The 3 m diameters spheres, given their small size compared to the predominate wavelengths investigated, not surprisingly are virtually invisible with only negligible perturbations to the far-field waveforms and resultant seismic source time functions. Similarly, the 11 m diameter basalt sphere has a larger, but still relatively minor impact on the wavefield. However, the 11 m diameter air-filled sphere has the largest impact on both waveforms and the estimated seismic moment of any of the investigated cases with a reduction of ~25% compared to the tamped moment. This significant reduction is likely due in large part to the cavity collapsing from the shock instead of being solely due to diffraction effects . Although the cylinders have the same diameters as the 3 m spheres, their length of interaction with the wavefield produces noticeable changes to the seismic waveforms and estimated source terms with reductions in the peak seismic moment on the order of 10%. Both the cylinders and 11 m diameter spheres generate strong shear waves that appear to emanate from body force sources.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2021

Imaging the Shallow Structure of the Yucca Flat at the Source Physics Experiment Phase II Site with Horizontal-to-Vertical Spectral Ratio Inversion and a Large- N Seismic Array

The Source Physics Experiment (SPE) is a series of chemical explosions at the Nevada National Security Site (NNSS) with the goal of understanding seismic-wave generation and propagation of underground explosions. To understand explosion source physics, accurate geophysical models of the SPE site are needed. Here, we utilize a large-N seismic array deployed at the SPE phase II site to generate a shallow subsurface model of shear-wave velocity. The deployment consists of 500 geophones and covers an area of, approximately, 2.5 × 2 km. The array is located in the Yucca Flat in the northeast corner of the NNSS, Nye County, Nevada. Using ambient-noise recordings throughout the large-N seismic array, we calculate horizontal-to-vertical spectral ratios (HVSRs) across the array. We obtain 2D seismic images of shear-wave velocities across the SPE phase II site for the shallow structure of the basin. In this work the results clearly image two significant seismic impedance interfaces at ~150–500 and ~350–600 m depth. The shallower interface relates to the contrast between Quaternary alluvium and Tertiary volcanic rocks. The deeper interface relates to the contrast between Tertiary volcanic rocks and the Paleozoic bedrock. The 2D subsurface models support and extend previous understanding of the structure of the SPE phase II site. This study shows that the HVSR method in conjunction with a large-N seismic array is a quick and effective method for investigating shallow structures.

58 GEOSCIENCES↗

DOE OSTI2022

Deformations, relaxation, and broken symmetries in liquids, solids, and glasses: A unified topological field theory

Here, we combine hydrodynamic and field theoretic methods to develop a general theory of phonons as Goldstone bosons in crystals, glasses, and liquids based on nonaffine displacements and the consequent Goldstone phase relaxation. We relate the conservation, or lack thereof, of specific higher-form currents with properties of the underlying deformation field—nonaffinity—which dictates how molecules move under an applied stress or deformation. In particular, the single-valuedness of the deformation field is associated with conservation of higher-form charges that count the number of topological defects. Our formalism predicts, from first principles, the presence of propagating shear waves above a critical wave vector in liquids, thus giving a formal derivation of the phenomenon in terms of fundamental symmetries. The same picture provides also a theoretical explanation of the corresponding “positive sound dispersion” phenomenon for longitudinal sound. Importantly, accordingly to our theory, the main collective relaxation timescale of a liquid or a glass (known as the α relaxation for the latter) is given by the phase relaxation time, which is not necessarily related to the Maxwell time. Finally, we build a nonequilibrium effective action using the in-in formalism defined on the Schwinger-Keldysh contour, that further supports the emerging picture. In summary, our work suggests that the fundamental difference between solids, fluids, and glasses has to be identified with the associated generalized higher-form global symmetries and their topological structure, and that the Burgers vector for the displacement fields serves as a suitable topological order parameter distinguishing the solid (ordered) phase and the amorphous ones (fluids, glasses).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗