Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatial variations”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Facile alignment estimation in carbon nanotube films using image processing

Whether a macroscopic assembly of carbon nanotubes can exhibit the one-dimensional properties expected from individual nanotubes critically depends on how well the nanotubes are aligned inside the assembly. Therefore, a simple and accurate method for assessing the degree of alignment is desired for the rapid characterization of carbon nanotube films and fibers. Here, we present an end-to-end solution for determining the global and local spatial orientation of carbon nanotubes in films within a short amount of time using a fast, precise, and economical approach based on an image processing method applied to scanning electron microscopy images. Further, we first use Laplacian edge enhancement filtering for improving the appearance of edge regions, which is followed by image partitioning into multiple blocks to capture the nanoscale orientation characteristics and total variation-based image decomposition of these image blocks. We then perform a 2D-fast Fourier transform on the image decomposed textural components of these edge-enhanced image blocks to determine the orientation distribution, which is utilized to estimate the 2D nematic order parameter. To show the effectiveness of our method, we corroborated our results against results obtained with other state-of-the-art image processing and experimental techniques.

2D-FFT↗
DOE OSTI2025

Molecular Structure and Thermodynamics of CO 2 and Water Adsorption on Mica

The adsorption of CO 2 and water on clay surfaces plays a key role in applications, such as gas storage in saline aquifers and depleted hydrocarbon reservoirs, but is not yet fully understood. Here, we study the adsorption of CO 2 and water vapor using Grand Canonical Monte Carlo and molecular dynamics simulations. At a bulk pressure of 100 bar, pure CO 2 adsorbs strongly on mica and forms extensive layers next to it. CO 2 adsorption is lowered substantially if introducing water vapor above mica and is largely eliminated when the relative humidity (RH) approaches about 60%. When pure water vapor is introduced above a mica surface, a subnanometer thick liquid water film develops on it to form apparent liquid–solid and liquid–vapor interfaces simultaneously. Using the identification of truly interfacial molecules (ITIM) analysis, we delineate how individual water layers develop in this film as RH increases. We highlight that the water film is spatially heterogeneous and the true liquid–vapor interface emerges only at an RH of 60–80%. Introducing 100 bar of CO 2 into the water vapor above the mica surface modulates water adsorption nonlinearly: at RH = 0.01%, the water adsorption is reduced by ∼30%; as RH increases, the reduction is weakened, and eventually, enhancement of water adsorption by about 7% occurs at RH = 90%. These variations are attributed to the interplay of film thinning by high-pressure CO 2 , competition of mica surface sites by CO 2 molecules, and energetic and entropic stabilization of interfacial water by CO 2 molecules.

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

How Do Ionic Liquids “Fold” Ionenes? Computational and Experimental Analysis of Imidazolium Polymers Based on Ether and Alkyl Chain Variations Dissolved in an Ionic Liquid

Molecular dynamics (MD) simulations and complementary experiments are used here to understand the inter- and intramolecular structure and conformational properties of ionenes (i.e., cationic polymers formed from condensation reactions) dissolved in an ionic liquid (IL), 1-ethyl-3-methylimidazolium bistriflimide [C 2 mim + ][Tf 2 N - ]. The simulated structural properties are benchmarked against experimental analyses of these same polymers, mainly using dynamic light scattering (DLS) experiments. Four different imidazolium ionenes are considered, corresponding to variations in the chemistry and structure of the repeat units including: (1) poly (decylimidazolium) (PD 10 ); (2) poly (tetraethyleneglycolimidazolium) (PE 10 ); (3) alternating copolymer (P(ED) 5 ); and (4) a block copolymer (PE 5 D 5 ). Detailed computational analyses of the polymer structure and conformational properties were performed, including the radius of gyration, end-to-end distance, torsional distributions, and site-site and spatial distribution functions. Overall, there is a competition between intramolecular associations between the imidazolium groups and the ether sites in the polymers versus the intermolecular interactions of the polymer imidazolium groups with the surrounding anion molecules. The polymer with only ether linkages (PE 10 ) results in the most intramolecular interactions, leading to significant coiling behavior and chain contraction. These strong interactions reduce chain flexibility, but they also result in much more linear chain configurations and alignment of the imidazolium groups. The polymer structure is strongly affected as the concentration of alkyl groups is increased, and we find that the polymer architecture (alternating versus block copolymer) also has an important influence. Furthermore, these detailed observations are important for understanding the connection between the molecular design of imidazolium polymers and their emergent structural properties, which could lead to unique opportunities for creating polymer composites.

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

X-ray nano-imaging of defects in thin film catalysts via cluster analysis

Functional properties of transition-metal oxides strongly depend on crystallographic defects; crystallographic lattice deviations can affect ionic diffusion and adsorbate binding energies. Scanning x-ray nanodiffraction enables imaging of local structural distortions across an extended spatial region of thin samples. Yet, localized lattice distortions remain challenging to detect and localize using nanodiffraction, due to their weak diffuse scattering. Here, in this study, we apply an unsupervised machine learning clustering algorithm to isolate the low-intensity diffuse scattering in as-grown and alkaline-treated thin epitaxially strained SrIrO 3 films. We pinpoint the defect locations, find additional strain variation in the morphology of electrochemically cycled SrIrO 3 , and interpret the defect type by analyzing the diffraction profile through clustering. Our findings demonstrate the use of a machine learning clustering algorithm for identifying and characterizing hard-to-find crystallographic defects in thin films of electrocatalysts and highlight the potential to study electrochemical reactions at defect sites in operando experiments.

42 ENGINEERING↗
DOE OSTI2025

Arctic shrub expansion generates regional variation in litter decomposition by altering litter quality and the decomposition environment

Abstract The expansion of deciduous shrubs into the graminoid‐dominated arctic tundra is expected to alter litter decomposition by changing litter quality and local abiotic and biotic conditions. However, it remains unclear how shrub expansion will affect litter decomposition at regional scales, where macroclimate is expected to be the dominant regulator of decomposition. To determine the relative influence of macroclimate and local controls on regional patterns of litter mass loss and nitrogen release, we conducted two hierarchical litter decomposition experiments across spatial scales. We decomposed leaf and root litter from a prominent graminoid ( Eriophorum vaginatum ) and three genera of deciduous shrubs ( Betula , Alnus and Salix ) for 1 year within replicated plots at five sites spanning a 160 km latitudinal gradient in northern Alaska. Using Eriophorum litter as a substrate, we found that macroclimate was the primary regulator of mass loss but had opposing effects on leaf and root litter. As summer temperature increased along the latitudinal gradient (11.9 to 13.9°C), leaf litter mass loss increased by 20% whereas root litter mass loss decreased by 33%. Leaf nitrogen release also increased with summer temperature. Conversely, root nitrogen release was controlled by the vegetation type of the decomposition environment. Using different shrub litters as substrates, we found that litter quality and its interaction with soil microclimate and macroclimate controlled decomposition. Overall, shrub root litters decomposed faster than Eriophorum root litter, losing 53% more mass and 190% more nitrogen across all sites and decomposition environments. For leaf litter, however, patterns varied by litter genus, with Salix losing more mass and Betula and Alnus losing less mass than Eriophorum . Our findings demonstrate that shrub expansion in the Arctic can regulate leaf and root litter decomposition at the regional scale through its effects on local controls, primarily litter quality. Ongoing increases in shrub cover are likely to accelerate the turnover of root litter carbon and nitrogen pools in tundra ecosystems. Therefore, including shrub‐related processes in Earth system models will improve our ability to predict regional‐scale litter decomposition and its effects on carbon and nutrient cycling in a warming Arctic. Read the free Plain Language Summary for this article on the Journal blog.

Vozzo, Justin T. [Department of Natural Resources ↗
DOE OSTI2023

Watersheds and stream networks viewed longitudinally: Example insights from novel spatial portrayals of watershed characteristics

Abstract Longitudinal depictions of watershed structure and characteristics, including topography, stream networks, wetlands, ground water levels, and land use, can provide watershed knowledge and understanding unavailable from standard plan view maps. Three case studies provide examples of knowledge gained by applying longitudinal views of stream networks, watershed hydrologic behavior, and land use distributions. Longitudinal views of mountain stream networks show extreme variability in the slope‐area relationships of low Strahler order streams, large discontinuities in drainage area (large parts of drainage area space are absent in networks), and large variations in network curvature. Longitudinal views of a groundwater‐dominated headwater watershed increase the inference available from limited groundwater observations and clearly reveal how groundwater connections affect the permanence of surface water features and the distribution of vadose zone storage in the landscape. Plotting land uses longitudinally illuminates and allows a quantitative analysis of how land uses are distributed relative to topographic position. Viewing watersheds and stream networks longitudinally can provide new insights into watershed forms and processes and motivate new questions and research.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Peridynamic modeling of cementitious materials for nuclear waste management

Radioactive and hazardous waste generated from fuel processing plants, nuclear reactors, and hospitals, requires effective management strategies. Cementitious materials are widely applied for these needs, serving as structural materials, reactive barriers, or waste forms. In these applications cracking poses a significant risk to performance, driven by inconsistent shrinkage of the components and the varying strength and permeability of their interfaces. Traditional modeling approaches face challenges in representing the complex fracture behavior of cementitious materials due to the reliance on spatial derivatives and difficulties with mesh generation. Here, to overcome these limitations, we employ peridynamics, a novel continuum mechanics formulation that uses integrals to describe mechanical equilibrium, avoiding discontinuities associated with traditional methods. Through incorporation of a bi-linear softening model and quasistatics, an experimentally validated model for Portland cement concrete samples was created. Mechanical parameters, including compressive strength and elastic modulus were validated and variation due to aggregate packing was evaluated. Additionally, sensitivity analysis of the peridynamic parameters to the bulk material properties was established. The results lay the groundwork for evaluating the impact of unique conditions of cementitious waste forms that can be assessed to improve the reliability of waste management strategies.

Aggregates↗
DOE OSTI2022

Seeking Repeating Anthropogenic Seismic Sources: Implications for Seismic Velocity Monitoring at Fault Zones

Abstract Seismic velocities in rocks are highly sensitive to changes in permanent deformation and fluid content. The temporal variation of seismic velocity during the preparation phase of earthquakes has been well documented in laboratories but rarely observed in nature. It has been recently found that some anthropogenic, high‐frequency (>1 Hz) seismic sources are powerful enough to generate body waves that travel down to a few kilometers and can be used to monitor fault zones at seismogenic depth. Anthropogenic seismic sources typically have fixed spatial distribution and provide new perspectives for velocity monitoring. In this work, we propose a systematic workflow to seek such powerful seismic sources in a rapid and straightforward manner. We tackle the problem from a statistical point of view, considering that persistent, powerful seismic sources yield highly coherent correlation functions (CFs) between pairs of seismic sensors. The algorithm is tested in California and Japan. Multiple sites close to fault zones show high‐frequency CFs stable for an extended period of time. These findings have great potential for monitoring fault zones, including the San Jacinto Fault and the Ridgecrest area in Southern California, Napa in Northern California, and faults in central Japan. However, extra steps, such as beamforming or polarization analysis, are required to determine the dominant seismic sources and study the source characteristics, which are crucial to interpreting the velocity monitoring results. Train tremors identified by the present approach have been successfully used for seismic velocity monitoring of the San Jacinto Fault in previous studies.

58 GEOSCIENCES↗
DOE OSTI2022

Neural simulation-based inference approach for characterizing the Galactic Center γ-ray excess

The nature of the Fermi γ-ray Galactic Center Excess (GCE) has remained a persistent mystery for over a decade. Although the excess is broadly compatible with emission expected due to dark matter annihilation, an explanation in terms of a population of unresolved astrophysical point sources e.g., millisecond pulsars, remains viable. The effort to uncover the origin of the GCE is hampered in particular by an incomplete understanding of diffuse emission of Galactic origin. This can lead to spurious features that make it difficult to robustly differentiate smooth emission, as expected for a dark matter origin, from more “clumpy” emission expected from a population of relatively bright, unresolved point sources. Here, we use recent advancements in the field of simulation-based inference, in particular density estimation techniques using normalizing flows, in order to characterize the contribution of modeled components, including unresolved point source populations, to the GCE. Compared to traditional techniques based on the statistical distribution of photon counts, our machine learning-based method is able to utilize more of the information contained in a given model of the Galactic Center emission, and in particular can perform posterior parameter estimation while accounting for pixel-to-pixel spatial correlations in the γ-ray map. This makes the method demonstrably more resilient to certain forms of model misspecification. On application to Fermi data, the method generically attributes a smaller fraction of the GCE flux to unresolved point sources when compared to traditional approaches. We nevertheless infer such a contribution to make up a non-negligible fraction of the GCE across all analysis variations considered, with at least $38^{+9}_{-19}$% of the excess attributed to unresolved point sources in our baseline analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Removing the Effects of Tropical Dynamics from North Pacific Climate Variability

Teleconnections from the Tropics energize variations of the North Pacific climate, but detailed diagnosis of this relationship has proven difficult. Simple univariate methods, such as regression on El Niño-Southern Oscillation (ENSO) indices, may be inadequate since the key dynamical processes involved -- including ENSO diversity in the Tropics, re-emergence of mixed layer thermal anomalies, and oceanic Rossby wave propagation in the North Pacific -- have a variety of overlapping spatial and temporal scales. Here we use a multivariate Linear Inverse Model to quantify tropical and extra-tropical multi-scale dynamical contributions to North Pacific variability, in both observations and CMIP6 models. In observations, we find that the Tropics are responsible for almost half of the seasonal variance, and almost three quarters of the decadal variance, along the North American coast and within the subtropical front region northwest of Hawaii. SST anomalies that are generated by local dynamics within the Northeast Pacific have much shorter time scales, consistent with transient weather forcing by Aleutian low anomalies. Variability within the Kuroshio-Oyashio Extension (KOE) region is considerably less impacted by the Tropics, on all time scales. Consequently, without tropical forcing the dominant pattern of North Pacific variability would be a KOE pattern, rather than the Pacific Decadal Oscillation (PDO). In contrast to observations, most CMIP6 historical simulations produce North Pacific variability that maximizes in the KOE region, with amplitude significantly higher than observed. Correspondingly, the simulated North Pacific in all CMIP6 models is shown to be relatively insensitive to the Tropics, with a dominant spatial pattern generally resembling the KOE pattern, not the PDO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Ambient field seismology in critical zone hydrological sciences

Passive ambient noise monitoring is an emerging tool in environmental seismology, leveraging the ambient seismic field to assess temporal variations in shallow subsurface properties. This review focuses on the potential and challenges of using scattered coda waves from noise correlation functions to monitor critical zone dynamics. The sensitivity of seismic velocities to various environmental factors, including precipitation, snowmelt, atmospheric pressure, and groundwater fluctuations, underscores the method’s versatility. While coda waves excel in detecting subtle changes due to their scattered nature, ballistic waves provide higher spatial resolution, albeit with challenges in source stability. Advances in seismic sensing, including distributed acoustic sensing and low-cost geophone networks, have enabled high-resolution monitoring of hydrological processes, subsurface deformation, and seismic hazards. Integrating seismic data with hydrological models provides insights into water storage, pore pressure changes, and soil moisture dynamics. However, limitations in spatial resolution, calibration with ground truth data, and coupled effects between environmental factors remain key challenges. This review emphasizes the importance of interdisciplinary approaches in refining methodologies, enhancing sensor deployments, and addressing data gaps. Passive seismic monitoring offers opportunities to understand critical zone processes and their broader impacts on seismic hazards and environmental sustainability.

58 GEOSCIENCES↗
DOE OSTI2022

Understanding Synthesis and Structural Variation of Nanomaterials Through In Situ/Operando XAS and SAXS

Nanostructured materials with high surface area and low coordinated atoms present distinct intrinsic properties from their bulk counterparts. However, nanomaterials' nucleation/growth mechanism during the synthesis process and the changes of the nanomaterials in the working state are still not thoroughly studied. As two indispensable methods, X-ray absorption spectroscopy (XAS) provides nanomaterials' electronic structure and coordination environment, while small-angle X-ray scattering (SAXS) offers structural properties and morphology information. Here, a combination of in situ/operando XAS and SAXS provides high temporal and spatial resolution to monitor the evolution of nanomaterials. This review gives a brief introduction to in situ/operando SAXS/XAS cells. In addition, the application of in situ/operando XAS and SAXS in preparing nanomaterials and studying changes of working nanomaterials are summarized.

25 ENERGY STORAGE↗
DOE OSTI2021

Along the Channel Gradients Impact on the Spatioactivity of Gas Diffusion Electrodes at High Conversions during CO 2 Electroreduction

Results of a 2-D transport model for a gas diffusion electrode performing CO 2 reduction to CO with a flowing catholyte are presented, including the concentration gradients along the flow cell, spatial distribution of the current density and local pH in the catalyst layer. The model predicts that both the concentration of CO 2 and the buffer electrolyte gradually diminish along the channels for a parallel flow of gas and electrolyte as a result of electrochemical conversion and nonelectrochemical consumption. At high single-pass conversions, significant concentration gradients exist along the flow channels leading to large local variations in the current density (>150 mA/cm 2 ), which becomes prominent when compared to ohmic losses. In addition, concentration overpotentials change dramatically with CO 2 flow rate, which results in significant differences in outlet concentrations at high conversions. The outlet concentration of CO attains a maximum of 80% along with 5% CO 2 and 15% H 2 , although the maximum single-pass conversion is limited to below 60% due to homogeneous consumption by the electrolyte. Fundamental and practical implications of our findings on electrochemical CO 2 reduction are discussed with a focus on the trade-off between high current density operation and high single-pass conversion efficiency.

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

Dissecting Nearby Galaxies with piXedfit. I. Spatially Resolved Properties of Stars, Dust, and Gas as Revealed by Panchromatic SED Fitting

We study spatially resolved properties (on spatial scales of ~1–2 kpc out to at least 3 effective radii) of the stars, dust, and gas in 10 nearby spiral galaxies. The properties of the stellar population and dust are derived by fitting the spatially resolved spectral energy distribution (SED) with more than 20 photometric bands ranging from far-ultraviolet to far-infrared. Our newly developed software piXedfit performs point-spread function matching of images, pixel binning, and models the stellar light, dust attenuation, dust emission, and emission from a dusty torus heated by an active galactic nucleus simultaneously through the energy-balance approach. With this self-consistent analysis, we present the spatially resolved version of the IRX–βrelation, finding that it is consistent with the relationship from the integrated photometry. We show that the old stellar populations contribute to the dust heating, which causes an overestimation of the star formation rate (SFR) derived from the total ultraviolet and infrared luminosities on kiloparsec scales. With archival high-resolution maps of atomic and molecular gas, we study the radial variation of the properties of the stellar populations (including stellar mass, age, metallicity, and SFR), dust (including dust mass, dust temperature, and abundance of polycyclic aromatic hydrocarbon), and gas, as well as dust-to-stellar mass and dust-to-gas mass ratios. We observe a depletion of the molecular gas mass fraction in the central region of the majority of the galaxies, suggesting that the lack of available fuel is an important factor in suppressing the specific SFR at the center.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Spatial and temporal characterization of municipal solid waste based on resource recovery pathways

This study presents a two-year, quarterly assessment of MSW across four source sectors (residential, schools, restaurants, and grocery stores) from sixteen sites across five U.S. states. MSW was manually sorted into 27 categories and aggregated into pathway fractions: high-moisture (HM) organics, low-moisture (LM) organics, recyclable (RC) materials, and residuals for disposal. Organics represented 89 % of the MSW stream. The largest fraction was HM organics consisting of food waste (31 %) and yard waste (3 %), with large coefficient of variations (CV), 79 and 278 %, respectively, reflecting high seasonal and site variability that varied significantly (p < 0.01) across sampling periods. The HM fraction showed properties favorable for anaerobic digestion, with moisture content ranging from 56 to 95 % and volatile solids ranges of 86-95 %. In contrast, the LM and RC fractions remained more stable (plastics CV = 41 %; paper CV = 53 %) with heating values up to 26.9 MJ/kg across sources, reflecting suitability for gasification. Microstructural analysis revealed less porosity in residential waste sampled at the landfill, which can influence preprocessing efficiency and microbial accessibility. Pathway informed allocations showed that 35 % of MSW is suitable for anaerobic digestion, 36 % for gasification, and 18 % for recycling, leaving 11 % requiring landfill disposal. These results provide quantitative evidence to determine feedstock allocation, waste-to-energy system design, and the development of data-driven sustainability and resource recovery strategies within a circular bioeconomy.

09 BIOMASS FUELS↗
DOE OSTI2020

Microscale Observation via High-Speed X-ray Diffraction of Alloy 718 During In Situ Laser Melting

The laser melting process is accompanied by rapid evolution in temperature, phase, structure, and strain because of its high heating and cooling rates. In this study, the evolution of grains within a thin solid plate of Ni alloy 718 during laser processing was probed with in situ high-energy x-ray diffraction experiments. The high temporal and spatial resolution available in the measurement allowed us to study the rapid evolution of the melted region beneath the surface of the sample. The characterization of the evolution of secondary phases, i.e., Laves and carbide, was captured despite the weak diffracted peaks caused by small volume fractions. Thermal history was estimated based on changes in the lattice spacing from the thermal contraction upon cooling. In conclusion, the temporal variation in 2θ with azimuthal direction revealed the evolution in anisotropy of lattice spacing and thus of the mechanical state during laser processing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sulphur variations in annually layered stalagmites using benchtop micro-XRF

Variation of sulphur in annually laminated stalagmites can be used to infer the impact of past volcanic activities, anthropogenic pollution, and climate change due to increased bushfire activity. The synchrotron radiation micro-X-Ray fluorescence (SR-XRF) microprobe is a powerful tool to analyse and image sulphur recorded in stalagmites with micrometre resolution. However, access to SR-XRF beamlines can be limited, so researchers must select the most promising stalagmites for imaging. Benchtop micro-XRF is an effective tool for trace elemental analysis of speleothem samples and is a candidate for routine laboratory measurement of sulphur along stalagmite laminae and screening for SR-XRF. Here, this study describes a protocol using matrix-matched standards to measure annual variations of sulphur at trace to percent level along the laminae of two Western Australian stalagmites, one of which already having been analysed using SR-XRF. Parameters that affected quantitation include X-ray tube voltage and current, spot size of the X-ray beam and stalagmite surface roughness and porosity. The use of a 20 μm X-ray spot size provides sub-annual spatial resolution that can be completed in an overnight scan. The features in a 1000 point micro-XRF analysis of sulphur along a 20 mm transect show good consistency with SR-XRF microprobe data. Micro-XRF mapping was also performed to produce chemical images on the stalagmite and compared with Raman and X-ray diffraction (XRD) to confirm that the stalagmite is exclusively calcite, with no aragonite, and that the source of sulphur in the samples was gypsum and anhydrite. Regions of very high sulphur in the micro-XRF maps were found to be artefacts due to diffraction of the incident beam but these could be efficiently removed by using a multiple point statistics approach to produce a clean image suitable for analysis of the laminae. This work shows the potential of micro-XRF for routine analysis of sulphur in stalagmites, and to streamline sample characterisation before SR-XRF imaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗