Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ITER design”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Report of METL Mutual Inductance Level Sensor Development for Use in Liquid Metals – FY2023

A robust mutual inductance level sensor (MILS) for use in high temperature liquid metals has been developed at Argonne National Laboratory (ANL). The current design utilizes mineral insulated cables with a 300-series stainless steel sheath, magnesium oxide insulation, and a single copper conductor. Two coils are wrapped on a common core made of a 300-series stainless steel tube. One coil is energized with an alternating current (AC) source, and this electromagnetically couples with the second coil to generate an induced voltage. The voltage is measured to determine the mutual induction between the coils, and the mutual inductance measurement can be used to determine the level of a nearby liquid metal volume. The MILS can be located in an isolating thimble that is fully sealed to a vessel containing liquid metal for ease of maintenance and replacement. When the primary coil is energized with the AC source, the coil not only electromagnetic (EM) couples with the second coil, but also the surrounding liquid metal. The EM coupling with the surrounding liquid metal reduces the EM coupling with the second coil, producing an inverse-linear relationship between liquid metal level and secondary coil voltage. This has all been demonstrated at the Mechanisms Engineering Test Loop (METL) liquid sodium facility at ANL. The most current iteration of the MILS system is the MILS-MK-II. This sensor system has been commissioned in a non-sodium test stand where calibrations were performed using a sodium analog. The calibrations proved to be highly linear and repeatable. The MILS-MK-II has been installed in the METL expansion tank where high temperature sodium tests have been performed. The MILS-MK-II has been calibrated against a known standard at temperature of 300°C, and the calibrations have proved to be highly linear and repeatable. The calibrated MILS-MK-II has been in operation in the METL facility for several thousands of hours at temperatures around 300°C. The experimental data has been used to develop and validate electromagnetic finite element models in COMSOL Multiphysics, and now these models can be used to advance the development of the sensor system. Next steps will include temperature compensation to allow for operation at various temperatures. Additionally, efforts to incorporate a self-calibration methodology are underway.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Design, development and commissioning of a multi-alkali semiconductor photocathode deposition system for the IUAC Delhi light source photoinjector

A fourth-generation light source, called Delhi Light Source (DLS) based on photocathode-based RF gun has been commissioned at Inter-University Accelerator Centre, New Delhi. Presently, the electron beam is being generated from copper photocathode and the beam is being used for scheduled experiments. Soon, the semiconductor photocathode will be used to produce higher beam current. Here, to develop the semiconductor photocathode, a dedicated photocathode deposition facility was developed in collaboration with Brookhaven national Laboratory (BNL) and has been successfully commissioned and becomes operational at IUAC. This deposition facility is an integrated system with the electron gun and is a unique system as it is capable of producing, preserving (without residual gas poisoning) and in-vacuum transfer of the deposited photocathodes from the deposition chamber up to the RF electron gun. The system is designed to operate under ultra-high vacuum (UHV) and is equipped with load-lock chambers, substrate heating assembly, thickness monitoring via a quartz crystal microbalance (QCM), and an in-situ setup for quantum efficiency (QE) measurements. After testing of all the subsystems and a detailed calibration, the first deposition of a cesium telluride (Cs 2 Te) photocathode was successfully performed on a copper (Cu) substrate. This successful commissioning and initial deposition mark a significant step toward the indigenous photocathode development and lays the groundwork for further research into advanced photo emissive materials at IUAC. This paper will discuss the salient features, installation, commissioning, first semiconductor photocathode deposition and its results.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

From natural language to control signals: a conceptual framework for semantic channel finding in complex experimental infrastructure

Modern experimental platforms such as particle accelerators, fusion devices, telescopes, and industrial process control systems expose tens to hundreds of thousands of control and diagnostic channels, accumulated over decades of hardware evolution. Operators and AI systems alike depend on informal expert knowledge, inconsistent naming conventions, and scattered documentation to locate the signals required for monitoring, troubleshooting, and automated control, creating a persistent bottleneck for reliability, scalability, and emerging language-model-driven interfaces. We formalize semantic channel finding, the task of mapping natural-language intent to concrete control-system signals, as a general problem in complex experimental infrastructure, and introduce a four-paradigm conceptual framework to guide architecture selection based on facility-specific data regimes. The paradigms span (i) direct in-context lookup over small, curated channel dictionaries, (ii) constrained hierarchical navigation through structured trees, (iii) interactive agent exploration using iterative reasoning and tool-based database queries, and (iv) ontology-grounded semantic search that decouples channel meaning from facility-specific naming conventions. We demonstrate the practical feasibility of each paradigm through proof-of-concept implementations at four operational facilities spanning two orders of magnitude in scale: from compact free-electron lasers to large synchrotron light sources, operating under diverse control-system architectures ranging from clean hierarchical naming schemes to legacy environments with decades of heterogeneous conventions. Where evaluated against expert-curated operational queries, these instantiations achieve 90%–97% accuracy, validating the framework’s applicability across real-world deployment scenarios. To accelerate adoption across the broader scientific and industrial control-system community, we release open-source, plug-and-play implementations of all three interactive paradigms-direct lookup, hierarchical navigation, and middle-layer exploration-within the Osprey framework, together with tools for channel database generation, interactive testing, and minimal-configuration deployment. This work establishes semantic channel finding as a foundational capability for human-centric and agentic AI interfaces at large-scale facilities, providing both a systematic framework for architecture design and practical resources to enable adoption without building custom infrastructure from scratch.

channel finding↗
DOE OSTI2023

Impact of Glass Irradiation on Laser-Induced Breakdown Spectroscopy Data Analysis

Increased absorption of optical materials arising from exposure to ionizing radiation must be accounted for to accurately analyze laser-induced breakdown spectroscopy (LIBS) data retrieved from high-radiation environments. We evaluate this effect on two examples that mimic the diagnostics placed within novel nuclear reactor designs. The analysis is performed on LIBS data measured with 1% Xe gas in an ambient He environment and 1% Eu in a molten LiCl-KCl matrix, along with the measured optical absorption from the gamma- and neutron-irradiated low-OH fused silica and sapphire glasses. Significant changes in the number of laser shots required to reach a 3σ detection level are observed for the Eu data, increasing by two orders of magnitude after exposure to a 1.7 × 1017 n/cm2 neutron fluence. For all cases examined, the spectral dependence of absorption results in the introduction of systematic errors. Moreover, if lines from different spectral regions are used to create Boltzmann plots, this attenuation leads to statistically significant changes in the temperatures calculated from the Xe II lines and Eu II lines, lowering them from 8000 ± 610 K to 6900 ± 810 K and from 15,800 ± 400 K to 7200 ± 800 K, respectively, for exposure to the 1.7 × 1017 n/cm2 fluence. The temperature range required for a 95% confidence interval for the calculated temperature is also broadened. In the case of measuring the Xe spectrum, these effects may be mitigated using only the longer-wavelength spectral region, where radiation attenuation is relatively small, or through analysis using the iterative Saha–Boltzmann method.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Experimental evaluation of the vapor box divertor concept with an open vapor box module in Magnum-PSI

A promising approach to handle the intense plasma heat flux in the divertor region of a tokamak is the Vapor Box Divertor (VBD). Here plasma-lithium interaction creates a dense lithium vapor cloud which interacts with the incoming plasma, effectively shielding the tungsten surface beneath, therefore preventing overheating and sputtering of the tungsten and increasing the component’s lifetime. Two key steps must be addressed to validate this concept: investigating plasma-Li interactions and the transport mechanism of the latter in the presence of plasma. To explore this, a Vapor Box Module (VBM) has been designed for use with the linear plasma generator Magnum-PSI. The VBM consists of a series of three cylindrical boxes, with Li being evaporated at a controlled temperature in the central box (CB). Divertor-like plasma enters the VBM from an upstream aperture, interacts with the Li vapor cloud and exits through a downstream aperture ultimately impacting a target equipped with a calorimetry system. Optical diagnostics Thomson scattering, Filtered fast camera Imaging and optical emission spectroscopy provided information on plasma parameters in terms of electron density n e , temperature T e and plasma composition before and after this interaction. A significant reduction in plasma power was observed upon the establishment of lithium vaporization determined via cooling water calorimetry and embedded thermocouples at the downstream target. This resulted in a drop of the target temperature from ∼800 °C to ∼350 °C (57%) at the highest applied power (11.1 MW m −2 ) and CB temperature (∼700 °C). Lithium condensation at the side boxes of the VBM in combination with strong plasma momentum transfer effectively prevented upstream migration of lithium, while enhancing Li transport toward the target. The achievement of the two main goals of reducing plasma power and confining the Li in the VBM, are consistent with earlier published preliminary SOLPS-ITER simulations. This study shows that the presence of lithium in a VBD-like configuration can result in a strong reduction of power to the target surface, while at the same time the lithium vapor is effectively confined by the VBM, aided by the incoming plasma, preventing its escape from the VBM geometry. This represents a valuable step toward validating the feasibility of the VBD configuration in future fusion reactors.

Magnum-PSI↗
DOE OSTI2023

PipeSight: A High-Performance Computing Platform for Pipeline Integrity Management

The Phase I feasibility study completed as part of this project has led to a number of innovative technologies being developed and has laid the foundation for a successful Phase II effort to commercialize a platform for managing the integrity of pipelines for the damage mechanisms of the new, hybrid-energy based economy. To ground the development efforts and direction of the project, an extensive market research and customer discovery effort was undertaken early in Phase I. Through this effort, a number of pipeline owners and operators were interviewed, and the following key findings were discovered about the pipeline industry: • Small pipeline operators do not have the central engineering groups necessary to perform their own independent analysis of inspection data, but instead rely on summarized tally sheets provided to them by inspection service providers. • The time it takes to go from an inspection to a completed engineering assessment, even for small segments of pipeline, can take anywhere from 30-120 days. During this delay, critical threats can (and have been known to) cause failures. • Uncertainty is often not accounted for in the assessment of pipeline integrity. The tally sheets provided by third-party service providers are almost always deterministic in nature, identifying threats that present a concern only to the current (not the future) integrity of the pipeline. • It is uncommon to apply the latest technologies to perform advanced assessments of damaged pipelines. There is a desire to use more advanced analysis capabilities to assess threats. Many pipeline operators indicated that they would often excavate a pipeline to perform an inspection and find that the damage was not as bad as they anticipated, thus using limited resources unnecessarily. Companies are not consistent in their use of inspection data to determine corrosion rates, and those that do only calculate deterministic corrosion rates. • The industry has prominently relied on time-based inspections but has recently started to transition to risk-based inspections. However, there appears to be no uniform guidance on how to do so while properly accounting for all sources of uncertainty. • Companies are not storing inspection data in a manner that allows for the ready determination of temporal trends. • Predictive maintenance principles and practices are beginning to be used by early adopters • Some pipelines are being re-purposed to transport different process fluids than they were designed for, e.g., H 2 and CO 2 rich process streams to serve the new hybrid-energy based economy, which are presenting new integrity concerns for the existing pipeline network that crisscrosses the United States. As a result of these discoveries, we were able to target the development efforts in Phase I to best serve the needs of the industry. In Phase I, we developed a way to correlate multiple large-scale scans of the pipeline to determine a probabilistic corrosion rate that accounts for all sources of error and uncertainty in the inspection process. This probabilistic corrosion rate can be used to predict the future thickness distribution of the pipe wall. We demonstrate how this analysis may be performed in an analytical fashion and has been implemented in such a manner that it can be readily distributed using GPU computing through integration of the Kokkos programming model. We also make a very novel extension of the analytical corrosion rate model to Bayesian Networks (an explainable AI technique) that can account for non-parametric distributions of corrosion rates. With the predictions made above for the probabilistic corrosion rate and corresponding future distribution of the pipe wall thickness, we can assess the integrity of the pipeline through the use of a probabilistic engineering assessment. We developed a novel screening data analysis approach that can rapidly identify ‘hotspots’ (local thin areas) where the integrity of the pipeline is a concern. Once more, we implemented this screening approach in C++ to leverage GPU computing via the Kokkos programming model. After the critical hotspots are identified, we developed a program that can automatically generate an advanced finite element model of the damaged regions. Since the number of damaged regions that require advanced analysis can number in the thousands, we integrated an open-source container-native workflow engine for orchestrating parallel jobs on the cloud. Initially, these advanced numerical models were only designed to account for loading due to internal pressure. However, in a slight pivot from the initial Phase I proposal, we developed a complete pipe stress analysis program (called Simflex) which can simulate the complete pipeline and its response to thermal expansion, pressure, thermal bowing, weight, wind, earthquake, support displacement, support friction and external forces. This pipe stress analysis program was written generically, to handle any piping system, but contains the features needed to model long pipelines (i.e., it incorporates a model for soil mechanics and can account for the nonlinear boundary conditions necessary to simulate long underground pipelines). This pipe stress analysis program can simulate any segment of the pipeline (simple or complex) under any set of conditions and loads, to determine the supplemental loads (axial forces and bending moments) at the location of damage. This enables the most accurate state of stress to be accounted for in the pipeline, which can prove critical when evaluating the integrity of a damaged region. In the process of developing the technologies to perform the integrity assessment of the pipeline, we also extended one of the industry standard approaches for performing the assessment of local thin areas that extend more in the circumferential direction than the longitudinal direction of the pipeline. This approach was presented to the API 579-1/AS ME FFS-1 steering committee in November 2021 for consideration in the next edition of the industry standard for Fitness-For-Service (expected to be released in 2023). To help pipeline operators make decisions with the results on any integrity assessment, we developed a new approach to the life-cycle management of pipelines which uses a Bayesian Decision Network. The network is designed to help pipeline operators plan and prioritize inspection activities and ultimately make smarter, more cost-effective decisions. The Bayesian approach accounts for all sources of uncertainty and carries them through to the final optimal decisions, providing a probabilistic framework for optimizing inspection intervals. The proof-of-concept networks developed in the feasibility study are complete, verified, and are focused on a subset of the pipeline. To expand this novel approach to the scale necessary for an entire network of pipelines in Phase II, we will leverage the DOE-funded Bengi solver for industrial-scale decision making with Bayesian Networks [22]. Once implemented, we will be able to provide the pipeline industry with a much-needed tool for optimal inspection planning using truly explainable artificial intelligence (XAI). To handle all of these advanced capabilities into a cloud-based platform, the architecture of the Equity Engineering Cloud (EEC) was extended to include Argo Workflows, a framework capable of distributing and managing a massive number of jobs that consume their own resources, such that thousands of serial finite element simulations can be run in parallel. As part of this substantial undertaking, we also integrated Argo Continuous Delivery (CD) into the EEC, to aid with the rapid prototyping and iterations that will be imperative to the success of the PipeSight platform’s Agile development process in Phase II. As part of the pipe stress analysis program, we also developed a custom visualizer that leverages the DOE-funded VTK visualization library. We added custom contouring capabilities and a means for interacting visually with both the inputs and outputs of the pipe stress analysis program. We also developed routines for automating the post-processing of the finite element simulations to determine if any failure criteria are met and to visualize the deformations, stresses and strains in ParaView using the exodus II file format (a subset of netCDF).

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Modeling DOE Standard Canister Configurations with Updated Surface Chemistry

Road-ready and final disposition packaging configurations for the advanced test reactor (ATR) fuel currently specifies storage within helium backfilled DOE sealed standard canisters. The aluminum cladding of the ATR fuel contains an oxyhydroxide layer of boehmite/bayerite that generates hydrogen when subjected to irradiation. Understanding the effect of this hydrogen buildup over time to important for long term storage considerations. Previous modeling efforts have built a coupled CFD-chemical model to simulate the temperature gas phase concentrations within the DOE sealed standard canisters. These models have coupled the temperature conditions to both the gas phase radiolysis chemistry and the surface chemistry associated with the oxyhydroxide layer. Continued experimental work has identified trends for aluminum surrogate samples with oxyhydroxide layers for tests done at higher dose rates. At low initial doses a fast generation rate of hydrogen occurs which starts to roll over to a lower generation rate as the total dose applied in increased. A previous study created a small-scale chemical model was built to replicate a mini-canister surrogate system at SRNL as well as for the smaller capsule tests performed at INL. The previous iteration of the model utilized step function for its G-values for the hydrogen generation over a 200-year period. This model makes an update to the surface chemistry to account for theorized surface chemistry reactions allowing for oxygen to remain bounded to the oxyhydroxide layer. As additional experiments continue, the kinetic fits may be adjusted to adapt to new data. Three canister configurations are modeled – the base 18-inch, 15-foot DOR standard canister with 30 ATR fuel elements, an 18-inch, 10-foot DOE standard canister with 32 ATR fuel elements and a 24-inch, 10-foot DOE standard canister with 40 ATR fuel elements. In previous reports, the primary sensitivity of the canister conditions was identified as the decay heat of the fuel and the dried condition of the fuel. Only these parameters are studied in the report. For the nominal case with dried fuel, the hydrogen generation is only about 2%, so it is even less than the flammability limit if it were exposed to oxygen. For the densely packed 18-inch case the total hydrogen concentration is only around 4% with the nominal decay heat. For the densely packed 24-inch case the total hydrogen concentration is only 2.5%, roughly 20% higher than the original packing design, and still under flammability conditions if exposed to oxygen.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

DOE Bioenergy Technologies Office (BETO) 2023 Project Peer Review: WBS 2.4.1.100 Bench Scale Research & Development

Bench Scale Integration develops and optimizes fermentation processes to produce bio-based fuels and chemicals for commercial scale-up. The project uses fermentation science to achieve high titers and production rates by, for example, manipulating how the microorganisms are fed biomass sugars and nutrients, modifying fermentation conditions (pH, temperature, aeration) or developing online control strategies for better fermentation operations and high titer, rates, and yield (TRY). For this period of performance, we continued our development of a commercial-ready 2,3-butanediol (BDO) fermentation from biomass sugars utilizing NREL's proprietary Zymomonas mobilis microorganism. The engineered Z. mobilis can use all the main sugars in corn stover biomass, which are glucose, xylose, and arabinose. BDO is a versatile, low-carbon chemical which can be catalytically upgraded to a variety of hydrocarbon fuels and chemicals. The project had three goals during this review period; evaluate the technical feasibility of using whole slurry pretreated corn stover to achieve the techno-economic analysis (TEA) performance goals, continue optimizing a liquor-based fed-batch fermentation process for high titer, and develop strategies to enable scale-up. After evaluating different iterations of a whole slurry fermentation that did not meet the TEA goals, a Go/No-Go decision was made to pivot to liquor-only with new TEA performance targets, the main one being 140 g/L titer. We successfully met this goal, producing 141 g/L BDO at 1 g/L-hr productivity and 84% process yield. This titer and productivity attracted industrial interest to scale the fermentation resulting in a Technology Commercialization Fund project award in 2022. We used an NREL developed near-infra red (NIR) spectroscopy method for rapid analysis which allowed for changes to aeration levels and sugar feeding during the fermentation to maximize BDO production. The NIR analysis can be done using a hand-held spectrometer, essentially taking the analysis on to the plant floor during scale-up and preliminary work shows the feasibility of using an online probe for continuous monitoring and control. The other scale-up tool is mapping oxygen transfer coefficient (kLa) and oxygen transfer rate (OTR) in various vessels to find conditions that match the optimized 500 mL vessels. Showing a correlation to the mapping work, which is done with water and a dissolved oxygen probe, can reduce the risk of failed fermentations during scale-up. Our end-of-project goal is to meet the design target BDO titer (140-150 g/L) from DMR corn stover liquor at 1000L or larger scale to demonstrate BDO process design case and transfer the technology to industrial fermentation stakeholders for commercialization.

2 3-butanediol↗
DOE OSTI2021

Comprehensive deletion landscape of CRISPR-Cas9 identifies minimal RNA-guided DNA-binding modules

Proteins evolve through the modular rearrangement of elements known as domains. Extant, multidomain proteins are hypothesized to be the result of domain accretion, but there has been limited experimental validation of this idea. Here, we introduce a technique for genetic minimization by iterative size-exclusion and recombination (MISER) for comprehensively making all possible deletions of a protein. Using MISER, we generate a deletion landscape for the CRISPR protein Cas9. We find that the catalytically-dead Streptococcus pyogenes Cas9 can tolerate large single deletions in the REC2, REC3, HNH, and RuvC domains, while still functioning in vitro and in vivo, and that these deletions can be stacked together to engineer minimal, DNA-binding effector proteins. In total, our results demonstrate that extant proteins retain significant modularity from the accretion process and, as genetic size is a major limitation for viral delivery systems, establish a general technique to improve genome editing and gene therapy-based therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A thermodynamic model of integrated liquid-to-liquid thermoelectric heat pump systems

Thermoelectric (TE) heat pumps (TEHPs) are advantageous for heating and cooling in various applications because of their modularity and simple design. A TEHP system includes the TE modules with p- and n-type materials bonded to substrates, plus heat exchangers, thermal interfaces to the heat exchangers, and heat transfer fluids. Although modeling an individual TE module has been extensively studied, limited studies have reported performance at the larger system-level. Furthermore, no prior study has addressed the impact of temperature-dependent TE material properties (e.g., electric resistivity, thermal conductivity, and Seebeck coefficient) on overall heat-pump-system-level performance. Here, this work presents a mathematical model for TEHP system performance based on Goldsmid's approach for TE material performance, “effective” TE material properties, Gnielinski's correlation for convective heat transfer, and thermal balance theory for a heat exchange network. This combined approach provides an accurate model of the liquid-to-liquid TEHP system. Three different approaches—one empirical, one based on the manufacturer's specifications, and one drawn from the literature—were then used to determine values for TE material properties. The first two methods treated properties as constants, while the last approach treated properties as surface-temperature-based functions. Finally, experimental TEHP data was used to validate the models, all with relative absolute deviations of approximately 10% when predicting heating capacity and 10%–25% when forecasting cooling capacity up to a 30 K surface temperature lift. The results demonstrated that, at the TEHP system level, the TE material properties could be treated as constants, avoiding solver iterations and reducing the performance uncertainty by up to 95%.

42 ENGINEERING↗
DOE OSTI2025

Neural network-based classification and regression of magnetohydrodynamic modes in tokamaks

We present a machine learning-based magnetohydrodynamic (MHD) classifier and regressor that utilizes real or complex-valued 3D magnetic sensor array data to determine neoclassical tearing mode (NTM) onset times in tokamaks with millisecond accuracy. The input dataset consists of poloidal profiles of complex Fourier amplitudes with an n = 1 toroidal mode number from 144 human-labeled ITER Baseline Scenario discharges in the DIII-D tokamak, spanning both tearing-dominated and sawtooth-dominated regimes. Since m, n = 2,1 NTMs frequently emerge alongside sawteeth at the same frequency in this scenario, the focus is on isolating the m = 1 and m = 2 components of the n = 1 MHD mode near the tearing onset. To improve model regularization and prediction stability, singular value decomposition was applied to balance the sawtooth and tearing datasets. The enriched datasets facilitated training neural networks that learn the key distinguishing features of sawtooth and tearing modes in the poloidal profiles of their magnetic amplitude and phase. When the modes occur independently, the networks achieve perfect classification due to the modes’ distinct characteristics and low measurement noise. In the more experimentally relevant case where both modes coexist, the networks maintain exceptional performance across key metrics. Tests on synthetic data with known ground truth demonstrate the superior accuracy of the neural network trained on complex-valued input compared to models using real amplitude, phase, or pseudo-complex data, achieving both a mean time delay and standard deviation below 1 ms. Notably, standard linear regression methods fitting the dominant singular modes to the data closely match the neural network’s performance. Applying these methods across a broad range of H-mode scenarios will enable future studies to systematically identify dominant NTM triggers as scenario-specific variables, paving the way for more effective tearing mode avoidance strategies in future fusion reactor designs.

machine learning↗
DOE OSTI2024

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

AC loss and contact resistance in REBCO CORC ® , Roebel, and stacked tape cables

Abstract Many high-temperature superconductor (HTS) applications require superconducting cables with high currents while operating in an alternating magnetic field. HTS cables should be composed of numerous superconducting tapes to achieve the required current capacity. Alternating current and magnetic fields cause AC losses in such cables and can provoke conductor instability. AC losses and contact resistances were measured of several cable designs based on commercially available REBCO tapes at the University of Twente. The AC loss was measured under identical conditions for eight REBCO conductors manufactured according to three types of cabling methods—CORC ® (Conductor on Round Core), Roebel, and stacked tape, including a full-size REBCO CICC (cable in conduit conductor). The measurements were done at T = 4.2 K without transport current in a sinusoidal AC magnetic field of 0.4 T amplitude and frequencies from 5 to 55 mHz. The AC loss was measured simultaneously by calibrated gas flow calorimeter utilizing the helium boil-off method and by the magnetization method using pick-up coils. Also, the AC loss of two CORC® conductors and a Roebel cable was measured at 77 K. Each conductor was measured with and without background field of 1 T. The measured AC coupling loss in the CORC ® and Roebel conductors is negligible at 4.2 K for the applied conditions while at 77 K coupling loss was observed for all conductors. The absence of coupling loss at 4.2 K can be explained by shielding of the conductor interior; this is confirmed with measurement and calculation of the penetration field of CORC ® and Roebel cables. The inter-tape contact resistance was measured for CORC ® and stacked tape samples at 4.2 and 77 K. It was demonstrated that a short heat treatment of CORC ® conductor with solder-coated tapes activates tape-to-tape soldering and decreases the contact resistance. The reduction of contact resistance by two orders in magnitude to tens of nΩm is comparable with the interstrand contact resistance in ITER Nb 3 Sn type conductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

First-principles determination of the phonon-point defect scattering and thermal transport due to fission products in ThO 2

This work presents the first-principles calculations of the lattice thermal conductivity degradation due to point defects in thorium dioxide using an iterative solution of the Pierels-Boltzmann transport equation. We have used the non perturbative Green’s function methodology to compute the phonon- point defect scattering rates that consider the local distortion around the point defect, including the mass difference changes, interatomic force constants, and structural relaxation near the point defects. The point defects considered in this work include the vacancy of thorium ($V$ $TH$ ) and oxygen ($V$ $O$ ), substitutions of helium ($He$ $Th$ ), krypton ($Kr$ $Th$ ), zirconium ($Zr$ $Th$ ), iodine ($I$ $Th$ ) and xenon ($Xe$ $Th$ ) in the thorium site, and the three different configurations of the Schottky defects. The results of the phonon-defect scattering rate reveal that, among all the considered intrinsic defects, the thorium vacancy and helium substitution in the thorium site scatter the phonon most due to the substantial changes in the force constant and structural distortions. The scattering of phonons due to the substitutional defects unveils that the zirconium atom scatters phonons the least, followed by xenon, iodine, krypton, and helium. This is contrary to the intuition that the scattering strength follows $He$ $Th$ > $Kr$ $Th$ > $Zr$ $Th$ > $I$ $Th$ > $Xe$ $Th$ based on the mass difference. This striking difference in the zirconium phonon scattering is due to the local chemical environment changes. Zirconium is an electropositive element with valency similar to thorium and, therefore, can bond with the oxygen atoms, thus creating less force constant variance compared to iodine, an electronegative element noble gas helium, xenon, and krypton. In conclusion, these results can serve as the benchmark for the analytical models and help the engineering-scale modeling effort for nuclear design.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Structural Analysis of the US ITER ECH Transmission Line System

The electron cyclotron heating (ECH) and current drive system is one of the main plasma heating systems for ITER. It uses high-power microwave beams with the power deposition location steerable across the plasma cross section. Microwave power is conveyed via transmission lines (TLs) that run from the gyrotrons in the radio frequency building through the assembly hall and tokamak building to the ECH launchers within the tokamak vacuum vessel. The ECH system includes a vast array of interconnected TL waveguides, in-line components, and support structures.Finite element (FE) modeling provides an essential means of simulating the system, applying loads and determining deflections, rotations, forces, moments, and stresses in order to evaluate various structural and microwave transmission performance metrics. A representative FE model of the overall ECH TL system is developed in ANSYS®. This top-level model defines the centerline of the waveguide system. Waveguide segments are represented by line elements (beams and pipes) with equivalent section properties, and support structures are represented by boundary conditions. A systematic approach is used to model each ECH component with lumped masses and structurally equivalent stiffness matrices or ANSYS superelements.The top-level TL FE model is used to evaluate the various loads (thermal, vacuum, seismic, etc.) and operating scenarios. The top-level model directly calculates stresses in the straight aluminum waveguide segments. The model provides the forces and moments acting on the in-line components for detailed submodel assessments. Displacement results from the top-level analysis feed into a separate microwave performance model to help determine operational efficiency. All TL performance and thermal-structural requirements are met, as specified by the applicable codes and standards, and successfully documented in numerous technical reports and demonstrated at the final design review.

Shanmugasundaram, Aravind↗
DOE OSTI2022

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↗
DOE OSTI2026

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗
DOE OSTI2020

Innovations in Building Energy Modeling: Research and Development Opportunities for Emerging Technologies

Building energy modeling (BEM) is a multipurpose tool for building energy efficiency (EE). The U.S. Department of Energy Building Technologies Office (BTO) seeks to expand the use and effectiveness of BEM in the design and operation of commercial and residential buildings with the goal of achieving lasting reductions in total and peak energy use. This report identifies gaps and outlines recommended initiatives to achieve this goal, based on a combination of technical analysis and stakeholder input. In addition to BTO, this report can benefit BEM professionals (architects, mechanical engineers, energy consultants, building auditors, equipment manufacturers, and BEM software vendors) and BEM clients (building owners and operators, EE program administrators, EE service providers, policymakers, and policy and code jurisdictions such as states and cities). This report was developed in two phases. In the first, BTO worked with a team from Navigant Consulting (now Guidehouse) to characterize objectives, opportunities, and current activities; identify gaps and barriers; and define initiatives. To collect input, Navigant conducted telephone interviews and workshops with industry experts. The initial phase produced a draft report, which was released for public review in 2016 and yielded over 400 comments. Based on these comments, BTO compiled a second draft report that addressed many of those comments while acknowledging changes that had occurred both at BTO and in the industry. Unlike the first draft report, the second focused much more heavily on BTO’s own role, portfolio, and activities. BTO is a direct player in the BEM field - it funds the development of several significant software packages that are embedded in commercial products - and transparency about its goals and future plans is requisite. BTO recognizes that a great number of other public and private organizations contribute to the BEM enterprise. With the second draft report, BTO did not attempt to produce a blueprint for the industry as a whole, but rather a working document BTO can use to iteratively solicit stakeholder input and synthesize it into a program. BTO released the second draft report for public review in 2019. The second round of review generated 83 pages of feedback and comments - almost exactly the length of the draft report itself - a significant portion of which was collected and synthesized by IBPSA-USA Advocacy Committee. This final report incorporates this feedback. This report does not address the use of BEM in support of building-based grid services, a recent BTO initiative called Grid-interactive Efficient Buildings (GEB). In 2019, BTO published a report that specifically addresses the role of BEM - and other “integration” technologies - in GEB.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗