Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analysis codes”

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 847 records · Page 47

Status of the NEAMS and ARC fast reactor tools integration to the NEAMS Workbench

The Workbench initiative was launched in FY-2017 within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to facilitate the transition from conventional tools to high-fidelity tools. The NEAMS Workbench provides a common user interface for model creation, real-time validation, execution, output processing, and visualization for integrated codes. The integration of the Argonne Reactor Computation (ARC) suite of codes into the NEAMS Workbench through the PyARC module was initiated in FY-2017. The ARC codes, which focus on fast reactor multiphysics analysis, contain both legacy codes like DIF3D and REBUS-3 that were developed with over 30 years of experience, and newer NEAMS additions like MC 2 -3, PERSENT, and PROTEUS. Recent work expended this integration to other tools to support the U.S. fast reactor community such as DASSH for sub-channel thermal-hydraulics, Griffin for high-fidelity deterministic transport calculations, and OpenMC for Monte Carlo simulations (with Shift integration planned for FY-2023). The integration of the “extended ARC” suite of codes into the NEAMS Workbench interface relies on the PyARC and PyGriffin modules to handles the pre- and post-processing of these codes input, and the runtime environment. The PyARC module together with the NEAMS Workbench interface are both released under Open Source Software licenses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Independent Analyses of Antares R1 Core Design

This report summarizes the reactor analysis work performed by Oak Ridge National Laboratory (ORNL) for Antares Nuclear Incorporated's R1 Mark-1 heat pipe reactor design. The work was performed under a collaboration through the Department of Energy GAIN Nuclear Energy Voucher program. The purpose of this work is to perform an independent reactor analysis of the Antares reactor design and, where possible, compare ORNL's results with the results obtained by Antares. The overarching goal is to provide Antares with an independent review and calculation of their design, thereby contributing to Antares' mission to further their reactor design concept. A report with all ORNL technical results and proprietary details was provided to Antares Nuclear Incorporated separately; the present report provides a summary of the work completed, and all proprietary information is omitted. The independent reactor analysis was performed using both the SCALE code system for neutronics analysis and Flownex for thermal hydraulics analysis.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Advanced modeling and simulation of research reactors using dynamic mode decomposition

Full text of publication follows. Due to the ever-increasing safety requirements, the current trend of nuclear reactor analysis is shifting towards high-fidelity multi-physics models, which have a very high computational cost and modelling complexity. As the cost of even a single model run makes it impossible to analyse the behaviour and performance of these models on large-scale commercial plants, it has become even more significant to provide suitable benchmarks to validate and test them extensively. In this sense, research reactors offer a promising solution for the initial validation of high-fidelity models, as they are significantly smaller than commercial reactors and their characteristics are well known. In particular, the reactors of the TRIGA family have been used to assess and validate models and methods for Generation-IV designs, as they have some similar features (such as the dominance of natural convection as cooling mechanism and the difficulties in performing sub-channel analysis using standard codes). Still, the computational requirements of high-fidelity models make them unsuitable for real-time analysis, even following their assessment on research reactors. In this sense, Model Order Reduction (MOR) techniques give an additional strategy to reduce the computational cost of high-fidelity models (whilst preserving sufficient accuracy). In particular, this work focuses on Dynamic Mode Decomposition (DMD), a non-intrusive MOR technique that aims at representing models with explicit temporal dynamics by extracting the time-varying characteristics and the governing structures based only on a set of available data, thus without needing any underlying knowledge of the governing equations. In addition, DMD also computes a low-dimensional surrogate of the dynamic matrix of the system, making it suited for stability analysis and real-time evaluations. This work focuses on the application and validation of the DMD method on the Computational Fluid-Dynamics (CFD) model TRIGA Mark II reactor, also discussing in detail the potentiality of this algorithm as an advanced modelling tool for nuclear reactor analysis. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Data and code for Daily and Multi-Day Extreme Rainfall Analysis Under Future Climates Using Stochastic Storm Transposition and NEX-GDDP-CMIP6 Over CONUS

This data package provides inputs, codes, and outputs for a comprehensive analysis of projected changes in extreme precipitation across 10 regions of the continental United States, using 34 downscaled Earth System Models (ESMs) from the NASA Earth Exchange Global Daily Downscaled Projections, Coupled Model Intercomparison Project Phase 6 (NEX-GDDP-CMIP6) dataset. These models are part of the Coupled Model Intercomparison Project Phase 6 (CMIP6), a coordinated climate modeling framework widely used to assess climate change impacts. The analysis applies a stochastic storm transposition method to quantify changes in extreme rainfall under two Shared Socioeconomic Pathway (SSP) climate scenarios—SSP2-4.5 (moderate emissions) and SSP5-8.5 (high emissions)—compared to historical conditions (1995–2014 vs. 2081–2100). The dataset includes rainfall depth estimates for extreme events with return periods from 2 to 500 years across multiple storm durations (1, 3, and 5 days) for each of the 10 U.S. regions. Weighted ensemble statistics are derived from individual ESM performance against historical precipitation patterns, enabling robust uncertainty quantification through both sign-based and permutation-test-based model agreement assessments. Key analyses address: (1) relative changes in extreme precipitation for each climate scenario, (2) differences between SSP scenarios (SSP5-8.5 vs. SSP2-4.5), (3) contrasts between rare and frequent events, and (4) variations between multi-day and daily storm durations. The workflow produces ensemble statistics—median, 5th, 25th, 75th, and 95th percentiles—along with model agreement metrics that identify regions and event types with robust climate change signals. The dataset includes: processed rainfall depth outputs (netCDF format) from the RainyDay Python package, ESM weights from historical performance evaluation using DayMet observations, ensemble statistics across all storm dimensions, and figures summarizing key findings.

54 ENVIRONMENTAL SCIENCES↗

The STAT7 Code for Statistical Propagation of Uncertainties In Steady-State Thermal Hydraulics Analysis of Plate-Fueled Reactors

The STAT code was written to automate many of the steady-state thermal hydraulic safety calculations for the MIT research reactor, both for conversion of the reactor from high enrichment uranium fuel to low enrichment uranium fuel and for future fuel re-loads after the conversion. A Monte-Carlo statistical propagation approach is used to treat uncertainties in important parameters in the analysis. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux or departure from nucleate boiling or onset of flow instability; but additional margin is obtained by basing the limiting safety settings on avoiding onset of nucleate boiling. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one stripe of a fuel element. The stripes run the length of the fuel, from the bottom to the top. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. By running STAT7 multiple times, full core analysis has been performed by analyzing the margin to ONB for each axial node of each stripe of each plate of each element in the core.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development of a LOCA Experimental Benchmark for BISON - Milestone L2:FMC.P15.08

Since 2013 CASL has been planning and working towards a PWR LOCA Challenge Problem. Since BISON is the CASL fuel rod simulation platform, it plays a central role in that activity. The goal of the FY-17 CASL milestone reported here was to develop an Experimental Benchmark for LOCA analysis using BISON. Activities encompassed both code development and validation. This report provides a current snapshot of BISON’s capability for LOCA behavior including: 1) a summary of code extensions to facilitate accident analysis, 2) a series of separate effects tests, and 3) initial code validation to complex integral rod LOCA behavior. Conclusions from each of these main activities are summarized below. The key material and behavior models required to address transient high-temperature phenomena occurring during LOCAs in a standard PWR have now been implemented in BISON. These apply specifically to UO2 fuel, Zircaloy cladding and water coolant. During FY-17 important new capability to address axial UO2 fuel relocation (for 1.5D geometry) and account for oxidation energy deposition in cladding, were included. Planned future development efforts include extending the axial fuel relocation model from 1.5D to 2D/3D and improving BISON’s fission gas release model to include transient gas release associated with the high burnup structure (HBS) in high burnup fuel. Additionally, investigation of potential anisotropic creep behavior in Zircaloy cladding is planned, assuming sufficient experimental data are available to support this effort. A substantial number of separate effects validation cases (42 tests from 3 experimental series) have been completed to compare BISON predictions to measured ballooning and burst behavior for Zircaloy cladding. Such experiments include a wide variety of pressures, temperatures and loading rates. In general, BISON predictions of burst temperature, pressure and burst time are very reasonable. For one experimental series, however, involving both very high temperatures and strain rates, BISON systematically over-predicts the cladding hoop strain. Investigation of this discrepancy is an important and planned activity. With another experiment an effort was made to investigate 3D cladding response due to an azimuthal temperature variation. Results indicated 3D effects are potentially important in fuel rod analysis during LOCAs and will be further investigated in the future. BISON validation to a series of integral fuel rod experiments has also been completed. These experiments involve all fuel and cladding phenomena relevant to LOCA conditions, and can include complexities associated with irradiated fuel relative to fresh fuel. Such experiments also generally include complex thermal-hydraulic boundary conditions. Four experiments (6 rods) have been considered to date including simulated fuel (ZrO2) and both fresh and high-burnup UO2. Test rods ranged from rodlets to full length commercial PWR fuel rods. As with the separate effects experiments, BISON predictions of burst temperature, pressure and burst time are generally very reasonable. Comparisons to cladding peak strain and rod outer diameter axial profiles are less satisfactory, and identify material models and possibly modeling approximations (e.g., 2D-RZ vs 3D geometries) requiring additional investigation. Validation of BISON for integral rod LOCA behavior is by no means complete, with additional cases planned.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SAM Developments for MSR Modeling (FY21)

An advanced system analysis tool, SAM, is under development for advanced non-LWR reactor safety analysis, including molten salt reactors (MSR). To support the development and utilization of the SAM code for MSR safety analysis and licensing, continuous efforts have been devoted to enhancing code capabilities and updating reference models for the MSRs. This report documents the FY21 progress in SAM code development, capability enhancements, and reference model development to support transient safety analysis of MSRs, including code enhancements in reactor kinetics and reactivity feedback modeling for liquid fuel reactors, updates of the Molten Salt Reactor Experiment (MSRE) primary system model, developments and updates of the molten salt fast reactor (MSFR) model based on EVOL design, and implementation of a drift flux model for modeling gas transport in MSR systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Radiation shielding analysis of the concrete hutch of a long beamline of the Advanced Photon Source Upgrade project

The radiation shielding analysis of the external concrete enclosure of a feature long beamline of the Advanced Photon Source Upgrade project is described in this work. The beamline has 2 undulators that can operate in two different configurations and the source terms are estimated for the worst configuration using the SRW code in the OASYS package as well as the STAC8 code. The shielding analysis is carried out using the STAC8 code, and for the case of pink beam, with the FLUKA code also. Effective dose rates on contact with the outside surface of the concrete walls for pink, reflected full spectrum and monochromatic beams are calculated. For the monochromatic beams, 2 sets of discrete bandwidths (BW) from the Double Crystal Monochromator (DCM) and the full reflectance of the Double Multilayer Monochromator (DMM) are used. Further, due to the size of the direct beam and the energy domain of interest in the direct and scattered fields, the importance of the equivalent dose to the lens of the eye is discussed with respect to the shielding calculations of such facilities. The dose rates outside the concrete hutch in the lateral directions calculated using STAC8 is about 50% higher than the FLUKA results for the case of pink beam but may be considered as a reasonable agreement due to the differences in the calculational methodologies. Those dose rates in the lateral and forward directions are higher for the pink beam to allow full occupancy from the radiation protection point of view. For the full spectrum after 3 reflections, the effective dose rates are less than 0.5 Sv/h in the lateral direction but is higher in the forward region below about 22° with respect to the beam direction. For the DCM beams, except for the scenario with one mirror and 0.1%BW, all dose rates are less than 0.5 Sv/h outside the concrete hutch in the lateral directions. Here, the higher harmonics are what contributes to the dose outside the shield and the bandwidths assumed to derive the photon flux thus becomes important. The dose rates with the DMM are generally lower compared to the 0.1% BW DCM beam but are higher than the results obtained with XOP calculated bandwidths.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Promoter analysis reveals globally differential regulation of human long non-coding RNA and protein-coding genes

Transcriptional regulation of protein-coding genes is increasingly well-understood on a global scale, yet no comparable information exists for long non-coding RNA (lncRNA) genes, which were recently recognized to be as numerous as protein-coding genes in mammalian genomes. We performed a genome-wide comparative analysis of the promoters of human lncRNA and protein-coding genes, finding global differences in specific genetic and epigenetic features relevant to transcriptional regulation. These two groups of genes are hence subject to separate transcriptional regulatory programs, including distinct transcription factor (TF) proteins that significantly favor lncRNA, rather than coding-gene, promoters. We report a specific signature of promoter-proximal transcriptional regulation of lncRNA genes, including several distinct transcription factor binding sites (TFBS). Experimental DNase I hypersensitive site profiles are consistent with active configurations of these lncRNA TFBS sets in diverse human cell types. TFBS ChIP-seq datasets confirm the binding events that we predicted using computational approaches for a subset of factors. For several TFs known to be directly regulated by lncRNAs, we find that their putative TFBSs are enriched at lncRNA promoters, suggesting that the TFs and the lncRNAs may participate in a bidirectional feedback loop regulatory network. Accordingly, cells may be able to modulate lncRNA expression levels independently of mRNA levels via distinct regulatory pathways. Our results also raise the possibility that, given the historical reliance on protein-coding gene catalogs to define the chromatin states of active promoters, a revision of these chromatin signature profiles to incorporate expressed lncRNA genes is warranted in the future.

59 BASIC BIOLOGICAL SCIENCES↗

Multi-Physics Simulations of Heat Pipe Micro Reactor

System Analysis Module (SAM) is under development at Argonne National Laboratory as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. This report presents the development of new capabilities for modeling the heat pipe type reactor systems. The need for power at remote locations away from a reliable electrical grid is an important niche for nuclear energy. Heat pipe-cooled fast-spectrum nuclear reactors are well suited for these applications. The key feature of the heat pipe reactors is the use of heat pipes for heat removal from the reactor core. The heat pipe makes use of the phase change of the working fluid and transports a large amount of heat from the evaporator to the condensation end with very small temperature drops. In contrast to the traditional nuclear reactor system that makes use of pumped loop for extracting the thermal power, the heat pipe reactors make use of hundreds of heat pipes for removing the thermal power (including the decay heat) passively. This could potentially significantly improve the reliability and safety of the reactor systems. The essential part in the analysis of a heat pipe type reactor is the modeling of heat transport inside the heat pipe. The capability of SAM was extended to enable the modeling of the conventional heat pipe. The full suite of non-LWR codes for confirmatory analysis at NRC is known as the Comprehensive Reactor Analysis Bundle (CRAB). It makes use of existing NRC codes, and integrates them with several codes developed through the DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this report, the heat pipe micro reactor is simulated using three submodules under CRAB. The MAMMOTH module is used to simulate the reactor kinetics behavior of the micro reactor; the SAM module is used to simulate the heat conduction in the reactor core and heat removal through the heat pipe heat exchangers and Reactor Cavity Cooling System (RCCS); and the MOOSE Tensor Mechanics module is used to simulate the thermal expansion of the reactor cores. The different sub-models are coupled together using MOOSE’s MultiApp system and executed using the BlueCrab application. The new capability has been demonstrated by a steady-state operation analysis, a failure of a single central heat pipe transient analysis, and a loss of heat sink transient analysis. The fully coupled model is shown to work well.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Atomic and Molecular Physics Effects on Ion Flow into Hot Surfaces

Research during this project at the University of Wisconsin Fusion Technology Institute (UW FTI) on ion and neutral flow through an arbitrary, monotonic potential difference created by nearly transparent electrodes accomplished the following: (1) developed and implemented an integral equation approach for atomic physics effects in helium plasmas; (2) extended the analysis to coupled integral equations that treat atomic and molecular deuterium ions and neutrals; (3) implemented the key deuterium and helium atomic and molecular cross sections; (4) added negative ion production and related cross sections; and (5) benchmarked the code against experimental results. The analysis and codes treat the species D 0 , D 2 0 , D + , D 2 + , D 3 + , D - and, separately at present, He 0 and He + . Extensions enhanced the analysis and related computer codes to include He ++ ions plus planar and cylindrical geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of an object-oriented ORIGEN for advanced nuclear fuel modeling applications

The ORIGEN package serves as the core depletion and decay calculation module within the SCALE code system. A recent major re-factor to the ORIGEN code architecture as part of an overall modernization of the SCALE code system has both greatly enhanced its maintainability as well as afforded several new capabilities useful for incorporating depletion analysis into other code frameworks. This paper will present an overview of the improved ORIGEN code architecture (including the methods and data structures introduced) as well as current and potential future applications utilizing the new ORIGEN framework. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Summary of improvements in modeling for fidelity and performance in MAMBA during FY18

This milestone report summarizes the FMC work on improvements to the MAMBA modeling methodology. The work was performed in close collaboration with the PHI as well as THM focus area and it was undertaken to improve both the physical fidelity and computational performance of MAMBA. We first describe the path towards identifying and then resolving issues in the MAMBA code and the underlying physics models that were still remaining at the beginning of FY18, after the code re-write in the previous year. These were primarily related to boiling and species transport, especially boron, within the CRUD. The code problems originated in the MAMBA3D code, from which the development of the new re-written code started, but they were only identified during testing of the new MAMBA code. We purposefully discuss a number of alternative explanations and solutions to the main issue identified, since the two-phase character of fluid transport in the CRUD and the related mass transport are still sources of model uncertainty that we continue to investigate. The hands-on code development referenced in this report was done under PHI and this milestone mainly provided support in the form of analysis, literature review, code review and feedback to the code developers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

High conversion Th-U{sup 233} fuel assembly for current generation of PWRs

This paper presents a preliminary design of a high conversion Th-U{sup 233} fuel assembly applicable for current generation of Pressurized Water Reactor (PWRs). The considered fuel assembly has a typical 17 x 17 PWR lattice. However in order to increase the conversion of Th{sup 232} to U{sup 233}, the assembly was subdivided into the two regions called seed and blanket. The central seed region has a higher than blanket U{sup 233} content and acts as a neutron source for the peripheral blanket region. The latest acts as a U{sup 233} breeder. While the seed fuel pins have a standard dimensions the blanket fuel radius was increased in order to reduce the moderation and to facilitate the resonance neutron absorption in blanket Th{sup 232}. The U{sup 233} content in the seed and blanket regions was optimized to achieve maximal initial to discharged fissile inventory ratio (FIR) taking into account the target fuel cycle length of 12 months with 3-batch reloading scheme. In this study the neutronic calculations were performed on the fuel assembly level using Helios deterministic lattice transport code. The fuel cycle length and the core k{sub eff} were estimated by applying the Non Linear Reactivity Model. The applicability of the HELIOS code for the analysis of the Th-based high conversion designs was confirmed with the help of continuous-energy Monte-Carlo code SERPENT. The results of optimization studies show that for the heterogeneous seed and blanket (SB) fuel assembly the FIR of about 0.95 can be achieved. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Identification of RESNET HERS Index Values Corresponding to Minimal Compliance with the IECC Performance Path

Pacific Northwest National Laboratory analyzed the relationship between the Residential Energy Services Network (RESNET) Home Energy Rating System (HERS) Index and the traditional simulation-based Performance Path used in the International Energy Conservation Code (IECC). The analysis evaluates, for a single-family residence with various characteristics, the ranges of HERS Index values that would imply compliance with the 2012 IECC Performance Path. Several building characteristics considered likely to result in quantifiable differences in the outcomes of the two approaches, or otherwise believed to be of interest to code developers and policy makers, are considered in the analysis.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Thermomechanical analysis and modeling of involute-shaped fuel plates using the Cheverton–Kelley experiments for the High Flux Isotope Reactor

Three research reactors with involute-shaped fuel plates are pursuing conversion from highly enriched uranium to low-enriched uranium fuel. Various core design and safety evaluation studies are essential to assess the feasibility of the conversion. The use of 3D computational multiphysics codes is being explored in these analyses and therefore they must undergo a thorough evaluation and quality assurance process due to their potential impact on nuclear safety. Here, the Cheverton and Kelley physical tests performed in the late 1960s to investigate the deflections of HFIR’s outer plate under uniform pressure and temperature fields are simulated by employing commercially available computational codes, with the goals to (1) verify and validate the models and numerical solvers implemented in the codes for thermomechanical analysis of involute reactor plates and (2) to develop a benchmark computational test to evaluate future versions of existing software or newly developed computational codes. The results of the simulations showed good agreement with each other as well as against the Cheverton–Kelley experimental data. Some minor deviations were observed for a few multiphysics cases and their potential origins and impact on the analysis results is investigated in the paper. The validated models increase the confidence in using multiphysics codes to evaluate existing or new LEU designs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Three-dimensional analysis using polygons in particle and heavy ion transport code system (PHITS) - 25504

Recently a function to treat tetrahedral-mesh geometry was introduced in Particle and Heavy Ion Transport code System (PHITS). Owing to this function, three-dimensional analysis using polygons became possible. The work flow of the three-dimensional analysis using a particular polygon object was demonstrated, converting the data into tetrahedral mesh with the tetrahedral-mesh generator TetGen, computing the radiation transport with PHITS, and finally analyzing with three-dimensional data analysis application ParaView. (author)

61 RADIATION PROTECTION AND DOSIMETRY↗

Flow-induced vibration analysis of a helical coil steam generator experiment using large eddy simulation

This paper describes a numerical study of flow-induced vibration in a helical coil steam generator experiment conducted at Argonne National Laboratory in the 1980s. In the experiment, a half-scale sector model of a steam generator helical coil tube bank was subjected to still and flowing air and water, and the vibrational characteristics were recorded. The research detailed in this document utilizes the multi-physics simulation toolkit SHARP developed at Argonne National Laboratory, in cooperation with Lawrence Livermore National Laboratory, to simulate the experiment. SHARP uses the spectral element code Nek5000 for fluid dynamics analysis and the finite element code DIABLO for structural analysis. The flow around the coil tubes is modeled in Nek5000 by using a large eddy simulation turbulence model. Transient pressure data on the tube surfaces is sampled and transferred to DIABLO for the structural simulation. The structural response is simulated in DIABLO via an implicit time-marching algorithm and a combination of continuum elements and structural shells. Tube vibration data (acceleration and frequency) are sampled and compared with the experimental data. Currently, only one-way coupling is used, which means that pressure loads from the fluid simulation are transferred to the structural simulation but the resulting structural displacements are not fed back to the fluid simulation

CFD↗