Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reaction Rates”

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 343 records · Page 19

Entropic Control of HD Exchange Rates over Dilute Pd-in-Au Alloy Nanoparticle Catalysts

Dilute Pd-in-Au alloy catalysts are promising materials for selective hydrogenation catalysis. Extensive surface science studies have contributed mechanistic insight on the energetic aspect of hydrogen dissociation, migration and recombination on dilute alloy systems. Yet, translating these fundamental concepts to the kinetics and free energy of hydrogen dissociation on nanoparticle catalysts operating at ambient pressures and temperatures remains challenging. Here, the effect of the Pd concentration and Pd ensemble size on the catalytic activity, apparent activation energy and rate limiting process is addressed by combining experiment and theory. Experiments in a flow reactor show that a compositional change from 4 to 8 atm% Pd of the Pd-in-Au alloy catalyst leads to strong increase in activity, even exceeding the activity per Pd atom of monometallic Pd under the same conditions, albeit with an increase in apparent activation energy. First-principles calculations show that the rate and apparent activation enthalpy for HD exchange increase when increasing the Pd ensemble size from single Pd atoms to Pd trimers in a Au surface, suggesting that the ensemble size distribution shifts from mainly single Pd atoms on the 4 atm% Pd alloy to larger Pd ensembles of at least three atoms for the 8 atm% Pd/Au catalyst. The DFT studies also indicated that the rate-controlling process is different: H 2 (D 2 ) dissociation determines the rate for single atoms whereas recombination of adsorbed H and D determines the rate on Pd trimers. Both experiment and theory suggest that the increased reaction rate with increasing Pd content and ensemble size stems from an entropic driving force. Lastly, our results support hydrogen migration between Pd sites via Au and indicate that the dilute alloy design prevents the formation of subsurface hydrogen, which is crucial in achieving high selectivity in hydrogenation catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Steam-created grain boundaries for methane C–H activation in palladium catalysts

Defects may display high reactivity because the specific arrangement of atoms differs from crystalline surfaces. We demonstrate that high-temperature steam pretreatment of palladium catalysts provides a 12-fold increase in the mass-specific reaction rate for carbon-hydrogen (C–H) activation in methane oxidation compared with conventional pretreatments. Through a combination of experimental and theoretical methods, we demonstrate that an increase in the grain boundary density through crystal twinning is achieved during the steam pretreatment and oxidation and is responsible for the increased reactivity. The grain boundaries are highly stable during reaction and show specific rates at least two orders of magnitude higher than other sites on the palladium on alumina (Pd/Al 2 O 3 ) catalysts. Theoretical calculations show that strain introduced by the defective structure can enhance C–H bond activation. Introduction of grain boundaries through laser ablation led to further rate increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Proton branching ratios of 23 Mg levels

Background: The anomalous 22 Ne abundance measured in certain presolar graphite grains is thought to arise from the decay of 22 Na that was synthesized at high temperatures in core-collapse supernovae. To better interpret this abundance anomaly, the primary destruction mechanism of 22 Na, the 22 Na(p,γ) 23 Mg reaction, must be better understood. Purpose: Determine proton branching ratios of several 23 Mg excited states that play a role in the high-temperature 22 Na(p,γ) 23 Mg reaction rate. Methods: Particle decays of 23 Mg excited states populated with the previously reported 24 Mg(p,d) 23 Mg transfer reaction measurement [Kwag et al., Eur. Phys. J. A 56, 108 (2020)] were analyzed to extract proton branching ratios. The reaction was studied using a 31-MeV proton beam from the Holifield Radioactive Ion Beam Facility of Oak Ridge National Laboratory and 24 Mg solid targets. Results: Proton branching ratios of several 23 Mg excited states in the energy range of Ex = 8.044 - 9.642 MeV were experimentally determined for the first time for the p0 and p1'(p1+p2+p3) decay channels. Conclusions: These new branching ratios for 23 Mg levels can provide an experimental foundation for an improved high-temperature rate of the 22 Na(p,γ) 23 Mg reaction needed to understand production of anomalously high 22 Ne abundance in core-collapse supernovae.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multiple redox mechanisms for water-gas shift reaction on Fe 3 O 4 (1 1 1) surface: A density functional theory and mean-field microkinetic modeling study

In this work, we investigated redox mechanisms for the high temperature (HT) water–gas shift reaction (WGSR) on the Fe 3 O 4 (1 1 1) surface using density functional theory (DFT) calculations and mean-field microkinetic modeling (MKM). The redox pathways branched into three Langmuir-Hinshelwood processes (LH1, LH2 and LH3) and one Mars-van-Krevelen (MVK) process (in the presence of oxygen vacancy) from multidentate binding of CO and CO 2 over four available reactive sites (Fe oct2 , Fe tet1 , Fe bridge, and O1). We found that the LH1 and LH2 processes have CO adsorption at a single iron site (Fe oct2 or Fe tet1 ), while the LH3 and MVK processes have stronger chemisorption of CO or CO 2 by both the Fe oct2 site and the O 1 site. From the mechanistic study of these reaction path, we recognized that availability of O 1 sites was key to proceed to either CO oxidation by the single site (LH1, LH2) or by the dual site (LH3, MVK). We observed that the reaction energetics were significantly different in the CO oxidation steps, where the single or dual site results in markedly different apparent activation energy (73–281 kJ/mol) and reaction rates (10 -3 –10 -10 mol∙m -2 ∙s -1 ) among the four reaction mechanisms. The utilization of mean field MKM with DFT reaction energetics helps explain the experimental debates for the catalytic reaction details as well as providing a possible direction to engineer the catalyst with higher activity.

42 ENGINEERING↗

Operability thresholds for thermally damaged EBW detonators

Operability thresholds that differentiate between functional RP-87 exploding bridge wire (EBW) detonators and nonfunctional RP-87 EBW detonators (duds) were determined by measuring the time delay between initiation and early wall movement (function time). The detonators were inserted into an externally heated hollow cylinder of aluminum and fired with current flow from a charged capacitor using an exploding bridge wire (EBW initiated). Functioning detonators responded like unheated pristine detonators when the function time was 4 μs or less. The operability thresholds of the detonators were characterized with a simple decomposition cookoff model calibrated using a modified version of the Sandia Instrumented Thermal Ignition (SITI) experiment. These thresholds are based on the calculated state of the PETN when the detonators fire. The operability threshold is proportional to the positive temperature difference (ΔT) between the maximum temperature within the PETN and the onset of decomposition (~406 K). Furthermore, the temperature difference alone was not sufficient to define the operability threshold. The operability threshold was also proportional to the time that the PETN had been at elevated temperatures. That is, failure was proportional to both temperature and reaction rate. The reacted gas fraction is used in the current work for the reaction correlation. Melting of PETN also had a significant effect on the operability threshold. Detonator failure occurred when the maximum temperature exceeded the nominal melting point of PETN (414 K) for 45±5 s or more.

42 ENGINEERING↗

Constraining (n,γ) cross sections via surrogate measurements with hyperion

Indirect measurement techniques, such as the surrogate method, are necessary to constrain cross sections and reaction rates that are inaccessible for direct measurement. The surrogate method uses an alternate reaction channel to populate a nucleus of interest, and combines resulting experimental data with theory to constrain the desired (n,γ) cross section. Individual γ-decay probabilities are experimentally extracted using γ-ray spectroscopy. This project utilizes surrogate analysis of 64,70 Zn(p,d) and 64,70 Zn(p,t) reactions measured with the Hyperion array at Texas A&M University's Cyclotron Facility to probe the relevant neutron capture reactions 61,62,67,68 Zn(n,γ). As an example, preliminary γ-decay probabilities are presented for 69 Zn.

Neutron capture↗

Nuclear spin conservation enables state-to-state control of ultracold molecular reactions

Quantum-state control of reactive systems has enabled microscopic probes of underlying interaction potentials and the alteration of reaction rates using quantum statistics. However, extending such control to the quantum states of reaction outcomes remains challenging. Here in this work, we realize this goal by utilizing the conservation of nuclear spins throughout the reaction. Using resonance-enhanced multiphoton ionization spectroscopy to investigate the products formed in bimolecular reactions between ultracold KRb molecules we find that the system retains a near-perfect memory of the reactants’ nuclear spins, manifested as a strong parity preference for the rotational states of the products. We leverage this effect to alter the occupation of these product states by changing the coherent superposition of initial nuclear spin states with an external magnetic field. In this way, we are able to control both the inputs and outputs of a reaction with quantum-state resolution. The techniques demonstrated here open up the possibilities to study quantum entanglement between reaction products, and ultracold reaction dynamics at the state-to-state level.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Lasers for the observation of multiple order nuclear reactions

Nuclear reaction rates become nonlinear with respect to flux (cm −2 s −1 ) in extreme environments such as those found during stellar nucleosynthesis and terrestrial nuclear detonations. To observe these effects directly in the laboratory, extremely high particle fluences (cm −2 ) are necessary but not sufficient. Reactor-based neutron sources, such as the Institut Laue-Langevin’s high-flux neutron reactor, were previously the closest to meeting this challenge, albeit over ∼hour time scales. In ultra-high flux environments, where multiple reactions occur on picosecond time scales, nuclei are unable to return to their ground states between reactions; consequently, reactions take place on excited nuclei. To accurately model high-flux environments, data on the cross-sections of excited nuclear states are required, which differ significantly from those of ground states due to spin/parity effects. In order to replicate these effects in the laboratory, short high-fluence pulses on the order of the lifetime of a typical nuclear excited state (generally ≲1 ns) are required. Particle beams generated by high-intensity lasers are uniquely positioned to meet this need with the potential to produce fluences of 10 17 protons/cm 2 and 10 22 neutrons/cm 2 over a few pico-seconds or less. In addition to providing a quantitative analysis of the rates of multiple rapid reactions in general, the present work examines a number of laser-based experiments that could be conducted in the near future to observe multiple rapid reactions for laboratory-based astrophysics and the measurement of exotic cross-sections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Effects of oligomerization and decomposition on the nanoparticle growth: a model study

Abstract. The rate at which freshly formed secondary aerosol particles grow is an important factor in determining their climate impacts. The growth rate of atmospheric nanoparticles may be affected by particle-phase oligomerization and decomposition of condensing organic molecules. We used the Model for Oligomerization and Decomposition in Nanoparticle Growth (MODNAG) to investigate the potential atmospheric significance of these effects. This was done by conducting multiple simulations with varying reaction-related parameters (volatilities of the involved compounds and reaction rates) using both artificial and ambient measured gas-phase concentrations of organic vapors to define the condensing vapors. While our study does not aim at providing information on any specific reaction, our results indicate that particle-phase reactions have significant potential to affect the nanoparticle growth. In simulations in which one-third of a volatility basis set bin was allowed to go through particle-phase reactions, the maximum increase in growth rates was 71 % and the decrease 26 % compared to the base case in which no particle-phase reactions were assumed to take place. These results highlight the importance of investigating and increasing our understanding of particle-phase reactions.

54 ENVIRONMENTAL SCIENCES↗

Bridging Scales to Model Reactive Diffusive Transport in Porous Media

Published by ECS. Two novel scale-bridging algorithms to model reaction-diffusion transport in porous media are presented. The algorithms are based on direct numerical simulations and couple the information of a micro-scale model, which accounts for the large field of view provided by micro X-ray computed tomography (X-ray CT), and a nano-scale model, which locally resolves transport in the fine structure extracted from nano X-ray CT. The micro-scale model is discretized in the through-plane direction into a 1D grid, where effective properties and internal boundaries are determined based on the results from the nano-scale model. The validated algorithms are used to examine transport of oxygen in precious group metal-free electrodes considering both zero- and first-order kinetics. Unlike conventional methods, the results show that the effective diffusivity is not a passive property but increases in regions where the reaction-rate coefficient is large. The proposed algorithms account for the multiscale coupling of reaction-diffusion transport and material microstructure, thus improving the predictions compared to conventional methods.

42 ENGINEERING↗

Improved Performance of Cu(InGa)(SeS) 2 PV Modules Using the Reaction of Metal Precursors. Final Report

This project “Improved Performance of Cu(InGa)(SeS) 2 PV Modules using the Reaction of Metal Precursors” was a partnership led by the Institute of Energy Conversion (IEC) at the University of Delaware with Columbia University and the Molecular Foundry at the Lawrence Berkeley National Laboratory. The aim was to develop pathways to improve Cu(InGa)(SeS) 2 (CIGSS) thin film photovoltaic modules using processes compatible with low manufacturing cost. The CIGSS approach investigated was a two-step process including deposition of metal precursor films following by reaction in hydride gases utilizing IEC’s novel reactor. The process was similar to that under commercial development by the project’s industry partner Stion. When Stion went out of business mid-project the focus changed to a rapid thermal process considered more commercially viable. Approaches to improve the performance of solar cells using the reacted films focused on two material innovations. First, the overall Ga content was increased to increase the operating voltage, which is desirable for scale-up to commercial modules. Second, the processing and performance advantages arising from Ag alloying were investigated. Advanced characterization guided process and material development including control of relative composition gradients. Research on the formation of Cu-Ga-In metal precursors utilized sputtering deposition which is normally used in commercial applications. The work resulted in processes for deposition of precursor stacks with increased relative Ga content and effects of deposition parameters on morphology and phase composition were established. It was shown that the metal precursor films have comparable phase composition and morphology so subsequent reaction follows from the same starting point. The addition of Ag to the metal precursors gave more uniform morphology and improved adhesion of reacted films which enable higher reaction temperature for faster processing. A novel outcome was the discovery of a previously undocumented material phase in sputter-deposited and evaporated Ag-Cu-In-Ga thin films. Hydride gas reaction processes including time-temperature-concentration profiles were developed for different precursor compositions. This enables control of composition profiles to engineer through-film gradients for solar cell optimization with characterization and simulations used to correlate measured film composition profiles to measurements of devices. In particular, the gradient of sulfur at the front of the CIGSS film was found to be critical. The simulations guided process development leading to improved reproducibility of devices improved performance with higher Ga content and higher voltage. With Ag-alloyed precursors, the reaction pathways leading were determined. A significant finding was that Ag-alloying increases the reaction rate to completely convert precursor films to the final chalcopyrite which could enable reduced reaction time to benefit manufacturability. To maintain potential commercial viability, the process under investigation was refocused to a rapid thermal process that could potentially be incorporated into an in-line process for manufacturing. Precursors with different composition were capped with an extra selenium layer and reacted in hydrogen sulfide 5-15 minutes, compared to typically 2 hours in the previous multi-step batch process. Critical RTP parameters were identified to control the reaction. Further optimization would be needed for high efficiency solar cells but pathways to high quality devices with further optimization and improved heating uniformity were developed. The project also developed new optoelectronic characterization approaches with a focus on development and application of spatial- and time-resolved photoluminescence and a custom mapping photoluminescence microscope built. It was shown how critical electronic transport properties strongly depend on the chemical composition of the material and that a wide range of samples show inhomogeneity on a length scale larger than the grains in the films. Additionally, two-photon excitation capability was developed to distinguish bulk vs surface losses. The project advances the state-of-the -art for precursor reaction processes in several ways that could impact manufacturing. This includes validation of approaches to increase voltage and establishment of model-guided control to form optimal composition profiles. The application of process control approaches with knowledge of phase formation and reaction pathways can be critically valuable in designing a large-scale process.

14 SOLAR ENERGY↗

Modeling of Thermal Decomposition of TATB-Based Explosive for Safety Analysis

We investigate and model the cook-off behavior of LX-17 to understand the response of explosive systems in abnormal thermal environments. Decomposition has been explored via conventional ODTX (One-Dimensional Time-to-eXplosion), PODTX (ODTX with pressure-measurement), TGA (Thermo-Gravimetric Analysis), and DSC (Differential Scanning Calorimetry) experiments under isothermal and ramped temperature profiles. The data were used to fit reaction rate parameters for proposed schemes in an ALE3D computational model. This model includes chemical reactions, thermo- and hydro-dynamics, and material properties, including thermal expansion, compressibility, and strength. These parameterizations were carried out utilizing a Python evolutionary optimization method on LLNL’s high-performance computing clusters. Additional experiments are being developed to further characterize and monitor decomposition intermediates to improve the model. Once experimentally validated, this model will be scalable to several applications involving LX-17. Furthermore, the optimization methodology developed herein should be applicable to other high explosive materials.

Chemistry - Chemical explosives↗

Reactive Transport Modeling of Anthropogenic Carbon Mineralization in Stacked Columbia River Basalt Reservoirs

Numerical simulation of CO2 storage in basalts and related reactive lithologies requires modeling complex, coupled hydrologic and chemical processes, including multi-phase flow and transport, partitioning of CO2 into the aqueous phase, and chemical interactions with aqueous fluids and rock minerals. We conducted reactive transport simulations of the Wallula pilot-scale CO2 injection into the flow tops of the Grande Ronde Basalt using our PNNL STOMP-CO2 simulator with the ECKEChem reactive module. Our mineralization simulation of the ~1,000 tons of injected CO2 into the interflow zones was based on the hydrologic transport model we previously developed. For this work, the simulations considered geochemical reactions involving the basalt components, precipitates, formation brine, and injected CO2. In our benchmark case, carbonate minerals precipitated, resulting in ~20% of the CO2 being mineralized in 10 years. Increasing the reaction rate of a single primary mineral phase (clinopyroxene) by an order of magnitude resulted in a carbon mineralization reaction extent of ~90% over the same time interval. Based on these initial sensitivity analysis results, it is clear that a thorough understanding of primary mineral dissolution rates is required for accurately predicting long-term fate and transport of injected CO2 into basalt formations. Our reactive transport numerical simulations will be key components of commercial-scale CO2 storage operation permitting, de-risking, and optimization in mafic and ultramafic reservoirs.

Cao, Ruoshi↗

Hydrothermal corrosion behavior of CVD SiC in high temperature water

Here, the hydrothermal corrosion of polished and as-cut high purity chemical vapor deposited (CVD) SiC was studied in a constantly refreshing water loop. Light water reactor (LWR) conditions were simulated at 288, 320, and 350 °C with dissolved gas concentrations between 0.15 and 3 ppm H 2 or between 1 and 4 ppm O 2 . In hydrogenated water, the rate of material loss was low, calculated to be ~1.3 μm of recession after 5 years of service in 320 °C water. Moreover, there was no observed localized attack at any temperature. In oxygenated conditions, the corrosion rate was higher, with a calculated material loss >10 μm after 5 years of service in 1 ppm O 2 , 320 °C water. Mass loss significantly increased when grain fallout became significant (as early as 200h with 4 ppm O 2 at 350 °C or after 1000–2000h with 2 ppm O 2 at 288 °C). Grain fallout more than doubled the corrosion rate and a steady state corrosion rate in the grain fallout regime was not observed but expected to eventually occur once large grains begin to be removed. Polished specimens had lower mass loss than unpolished coupons. A kinetic analysis of the data in this work suggests that the corrosion rates are controlled by a single activation step in both oxygenated and deoxygenated conditions, with the reaction order with respect to oxygen being 1. A resulting reaction rate equation to predict corrosion of SiC (in mg/cm 2 s) in high purity water from 288 to 350 °C and up to 4 ppm O 2 was constructed: Rate = $\frac{0.1458}{1+SA}$ T (1.09(1–10 -3 T)[O 2 ]e– $\frac{1.275x10^4}{T}$ + 7.91x10 -6 e– $\frac{7.39x10^3}{T}$).

36 MATERIALS SCIENCE↗

Fission and the r -process nucleosynthesis of translead nuclei in neutron star mergers

In this report we study the impact of fission on the production and destruction of translead nuclei during the r-process nucleosynthesis occurring in neutron-star mergers. Abundance patterns and rates of nuclear energy production are obtained for different ejecta conditions using three sets of stellar reaction rates, one of which is based on microscopic and consistent calculations of nuclear masses, fission barriers, and collective inertias. We show that the accumulation of fissioning material during the r process can strongly affect the free neutron abundance after the r-process freeze-out. This leads to a significant impact on the abundances of heavy nuclei that undergo α decay or spontaneous fission, affecting the radioactive energy production by the ejecta at timescales relevant for kilonova emission.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A new way to unravel the $^{12}$C($\alpha $,$\gamma $)$^{16}$O cross section components using neural networks

The $^{12}$C($\alpha $,$\gamma $)$^{16}$O reaction rate is crucial in determining the carbon-to-oxygen abundance ratio in stellar nucleosynthesis. Measuring this reaction’s cross section at stellar energies is challenging due to its extremely small value, approximately 10 -17 barn at E c.m. = 300 keV. To address this, R-matrix calculations are employed to extrapolate data to lower energies, requiring a comprehensive understanding of each contribution to the cross section. The dominant contributions to the cross section at stellar energies arise from electric dipole (E1) and electric quadrupole (E2) transitions to the ground state of 16 O, along with a significant cascade contribution. Traditionally, these contributions have been separated using the γ-ray angular distribution. In this work, we propose a novel technique using the energy distribution of the 16 O recoils at the focal plane. This method involves a neural network trained on detailed Monte Carlo simulations of the energy distribution of recoils transported through the recoil mass separator ERNA. This approach enables the simultaneous determination of all three contributions with errors around 10% in the energy range E c.m. = 1.0–2.2 MeV. In conclusion, by employing this new technique, we aim to significantly improve the accuracy of determining the cross section of the $^{12}$C($\alpha $,$\gamma $)$^{16}$O reaction at astrophysical energies.

Astrochemistry↗

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗