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 1,027 records · Page 57

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↗

Recent Updates in ETOE-2 and MC 2 -3

The MC 2 -3 code is a multigroup cross section generation code for fast reactor applications, developed by Argonne National Laboratory (ANL) under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, first, the cross section libraries for ENDF/B-VII.1 and ENDF/B-VIII.0 data were generated using the updated ETOE-2 code. Preliminary verification and validation tests of the ENDF/B-VII.1 MC 2 -3 library were performed with the selected fast reactor benchmark problems. Processing and verifying the libraries suggested that significant efforts would be required for thoroughly verifying the ENDF/B-VII.1 MC 2 -3 library and successfully processing the ENDF/B-VIII.0 MC 2 -3 library. Secondly, the cross section generation capability of MC 2 -3 was updated with the intermediate group lattice calculation and the equivalent Dancoff-factor cell (EDC) method to significantly improve the performance of a twodimensional assembly calculation using the method of characteristics (MOC). This effort was made useful in implementing and verifying the EDC method in the Griffin cross section API. Finally, we analyzed a SFR problem, for which unusually large deviations in core eigenvalues from Monte Carlo solutions were reported, by generating multigroup cross sections with MC 2 -3, performing core calculations with DIF3D-VARIANT, and analyzing cross section and eigenvalue results against Monte Carlo solutions. The analysis confirmed that the MC 2 -3/DIF3D solutions were in good agreement with Monte Carlo solutions, providing the correct process of accurately generating broad-group cross sections with MC 2 -3

97 MATHEMATICS AND COMPUTING↗

Doppler-Free Saturation Spectroscopy (DFSS) for Magnetic and Electric Field Measurements in an FRC plasma

The Cooperative Research and Development Agreement (CRADA) between Oak Ridge National Laboratory (ORNL) and TAE Technologies (TAE) focused on a feasibility study of Doppler-free saturation spectroscopy (DFSS) as a diagnostic strategy to measure the internal magnetic and electric fields of TAE’s C-2W field-reversed configuration (FRC) fusion device. Non-perturbative measurements of magnetic and electric fields in the core of FRCs can provide critical information needed to constrain equilibrium models and for the continued advancement of the fusion research program at TAE. DFSS is a well-established, laser-based absorption spectroscopy technique used to suppress inhomogeneous (Doppler) broadening linewidths and reveal Zeeman/Stark split hyperfine structures of target emission lines. The collaboration involved the use of ORNL’s suite of experimentally validated quantum mechanical modelling codes to simulate theoretical atomic spectra for a variety of machine parameters and laser beam configurations. TAE’s expertise was heavily used to provide measurements/models of the required machine parameters needed as inputs to the simulations. Results of an all-encompassing error analysis demonstrate the favourable capabilities of the diagnostic and provided invaluable insights as to how DFSS can be best used to accurately measure internal B/E-fields in C-2W.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Multi-scale fission product release model with comparison to AGR data

TRistructural ISOtropic (TRISO) particle fuel is central to several advanced, high-temperature reactor designs. Each particle consists of a fuel kernel encapsulated by three layers of carbon and ceramics that prevent the release of fission products and ensure physical integrity. Despite outstanding retention properties, fission product release has been observed from intact particles. To better understand and quantify fission product release from TRISO particles, a multiscale, mechanistic model of fission product transport is being developed by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. Previous work focused on silver (Ag) transport and improved Ag release predictions. The work described in this report builds on this experience to better understand cesium (Cs) transport in silicon carbide (SiC), the main barrier to the release of fission products. Atomistic simulations provide bulk and grain boundary (GB) Cs diffusivities in SiC, which are used by phase field simulations in the mesoscale code Marmot to determine the temperature, microstructure, and irradiation-dependent Cs diffusivity at the mesoscale in SiC. This approach attributes the different temperature regimes experimentally observed for Cs diffusivities in SiC to a transition from bulk-dominated diffusivity at high temperatures to a GB-dominated regime at low temperatures, providing new insight. The multiscale, mechanistic effective diffusivity is then implemented in the fuel performance code BISON and further validated by comparing Cs release predictions from Advanced Gas Reactor (AGR)-1 and AGR-2 post-irradiation measurements. The new model improves BISON’s predictability. This document also reports improvements made on Ag transport modeling by accounting for different GB types having different diffusivities. Moreover, this report details preliminary efforts to model palladium (Pd) attack of the SiC at the mesoscale using a phase field approach. Pd attack and its impact on accelerated Ag transport remains a misunderstood phenomenon, and we use the model to demonstrate that the formation of lamellae that has been observed in experiments can be explained by the reaction of Pd with SiC to form alternating layers of graphite and Pd 2 Si. This effort aims to improve our understanding of the reaction and eventually provide a model for BISON to account for Pd penetration and its effects on fission product release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Basic research in wake vortex alleviation using a variable twist wing

The variable twist wing concept was used to investigate the relative effects of lift and turbulence distribution on the rolled up vortex wake. Several methods of reducing the vortex strength behind an aircraft were identified. These involve the redistribution of lift spanwise on the wing and drag distribution along the wing. Initial attempts to use the variable twist wing velocity data to validate the WAKE computer code have shown a strong correlation, although the vorticity levels were not exactly matched.

Morris, D. J.↗

Transonic swept wings studied by the lifting-line theory

Transonic swept wings are analyzed as a lifting-line problem under a small-disturbance approximation. Basic concepts and principal results of the asymptotic theory are discussed. The study focuses on straight oblique wings and V-shaped swept wings, of which the local centerline curvature can be equated to zero. The three-dimensional (3-D) perturbation of the nonlinear component flow admits a similarity flow structure but requires that all wing sections are generated from a single airfoil profile; the reduced 2-D problems in this case are solved only once for all span stations. Examples of solutions involving high subcritical and slightly supercritical component flows are demonstrated and compared with surface pressure data from 3-D computer codes based on the full-potential equation (FLO 22). Except in the neighborhood of leading edges, where the small-disturbance assumption breaks down, and in the vicinities of wing tips and the symmetry plane, where neither the theory nor the 3-D codes may claim full validity, reasonable agreement is consistently found. The explicit results from the upwash analysis, along with the similarity flow structure, provides a rational approach to the control of 3-D effects in transonic aerodynamic design studies.

Cheng, H. K.↗

Addressing the hypersonic simulation problem

The hypersonic simulation problem is solved by the meshing of bench-marked experimental data with results from validated computational fluid dynamic codes. The example presented is for the reentry of the Space Shuttle Orbiter. Mach number effects were assessed by parametrically varying free-stream Mach number and angle of attack in a series of inviscid, perfect-gas computations carried out on a modified Orbiter geometry. Real-gas effects were determined by making calculations at specific points of the reentry trajectory using equilibrium air thermodynamics and comparing with corresponding perfect-gas computations. Viscous computations were also made for both the basic Orbiter body and control surfaces. A methodology is developed that permits the extrapolation of wind tunnel data to flight situations, providing a solution to the hypersonic simulation problem.

Griffith, B. J.↗

Simplified curve fits for the transport properties of equilibrium air

New, improved curve fits for the transport properties of equilibruim air have been developed. The curve fits are for viscosity and Prandtl number as functions of temperature and density, and viscosity and thermal conductivity as functions of internal energy and density. The curve fits were constructed using grabau-type transition functions to model the tranport properties of Peng and Pindroh. The resulting curve fits are sufficiently accurate and self-contained so that they can be readily incorporated into new or existing computational fluid dynamics codes. The range of validity of the new curve fits are temperatures up to 15,000 K densities from 10 to the -5 to 10 amagats (rho/rho sub o).

Srinivasan, S.↗

Direct simulations of chemically reacting turbulent mixing layers, part 2

The results of direct numerical simulations of chemically reacting turbulent mixing layers are presented. This is an extension of earlier work to a more detailed study of previous three dimensional simulations of cold reacting flows plus the development, validation, and use of codes to simulate chemically reacting shear layers with heat release. Additional analysis of earlier simulations showed good agreement with self similarity theory and laboratory data. Simulations with a two dimensional code including the effects of heat release showed that the rate of chemical product formation, the thickness of the mixing layer, and the amount of mass entrained into the layer all decrease with increasing rates of heat release. Subsequent three dimensional simulations showed similar behavior, in agreement with laboratory observations. Baroclinic torques and thermal expansion in the mixing layer were found to produce changes in the flame vortex structure that act to diffuse the pairing vortices, resulting in a net reduction in vorticity. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layers were shown to result from vorticity generation by baroclinic torques.

Metcalfe, Ralph W.↗