Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fuel lattices”

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 55 records · Page 3

The Impact of Micro Porous Layer on Liquid Water Evolution inside PEMFC using Lattice Boltzmann Method

Proton exchange membrane fuel cells (PEMFCs) require a gas diffusion layer to aid in fuel transport to the catalyst sites. A microporous layer (MPL) is often added to the GDL to improve liquid saturation inside gas diffusion media. In this work, the lattice Boltzmann method was applied to four GDL samples with the addition of an MPL. Three injection orientations were used to study liquid evolution through the samples. Each orientation used four different injection pressures, ranging from 5,000 Pa to 8,000 Pa. Saturation data for GDL samples with and without an MPL were compared. Results showed that when adding an MPL, liquid tends to distribute laterally under the MPL until pressure is reached to allow liquid to travel through the cracks of the MPL surface and into the GDL geometry. A more uniform saturation distribution across the sample is seen when comparing both types of GDLs. The outcomes of this work will help research that requires knowledge of the internal liquid transport through gas diffusion media for PEMFC application.

Sepe, M.↗

Impact of Irradiation-induced Lattice Defects and Dopants on Thermal Transport in Oxide Nuclear Fuels

Advanced nuclear reactor concepts aim to use fuels that must withstand unprecedented temperature and radiation extremes. In these fuels, thermal energy transport under irradiation is directly related to fuel longevity, reactor safety, and is arguably one of the most important performance metrics. Here we provide a comprehensive, first-principles-informed treatment of phonon mediated thermal transport in a defect-bearing actinide oxide with direct comparison to experimental measurements. Pristine and proton irradiated thorium dioxide was chosen as a model system to treat the complexity of thermal transport in the presence of lattice defects. A thermal transport model is implemented using the linearized Boltzmann transport equation (LBTE) with input from first principles calculations and defect evolution models. The output of the LBTE is compared directly to mesoscopic measurements of thermal conductivity on length scales commensurate with defect accumulation. Parametric measurements of conductivity with irradiation dose and temperature reveal a monotonic decrease in conductivity with irradiation dose, in the range of 0.001 dpa to 0.1 dpa. Besides radiation-induced defects, doping the thorium dioxide crystal with small quantities of uranium atoms also results in a drastic reduction in thermal conductivity. Temperature-dependent measurements of thermal conductivity in uranium-doped thorium dioxide single crystals suggests a resonant scattering mechanism is responsible for the observed reduction. This comprehensive, atomistic- to meso-scale treatment provides the necessary basis to investigate thermal transport under irradiation in more complex systems that exhibit strong electron correlation. This is a poster presentation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NCSP Integral Experiments at Sandia in FY21 [Slides]

This presentation provides details regarding integral experiments at Sandia National Laboratory for fiscal year 2021. The experiments discussed are as follows: IER 230: Characterize the Thermal Capabilities of the 7uPCX; IER 304: Temperature Dependent Critical Benchmarks; IER 305: Critical Experiments with UO 2 Rods and Molybdenum Foils; IER 306: Critical Experiments with UO 2 Rods and Rhodium Foils ; IER 441: Epithermal HEX Lattices with SNL 7uPCX Fuel for Testing Nuclear Data; IER 452: Inversion Point of the Isothermal Reactivity Coefficient; and IER 523: Critical Experiments with ACRR UO 2 -BeO Fuel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Origin of High Activity and Durability of Twisty Nanowire Alloy Catalysts under Oxygen Reduction and Fuel Cell Operating Conditions

The ability to control the surface composition and morphology of alloy catalysts is critical for achieving high activity and durability of catalysts for oxygen reduction reaction (ORR) and fuel cells. Here, this report describes an efficient surfactant-free synthesis route for producing a twisty nanowire (TNW) shaped platinum–iron (PtFe) alloy catalyst (denoted as PtFe TNWs) with controllable bimetallic compositions. PtFe TNWs with an optimal initial composition of ~24% Pt are shown to exhibit the highest mass activity (3.4 A/mg Pt , ~20 times higher than that of commercial Pt catalyst) and the highest durability (<2% loss of activity after 40000 cycles and <30% loss after 120000 cycles) among all PtFe-based nanocatalysts under ORR or fuel cell operating conditions reported so far. Using ex situ and in situ synchrotron X-ray diffraction coupled with atomic pair distribution function (PDF) analysis and 3D modeling, the PtFe TNWs are shown to exhibit mixed face-centered cubic (fcc)-body-centered cubic (bcc) alloy structure and a significant lattice strain. A striking finding is that the activity strongly depends on the composition of the as-synthesized catalysts and this dependence remains unchanged despite the evolution of the composition of the different catalysts and their lattice constants under ORR or fuel cell operating conditions. Notably, dealloying under fuel cell operating condition starts at phase-segregated domain sites leading to a final fcc alloy structure with subtle differences in surface morphology. Due to a subsequent realloying and the morphology of TNWs, the surface lattice strain observed with the as-synthesized catalysts is largely preserved. This strain and the particular facets exhibited by the TNWs are believed to be responsible for the observed activity and durability enhancements. These findings provide new insights into the correlation between the structure, activity, and durability of nanoalloy catalysts and are expected to energize the ongoing effort to develop highly active and durable low-Pt-content nanowire catalysts by controlling their alloy structure and morphology.

25 ENERGY STORAGE↗

High-burnup isotopics bias and uncertainty estimations for LANCR02

Calculated nuclide concentrations are compared against destructive radiochemical assays of modern, high-burnup BWR fuel to validate the LANCR02 lattice physics code. Good agreement is observed between different lattice physics methodologies - including the LANCR02 implementation of Method of Characteristics neutron and gamma transport solver and the TGBLA06 methodology - and against measurements. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SCALE Code System

The SCALE Code System is a widely used modeling and simulation suite for nuclear safety analysis and design that is developed, maintained, tested, and managed by the Reactor and Nuclear Systems Division (RNSD) of Oak Ridge National Laboratory (ORNL). SCALE provides a comprehensive, verified and validated, user-friendly tool set for criticality safety, reactor and lattice physics, radiation shielding, spent fuel and radioactive source term characterization, and sensitivity and uncertainty analysis. Since 1980, regulators, licensees, and research institutions around the world have used SCALE for safety analysis and design. SCALE provides an integrated framework with dozens of computational modules, including three deterministic and three Monte Carlo radiation transport solvers that are selected based on the desired solution strategy. SCALE includes current nuclear data libraries and problem-dependent processing tools for continuous-energy (CE) and multigroup (MG) neutronics and coupled neutron-gamma calculations, as well as activation, depletion, and decay calculations. SCALE includes unique capabilities for automated variance reduction for shielding calculations, as well as sensitivity and uncertainty analysis. SCALE’s graphical user interfaces assist with accurate system modeling, visualization of nuclear data, and convenient access to desired results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Linear source method of characteristics in HELIOS2

The Studsvik Scandpower code HELIOS2 is a fast lattice transport code that has gained wide use for reactor fuel cycle analysis in commercial, regulatory, and research environments. It contains two transport solvers - Collision Probabilities (CP) and Method of Characteristics (MoC) - with a general geometry ability and generalized depletion for a variety of fuel types, experimental facilities, and isotope production applications. The HELIOS2 lattice code has been extended with the development of a linear source MoC option. The linear source equations are derived within the scope of the HELIOS2 general geometry and non-cyclic ray tracing. The implementation was tested for fine and coarse spatial mesh models of the 2D C5G7 MOX computational benchmark, which provides a pure and extreme test of a transport solver. The linear source yields at least comparable accuracy with reduced computational resource requirements relative to the flat-source approximation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Reduced-order modeling of neutron transport separated in energy by Minimax Proper Generalized Decomposition

In this article, we demonstrate a Petrov-Galerkin Proper Generalized Decomposition (PGD) known as Minimax PGD for modeling neutron transport separated in energy. To compare the Minimax with the classical Galerkin PGD, we assess both on a model problem of UO{sub 2} or Mixed Oxide (MOX) fuel pins in an infinite lattice with 3 industry-standard energy meshes. We find the Minimax PGD achieves a superior decomposition to Galerkin PGD, both with and without update of the energy modes. This suggests Minimax PGD may be more computationally efficient, provided this reduction in modes (to achieve a given accuracy) outweighs the cost of solving the necessary adjoint problems. In either case, we note that PGD offers an a priori Reduced-Order Model (ROM) which may be dramatically cheaper to solve than the full-order model, especially in problems with fine to ultrafine energy meshes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Lattice Physics Calculations Using the Embedded Self-Shielding Method in Polaris, Part I: Methods and Implementation

Polaris is a 2-dimensional multigroup lattice physics capability in the SCALE code system for the analysis of light water reactor fuel designs. The goal of light water reactor lattice physics codes is to generate few-group homogenized cross sections for downstream full-core nodal diffusion calculations. Additionally, lattice physics calculations contain three primary components: the cross section processing calculation, the 2D transport calculation, and the depletion calculation. This paper summarizes the calculational methods and their implementation into Polaris, with an emphasis on implementation of the embedded self-shielding method. The accuracy of the embedded self-shielding method depends on the procedure used to generate self-shielding factors on the multigroup library. Numerical benchmarks calculations reveal that the accuracy of Polaris eigenvalue predictions is enhanced by (1) using heterogeneous unit cell models to generate the self-shielding factors on the library and (2) using self-shielding factors for within-group scattering cross section.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Ceramic Composite Inert Matrix Fuel Forms in High-Temperature Gas-Cooled Microreactors

Here, this work optimizes micro-prismatic high-temperature gas reactor (HTGR) designs to reduce the energy-normalized mass of spent nuclear fuel (SNF) and high-level waste (HLW) produced. The optimization was performed for the current graphite moderator and an inert matrix fuel (IMF) concept employing different composite moderators in a prismatic design architecture. The fuel matrix is magnesium oxide (MgO) with entrained tristructural-isotropic (TRISO) fuel. The moderator materials, including beryllium oxide (MgO-BeO) and beryllium (MgO-Be) at 40 vol % loading and yttrium hydride (MgO-YH x=1.9 ) and zirconium hydride (MgO-ZrH x=1.9 ) at 15 vol % loading, were entrained within the MgO host matrix. A generic graphite micro-prismatic HTGR is used as the baseline point design where the external dimensions are held constant. The composite moderator designs use 19.9% enriched uranium nitride TRISO fuel and hexagonal assemblies. For each IMF concept, an optimization study was performed to maximize the discharge burnup of the fuel by varying the TRISO packing fraction and the lattice pitch of the assemblies. The mass of SNF and HLW, other waste metrics, fuel cost, environmental impact metrics, and the activity of the SNF and HLW at 100 years and 100 000 years were calculated for the optimized IMF and graphite reference designs. The IMF results were subsequently compared to those of the graphite reference and the values for a light water reactor (LWR) and a small modular LWR. For the SNF and HLW, all the IMF concepts and the graphite reference produced less waste compared to the traditional LWR designs. However, the IMF concepts outperformed the graphite reference regarding the mass of SNF and HLW. For the other waste metrics, the IMF concepts showed reductions in fuel cost with improved environmental metrics relative to the graphite reference. Overall, the IMF concepts significantly reduced the SNF and HLW produced per unit of energy generated compared to traditional LWR designs.

TRISO↗

Microstructural and crystallographic effects of sol-gel synthesized Ti-doped UO 2 sintered under reducing conditions

Titanium (Ti)-doped UO 2 microspheres of three different Ti concentrations (1000, 2000, and 4000 wppm) were synthesized using an internal gelation process. The microspheres were pressed into pellets, and a two-step heat treatment was applied to form monolithic cylindrical pellets with high densities (≥95%TD). Microstructure of these samples consisted of equiaxed grains with >300% increase in average grain size compared to the undoped UO 2 pellets. Secondary Ti-rich chemical phases corresponding to a liquid eutectic formed during sintering were observed at grain boundaries of UO 2 for samples doped with 4000 wppm Ti. Furthermore, these Ti-rich chemical phases were not observed in 1000 or 2000 wppm Ti samples at microscale using electron microscopy investigations. The 0.02–0.04% lower lattice parameter values for the Ti-doped UO 2 samples compared to the undoped UO 2 confirms the incorporation of Ti into the UO 2 lattice.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Cell Dancoff-based embedded self-shielding capability for doubly heterogeneous particulate fuels in SCALE/Polaris

Polaris is a two-dimensional (2D) transport lattice physics capability in the SCALE code system for light-water reactor (LWR) analysis. SCALE/Polaris is being extended to treat doubly heterogeneous particulate fuels in LWRs with accident-tolerant fuels and prismatic high-temperature gas-cooled reactor (HTGR) fuel. Recently, a Dancoff-based Wigner-Seitz approximation capability based on the embedded self-shielding method (DWA-ESSM) was implemented into Polaris for efficient computation. A new double heterogeneity treatment capability based on DWA-ESSM in conjunction with the Hebert's collision probability method for double heterogeneity was developed and implemented for Polaris. The new capability was verified by performing benchmark calculations for the prismatic HTGR fuel compact problems with various design parameters through a code-to-code comparison between the Polaris and continuous-energy Monte Carlo results. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

SCALE 6.2 Lattice Physics Performance Assessment

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Computational Investigation of Fuel Dispersal Phenomena during Large-Break Loss of Coolant Accident in Light-Water Reactors

In the event of cladding rupture, which could occur in light water reactor fuel assemblies during a loss-of-coolant accident (LOCA), fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod. This expulsion of fragmented fuel particles, referred as fuel dispersal, is the subject of investigation to evaluate the safety implications of increasing fuel burnup in light water reactors, with a specific focus on fuel fragmentation, relocation, and dispersal. Particle trajectories and the resulting mass distribution of the settled particles within the reactor pressure vessel can pose a long-term cooling challenge for the reactor core. The fuel dispersal phenomenon is significantly influenced by the ejection characteristics of the fuel fragments, as well as the size and shape of the cladding rupture and fuel rod depressurization history during LOCA transients. In this study, the transport of fuel particles within a scaled 5 × 5 lattice of a pressurized water reactor rod bundle geometry is modeled through a two-fluid Eulerian framework that treats the gas and solid phases as interpenetrating continua. The required boundary conditions are evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considers solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum and energy exchange between the gas and solid phases. The simulation results provide the volume fraction of the solids settled on the bottom surface of the fuel bundle, quantifying the deposition within the bundle geometry.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Implications of point defect accumulation on UO 2 thermal conductivity and fission gas release under accelerated fuel irradiation

Evaluation of thermal properties is a crucial factor for nuclear fuel performance. During reactor operation, the accumulation of fission products and irradiation-induced lattice defects are responsible for degradation in thermal conductivity. Consequently, it affects fuel temperature and fission gas release (FGR) among other Multiphysics processes important for economics and safety analysis. We analyze the implications of point defects (PD) accumulation described using a rate theory (RT) Model on lattice thermal conductivity of UO 2 . Here, we demonstrate that fission rate-dependent point defect concentrations have the largest impact on in-pile thermal conductivity in the periphery of light water reactor fuels below a temperature threshold governed by the migration barrier of defects. Our analysis provides a mechanistic description of this phenomena which current fuel performance codes treat empirically. The reduction of thermal conductivity in the low -temperature rim region acts as additional thermal resistance and leads to a temperature notably larger than suggested by Lucuta thermal conductivity correlation. These effects are anticipated to have notable impacts when fuels are exposed to accelerated radiation. The impact of such point defect-informed treatment of thermal conductivity on fuel performance is evaluated by a detailed analysis of fission gas behavior and its release. We consider several models capturing different stages of fission gas bubble evolution and fission gas release (FGR). Finally, a new fission rate-dependent correction to the Lucuta correlation is proposed. The results show a significant reduction in thermal conductivity at the fuels’ periphery and an increase in fuel centerline temperature specifically at low burnups. Ultimately a modified LC shows a higher FGR compared to the original LC, while the acceleration process results in a reduction in overall FGR.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Engineering of a Coupled Nanocomposite as a High-Performance Protonic Ceramic Fuel Cell Cathode

The lack of high-performance cathode catalysts is a salient issue that bedeviled the commercialization of protonic ceramic fuel cells (PCFCs). Here, in this work, we report a remarkable electrocatalytic activity and stability enhancement of cathode electrodes by engineering a coupled nanocomposite. The as-prepared Pr 0.3 (Ba 0.5 Sr 0.5 ) 0.7 Co 0.8 Fe 0.2 O 3−δ nanocomposite possesses a bulk cubic phase on which homogeneous and intimate orthorhombic PrCo 0.5 Fe 0.5 O 3−δ nanoparticles are uniformly decorated. X-ray diffraction and Raman spectroscopy reveal the excellent thermal stability of the nanocomposite. It achieves a high peak power density of 1.02 W cm –2 based on protonic electrolytes at 600 °C. No noticeable structural degradation is observed over ∼210 h at 550 °C according to scanning electron microscopy analysis. This work demonstrates an effective strategy to boost the performance of perovskite oxides for PCFCs via nanocomposite engineering. It may apply to other catalyst designs and discoveries, such as for batteries, electrolyzers, and membrane reactors.

08 - HYDROGEN↗

Electrochemical Degradation of Pt 3 Co Nanoparticles Investigated by Off-Lattice Kinetic Monte Carlo Simulations with Machine-Learned Potentials

In fuel cell applications, the durability of catalysts is critical for large-scale industrial implementation. However, limited synthesis controllability and spectroscopic resolution impede a comprehensive understanding of degradation mechanisms at the atomic level. In this study, we develop a machine-learned potential (MLP) to simulate the degradation processes for Pt3Co nanoparticles. The precision of MLP is determined to be comparable to that of density functional theory calculations. Using off-lattice kinetic Monte Carlo simulations with MLP, we successfully replicate established experimental trends and offer a logical resolution to ongoing debates regarding atomic orderings. Based on the simulation results, we suggest design principles for Pt3Co nanoparticles that combine high activity and durability. Finally, we validate the wide applicability of our method by successfully applying it to Pt3Ni and Pt3Co0.5Ni0.5 nanoparticles. The research serves as a guideline for developing MLPs for alloy electrochemical catalysts and lays the foundation for designing more durable and active fuel-cell catalysts.

36 MATERIALS SCIENCE↗