Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic kinetics”

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 361 records · Page 20

Atomic structure of different surface terminations of polycrystalline ZnPd

The intermetallic compound ZnPd has been found to have desirable characteristics as a catalyst for the steam reforming of methanol. The understanding of the surface structure of ZnPd is important to optimize its catalytic behavior. However, due to the lack of bulk single-crystal samples and the complexity of characterizing surface properties in the available polycrystalline samples using common experimental techniques, all previous surface science studies of this compound have been performed on surface alloy samples formed through thin-film deposition. In this study, we present findings on the chemical and atomic structure of the surfaces of bulk polycrystalline ZnPd studied by a variety of complementary experimental techniques, including scanning tunneling microscopy (STM), x-ray photoelectron spectroscopy (XPS), low energy electron microscopy (LEEM), photoemission electron microscopy (PEEM), and microspot low-energy electron diffraction ( μ -LEED). These experimental techniques, combined with density functional theory (DFT)-based thermodynamic calculations of surface free energy and detachment kinetics at the step edges, confirm that surfaces terminated by atomic layers composed of both Zn and Pd atoms are more stable than those terminated by only Zn or Pd layers. DFT calculations also demonstrate that the primary contribution to the tunneling current arises from Pd atoms, in agreement with the STM results. The formation of intermetallics at surfaces may contribute to the superior catalyst properties of ZnPd over Zn or Pd elemental counterparts. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗

Radiation-Induced Oxidation Reactions of 2-Selenouracil in Aqueous Solutions: Comparison with Sulfur Analog of Uracil

One-electron oxidation of 2-selenouracil (2-SeU) by hydroxyl (●OH) and azide (●N3) radicals leads to various primary reactive intermediates. Their optical absorption spectra and kinetic characteristics were studied by pulse radiolysis with UV-vis spectrophotometric and conductivity detection and by the density functional theory (DFT) method. The transient absorption spectra recorded in the reactions of ●OH with 2-SeU are dominated by an absorption band with an λmax = 440 nm, the intensity of which depends on the concentration of 2-SeU and pH. Based on the combination of conductometric and DFT studies, the transient absorption band observed both at low and high concentrations of 2-SeU was assigned to the dimeric 2c-3e Se-Se-bonded radical in neutral form (2●). The dimeric radical (2●) is formed in the reaction of a selenyl-type radical (6●) with 2-SeU, and both radicals are in equilibrium with Keq = 1.3 × 104 M−1 at pH 4 (below the pKa of 2-SeU). Similar equilibrium with Keq = 4.4 × 103 M−1 was determined for pH 10 (above the pKa of 2-SeU), which admittedly involves the same radical (6●) but with a dimeric 2c-3e Se-Se bonded radical in anionic form (2●−). In turn, at the lowest concentration of 2-SeU (0.05 mM) and pH 10, the transient absorption spectrum is dominated by an absorption band with an λmax = 390 nm, which was assigned to the ●OH adduct to the double bond at C5 carbon atom (3●) based on DFT calculations. Similar spectral and kinetic features were also observed during the ●N3-induced oxidation of 2-SeU. In principle, our results mostly revealed similarities in one-electron oxidation pathways of 2-SeU and 2-thiouracil (2-TU). The major difference concerns the stability of dimeric radicals with a 2c-3e chalcogen-chalcogen bond in favor of 2-SeU.

2-selenouracil↗

Plasma-Coupled Flow Reactor Studies of Low-Temperature Plasma Assisted Kinetics of Methanol Blended with CO2

Ignition technologies based on low-temperature plasmas (LTP) have the potential to operate next generation engines at elevated pressures and increased dilution limits promoting higher efficiencies. From a practical standpoint, research on LTP igniters has shown to enhance combustion and ignition, improve flame stability, and extend the dilution limits of combustion. All of the aforementioned gains are complemented with higher ignition efficiencies. However, the biggest challenge with incorporating this technology into engines is the knowledge gap of how exactly plasma chemistry effects can enhance the basic combustion phenomena. This coupled with the lack of validated kinetic mechanisms for plasma-combustion chemistry is the biggest obstruction to recognize efficient ignition, especially for application relevant fuels and biofuels. In order to comprehensively evaluate the effects of LTP on an oxygenated fuel specific system, this present study examines the kinetics of methanol plasma-assisted pyrolysis and oxidation using a custom-built plasma flow reactor (PFR). Experimental regimes are also further extended to understand the effects of adding CO2 to the mixture and its consequence on reaction kinetics. The PFR is installed with a dielectric-barrier discharge (DBD) configuration to induce LTP into the fuel mixture. Non-equilibrium plasmas are generated by high-voltage pulses (nearing 20 KV) administered by a plasma pulser at high-pulse repetition rates (up to 10 kHz). In order to better understand and isolate plasma chemistry and its effect on neutral chemistry, the experiments were carried out by heavily diluting reactive mixtures in nitrogen at near isothermal conditions. This suppresses the effect of exothermic reactions on chemistry allowing stable intermediates and products to be detected and quantified using ex-situ GC/MS diagnostic methods. Experiments were carried out at 0.5 atm pressure and over a wide range of temperatures from 523 K to 1203 K. Experimental results depicted the enhancement in intermediates production as well as overall lower temperatures required for complete fuel consumption in the plasma specific cases as opposed to their pure thermal counterpart. Formation of oxygenated and nitrile compounds specific to plasma assisted pyrolysis cases illustrated the efficacy of LTP to introduce new reaction pathways accelerating fuel decomposition. Increase in reactivity at lower temperatures is sought to be an effect induced by plasma chemistry postulating the acceleration in intermediates production. The onset of thermal ignition in plasma assisted oxidation is seen 200 K earlier than thermal oxidation highlighting efficient fuel conversion to final byproducts. Enhanced collisional processes afforded by LTP is seen to perturb reaction pathways between oxygenated fuel radicals and N-atoms to alter overall chemical reactivity. New insights into kinetics are gained from this study highlighting governing plasma assisted combustion (PAC) pathways for oxygenated fuel reaction chemistry. The results from this study can be used to develop future mechanisms specific to plasma chemistry which can bridge the mechanistic knowledge gap for LTP ignition so that future engines can adopt LTPs in their design for efficient combustion.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Mechanism of Nitrogen-Carbon Bond Formation from Iron(IV) Disilylhydrazido Intermediates during N 2 Reduction

We recently reported a reaction sequence that activates C–H bonds in simple arenes as well as the N–N triple bond in N 2 , delivering the aryl group to N 2 to form a new N–C bond ( Nature 2020 , 584, 221). This enables the transformation of abundant feedstocks (arenes and N 2 ) into N-containing organic compounds. The key N–C bond forming step occurs upon partial silylation of N 2 . However, the pathway through which reduction, silylation, and migration occurred was unknown. Here, in this work, we describe synthetic, structural, magnetic, spectroscopic, kinetic, and computational studies that elucidate the steps of this transformation. N 2 must be silylated twice at the distal N atom before aryl migration can occur, and sequential silyl radical and silyl cation addition is a kinetically competent pathway to a formally iron(IV)–NN(SiMe 3 ) 2 intermediate that can be isolated at low temperature. Kinetic studies show its first-order conversion to the migrated product, and DFT calculations indicate a concerted transition state for migration. The electronic structure of the formally iron(IV) intermediate is examined using DFT and CASSCF calculations, which reveal contributions from iron(II) and iron(III) resonance forms with oxidized NNSi 2 ligands. The depletion of electron density from the Fe-coordinated N atom makes it electrophilic enough to accept the incoming aryl group. This new pathway for the N–C bond formation offers a method for functionalizing N 2 using organometallic chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laser flash photolysis studies of atmospheric free radical chemistry using optical diagnostic techniques

Some recent studies carried out in our laboratory are described where laser flash photolytic production of reactant free radicals has been combined with reactant and/or product detection using time-resolved optical techniques to investigate the kinetics and mechanisms of important atmospheric chemical reactions. Discussed are (1) a study of the radical-radical reaction O + BrO yields Br + O2 where two photolysis lasers are employed to prepare the reaction mixture and where the reactants O and BrO are monitored simultaneously using atomic resonance fluorescence to detect O and multipass UV absorption to detect BrO; (2) a study of the reaction of atomic chlorine with dimethylsulfide (CH3SCH3) where atomic resonance fluorescence detection of Cl is employed to elucidate the kinetics and tunable diode laser absorption spectroscopy is employed to investigate the HCl product yield; and (3) a study of the aqueous phase chemistry of Cl2(-) radicals where longpath UV absorption spectroscopy is employed to investigate the kinetics of the Cl2(-) + H2O reaction.

Wine, Paul H.↗

Best of both worlds: Enforcing detailed balance in machine learning models of transition rates

The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.

36 MATERIALS SCIENCE↗

Unraveling Nonequilibrium Generation of Atomic Nitrogen and Hydrogen in Plasma-Aided Ammonia Synthesis

In this study, nonequilibrium generation of atomic nitrogen and hydrogen governing ammonia production both in the gas phase and on the catalyst surface is critical to plasma-aided ammonia synthesis. Here, this work studies the nonequilibrium generation of atomic nitrogen and hydrogen with a focus on the kinetic role of vibrational energy transfer of hydrogen molecules in plasma-aided ammonia synthesis. By combining two-photon absorption laser-induced fluorescence measurements and plasma kinetic modeling, we found that plasma not only generates ammonia but also produces critical H, N, and NH radicals via both electron impact and vibrational energy transfer excitations. The vibrational energy transfer from the excited hydrogen H 2 (v = 1) to higher vibrational levels H 2 (v = 2–3) via the V–V exchange (H 2 (v)–H 2 (v)) and V–V' exchange (N 2 (v)–H 2 (v)) significantly enhances the H and NH production and then promotes the coupling between N and NH for ammonia synthesis both in the gas phase and on the catalyst surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enantioselective C–H Amination Catalyzed by Nickel Iminyl Complexes Supported by Anionic Bisoxazoline (BOX) Ligands

The trityl-substituted bisoxazoline ( TrH BOX) was prepared as a chiral analogue to a previously reported nickel dipyrrin system capable of ring-closing amination catalysis. Ligand metalation with divalent NiI 2 (py) 4 followed by potassium graphite reduction afforded the monovalent (TrH BOX)Ni(py) (4). Slow addition of 1.4 equiv of a benzene solution of 1-adamantylazide to 4 generated the tetrazido ( TrH BOX)Ni(κ 2 -N 4 Ad 2 ) (5) and terminal iminyl adduct ( TrH BOX)Ni(NAd) (6). Investigation of 6 via single-crystal X-ray crystallography, NMR and EPR spectroscopies, and computations revealed a Ni(II)-iminyl radical formulation, similar to its dipyrrinato congener. Complex 4 exhibits enantioselective intramolecular C–H bond amination to afford N-heterocyclic products from 4-aryl-2-methyl-2-azidopentanes. Catalytic C–H amination occurs under mild conditions (5 mol % catalyst, 60 °C) and provides pyrrolidine products in decent yield (29%–87%) with moderate ee (up to 73%). Substrates with a 3,5-dialkyl substitution on the 4-aryl position maximized the observed enantioselectivity. Kinetic studies to probe the reaction mechanism were conducted using 1H and 19F NMR spectroscopies. A small, intermolecular kinetic isotope effect (1.35 ± 0.03) suggests an H-atom abstraction step with an asymmetric transition state while the reaction rate is measured to be first order in catalyst and zeroth order in substrate concentrations. Enantiospecific deuterium labeling studies show that the enantioselectivity is dictated by both the H-atom abstraction and radical recombination steps due to the comparable rate between radical rotation and C–N bond formation. Furthermore, the competing elements of the two-step reaction where H-removal from the pro-R configuration is preferred while the preferential radical capture occurs with the Si face of the carboradical likely lead to the diminished ee observed, as corroborated by theoretical calculations. Furthermore, based on these enantio-determining steps, catalytic enantioselective synthesis of 2,5-bis-tertiary pyrrolidines is demonstrated with good yield (50–78%) and moderate ee (up to 79%).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An assessment of potential degradation products in the gas-phase reactions of alternative fluorocarbons in the troposphere

Tropospheric chemical transformations of alternative hydrofluorocarbons (HCF's) and hydrochlorofluorocarbons (HCFC's) are governed by hydroxyl radical initiated oxidation processes, which are likely to be analogous to those known for alkanes and chloroalkanes. A schematic diagram is used to illustrate plausible reaction mechanisms for their atmospheric degradation, where R, R', and R'' denote the F- and/or Cl-substituted alkyl groups derived from HCF's and HCFC's subsequent th the initial H atom abstraction by HO radicals. At present, virtually no kinetic data exist for the majority of these reactions, particularly for those involving RO. Potential degradation intermediates and final products include a large variety of fluorine- and/or chlorine-containing carbonyls, acids, peroxy acids, alcohols, hydrogen peroxides, nitrates and peroxy nitrates, as summarized in the attached table. Probably atmospheric lifetimes of these compounds were also estimated. For some carbonyl and nitrate products shown in this table, there seem to be no significant gas-phase removal mechanisms. Further chemical kinetics and photochemical data are needed to quantitatively assess the atmospheric fate of HCF's and HCFC's, and of the degradation products postulated in this report.

Niki, Hiromi↗

Structures Self-Assembled Through Directional Solidification

Nanotechnology has created a demand for new fabrication methods with an emphasis on simple, low-cost techniques. Directional solidification of eutectics (DSE) is an unconventional approach in comparison to low-temperature biomimetic approaches. A technical challenge for DSE is producing microstructural architectures on the nanometer scale. In both processes, the driving force is the minimization of Gibb's free energy. Selfassembly by biomimetic approaches depends on weak interaction forces between organic molecules to define the architectural structure. The architectural structure for solidification depends on strong chemical bonding between atoms. Constituents partition into atomic-level arrangements at the liquid-solid interface to form polyphase structures, and this atomic-level arrangement at the liquid-solid interface is controlled by atomic diffusion and total undercooling due to composition (diffusion), kinetics, and curvature of the boundary phases. Judicious selection of the materials system and control of the total undercooling are the keys to producing structures on the nanometer scale. The silicon-titanium silicide (Si-TiSi2) eutectic forms a rod structure under isothermal cooling conditions. At the NASA Glenn Research Center, directional solidification was employed along with a thermal gradient to promote uniform rods oriented with the thermal gradient. The preceding photomicrograph shows the typical transverse microstructure of a solidified Si-TiSi2 eutectic composition. The dark and light gray regions are Si and TiSi2, respectively. Preferred rod orientation along the thermal gradient was poor. The ordered TiSi2 rods have a narrow distribution in diameter of 2 to 3 m, as shown. The rod diameter showed a weak dependence on process conditions. Anisotropic etch behavior between different phases provides the opportunity to fabricate structures with high aspect ratios. The photomicrographs show the resulting microstructure after a wet chemical etch and a dry plasma etch. The wet chemical etches the silicon away, exposing the TiSi2 rods, whereas plasma etching preferentially etches the Si-TiSi2 interface to form a crater. The porous architectures are applicable to fabricating microdevices or creating templates for part fabrication. The porous rod structure can serve as a platform for fabricating microplasma devices for propulsion or microheat exchangers and for fabricating microfilters for miniatured chemical reactors. Although more work is required, self-assembly from DSE can have a role in microdevice fabrication.

Dynys, Frederick W.↗

Atomic-Scale Structure and Non-Stoichiometry of Meteoritic Hibonite: A Transmission Electron Microscope Study

Hibonite (CaAl_12O_19) is a common refractory mineral in Ca-Al-rich inclusions (CAIs) in primitive meteorites. Transmission electron microscope (TEM) studies have identified enigmatic planar defects in different occurrences of hibonite in the Allende meteorite that give rise to strong streaking along c* in electron diffraction patterns. Atomic resolution high-angle annular dark-field (HAADF) imaging and energy-dispersive X-ray (EDX) analyses were used to determine the nature and origin of these planar features. HAADF images of hibonite grains reveal lamellar intergrowths of common 1.6 nm spacing, and less commonly 2.0 and 2.5 nm spacings, interspersed in stoichiometric hibonite showing 1.1 nm (002) spacing. Stoichiometric hibonite consists of alternating Ca-containing (“R”) and spinel-structured (“S”) blocks stacked in a sequence RS. In contrast, the 1.6 nm layers result from a doubled S block such that the stacking sequence is RSS, while in the widest defect observed, the stacking sequence is RSSSS. These intergrowths are epitaxial and have coherent, low-strain boundaries with the host hibonite Meteoritic hibonite shows common Ti and Mg substitution for Al in its structure. Atomic-resolution EDX maps of hibonite grains in the Allende CAI confirm the preferred site occupancy of Mg on tetragonal M3 sites in S blocks and of Ti on trigonal bipyramidal M2 and octahedral M4 sites in R blocks. Mg is highly concentrated, but Ti is absent in the planar defects where wider S blocks show Al-rich compositions compared to stoichiometric MgAl_2O_4 spinel. Therefore, Mg likely played the major role in the formation and metastability of planar defects in hibonite. Electron energy loss spectroscopy data from the Ti L_2,3 edge show the presence of mixed Ti oxidation states with ~ 15–20% of Ti as Ti^(3+) in hibonite, suggesting a direct substitution of Ti^(3+) ↔ Al^(3+) in hibonite. The remaining ~ 80–85% of Ti is present as Ti^(4+) and corresponding EDX analyses are consistent with the well-known coupled substitution 2Al^(3+) ↔ Ti^(4+) + Mg^(2+) being the major mechanism for Ti and Mg substitution in hibonite. The formation of planar defects in hibonite occurred during high-temperature nebular condensation or melting/crystallization processes. The occurrence of non-stoichiometric hibonite in the Allende CAI deviates from the mineral formation sequence predicted from equilibrium condensation models. Overall, our atomic resolution TEM observations signify non-equilibrium, kinetic-controlled crystal growth during the high-temperature formation of refractory solids in the early solar nebula.

hibonite↗

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fluorine negative ion detachment kinetics

A study of the rate of F(-) detachment by O and H atoms via the reactions F(-) + O yields FO + e and F(-) + H yields FH+ e was undertaken using a drift tube to produce F(-) ions at various drift velocities and therefore different ion temperatures. Preliminary mobility measurements of F(-) ions in Ar were made, indicating that ion temperatures in the 300 K to 5000 K range could be achieved; however due to numerous difficulties experienced in obtaining a reliable F(-) ion source, the study could not be completed.

Burke, R. R.↗

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstractions by $CH_3\dot{O}$ radical

H-atom abstraction by methoxy radical ($CH_3\dot{O}$) plays an important role in capturing the kinetics of reactions between gasoline components and alcohols. This study focuses on determining the reaction rates and thermodynamic properties of methoxy radical reactions with five gasoline fuel components: n-heptane, iso-octane, 1-hexene, cyclopentane and toluene. Electronic structure calculations were per-formed for all the stationary points with M06-2X/6 -311 ++ g(d,p) method. G3 composite method with atomization method is used for determining Δ f H 0 of all the closed shell and radical species, using which the necessary thermodynamic data of all the species was determined. Coupled cluster theory QCISD(T)/cc-pVXZ (where X = D and T) and Møller-Plesset perturbation theory MP2/cc-pVXZ (where X = D, T and Q) were used to calculate single point energies. Subsequently, rate constants for all hydrogen atom ab-straction channels have been performed using conventional transition state theory with unsymmetric tunneling corrections. A systematic comparison of rates for abstraction from different sites within the same species and same site from different species is done in order to get insights into this reaction class. Here, the computed thermodynamic properties and rate constants were incorporated into a recent gasoline mechanism to investigate the impact of the calculations performed in this work. A shift in predicted NTC (negative temperature coefficient) behavior and significant reduction in model reactivity is observed upon incorporating the rates calculated herein.

33 ADVANCED PROPULSION SYSTEMS↗

A quantitative phase-field model for gas bubble evolution in UO2

Due to the large formation energy of vacancies and noble gas atoms in the form of interstitials or substitutional atoms in nuclear fuel (UO2), the thermodynamic equilibrium concentrations of these species are very low in the nuclear fuel matrix even at very high temperatures, which imposes difficulties upon the quantitative study of bubble evolution via the phase-field method. In this study, a quantitative phase-field model is proposed to deal with this problem. The system’s free energy density is derived according to the principles of thermodynamics, with consideration of the elastic effect and with the use of real material parameters from experiments. The model is useful for the study of the kinetics of gas bubble growth with very dilute concentrations of vacancy and gas atoms in the matrix. This model is applied to study single bubble growth and multiple bubble growth under various concentrations of vacancy and gas atoms and at various temperatures. The elastic effect and the effects of the generation rate of vacancies and gas atoms on bubble growth are analyzed.

low defect concentration, nuclear fuel, gas bubble↗