Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling constant”

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 1,513 records · Page 84

Vivid Motor Imagery as an Adaptation Method for Head Turns on a Short-Arm Centrifuge

Artificial gravity (AG) has been proposed as a potential countermeasure to the debilitating physiological effects of long duration space flight. The most economical means of implementing AG may be through the use of a short-radius (2m or less) centrifuge. For such a device to produce gravitational forces comparable to those on earth requires rotation rates in excess of 20 revolutions per minute (rpm). Head turns made out of the plane of rotation at these rates, as may be necessary if exercise is combined with AG, result in cross-coupled stimuli (CCS) that cause adverse side effects including motion sickness, illusory sensations of motion, and inappropriate eye movements. Recent studies indicate that people can adapt to CCS and reduce these side effects by making multiple head turns during centrifuge sessions conducted over consecutive days. However, about 25% of the volunteers for these studies have difficulty tolerating the CCS adaptation paradigm and often drop out due to motion sickness symptoms. The goal of this investigation was to determine whether vivid motor imagery could be used as a pseudostimulus for adapting subjects to this unique environment. Twenty four healthy human subjects (14 males, 10 females), ranging in age from 21 to 48 years (mean 33, sd 7 years) took part in this study. The experimental stimuli were produced using the NASA JSC short-arm centrifuge (SAC). Subjects were oriented supinely on this device with the nose pointed toward the ceiling and head centered on the axis of rotation. Thus, centrifuge rotation was in the body roll plane. After ramp-up the SAC rotated clockwise at a constant rate of 23 rpm, producing a centrifugal force of approximately 1 g at the feet. Semicircular canal CCS were produced by having subjects make yaw head turns from the nose up (NU) position to the right ear down (RED) position and from RED to NU. Each head turn was completed in about one second, and a 30 second recovery period separated consecutive head movements. Participants were randomly assigned to one of three groups (n=8 per group): physical adapters (PA), mental adapters (MA), or a control group (CG). Each subject participated in a one hour test session on each of three consecutive days. Each test session consisted of an initial (preadaptation) period during which the subject performed six CCS maneuvers in the dark, followed by an adaptation period with internal lighting on the centrifuge, and a final (postadaptation) period during which six more CCS maneuvers were performed in the dark. For the PA group, the adaptation period consisted of performing 30 additional CCS maneuvers in the light. For the MA and CG group the centrifuge was ramped down to 0 rpm after the pre-adaptation period and ramped back up to 23 rpm before the post-adaptation period. For the both of these groups, the adaptation period consisted of making 30 CCS maneuvers in the light with the centrifuge stationary (so no cross-coupling occurred). MA group subjects were instructed to vividly imagine the provocative sensations produced by the preadaptation CCS maneuvers in terms of magnitude, duration, and direction of illusory body tilt, as well as any accompanying levels of motion sickness. CG group subjects were asked to answer low imagery content questions (trivial pursuit) during each adaptation period head turn. During the 30 second recovery following each head turn, psychophysical data were collected including self reports of motion sickness, magnitude and direction estimates of illusory body tilt, and the overall duration of these sensations. A multilevel mixed effects linear regression analysis performed on all response variables indicated that all three groups experienced some psychophysical adaptation across the three test sessions. For illusory tilt magnitude, the PA group exhibited the most overall adaptation, followed by the MA group, and the CG group. The slopes of these adaptation trajectories by group over day were significantly diffent from one another. For the perceived duration of sensations, the CG group again exhibited the least amount of adaptation. However, the rates of adaptation of the PA and the MA groups were indistinguishable, suggesting that the imagined pseudostimulus appeared to be just as effective a means of adaptation as the actual stimulus. The MA group's rate of adaptation to motion sickness symptoms was also comparable to the PA group. The use of vivid motor imagery may be an effective method for adapting to the illusory sensations and motion sickness symptoms produced by cross-coupled stimuli. For space-based AG applications, this technique may prove quite useful in retaining astronauts considered highly susceptible to motion sickness as it reduces the number of actual CCS required to attain adaptation.

Newby, N. J.↗

Axion dark matter from frictional misalignment

We study the impact of sphaleron-induced thermal friction on the axion dark-matter abundance due to the interaction of an axion-like particle (ALP) with a dark non-abelian gauge sector in a secluded thermal bath. Thermal friction can either enhance the axion relic density by delaying the onset of oscillations or suppress it by damping them. We derive an analytical formula for the frictional adiabatic invariant, which remains constant along the axion evolution and which allows us to compute the axion relic density in a general set-up. Even in the most minimal scenario, in which a single gauge group is responsible for both the generation of the ALP mass and the friction force, we find that the resulting dark-matter abundance from the misalignment mechanism deviates from the standard scenario for axion masses m a ≳ 100 eV. We also generalize our analysis to the case where the gauge field that induces friction and the gauge sector responsible for the ALP mass are distinct and their couplings to the axion have a large hierarchy as can be justified by means of alignment or clockwork scenarios. We find that it is easy to open up the ALP parameter space where the resulting axion abundance matches the observed dark-matter relic density both in the traditionally over- and underabundant regimes. This conclusion also holds for the QCD axion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Reduced global plant respiration due to the acclimation of leaf dark respiration coupled with photosynthesis

Summary Leaf dark respiration ( R d ) acclimates to environmental changes. However, the magnitude, controls and time scales of acclimation remain unclear and are inconsistently treated in ecosystem models. We hypothesized that R d and Rubisco carboxylation capacity ( V cmax ) at 25°C ( R d,25 , V cmax,25 ) are coordinated so that R d,25 variations support V cmax,25 at a level allowing full light use, with V cmax,25 reflecting daytime conditions (for photosynthesis), and R d,25 / V cmax,25 reflecting night‐time conditions (for starch degradation and sucrose export). We tested this hypothesis temporally using a 5‐yr warming experiment, and spatially using an extensive field‐measurement data set. We compared the results to three published alternatives: R d,25 declines linearly with daily average prior temperature; R d at average prior night temperatures tends towards a constant value; and R d,25 / V cmax,25 is constant. Our hypothesis accounted for more variation in observed R d,25 over time ( R 2 = 0.74) and space ( R 2 = 0.68) than the alternatives. Night‐time temperature dominated the seasonal time‐course of R d , with an apparent response time scale of c. 2 wk. V cmax dominated the spatial patterns. Our acclimation hypothesis results in a smaller increase in global R d in response to rising CO 2 and warming than is projected by the two of three alternative hypotheses, and by current models.

Plant Sciences↗

Exact analytical solution to a transient conjugate heat-transfer problem

An exact analytical solution is found for laminar, constant-property, slug flow over a thin plate which is also convectively cooled from below. The solution is found by means of two successive Laplace transformations when a transient in the plate and the fluid is initiated by a step change in the fluid inlet temperature. The exact solution yields the transient fluid temperature, surface heat flux, and surface temperature distributions. The results of the exact transient solution for the surface heat flux are compared to the quasi-steady values, and a criterion for the validity of the quasi-steady results is found. Also the effect of the plate coupling parameter on the surface heat flux are investigated.

Sucec, J.↗

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗

Mesoscopic analyses of the impact of morphology and operating conditions on the transport resistances in a proton-exchange-membrane fuel-cell catalyst layer

Exploring the origins of local transport resistance and characterizing the oxygen transport resistances in the catalyst layer (RCL) are critical for cost reduction. In this paper, a comprehensive mesoscopic model for simulating coupled transport processes of oxygen and water vapor for different structural parameters under different operating conditions in reconstructed microstructures is proposed. The local transport resistance is calculated after achieving the limiting current density and the effective diffusivity of oxygen. The results demonstrate that RCL increases greatly with decreasing platinum loading (LPt) and the transport resistances in other components of the cell dominate for high-loadings. Both the reduced oxygen permeation coefficient in the ionomer thin-film and the adsorption resistance account for the origins of local transport resistance. The local transport resistance increases with the bare carbon ratio for a constant LPt and Pt/C ratio due to the decreased effective ionomer surface, and increases with the I/C ratio due to the increased ionomer thickness and decreased Knudsen diffusivity. Due to the presence of liquid water, a slight decrease followed by an increase of the local transport resistance versus relative humidity is obtained. The contribution of ionomer thin-films to RCL is more sensitive to liquid saturation compared with that of pores. This journal is

Mu, YT↗

Nonadiabatic transition paths from quantum jump trajectories

We present a means of studying rare reactive pathways in open quantum systems using transition path theory and ensembles of quantum jump trajectories. This approach allows for the elucidation of reactive paths for dissipative, nonadiabatic dynamics when the system is embedded in a Markovian environment. Here, we detail the dominant pathways and rates of thermally activated processes and the relaxation pathways and photoyields following vertical excitation in a minimal model of a conical intersection. We find that the geometry of the conical intersection affects the electronic character of the transition state as defined through a generalization of a committor function for a thermal barrier crossing event. Similarly, the geometry changes the mechanism of relaxation following a vertical excitation. Relaxation in models resulting from small diabatic coupling proceeds through pathways dominated by pure dephasing, while those with large diabatic coupling proceed through pathways limited by dissipation. The perspective introduced here for the nonadiabatic dynamics of open quantum systems generalizes classical notions of reactive paths to fundamentally quantum mechanical processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cr silicate as a prototype for engineering magnetic phases in air-stable two-dimensional transition-metal silicates

Abstract Identifying environmentally inert, ferromagnetic two-dimensional (2D) materials with high Curie temperatures ( T c ) down to the single layer limit has been an obstacle to fundamental studies of 2D magnetism and application of 2D heterostructures to spin-polarized devices. To address this challenge, the growth, structure and magnetic properties of a 2D Cr-silicate single layer on Pt(111) was investigated experimentally and theoretically. The layer was grown by sequentially depositing SiO and Cr followed by annealing in O 2 . Scanning tunneling microscopy (STM), low-energy electron diffraction (LEED), and low energy electron microscopy all indicated a well-ordered layer that uniformly covered the surface, with STM and LEED indicating that the silicate relaxed to its favored lattice constant. Further experimental characterizations demonstrated that the Cr was nominally 3+ but with a lower electron density than typical trivalent Cr compounds. Comparison with theory identified a Cr 2 Si 2 O 9 structure that resembles a single layer of a dehydrogenated dioctahedral silicate. Magnetic circular dichroism in x-ray absorption spectroscopy revealed a ferromagnetically ordered state up to at least 80 K. Theoretical analysis revealed that the Cr in a dehydrogenated Cr-silicate/Pt(111) is more oxidized than Cr in freestanding Cr 2 Si 2 O 9 H 4 layers. This greater oxidation was found to enhance ferromagnetic coupling and suggests that the magnetism may be tuned by doping. The 2D Cr-silicate is the first member of a broad series of possible layered first-row transition metal silicates with magnetic order; thus, this paper introduces a new platform for investigating 2D ferromagnetism and the development of magnetoelectronic and spintronic devices by stacking 2D atomic layers.

2D materials↗

A Three-Wave Model of the Stratosphere with Coupled Dynamics, Radiation and Photochemistry

A zonal mean chemistry transport model (2-D CTM) coupled with a semi-spectral dynamical model is used to simulate the distributions of trace gases in the present day atmosphere. The zonal-mean and eddy equations for the velocity and the geopotential height are solved in the semi-spectral dynamical model. The residual mean circulation is derived from these dynamical variables and used to advect the chemical species in the 2- D CTM. Based on a linearized wave transport equation, the eddy diffusion coefficients for chemical tracers are expressed in terms of the amplitude, frequency and growth rate of dynamical waves; local chemical loss rates; and a time constant parameterizing small scale mixing. The contributions to eddy flux are from the time varying wave amplitude (transient eddy), chemical reactions (chemical eddy) and small scale mixing. In spite of the high truncation in the dynamical module (only three longest waves are resolved), the model has simulated many observed characteristics of stratospheric dynamics and distribution of chemical species including ozone. Compared with the values commonly used in 2-D CTMs, the eddy diffusion coefficients for chemical species calculated in this model are smaller, especially in the subtropics. It is also found that the chemical eddy diffusion has only a small effects in determining the distribution of most slow species, including ozone in the stratosphere.

Shia, Run-Lie↗

Cloud Condensation Nuclei in the Early Martian Atmosphere

Explanations for how early Mars supported warmer and wetter conditions 3.5 to 4 billion years ago have long been debated. A water cloud greenhouse has been proposed as a mechanism of warming early Mars, but modeling studies of these environments produce mixed results. Studies show that both H 2 O and CO 2 clouds are capable of influencing atmospheric and surface temperatures, but the amount of warming induced is very sensitive to the cloud physics assumptions made. Producing a strong water cloud greenhouse effect is dependent on forming clouds as high in the atmosphere as possible to maximize warming. However, the majority of studies of H 2 O and CO 2 clouds on early Mars assume that a constant number of seed nuclei are available at all altitudes for clouds to condense onto. Here we test the validity of this assumption in massive CO 2 atmospheres. Explaining how dust is lofted into the middle atmosphere to serve as seed nuclei for cloud formation is a challenge even for present day Mars. This problem may change in an early Mars environment where having a more massive atmosphere could result in differences in the general circulation, vertical transport, and settling time scales of dust. Here, we begin a first-ever study of the early Mars dust cycle and how it might couple to the water cycle. We present early Mars simulations in which dust is injected at the surface and can be transported throughout the atmosphere. We explore the resulting vertical and horizontal distributions of dust particles and their sizes with an emphasis on the maximum height of the dust. We aim to investigate what the distribution of cloud condensation nuclei may have been on early Mars and how this might impact cloud formation.

Kathryn Steakley↗

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗

Probing Electrolyte Influence on CO 2 Reduction in Aprotic Solvents

Selective CO 2 capture and electrochemical conversion are important tools in the fight against climate change. Industrially, CO 2 is captured using a variety of aprotic solvents due to their high CO 2 solubility. However, most research efforts on electrochemical CO 2 conversion use aqueous media and are plagued by competing hydrogen evolution reaction (HER) from water breakdown. Fortunately, aprotic solvents can circumvent HER, making it important to develop strategies that enable integrated CO 2 capture and conversion. However, the influence of ion solvation and solvent selection within nonaqueous electrolytes for efficient and selective CO 2 reduction is unclear. In this work, we show that the bulk solvation behavior within the nonaqueous electrolyte can control the CO2 reduction reaction and product distribution occurring at the catalyst-electrolyte interface. We study different tetrabutylammonium (TBA) salts in two electrolyte systems with glyme ethers (e.g., 1,2 dimethoxyethane or DME) and dimethyl sulfoxide (DMSO) as a low and high dielectric constant medium, respectively. Using spectroscopic tools, we quantify the fraction of ion pairs that forms within the electrolyte. Also, we show how ion pair formation is prevalent in DME and is dependent on the anion type. More importantly, we show that as ion pair formation decreases within the electrolyte, CO 2 current densities increase, and a higher CO Faradaic efficiency is observed at low overpotentials. Meanwhile, in an electrolyte medium where the ion pair fraction does not change with the anion type (such as in DMSO), a smaller influence of solvation is observed on CO 2 current densities and product distribution. By directly coupling bulk solvation to interfacial reactions and product distribution, we showcase the importance and utility of controlling the reaction microenvironment in tuning the electrocatalytic reaction pathways. Insights gained from this work will enable novel electrolyte designs for efficient and selective CO 2 conversion to desired fuels and chemicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transonic aeroelastic analysis of launch vehicle configurations

A numerical study of the aeroelastic stability of typical launch vehicle configurations in transonic flight is performed. Recent computational fluid dynamics techniques are used to simulate the transonic aerodynamic flow fields, as opposed to relying on experimental data for the unsteady aerodynamic pressures. The flow solver is coupled to an appropriate structural representation of the vehicle. The aerodynamic formulation is based on the thin layer approximation to the Reynolds-Averaged Navier-Stokes equations, where the account for turbulent mixing is done by the two-layer Baldwin and Lomax algebraic eddy viscosity model. The structural-dynamic equations are developed considering free-free flexural vibration of an elongated beam with variable properties and are cast in modal form. Aeroelastic analyses are performed by integrating simultaneously in the two sets of equations. By tracing the growth or decay of a perturbed oscillation, the aeroelastic stability of a given constant configuration can be ascertained. The method is described in detail, and results that indicate its application are presented. Applications include some validation cases for the algorithm developed, as well as the study of configurations known to have presented flutter programs in the past.

Filgueirasdeazevedo, Joao Luiz↗

An accurate ab initio quartic force field for ammonia

The quartic force field of ammonia is computed using basis sets of spdf/spd and spdfg/spdf quality and an augmented coupled cluster method. After correcting for Fermi resonance, the computed fundamentals and nu 4 overtones agree on average to better than 3/cm with the experimental ones except for nu 2. The discrepancy for nu 2 is principally due to higher-order anharmonicity effects. The computed omega 1, omega 3, and omega 4 confirm the recent experimental determination by Lehmann and Coy (1988) but are associated with smaller error bars. The discrepancy between the computed and experimental omega 2 is far outside the expected error range, which is also attributed to higher-order anharmonicity effects not accounted for in the experimental determination. Spectroscopic constants are predicted for a number of symmetric and asymmetric top isotopomers of NH3.

Martin, J. M. L.↗

Erosion of focus rings in capacitively coupled plasma etching reactors

In plasma etching reactors, the structure surrounding the wafer, often called a focus ring (FR), plays an important role in maintaining uniform fluxes of reactants across the wafer. The FR is typically made of dielectric materials. During the plasma etching process, the sheath that forms over the wafer to accelerate ions anisotropically into the surface extends over the FR. Electrical charging of the FR modifies the sheath relative to that over the wafer. On the one hand, one wants the sheath to be uniform across the wafer-FR boundary to enable uniform fluxes to the edge of the wafer. On the other hand, maintaining a high voltage sheath over the FR will erode the FR, which is undesirable as the FR is a consumable component that must be periodically replaced in high volume manufacturing. In this work, we computationally investigated the consequences of dielectric constant ε of the FR materials on erosion of the FR. The series capacitance of the FR and its underlying structure is typically smaller than that of the wafer and its underlying structure. As a result, the FR charges quickly relative to the wafer, which then reduces the voltage across the sheath on top of the FR. The ion energy and angular distributions (IEADs) striking the FR are, therefore, generally lower in energy with a broader angular distribution. With ε = 2, the ion energies striking the middle to the outer edge of the FR are 30–180 eV, whereas for ε = 100, the ion energies are 120–380 eV. At the transition between the wafer and the FR, there is a skew in the IEAD as large as 15° that results from the difference in sheath thickness above the wafer and the FR. This skew and the erosion rate across the FR are functions of the dielectric constant of the FR material. With low ε, the FR charges quickly, less plasma is produced above the FR, and there is less voltage across the sheath that results in less FR erosion. Increasing ε of the FR produces a higher sheath voltage as well as higher ion fluxes over the FR, which increases erosion, while the skew at the edge of the wafer is less severe. Further, the material of the subsurface portion of the FR, which dominates its capacitance, is an important consideration in the design of the substrate assembly.

36 MATERIALS SCIENCE↗

Transition-Edge Hot-Electron Microbolometers for Millimeter and Submillimeter Astrophysics

The millimeter and the submillimeter wavelengths of the electromagnetic spectrum hold a wealth of information about the evolution of the universe. In particular, cosmic microwave background (CMB) radiation and its polarization carry the oldest information in the universe, and provide the best test of the inflationary paradigm available to astronomy today. Detecting gravity waves through their imprint on the CMB polarization would have extraordinary repercussions for cosmology and physics. A transition-edge hot-electron micro - bolometer (THM) consists of a superconducting bilayer transition-edge sensor (TES) with a thin-film absorber. Unlike traditional monolithic bolometers that make use of micromachined structures, the THM em ploys the decoupling between electrons and phonons at millikelvin temperatures to provide thermal isolation. The devices are fabricated photolithographically and are easily integrated with antennas via microstrip transmission lines, and with SQUID (superconducting quantum interference device) readouts. The small volume of the absorber and TES produces a short thermal time constant that facilitates rapid sky scanning. The THM consists of a thin-film metal absorber overlapping a superconducting TES. The absorber forms the termination of a superconducting microstripline that carries RF power from an antenna. The purpose of forming a separate absorber and TES is to allow flexibility in the optimization of the two components. In particular, the absorbing film's impedance can be chosen to match the antenna, while the TES impedance can be chosen to match to the readout SQUID amplifier. This scheme combines the advantages of the TES with the advantages of planar millimeter-wave transmission line circuits. Antenna-coupling to the detectors via planar transmission lines allows the detector dimensions to be much smaller than a wavelength, so the technique can be extended across the entire microwave, millimeter, and submillimeter wavelength ranges. The circuits are fabricated using standard microlithographic techniques and are compatible with uniform, large array formats. Unlike traditional monolithic bolometers that make use of micromachined structures, the THM employs the decoupling between electrons and phonons at millikelvin temperatures to provide thermal isolation. There is no fragile membrane in the structure for thermal isolation, which improves the fabrication yield.

Hsieh, Wen-Ting↗

Electrochemical Corrosion of SIMFUEL: Effects of dissolved H2 and noble metal particles

This report describes the results from electrochemical corrosion tests that were conducted using simulated spent fuel (SIMFUEL) materials comprised of UO2 and surrogate fission products. Two SIMFUEL compositions were tested to quantify the effect of noble metal inclusions and dissolved H2 on the UO2 dissolution rate. One material consisted of UO2 with added lanthanide oxides (UO2 N) and the other consisted of UO2 with added lanthanide oxides and noble metals (UO2-H) at concentrations to simulate high burnup fuel. The electrochemical corrosion tests on the SIMFUEL materials were conducted in aqueous electrolyte solutions (pH 10) that were either saturated with air or purged with a H2/Ar gas mixture to maintain a constant dissolved H2 concentration. The results from electrochemical corrosion tests on SIMFUEL can be applied to qualitatively understand and also quantify the separate effects of water chemistry and fuel composition on the degradation behavior of the UO2 matrix. Open circuit potential (OCP) measurements on a SIMFUEL material of known composition (i.e., known fraction of NMPs at the fuel surface) immersed in a known solution chemistry (i.e., pH, Eh, [O2], [H2]) at known temperature provide qualitative insight into the degradation behavior of the UO2 matrix. If the OCP is above the threshold potential at which the oxidative dissolution of U(IV) to U(VI) occurs, the SNF is expected to degrade by oxidative processes under the experimental conditions. The net currents that are measured during potentiostatic tests, during which the surface potential of the SIMFUEL material is fixed by a potentiostat, can be used to quantify the UO2 degradation rate and optimize the rate constant values used in the Fuel Matrix Degradation Model (FMDM) for half reactions that occur on the SNF surface. Specifically, electrochemical measurements enable the estimation of the total anodic current at ECORR–which is the surface potential at which the total anodic and total cathodic currents are equal–so that the rate constant values for key reactions can be calculated. The open circuit potential measurements, potentiostatic tests, surface property measurements (scanning electron microscopy images and electrochemical impedence spectroscopy plots), and solution elemental composition analyses are being performed to update and optimize the FMDM. A case study is presented herein to show how electrochemical corrosion test results and accompanying characterization results for the UO2 N material in air-saturated solution can be used to validate the SNF surface reaction module of the FMDM. Future electrochemical tests will be conducted to provide quantitative information on the effects of NM content (burnup), H2 concentration, water chemistry, temperature, and galvanic couples with cladding and EBS alloys on UO2 degradation that can be used to improve the accuracy and functionality of the FMDM.

Thomas, Sara↗