Search NASA⌕ Search

SEARCH · Search NASA

Results for “Activity coefficient”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Advances in Studies of Electrode Kinetics and Mass Transport in AMTEC Cells (abstract)

Previous work reported from JPL has included characterization of electrode kinetics and alkali atom transport from electrodes including Mo, W, WRh(sub x), WPt(sub x)(Mn), in sodium AMTEC cells and vapor exposure cells, and Mo in potassium vapor exposure cells. These studies were generally performed in cells with small area electrodes (about 1 to 5 cm(sup 2)), and device geometry had little effect on transport. Alkali diffusion coefficients through these electrodes have been characterized, and approximate surface diffusion coefficients derived in cases of activated transport. A basic model of electrode kinetic at the alkali metal vapor/porous metal electrode/alkali beta'-alumina solid electrolyte three phase boundary has been proposed which accounts for electrochemical reaction rates with a collision frequency near the three phase boundary and tunneling from the porous electrode partially covered with adsorbed alkali metal atoms. The small electrode effect in AMTEC cells has been discussed in several papers, but quantitative investigations have described only the overall effect and the important contribution of electrolyte resistance. The quantitative characterization of transport losses in cells with large area electrodes has been limited to simulations of large area electrode effects, or characterization of transport losses from large area electrodes with significant longitudinal temperature gradients. This paper describes new investigations of electrochemical kinetics and transport, particularily with WPt(sub 3.5) electrodes, including the influence of electrode size on the mass transport loss in the AMTEC cell. These electrodes possess excellent sodium transport properties making verification of device limitations on transport much more readily attained.

electrode kinetics alkali atom transport transport↗
DOE OSTI2026

Evaluating the limitations of Bayesian metabolic control analysis

Bayesian Metabolic Control Analysis (BMCA) is a promising framework for inferring metabolic control coefficients in data-limited scenarios, combining Bayesian inference with linear-logarithmic (lin-log) rate laws. These metabolic control coefficients quantify how changes in enzyme activities affect steady-state fluxes and metabolite concentrations across a metabolic network. However, its predictive accuracy and limitations remain underexplored. This study systematically evaluates BMCA’s ability to infer elasticity values, flux control coefficients (FCC), and concentration control coefficients (CCC) under varying data availability conditions using three synthetic metabolic network models. We demonstrate that BMCA predictions are highly dependent on the inclusion of flux and enzyme concentration data, with the omission of these datasets leading to severe inaccuracies. In our synthetic, enzyme-perturbation datasets, external metabolite concentrations had minimal impact and, in some cases, their exclusion improved predictions; when external-nutrient perturbations were introduced and those concentrations were observed, gains were at most modest. Additionally, we find that posterior estimation with both ADVI and HMC can underestimate large-magnitude elasticities in our synthetic settings, with ADVI showing somewhat higher variance under strong up-regulation; thus, recovering |elasticity| ≳ 1.5 remains challenging regardless of the inference engine. ADVI also fails to accurately infer allosteric interactions, even when regulatory effects are strong. While BMCA maintains reasonable accuracy in partially recovering the rankings of the highest FCC values, its estimates of absolute values remain constrained by prior assumptions and data limitations. Our findings reveal the BMCA algorithm’s strengths and weaknesses, providing guidance on its application in metabolic engineering, and highlighting the need for methodological refinements to enhance its predictive capabilities.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1983

X-ray photoelectron spectroscopy and friction studies of nickel-zinc and manganese-zinc ferrites in contact with metals

X-ray photoelectron spectroscopy analysis and sliding friction experiments were conducted with hot-pressed, polycrystalline Ni-Zn and Mn-Zn ferrites in sliding contact with various transition metals at room temperature in a vacuum of 30 nPa. The results indicate that the coefficients of friction for Ni-Zn and Mn-Zn ferrites in contact with metals are related to the relative chemical activity in these metals: the more active the metal, the higher is the coefficient of friction. The coefficients of friction for the ferrites correlate with the free energy of formation of the lowest metal oxide. The interfacial bond can be regarded as a chemical bond between the metal atoms and the oxygen anions in the ferrite surfaces. The adsorption of oxygen on clean metal and ferrite surfaces increases the coefficients of friction for the Ni-Zn and Mn-Zn ferrite-metal interfaces.

Miyoshi, K.↗
DOE OSTI2019

Static Polystyrene Gate Charge Density Modulation of Dinaphthothienothiophene with Tetrafluorotetracyanoquinodimethane Layer Doping: Evidence from Conductivity and Seebeck Coefficient Measurements and Correlations

We relate the static charging of organic field effect transistor (OFET) gate dielectric polymers, the resulting change in threshold voltage shift, and the modulation of thermoelectric properties of the transistor semiconductor. We utilize a bottom gate OFET structure with cross-linkable polystyrene and plain polystyrene (XLPS/PS) bilayer dielectrics and dinaphthothienothiophene (DNTT) organic molecular semiconductor with tetrafluorotetracyanoquinodimethane (F4TCNQ) doping layer. Conductivity and Seebeck coefficient are measured before and after charging. We find that the conductivity increases significantly, while the Seebeck coefficient decreases correspondingly in a trend of S ∝ ln(σ) after the dielectrics were negatively charged, resulting in an average of fivefold increase in power factor at the maximum charging level. Additionally, temperature-dependent conductivity measurements show that the activation energy decreases with the increasing conductivity. We calculate the Seebeck coefficient based on the activation energy and find the result is consistent with the Mott mobility edge model. This work confirms the charging mechanism in our previous published work, including the location of static charges in the bulk of the dielectric, and demonstrates modulation and enhancement of the thermoelectric properties of organic thermoelectric materials by adjacent chargeable dielectrics without changing the molecular microstructure or applying external gate voltage during operation.

static gate, Seebeck coefficient, organic semicond↗
NASA NTRS1981

Global ozone - Solar activity relationship from satellite measurements

Global ozone variations are evaluated by analyzing satellite measurements, using the technique of integrating out dynamical effects causing large variations in total ozone at individual locations. Analysis of Nimbus 4 reprocessed BUV data was extended to cover the period from April 1970 through December 1975. Temporal and latitudinal variation of biases and yearly change in global biases between the satellite and ground based ozone measurements are indicated. The bias-corrected 6-month-running filtered global means are shown, and the correlation coefficient between ozone and the solar activity index remains high at 0.94. It is concluded that the results give a strong indication of an ozone-solar activity relationship that is caused by solar UV variability.

Keating, G. M.↗
NASA NTRS2014

The Role of Temperature and Humidity on Seasonal Influenza in Tropical Areas: Guatemala, El Salvador and Panama, 2008-2013

Background: The role of meteorological factors on influenza transmission in the tropics is less defined than in the temperate regions. We assessed the association between influenza activity and temperature, specific humidity and rainfall in 6 study areas that included 11 departments or provinces within 3 tropical Central American countries: Guatemala, El Salvador and Panama. Method/ Findings: Logistic regression was used to model the weekly proportion of laboratory-confirmed influenza positive samples during 2008 to 2013 (excluding pandemic year 2009). Meteorological data was obtained from the Tropical Rainfall Measuring Mission satellite and the Global Land Data Assimilation System. We found that specific humidity was positively associated with influenza activity in El Salvador (Odds Ratio (OR) and 95% Confidence Interval of 1.18 (1.07-1.31) and 1.32 (1.08-1.63)) and Panama (OR = 1.44 (1.08-1.93) and 1.97 (1.34-2.93)), but negatively associated with influenza activity in Guatemala (OR = 0.72 (0.6-0.86) and 0.79 (0.69-0.91)). Temperature was negatively associated with influenza in El Salvador's west-central departments (OR = 0.80 (0.7-0.91)) whilst rainfall was positively associated with influenza in Guatemala's central departments (OR = 1.05 (1.01-1.09)) and Panama province (OR = 1.10 (1.05-1.14)). In 4 out of the 6 locations, specific humidity had the highest contribution to the model as compared to temperature and rainfall. The model performed best in estimating 2013 influenza activity in Panama and west-central El Salvador departments (correlation coefficients: 0.5-0.9). Conclusions/Significance: The findings highlighted the association between influenza activity and specific humidity in these 3 tropical countries. Positive association with humidity was found in El Salvador and Panama. Negative association was found in the more subtropical Guatemala, similar to temperate regions. Of all the study locations, Guatemala had annual mean temperature and specific humidity that were lower than the others.

Central America↗
DOE OSTI2025

Cations Enhance Hydride Transfer to Noncatalytic Metals in Concentrated Alkaline Electrolytes

Alkali metal cations are known to influence the kinetics of the hydrogen evolution reaction (HER), acting as either promoters or inhibitors to the rate-determining step depending on the metal surface, local pH, and cation concentration. Despite the importance of concentrated electrolytes for commercial electrochemical cells, the impact of cations on the HER in concentrated alkaline environments (>1 M) and in mixed cation systems remains poorly understood. Here, this study quantifies the HER kinetics at polycrystalline metal surfaces (Pt, Au, Cu, and Fe) and the equilibrium solvation environment in pure and mixed alkali metal hydroxide electrolytes at concentrations up to 3.0 M. Kinetic analyses of Au, Cu, and Fe revealed a positive cation-concentration-effect that was primarily driven by changes to the charge transfer coefficient. Multinuclear NMR spectroscopy examined the solvation of H 2 O/OH – species and the alkali cations as a function of (mixed) alkali cation concentration(s), and demonstrated rapid exchange between solvent, hydroxide, and solvated cations. Together, these findings support models where HER kinetics on noncatalytic metal surfaces in strongly alkaline conditions are primarily governed by the average polarization and polarizability of the metal/solution interface and that increasing cation activity continues to increase the transfer coefficient at metal-hydroxide concentrations up to 3.0 M. Electrolytes and additives which can outcompete weakly hydrated cations and disrupt the interfacial water structure are expected to suppress parasitic HER at electrodes for energy storage and electroplating.

Cations↗
DOE OSTI2023

The P( 4 S) + NH( 3 Σ – ) and N( 4 S) + PH( 3 Σ – )reactions as sources of interstellar phosphorus nitride

Phosphorus nitride (PN) is believed to be one of the major reservoirs of phosphorus in the interstellar medium (ISM). For this reason, understanding which reactions produce PN in space and predicting their rate coefficients is important for modelling the relative abundances of P-bearing species and clarifying the role of phosphorus in astrochemistry. Here, in this work, we explore the potential energy surfaces of the $\textrm{P}(^4\textrm{S}) + \textrm{NH}(^3\Sigma^-)$ and $\textrm{N}(^4\textrm{S}) + \textrm{PH}(^3\Sigma^-)$ reactions and the formation of $\textrm{H}(^2\textrm{S}) + \textrm{PN}(^1\Sigma^+)$ through high accuracy ab initio calculations and the variable reaction coordinate transition state theory (VRC-TST). We found that both reactions proceed without an activation barrier and with similar rate coefficients that can be described by a modified Arrhenius equation ( $k(T)=\alpha\!\left( T/300 \right)^{\beta} \exp\!{(\!-\!\gamma/T)})$ with $\alpha=0.93\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=0.24\, \rm K$ for the $\textrm{P} + \textrm{NH} \longrightarrow \textrm{H} + \textrm{PN}$ reaction and $\alpha=0.88\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=1.01\, \rm K$ for the $\textrm{N} + \textrm{PH} \longrightarrow \textrm{H} + \textrm{PN}$ one. Both reactions are expected to be relevant for modelling PN abundances even in the cold environments of the ISM. Given the abundance of hydrogen in space, we have also predicted rate coefficients for the destruction of PN via H + PN collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Comprehensive Dynamics in a Polyelectrolyte Complex Coacervate

The linear viscoelastic response, LVR, of a hydrated polyelectrolyte complex coacervate, PEC, was evaluated over a range of frequencies, temperatures, and salt concentrations. The PEC was a nearly stoichiometric blend of a quaternary ammonium poly([3-(methacrylamido)propyl]trimethylammonium chloride), PMAPTAC, and poly(2-acrylamido-2-methyl-1-propanesulfonic acid sodium salt), PAMPS, an aliphatic sulfonate, selected because they remain fully charged over the conditions of use. Narrow molecular weight distribution polyelectrolytes were prepared using fractionation techniques. A partially deuterated version of PMAPTAC was incorporated to determine its coil radius of gyration, R g , within PECs using small-angle neutron scattering. Chain dimensions were determined to be Gaussian with a Kuhn length of 2.37 nm, which remained constant from 25 to 65 °C. The LVR for a series of matched molecular weight PECs, mostly above the entanglement threshold, exhibited crossovers of modulus versus frequency classically attributed to the reptation time, relaxation between entanglements, and the relaxation of a Kuhn length of units (the “monomer” time). The scaling for zero shear viscosity, η 0 , versus chain length, N, was η 0 ~ N 3.1 , in agreement with “sticky reptation” theory. The lifetime and activation energy, E p , of a pair between polyanion and polycation repeat units, Pol + Pol – , were determined from diffusion coefficients of salt ions within the PEC. The activation energy for LVR of salt-free PECs was 2E p , showing that the key mechanism limiting the dynamics of undoped PECs is pair exchange. An FTIR technique was used to distinguish whether SCN– acts as a counterion or a co-ion within PECs. Doping of PECs with NaSCN breaks Pol + Pol – pairing efficiently, which decreases effective cross-linking and decreases viscosity. Finally, an equation was derived that quantitatively predicts this effect.

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

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Non-Additive Voltametric Currents From a Mixture of Two, Three and Four Redox-Active Compounds and Electroanalytical Implications

We have published recently the effect of dissimilar diffusion coefficients on the size of the voltammetric waves from a mixture of two redox-active compounds. Similarly, at the potential range where three redox-active species, decamethylferrocene (dMeFc), ferrocene (Fc) and N-methylphenothiazine (MePTZ), are oxidized simultaneously with rates controlled by linear diffusion, electrogenerated radicals diffusing outwards from the electrode react with the original species diffusing towards the electrode from the bulk; thus, Fc(+) reacts with dMeFc producing Fc and dMeFc(+), while MePTZ(+) reacts both with dMeFc producing MePTZ and dMeFc(+), and with Fc producing MePTZ and Fc(+). These reactions replace dMeFc with Fc at the second plateau, and both dMeFc and Fc with MePTZ at the third plateau. Since the diffusion coefficients of the three species are not equal, the mass-transfer limited currents of the second and the third oxidation wave plateaus change by approx. 10%. Numerical simulations of the experimental voltamograms support this mechanism. Similar results were also obtained for a mixture of four redoxactive compounds. The implications of this non-additive nature of currents on: (a) the use of internal voltammetric standards for quantitative analysis of a mixture of redox-active compounds; and, (b) the half wave potentials (E1/2) of the 2nd, 3rd and 4th waves for qualitative analysis, will be discussed.

Dass, Amala↗
DOE OSTI2024

Insights into the soft brittle-to-ductile transition from discrete dislocation dynamics

The Brittle-to-ductile transition (BDT) in body centered cubic metals exhibits a soft transition wherein the fracture toughness gradually rises to before the onset of ductility. The resultant brittle-to-ductile transition temperature can be described with an Arrhenius relationship whose activation energy is related to plasticity in the material. To provide further insight into the nature of the BDTT, in this work we utilized a discrete dislo- cation dynamics model with a crack to simulate the BDT and how it depends on the thermally activated nature of plasticity. The interrelationship between the BDT activation energy and the dislocation mobility parameters were determined via the calculation of first order sensitivity coefficients. This analysis allows us to demonstrate that the activation energy for the BDT is directly related to the activation energy for plasticity through an effective stress that defines this relationship. This effective stress physically is the average stress on the dislocations that move out of the crack. Lastly, we are able to show that this effective stress is dictated by the low temperature fracture toughness or cleave energy of the material and the source position, the latter of which can be affected by processing. Collectively, these results provide new insight into what controls the thermal activation of the BDT and what are the important parameters to control it.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Altimeter and gravity data analysis

The studies carried out under this grant fell into two broad areas. The first area was the analysis of surface gravity data with the ultimate goal of providing normal equations that could be used in combination with normal equations from the analysis of satellite orbit perturbations to obtain an optimal estimate of the gravitational potential coefficients of the Earth. The second main research activity was the estimation of gravity anomalies in ocean areas from satellite altimeter data. Such anomalies could enable the improved calibration of potential coefficient models derived solely from the analysis of orbital perturbation information. The studies in these two areas are discussed.

Rapp, Richard H.↗
DOE OSTI2025

Speciation and diffusive dynamics in hydrated grain boundaries of complex oxide Gd2Ti2O7

Abstract Grain boundaries in polycrystalline materials significantly affect their properties, such as ionic transport, corrosion, and chemical durability. The pyrochlore compound (Gd 2 Ti 2 O 7 ) is employed as a model for complex oxides and is known for its diverse applications, including nuclear waste immobilization. Density functional theory-based first-principles molecular dynamics simulations were performed at different temperatures on the hydrated grain boundary system. The results show extensive transformations within the grain boundaries among hydrous water species (OH − , H 2 O, and H 3 O + ). The temperature dependence of self-diffusion coefficients follows Arrhenius behavior, with an activation energy of 35.9 kJ/mol for hydrogen and 46.3 kJ/mol for oxygen. The lifetime of OH − is about three to four times longer than that of H 2 O at temperatures from 800 to 2100 K, suggesting the greater stability of OH − over H 2 O, a unique characteristic of the grain boundaries. The estimated lifetime of the hydrous species decreases as the temperature increases, with an activation energy of 9.9 kJ/mol for OH − and 13.4 kJ/mol for H 2 O. While Gd 3 + is more mobile than Ti 4+ , both the Gd 3 + and Ti 4+ cations are orders of magnitude less mobile than the water species. The results suggest that water species are much more mobile within grain boundaries than in the bulk crystal and have the potential to penetrate deep into polycrystalline materials through grain boundaries, leading to grain boundary degradation and dissolution. The different mobilities of cations in complex oxides can lead to leaching of certain cations and incongruent dissolution during the chemical weathering of Earth and industrial materials.

B. Ghosh, Dipta↗
DOE OSTI2021

Characterization of the indoor near-field aerosol transmission in a model commercial office building

To evaluate the exposure potential of infectious aerosols containing SARS-CoV-2 in an office building setting, synthetic test aerosols were used to experimental study airborne particle transmission in a multizone small office test building at the Oak Ridge National Laboratory. Nine measurement points in a single zone using active aerosol impactors report that the coefficient of variation of the time-averaged concentration is <10% in two campaigns and < 15% in one campaign, so a nearly well-mixed condition was noted. To understand the effect of HVAC system operation on the dynamic concentration of aerosols in office spaces, an aerosol transport model that includes factors such as outside air (OA) ratio, filtration, return air fraction, transport loss in air ducts, and particle deposition has been developed. The results of model fitting demonstrate strong agreement with experimental data. Overall, our investigation finds the return air fraction effects outweigh other mechanisms for the aerosol recirculation in this study, and the impact of air change rate (ACR) is more important than the small particle deposition for aerosol removal. Because ACR dominates the aerosol transport, the full model can be simplified to just one factor, the ACR, while maintaining an acceptable representation of the experimental data.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Diffusion in undoped and Cr-doped amorphous UO 2

UO 2 fuel pellets are often doped with chromium oxide to obtain favourable properties such as higher density, improved thermal stability, large grain sizes, improved pellet-clad interaction margins, and increased fission gas retention during transients. Chromium has a low solubility limit in UO 2 , with past experimental work reporting solubility limits ranging between 0.004 to 0.06 wt.% Cr. Due to its low solubility, segregation of Cr ions to the grain boundary may occur. Further, the complexity of these boundaries may be high as observed in other ceramics resulting in disordered or amorphous regions along the boundary, affecting a range of material and operational properties of the fuel pellet. To assess these disordered regions, in this work we study amorphous undoped and Cr doped UO 2 systems (containing 10–50 at.% Cr 3+ ) that have been modelled using classical molecular dynamics methods incorporating Cr 3+ into the well-used CRG potential library. Diffusion coefficients, pre-exponential factors, and activation energies for diffusion were computed for oxygen ions, assessing the impact of structure and extrinsic species on migration. Oxygen diffusion was observed to be much faster in the undoped amorphous system compared to its crystalline counterpart. Oxygen diffusion in doped systems decreased with increasing Cr concentration, highlighting the importance of additives to retain fission products and other migratory species.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗