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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Shutdown-induced tensile stress in monolithic miniplates as a possible cause of plate pillowing at very high burnup

Post-irradiation examination of Reduced Enrichment for Research and Test Reactors (RERTR)-12 miniplates showed that in-reactor pillowing occurred in at least 4 plates, rendering performance of these plates unacceptable. To address in-reactor failures, efforts are underway to define the mechanisms responsible for in-reactor pillowing, and to suggest improvements to the fuel plate design and operational conditions. To achieve these objectives, the mechanical response of monolithic fuel to fission and thermally-induced stresses was modeled using a commercial finite element analysis code. Calculations of stresses and deformations in monolithic miniplates during irradiation and after the shutdown revealed that the tensile stress generated in the fuel increased from 2 MPa to 100 MPa at shutdown. The increase in tensile stress at shutdown possibly explains in-reactor pillowing of several RERTR-12 miniplates irradiated to the peak local burnup of up to 1.11x1022 fissions/cm3 . This paper presents the modeling approach and calculation results, and compares results with post-irradiation examinations and mechanical testing of irradiated fuel. The implications for the safe use of the monolithic fuel in research reactors are discussed, including the influence of fuel burnup and power on the magnitude of the shutdown-induced tensile stress.

rertr↗
DOE OSTI2015

Experimental study of forced convection heat transfer during upward and downward flow of helium at high pressure and high temperature

Fundamental high pressure/high temperature forced convection experiments have been conducted in support of the development of a Very High Temperature Reactor (VHTR) with a prismatic core. The experiments utilize a high temperature/high pressure gas flow test facility constructed for forced convection and natural circulation experiments. The test section has a single 16.8 mm ID flow channel in a 2.7 m long, 108 mm OD graphite column with four 2.3kW electric heater rods placed symmetrically around the flow channel. This experimental study presents the role of buoyancy forces in enhancing or reducing convection heat transfer for helium at high pressures up to 70 bar and high temperatures up to 873 degrees K. Wall temperatures have been compared among 10 cases covering the inlet Re numbers ranging from 500 to 3,000. Downward flows display higher and lower wall temperatures in the upstream and downstream regions, respectively, than the upward flow cases due to the influence of buoyancy forces. In the entrance region, convection heat transfer is reduced due to buoyancy leading to higher wall temperatures, while in the downstream region, buoyancyinduced mixing causes higher convection heat transfer and lower wall temperatures. However, their influences are reduced as the Reynolds number increases. This experimental study is of specific interest to VHTR design and validation of safety analysis codes.

Buoyancy↗
DOE OSTI2015

SOFAST: Sandia Optical Fringe Analysis Slope Tool

SOFAST is used to characterize the surface slope of reflective mirrors for solar applications. SOFAST uses a large monitor or projections screen to display fringe patterns, and a machine vision camera to image the reflection of these patterns in the subject mirror. From these images, a detailed map of surface normals can be generated and compared to design or fitted mirror shapes. SOFAST uses standard Fringe Reflection (Deflectometry) approaches to measure the mirror surface normals. SOFAST uses an extrinsic analysis of key points on the facet to locate the camera and monitor relative to the facet coordinate system. It then refines this position based on the measured surface slope and integrated shape of the mirror facet. The facet is placed into a reference frame such that key points on the facet match the design facet in orientation and position. This is key to evaluating a facet as suitable for a specific solar application. SOFAST reports the measurements of the facet as detailed surface normal location in a format suitable for ray tracing optical analysis codes. SOFAST also reports summary information as to the facet fitted shape (monomial) and error parameters. Useful plots of the error distribution are also presented.

Andraka, Charles E.↗
DOE OSTI2016

IMPLEMENTATION AND VALIDATION OF A FULLY IMPLICIT ACCUMULATOR MODEL IN RELAP-7

This paper presents the implementation and validation of an accumulator model in RELAP-7 under the framework of preconditioned Jacobian free Newton Krylov (JFNK) method, based on the similar model used in RELAP5. RELAP-7 is a new nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). RELAP-7 is a fully implicit system code. The JFNK and preconditioning methods used in RELAP-7 is briefly discussed. The slightly modified accumulator model is summarized for completeness. The implemented model was validated with LOFT L3-1 test and benchmarked with RELAP5 results. RELAP-7 and RELAP5 had almost identical results for the accumulator gas pressure and water level, although there were some minor difference in other parameters such as accumulator gas temperature and tank wall temperature. One advantage of the JFNK method is its easiness to maintain and modify models due to fully separation of numerical methods from physical models. It would be straightforward to extend the current RELAP-7 accumulator model to simulate the advanced accumulator design.

Accumulator↗
DOE OSTI2015

A STRONGLY COUPLED REACTOR CORE ISOLATION COOLING SYSTEM MODEL FOR EXTENDED STATION BLACK-OUT ANALYSES

The reactor core isolation cooling (RCIC) system in a boiling water reactor (BWR) provides makeup cooling water to the reactor pressure vessel (RPV) when the main steam lines are isolated and the normal supply of water to the reactor vessel is lost. The RCIC system operates independently of AC power, service air, or external cooling water systems. The only required external energy source is from the battery to maintain the logic circuits to control the opening and/or closure of valves in the RCIC systems in order to control the RPV water level by shutting down the RCIC pump to avoid overfilling the RPV and flooding the steam line to the RCIC turbine. It is generally considered in almost all the existing station black-out accidents (SBO) analyses that loss of the DC power would result in overfilling the steam line and allowing liquid water to flow into the RCIC turbine, where it is assumed that the turbine would then be disabled. This behavior, however, was not observed in the Fukushima Daiichi accidents, where the Unit 2 RCIC functioned without DC power for nearly three days. Therefore, more detailed mechanistic models for RCIC system components are needed to understand the extended SBO for BWRs. As part of the effort to develop the next generation reactor system safety analysis code RELAP-7, we have developed a strongly coupled RCIC system model, which consists of a turbine model, a pump model, a check valve model, a wet well model, and their coupling models. Unlike the traditional SBO simulations where mass flow rates are typically given in the input file through time dependent functions, the real mass flow rates through the turbine and the pump loops in our model are dynamically calculated according to conservation laws and turbine/pump operation curves. A simplified SBO demonstration RELAP-7 model with this RCIC model has been successfully developed. The demonstration model includes the major components for the primary system of a BWR, as well as the safety system components such as the safety relief valve (SRV), the RCIC system, the wet well, and the dry well. The results show reasonable system behaviors while exhibiting rich dynamics such as variable flow rates through RCIC turbine and pump during the SBO transient. The model has the potential to resolve the Fukushima RCIC mystery after adding the off-design two-phase turbine operation model and other additional improvements.

RCIC↗
DOE OSTI2017

SOLVING CONJUGATE HEAT TRANSFER PROBLEM IN TWO-PHASE FLOW USING NEWTON-KRYLOV METHOD

In this paper, we successfully demonstrated an application of the Jacobian-free Newton Krylov (JFNK) method in solving the conjugated two-phase channel flow and heat conduction problem, in a fully coupled and fully implicit manner. Such a problem represents the most fundamental piece of safety analysis codes, for example the heat transfer between fuel rods and reactor coolant during postulated accidental scenarios. The two-phase flow is modeled with the widely accepted one-dimensional two-fluid two-phase flow model, and the heat conduction is modeled with two-dimensional heat conduction model. The coupling between these two physics is realized through the wall heat transfer term in fluid equations and heat flux boundary conditions of heat conduction model. For spatial discretization, both the first-order and second- order upwind schemes are considered for flow equations, while the standard central differencing method is used for the heat conduction model. For time integration, both the first-order backward Euler (BDF1) and second-order BDF2 schemes are considered. Several treatments were also introduced to enhance the robustness of the JFNK method in solving such a problem that is discontinuous in nature because of the discrete closure correlations. These treatments include the re-introduction of a transition region between the slug and annular mist flow regimes, and a simplification in the interfacial drag term. Numerical results have been extensively validated with experimental data in both vertical pipe and rod bundle geometries with a wide range of thermal-hydraulics conditions, including pressure, inlet subcooling, wall heat flux, outlet void fraction, and etc. Further benchmark with the RELAP5-3D code confirms the correct implementation of wall heat transfer models.

42 ENGINEERING↗
DOE OSTI2016

An Equation of State for Liquid Pb83Li17

Thermal hydraulic and accident analysis codes such as RELAP5-3D and MELCOR for Fusion rely on an equation of state to specify all the thermodynamic properties of fusion-relevant working fluids such as PbLi. The existing liquid metal fluid properties in both RELAP5-3D and MELCOR are based on a five parameter "soft sphere" equation of state for which parameter sets that approximately reproduce experiment data were available in the literature. The PbLi parameters were based on a linear, mass-weighted average of the parameters for pure Pb and pure Li, which does not always result in adequate agreement between the equation and experimental data for PbLi, much of which was not available at the time the original equations were implemented. In order to address this shortcoming, we describe here a modification and non-linear least squares fitting of the parameters to closely match available experiment data for the density, sound speed, and specific heat capacity of liquid PbLi over a range of temperatures. The new equation has been implemented in both RELAP5-3D and MELCOR, and we make comparisons with both codes to experiment data to verify the implementation. Transport properties such as the thermal conductivity and viscosity, though not derived from the equation of state, have also been updated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Flexible and Scalable Data-Acquisition Using the artdaq Toolkit

The Real-Time Systems Engineering Department of the Scientific Computing Division at Fermilab is developing a flexible, scalable, and powerful data-acquisition (DAQ) toolkit which serves the needs of experiments from bench-top hardware tests to large high-energy physics experiments. The toolkit provides data transport and event building capabilities with the option for experimenters to inject art analysis code at key points in the DAQ for filtering or monitoring. The toolkit also provides configuration management, run control, and low-level hardware communication utilities. Firmware blocks for several commercial data acquisition boards are provided, allowing experimenters to approach the DAQ from a high level. A fully-functional DAQ "solution" of the toolkit is provided in otsdaq, sacrificing some flexibility in favor of being ready-to-use. artdaq is being used for several current and upcoming experiments, and will continue to be refined and expanded for use in the next generation of neutrino and muon experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

RISK MONITORING CAPABILITES FROM DYNAMIC PRA DATA

Risk Monitors (RMs) are software tools that are employed in nuclear power plants to monitor the safety properties of the plant itself during operation. These safety properties are generated using classical PRA methods which are typically based on static Boolean structures such as Fault-Trees (FTs) and Event-Trees (ETs). The outcomes of these methods are Core Damage Frequency (CDF) and risk importance values of plant components (e.g., valves and pumps). Examples of Risk Importance Measures (RIMs) commonly used are: Risk Achievement Worth (RAW), Risk Reduction Worth (RRW), Birnbaum (B) and Fussell-Vesely (FV) [1,2]. In a RM these values are updated in real time once plant components status is updated (e.g., under maintenance, off-line, testing). In addition, RMs can be used to optimize maintenance schedules so those safety features of the plant are below plant standard values. This paper presents an approach to create temporal profiles of RIMs (i.e., RM capabilities) given a temporal change of the plant configuration. Instead of employing data generated by ETs and FTs, the proposed approach will employ data generated by Dynamic PRA methods. In contrast to classical PRA methods, Dynamic PRA methods couple stochastic methods with safety analysis codes to determine risk associate to complex systems such as nuclear plants. Compared to classical PRA methods, they can evaluate with higher resolution the safety impact of timing and sequencing of events on the accident progression. We will describe how classical RIMs can be generated by Dynamic PRA data and then we will show how the temporal profiles of the RIMs can be generated given a temporal profile of the system configuration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Operation Data from the UW Natural Circulation FLiBe Flow Loop

Operating measurements from three multi-week test campaigns of the natural circulation FLiBe loop. This system is investigating thermal hydraulic behavior of the molten salt BeF_2 -LiF (33 - 67 mol %). The system behaves in a stable fashion, but shows unusual local transients, such as flow break-down in the riser and thermal jumps at the cooler exit. Heat transfer shows promising trends that FLiBe will behave as a normal heat transfer fluid, as long as salt purity can be maintained. The most recent test shows heat transfer degradation that is likely a result of oxides and impurities. The data archive contains as-built dimensions, Matlab analysis codes, and the raw data files.

42 ENGINEERING↗
DOE OSTI2019

Parfume (particle Fuel Model)

PARFUME is an advanced gas-cooled reactor fuel performance modeling and analysis code. It has been developed as an integrated mechanistic code that evaluates the thermal, mechanical, and physio-chemical behavior of fuel particles during irradiation and postulated accident conditions to determine the failure probability of a population of fuel particles given the particle-to-particle statistical variations in physical dimensions and material properties that arise from the fuel fabrication process, accounting for all viable mechanisms that can lead to particle failure. The code also determines the diffusion of fission products from the fuel through the particle coating layers, and through the fuel matrix to the coolant boundary. The subsequent release of fission products is calculated at the fuel element level (release of fission products from the compact or pebble). PARFUME calculates the release fraction as the ratio of the number of atoms released from the fuel element to the amount.

Gougar, Hans D.↗
DOE OSTI2019

NOBLE METAL MASS TRANSPORT MODEL FOR MOLTEN SALT REACTOR ANALYSIS IN VERA-CS

Molten Salt Reactor designs require robust multiphysics modelling capabilities to further the development of this next generation technology. Due to the unique liquid fuel design, traditional reactor physics analysis codes do not capture all of the important physics of these reactors. This work builds upon a general species transport model implemented into the multiphysics core simulator suite VERA-CS, by adding source and sink models for non-soluble fission products, specifically the noble metals. Noble metal fission products are born in the fuel-salt, but do not form stable fluorides, and instead plate out on various surfaces in the reactor fuel loop. A concentration gradient diffusion-convection based mass transfer model is implemented into the general species transport model within CTF-the subchannel thermal hydraulics code in VERA-CS-and verified. Two noble metal decay chains are analyzed, and their steady state bulk and surface concentrations are calculated for a simple flow loop in CTF that is representative of the Oak Ridge National Laboratory Molten Salt Reactor Experiment (MSRE). Non-soluble species 99Mo to 99Tc captures wall plate out of decaying species, and non-soluble 135Te to soluble 135I and gaseous 135Xe captures wall loss of decaying species. Results are compared with MSRE data and the possible multiphysics effects of this phenomena are identified. Future work to enhance the robustness of the physical model and couple it with neutronics to evaluate these effects of Noble Metal mass transport is discussed.

Walker, Samuel↗
DOE OSTI2021

Analysis of PNAR Spent Fuel Safeguards Measurements using the ORIGEN Data Analysis Module

This paper summarizes the analysis of the Passive Neutron Albedo Reactivity (PNAR) measurements using the ORIGEN data analysis module for 23 boiling water reactor spent fuel assemblies that were performed in Finland under an international collaboration on spent fuel safeguards verification methods. PNAR is part of the proposed integrated nondestructive system to be used for safeguards verifications at the planned Finnish encapsulation facility. Besides measuring passive neutron and gamma emission rates from an assembly like a Fork detector, PNAR also measures the PNAR ratio, which is expected to correlate with the fissile content in the assembly. The emission rates and PNAR ratio can be used to verify the operator declarations and the fissile content of an assembly, respectively. The analysis was performed using the ORIGEN Data Analysis Module, which was originally developed for predicting Fork detector neutron and gamma signals for spent fuel measurements and has been integrated into the Integrated Review and Analysis Package developed by Euratom and the IAEA. The Module includes the ORIGEN burnup analysis code and integrates detector response functions pre-generated using MCNP to predict detector signals in several seconds per assembly. In this study, new response functions specific to PNAR measurements were generated for ORIGEN Module. The study also analyzes impacts of using detailed fuel design and operation information vs. standard safeguards information on results calculated by ORIGEN Module. Using detailed information reduced the standard deviation of relative differences between calculated and measured neutron count rates among the 23 assemblies from ~10% to ~4%. The results obtained using standard safeguards information for these PNAR measurements were similar to those obtained for the Fork detector. A clear trend was found between the calculated net neutron multiplications and the measured PNAR ratios of the 23 assemblies. This paper describes how ORIGEN Module calculates the expected PNAR neutron and gamma signals and PNAR ratio and how they compare with corresponding measured values.

Ilas, Germina↗
DOE OSTI2022

Design of a Test Article and Test Facility to Investigate Flow-Induced Vibration of a Helical Coil Steam Generator

NuScale Power LLC (NuScale) is developing a small modular reactor (SMR) with unique design features and arrangements that require testing to provide quality data for safety analysis codes and design validation, as well as for technology maturation. A multi-phase testing program is underway at NuScale to demonstrate that its unique helical coil steam generator design is not susceptible to flow-induced vibration during operation. The phases of the testing program include modal testing of individual tubes during test article fabrication and modal and flow testing after test article completion. NuScale has contracted SIET S.p.A. in Piacenza, Italy for the design and fabrication of the test article and test facility. In-fabrication testing was performed at the site of test article build in Cremosano, Italy. Modal testing and flow testing of the completed test article will be performed at SIET’s facilities in Piacenza. This paper discusses the design of the test article and test facility, including supporting systems, required to generate the quality data required to support NuScale’s comprehensive vibration assessment program as part of the overall licensing effort for NuScale’s first-of-a-kind small modular light water reactor.

Flow-induced vibration, helical coil, steam genera↗
DOE OSTI2021

Radial Compressor Design and Off-Design for Trans-critical CO2 Operating Conditions

Presentation for the The 7th International Supercritical CO2 Power Cycles Symposium, February 21 – 24, 2022, San Antonio, Texas. A multi-staged intercooled centrifugal compressor design is performed for an indirect-fired sCO2 cycle and the design point data are used to generate the performance maps using a 2D mean-line analysis code. Fluid and turbomachinery similarity laws are applied to find an appropriate surrogate ideal gas fluid to be used in the design and performance map generation codes. The performance maps are then scaled back for the sCO2 fluid case using the affinity laws.

Uysal, Selcuk Can↗
DOE OSTI2021

Radial Compressor Design and Off-Design for Trans-critical CO2 Operating Conditions

A multi-staged intercooled centrifugal compressor design is performed for an indirect-fired sCO2 cycle and the design point data are used to generate the performance maps using a 2D mean-line analysis code. Fluid and turbomachinery similarity laws are applied to find an appropriate surrogate ideal gas fluid to be used in the design and performance map generation codes. The performance maps are then scaled back for the sCO2 fluid case using the affinity laws. The developed methodology is validated using the experimental data for a compressor operating at similar inflow conditions.

Uysal, Selcuk Can↗
DOE OSTI2022

Assessment of Thermal Stratification in Versatile Test Reactor Transients

The Versatile Test Reactor (VTR) is a fast-spectrum test reactor currently being developed in the United States under the direction of the US Department of Energy. The conceptual design of the 300 MWth pool-type sodium-cooled fast reactor (SFR) has been led by the US National Laboratories in collaboration with General Electric-Hitachi and Bechtel National Inc. Safety performance analysis for the VTR conceptual design is being performed with the systems thermal-hydraulics (SYSTH) module of the SAS4A/SASSYS-1 liquid-metal reactor safety analysis code system. Since the current model of the VTR employs a simple perfect mixing model for large plena like the hot pool, it is not able to predict temperature variations that may develop during the transient. Prior work simulating the response of SFRs to postulated events like the Protected Station Blackout (PSBO) has shown that the phenomenon of thermal stratification, where stable thermal layers accumulate in the hot pool, may delay the transition to natural circulation and thus impact the predicted transient progression. Thus, an effort has begun to model this transient by integrating a Computational Fluid Dynamics (CFD) model of the hot pool into the SAS4A/SASSYS-1 model of the Primary Heat Transport System during the simulation of the PSBO event. A three-dimensional Volume-Of-Fluid CFD model of the VTR hot pool has been developed for the co-simulation of SAS4A/SASSYS-1 with CFD. In this work, the standalone SAS4A/SASSYS-1 calculation and the standalone CFD calculation based on the SAS4A/SASSYS-1 calculation result were produced on the Idaho National Laboratory High Performance Computing cluster, SAWTOOTH. At this time, only the standalone CFD and SAS4A/SASYS-1 simulation results are provided.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2013

Calibration of a fuel relocation model in BISON

We demonstrate parameter calibration in the context of the BISON nuclear fuels performance analysis code. Specifically, we present the calibration of a parameter governing fuel relocation: the power level at which the relocation model is activated. This relocation activation parameter is a critical value in obtaining reasonable comparison with fuel centerline temperature measurements. It also is the subject of some debate in terms of the optimal values. We show that the optimal value does vary across the calibration to individual rods. We also demonstrate an aggregated calibration, where we calibrate to observations from six rods. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗