Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary effects”

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 469 records · Page 26

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flux-Closure Domain Structures in Ferroelectric K 0.5 Na 0.5 NbO 3 Thin Films

Topological domain structures in ferroelectric materials have garnered increasing attention due to their intriguing physical properties and promising applications. While most existing topological structures in ferroelectric perovskite oxides originate from tetragonal or rhombohedral bulk phases, much less is understood about their counterparts in orthorhombic ferroelectrics. Here, in this work, we employ ferroelectric K 0.5 Na 0.5 NbO 3 (KNN) thin films as a model system and leverage phase-field simulations to theoretically predict the static structures and dynamic behaviors of three types of flux-closure domain configurations: in-plane (Type-I), out-of-plane (Type-II), and superdomain (Type-III) flux-closure structures. We systematically investigate the effects of finite size, misfit strains, and electrical boundary conditions on the formation and switching of these topological structures. For the Type-I structure, size reduction or small misfit strain facilitates a transition of the flux-closure pattern to polar vortices. Type-II structures emerge under open-circuit electrical boundary conditions of the film, forming at the junctions of specific domain walls with the film surface or the film–substrate interface. The formation mechanisms of these two flux-closure structures are rationalized from an energy perspective. We demonstrated switching capabilities of Type-II and Type-III structures by obtaining polarization–electric field hysteresis loops. Our simulations also reveal a reversible electric-field-induced transition between the orthorhombic and rhombohedral ferroelectric phases with a checkerboard domain pattern, during which the integrity of the flux-closure structure is preserved. These findings provide theoretical insights and practical guidance for identifying and manipulating topological structures in low-symmetry ferroelectrics, paving the way for developing energy-efficient microelectronic devices based on topological structures.

P-E loop↗

In situ energetics modulation enables high-efficiency and stable inverted perovskite solar cells

In contrast to conventional (n-i-p) perovskite solar cells (PSCs), inverted (p-i-n) PSCs offer enhanced stability and integrability with tandem solar cell architectures, which have garnered increasing interest. However, p-i-n cells suffer from energy level misalignment with transport layers, imbalanced transport of photo-generated electrons and holes, and significant defects with the perovskite films. Here we introduce tris(2,4,6-trimethyl-3-(pyridin-3-yl)phenyl)borane (3TPYMB), a nonionic n-type molecule that, through hydrogen bonding and Lewis acid-base reactions with perovskite surfaces or grain boundaries, enables in situ modulation of perovskite energetics, effectively mitigating the key challenges of p-i-n PSCs. The p-i-n PSCs incorporating 3TPYMB achieve a certified quasi-steady-state power conversion efficiency of 24.55 +- 0.33%, with a reverse scan efficiency of 25.58%. Finally, they also exhibit exceptional stability, with unencapsulated devices retaining 97.8% of their initial efficiency after 1,800 h of continuous operation at maximum power point under N 2 atmosphere, 1 sun illumination and 60 °C conditions.

14 SOLAR ENERGY↗

Atomic forces from Dirac-Kohn-Sham equations: Implementation in flexible (APW+lo/LAPW)+LO basis set

The atomic forces formulation based on the Dirac-Kohn-Sham equation and the flexible (APW+lo/LAPW)+LO basis set is presented here. The formulation was implemented in the code FlapwMBPT and allows a user to easily switch between different basis functions of the augmentation type (APW or LAPW) and between different kind of local orbitals. Similar to the work (Phys.Rev.B 91 (2015) 035105), the implementation takes into account small discontinuities of the wave functions, density, and potential at the muffin-tin sphere boundaries. Applications to the materials with strong relativistic effects, such as $\alpha$-Uranium, PuCoGa$_{5}$, and FePt, demonstrate robustness of the method. The comparison of the calculated forces with the ones obtained by numerical differentiation of free energy shows close agreement with deviations about 0.1% or less.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Anisotropic weak antilocalization in thin films of the Weyl semimetal TaAs

Device applications of topological semimetals await the development of epitaxial films in the ultrathin limit. Weak antilocalization (WAL) has been extensively utilized in the understanding of surface states in topological insulators and shows promise for use in elucidating the properties of thin film topological semimetals. Here, we report insights from WAL in the surface state and interface transport properties of our recently synthesized single-crystal-like thin films of the Weyl semimetal TaAs(001) grown on GaAs(001). We observe robust, anisotropic WAL in the magnetoconductance from 2 to 20 K in films from 10 to 200 nm thick. We link the anisotropic WAL magnetoconductance to anisotropic mobility stemming from film topography. We conclude that WAL in the films likely originates from the antilocalization of topological surface states. The WAL magnetoconductance is impacted by the film thickness and topography, solidifying the useful role of WAL in the study of topological semimetal/semiconductor heterointerfaces. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Transport, correlations, and chaos in a classical disordered anharmonic chain

Here, we explore transport properties in a disordered nonlinear chain of classical harmonic oscillators, and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size $\textit{N}$, disorder strength Δ, and temperature $\textit{T}$. The conductivity $κ_N$, obtained for finite length ($\textit{N}$), saturates to a value $κ_∞$ > 0 in the large $\textit{N}$ limit, for all values of disorder strength Δ and temperature $\textit{T}$ > 0. We show evidence that for any Δ > 0 the conductivity goes to zero faster than any power of $\textit{T}$ in the ($\textit{T}$/Δ) → 0 limit, and find that the form $κ_∞$ ~ $e^{{–B|\text{ln}(CΔ/T)|}^{3}}$ fits our data. This form has earlier been suggested by a theory based on the dynamics of multioscillator chaotic islands. The finite-size effect can be $κ_N < κ_∞$ due to boundary resistance when the bulk conductivity is high (the weak disorder case), or $κ_N > κ_∞$ due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the out-of-time-ordered commutator.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The Making of a Metric: Co-Producing Decision-Relevant Climate Science

Developing decision-relevant science for adaptation requires the identification of climatic parameters that are both actionable for practitioners as well as tractable for modelers. In many sectors, these decision-relevant climatic metrics and the approaches that enable their identification remain largely unknown. “Co-production” of science with scientists and decision-makers is one potential way to identify these metrics, but there is little research describing specific and successful co-production approaches. This paper examines the negotiations and outcomes from Project Hyperion, wherein scientists and water managers jointly developed decision-relevant climatic metrics for adaptive water management. We identify successful co-production strategies by analyzing the project’s numerous back-and-forth engagements and tracing the evolution of the science during these engagements. We found that effective mediation between scientists and managers needed dedicated “boundary spanners” with significant modeling expertise. Translating practitioners’ information needs into tractable climatic metrics required direct and indirect methods of eliciting knowledge. We identified four indirect methods that were particularly salient for extracting tacitly held knowledge and enabling shared learning: developing a hierarchical framework linking management issues with metrics, starting discussions from the planning challenges, collaboratively exploring the planning relevance of new scientific capabilities, and using analogies of other “good” metrics. The decision-relevant metrics we developed provide insights into advancing adaptation-relevant climate science in the water sector. The co-production strategies we identified can be used to design and implement productive scientist–decision-maker interactions. Overall, the approaches and metrics we developed can help climate science to expand in new and more use-inspired directions.

Planning↗

Free fermion entanglement with a semitransparent interface: the effect of graybody factors on entanglement islands

We study the entanglement entropy of free fermions in 2d in the presence of a partially transmitting interface that splits Minkowski space into two half-spaces. We focus on the case of a single interval that straddles the defect, and compute its entanglement entropy in three limits: Perturbing away from the fully transmitting and fully reflecting cases, and perturbing in the amount of asymmetry of the interval about the defect. Using these results within the setup of the Poincaré patch of AdS _2 2 statically coupled to a zero temperature flat space bath, we calculate the effect of a partially transmitting AdS _2 2 boundary on the location of the entanglement island region. The partially transmitting boundary is a toy model for black hole graybody factors. Our results indicate that the entanglement island region behaves in a monotonic fashion as a function of the transmission/reflection coefficient at the interface.

Kruthoff, Jorrit↗

Uncertainties in Local Intense Precipitation Flood Modeling

In the Task 1 report for the U.S. Nuclear Regulatory Commission (NRC) funded Local Intense Precipitation (LIP) Probabilistic Flood Hazard Assessment (PFHA) Pilot Study, the Pacific Northwest National Laboratory (PNNL) reviewed hydrologic and hydraulic modeling approaches for estimation of LIP flooding and their implementation in readily available simulation software packages (Prasad and Yuan 2020). NRC’s interest is in developing guidance for performing LIP PFHA at nuclear power plant (NPP) sites. This report describes the work done under Task 2 of the study that focused on describing uncertainties associated with LIP flood modeling. Two main types of uncertainties exist in LIP flood modeling: (1) aleatory uncertainties that arise from the inherent natural variability of the hydrometeorological system and (2) epistemic uncertainties that arise from the analysts’ incomplete knowledge of the hydrometeorological and site configuration. Principal sources for aleatory variability in LIP flood simulations at NPP sites include precipitation (i.e., magnitude, duration, and temporal distribution of LIP events), initial conditions (e.g., soil moisture content, stormwater drainage discharge, surface storage), boundary conditions (e.g., upstream discharge, downstream water surface elevations), and long-term temporal trends (e.g., climate change). This report describes approaches used for estimating the aleatory variability in precipitation including precipitation-frequency analysis, numerical weather prediction, and stochastic weather generation. Data sources that can be used to estimate aleatory variability in initial and boundary conditions are also described. Approaches to include effects of long-term trends like those from climate change into estimation of aleatory variability are summarized. Sources of epistemic uncertainty in LIP flood simulations at NPP sites include process representation (e.g., multiple approaches to represent runoff generation, stormwater drainage, and hydraulic routing), site configurations (e.g., site layout, flow features, status of temporary flood protection, blockage of drains), model resolution, and long-term temporal trends (e.g., known/planned site alterations, land-use changes at and in the vicinity of the site). The report describes alternative process representations (methods and models), particularly those implemented in the LIP flood simulation software packages reviewed in Task 1 report, and lists the associated model parameters. Approaches for estimating model parameters when surface and subsurface water exchanges occur are also described.

42 ENGINEERING↗

AUSTENITIC STAINLESS STEEL WELD EMBRITTLEMENT BY HYDROGEN AND TRITIUM

-Commercialization of fusion power requires large scale, local re processing of D/T fuel, with throughputs ~1/6 of current global tritium supply per plant per day -Tritium is a radioactive material that has special needs in handling, processing, and storage -Tritium’s effects on the safety performance of primary pressure boundary materials (likely austenitic stainless steel) will be important for regulatory acceptance

KRENTZ, TIMOTHY M.↗

Surface impacts and associated mechanisms of a moisture intrusion into the Arctic observed in mid-April 2020 during MOSAiC

Distinct events of warm and moist air intrusions (WAIs) from mid-latitudes have pronounced impacts on the Arctic climate system. We present a detailed analysis of a record-breaking WAI observed during the MOSAiC expedition in mid-April 2020. By combining Eulerian with Lagrangian frameworks and using simulations across different scales, we investigate aspects of air mass transformations via cloud processes and quantify related surface impacts. The WAI is characterized by two distinct pathways, Siberian and Atlantic. A moist static energy transport across the Arctic Circle above the climatological 90th percentile is found. Observations at research vessel Polarstern show a transition from radiatively clear to cloudy state with significant precipitation and a positive surface energy balance (SEB), i.e., surface warming. WAI air parcels reach Polarstern first near the tropopause, and only 1–2 days later at lower altitudes. In the 5 days prior to the event, latent heat release during cloud formation triggers maximum diabatic heating rates in excess of 20 K d -1 . For some poleward drifting air parcels, this facilitates strong ascent by up to 9 km. Based on model experiments, we explore the role of two key cloud-determining factors. First, we test the role moisture availability by reducing lateral moisture inflow during the WAI by 30%. This does not significantly affect the liquid water path, and therefore the SEB, in the central Arctic. The cause are counteracting mechanisms of cloud formation and precipitation along the trajectory. Second, we test the impact of increasing Cloud Condensation Nuclei concentrations from 10 to 1,000 cm -3 (pristine Arctic to highly polluted), which enhances cloud water content. Resulting stronger longwave cooling at cloud top makes entrainment more efficient and deepens the atmospheric boundary layer. Finally, we show the strongly positive effect of the WAI on the SEB. This is mainly driven by turbulent heat fluxes over the ocean, but radiation over sea ice. The WAI also contributes a large fraction to precipitation in the Arctic, reaching 30% of total precipitation in a 9-day period at the MOSAiC site. However, measured precipitation varies substantially between different platforms. Therefore, estimates of total precipitation are subject to considerable observational uncertainty.

54 ENVIRONMENTAL SCIENCES↗

Observations of wind farm wake recovery at an operating wind farm

Abstract. The interplay of momentum surrounding wind farms significantly influences wake recovery, affecting the speed at which wakes return to their freestream velocities. Under stable atmospheric conditions, wind farm wakes can extend over considerable distances, leading to sustained vertical momentum flux downstream, with variations observed throughout the diurnal cycle. Particularly in regions such as the US Great Plains, stable conditions can induce low-level jets (LLJs), impacting wind farm performance and power output. This study examines the implications of wake recovery using long-term observations of vertical momentum flux profiles across diverse atmospheric conditions. In these observations, several key findings were observed, such as (a) LLJ heights being altered downstream of a wind farm, especially when the LLJs are below 250 m above ground level; (b) a notable impact of LLJ height on wake recovery being observed using momentum flux profiles at upwind and downwind locations, wherein LLJs between 250 and 500 m above ground level resulted in larger momentum transfer within the wake (i.e., smaller velocity deficit) compared to LLJs below 250 m above ground level; (c) the largest momentum flux variability being observed during stable atmospheric conditions, with non-negligible variability observed during neutral and unstable atmospheric conditions; (d) detection of wake effects almost always being observed throughout the atmospheric boundary layer height; and finally (e) enhancement of wake recovery being observed in the presence of propagating gravity waves. These insights deepen our understanding of the intricate dynamics governing wake recovery in wind farms, advancing efforts to model and predict their behavior across varying atmospheric contexts. In addition, the performance of large-eddy-simulation-based semi-empirical internal boundary layer height model estimates incorporating real-world atmospheric and turbine inputs was evaluated using observations during LLJ conditions.

17 WIND ENERGY↗

Progress on the MARVEL Cybersecurity by Design Model-Based Systems Engineering Project

Formal model-based systems engineering (MBSE) combines a model, systems thinking, and systems engineering to visually depict the boundaries, context, and behavior of interconnected systems, facilitating effective design, development, and utilization of engineered systems throughout the systems engineering lifecycle. Although nuclear reactor vendors employ these tools to integrate functionality, performance, and safety, they are not yet addressing digital risk concerns introduced by use of operational technology, such as digital instrumentation and control systems. To accomplish this objective, the Microreactor Applications Research Validation and EvaLuation (MARVEL) microreactor was used as an MBSE case study. This real-world application provides a first-of-a-kind opportunity to demonstrate the benefits of integrating digital risk and cybersecurity into the MBSE design process of a nuclear reactor. This paper provides an update of the ongoing MARVEL Cyber MBSE project as it specifically relates to the integration of digital risk management and cybersecurity by design.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Grounding our Understanding of the Impacts of Boreal Forest Expansion on Shallow Cumulus Clouds with a Simple Modeling Framework

Abstract The expansion of the boreal forest poleward is a potentially important driver of feedbacks between the land surface and Arctic climate. A growing body of work has highlighted the importance of differences in evaporative resistance between different possible future Arctic land covers, which in turn alters humidity and cloudiness in the boundary layer, for these feedbacks. While thus far this problem has been studied primarily with complex Earth system models, we turn to a locally focused, idealized model capable of diagnosing and testing the sensitivity of first-order processes connecting vegetation, the atmospheric boundary layer, and low clouds in this critical region. This allows us to benchmark the mechanisms and results at the center of predictions from larger-scale simulations. A surface dominated by broadleaf trees, characterized by higher albedo and lower surface evaporative resistance, drives cooling and moistening of the boundary layer relative to a surface of needleleaf trees, characterized by lower albedo and higher surface evaporative resistance. Differences in evaporative resistance between these hypothetical Arctic vegetation covers are of equal importance to changes in albedo for the initial response of the boundary layer to boreal expansion, even with our idealized approach. However, compensation between the elevation of the lifting condensation level (LCL) and more rapid growth of the mixed layer over higher evaporative resistance surfaces can minimize changes in the favorability of shallow clouds over different land cover types under some conditions. We then perform two tests on the sensitivity of this compensating effect, to changes in water availability, represented first by a reduction in boundary layer humidity and then by both a reduction in humidity and soil moisture available to our vegetation surface. Finally, given the importance of this potential LCL–mixed-layer height compensation in our idealized modeling results, we look to determine its relevance in observational data from a field campaign in boreal Finland. These observations do confirm that such a coupling plays an important role in cumulus-topped boundary layers over a needleleaf forest surface. While our results confirm some underlying mechanisms at the center of prior work with Earth system models, they also provide motivation for future work to constrain the impact of boreal forest expansion. This will include both large eddy simulations to examine the impact of processes and feedbacks not resolved by a mixed-layer model, as well as a more systematic evaluation and comparison of relevant observations at the site in Finland and sites from prior boreal field campaigns. Significance Statement Clouds and vegetation are both important components of the climate system that interact across a range of scales. These interactions are central to understanding how changes at the land surface feedback on climate. For example, if a forest expands or recedes, diagnosing how that will impact clouds will determine whether you predict warming or cooling temperatures from that shift in the forest area. These predictions are often made with complex Earth system models, but we look to a more idealized representation of the land–atmosphere system to diagnose how shallow clouds should respond to changes in surface properties with different scenarios of boreal forest expansion at a more foundational level. This both grounds our understanding of previous analysis and provides helpful direction for future studies of this relevant and impactful land cover change.

Meteorology & Atmospheric Sciences↗

On the efficacy of post-build thermomechanical treatments to improve properties of Zirconium fabricated using ultrasonic additive manufacturing

Here, hot-isostatic pressing has been applied to Zirconium plate fabricated using the ultrasonic additive manufacturing (UAM) technique to enhance interfacial bond quality. Specimens heated to 800 C for 1 h at 100 MPa pressure showed grain growth across many prior foil-to-foil interfaces, thereby increasing foil adhesion. In addition to the material softening induced by a loss of Hall-Petch strengthening, premature failure of specimens loaded parallel to the build direction was observed. Premature specimen failure was attributed to the local delamination at prior foil-to-foil boundaries where grain growth was pinned by Ti impurities introduced during the UAM process via the interaction between a Ti-alloy buffer foil and the welded Zr-foils underneath. In addition, the presence of Ti along select foil interfaces resulted in the nucleation and growth of secondary (Zr,Ti)(Fe,Cr) laves phases during higher-temperature thermomechanical processing. Using a combination of micro-scale X-ray computed tomography, fractography, and in-situ digital image correlation, the effect of defects along prior foil-to-foil boundaries was revealed as delamination-assisted plasticity accelerated specimen failure in preferred tensile orientations. These findings underscore the importance of impurity control when optimizing weld quality of higher-strength material systems using ultrasonic welding.

36 MATERIALS SCIENCE↗

Regular and in-plane skyrmions and antiskyrmions from boundary instabilities

We formulate a theory of skyrmion and antiskyrmion generation using magnetic field and charge current pulses. We show that the topological defect can be created at an edge of a system with Dzyaloshinskii-Moriya interaction (DMI) as well as at a boundary between regions with different DMI. We consider both perpendicular and in-plane (also known as magnetic bimerons) versions of skyrmions and antiskyrmions. We show that the magnetization twist in the vicinity of an edge or a boundary is described by a kink solution, the presence of which can instigate the generation of topological defects. We study the collective excitations of magnetization analytically and numerically, and demonstrate that under application of magnetic field and charge current pulses the magnon modes localized near boundaries can develop instabilities leading to the formation of skyrmions or antiskyrmions. As a result, due to the skyrmion and antiskyrmion Hall effects, a properly chosen current direction can push the topological defects away from the boundary, thus facilitating their generation.

36 MATERIALS SCIENCE↗

Field-induced spin level crossings within a quasi-$XY$ antiferromagnetic state in $Ba$ 2 $FeSi$ 2 $O$ 7

Here, we present a high-field study of the strongly anisotropic easy-plane square lattice $S$ = 2 quantum magnet $Ba$ 2 $FeSi$ 2 $O$ 7 . This compound is a rare high-spin antiferromagnetic system with very strong easy-plane anisotropy, such that the interplay between spin level crossings and antiferromagnetic order can be studied. We observe a magnetic field-induced spin level crossing occurring within an ordered state. This spin level crossing appears to preserve the magnetic symmetry while producing a nonmonotonic dependence of the order parameter magnitude. The resulting temperature–magnetic field phase diagram exhibits two dome-shaped regions of magnetic order overlapping around 30 T. The ground state of the lower-field dome is predominantly a linear combination of |$S^z$ = 0$\rangle$ and |$S^z$ = 1$\rangle$ states, while the ground state of the higher-field dome can be approximated by a linear combination of |$S^z$ = 1$\rangle$ and |$S^z$ = 2$\rangle$ states. At 30 T, where the spin levels cross, the magnetization exhibits a slanted plateau, the magnetocaloric effect shows a broad hump, and the electric polarization shows a weak slope change. We determined the detailed magnetic phase boundaries and the spin level crossings using measurements of magnetization, electric polarization, and the magnetocaloric effect in pulsed magnetic fields to 60 T. Furthermore, we calculate these properties using a mean-field theory based on direct products of SU(5) coherent states and find good agreement. Finally, we measure and calculate the magnetically induced electric polarization that reflects magnetic ordering and spin level crossings. This multiferroic behavior provides another avenue for detecting phase boundaries and symmetry changes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multiscale, mechanistic calculation of the effective silver diffusion coefficient in polycrystalline silicon carbide: application to silver release in AGR-1 TRISO particles

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) fuel particles serves as a barrier to prevent the escape of fission from the fuel kernel. The release of silver (Ag) is a concern due to the long half-life of the 110mAg isotope. In this study, the effective diffusion coefficient of the fission product Ag through the grain boundary (GB) network is calculated using a combination of atomistic and phase-field methods. Atomistic calculations of Ag diffusivity in SiC bulk and GBs are leveraged to develop a mesoscale effective Ag diffusion coefficient (Deff) in SiC. Since GBs serve as pathways for Ag diffusion, Deff is defined as a function of temperature, microstructure variables, and fluence. Deff is implemented in the fuel performance code Bison to predict Ag release from AGR-1 TRISO fuel particles. We hereby quantify the impact of SiC grain size and irradiation on Ag release and improve Bison's predictions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗