Search NASA⌕ Search

SEARCH · Search NASA

Results for “PIN”

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 523 records · Page 29

Thermodynamic approach for enhancing superconducting critical current performance

Abstract The addition of artificial pinning centers has led to an impressive increase in the critical current density ( J c ) of superconductors, enabling record-breaking all-superconducting magnets and other applications. The J c of superconductors has reached ~0.2–0.3 J d , where J d is the depairing current density, and the numerical factor depends on the pinning optimization. By modifying λ and/or ξ, the penetration depth and coherence length, respectively, we can increase J d . For (Y 0.77 Gd 0.23 )Ba 2 Cu 3 O y ((Y,Gd)123), we can achieve this by controlling the carrier density, which is related to λ and ξ. We can also tune λ and ξ by controlling the chemical pressure in Fe-based superconductors, i.e., BaFe 2 (As 1− x P x ) 2 films. The variation in λ and ξ leads to an intrinsic improvement in J c via J d , allowing extremely high values of J c of 130 MA/cm 2 and 8.0 MA/cm 2 at 4.2 K, consistent with an enhancement in J d of a factor of 2 for both incoherent nanoparticle-doped (Y,Gd)123 coated conductors (CCs) and BaFe 2 (As 1− x P x ) 2 films, showing that this new material design is useful for achieving high critical current densities in a wide array of superconductors. The remarkably high vortex-pinning force in combination with this thermodynamic and pinning optimization route for the (Y,Gd)123 CCs reached ~3.17 TN/m 3 at 4.2 K and 18 T ( H || c ), the highest values ever reported for any superconductor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Coupled Monte Carlo Transport and Conjugate Heat Transfer for Wire-Wrapped Bundles Within the MOOSE Framework

Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics (CFD) to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. We first provide a brief introduction to Cardinal's software design, data mapping, and coupling strategy to highlight our approach to overcoming common challenges in high-fidelity multiphysics simulations. Here we then present two Cardinal simulations for hexagonal pin bundles. The first is a validation of Cardinal's conjugate heat transfer coupling of NekRS's Reynolds-Averaged Navier Stokes model with MOOSE's heat conduction physics for a bare seven-pin Freon-12 bundle flow experiment. Predictions for pin surface temperatures under three different heating modes agree reasonably well with experimental data and similar CFD modeling from the literature. The second simulation is a multiphysics coupling of OpenMC, NekRS, and BISON for a reduced-scale, seven-pin wire-wrapped version of an Advanced Burner Reactor bundle. Wire wraps are approximated using a momentum source model, and coupled predictions are provided for velocity, temperature, and power distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The nontrivial effects of annealing on superconducting properties of Nb single crystals

The effect of annealing on the superconducting properties of niobium single crystals was studied using optical, magnetic, and scanning tunneling microscopy (STM) methods. Pieces of the same crystal boule were studied before and after the annealing at 800 ${^\circ}\textrm{C}$, 1400 ${^\circ}\textrm{C}$, and near the melting point of niobium (2477 ${^\circ}\textrm{C}$). The initial samples had a high hydrogen content and low-temperature imaging revealed large hydrides (hundreds of micrometers) appearing below 190 K. The formation of these large precipitates is already completely suppressed by annealing at 800 ${^\circ}\textrm{C}$. However, the overall superconducting properties of the annealed samples did not improve and, in fact, worsened. In particular, the superconducting transition temperature decreased, the upper critical field increased, and the pinning strength increased. In the STM study, the sample was annealed initially at 400 ${^\circ}\textrm{C}$, measured, annealed at 1700 ${^\circ}\textrm{C}$, and measured again. The STM revealed a ‘dirty’ superconducting gap with a significant spatial variation in tunneling conductance after annealing at 400 ${^\circ}\textrm{C}$. The clean gap was recovered after annealing at 1700 ${^\circ}\textrm{C}$. This is likely due to oxygen redistribution near the surface, which is always covered by oxide layers in as-grown crystals. Our results indicate that vacuum annealing at least up to 1400 ${^\circ}\textrm{C}$, while removing a large percentage of hydrogen, introduces additional nanosized defects, likely hydride precipitates, that act as efficient pair-breaking and pinning centers. The dimensionless scattering rate is estimated to have increased from $\Gamma = 0.2$ to about $\Gamma = 0.4$ after annealing at 1400 ${^\circ}\textrm{C}$. These results on single crystals differ drastically from those obtained in polycrystalline bulk niobium (i.e. cut from superconducting radio-frequency cavities), where annealing is known to have a significant positive effect that is attributed to the improvement of the crystalline structure masking the more subtle influence of the hydrides.

43 PARTICLE ACCELERATORS↗

The effect of non-linear mutual friction on pulsar glitch sizes and rise times

ABSTRACT Observations of pulsar glitches have the potential to provide constraints on the dynamics of the high density interior of neutron stars. However, to do so, realistic glitch models must be constructed and compared to the data. We take a step towards this goal by testing non-linear models for the mutual friction force, which is responsible for the exchange of angular momentum between the neutron superfluid and the observable normal component in a glitch. In particular, we consider a non-linear dependence of the drag force on the relative velocity between superfluid vortices and the normal component, in which the contributions of both kelvin and phonon excitations are included. This non-linear model produces qualitatively new features, and is able to reproduce the observed bimodal distribution of glitch sizes in the pulsar population. The model also suggests that the differences in size distributions in individual pulsars may be due to the glitches being triggered in regions with different pinning strengths, as stronger pinning leads to higher vortex velocities and a qualitatively different mutual friction coupling with respect to the weak pinning case. Glitches in pulsars that appear to glitch quasi-periodically with similar sizes may thus be due to the same mechanisms as smaller events in pulsars that have no preferred glitch size, but simply originate in stronger pinning regions, possibly in the core of the star.

79 ASTRONOMY AND ASTROPHYSICS↗

Electronic noise of a single skyrmion

To enable the practical use of skyrmion-based devices, it is essential to achieve a balance between energy efficiency and thermal stability while also ensuring reliable electrical detection against noise. Understanding how a skyrmion interacts with material disorder and external perturbations is thus essential. Here, we investigate the electronic noise of a single skyrmion under the influence of thermal fluctuations and spin currents in a magnetic thin film. We detect the thermally induced noise with a 1/ƒ $\gamma$ signature in the strong pinning regime but a random telegraph noise in the intermediate pinning regime. Both the thermally dominated and current-induced telegraphlike signals are detected in the weak pinning regime. Our results provide a comprehensive electronic noise picture of a single skyrmion, demonstrating the potential of noise fluctuation as a valuable tool for characterizing the pinning condition of a skyrmion. In conclusion, these insights could also aid in the development of low-noise and reliable skyrmion-based devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Anomalies in the electromagnetic response and peculiarities of the domain structure in the magnetic Weyl semimetal Co 3⁢ Sn 2 ⁢S 2

In this paper, we present a comprehensive study of magnetic properties of Co 3 ⁢Sn 2 ⁢S 2 (CoSnS) crystals using macroscopic characterization and observations of the domain structure (DS) in a wide range of easy- and hard-axis fields (𝐻) and temperatures (𝑇). From macroscopic measurements, we derive temperature dependences of the basic magnetic parameters of the material [saturation magnetization 𝑀 𝑠 ⁡(𝑇), uniaxial anisotropy constant 𝐾⁡(𝑇), and exchange stiffness 𝐴⁡(𝑇)] and, based on them, calculate the thermodynamic magnetization 𝑀⁡(𝐻,𝑇) curves and DS parameters for pinning-free samples. We find important qualitative and quantitative differences between the theory and experiment defined by the efficient pinning of domain walls (DWs) in our samples. It turns out that, already at temperatures close to the ferromagnetic transition 𝑇 𝑐 , the DW dynamics falls into the heavy creep regime, resulting in a strong deferral of the DS adjustment into the equilibrium state. Here, we suggest that specific magnetic defects, responsible for pinning and causing field-tuned reproducible domain nucleation patterns observed in our high-quality samples, are pointlike. Recently found magnetic polarons formed around vacancies and impurities, which carry strong diamagnetic moments of spin-orbit nature and are intrinsic for CoSnS and expected in other Weyl semimetals, are a possible source of pinning yielding peculiar 𝑀⁡(𝐻,𝑇) response of CoSnS.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Continuous and discontinuous transitions in the depinning of two-dimensional dusty plasmas on a one-dimensional periodic substrate

We report Langevin dynamical simulations are performed to study the depinning dynamics of two-dimensional dusty plasmas on a one-dimensional periodic substrate. From the diagnostics of the sixfold coordinated particles P 6 and the collective drift velocity V x , three different states appear, which are the pinning, disordered plastic flow, and moving ordered states. It is found that the depth of the substrate is able to modulate the properties of the depinning phase transition, based on the results of P 6 and V x , as well as the observation of hysteresis of V x while increasing and decreasing the driving force monotonically. When the depth of the substrate is shallow, there are two continuous phase transitions. When the potential well depth slightly increases, the phase transition from the pinned to the disordered plastic flow states is continuous; however, the phase transition from the disordered plastic flow to the moving ordered states is discontinuous. When the substrate is even deeper, the phase transition from the pinned to the disordered plastic flow states changes to discontinuous. When the depth of the substrate further increases, as the driving force increases, the pinned state changes to the moving ordered state directly, so that the disordered plastic flow state disappears completely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advancing Temperature Monitoring of the Bottom Anode in a Direct Current Electric Arc Furnace Operations with Distributed Optical Fiber Sensors.

The bottom anode in the Direct Current Electric Arc Furnace (DC EAF) is critical for completing the electrical circuit necessary for sustaining the arc within the furnace. For pin-type bottom anodes, monitoring of the temperature of select pins instrumented with thermocouples is performed to track bottom wear in the EAF and inform the operator when the furnace should be removed from service. Furthermore, this work presents the results from a plant trial using distributed temperature monitoring of bottom anode pins in a 165-ton DC EAF over a two-month service period utilizing two optical fiber sensing techniques: fiber Bragg grating (FBG) and Rayleigh backscattering (RBS). The early detection of temperature anomalies along the length of the anode pin through distributed sensing enhances operational safety, providing a robust alternative to traditional thermocouples.

42 ENGINEERING↗

Load Profile Inpainting for Missing Load Data Restoration and Baseline Estimation

This paper introduces a Generative Adversarial Nets (GAN) based, Load Profile Inpainting Network (Load-PIN) for restoring missing load data segments and estimating the baseline for a demand response event. The inputs are time series load data before and after the inpainting period together with explanatory variables (e.g., weather data). Here, we propose a Generator structure consisting of a coarse network and a fine-tuning network. The coarse network provides an initial estimation of the data segment in the inpainting period. The fine-tuning network consists of self-attention blocks and gated convolution layers for adjusting the initial estimations. Loss functions are specially designed for the fine-tuning and the discriminator networks to enhance both the point-to-point accuracy and realisticness of the results. We test the Load-PIN on three real-world data sets for two applications: patching missing data and deriving baselines of conservation voltage reduction (CVR) events. We benchmark the performance of Load-PIN with five existing deep-learning methods. Our simulation results show that, compared with the state-of-the-art methods, Load-PIN can handle varying-length missing data events and achieve 15-30% accuracy improvement.

14 SOLAR ENERGY↗

Improved charge extraction in inverted perovskite solar cells with dual-site-binding ligands

Inverted (pin) perovskite solar cells (PSCs) afford improved operating stability in comparison to their nip counterparts but have lagged in power conversion efficiency (PCE). The energetic losses responsible for this PCE deficit in pin PSCs occur primarily at the interfaces between the perovskite and the charge-transport layers. Additive and surface treatments that use passivating ligands usually bind to a single active binding site: This dense packing of electrically resistive passivants perpendicular to the surface may limit the fill factor in pin PSCs. We identified ligands that bind two neighboring lead(II) ion (Pb 2+ ) defect sites in a planar ligand orientation on the perovskite. We fabricated pin PSCs and report a certified quasi–steady state PCE of 26.15 and 24.74% for 0.05– and 1.04–square centimeter illuminated areas, respectively. The devices retain 95% of their initial PCE after 1200 hours of continuous 1 sun maximum power point operation at 65°C.

Science & Technology - Other Topics↗

Demonstration of MOOSE-based Griffin reactor physics, code for heterogeneous lead-cooled fast reactor analysis

The MOOSE-based reactor physics code Griffin was assessed on a heterogeneous pin-resolved model of a prototype lead-cooled fast reactor assembly. This model was developed in preparation for future use in MOOSE-based multiphysics calculations for computing hot channel factors. Heterogeneous multigroup cross sections were prepared using the fast reactor multi-group cross section processing code MC{sup 2}-3 using a two-step method. Griffin simulations were performed using the DFEM-SN solver on 576 cores on Argonne's LCRC cluster. Diffusion-based acceleration methods were applied (NDA and CMFD). Reference solutions were generated with continuous energy MCNP and the hybrid MOC/finite element solver PROTEUS-MOC for code-to-code comparison. Space-angle convergence studies were conducted to observe convergence in k-eigenvalue and axial pin power distributions. The fully resolved Griffin calculation was within 68 pcm of the MCNP eigenvalue and exhibited max 1.2% relative error in the axial pin power distribution. Griffin produced nearly identical results to PROTEUS-MOC when using the same 9-group multigroup cross-section set. Griffin demonstrated favorable scaling in wall-clock time and memory usage when using diffusion-based acceleration methods. Griffin is capable of simulating the pin-resolved heterogeneous LFR assembly with good accuracy and performance, and is suitable for future use in coupled high-fidelity hot channel factor simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Recent progress on the APOLLO3 calculation scheme for the FLUOLE2 experiment

Studies of the neutron fluence on reactor pressure vessel and the correlated reactor lifetime issues require very good knowledge of the neutron source distribution in the reactor core and especially in the peripheral region. The FLUOLE2 experiment offers an appropriate starting configuration to investigate the possibility of using the APOLLO3 deterministic code with a two-step calculation scheme. It enables us to calculate a pin-by-pin neutron source distribution in the core, then to carry out the neutron propagation from the core to neutron detectors. In this paper, we focus on a calculation scheme dedicated to the determination of core neutron source distribution using APOLLO3 with the SHEM281 energy mesh, a Method of Characteristics (MOC) to solve the transport equation for the lattice computation step and a 20-groups 3D computation using the MINARET SN solver of APOLLO3. Reactivity and pin-by-pin reaction rate results have been analyzed. Results obtained with both lattice and core computations show a good agreement with those produced through a reference calculation using the TRIPOLI-4 Monte Carlo transport code. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Application of Gaussian Bayes classifier to differentiate chlorine-based chemical agents

The Portable Isotopic Neutron Spectroscopy (PINS) is a commercialized system developed by Idaho National Laboratory (INL) to examine chemical warfare agents (CWA) non-destructively, utilizing Prompt Gamma Neutron Activation Analysis (PGNAA) techniques. The PINS system takes advantage of a high-resolution gamma-ray spectrum from a mechanically-cooled high-purity germanium (HPGe) detector. One of the difficult technical challenges is to discriminate the chlorine-based chemical agents. Especially, CN, CNB, CNS and CG have similar chemical compositions to make it hard to discriminate them with a higher confidence. Current identification algorithms for PINS systems with 252-Cf sources have been improved and updated continuously as more field data became available, and new algorithms was studied to complement the current algorithms by adopting the Gaussian Bayes classifier. These new algorithms were intended to be applied to a subset of chlorine-based chemical agents, and their main goal is discriminate CN, CNB, CNS and CG with their ratios of the chlorine neutron inelastic 1763keV peak to the chlorine thermal neutron capture 1959keV peak, which is referred to as the “Cl i/c” or “clic” ratio in this study. The Cl i/c ratios were assumed to follow Gaussian distributions with the means and the standard deviations unique to their corresponding chemical agents. The prior probabilities of these four chemical agents were optimized with a collection of field data to achieve the best performance in terms of precision or positive predictive value (PPV). Finally, their posterior probabilities as functions of the Cl i/c ratio were implemented in the current version of PINS analysis software in order to be tested with more field data.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Final Fugitive Dust Control NOV Causal Report

On June 30, 2020, a Notice of Violation (NOV) was issued by the City of Albuquerque (COA) Environmental Health Department, Air Quality Program. The NOV identified two violations of New Mexico Administrative Code (NMAC) 20.11.20, Fugitive Dust Control, stemming from an August 30, 2019 inspection of the construction site at Sandia/New Mexico (SNL/NM) Building 812. After the August 30 inspection, a Post-Instruction Notification (PIN) was issued to the SNL Construction Facilities Manager. The PIN was acknowledged by Department 4722 and sent to the COA on September 13, 2020 via email. The PIN "Comply by" date was September 13, 2019 was transmitted by National Nuclear Security Administration/ Sandia Field Office (NNSA/SFO) to COA offices on November 12, 2019. The PIN response was 48 business days past due for various reasons which were explored during the causal analysis.

54 ENVIRONMENTAL SCIENCES↗

Tomographic Image Reconstruction for the Parallel-Slit Ring Collimator Fast Neutron Emission Tomography System

For the past three years, Oak Ridge National Laboratory (ORNL) has been developing a passive fast-neutron emission tomography capability. The goal of this development is the ability to quantify the neutron source strength of individual fuel pins (rods) in spent nuclear fuel assemblies. Such a system could be used to measure the burnup of each fuel pin in a spent fuel assembly in order to take burnup credit when loading dry storage casks or to count individual fuel pins in spent fuel assemblies for safeguards purposes. At present, a laboratory prototype imager is under construction. The purpose of this prototype is to demonstrate imaging capability sufficient to resolve individual fuel pins in spent fuel assemblies. This report documents the development of the iterative reconstruction code used to perform tomographic image reconstruction, the imager response calculation used by the reconstruction code, and the results of reconstructions of simulated tomographic imaging measurements for the prototype imager design.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Assessment of the Griffin Reactor Multiphysics Application Using the Empire Micro Reactor Design Concept

In late 2019, INL and ANL agreed to jointly develop the reactor physics code named Griffin based on the integration of the two code suites, MAMMOTH/Rattlesnake (INL) and MC2 - 3/PROTEUS (ANL). Griffin is being developed based on the MOOSE framework and MOOSE quality assurance procedures. This decision was made to be able to allow DOE-NE to efficiently invest funding to this area and to provide effective and timely support for existing and potential users; the latter includes industry and government organizations who are developing various types of advanced reactors in the near and long term. Since MAMMOTH/Rattlesnake has been developed based on the MOOSE framework, the INL/ANL Griffin development team agreed to build Griffin beginning with a merger of MAMMOTH and Rattlesnake into a single code and moving forward by implementing capabilities from the PROTEUS suite into Griffin. Moving forward, both ANL and INL efforts are equally invested in the Griffin project, with management support, to provide an advanced reactor multiphysics tool to assist in reactor design, optimization, and safety analysis. Much work remains in moving Griffin forward to migrate PROTEUS capabilities and to optimize performance to meet user needs. The main objective of this work is to assess the current status of Griffin capabilities in terms of performance and accuracy, to determine priorities for PROTEUS migration, and to identify capabilities and features to improve for supporting the code integration effort. For this assessment, the Empire micro reactor problem that was developed in the ARPA-E MEITNER program was selected as an advance reactor concept of interest to the technical community. The Empire reactor problem was expanded from its original incomplete specification to be a small heat-pipe-cooled micro reactor core with ~113 cm radius and 70 cm in height, composed of 18 fuel assemblies, 12 control drums, and beryllium radial and axial reflectors. In the current model, using 5 cm axial reflectors specified in the original Empire assembly model, more than 10% of neutrons leak axially and through the empty center safety hole, as well as through heat pipe channels in fuel assembly elements that extend through the top reflector region. Several calculation models of the core were defined for systematic assessment, including 2-D and 3-D fuel assemblies and whole cores with cylindrical boundaries. Cross sections were generated using Serpent 2, and meshes were produced using the Argonne mesh tool or the INL neutronics meshing tools combined with CUBIT. Cross sections and meshes were converted to the ISOXML and Exodus formats, respectively, so that Griffin and PROTEUS could use consistent data for solving the reactor problems. With the prepared cross sections and meshes, PROTEUS was run first to ensure that all input data were correctly generated and input options in terms of angle, mesh, and energy group were accurately determined. Comparisons against Serpent 2 solutions were made in terms of eigenvalue and pin power. The same calculations and comparisons were then conducted using Griffin. For the fuel assembly and whole core problems, the PROTEUS eigenvalues agreed well with reference Serpent 2 solutions within 100 and 30 pcm, respectively, and pin power differences relative to Serpent 2 were overall less than 2.2% and RMS 0.8% for the whole core models. This indicated that all input data were properly prepared. Using the same data, Griffin was run selecting the SAAF-CFEM SN solver with Legendre-Gaussian quadrature and NDA and DSA for acceleration. It was found that the SAAF-CFEM solver of Griffin required finer meshes to achieve eigenvalue and pin power solutions in good agreement with Serpent 2, consequently requiring more memory requirement and longer computation time. On the other hand, the SPH-Diffusion 2-D core calculations performed using Griffin were able to recover the exact eigenvalue from the reference Serpent 2 solutions, resulting in a pin-power distribution with an RMS of 0.6% and maximum absolute difference of less than 1.4%. The runtimes for SPH-Diffusion for the 2-D core were less than 3 minutes on 40 cores. During this evolution of this evaluation, many updates were made in Griffin by the Griffin development team of INL (focusing on software updates) and ANL (reviewing and supporting software updates) to complete this assessment. Observations from the code assessment are presented in the conclusion section of this report, followed by a discussion of recommendations for future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Analysis of Approximations in Modeling of BWR Bundle Void Distributions

In boiling water reactors, complex heterogeneous bundle designs, control blades adjacent to the corner of bundles, and the presence of boiling can lead to complex internal void distributions. A few approximations exist to model these void distributions. They could be modeled using a 1D axial solver in which each axial node is assumed to be at an average void, or each pin cell could be modeled with its own void concentration. In the latter case, the void could be discretized in pin-centered or coolant-centered channels. The goal of this project was to quantify the effect of using the different approximations for modeling internal void distributions on neutronics calculations. Using 3D void distributions calculated with CTF, Monte Carlo Neutral Particle (MCNP) transport code models were created for GE-9 and GE-14 lattices. For each model, the internal void distribution from CTF at a given axial node was selected, and a lattice calculation was carried out with MCNP. Comparisons between models using a lattice-averaged void, or using a void distribution in coolant-centered channels, showed large differences in reactivity which in some cases were well above 1,000 pcm, and it also showed differences in normalized fission rates greater than 20%. It was also found that using a lattice average void can lead to a significant difference in the estimation of the worth of a control blade. The differences found when comparing results from models using pin-centered and coolant-centered channels were up to 200 pcm in reactivity and up to 1.4% in the normalized fission rates. In addition to these two sets of comparisons, MCNP models were set up so that each subchannel had a saturated liquid component around the fuel pins and a saturated vapor component in the center to approximate annular flow. In comparison to the models using coolant-centered subchannels, up to 1–3% differences in normalized fission rates could be found.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗