Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical Dynamics Initiative”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Comparison of two pressure–temperature equilibration methods

We compare and contrast the traditionally used method of solving the pressure–temperature equilibration problem in hydrodynamics, where specific internal energy and density are considered independent variables, with a different method where pressure and temperature are independent variables. With the goal of examining the robustness of the two methods as the number of components increases, we examine 2-, 4-, 6-, and 8-component systems. After equilibrating more than 10 4 initial conditions for each system using both methods, we demonstrate that the latter method constrains the search space by lowering its dimensionality and forces a better initial guess, resulting in a higher probability of convergence to solution with fewer, cheaper iterations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Equilibrium-gated pattern formation: How molecular dissociation thermodynamics drive emergent behavior in dissipative polymeric systems

Emergent patterns in biological systems arise through dissipative processes that balance reaction and transport phenomena, producing highly functional properties from self-regulating mechanisms. Synthetic fabrication, by contrast, often relies on user-controlled, multistep methods that lack the self-organizing capabilities of natural systems. Inspired by nature, we sought chemical systems that integrate strongly coupled reaction and transport phenomena, identifying frontal ring-opening metathesis polymerization (FROMP) as a method capable of creating diverse forms and functions through reactive processing. By employing discrete molecular initiators, FROMP allows precise control of key reaction steps—inhibition, initiation, and propagation. Using an integrated computational and experimental framework, we uncover how near-equilibrium inhibition dynamics, coupled with far-from-equilibrium reaction kinetics, drive pattern formation in frontally polymerized synthetic materials. We propose the concept of equilibrium-gated pattern formation, demonstrating how initiator chemistry can be tuned to achieve programmable macroscale properties. Our study reveals a surprising insight: Emergent behavior in FROMP systems arises from the inhibition-dominated regime of resin composition, expanding prior observations that such behavior is confined to a narrow compositional space near the boundary between front quenching and uniform front propagation. We identify a broader compositional window, far from the quenching regime, where emergent behavior reliably manifests. This expanded design space significantly enhances the operational flexibility of reactive systems and their capacity for self-organization. Furthermore, these insights provide a roadmap for designing bioinspired materials with self-organizing capabilities, unlocking possibilities in synthetic manufacturing.

Architected Polymers↗
DOE OSTI2023

Atomistic simulations for investigation of substrate effects on lipid in-source fragmentation in secondary ion mass spectrometry

In beam-based ionization methods, the substrate plays an important role on the desorption mechanism of molecules from surfaces. Both the specific orientation that a molecule adopts at a surface and the strength of the molecule-surface interaction can greatly influence desorption processes, which in turn will affect the ion yield and the degree of in-source fragmentation of a molecule. In the beam-based method of secondary ion mass spectrometry (SIMS), in-source fragmentation can be significant and molecule specific due to the hard ionization method of using a primary ion beam for molecule desorption. To investigate the role of the substrate on orientation and in-source fragmentation, we have used atomistic simulations—molecular dynamics in combination with density functional theory calculations—to explore the desorption of a sphingolipid (palmitoylsphingomyelin) from a model surface (gold). We then compare SIMS data from this model system to our modeling findings. Using this approach, we found that the combined adsorption and binding energy of certain bonds associated with the headgroup fragments (C3H8N+, C5H12N+, C5H14NO+, and C5H15PNO4+) was a good predictor for fragment intensities (as indicated by relative ion yields). This is the first example where atomistic simulations have been applied in beam-based ionization of lipids, and it presents a new approach to study biointerfacial lipid ordering effects on SIMS imaging.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Dynamic System Scaling Application to Accelerated Nuclear Fuel Testing

The development of nuclear fuel and materials requires a continuous effort to investigate the acute and prolonged effects of irradiation, thermal-material stress, chemical change, or other conceivable damage mechanics acquired during normal operation or accident scenarios throughout its lifetime. As in-core fuel property measurement techniques advance to support in the real-time, non-invasive, and enhanced accuracy realm, it is the future of fuel development to pursue to higher degrees of control, predictability of integral test behavior via separate effects test (SET), and shorter test time intervals. The fuel development life cycle from initial concept to commercial licensing is approximated to be 20 years and current literature suggests by optimizing fuel performance codes with SETs, the process could possibly be compressed to 5 to 10 years. Recently, a research group in the Idaho National Laboratory (INL) is testing reduced scale fuel rods and increased power density to accelerate evolution of fuel phenomena in metallic fuels. In support of nuclear fuel rod development, compressing fuel test process, and accelerating fuel phenomena, it was the purpose of this study to investigate nuclear fuel performance phenomena via literature review and effectively scale the initial conditions, boundary conditions, and geometric properties to describe the time-dependent response including to fuel burnup, thermal and mechanical stress, transmutations and inter-diffusion, and other relevant observed phenomena. The study was based on BISON simulations of historic EBR-II metallic fuel experiments and Dynamical System Scaling (DSS) method are utilized to assess effects of scaling activity including fuel phenomena acceleration and calculations of time-dependent distortions. The research successfully scaled metallic fuel phenomena, accelerated fuel testing, and assessed the distortions for each scaled case derived.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

A Scaling Analysis of Post-Detonation Mixing With Detailed Chemical Kinetics

This report documents a scaling study of hemispherical pentaerythritol tetranitrate explosives detonated in ambient air. The focus is on the mixing layer and evolution of the chemical species. The HyBurn code developed at the University of Florida is used to simulate center initiated 12, 38, 120, and 1200 mm diameter charges. The results show blast structures commonly found in the literature. Azimuthal averaging is used to radially average the flow variables to get mean values. The charges exhibit self-similar behavior at early times in a primarily gas dynamic regime, but discrepancies between the scales form at later times due to the finite rate chemistry. The turbulent mixing layer is calculated and used to define a region to conditionally average by mixture fraction. The results of the conditional averages imply a transition from a kinetic-limited regime to a mixing-limited regime with increasing scale. The scaling of the total heat release rate confirms an early time gas dynamic regime, and a later time regime where mixing and chemistry becomes impactful.

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

MULTI-DIMENSIONAL MODELING OF THE CFR ENGINE FOR THE INVESTIGATION OF SI NATURAL GAS COMBUSTION AND CONTROLLED END-GAS AUTOIGNITION

Engine knock and misfire are barriers to pathways leading to high-efficiency Spark-Ignited (SI) Natural Gas engines. The general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by low emission limits and the wide range of chemical reactivity in pipeline quality natural gas. Depending on the region and the source of the natural gas, its reactivity, described by its methane number (analogous to the octane number for liquid SI fuels) can span from 65 - 95. In order to realize diesel-like efficiencies, SI natural gas engines must be designed to operate at high BMEP near knock limits over a wide range of fuel reactivity. This requires a deep understanding regarding the combustion-engine interactions pertaining to flame propagation and end-gas autoignition (EGAI). However, EGAI, if controlled, provides an opportunity to increase SI natural gas engine efficiency by increasing combustion rate and the total burned fuel, mitigating the effects of the slow flame speeds of natural gas fuels which generally reduce BMEP and increase unburned hydrocarbon emissions. For this reason, in order to study EGAI phenomenon, the present work highlights multi-dimensional computational fluid dynamics (CFD) models of the Cooperative Fuel Research (CFR) engine. The CFR engine models are used to investigate fuel-engine interactions that lead to EGAI with natural gas, including effects of fuel reactivity, engine operating parameters, and exhaust gas recirculation (EGR). A Three-Pressure Analysis, performed with GT-Power, was used to estimate initial and boundary conditions for the three-dimensional CFD model. CONVERGE CFD v2.4 was used for the three-dimensional CFD modeling where the level set G-Equation model and SAGE detailed chemical kinetics solver were used. An assessment of the different modeling approaches is also provided to evaluate their limitations, advantages and disadvantages, and for which situations they are most applicable. Model validation was performed with experimental data taken with a CFR engine over varying compression ratio, CA50, EGR fraction, and IMEP and shows good agreement in Peak Cylinder Pressure (PCP), PCP crank angle, and the location of the 10%, 50%, and 90% mass fraction burned (CA10, CA50, and CA90, respectively). The models can predict the onset crank angle and pressure rise rate for light, medium, and heavy EGAI under a variety of fuel reactivities and engine operating conditions.

Bestel, Diego↗
DOE OSTI2025

Why are Lead Iodide‐Based Perovskite Precursor Inks Yellow?

A challenge faced by metal halide perovskite (MHP) photovoltaics is scaling up solution deposition processes to realize rapid and inexpensive manufacturing. The challenge lies in completely understanding and controlling solution speciation, nucleation, and self-assembly of iodoplumbate complexes during solvent evaporation as the liquid transforms into gels and solids. An accurate description of solution species, at all points in the transformation, is a prerequisite to design robust and reliable processes. Here, in this study, the common assumption that initial monoplumbate solution species typically invoked (e.g., [PbI 6 ] 4− ) are certainly not the origin of optical absorbance at >400 nm wavelengths is disproved, as are many large particles of common “intermediate” iodoplumbate phases with face- or edge-sharing connectivity. Instead, a new perspective is offered, involving (partially) corner-sharing iodo(poly)plumbates (>1 Pb 2+ per complex) that experience highly dynamic chemical environments. It is outlined how the MHP field would benefit by elucidating these phenomena. Future work is required to determine the size and kinetic behavior of polyplumbate species, and contextualize these findings in relation to broader trends in materials chemistry beyond MHPs. Ultimately, a complete explanation for the solution speciation, optical absorbance signatures, and the color of MHP precursor inks remains an open challenge to the community.

14 SOLAR ENERGY↗
DOE OSTI2021

Modeling of thermal and kinetic processes in non-equilibrium plasma ignition applied to a lean combustion engine

In recent years novel ignition systems have been developed to enable stable and efficient engine operations with lean mixtures. Among them, radio-frequency corona ignition systems create discharges that involve a much wider region compared to traditional spark, and produce non equilibrium plasma with high levels of active radicals and excited species. These devices considerably increase the early flame growth speed and extend stable operating limits. With the aim of expanding the knowledge on high efficiency lean-burn SI engines, this paper investigates and compares the combustion development generated by spark and corona ignitions through computational fluid dynamics, within the Reynolds-Averaged Navier-Stokes framework for turbulence modeling. In order to simultaneously take thermal and chemical effects into account, the Perfectly Stirred Reactor combustion model is used. Experimental data are also collected for validation in an optical access engine, for different mixture levels, from stoichiometric to very lean. Furthermore, the faster burn rate generated by the corona system in the initial stage of the combustion is well predicted by the simulations, in all the relative air-fuel ratio conditions. Remarkably, as the mixture becomes lean, simulations are able to capture the non-linear transition from fast to slow kernel growth, before a self-sustainable flame propagation is established. This correlates very well with the measured engine cyclic variability and the corresponding steep change in the duration of the flame kernel formation. Ultimately, this study highlights the important role of the atomic oxygen, as active radical, in promoting and enhancing the combustion initiated by a corona discharge, in addition to the volumetric ignition effect. By contrast, the validated simulations allow to explain that the high-temperature thermal plasma generated in a traditional spark discharge is insensitive to kinetic aspects.

42 ENGINEERING↗
DOE OSTI2021

Interfacial Response and Structural Adaptation of Structured Polyelectrolyte Thin Films

Ionizable block copolymers with distinctive block characteristics display the diversity crucial for the design of macromolecules for targeted applications. In contrast to van der Waals copolymers, their interfaces, which are critical to their function, consist of nanodomains, each of a different nature and thus unique interfacial behavior. Here, the interfacial response of a symmetric block copolymer with a sulfonated polystyrene polyelectrolyte center, tethered to polyethylene-r-propylene and terminated by poly(t-butyl styrene) is probed as polymer films are exposed to three polar solvents, water, propanol, and tetrahydrofuran (THF), using molecular dynamics simulations. Each of the solvents captures a distinctive interaction with the individual blocks. We find that at the film boundary, the interfacial response is initially dominated by that of the hydrophobic blocks to all solvents. At later times, the solvent distribution among the blocks, where water molecules associate predominantly with the sulfonated groups and propanol and THF reside at multiple different sites, determines the chemical composition and the polymer conformation at the interface. Altogether, these simulations provide the first direct molecular insight into the interfacial response of ionizable copolymers.

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

Modeling Approach for the Aluminum-clad Dry Storage Pilot using HFIR Fuel

To confirm that the dry storage of aluminum-clad research reactor spent nuclear fuel (ASNF) will remain within the safety envelope after applied drying schemes and that the resulting evolution of the gas space composition, temperature, and pressure conditions are understood, a dry storage pilot project is being established. The pilot will incorporate an instrumented lid for discrete interval or for on-demand gas composition and temperature monitoring of two DOE Standard Canisters (DSCs) loaded with three High Flux Isotope Reactor (HFIR) inner cores per DSC. Each DSC would be subjected to a separate alternative candidate drying scheme. Canisters will undergo 1 to 5 years of monitoring, including internal temperature and gas sampling to track pressure and composition changes. This report outlines the approach for modeling the ASNF-in-canister behavior in terms of evolving gas space conditions for the ASNF dry storage pilot using HFIR fuel. The ASNF has an adherent surface oxyhydroxide layer comprised of boehmite/bayerite that generates hydrogen when subjected to irradiation. Three-dimensional multi-physics computational fluid dynamics simulations will be executed to compute the thermal field within the DSC and provide inputs to a chemical model employed to compute pressure buildup as hydrogen is generated in the system. Implemented in Cantera, the chemical model solves gas phase and aluminum oxyhydroxide surface-mediated radiolysis reactions. Gas phase reactions are sourced from Wittman and Hanson (2015), whereas surface-mediated reactions are incorporated by fitting experimental data using an optimization algorithm (Abboud, 2023). Water radiolysis reactions from Wren and Ball (2001) are adopted with modifications as described in Abboud (2023c). Understanding the effect of the hydrogen buildup over time is important for long-term storage safety considerations. Modeling results will include the canister pressure, temperature, and composition evolution from the initial helium backfill with the addition of radiolytically-evolved chemical species (e.g., hydrogen and oxygen). The specific HFIR cores for the pilot program have not yet been selected, and the overall design is still in development. The CFD-chemical model used for this work will be based on prior models with necessary updates to allow for improved accuracy and efficiency. The experimental data obtained from the HFIR demonstration will be used to improve and validate the computational models to predict the ASNF-in-canister behavior.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Liquid-induced topological transformations of cellular microstructures

The fundamental topology of cellular structures—the location, number and connectivity of nodes and compartments—can profoundly affect their acoustic, electrical5, chemical, mechanical and optical properties, as well as heat, fluid and particle transport. Approaches that harness swelling, electromagnetic actuation and mechanical instabilities in cellular materials have enabled a variety of interesting wall deformations and compartment shape alterations, but the resulting structures generally preserve the defining connectivity features of the initial topology. Achieving topological transformation presents a distinct challenge for existing strategies: it requires complex reorganization, repacking, and coordinated bending, stretching and folding, particularly around each node, where elastic resistance is highest owing to connectivity. Here we introduce a two-tiered dynamic strategy that achieves systematic reversible transformations of the fundamental topology of cellular microstructures, which can be applied to a wide range of materials and geometries. Our approach requires only exposing the structure to a selected liquid that is able to first infiltrate and plasticize the material at the molecular scale, and then, upon evaporation, form a network of localized capillary forces at the architectural scale that ‘zip’ the edges of the softened lattice into a new topological structure, which subsequently restiffens and remains kinetically trapped. Reversibility is induced by applying a mixture of liquids that act separately at the molecular and architectural scales (thus offering modular temporal control over the softening–evaporation–stiffening sequence) to restore the original topology or provide access to intermediate modes. Guided by a generalized theoretical model that connects cellular geometries, material stiffness and capillary forces, we demonstrate programmed reversible topological transformations of various lattice geometries and responsive materials that undergo fast global or localized deformations. In conclusion, we then harness dynamic topologies to develop active surfaces with information encryption, selective particle trapping and bubble release, as well as tunable mechanical, chemical and acoustic properties.

42 ENGINEERING↗
DOE OSTI2023

Structure and thermodynamics of water adsorption in NU-1500-Cr

Abstract Metal-organic frameworks (MOFs) are a class of materials with diverse chemical and structural properties, and have been shown to effectively adsorb various types of guest molecules. The mechanism of water adsorption in NU-1500-Cr, a high-performance atmospheric water harvesting MOF, is investigated using a combination of molecular dynamics simulations and infrared spectroscopy. Calculations of thermodynamic and dynamical properties of water as a function of relative humidity allow for following the adsorption process from the initial hydration stage to complete filling of the MOF pores. Initial hydration begins at the water molecules that saturate the open Cr 3+ sites of the framework, which is then followed by the formation of water chains that extend along the channels connecting the hexagonal pores of the framework. Water present in these channels gradually coalesces and fills the hexagonal pores sequentially after the channels are completely hydrated. The development of hydrogen-bond networks inside the MOF pores as a function of relative humidity is characterized at the molecular level using experimental and computational infrared spectroscopy. A detailed analysis of the OH-stretch vibrational band indicates that the low-frequency tail stems from strongly polarized hydrogen-bonded water molecules, suggesting the presence of some structural disorder in the experimental samples. Strategies for designing efficient water harvesting MOFs are also proposed based on the mechanism of water adsorption in NU-1500-Cr.

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

Light-Mediated Synthesis and Reprocessing of Dynamic Bottlebrush Elastomers under Ambient Conditions

We introduce a novel grafting-through polymerization strategy to synthesize dynamic bottlebrush polymers and elastomers in one step using light to construct a disulfide-containing backbone. The key starting material—α-lipoic acid (LA)—is commercially available, inexpensive, and biocompatible. When installed on the chain end(s) of poly(dimethylsiloxane) (PDMS), the cyclic disulfide unit derived from LA polymerizes under ultraviolet (UV) light in ambient conditions. Significantly, no additives such as initiator, solvent, or catalyst are required for efficient gelation. Formulations that include bis-LA-functionalized cross-linker yield bottlebrush elastomers with high gel fractions (83–98%) and tunable, supersoft shear moduli in the ~20–200 kPa range. An added advantage of these materials is the dynamic disulfide bonds along each bottlebrush backbone, which allow for light-mediated self-healing and on-demand chemical degradation. These results highlight the potential of simple and scalable synthetic routes to generate unique bottlebrush polymers and elastomers based on PDMS.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Spin-Dependent Stereochemistry: A Nonadiabatic Quantum Dynamics Case Study of S + H 2 → SH + H Reaction

We study the spin-dependent stereodynamics of the S + H 2 → SH + H reaction by using full-dimensional quantum dynamics calculations with zero total nuclear angular momentum along the triplet 3 A" states and singlet 1 A' states. We find that the interplay between the electronic spin direction and the molecular geometry has a measurable influence on the singlet–triplet intersystem crossing reaction probabilities. Our results show that for some incident scattering angles in the body-fixed frame, the relative difference in intersystem crossing reaction probabilities (as determined between spin up and spin down initial states) can be as large as 15%. Furthermore, our findings are an ab initio demonstration of spin-dependent nonadiabatic dynamics, which we hope will shine light as far as understanding the chiral-induced spin selectivity effect.

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

Rapid assessment of structural and compositional changes during early stages of zirconium alloy oxidation

A multimodal chemical imaging approach has been developed and applied to detail the dynamic, atomic-scale changes associated with oxidation of a zirconium alloy (Zircaloy-4). Scanning transmission electron microscopy, a gas-phase reactor chamber attached to an atom probe tomography instrument, and synchrotron-based X-ray absorption near-edge spectroscopy were employed to reveal morphology, composition, crystal, and electronic structure changes that occur during initial stages of oxidation at 300 °C. Oxidation was carried out in 10 mbar O 2 gas for short exposure times of 1 and 5 min. A multilayered oxide film with a cubic ZrO adjacent to the oxide/metal interface, a nanoscopic transition region with a graded composition of ZrO 2-x (where 0 < x < 1), and tetragonal ZrO 2 in the outermost oxide were formed. Partitioning of the major alloying element (tin) to the oxide/metal interface and heterogeneously within the oxide accompanied the development of the layered oxide. Our work provides a rapid, high-throughput approach for detailed characterisation of initial stages of zirconium alloy oxidation at an accelerated time scale, with implications for several other alloy systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Subnanosecond phase transition dynamics in laser-shocked iron

Iron is one of the most studied chemical elements due to its sociotechnological and planetary importance; hence, understanding its structural transition dynamics is of vital interest. By combining a short pulse optical laser and an ultrashort free electron laser pulse, we have observed the subnanosecond structural dynamics of iron from high-quality x-ray diffraction data measured at 50-ps intervals up to 2500 ps. We unequivocally identify a three-wave structure during the initial compression and a two-wave structure during the decaying shock, involving all of the known structural types of iron (α-, γ-, and ε-phase). In the final stage, negative lattice pressures are generated by the propagation of rarefaction waves, leading to the formation of expanded phases and the recovery of γ-phase. Our observations demonstrate the unique capability of measuring the atomistic evolution during the entire lattice compression and release processes at unprecedented time and strain rate.

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

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗