Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “code validation”

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 343 records · Page 19

A methodology for domain overlapping coupling of thermal-hydraulic systems

Multi-scale coupling has increasingly drawn attention as a promising approach for modeling thermal systems. Traditional system codes provide fast and robust predictions at the plant scale, while high-fidelity computational fluid dynamics (CFD)-based tools resolve localized flow and heat transfer phenomena with greater accuracy. By combining these complementary strengths, co-simulations enable multi-scale analysis that would otherwise be computationally prohibitive for a standalone CFD code. Here, this work introduces a robust and problem-agnostic domain overlapping (DO) coupling between the system thermal-hydraulic (STH) code System Analysis Module (SAM) and the coarse-mesh CFD code Pronghorn. Both applications belong to the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite, a code suite in active development at the Idaho National Laboratory (INL), tailored for multi-physics analysis of advanced reactors. Unlike previous approaches, BlueCRAB supports an agnostic interface between codes based on different fidelity, while its coupling formulation can address arbitrary flow geometries with multiple inlets and outlets in coupled components. The implemented method leads to consistent pressure drops, enthalpies, and scalar concentrations between coupled SAM and Pronghorn simulations. The methodology is demonstrated through two verification tests, which ensure the numerical consistency and conservation across the codes, and through one validation test against experimental data. The proposed problems explore different physical aspects inherent to thermal systems, with particular attention given to nuclear reactor analysis. These include buoyancy-driven flows, complex flow patterns, and setups with multiple inlets and outlets, representing challenges in advanced reactor applications.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

End-To-End Decentralized Transmission Line Protection in IBR-Dominated Weak Grids Using Interpretable Data-Driven Methods

Traditional transmission line protection relies on predictable synchronous-based fault signatures, which frequently fail under the non-standard, current-limited fault characteristics of Inverter-Based Resources (IBRs). This study investigates how to achieve secure, communication-free fault isolation in IBR-dominated weak grids without relying on opaque, computationally heavy "black-box" machine learning algorithms. To address this, we propose a novel, standalone, and inherently interpretable data-driven protection framework. Unlike centralized methods requiring multi-terminal communication, this decentralized approach relies solely on local measurements using a hierarchical linear-kernel Support Vector Machine (SVM). The methodology decomposes the protection task into four sequential stages that mimic traditional protection elements: fault detection and fault direction identification, fault type classification, zone classification, and location estimation. This multi-stage architecture allows for specialized feature engineering at each stage, combining high computational efficiency with logic traceability. The framework's end-to-end performance was validated via C-code and PSCAD/EMTDC co-simulation, utilizing a real-world utility network and an OEM black-box IBR model. The proposed relay achieves 97.2% overall accuracy and provides a reliable trip decision within a 2.5-cycle window. The results confirm 100% accuracy in fundamental fault detection, reliable zone selectivity across low to moderate fault resistances, and robust security against non-fault transients, proving its immediate viability for integration into commercial numerical relays.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A multi-dimensional dataset for two-phase instability in low pressure natural circulation based on direct transient local measurement

Thermal hydraulic instability can challenge the control and the safety of nuclear reactors. This experimental research captures and investigates the instability in a low-pressure (1–10 atmospheres) natural circulation loop with an adiabatic chimney installed above the heated section. Direct transient local measurement of two-phase flow is performed at multiple axial locations along the flow channel with thermocouples and conductivity probes traversed by customized high-precision mechanisms. Based on the periodicity of oscillations, a post-processing method is proposed to align the data captured non-simultaneously to reconstruct the instability behavior in one representative period. The averaging operation over multiple oscillation cycles is performed to reduce the random errors. Cyclic phenomena are clearly demonstrated by the post-processed data, and rapid void fraction change induced by the flashing is observed to be significant in the driving mechanism for low pressure oscillations. Parametric effects of inlet subcooling, heat flux, inlet flow restriction and system pressure are studied within the experiment range, and qualitative results are explained and compared with existing literature. The ratio between the oscillation period and the travelling time of the fluid particle through the test section is compared with the theoretical value for the density wave oscillation. The new dataset contains directly measured transient multi-dimensional two-phase data along the flow channel for natural circulation under density wave oscillations. The generated dataset is valuable for the validation of system analysis codes and CFD codes in the simulation of natural circulation phenomena and dynamic two-phase flow.

42 ENGINEERING↗

Microstructure-based prediction of thermal aging strength reduction factors for grade 91 ferritic-martensitic steel

This study aimed to develop a microstructure-based mechanistic approach to address the long-term thermal aging effect on yield stress and ultimate tensile strength and to provide a physical basis for developing thermal aging factors for G91 for a design life of 60 years. Several heats of G91 steel were examined. Controlled aging experiments were conducted on two heats at temperatures of 550, 600, and 650 °C for times up to ~64,000 h. Specimens were also taken from archived creep-tested specimens of G91 and from a tube removed from the Kingston coal-fired power plant to obtain data that cover a wide range of temperatures and for aging times up to 155,000 h. Thermal aging caused significant subgrain recovery, coarsening of M 23 C 6 carbides at subgrain boundaries and MX carbonitrides within subgrains. The intermetallic Laves phase forms during aging and grows rapidly. The growth rate of the subgrain width and MX mean size during thermal aging were described by kinetic models. Thermal aging results in the reductions in the yield stress and the ultimate tensile strength of G91. The effects of thermal aging on the reductions in yield stress and ultimate tensile strength were well described by the microstructure-strength model that considers three superimposed strengthening mechanisms, namely, sub-boundary strengthening, MX precipitation hardening, and Mo solid solution strengthening. Furthermore, the model was independently validated by the ASME Code values and is being used by the ASME to develop the aging-induced strength reduction factors for G91 steel for a design life of 60 years.

36 MATERIALS SCIENCE↗

Scaling methodologies and similarity analysis for thermal hydraulics test facility development for water-cooled small modular reactor

Small modular reactors (SMRs) represent a promising option for providing clean and sustainable energy due to their potential for enhanced safety, reduced capital costs, and increased siting flexibility. However, new reactor systems require the development and operation of representative scaled-down test facilities to support the verification and validation of system computer codes and models. Here this study reviews the research on scaling methodologies and similarity principles pivotal in developing non-nuclear integral effects test and separate effects test facilities for water-cooled SMRs. The study focuses on a review of the scaling methods, similarity approaches, and possible challenges posed by the unique and compact design features of integral-pressurized water reactor-type SMRs, and their representative test facilities. This study also reviews previous research related to scaling and similarity methodologies and provides insights into design considerations for achieving prototypic conditions in test facilities. The findings and recommendations emphasize the broader impact of appropriate scaling and similarity principles to ensure meaningful and transferable results from non-nuclear test facilities to accelerate the safe and efficient deployment of next-generation water-cooled SMRs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Resolving ECRH deposition broadening due to edge turbulence in DIII-D

Microwave heat pulse propagation experiments have demonstrated a correlation between millimeter-scale turbulence and deposition profile broadening of electron cyclotron (EC) waves on the DIII-D tokamak. In a set of discharges in DIII-D, a variation in edge density fluctuations on the mm-scale is associated with 40%–150% broader deposition profiles, expressed in terms of normalized minor radius, as compared with equilibrium ray tracing. The 1D power profile is determined from transport analysis of the electron temperature response to EC power modulation using perturbative analysis with a square wave power modulation at 20–70 Hz, producing a series of Fourier harmonics that are fit collectively to resolve transport. Fitting an integrated heat flux expressed in the Fourier basis of the modulation to diffusive, convective, and coupled transport terms in a linear model can resolve the broadened EC deposition width from the power perturbation to resolve a broadening in each case. The best fit degree of beam broadening observed scales approximately linearly with the Doppler backscattering measured fluctuation level in the steep gradient region. Quantifying the effect of edge fluctuation broadening on EC current drive power needs of future devices will require 3D full-wave codes that can be validated on the current generation of machines. These DIII-D experiments provide a quantitative measure of fluctuation effects and a dataset to benchmark full-wave simulations that can model and eventually predict nonlinear effects neglected by 1D equilibrium beam and ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydro-scaling of direct-drive cylindrical implosions at the OMEGA and the National Ignition Facility

Deceleration-phase Rayleigh–Taylor instability (RTI) growth during inertial confinement fusion capsule implosions significantly affects the performance as it mixes cold ablator material into the fuel. Precise measurements of such instability growth are essential for both validating the existing simulation codes and improving our predictive capability. RTI measurements on the inner surface of a spherical shell are limited and are often inferred indirectly at limited convergence. In contrast, cylindrical implosions allow for direct diagnostic access to the converging interface by imaging down the cylinder axis while retaining the effects of convergence. We have performed direct-drive cylindrical implosion experiments at both the OMEGA and the NIF laser facilities using scaled targets. RTI growth is demonstrated to be scale-invariant between the cylindrical targets at OMEGA and similar targets at the NIF, which are scaled up by a factor of three in the radial dimension. Single-mode (m = 0) instability growth factors of ~14 are measured at a convergence ratio (CR) ~ 2.5 with nearly identical mode growth at both scales. The measurements are in agreement with xRAGE radiation-hydrodynamics simulations. In addition, we have developed the Bayesian-inference-engine method to account for the variations in the target alignment, magnification, and the parallax effect in the measurement, allowing a more precise comparison between the experimental data and the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental investigation of beam-ion losses induced by magnetic perturbations using the light ion beam probe technique in the ASDEX Upgrade tokamak

The impact of externally applied magnetic perturbations (MPs) on fast-ion losses has been investigated by means of the light ion beam probe (LIBP) technique in the ASDEX Upgrade tokamak. The LIBP technique allows to experimentally infer the fast-ion orbit displacement induced by MPs via first-orbit losses using scintillator based fast-ion loss detector (FILD) measurements. The fast-ion orbit displacement against different applied MP spectra has been studied. These shots were conducted in ELM mitigated H-mode plasmas. A rigid rotation of the MP coils was applied with a frequency of 1 Hz, with an n = 2 configuration and changing the differential phase between the upper and lower set of coils (ΔΦ ul ) on a shot-to-shot basis. Beam sources Q7 (tangential) and Q8 (radial) were used to probe different fast-ion orbits with FILD1. The measured fast-ion orbit displacement ranges from 3 to 20 mm approximately, and no qualitative difference is observed between ions from beam sources Q7 and Q8. The minimum is found for a ΔΦ ul ~ 50°, which is shifted with respect to the minimum of the plasma boundary displacement, found at ΔΦ ul ~ 0°. A first attempt to validate the orbit following code ASCOT—including the plasma response calculated with the MARS-F code—against these experimental measurements is performed. Finally, while the dependence of the first-orbit fast-ion displacement with ΔΦ ul does not match the experimental measurements, these simulations do capture other features such as the order of magnitude of the orbit displacement and the importance of the toroidal spectrum of the applied perturbation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Determining the 9 $\mathrm{Be}$($n, γ$) 10 $\mathrm{Be}$ integral cross section at fission neutron energies

The 9 B neutron capture cross section has significant implications for Be materials in the nuclear industry as well as the α process in stellar nucleosynthesis. While the cross section is well constrained at thermal neutron energies, there is a lack of experimental data at higher neutron energies, and the evaluated nuclear data libraries can differ by up to two orders of magnitude. We calculate the 9 Be(n, γ) 10 Be integral cross section at fission neutron energies in an effort to resolve disagreements amongst the nuclear data libraries. Foil irradiation experiments were performed using the Flattop critical assembly at the National Criticality Experiments Research Center with either the highly enriched U or Pu cores, with target foil stacks placed at multiple locations to exploit different neutron energy profiles. Accelerator mass spectrometry was used to measure the 10 Be/ 9 Be ratio in irradiated Be foils, while all other activation products were quantified through gamma spectrometry. The experiments were simulated using the Monte Carlo N-Particle radiation transport code and combined with experimental results to determine the total neutron fluence, while the staysl-pnnl suite and fispact-ii code were used to validate the model and assess the systematic uncertainty. The new 9 Be(n, γ) 10 Be integral cross sections calculated in this work are 26.5 ± 2.2µb at 0.59 ± 0.07 MeV, 24 ± 3 µb at 0.98 ± 0.14 MeV, 21.7 ± 1.3 µb at 1.26 ± 0.11 MeV, 21.8 ± 1.4 µb at 1.32 ± 0.11 MeV, and 18.6 ± 1.1 µb at 1.46 ± 0.13 MeV. These results do not agree with integral cross sections from any of the nuclear data library evaluations. Discrepancies between the new integral cross sections reported here and the nuclear data libraries suggest a more complex cross-section structure in the MeV range which allows for more resonance contributions, and more work is needed to further constrain the evaluated cross sections.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

An iron evaluation story: from TALYS model parameters to validation on the ASPIS benchmark with the Monte-Carlo code TRIPOLI-4

Iron is a widespread structural material in nuclear applications. However some benchmark studies show that for many libraries, including TENDL-2019, producing good iron evaluations is still an open challenge. This work aims at improving TENDL-2019 {sup 56}Fe and {sup 54}Fe evaluations through feedbacks from experimental data and Monte Carlo results of the iron benchmark ASPIS. Part of the T6 code infrastructure was used to produce the evaluations. In the new evaluations, particular attention to resonance parameters and formalism was paid to address the importance of elastic scattering. The total inelastic scattering and inelastic continuum cross sections were shaped by modifying the input parameters of TALYS, one of the T6 codes. Additionally, it was necessary to normalise some reaction channels to those of JENDL-4.0. With the new files ASPIS results, obtained with the Monte Carlo code TRIPOLI-4, improved noticeably compared to TENDL 2019. The evaluations and parameter setting modifications are to be included into TENDL-2021. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Study of reference burnup steps optimization in fuel segment data file generation for NEXUS/ANC9 code system

For any two-step core design code system, the cross-section files are the primary factor to determine the accuracy of the system prediction. For a once-through cross-section system to cover all potential applicable conditions, the system may need to perform tens of thousands state points of lattice calculations. These calculations take a significant amount of CPU times and computing power, especially when a fine or ultra-fine energy group library is used. As computer power increases, so does the calculational complexity. Therefore, it is important to make the lattice calculations effective and efficient. Based on the Westinghouse core design code system NEXUS/ANC9, this work focuses on the reference burnup steps used by the lattice code when performing the calculations. Following the fundamental cross-section methodology and analyzing the contribution of each individual terms, this study provides an applicable solution of the optimized reference burnup steps. The test results show that, with the existing cross-section methodology, it is possible to significantly reduce the lattice calculation cases and cross-section file generation time without sacrificing the accuracy of system prediction. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Critical Unresolved Region Integral Experiment Execution

Integral Experiments are a key aspect of validating nuclear data behavior and simulation capabilities. Further, integral experiments are the earliest form of validation of equations and codes. They have been used extensively since the 1940s. Today, simulation and predictive capabilities have improved greatly from 75 years ago, but nonetheless, integral experiments are still needed. As simulation capabilities improve, their uncertainties get much smaller. The current focus of many new experiments is understanding the intermediate energy range. The intermediate energy range does not have a wide application space and it can be challenging to design experiments that are sensitive to neutron energies in this region. The Zeus experiments evaluated highly enriched uranium (HEU) in the intermediate energy region. Building on the success and knowledge gained from the Zeus experiments, the Critical Unresolved Region Integral Experiment (CURIE) experiments were designed to evaluate HEU in the narrower unresolved resonance region (URR). These experiments were executed during June and July 2020 at the National Criticality Experiments Research Center (NCERC). The National Criticality Experiments Research Center (NCERC), operated by Los Alamos National Laboratory, is the only general purpose critical experiments facility in the United States of America. NCERC regularly designs and executes critical experiments and other measurements useful to a wide variety of fields including nuclear criticality safety, commercial nuclear energy, nonproliferation, and nuclear data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spatio-temporal focal spot characterization and modeling of the NIF ARC kilojoule picosecond laser

The advanced radiographic capability (ARC) laser system, part of theNational Ignition Facility (NIF) at Lawrence Livermore National Laboratory, is a short-pulse laser capability integrated into the NIF. The ARC is designed to provide adjustable pulse lengths of ~1–38 ps in four independent beamlets, each with energies up to 1 kJ (depending on pulse duration). A detailed model of the ARC lasers has been developed that predicts the time- and space-resolved focal spots on target for each shot. Measurements made to characterize static and dynamic wavefront characteristics of the ARC are important inputs to the code. Modeling has been validated with measurements of the time-integrated focal spot at the target chamber center (TCC) at low power, and the space-integrated pulse duration at high power, using currently available diagnostics. In this work, these simulations indicate that each of the four ARC beamlets achieves a peak intensity on target of up to a few 10 18 W=cm 2 .

47 OTHER INSTRUMENTATION↗

Towards Auto-Generated Data Systems

After decades of progress, database management systems (DBMSs) are now the backbones of many data applications that we interact with on a daily basis. Yet, with the emergence of new data types and hardware, building and optimizing new data systems remain as difficult as the heyday of relational databases. In this paper, we summarize our work towards automating the building and optimization of data systems. Drawing from our own experience, we further argue that any automation technique must address three aspects: user specification, code generation, and result validation. We conclude by discussing a case study using videos data processing, along with opportunities for future research towards designing data systems that are automatically generated.

Computer Science↗

Transverse BBU Suppression with Feedback for Energy-Recovery Linacs

The Electron–Ion Collider (EIC) achieves its design luminosity by cooling the ion beam with a high-current electron beam generated in an energy-recovery linac (ERL). The baseline ERL lattice employs a BNL five-cell cavity that is frequency-scaled to 197 MHz, 591 MHz, and 1.773 GHz, raising concerns about multibunch beam-breakup (BBU) instabilities. Threshold currents for each frequency option are established with two independent BBU tracking codes, providing cross-validated operating margins that guide cavity selection. To further increase the current limit, we incorporate a digital transverse feedback (FB) system that targets the dominant higher-order modes (HOMs). Simulations show the FB raises the BBU threshold by roughly an order of magnitude: power-spectral-density analysis of the beam centroid at the linac exit confirms strong suppression of the aliased HOM peak, although some spectral growth appears at secondary frequencies. This combined study quantifies baseline BBU limits, demonstrates effective active mitigation, and charts a practical path toward robust, high-current operation of the EIC ERL.

Accelerator Physics↗

Collaborative Research: Improved Efficiency and Coupling of the Radiation Code in the ACME Earth System Model (Final Report)

The complexity of radiative transfer, its importance to the exchange of energy in the climate system, and its high computational cost establishes importance of an accurate and efficient radiative transfer parameterization for climate simulation. Previous work by the proposing team at AER led to the development of the radiation code, RRTMG, which has been widely accepted by the global modeling community as a fast and accurate advancement over the previous generation of radiation codes. It has been in use in the NCAR CESM for many years, and it has been implemented in the initial version the DOE E3SM model. However, its computational cost remains high relative to other components in part due to its complexity and to its inefficient use of modern optimization strategies, and this project helped address this limitation for the code’s application in E3SM. Under other funding, the Investigators of this project led an effort to develop a high-performance broadband radiation code, called RTE+RRTMGP, which is a completely restructured code that will take advantage of modern computational capabilities to enhance its performance while retaining the strengths and accuracy of the original code. Designed to perform over a range of computer architectures, RTE+RRTMGP makes extensive use of Fortran 2003 features to improve both its efficiency and its use of memory. RTE+RRTMGP is expected to be adopted in the next generation of RTE+RRTMGP. This project allowed for advancements to the code’s computational capabilities and adding new and enhanced features. This included revising RTE+RRTMGP to run on GPU processors, analysis of and improvements to the code’s timing, modifying the gas optics in RRTMGP to increase the code’s accuracy, extensive validation of computed fluxes, heating rates and forcings, generating a cloud optical property and vertical sampling capabilities for RTE+RRTMGP, implementing a fast and accurate longwave scattering capability, and groundwork for the inclusion of a capability to specify solar variability. Many of the accomplishment in this project necessitated significant collaboration with the E3SM development team. The result of this project was optimization of a key physical component (radiative transfer calculations) of E3SM, directly supporting E3SM’s overarching global modeling objectives. More broadly, this project provided overall advancements in the use of radiative transfer calculations in atmospheric modeling and simulation, particularly for climate.

58 GEOSCIENCES↗

Analysis of Space-Conditioning Loads in Commercial Buildings

Space conditioning end-uses, which include heating, cooling, and ventilation, represent a significant fraction of commercial building energy use, with a wide variety of heating and cooling technology options available in the market. In the interest of improving the overall efficiency of heating, ventilation and air-conditioning (HVAC) technologies, governments, utilities and private sector entities have implemented a variety of market transformation policies that aim to influence consumer purchase decisions. To evaluate the costs and benefits of such programs, analysts typically postulate a hypothetical default equipment choice, and compare it to one that provides comparable service with lower energy and/or power use. The corresponding reduced operating cost provides a benefit that offsets the potential higher cost of improved efficiency. Typically, life-cycle cost or cash-flow analyses are used to quantify the net economic benefit. These analyses require the capability to assess how a given equipment design would perform across a broad range of characteristics, both of the building and of the local weather. While these assessments can be performed using customized building simulations, it is generally not practical to develop and validate detailed building simulation code to cover all the potential variations of equipment design and installation. An alternative, and somewhat simpler, approach is to solely use detailed building simulations to generate time series of heating and cooling loads in commercial buildings. These loads can then be used as input to more detailed, stand-alone engineering models that simulate HVAC system performance under different equipment designs. This approach was used to evaluate a range of high-efficiency commercial packaged air conditioner design options for the Department of Energy’s Appliance and Equipment Standards Program (DOE-EERE 2015). While there may be some loss of precision relative to full simulation, the accuracy of this approach is sufficient for practical applications of cost-benefit analysis. This report describes the development of a database of commercial building heating and cooling loads, generated using the EnergyPlus software package, a whole building energy use model supported by the Department of Energy (DOE-EERE 2020a). EnergyPlus takes as input a set of configuration files that describe the building itself (size, zoning, envelope characteristics, etc.) and the various systems within it (HVAC, lighting, water heating, etc.). This analysis uses a publicly available collection of commercial reference buildings (CRB), comprised of sixteen building types and three vintages (DOE-EERE 2020b; Deru et al 2011). Each building is simulated in eighteen different locations, covering a wide range of climatic conditions. The prototype building description files assign the type of HVAC equipment used, and capacities across climate zones.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗