Search NASA⌕ Search

SEARCH · Search NASA

Results for “scaling reactions”

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

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗
DOE OSTI2023

Computational Discovery of Active and Selective Metal‐Nitrogen‐Graphene Catalysts for Electrooxidation of Water to H 2 O 2

A direct electrosynthesis of H 2 O 2 from either O 2 or H 2 O is an attractive strategy to replace the energy-intensive industrial anthraquinone process. Two-electron water oxidation reaction (2e-WOR) offers several advantages over the oxygen reduction reaction such as better mass transfer due to the absence of gas-phase reactants. However, 2e-WOR is a more challenging and less studied process with only a handful of metal oxides exhibiting reasonable activity/selectivity properties. Herein, we employ density-functional-theory calculations to screen a variety of metal-nitrogen-graphene structures for 2e-WOR. As a consequence of scaling between the adsorption energies of reaction intermediates, we determine a linear relation between selectivities for the first and second reaction steps of 2e-WOR, viz. that if selectivity toward adsorbed OH is improved, then selectivity toward H 2 O 2 at the subsequent step is decreased. We also find that selectivity and activity are linearly scaled in such a way that a higher activity (i.e., a lower overpotential) leads to a lower selectivity for the H 2 O 2 formation step. Based on the obtained results several chemistries, e.g., containing NiN x -C moieties, are predicted to rival the best-performing metal oxides such as ZnO and CaSnO 3 in terms of combination of their activity/selectivity characteristics for 2e-WOR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

A numerical procedure for analysis of finite rate reacting flows

Combustion processes in rocket propulsion systems are characterized by the existence of multiple, vastly differing time and length scales, as well as flow-speeds at wide variation of Mach numbers. The chemical kinetics processes in the highly active reaction zone are characterized by much smaller scales compared to fluid convective and diffusive time scales. An operator splitting procedure for transient finite rate chemistry problems has been developed using a pressure based method, which can be applied to all speed flows without difficulties. The splitting of chemical kinetics terms formed the fluid-mechanical terms of the species equation ameliorated the difficulties associated with the disparate time scales and stiffness in the set of equations which describes highly exothermic combustion. A combined efficient ordinary differential equations (ODE) solver was used to integrate the effective chemical source terms over the residence time at each grid cell. One and two dimensional reacting flow situations were carried out to demonstrate and verify the current procedure. Different chemical kinetics with different degrees of nonlinearity have also been incorporated to test the robustness and generality of the proposed method.

Shang, H. M.↗
NASA NTRS2020

Molecular Simulations of Ultrafast Radiation Induced Melting at Metal-Semiconductor Interfaces

Understanding radiation induced ultrafast melting at material interfaces is essential in designing robust electronic devices for aviation/space applications and in laser machining. While it is difficult to achieve the spatial and temporal resolution required to quantify the phenomenon experimentally, simulations can provide the detailed mechanisms of the structural changes that happen during phase transition. In this work, we use molecular simulations to study the effect of radiation damage on silicon carbide (SiC) - tungsten (W) interfaces which is of interest in high power electronics. A multi-scale approach is involved wherein the reactions at the interfaces are quantified using ab-initio molecular dynamics (MD) simulations and classical MD simulations are employed to understand the structural and diffusional changes across the material interface. Finally, coarse-grained Lennard-Jones type models are used to study the larger scale mechanisms and structures obtained due to the induced damages. We show that the response of the material to radiation damage depends on factors such as energy of the incident radiation, thermal properties, and molecular structure of the material.

Ravichandran, Ashwin↗
DOE OSTI2022

Support Effect and Surface Reconstruction in In 2 O 3 / m- ZrO 2 Catalyzed CO 2 Hydrogenation

Here, we investigate the chemical and structural dynamics at the interface of In 2 O 3 /m-ZrO 2 and their consequences on the CO 2 hydrogenation reaction (CO 2 HR) under reaction conditions. While acting to enrich CO 2 , monoclinic zirconia (m-ZrO 2 ) was also found to serve as a chemical and structural modifier of In 2 O 3 that directly governs the outcome of the CO 2 HR. These modifying effects include the following: (1) Under reaction conditions (above 623 K), partially reduced In 2 O 3 , i.e., InO x (0 < x < 1.5), was found to migrate in and out of the subsurface of m-ZrO 2 in a semireversible manner, where m-ZrO 2 accommodates and stabilizes InO x by serving as a reservoir. The decreased concentration of surface InO x under elevated temperatures coincides with significantly decreased selectivity toward methanol and a sharp increase of the reverse water–gas shift reaction. The reconstruction-induced variation of InO x concentration appears to be one of the most important factors contributing to the altered catalytic performance of CO 2 HR at different reaction conditions. (2) The strong interactions and reactions between m-ZrO 2 and In 2 O 3 result in the activation of a pool of In–O bonds at the In 2 O 3 /m-ZrO 2 interface to form oxygen vacancies. On the other hand, the high dispersity of In 2 O 3 nanostructures onto m-ZrO 2 prevents their over-reduction under catalytically relevant conditions (up to 673 K), when bare In 2 O 3 is unavoidably reduced into the metallic phase (In 0 ). The relationship between the extent of reduction of In 2 O 3 and catalytic performance (CO 2 conversion, CH 3 OH selectivity, or yield of CH 3 OH) suggests the presence of an optimum coverage of surface InO x and oxygen vacancies under reaction conditions. The conventional model that links catalytic performance solely to the coverage of oxygen vacancies appears invalid in the present case. In situ analysis also allows the observation of surface reaction intermediates and their interconversions, including the reduction of CO 3 * into formate, a precursor for the formation of methanol and CO. The combinative ex situ and in situ study sheds light on the reaction mechanism of the CO 2 HR on In 2 O 3 /m-ZrO 2 -based catalysts. Our findings on the large-scale surface reconstructions, support effect, and the reaction mechanism of In 2 O 3 /m-ZrO 2 for CO 2 HR may apply to other related metal oxide catalyzed CO 2 reduction reactions.

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

A new phenomenological model to describe root-soil interactions based on percolation theory

In his paper on net primary productivity of terrestrial communities predicted from climatological data, Rosenzweig (1968) argued that variability in productivity is well accounted for by (evapo)-transpiration, and that water from transpiration is, on global scales, the most variable component in the photosynthesis reaction. The goal of this paper is to investigate whether variability in plant growth on local scales and within species is primarily related to transpiration under several scenarios including different terrain curvature, slope aspect, soil characteristics, and climate ranges. Here, we test the hypothesis that this relationship exists because root growth into the surface soil layers (0–2 m) tends to follow paths with minima in resistance, which in turn maximizes water flow and nutrient delivery rates that regulate growth. The set of all connected paths with individual pore-to-pore flow resistances less than a critical, percolating, value forms a cluster with mass fractal dimensionality, d f . We propose that roots follow paths through the 2D percolation cluster, defining the set of all optimal flow paths, making the 2D value of d f from percolation relevant to root fractal dimensionality. The tortuosity of such optimal paths as defined in percolation theory should then relate root length to root radial extent, linking the parameters of root tortuosity and plant productivity. Our analysis of large data sets across species implies that root radial extent and tree height are both proportional to cumulative transpiration until trees approached maximum height, and their growth rates are proportional to the transpiration rate, not to the moisture content. Local variations in tree height as functions of the variables investigated appear generally consistent with deduced variations in transpiration. Here this correlation is investigated more closely in the context of studies addressing individual tree species.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Atomic-scale 3D structural dynamics and functional degradation of Pt alloy nanocatalysts during the oxygen reduction reaction

Pt-based electrocatalysts are the primary choice for fuel cells due to their superior oxygen reduction reaction (ORR) activity. To enhance ORR performance and durability, extensive studies have investigated transition metal alloying, doping, and shape control to optimize the three key governing factors for ORR: geometry, local chemistry, and strain of their surface and subsurface. However, systematic optimization remains incomplete, as it requires an atomic-scale understanding of these factors and their dynamics over potential cycling, as well as their relationship to ORR activity. Here, we implement neural network-assisted atomic electron tomography to measure the 3D atomic structural dynamics and their effects on the functional degradation of PtNi alloy catalysts. Our results reveal that PtNi catalysts undergo shape changes, surface alloying, and strain relaxation during cycling, which can be effectively mitigated by Ga doping. By combining geometry, local chemistry, and strain analysis, we calculated the changes in ORR activity over thousands of cycles and observed that Ga doping leads to higher initial activity and greater stability. These findings offer a pathway to understanding 3D atomic structural dynamics and their relation to ORR activity during cycling, paving the way for the systematic design of durable, high-efficiency nanocatalysts.

Jeong, Chaehwa↗
DOE OSTI2022

Molecular-Modified Photocathodes for Applications in Artificial Photosynthesis and Solar-to-Fuel Technologies

Nature offers inspiration for developing technologies that integrate the capture, conversion, and storage of solar energy. In this review article, we highlight principles of natural photosynthesis and artificial photosynthesis, drawing comparisons between solar energy transduction in biology and emerging solar-to-fuel technologies. Key features of the biological approach include use of earth-abundant elements and molecular interfaces for driving photoinduced charge separation reactions that power chemical transformations at global scales. Finally, for the artificial systems described in this review, emphasis is placed on advancements involving hybrid photocathodes that power fuel-forming reactions using molecular catalysts interfaced with visible-light-absorbing semiconductors.

14 SOLAR ENERGY↗
DOE OSTI2021

Advancing Catalytic Fast Pyrolysis Through Integrated Experimentation and Multi-Scale Computational Modeling

This webinar will highlight recent results from a multi-disciplinary research effort in which integrated reaction testing was coupled with particle- and reactor-scale computational modeling to advance catalytic fast pyrolysis (CFP) for the production of renewable hydrocarbon fuels. Data will be presented from a series of ex situ CFP experiments in which a fixed bed of Pt/TiO2 was utilized with co-fed H2 to upgrade woody biomass pyrolysis vapors. Further discussion will include the application of these data towards the development of (1) a multiscale simulation framework to de-couple apparent kinetics from both intraparticle and reactor-scale transport phenomena and (2) a finite element computational model to understand and predict thermal excursions during catalyst regeneration. Throughout the presentation, the speakers will emphasize synergistic outcomes derived from the collaborative approach and highlight ongoing research efforts to accelerate technology maturation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRSDate not supplied

Simulation of etch pit formation through active sites in carbon fiber micro-structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation as a result of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gasses such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various reaction mechanisms such as adsorption, desorption, Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗
NASA NTRSDate not supplied

Simulation of Etch Pit Formation in DSMC Through Active Sites in Carbon Fiber Micro-Structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation because of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gases such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various detailed surface reaction mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗
NASA NTRSDate not supplied

DSMC Simulation of Etch Pit Formation Through Active Sites in Carbon Fiber Micro-Structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation because of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gases such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various detailed surface reaction mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗
DOE OSTI2024

Inferring fusion nuclear burnwidths with low gain photomultiplier impulse response functions

When an inertial confinement fusion implosion is compressed, it maintains thermonuclear density and temperatures for a very short time scale, about 100 ps. The Gamma Reaction History diagnostic measures the time evolution of the fusion burn, but its temporal resolution is limited by the use of a photomultiplier tube (PMT) to amplify the photon signal. Multichannel plate-based PMTs have a fast (~120 ps) full-width at half-max impulse response function (IRF), but the time scale is similar to the incoming physics signal. An analysis routine is used to remove the effect of the PMT IRF and infer the incident fusion burnwidth. With the National Ignition Facility achieving ignition and creating much brighter signals, the PMTs are run at gains three orders of magnitude lower than nominal operation. Calibration at these settings shows the PMT IRFs get ~15% wider. Taking the gain-dependent IRF can affect the inferred nuclear burnwidths by up to ~15%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ruling out Color Transparency in Quasielastic 12 C(e,e'p) up to Q 2 of 14.2 (GeV/c) 2

Quasielastic $^{12}$C$(e,e'p)$ scattering was measured at space-like 4-momentum transfer squared $Q^2$ = 8, 9.4, 11.4, and 14.2 (GeV/c)$^2$, the highest ever achieved to date. Nuclear transparency for this reaction was extracted by comparing the measured yield to that expected from a plane-wave impulse approximation calculation without any final state interactions. Additionally, the measured transparency was consistent with no $Q^2$ dependence, up to proton momenta of 8.5 GeV/c, ruling out the quantum chromodynamics effect of color transparency at the measured $Q^2$ scales in exclusive $(e,e'p)$ reactions. Furthermore, these results impose strict constraints on models of color transparency for protons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRSDate not supplied

Building and Breaking Carbon Composites with REACTER

Carbon composites have become indispensable for aerospace and other high-performance applications, and a detailed picture of their morphology and failure mechanisms remains difficult to obtain through experiment. REACTER is a versatile computational modeling tool for atomistic molecular dynamics designed to model chemical reactions at the speed and length scales of classical force fields.1 In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites, carbon-fiber reinforced polymers (CFRP) and carbon nanotube (CNT) composites. Carbon fiber core morphologies were created by the method of Desai et al.,2 but using the advanced reaction constraints framework of REACTER, their proposed multistep procedure was reduced to a single uninterrupted molecular dynamics simulation. The carbon fiber filler was embedded into a polymer matrix by simulated in situ polymerization of several thermosetting resins, including bismaleimide and polyarylacetylene, to obtain the final CFRP model. To generate the second class of carbon composite, CNT networks were grown dynamically using the new ‘create atoms’ feature of REACTER, and similarly infiltrated with resin to obtain CNT composites. The resulting models were compared directly to experiment using simulated high resolution transmission electron microscopy and x-ray diffraction. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a reaction constraint on the total potential energy of the involved atoms.

Molecular Dynamics↗
DOE OSTI2021

Scaling of Transition State Vibrational Frequencies and Application of d -Band Theory to the Brønsted–Evans–Polanyi Relationship on Surfaces

Semiempirical energy relations provide a means of estimating thermodynamic properties. Specifically, linear scaling relationships (LSRs) and Brønsted–Evans–Polanyi (BEP) relationships correlate adsorption energies between adsorbates across surfaces and reaction energies with activation barriers, respectively. Although vibrational scaling relations (VSRs) exist between adsorbates at identical sites, scaling between vibrational frequencies of adsorbed local minima and transition states is lacking. Here, we present density functional theory calculations for AH X (A = C, N, O) diffusions on transition metal surfaces and reveal linear scaling between frequencies of local minima and the transition state between those minima. Using d-band theory and linear muffin tin orbital theory (LMTO), we derive the slopes of these transition state vibrational scaling relations (TSVSRs) and, in so doing, provide a rigorous theory extending the original BEP relations developed for solution chemistry to surface chemistry. Furthermore, with a single reference DFT calculation, we predict the slopes and quantify uncertainty in the predictions.

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

Assessment of enzyme active site positioning and tests of catalytic mechanisms through X-ray–derived conformational ensembles

How enzymes achieve their enormous rate enhancements remains a central question in biology, and our understanding to date has impacted drug development, influenced enzyme design, and deepened our appreciation of evolutionary processes. While enzymes position catalytic and reactant groups in active sites, physics requires that atoms undergo constant motion. Numerous proposals have invoked positioning or motions as central for enzyme function, but a scarcity of experimental data has limited our understanding of positioning and motion, their relative importance, and their changes through the enzyme’s reaction cycle. To examine positioning and motions and test catalytic proposals, we collected “room temperature” X-ray crystallography data for Pseudomonas putida ketosteroid isomerase (KSI), and we obtained conformational ensembles for this and a homologous KSI from multiple PDB crystal structures. Here, ensemble analyses indicated limited change through KSI’s reaction cycle. Active site positioning was on the 1- to 1.5-Å scale, and was not exceptional compared to noncatalytic groups. The KSI ensembles provided evidence against catalytic proposals invoking oxyanion hole geometric discrimination between the ground state and transition state or highly precise general base positioning. Instead, increasing or decreasing positioning of KSI’s general base reduced catalysis, suggesting optimized Ångstrom-scale conformational heterogeneity that allows KSI to efficiently catalyze multiple reaction steps. Ensemble analyses of surrounding groups for WT and mutant KSIs provided insights into the forces and interactions that allow and limit active-site motions. Most generally, this ensemble perspective extends traditional structure–function relationships, providing the basis for a new era of “ensemble–function” interrogation of enzymes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Effect of Heterogeneities on the Reaction-zone of a Propagating Detonation Wave

Reactive burn models for propagating detonation waves in an explosive are based on the ZND (Zel'dovich-Von Neumann-Doring) theory. A key property of the ZND theory, which leads to self-sustaining detonation waves, is a sonic point relative to the front in the reaction zone. Moreover, burn models assume the explosive is a homogeneous material. This assumption is used for plastic-bonded explosives, even though they are heterogeneous materials. The heterogeneities are accounted for only with an empirical burn rate. Furthermore, the reaction-zone width can be less than the length scale of heterogeneities. This raises the question of the effect of heterogeneities on the reaction zone. In particular, whether the reaction zone can be steady and have a well defined sonic point. To examine this question we focus on PBX 9501 since there is data on the grain-binder heterogeneities and on the reaction-zone velocity time history of a steady planar detonation wave. Using pseudo-mesoscale 1-D simulations, we examine the effect of a detonation wave transversing explosive and binder segments of random length. The simulations show that the reaction zone is quasi-steady with only short wavelength small amplitude variations in the detonation wave pressure and speed. The main effect of the heterogeneities is to introduce noise in the pressure and velocity fields that propagates into the Taylor wave behind the detonation. This is consistent with the scatter in the reaction-zone measurements for PBX 9501.

42 ENGINEERING↗