Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion couples”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗
NASA NTRS2008

Dead Zone Accretion Flows in Protostellar Disks

Planets form inside protostellar disks in a dead zone where the electrical resistivity of the gas is too high for magnetic forces to drive turbulence. We show that much of the dead zone nevertheless is active and flows toward the star while smooth, large-scale magnetic fields transfer the orbital angular momentum radially outward. Stellar X-ray and radionuclide ionization sustain a weak coupling of the dead zone gas to the magnetic fields, despite the rapid recombination of free charges on dust grains. Net radial magnetic fields are generated in the magnetorotational turbulence in the electrically conducting top and bottom surface layers of the disk, and reach the midplane by ohmic diffusion. A toroidal component to the fields is produced near the midplane by the orbital shear. The process is similar to the magnetization of the solar tachocline. The result is a laminar, magnetically driven accretion flow in the region where the planets form.

solar system↗
DOE OSTI2021

Redox-stable symmetrical solid oxide fuel cells with exceptionally high performance enabled by electrode/electrolyte diffuse interface

Here, in this study, we report a high performance and redox-stable symmetrical solid oxide fuel cell (SOFC) based on (Ba 0.5 Sr 0.5 ) (Mo 0.1 Fe 0.9 )O 3-δ (BSMF) electrode and La 0.8 Sr 0.2 Ga 0.8 Mg 0.2 O 3-δ (LSGM) electrolyte. BSMF is able to operate both as anode and cathode. Excellent electrocatalytic activity has been achieved on BSMF towards hydrogen oxidation and oxygen reduction. Due to its closely matched lattice parameter to LSGM electrolyte, a unique diffuse interface is formed between BSMF and LSGM. Compared to a clean interface, e.g. BSMF/gadolinium doped ceria interface, this diffuse interface promotes the performance of BSMF electrode 1–1.8 times in 600–800 °C. Polarization resistance of the BSMF/LSGM specimen is as low as 0.047 and 0.007 Ωcm 2 in humidified H 2 and in air at 800 °C, respectively. On the BSMF/LSGM/BSMF symmetrical cell, a maximum power density of 2.28 W/cm 2 is achieved at 800 °C, the highest among with redox-stable ceramic electrodes to the best of our knowledge. Redox stability of this cell is confirmed. The role of anode and cathode is reversed back and forth in different operation modes. No apparent degradation is observed through 4 cycles within a 110 h operation period. These findings demonstrate that (Ba 0.5 Sr 0.5 ) (M o0.1 Fe 0.9 )O 3-δ coupled with LSGM electrolyte is an excellent choice to build a high performance, redox-stable SOFC.

25 ENERGY STORAGE↗
DOE OSTI2021

Impacts of Anisotropic Porosity on Heat Transfer and Off-Gassing during Biomass Pyrolysis

The pore structure of biogenic materials imbues the ability to deliver water and nutrients through a plant from root to leaf. This anisotropic pore granularity can also play a significant role in processes such as biomass pyrolysis that are used to convert these materials into useful products like heat, fuel, and chemicals. Evolutions in modeling of biomass pyrolysis as well as imaging of pore structures allow for further insights into the concerted physics of phase change-induced off-gassing, heat transfer, and chemical reactions. In this work, we report a biomass single particle model which incorporates these physics to explore the impact of implementing anisotropic permeability and diffusivity on the conversion time and yields predicted for pyrolysis of oak and pine particles. Simulation results showed that anisotropic permeability impacts predicted conversion time more than 2 times when the Biot number is above 0.1 and pyrolysis numbers (Py 1 , Py 2 ) are less than 20. Pore structure significantly impacts predicted pyrolytic conversion time (>8 times) when the Biot number is above 1 and the pyrolysis number is below 1, i.e., the “conduction controlled” regime. Therefore, these nondimensional numbers reflect that when internal heat conduction limits pyrolysis performance, internal pyrolysis off-gassing further retards effective heat transfer rates as a closely coupled phenomenon. Overall, this study highlights physically meaningful opportunities to improve particle-scale pyrolysis modeling and experimental validation relevant to a variety of feedstock identities and preparations, guiding the future design of pyrolyzers for efficient biomass conversion.

09 BIOMASS FUELS↗
DOE OSTI2026

Model‐Based Interpretation of Solute Exports and Carbon Partitioning During Shale Weathering in a Mountainous Hillslope

The weathering of sedimentary rocks in high-elevation catchments influences freshwater quality and the global carbon cycle. While individual biogeochemical mechanisms involved in this process are relatively well understood, quantifying their contributions to solute export and carbon fluxes under natural, transient conditions remains challenging. Here, we implement a numerical multidimensional and multiphase model to simulate coupled hydrological and biogeochemical processes in a shale-underlain, snow-dominated hillslope in the Rocky Mountains, Colorado. The model captures the dynamic interplay between soil respiration, mineral weathering, and climate-driven hydrological forcing, reproducing observed soil CO 2 dynamics, groundwater chemistry, and subsurface flow. Our results reveal that seasonal snowmelt enhances carbonate weathering by promoting the infiltration of CO 2 -rich water to depth, while pyrite oxidation is primarily sensitive to low water saturation that facilitates O 2 diffusion through the regolith. Topography modulates the spatial distribution of shale weathering, as steeper slopes enhance lateral drainage, favoring the delivery of reactants to greater depths. While shale weathering at our site acts as a transient carbon sink, with silicates and carbonates buffering acidity and promoting atmospheric CO 2 consumption (1% of soil-derived CO 2 ), the exported dissolved inorganic carbon is predominantly geogenic (∼73%). Consequently, when accounting for long-term marine carbonate precipitation. The current weathering regime represents a net source of carbon to the atmosphere. The oxidation of pyrite and petrogenic organic carbon together release approximately 0.9 mol·m −2 ·yr −1 of CO 2 . Our findings highlight the role of topography, hydroclimate, and the coupling between acid-base reactions in shaping the carbon balance and the solute exports in mountainous critical zones.

carbon cycling↗
NASA NTRS1983

HCN formation on Jupiter - The coupled photochemistry of ammonia and Acetylene

HCN formation in the upper troposphere and lower stratosphere of Jupiter is presently modeled in terms of UV pyrolysis of the C2H5N isomer aziridine, which is a product of the NH2 and C2H3 radicals that originate from ammonia photolysis and the addition of H atoms to acetylene, respectively. The sensitivity of the HCN column density to the individual rate constants and the eddy diffusion coefficient profile is considered, along with the possibility that additional HCN-yielding pathways may exist. Both ammonia and phosphine are strongly depleted by photolysis.

Kaye, J. A.↗
DOE OSTI2023

Surfactants Used in Colloidal Synthesis Modulate Ni Nanoparticle Surface Evolution for Selective CO 2 Hydrogenation

Colloidal chemistry holds promise to prepare uniform and size-controllable pre-catalysts; however, it remains a challenge to unveil the atomic-level transition from pre-catalysts to active catalytic surfaces under the reaction conditions to enable the mechanistic design of catalysts. In this work, we report an ambient-pressure X-ray photoelectron spectroscopy study, coupled with in situ environmental transmission electron microscopy, infrared spectroscopy, and theoretical calculations, to elucidate the surface catalytic sites of colloidal Ni nanoparticles for CO 2 hydrogenation. We show that Ni nanoparticles with phosphine ligands exhibit a distinct surface evolution compared with amine-capped ones, owing to the diffusion of P under oxidative (air) or reductive (CO 2 + H 2 ) gaseous environments at elevated temperatures. The resulting NiP x surface leads to a substantially improved selectivity for CO production, in contrast to the metallic Ni, which favors CH 4 . The further elimination of surface metallic Ni sites by designing multi-step P incorporation achieves unit selectivity of CO in high-rate CO 2 hydrogenation.

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

Spin-orbital correlations from complex orbital order in MgV 2⁢ O 4

MgV 2 ⁢O 4 is a spinel based on magnetic V 3+ ions, which host both spin (𝑆 = 1) and orbital (𝑙 eff = 1) moments. Owing to the underlying pyrochlore coordination of the magnetic sites, the spins in MgV 2 ⁢O 4 only antiferromagnetically order once the frustrating interactions imposed by the 𝐹⁡𝑑⁢$\overline{3}$⁢𝑚 lattice are broken through an orbitally-driven structural distortion at T 𝑆 ≃ ⁢60 K. Consequently, a Néel transition occurs at T 𝑁 ≃ ⁢40 K. Low-temperature spatial ordering of the electronic orbitals is fundamental to both the structural and magnetic properties; however, considerable discussion on whether it can be described by complex or real orbital ordering is ambiguous. We apply neutron spectroscopy to resolve the nature of the orbital ground state and characterize hysteretic spin-orbital correlations using x-ray and neutron diffraction. Neutron spectroscopy finds multiple excitation bands and we parametrize these in terms of a multilevel (or excitonic) theory based on the orbitally degenerate ground state. Meaningful for the orbital ground state, we report an “optical-like” mode at high energies that we attribute to a crystal-field-like excitation from the spin-orbital 𝑗 eff = 2 ground-state manifold to an excited 𝑗 eff =1 energy level. We parametrize the magnetic excitations in terms of a Hamiltonian with spin-orbit coupling and local crystalline electric field distortions resulting from deviations from perfect octahedra surrounding the V 3+ ions. We suggest that this provides compelling evidence for complex orbital order in MgV 2⁢ O 4 . We then apply the consequences of this model to understand hysteretic effects in the magnetic diffuse scattering where we propose that MgV 2⁢ O 4 displays a high-temperature orbital memory of the low-temperature spin order.

Frustrated magnetism↗
DOE OSTI2021

Standoff Detection of Oil and Powder Mixtures at 12 Meters Using a Tunable Quantum Cascade Laser-Based System with a Close Focus Telescope and Uncooled Infrared Detector

We have designed and demonstrated a quantum cascade laser (QCL) based standoff system that utilizes an uncooled mercury cadmium telluride (MCT) detector with lock-in signal processing for chemical identification at a distance of 12.5 meters in indoor ambient light conditions. In the system, a tunable quad-QCL operating (1 MHz) in quasi-continuous wave mode between 8.45 and 10.03 μm (~1182 to 1000 cm –1 ) serves as the active mid-infrared source for remotely interrogating mineral, powder, and thin film oil samples including powder mixtures (6, 12.5, 25, and 50%) of crystalline quartz (SiO 2 ) in KBr. Light as reflected from a given sample is collected using a 10-inch (25.4 cm) Dall Kirkham telescope and coupled with ZnSe optics to an uncooled MCT detector. Furthermore, the mixture dependence of the highly transparent KBr and strongly absorbing quartz was found to fit a modified version of the Schatz reflectance model for compacted powder mixtures. All reflectance spectra reported are relative to an Au-coated diffuse reflector. A NIST traceable polystyrene standard reflector was also used to determine the QCL wavelength tuning range and calibration.

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

Comparison of Zirconium Redistribution in BISON EBR-II Models Using FIPD and IMIS Databases with Experimental Post Irradiation Examination

Metallic fuels have seen increased interest for future sodium fast reactors due to their material properties: high thermal conductivities and advantageous neutronic properties allow for greater fission densities. One drawback to typical metallic fuels is zirconium redistribution, which impacts this advantageous material and its neutronic properties. Unfortunately, the processes behind zirconium migration behavior are understood using first principles, so before these fuels are implemented in future fast reactors, characterization and fuel qualification regimes must be completed. These activities can be supported through the use of robust modeling using the most accurate empirical models currently available to fuel researchers around the world. The tool that allows researchers to model this complex coupled thermo-mechanical behavior and nuclear properties is BISON. Additionally, BISON model parameters need to be compared against PIE measurements. The current work utilizes two fuel pins from EBR-II experiment X441 to optimize various model parameters, including porosity correction factor, thermal conductivity, phase transition temperature, and diffusion coefficient multipliers, before implementing the final model for seven fuel pins with differing characteristics. To properly evaluate the BISON simulations, the results are compared to PIE metallography data for each fuel pin, to ensure the zirconium redistribution is properly reflected in the simulation results. Six out of seven analyzed fuel pins demonstrate good agreement between the metallography images and BISON results, showing alignment of the Zr-rich, Zr-depleted, and moderately Zr-enriched zones at various axial heights along the fuel pins. Further work is needed to refine the model parameters for general pin use.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1972

Silicon-germanium technology program of the Jet Propulsion Laboratory.

The outer planetary exploration missions studied by the Jet Propulsion Laboratory require silicon-germanium radioisotope thermoelectric generators (RTGs) in which the factors of safety are as low as is compatible with the reliable satisfaction of the power needs. Work on silicon germanium sublimation performed at the Jet Propulsion Laboratory is presented. Analytical modeling work on the solid-diffusion process involved in the steady-state (free) sublimation of silicon germanium is described. Analytical predictions are made of the sublimation suppression which can be achieved by using a cover gas. A series of accelerated (high-temperature) tests which were performed on simulated hardware (using four SiGe couples) to study long-term sublimation and reaction mechanisms are also discussed.

De Winter, F.↗
NASA NTRS1991

Modeling of InP metalorganic chemical vapor deposition

The growth of InP by metalorganic chemical vapor deposition (MOCVD) in a horizontal reactor is being modeled with a commercially available computational fluid dynamics modeling code. The mathematical treatment of the MOCVD process has four primary areas of concern: 1) transport phenomena, 2) chemistry, 3) boundary conditions, and 4) numerical solution methods. The transport processes involved in CVD are described by conservation of total mass, momentum, energy, and atomic species. Momentum conservation is described by a generalized form of the Navier-Stokes equation for a Newtonian fluid and laminar flow. The effect of Soret diffusion on the transport of particular chemical species and on the predicted deposition rate is examined. Both gas-phase and surface chemical reactions are employed in the model. Boundary conditions are specified at the inlet and walls of the reactor for temperature, fluid flow and chemical species. The coupled set of equations described above is solved by a finite difference method over a nonuniform rectilinear grid in both two and three dimensions. The results of the 2-D computational model is presented for gravity levels of zero- and one-g. The predicted growth rates at one-g are compared to measured growth rates on fused silica substrates.

Black, Linda R.↗
NASA NTRS2003

Modeling Charge Collection in Detector Arrays

A detector array charge collection model has been developed for use as an engineering tool to aid in the design of optical sensor missions for operation in the space radiation environment. This model is an enhancement of the prototype array charge collection model that was developed for the Next Generation Space Telescope (NGST) program. The primary enhancements were accounting for drift-assisted diffusion by Monte Carlo modeling techniques and implementing the modeling approaches in a windows-based code. The modeling is concerned with integrated charge collection within discrete pixels in the focal plane array (FPA), with high fidelity spatial resolution. It is applicable to all detector geometries including monolithc charge coupled devices (CCDs), Active Pixel Sensors (APS) and hybrid FPA geometries based on a detector array bump-bonded to a readout integrated circuit (ROIC).

Hardage, Donna↗
NASA NTRS2009

Solar Wind Ablation of Terrestrial Planet Atmospheres

Internal plasma sources usually arise in planetary magnetospheres as a product of stellar ablation processes. With the ignition of a new star and the onset of its ultraviolet and stellar wind emissions, much of the volatiles in the stellar system undergo a phase transition from gas to plasma. Condensation and accretion into a disk is replaced by radiation and stellar wind ablation of volatile materials from the system- Planets or smaller bodies that harbor intrinsic magnetic fields develop an apparent shield against direct stellar wind impact, but UV radiation still ionizes their gas phases, and the resulting internal plasmas serve to conduct currents to and from the central body along reconnected magnetic field linkages. Photoionization and thermalization of electrons warms the ionospheric topside, enhancing Jeans' escape of super-thermal particles, with ambipolar diffusion and acceleration. Moreover, observations and simulations of auroral processes at Earth indicate that solar wind energy dissipation is concentrated by the geomagnetic field by a factor of 10-100, enhancing heavy species plasma and gas escape from gravity, and providing more current carrying capacity. Thus internal plasmas enable coupling with the plasma, neutral gas and by extension, the entire body. The stellar wind is locally loaded and slowed to develop the required power. The internal source plasma is accelerated and heated, inflating the magnetosphere as it seeks escape, and is ultimately blown away in the stellar wind. Bodies with little sensible atmosphere may still produce an exosphere of sputtered matter when exposed to direct solar wind impact. Bodies with a magnetosphere and internal sources of plasma interact more strongly with the stellar wind owing to the magnetic linkage between the two created by reconnection.

Moore, Thomas Earle↗
DOE OSTI2025

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗
DOE OSTI2022

Vapor–Fed Electrolyzers for Carbon Dioxide Reduction Using Tandem Electrocatalysts: Cuprous Oxide Coupled with Nickel–Coordinated Nitrogen–Doped Carbon

We report ethylene is particularly attractive due to its major importance as a feedstock for various applications including the polymer industry. As such, catalyst and electrolyzer developments are crucial to achieve industrially relevant ethylene production and efficiency levels. In this work, a tandem electrocatalyst composed of copper nanocubes and nickel-coordinated nitrogen-doped carbon (NiNC) is presented, which is integrated into gas diffusion electrodes (GDEs) for direct conversion of vapor-fed CO 2 into ethylene. Evaluation of tandem GDEs in the vapor-fed flow electrolyzer shows significantly increased ethylene selectivity in terms of faradaic efficiency and C 2 H 4 /CO ratio compared to a non-tandem copper GDE. The enhancements are attributed to the increased local CO availability near the copper surface via effective CO 2 to CO conversion on neighboring NiNC. The experimental results are validated by 3D resolved continuum simulations, which show increased flux of higher-order prod-ucts with the added CO flux from NiNC. The practical viability of Cu/NiNC catalyst is further evaluated in a membrane electrode assembly electrolyzer, achieving 40% FE toward ethylene at 150 mA cm -2 and 3.2 V. These find-ings highlight the high selectivity and formation rate of ethylene achieved by successful device integration of the Cu/NiNC catalyst, demonstrating the potential for implementation in large-scale sustainable CO 2 electrolyzers.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Review of Recent Results in Global MHD Modeling: ISTP Project Scientist for Theory and Ground-Based Observations

Global MHD (magnetohydrodynamic) simulations have shown a remarkable ability to describe the global dynamics of geospace. The limitations of the physical approximations underlying MHD would seem to limit the effectiveness of these codes, since kinetic and hybrid effects should manifest themselves by cross-scale coupling from microscales to mesoscales to global scales. However three effects appear to allow the codes to operate much more successfully than one would at first believe. They are:(l) the globally self-consistent nature of the codes with very well defined exterior boundary conditions (the solar wind) which allows the proper intercommunication between magnetospheric regions on MHD scales, (2) the control by global dynamics of the boundary layer locations where micro and meso scale processes operate, and (3) the critical role of numerical diffusion and with a sufficiently high resolution grid, the use of an empirical resistivity term, which if set at a level where the major magnetosphere boundaries properly calibrate against their observed locations, appear to well represent the effects of kinetic and hybrid processes on the global dynamics. The effectiveness of the global MHD codes, which have been developed under the ISTP mission, in describing Wind, Polar and Geotail observations, as well as ground-based observations are described. Particular emphasis is placed upon the Polar imaging data which when combined with ground-based data and global MHD-based synthetic aurora and convection patterns provide a powerful tool in understanding the final link in the solar-terrestrial chain: coupling into the atmosphere and ionosphere.

Curtis, Steven↗
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↗