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 145 records · Page 8

RELAP-7 Development Updates

RELAP-7 is a nuclear systems safety analysis code being developed at the Idaho National Laboratory, and is the next generation tool in the RELAP reactor safety/systems analysis application series. RELAP-7 development began in 2011 to support the Risk Informed Safety Margins Characterization (RISMC) Pathway of the Light Water Reactor Sustainability (LWRS) program. The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical methods, and physical models in order to provide capabilities needed for the RISMC methodology and to support nuclear power safety analysis. The code is being developed based on Idaho National Laboratory’s modern scientific software development framework – MOOSE (the Multi-Physics Object-Oriented Simulation Environment). The initial development goal of the RELAP-7 approach focused primarily on the development of an implicit algorithm capable of strong (nonlinear) coupling of the dependent hydrodynamic variables contained in the 1-D/2-D flow models with the various 0-D system reactor components that compose various boiling water reactor (BWR) and pressurized water reactor nuclear power plants (NPPs). During Fiscal Year (FY) 2015, the RELAP-7 code has been further improved with expanded capability to support boiling water reactor (BWR) and pressurized water reactor NPPs analysis. The accumulator model has been developed. The code has also been coupled with other MOOSE-based applications such as neutronics code RattleSnake and fuel performance code BISON to perform multiphysics analysis. A major design requirement for the implicit algorithm in RELAP-7 is that it is capable of second-order discretization accuracy in both space and time, which eliminates the traditional first-order approximation errors. The second-order temporal is achieved by a second-order backward temporal difference, and the one-dimensional second-order accurate spatial discretization is achieved with the Galerkin approximation of Lagrange finite elements. During FY-2015, we have done numerical verification work to verify that the RELAP-7 code indeed achieves 2nd-order accuracy in both time and space for single phase models at the system level.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Assessment of the MACCS Code Applicability for Nearfield Consequence Analysis

Three codes are used for comparisons in this report to evaluate the adequacy of MACCS in the nearfield, AERMOD, ARCON96 and QUIC. Test cases were developed to give a broad range of weather conditions, building dimensions and plume buoyancy. Based on the comparisons of MACCS with AERMOD, ARCON96 and QUIC across the test cases, the following observations are made: MACCS calculations configured with point-source, ground-level, nonbuoyant plumes provide nearfield results that bound the centerline, ground-level air concentrations from AERMOD, ARCON96, and QUIC. MACCS calculations with ground-level, nonbuoyant plumes that include the effects of the building wake (area source) provide nearfield results that bound the results from AERMOD and QUIC and the results from ARCON96 at distances >250 m. If using a point-source is too conservative and it is desired to bound the results from all three codes, another alternative is to use area source parameters in MACCS that are less than the standard values, i.e., an area source intermediate between the standard recommendation and a point source. All these options provide results from MACCS that are bounding for the test cases evaluated. Based on these observations, it appears that MACCS is adequate for use in nearfield calculations, given the correct parameterization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An efficient numerical model for predicting residual stress and strain in parts manufactured by laser powder bed fusion

Abstract Computational modeling of additively manufactured structures plays an increasingly important role in product design and optimization. For laser powder bed fusion processes, the accurate modeling of stress and distortion requires large amount of computational cost due to very localized heat input and evolving complex geometries. The current study takes advantage of a graphics processing unit accelerated explicit finite element analysis code and approximated heat conduction analysis to predict the macroscopic thermo-mechanical behavior in laser selective melting. Adjacent layers and tracks were lumped to reduce the number of time steps and elements in the finite element model. The effects of track and layer grouping on prediction accuracy and solution efficiency are investigated to provide a guidance for a cost-effective simulation. Thin-wall builds from Inconel alloy 625 (IN625) powders were simulated by applying the developed modeling approach to get the detailed residual stress and distortion at a computational speed 50 times higher than conventional approach. Under repeated heating and cooling cycles, a high tensile stress was produced near surfaces of a build due to a larger shrinkage on surface than that in central area. It is also shown that horizontal stresses concentrate near the root and top layers of the IN625 build. The predicted residual elastic strain distribution was validated by the experimental measurement using x-ray synchrotron diffraction.

36 MATERIALS SCIENCE↗

Implementation and Verification of the Sesame Equation of State Database into Abaqus

Towards ability to perform advanced engineering calculations of shock driven processes requiring an accurate prediction of temperature inside the shock, the Sesame equation of state database has been implemented into Abaqus explicit, a finite element analysis code commonly used in the analysis community. The interface algorithm and supporting details are discussed, and results from a simple benchmark problem comparing Abaqus predictions to other previously validated calculations and published data are presented. Features of the Sesame-Abaqus interface are also discussed and a sample input file is provided for the interested user.

42 ENGINEERING↗

Validation of Numerical Tools for Calculating Reactivity Feedback in Sodium Fast Reactors Using SEFOR Experimental Data

The Southwest Experimental Fast Oxide Reactor (SEFOR) was an experimental sodium-cooled fast breeder reactor operated from 1969 to 1972 with experiments designed to measure Doppler reactivity feedback in a wide temperature range from around 350 °F to temperatures approaching the melting point of mixed oxide fuel of around 5000 °F, providing valuable data for code validations. Co-supported by the Department of Energy (DOE) Fast Reactor Program (FRP) and the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the SEFOR benchmark project focused on using the experimental data to validate numerical tools that are used in industry and academia to design and license sodium-cooled fast reactors (SFRs). By the end of FY-25, substantial progress was achieved in the SEFOR benchmark study. A variety of numerical tools commonly used for modeling SFRs were applied to develop models for SEFOR core configurations I-D, I-E, I-I, and I-J. These included Monte Carlo codes such as MCNP, Serpent, and Shift; deterministic codes such as the legacy Argonne Reactor Computation (ARC) suite and the high-fidelity NEAMS code Griffin; and the system analysis code SAS4A/SASSYS-1 (SAS). Using these models, both SEFOR zero-power experiments and power-ascending tests were successfully simulated. Comparisons were performed against experimental measurements of core criticalities, reflector worth, kinetics parameters (Λ/βeff), isothermal reactivity feedback (from 350 °F to 760 °F at zero power), and power-ascending reactivity feedback (as power increased from 0.4 MW to 17 MW). In general, these comparisons demonstrated very good agreement between numerical results and experimental data. In Fiscal Year 26 (FY-26), the SEFOR benchmark project will continue to address the modeling issues identified in FY-25. Effort will focus on the simulation of reactivity insertion transients in SEFOR core II using the ARC/SAS model. Future work will also focus on incorporating BISON into the SEFOR core modeling process to enable the first Multiphysics simulations of the isothermal tests based on the MOOSE framework.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Used Nuclear Fuel-Storage, Transportation & Disposal Analysis Resource and Data System (UNF-ST&DARDS)

This will be licensed via RSICC. A new, integrated data and analysis system has been designed to simplify and automate the performance of accurate and efficient evaluations for characterizing the input to the overall nuclear waste management system -UNF-Storage, Transportation & Disposal Analysis Resource and Data System (UNF-ST&DARDS). A relational database within UNF-ST&DARDS provides a standard means by which UNF-ST&DARDS can succinctly store and retrieve modeling and simulation (M&S) parameters for specific spent nuclear fuel analysis. A library of various analysis model templates provides the ability to communicate the various set of M&S parameters to the most appropriate M&S application. Interactive visualization capabilities facilitate data analysis and results interpretation. UNF-ST&DARDS current analysis capabilities include (1) assembly-specific depletion and decay, (2) and spent nuclear fuel cask-specific criticality and shielding. Currently, UNF-ST&DARDS uses SCALE nuclear analysis code system for performing nuclear analysis.

Banerjee, Kaushik↗

Verification and Validation Activities of Molten Salt Reactors Multiphysics Coupling Schemes at Idaho National Laboratory

This paper presents the latest verification and validation activities in molten salt reactor modeling and simulation performed at Idaho National Laboratory. Multiphysics solutions are obtained by coupling the neutronics code Griffin, the thermal hydraulics code Pronghorn, and the system analysis code SAM, under the MOOSE framework. We present various multiphysics coupling schemes with these codes for molten salt reactor problems and provide verification and validation results. First, we present verification test results of the Griffin-Pronghorn coupled scheme for the CNRS benchmark. Then validation test results are presented for the Griffin-SAM coupled scheme for the pump startup and coast down transients of the Molten Salt Reactor Experiment. Finally, the Griffin-Pronghorn-SAM coupled scheme is demonstrated for the Molten Salt Reactor Experiment reactivity insertion transient using a domain-overlapping coupling algorithm between Pronghorn and SAM. The results of these various coupling schemes demonstrate the ability to capture the effect of fuel flow and the various feedback mechanisms important to MSRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Data Analysis Framework for Earth System Simulation within an In-Situ Infrastructure

This paper presents a generic procedure to implement a scalable and high performance data analysis framework for large-scale scientific simulation within an in-situ infrastructure. It demonstrates a unique capability for global Earth system simulations using advanced computing technologies (i.e ., automated code analysis and instrumentation), in-situ infrastructure (i.e ., ADIOS) and big data analysis engines (i.e ., SciKit-learn). This paper also includes a useful case that analyzes a globe Earth System simulations with the integration of scalable in-situ infrastructure and advanced data processing package. Finally, the in-situ data analysis framework can provides new insights on scientific discoveries in multiscale modeling paradigms.

97 MATHEMATICS AND COMPUTING↗

Parametric Evaluation of SiC/SiC Composite Cladding with UO 2 Fuel for LWR Applications: Fuel Rod Interactions and Impact of Nonuniform Power Profile in Fuel Rod

SiC/SiC composites are candidates for accident tolerant fuel cladding in light water reactors. In the extreme nuclear reactor environment, SiC-based fuel cladding will be exposed to neutron damage, significant heat flux, and a corrosive environment. Here, to ensure reliable and safe operation of accident tolerant fuel cladding concepts such as SiC-based materials, it is important to assess thermo-mechanical performance under in-reactor conditions including irradiation and realistic temperature distributions. The effect of non-uniform dimensional changes caused by neutron irradiation with spatially varying temperatures, along with the closing of the fuel–cladding gap, on the stress development in the cladding over the course of irradiation were evaluated. The effect of non-uniform circumferential power profile in the fuel rod on the mechanical performance of the cladding is also evaluated. These analyses have been performed using the BISON fuel performance modeling code and the commercial finite element analysis code Abaqus. A constitutive model is constructed and solved numerically to predict the stress distribution in the cladding under normal operating conditions. The dependence of dimensions and thermophysical properties on irradiation dose and temperature has been incorporated into the models. Initial scoping results from parametric analyses provide time varying stress distributions in the cladding as well as the interaction of fuel rod with the cladding under different conditions of initial fuel rod-cladding gap and linear heat rate. It is found that a non-uniform circumferential power profile in the fuel rod may cause significant lateral bowing in the cladding, and motivates further analysis and evaluation.

Accident tolerant fuel↗

Delayed Gamma Signature Accounting for Field Line Transport of Debris Ions after a Nuclear Explosion

Delayed gammas (DGs) are produced by the decay of radioactive fission fragments after a nuclear explosion. Currently, the time and energy profile of the delayed gamma signal is propagated in DIORAMA to the GEO satellite sensors by assuming a point source at the burst location for the production of the delayed gammas. However, in reality the production of delayed gammas is spread over some spatial extent. First, the debris ions are confined within an expanding region by the magnetic field and/or surrounding plasma. This expanding confinement bubble eventually collapses, at which point the debris ions can be trapped on magnetic field lines and transported away from the original burst location. In this work the effects of debris ion transport on the associated time profile of delayed gammas at a GEO satellite is presented. The first step was the development of the Delayed Gamma Source Model (DGSM), which is based on the analytical Electron Source Model (ESM) code developed at LANL that models the spatial distribution of energetic electrons produced after a high-altitude nuclear explosion (HANE) for use as a source term in LANL artificial belt diffusion code. This model creates a spatial grid in latitude, longitude, and altitude of the distribution of radioactive debris ions. The output of this code is combined with a separate analysis code also described, which produces the DG time profile signature at a GEO sensor faceplate when accounting for the spatial spread of delayed gamma production during field line transport. Some simple case studies are presented to show how the DG signature at GEO calculated with this higher-fidelity model and a simple point-source model, as is currently assumed in DIORAMA, compare.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling and Simulation Development Pathways to Accelerating KP-FHR Licensing (Final Report)

This project assembles a strong U.S. industry and national laboratory team to complete scope of work. Kairos Power (KP), headquartered in Alameda, CA, is the leader of this effort and has built an internal team of highly competent engineers and managers with extensive combined experience in nuclear power, conventional power, product development, and licensing. INL, ANL, and LANL bring unique capabilities in advanced reactor R&D and licensing. The project funding source is the result of FOA No. 0001817, U. S. Industry Opportunities for Advanced Nuclear Technology Development. There has been on-going work and this Access Cooperative Research and Development Agreement (CRADA) will cover the remaining work scope of FOA 0001817. KP is implementing innovative strategies that can reduce the cost and accelerate the initial demonstration of the Kairos Power Fluoride-salt-cooled, High-temperature Reactor (KP-FHR) to meet the needs of the U.S. electricity market by 2030. Licensing of the KP-FHR could be significantly accelerated using advanced computing methods with sufficient predictive capabilities to be able to extrapolate potential response of the structure in different scenarios. However, currently used computational methods heavily rely on empirical fits and cannot be used for extrapolation. The scope focuses on improving the modeling capability in the NEAMS Grizzly structural mechanics code and the NEAMS BISON fuel performance code and the NEAMS SAM systems analysis code. This work leverages the expertise and know-how gathered in three DOE National Laboratories – INL, ANL, and LANL.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Picometer-Precision Atomic Position Tracking through Electron Microscopy

The modern aberration-corrected scanning electron microscopes have successfully achieved direct visualization of atomic columns with sub-angstrom resolution. With this significant progress, advanced image quantification and analysis is still at its early stages. In this work, we present the complete pathway for the metrology of atomic resolution STEM images. This includes: 1) tips for acquiring high-quality STEM images; 2) denoising and drift-correction for enhancing measurement accuracy; 3) obtaining initial atom positions; 4) indexing the atoms based on unit cell vectors; 4) quantifying the atom column positions with either 2D-gaussian single peak fitting or 5) multi-peak fitting routines for slightly overlapping atomic columns; 6) quantification of lattice distortion/strain within the crystal structures or at the defects/interfaces where the lattice periodicity is disrupted, and 7) some common methods to visualize and present the analysis. Furthermore, a simple self-developed free MATLAB app (EASY-STEM) with a graphical user interface (GUI) will be introduced that can help with the analysis of STEM images without the need for writing dedicated analysis code or software. The advanced data analysis methods presented here can be applied for the local quantification of defect relaxations, local structural distortions, local phase transformations, and non-centrosymmetry in a wide range of materials.

47 OTHER INSTRUMENTATION↗

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗

RELAP-7 Pressurizer Component Development Updates

RELAP-7 is a nuclear systems safety analysis code being developed at the Idaho National Laboratory (INL). RELAP-7 development began in 2011 to support the Risk Informed Safety Margins Characterization (RISMC) Pathway of the Light Water Reactor Sustainability (LWRS) program. The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical methods, and physical models in order to provide capabilities needed for the RISMC methodology and to support nuclear power safety analysis. The code is being developed based on Idaho National Laboratory’s modern scientific software development framework – MOOSE (the Multi-Physics Object-Oriented Simulation Environment). The initial development goal of the RELAP-7 approach focused primarily on the development of an implicit algorithm capable of strong (nonlinear) coupling of the dependent hydrodynamic variables contained in the 1-D/2-D flow models with the various 0-D system reactor components that compose various boiling water reactor (BWR) and pressurized water reactor nuclear power plants (NPPs). As part of the efforts to expand the capability for PWR simulation, an equilibrium single-region pressurizer model has been implemented in RELAP-7. The pressurizer component can simulate pressure and water level change through insurge, spray, and heating processes. Two simple tests – one for insurge process and another for outsurge process – have been reported to demonstrate and verify the functions of the pressurizer model. The typical single-phase PWR system model presented in the first RELAP-7 milestone report has been updated, as part of system level test for the new pressurizer model. The updated PWR system model with the pressurizer component can be used for more realistic transient simulations. The addition of the equilibrium single-region pressurizer model represents the first step of developing a suite of pressurizer models with different levels of fidelities. Non-equilibrium two-region two-zone model and non-equilibrium five-region three zone model will be implemented in the near future. Cross comparison and validation will be pursued when higher fidelity pressurizer models are completed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Co-Design Study Of Fusion Whole Device Modeling Using Code Coupling

Complex workflows consisting of multiple simulation and analysis codes running concurrently through in-memory coupling is becoming popular due to inherent advantages in online management of large-scale data, resilience, and the code development process. However, orchestrating such a multi-application workflow to efficiently utilize resources on a heterogeneous architecture is challenging. In this paper, we present our results with running the Fusion Whole Device Modeling benchmark workflow on Summit, a pre-exascale supercomputer at Oak Ridge National Laboratory. We explore various resource distribution and process placement mechanisms, including sharing compute nodes between processes from separate applications. We show that fine-grained process placement can have a significant impact towards efficient utilization of the compute power of a node on Summit, and conclude that sophisticated tools for performing co-design studies of multi-application workflows can play an important role towards efficient orchestration of such workflows.

Choi, Jong Youl↗

RELAP-7 Beta Release: Summary of Capabilities

RELAP-7 is a nuclear systems safety analysis code being developed at the Idaho National Laboratory (INL). Building upon the decades of software development at the INL, we began the development of RELAP-7 in 2011 to support the Risk Informed Safety Margins Characterization (RISMC) Pathway. As part of this development, the first lines of RELAP-7 code were committed to the software revision control repository on November 7th, 2011. The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical methods, and physical models in order to provide capabilities needed for the RISMC methodology and to support nuclear power safety analysis. RELAP-7 is built using the INL’s modern scientific software development framework, MOOSE (Multi-physics Object Oriented Simulation Environment). MOOSE provides improved numerical calculations (including higher-order integration in both space and time, yielding converged second-order accuracy). The RELAP-7 code structure is based on multiple physical component models such as pipes, junctions, pumps, etc. Each component can have options for different fluid models such as single- and two-phase flow. This component-based and physics-based software architecture allows RELAP-7 to adopt different physical models for different applications. A relatively new two-phase hydrodynamic model, termed the ''7-Equation model'' for two phasic pressures, velocities, energies, and volumetric fraction, is incorporated into RELAP-7 for liquid-gas (water-steam) flows. This new model allows second-order integration because it is well-posed, which will reduce the numerical error associated with traditional systems analysis codes. In this paper, we provide a RELAP-7 capability list describing analysis features, range of applicability, and reactor components that will be available for the December 15th, 2014 beta release of the software.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Requirements Description of DASSH-F

This report reviews the modeling and simulation capabilities of Argonne National Laboratory’s DASSH code that is used in present reactor analysis activities. These capabilities will be used to establish the set of verification tasks necessary to verify DASSH for use on commercial projects. A similar approach was taken for the PERSENT, REBUS and DIF3D software packages. The DASSH program is a thermal analysis code designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. DASSH is a follow-on development to the SE2-ANL software and SUPERENERGY-2 software that it is based upon. DASSH was designed to account for both neutron and gamma heating and is inherently connected to the GAMSOR part of the ARC suite of fast reactor analysis software. SE2-ANL is a developed piece of software from the 1980s while DASSH is a modern implementation with notable improvements in geometry handling. The most important upgrade of DASSH relative to SE2-ANL is that it can analyze multiple time points in a single run where SE2-ANL can only treat a single time point. This allows the user to understand the impact of and search the flow distribution for the entire operational period of a reactor design considering pressure drop, peak coolant and fuel temperatures, and thermal striping. DASSH has three input paths that have to be verified. The first input path builds the geometry and power distribution based upon the DIF3D model but ignores the gamma heating aspects of the problem. The second input path also builds the geometry from the DIF3D model but it takes the neutron and gamma heating distributions from GAMSOR. The third input path is to take the geometry and power distribution directly from user input (i.e. not coupled to DIF3D or GAMSOR). DASSH also has many built in correlations for material properties along with a user defined specification of the fuel, structure, and coolant properties. There are correlations for flow split, mixing, pressure drop, and heat transfer coefficients (subchannel rather than a direct methodology). In total, verification of DASSH will require an extensive testing to cover all possible user features of the software.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗