Search NASA⌕ Search

SEARCH · Search NASA

Results for “Heating elements”

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 901 records · Page 50

Metallic Fuels Handbook, Part 1: Alloys Based on U-Zr, Pu-Zr, U-Pu, or U-Pu-Zr, Including Those with Minor Actinides (Np, Am, Cm), Rare-earth Elements (La, Ce, Pr, Nd, Gd), and Y; and Part 2: Elements and Alloys not Based on U-Zr, Pu-Zr, U-Pu, or U-Pu-Zr

The Metallic Fuels Handbook summarizes currently available information about phases and phase diagrams, heat capacity, thermal expansion, and thermal conductivity of elements and alloys in the U-Pu-Zr-Np-Am-La-Ce-Pr-Nd system. Although many sections are reviews and updates of material in previous versions of the Handbook, this revision is the first to include alloys with four or more elements. In addition to presenting information about materials properties, the handbook attempts to provide information about how well each property is known and how much variation exists between measurements. Although it includes some results from models, its primary focus is experimental data. Because of its size, the Handbook has been split into two parts, which are presented as separate documents. This (the first) part contains introductory material and information about the U-Pu-Zr system (U-Zr, Pu-Zr, U-Pu, and U-Pu-Zr alloys, including those with minor actinides, rare-earth elements, and Y). The second part contains information about elements and the remaining alloys in the U-Np-Pu-Am-La-Ce-Pr-Nd-Zr system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Preflight reentry heat transfer analysis of Space Shuttle

Preflight predictions of the structural temperature distributions during entry are compared with data from the initial Shuttle flight. Finite element thermal analysis programming was used to model the heat flow on Shuttle structures and actual gas properties of air were employed in the analyses of aerodynamic heating. Laminar, separated, and turbulent heat fluxes were calculated for varying locations on the craft using velocity-attitude and angle-of-attack projections taken from the nominal STS-1 trajectory. Temperature time histories of the first flight are compared with laminar and turbulent flow assumptions and an unpredicted rapid cooling 1800 sec into entry is credited to inaccurate assumptions of structural heat dissipative properties or flow conditions in that time phase of the flight; additional discrepancies in descriptions of heating of the upper fuselage are attributed to a lack of knowledge of the complex flow patterns existing over that area of the Shuttle body.

Ko, W. L.↗

A Numerical Tool for the Design of Enhanced Geothermal Reservoirs based on Modeling of Thermo-Hydro-Mechanical Processes and Explicit Representation of 3D Discrete Fracture Networks

This project has developed a coupled thermo-hydro-mechanical numerical modeling tool for prediction of energy production from Enhanced Geothermal Reservoirs (EGS). The tool is unique because it can explicitly represent thousands of pre-existing natural fractures. The capabilities developed in this project were built on the already existing expertise of Itasca arising from thirty years of research and development in discrete element modeling of fractured rocks. A thermal module suitable for analysis of different heat transfer mechanisms in fractured rocks percolated by fluids was developed and coupled with the already existing hydro-mechanical calculation in 3DEC, a commercial software for analysis of jointed rock based on the discrete element method (DEM). The heat transfer module simulates heat conduction, heat advection by moving fluid in rock fractures, and convective heat exchange at the fracture walls between the rock and the moving fluid. Techniques for spatial discretization of heat transfer in EGS-type problems were evaluated, and a boundary layer meshing approach was developed. Also, given the significant differences in time scale of physical processes involved, temporal integration schemes for coupled advection-diffusion processes and coupled thermal-hydro-mechanical processes were developed by taking advantage of both implicit and explicit time integration and also optimizing the coupling intervals. These improvements were necessary for modeling EGS reservoirs. Also, the DFN capacities in 3DEC were extended to include a larger number of fractures. Research was carried out on approaches to simplify the DFN so that the model can meet computational demands, while preserving thermo-hydro-mechanical characteristics of the DFN. The developed modules were coupled with the already existing capabilities of 3DEC. The results were presented as a submitted conference paper (Riahi et al., 2018a). A series of benchmark and test examples showed the developed capabilities. Finally, the developed tool was tested for a large-scale reservoir model for longer periods of injection and production. The thermal draw-down curves were generated and presented. The project successfully met all the milestones laid out in the proposal. Results demonstrated the feasibility of developing a fast and reliable software for the EGS community to use for assessment and prediction of reservoir performance. Some of its applications include assessment of stimulated volume and produced temperature and power. It can also be used for improved engineering, design and optimization of EGS reservoirs as well as other unconventional resources such as shale gas. Finally, the 3D analysis presented in this study showed that simulating EGS reservoirs with tens of thousands of fractures is feasible, and confirmed that the developed tool can serve toward meeting the goals and objectives of Geothermal Technologies Office in overcoming technical barriers associated with commercial EGS development.

15 GEOTHERMAL ENERGY↗

On Heating the Sun's Corona by Magnetic Explosions: Feasibility in Active Regions and prospects for Quiet Regions and Coronal Holes

We build a case for the persistent strong coronal heating in active regions and the pervasive quasi-steady heating of the corona in quiet regions and coronal holes being driven in basically the same way as the intense transient heating in solar flares: by explosions of sheared magnetic fields in the cores of initially closed bipoles. We begin by summarizing the observational case for exploding sheared core fields being the drivers of a wide variety of flare events, with and without coronal mass ejections. We conclude that the arrangement of an event's flare heating, whether there is a coronal mass ejection, and the time and place of the ejection relative to the flare heating are all largely determined by four elements of the form and action of the magnetic field: (1) the arrangement of the impacted, interacting bipoles participating in the event, (2) which of these bipoles are active (have sheared core fields that explode) and which are passive (are heated by injection from impacted active bipoles), (3) which core field explodes first, and (4) which core-field explosions are confined within the closed field of their bipoles and which ejectively open their bipoles. We then apply this magnetic-configuration framework for flare heating to the strong coronal heating observed by the Yohkoh Soft X-ray Telescope in an active region with strongly sheared core fields observed by the MSFC vector magnetograph. All of the strong coronal heating is in continually microflaring sheared core fields or in extended loops rooted against the active core fields. Thus, the strong heating occurs in field configurations consistent with the heating being driven by frequent core-field explosions that are smaller but similar to those in confined flares and flaring arches. From analysis of the thermal and magnetic energetics of two selected core-field microflares and a bright extended loop, we find that (1) it is energetically feasible for the sheared core fields to drive all of the coronal heating in the active region via a staccato of magnetic microexplosions, (2) the microflares at the feet of the extended loop behave as the flares at the feet of flaring arches in that more coronal heating is driven within the active bipole than in the extended loop, (3) the filling factor of the X-ray plasma in the core field microflares and in the extended loop is approximately 0.1, and (4) to release enough magnetic energy for a typical microflare (10^27 - 10^28 erg), a microflaring strand of sheared core field need expand and/or untwist by only a few percent at most. Finally, we point out that (1) the field configurations for strong coronal heating in our example active region (i.e., neutral-line core fields, many embedded in the feet of extended loops) are present in abundance in the magnetic network in quiet regions and coronal holes, and (2) it is known that many network bipoles do microflare and that many produce detectable coronal heating. We therefore propose that exploding sheared core fields are the drivers of most of the heating and dynamics of the solar atmosphere, ranging from the largest and most powerful coronal mass ejections and flares, to the vigorous microflaring and coronal heating in active regions, to the multitude of fine-scale explosive events in the magnetic network. The low-lysing exploding core fields in the network drive microflares, spicules, global coronal heating, and ,consequently, the solar wind.

Moore, Ronald L.↗

Design and Engineering of LUPIN: A Test-Bed Radio-Frequency Ion Source for Enhanced Neutral Beam Injection on DIII-D

The Large, Uniform Plasma for Ionizing Neutrals (LUPIN) is a radio-frequency (RF) inductively coupled plasma (ICP) chamber for demonstrating plasma performance of an RF ICP positive ion source upgrade for the DIII-D neutral beam injection (NBI) system. LUPIN will be used to investigate ion source physics, including neutral gas dynamics, plasma density uniformity, interactions with Faraday shields, and power coupling to novel RF antenna designs. LUPIN has an RF generator capable of delivering 20 kW of power at 2 MHz, which is coupled into a cylindrical quartz vessel measuring 20 cm in length and 10 cm in radius. This configuration matches the power density requirements for a full-scale ion source. Target hydrogen and deuterium plasma densities exceeding 10 18 m -3 would relate to extracted ion current densities of 2100 A/m 2 for 10s. Vacuum conductance and gas flow calculations predict a maximum achievable neutral gas flow rate of 1675 Pa ⋅ L/s at 5 Pa of He, which mimics the gas flow of the DIII-D NBI system. Designs have been developed for an internal Faraday shield to mitigate heat flux and ion sputtering on the dielectric vessel. Thermomechanical finite element simulations demonstrated the Faraday shield design to be capable of withstanding anticipated heat loads from worst-case operation scenarios. Finally, results of upcoming experimental investigations on LUPIN will guide the design of a full-scale prototype for DIII-D integration.

Faraday shield↗

Design and Engineering of LUPIN: A Test-Bed Radiofrequency Ion Source for Enhanced Neutral Beam Injection on DIII-D

The Large, Uniform Plasma for Ionizing Neutrals (LUPIN) is an RF inductively coupled plasma (ICP) chamber for demonstrating plasma performance of an RF ICP positive ion source upgrade for the DIII-D Neutral Beam Injection (NBI) system. LUPIN will be used to investigate ion source physics, including neutral gas dynamics, plasma density uniformity, interactions with Faraday shields, and power coupling to novel RF antenna designs. LUPIN has an RF generator capable of delivering 20 kW of power at 2 MHz, which is coupled into a cylindrical quartz vessel measuring 20 cm in length and 10 cm in radius. This configuration matches the power density requirements for a full-scale ion source. Target hydrogen and deuterium plasma densities exceeding 10 18 m -3 would relate to extracted ion current densities of 2100 A/m 2 for 10 s. Vacuum conductance and gas flow calculations predict a maximum achievable neutral gas flow rate of 1675 Pa·L/s at 5 Pa of He, which mimics the gas flow of the DIII-D NBI system. Designs have been developed for an internal Faraday shield to mitigate heat flux and ion sputtering on the dielectric vessel. Thermomechanical finite element simulations demonstrated the Faraday shield design to be capable of withstanding anticipated heat loads from worst-case operation scenarios. Results of upcoming experimental investigations on LUPIN will guide the design of a full-scale prototype for DIII-D integration.

Faraday shield↗

Heat Emission and Temperature Non-uniformity in Rod Bundles with Different Spacing, Cooled by Heavy Liquid-Metal Coolant

Heat emission and temperature non-uniformity along the perimeter of rods in free lattices of fuel elements with heavy liquid-metal coolant are studied. The experimental data for an assembly of smooth rods, assemblies of fuel elements with spiral ribbing, and transverse lattices are analyzed. The studies show that there is almost no temperature non-uniformity in a lattice of smooth fuel elements along its perimeter, which is in contrast to the high general azimuthal temperature non-uniformity in an assembly with spacing by spiral coils. Rib spacing resulted in significant reduction of heat emission. In an assembly with transverse spacer lattices, the heat emission increases only in the region containing the lattices; in the spaces between lattices, it is approximately equal to the heat emission in a bundle with smooth fuel elements. Inside a spacer lattice, periodic temperature non-uniformities appear along the perimeter of the fuel elements; these non-uniformities are due to the fuel elements touching the elements of the lattice. A relation is recommended for calculating heat emission (Nusselt numbers) and temperature non-uniformity along the perimeter of the fuel elements for the studied methods of spacing the fuel elements.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Bi‐Continuous W‐Rich Refractory High Entropy Alloy‐Cu Composite: Toward Material Innovation of Nuclear Reactor Coolant System

Abstract Refractory high‐entropy alloys (RHEAs) are considered promising candidate materials for next‐generation nuclear reactors due to their superior mechanical strength, irradiation resistance, and thermal stability at high temperatures. However, the significant positive heat of mixing between refractory alloying elements and Cu, commonly used in cooling systems, poses challenges in forming composite structures. This study addresses the issue using a liquid metal dealloying (LMD) process. A precursor alloy (WTaVTi) with a directional dendrite‐interdendrite structure is fabricated and reacted with molten Cu at 1200 °C for 96 h. This approach produced a RHEA‐Cu composite with a stable interface between RHEA (W 31.5 Ta 30.9 V 21.4 Ti 14.3 ) and Cu, featuring a spontaneously formed W‐rich interlayer that enhances interfacial bonding. The composite showed excellent irradiation resistance, with 30% less swelling under α‐ion irradiation than pure W. It also exhibited low thermal conductivity at room temperature, but reached ≈120 W m −1 ·K −1 at ≈650 °C, surpassing pure W. This temperature‐dependent rise in κ, with a positive gradient of +0.075 W m −1 ·K − 2 , is attributed to decreasing diffuse mismatch at elevated temperatures. The large‐scale reaction and stable microstructure achieved through LMD process highlight its industrial potential. This work offers a strategy for developing high‐performance materials by combining RHEA's radiation resistance with Cu's thermal conductivity for extreme environments.

Chemistry↗

Plastic Laminate Antireflective Coatings for Millimeter-Wave Optics in BICEP Array

The BICEP/Keck series of experiments target the cosmic microwave background at degree-scale resolution from the South Pole. Over the next few years, the “Stage-3” BICEP Array (BA) telescope will improve the program’s frequency coverage and sensitivity to primordial B-mode polarization by an order of magnitude. The first receiver in the array, BA1, began observing at 30/40 GHz in early 2020. The next two receivers, BA2 and BA3, are currently being assembled and will map the southern sky at frequencies ranging from 95 to 150 GHz. Common to all BA receivers is a refractive, on-axis, cryogenic optical design that focuses microwave radiation onto a focal plane populated with antenna-coupled bolometers. High-performance antireflective coatings up to 760 mm in aperture are needed for each element in the optical chain, and must withstand repeated thermal cycles down to 4 K. Here, in this work, we present the design and fabrication of the 30/40 GHz anti-reflection coatings for the recently deployed BA1 receiver, with indices matched to its various polyethylene, nylon and alumina optical components. We describe an epoxy coating technique designed for alumina optics, which achieves better than 80% transmission at room temperature. For polyethylene optical elements, we present a new heat-compression approach that allows low-density polytetrafluoroethylene AR layers to reach sub-percent reflected power. We describe the planned use of these methods for the next BA cryostats, which may inform technological choices for future small-aperture telescopes of the CMB-S4 experiment.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermal Wear Sensing System: Proof of the Concept

A technique for wear sensing using a sacrificial layer with specific material properties has been developed successfully. This is achieved by recording the temperature evolution underneath the sacrificial organic layer during sliding tests, and by correlating the recorded temperature with a wear depth. Finite element modeling was used to calculate the contact temperature which spreads through the sacrificial layer under different sliding test conditions. Numerical simulation was used to determine the optimal parameters for the best correlation between temperature and wear depth. These parameters were used to select the appropriate sacrificial layer for wear sensing. Thus, an epoxy layer was implemented and experimental measurements were carried out to check the correlation between the measured temperature underneath the epoxy layer and wear depth. It was found that a low thermal conductivity of the sacrificial layer of less than 1 W/m K and a relative high thermal conductivity of the substrate, namely more than 4 W/m K are the optimal parameters for our thermal wear sensing system. The thickness of the sacrificial layer is also a crucial parameter and has to be taken into account when designing such a sensing system.

36 MATERIALS SCIENCE↗

Tailoring Microstructure and Properties of a Superelastic Ti–Ta Alloy by Incorporating Spark Plasma Sintering with Thermomechanical Processing

A powder metallurgy Ti–35 at.% Ta alloy was produced using a spark plasma sintering process. The sintered alloy exhibited Ta-rich regions at a sintering temperature of 1100 °C. Complete diffusion of elemental powders was achieved via homogenization heat treatment at 1500 °C for 12 h, but the ductility declined drastically. The brittle behavior of the homogenized sample was due to the high-oxygen content, continuous grain boundary α-phase and high-angle grain boundaries. Excellent mechanical properties were achieved when hot forging was used along with preheating the homogenized specimen at 1100 °C. The mechanisms underlying this phenomenon were suppression of the solid grain boundary α-phase, high fraction of low-angle grain boundaries and dynamic recrystallization. Tensile strength improved to the maximum value of 868 MPa. The strain recovery ratio increased to 90% in the first cycle during a cyclic loading–unloading test. Incomplete strain recovery of the first cycle was due to the presence of the martensite phase (α″).

36 MATERIALS SCIENCE↗

Microstructure characteristics of LPBF&HIP fabricated graded composite transition joint between ferritic steel and austenitic stainless steel

Graded composite transition joints (GCTJs) offer a promising alternative to conventional dissimilar metal welds (DMWs) by enabling smooth compositional and microstructural transitions. However, GCTJs fabricated solely through additive manufacturing (AM) face challenges such as heat accumulation, complex parameter control, and elemental segregation. In this study, we propose a novel approach that relies on AM to design a spatially graded structure in one alloy and then employs hot isostatic pressing (HIP) as a diffusion bonding method to join it with a second alloy. Here, this method combines the flexibility of AM with the powder net-shaping advantage of HIP. Specifically, a series of closely packed austenitic stainless steel 304 conical structures were printed using laser powder bed fusion (LPBF) and then combined with ferritic steel P91 powder via HIP. By using electron backscatter diffraction (EBSD), electron probe microanalysis (EPMA), and transmission electron microscopy (TEM) techniques, the microstructure characteristics of the GCTJ of 304&P91, especially the interdiffusion zone (IDZ), have been systematically investigated. The microstructure at the interface transitions from austenite-ferrite (A+F) to austenite-martensite-ferrite (A+M+F), and finally to martensite-ferrite (M+F) due to diffusion. Additionally, the diffusion width between 304 and P91 increases with the volume fraction of P91. This unique design also ensures a gradual transition in both hardness and thermal expansion coefficient from 304 to P91, thereby enabling a smooth gradient in functional properties. Overall, this study proposes a novel approach for fabricating GCTJs and contributes to advancing design concepts in the field of dissimilar metal joining.

Additively manufacturing (AM)↗

Thermomechanical analysis and modeling of a high-temperature light trapping planar cavity solar receiver

The development of durable particle-based high-temperature solar receivers is critical for advancing concentrating solar-thermal (CST) technologies to enable high-efficiency power generation and industrial process heat. Here, this study presents a computational framework to evaluate the thermomechanical performance of a proposed enclosed light-trapping planar cavity receiver designed for particle-based thermal energy systems. The receiver incorporates absorptive cavities and fluidized particle-bed channels to enhance heat capture and reduce thermal losses. Finite element analysis (FEA) is employed to assess stress, strain, and creep-fatigue behavior under concentrated solar flux using realistic thermal boundary conditions derived from coupled system models and experimental assembly parameters. The analysis investigates the influence of particle-to-wall heat transfer coefficients (HTC) ranging from 800 to 1800 W/m 2 .K on the thermomechanical response of six candidate high-temperature alloys: Alloy 740H, Alloy 282, Alloy 617, 316H, Alloy 230, and 800H. Results show that increasing HTC reduces thermal gradients, leading to lower stresses and strains and extended minimum predicted creep life. While all materials satisfy fatigue life requirements under the investigated conditions, significant differences in creep resistance are observed. Alloy 740H consistently exhibits the longest minimum predicted creep life and the most favorable durability margins, followed by Alloy 282, with the remaining materials showing reduced creep resistance under identical loading. The reported creep lives are conservative lower-bound estimates intended for comparative material evaluation. This framework highlights the critical roles of material selection and geometry optimization in improving mechanical durability and reliability of solar-thermal receivers, forming a foundation for future experimental validation and design optimization.

14 SOLAR ENERGY↗

Turning a pyroelectric infrared motion sensor into a high-accuracy presence detector by using a narrow semi-transparent chopper

Pyroelectric infrared (PIR) sensors are the most popular presence detectors. However, their applications are limited to motion detection only, as the pyroelectric sensing element is only sensitive to radiated heat power fluctuation. Here, this letter reports a chopped PIR (C-PIR) presence sensor, capable of detecting both the stationary and the moving occupants at high accuracy by introducing a narrow semi-transparent optical chopper to shutter the incident radiation periodically. More importantly, even though the narrow chopper can efficiently shutter the unidirectional infrared (IR) power radiated by stationary occupants, it has near-zero influence on the omnidirectional power radiated from the ambient environment. Therefore, the voltage difference generated by the C-PIR sensor between occupied and unoccupied scenarios can be directly used for identifying the presence of stationary occupants with high-accuracy. For moving occupant detection, the optical chopper can barely affect the variation of the radiated power and thus lead to the same high-accuracy detection as the on-board PIR. An approximate theoretical model is created to illustrate the physics and dynamics of IR radiation shuttered by the optical chopper. Parametric studies, experimental characterization, and empirical analysis were performed to obtain the optimal setting of the chopper in terms of its thickness and chopping frequency. Experimental results reveal that the C-PIR sensor can detect stationary occupants for up to 4.0 m and moving occupants for up to 8.0 m, while maintaining an accuracy of 100%. The field of view of the C-PIR was found at 110° horizontally and 90° vertically, consistent with that of the on-board PIR sensor.

47 OTHER INSTRUMENTATION↗

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Factors and Considerations for Modeling Loss of a GaN-based Inverter

The article investigates the impacts of four often-neglected factors on the loss model of a GaN-based full-bridge inverter: parasitic capacitance of the devices, dynamics of junction temperature (Tj) under time-varying power dissipation (Ploss), case temperature estimation, and detailed considerations of the passive components. Procedures to calculate the converter loss considering the above factors are proposed and implemented. A 4.5-kW hard-switching inverter prototype using gallium nitrite (GaN) high-electron-mobility transistors is used to experimentally demonstrate the impact of each factor on the converter loss model. Furthermore, it is found that the accuracy of a converter loss model is mainly affected by the passive components at the light load condition, whereas the thermal and loss models of the active components become the major factors as the output power increases. The results show that after considering the above factors, the converter loss discrepancy between calculation and measurement can be reduced from 30.6 W (28%) to 2.5 W (less than 3%) at heavy load (Po = 4.5 kW), while at the light load condition (Po = 500 W), it is reduced from 3.9 W (28%) to 2.6 W (16%). Furthermore, the difference between simulated and measured case temperature of the GaN devices is within 6 °C.

42 ENGINEERING↗

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↗