Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fuel performance 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 1,099 records · Page 61

Ignition and combustion characteristics of metallized propellants, phase 2

Secondary atomization and ignition characteristics of aluminum/hydrocarbon gel propellants were investigated. Models of gel droplet shell formation were applied to aluminum/liquid hydrocarbon propellants to examine the effects of solid loading and ultimate particle size on the minimum droplet diameter permitting secondary atomization. A one-dimensional model of a gel-fueled rocket combustion chamber was developed. A model for radiant heat transfer from hot aluminum oxide particles to the chamber walls is included. A two-dimensional, two-phase nozzle code was used to estimate nozzle two-phase losses and overall engine performance.

Mueller, D. C.↗

Structural Evaluation of a Space Shuttle Main Engine (SSME) High Pressure Fuel Turbopump Turbine Blade

Thermal and structural finite-element analyses were performed on the first high pressure fuel turbopump turbine blade of the space shuttle main engine (SSME). A two-dimensional (2-D) finite-element model of the blade and firtree disk attachment was analyzed using the general purpose MARC (finite-element) code. The loading history applied is a typical test stand engine cycle mission, which consists of a startup condition with two thermal spikes, a steady state and a shutdown transient. The blade material is a directionally solidified (DS) Mar-M 246 alloy, the blade rotor is forged with waspalloy material. Thermal responses under steady-state and transient conditions were calculated. The stresses and strains under the influence of mechanical and thermal loadings were also determined. The critical regions that exhibited high stresses and severe localized plastic deformation were the blade-rotor gaps.

Abdul-Aziz, Ali↗

Impact of Reactor Operating Parameters on Cask Reactivity in BWR Burnup Credit

This paper discusses the effect of reactor operating parameters used in fuel depletion calculations onspent fuel cask reactivity, with relevance for boiling-water reactor (BWR) burnup credit (BUC)applications. Assessments that used generic BWR fuel assembly and spent fuel cask configurationsare presented. The considered operating parameters, which were independently varied in thedepletion simulations for the assembly, included fuel temperature, bypass water density, specificpower, and operating history. Different operating history scenarios were considered for the assemblydepletion to determine the effect of relative power distribution during the irradiation cycles, as wellas the downtime between cycles. Depletion, decay, and criticality simulations were performed usingcomputer codes and associated nuclear data within the SCALE code system. Results quantifying thedependence of cask reactivity on the assembly depletion parameters are presented herein.

Ilas, Germina↗

MONTE CARLO ANALYSES OF THE YALINA THERMAL FACILITY WITH SERPENT STEREOLITHOGRAPHY GEOMETRY MODEL

This paper analyzes the YALINA Thermal subcritical assembly of Belarus using two different Monte Carlo transport programs, SERPENT and MCNP. The MCNP model is based on combinatorial geometry and universes hierarchy, while the SERPENT model is based on Stereolithography geometry. The latter consists of unstructured triangulated surfaces defined by the normal and vertices. This geometry format is used by 3D printers and it has been created by: the CUBIT software, MATLAB scripts, and C coding. All the Monte Carlo simulations have been performed using the ENDF/B-VII.0 nuclear data library. Both MCNP and SERPENT share the same geometry specifications, which describe the facility details without using any material homogenization. Three different configurations have been studied with different number of fuel rods. The three fuel configurations use 216, 245, or 280 fuel rods, respectively. The numerical simulations show that the agreement between SERPENT and MCNP results is within few tens of pcms.

CUBIT↗

Program ELM: A tool for rapid thermal-hydraulic analysis of solid-core nuclear rocket fuel elements

This report reviews the state of the art of thermal-hydraulic analysis codes and presents a new code, Program ELM, for analysis of fuel elements. ELM is a concise computational tool for modeling the steady-state thermal-hydraulics of propellant flow through fuel element coolant channels in a nuclear thermal rocket reactor with axial coolant passages. The program was developed as a tool to swiftly evaluate various heat transfer coefficient and friction factor correlations generated for turbulent pipe flow with heat addition which have been used in previous programs. Thus, a consistent comparison of these correlations was performed, as well as a comparison with data from the NRX reactor experiments from the Nuclear Engine for Rocket Vehicle Applications (NERVA) project. This report describes the ELM Program algorithm, input/output, and validation efforts and provides a listing of the code.

Walton, James T.↗

Multi-fluid transport code modeling of time-dependent recycling in ELMy H-mode

Simulations of a high-confinement-mode (H-mode) tokamak discharge with infrequent giant type-I ELMs are performed by the multi-fluid, multi-species, two-dimensional transport code UEDGE-MB, which incorporates the Macro-Blob approach for intermittent non-diffusive transport due to filamentary coherent structures observed during the Edge Localized Modes (ELMs) and simple time-dependent multi-parametric models for cross-field plasma transport coefficients and working gas inventory in material surfaces. Temporal evolutions of pedestal plasma profiles, divertor recycling, and wall inventory in a sequence of ELMs are studied and compared to the experimental time-dependent data. Short- and long-time-scale variations of the pedestal and divertor plasmas where the ELM is described as a sequence of macro-blobs are discussed. It is shown that the ELM recovery includes the phase of relatively dense and cold post-ELM divertor plasma evolving on a several ms scale, which is set by the transport properties of H-mode barrier. The global gas balance in the discharge is also analyzed. The calculated rates of working gas deposition during each ELM and wall outgassing between ELMs are compared to the ELM particle losses from the pedestal and neutral-beam-injection fueling rate, correspondingly. A sensitivity study of the pedestal and divertor plasmas to model assumptions for gas deposition and release on material surfaces is presented. The performed simulations show that the dynamics of pedestal particle inventory is dominated by the transient intense gas deposition into the wall during each ELM followed by continuous gas release between ELMs at roughly a constant rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Preliminary Design of the TEX-MOX Criticality Experimental Program using Optimization Algorithms Implemented in the IRSN PROMETHEE Tool

In the frame of its missions, IRSN is widely involved in design, evaluation and analyses of critical experiments in order to contribute to computer codes and nuclear data validation. The work described in this paper was carried out in the framework of the PRINCESS IRSN research project and is part of the IRSN/LLNL collaboration on the TEX program in the frame of the NCSP. The need for epithermal and intermediate energy range critical benchmarks is an established international criticality safety data need, underpinned by the subsequent use of data assimilation methods like GLLSM to reduce Nuclear Data uncertainties. The goal of the TEX program is to address this need by executing experiments with NCSP fissile assets, which can be used to create critical assemblies that span a wide range of fission energy spectra, from the thermal energy range (below 0.625 eV), through the intermediate energy range (0.625 eV to 100 keV), to the fast one (above 100 keV). The TEX program involves different materials of interest, such as TEX-Ta experiments using plutonium ZPPR plates proposed by LLNL. This paper presents an overview of the studies performed at IRSN using numerical optimization algorithms for the design of critical experiments involving Mixed U and Pu fuel plates, referred to as TEX-MOX program. This program ultimately aims to perform critical experiments to validate the MOX fuel configurations encountered in the MOX fuel fabrication process. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Monte Carlo estimation of an effect of uncertainties in initial data on solving the transport equation by means of the MCU code

A new version of the tally module of the MCU software package is developed in which the approach for taking directly into account the uncertainty in initial data is implemented that is recommended by the international standard on estimating the uncertainty in results of measuring (ISO 13005). The new module makes it possible to evaluate the effect of uncertainty in initial data (caused by technological tolerances in fabrication of structural members of the core) on neutronic characteristics of the reactor. The developed software is adapted to parallel computing with the use of multiprocessor computers, which significantly reduces the computation time: the parallelization coefficient is almost equal to 1. Testing is performed by examples of solving the problem on criticality for the Godiva benchmark experiment and also for the infinite lattice of fuel assemblies of the VVER-440, VVER-1000, and VVER-1200. The results of calculations of the uncertainty in neutronic characteristics (effective multiplication factor, fission reaction rate), which is caused by uncertainties in initial data due to technological tolerances, are compared (in the first case) to the published results obtained using the precision MCNP5 code and (in the second case) to those obtained by means of the RADAR engineering program. A good agreement of results is achieved for all cases.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

TRACE Code Validation of BWR SVEA Fuel Assembly Spray Cooling

Emergency spray cooling injection is a safety feature implemented in many boiling water reactor designs to reflood a reactor during an accident. Significant experimental work has qualified the efficacy of spray cooling, and ongoing computational modeling efforts strive to more accurately portray the phenomena involved. This work examines the physical phenomena pertaining to emergency spray cooling injection over SVEA-type fuel assemblies. BWR emergency core cooling system (ECCS) designs commonly incorporate a Low-Pressure Core Spray (LPCS) system that introduces coolant through spray nozzles directed to the top of the core, and are located in the upper plenum region. This design is effective for LOCA response, particularly during the initial refill and reflood stages of accident response where the goal is limiting the peak cladding temperature rise in the core. The U.S. NRC thermal-hydraulics code TRACE version 5.0 Patch 4 has been chosen to simulate the separate-effect tests performed by ASEA-ATOM. The computational model was evaluated by performing forward uncertainty quantification using Dakota as the analysis tool and code driver on 31 identified parameters from the input model and TRACE physical constitutive models. This paper discusses the sensitivity of the investigated parameters with respect to spray cooling simulation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Three-dimensional neutronics optimization of helium-cooled blanket for multi-functional experimental fusion-fission hybrid reactor (FDS-MFX)

Three-dimensional neutronics optimization calculations were performed to analyse the parameters of Tritium Breeding Ratio (TBR) and maximum average Power Density (PDmax) in a helium-cooled multi-functional experimental fusion-fission hybrid reactor named FDS (Fusion-Driven hybrid System)-MFX (Multi-Functional experimental) blanket. Three-stage tests will be carried out successively, in which the tritium breeding blanket, uranium-fueled blanket and spent-fuel-fueled blanket will be utilized respectively. In this contribution, the most significant and main goal of the FDS-MFX blanket is to achieve the PDmax of about 100 MW/m3 with self-sustaining tritium (TBR {>=} 1.05) based on the second-stage test with uranium-fueled blanket to check and validate the demonstrator reactor blanket relevant technologies based on the viable fusion and fission technologies. Four different enriched uranium materials were taken into account to evaluate PDmax in subcritical blanket: (i) natural uranium, (ii) 3.2% enriched uranium, (iii) 19.75% enriched uranium, and (iv) 64.4% enriched uranium carbide. These calculations and analyses were performed using a home-developed code VisualBUS and Hybrid Evaluated Nuclear Data Library (HENDL). The results showed that the performance of the blanket loaded with 64.4% enriched uranium was the most attractive and it could be promising to effectively obtain tritium self-sufficiency (TBR-1.05) and a high maximum average power density ({approx}100 MW/m{sup 3}) when the blanket was loaded with the mass of {sup 235}U about 1 ton. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

CFD Verification and Validation of Wire-Wrapped Pin Assemblies

TerraPower is participating in a cooperative project among industry, the national lab system, and academia to perform Verification and Validation (V and V) of Computational Fluid Dynamics (CFD) methods for predicting the flow and heat transfer within liquid-metal-cooled nuclear fuel assemblies with wire-wrapped fuel pins. This project, consisting of both experimental and numerical components, uses surrogate fluids and electrically heated fuel pins to substitute for the liquid metal coolant and the nuclear heat source. The experiments include both unheated assemblies to measure the flow velocities with Laser Doppler Velocimetry (LDV) and Particle Image Velocimetry (PIV), as well as heated assemblies to measure the temperatures of the fluid and the simulated fuel pins with thermocouples. Both sets of experiments also measure the pressure at multiple locations in the system. The numerical component of this work involves high-fidelity Large-Eddy Simulation (LES) modeling and industry-standard Reynolds-Averaged Navier-Stokes (RANS) modeling of the experiments. Furthermore, the experiments and simulations will be conducted on two duct geometries: the first is a non-deformed hexagonal duct that has flat internal duct faces that represents a beginning of life geometry. The other is a deformed duct where the internal hexagonal flow area enlarges then decreases along the axial height to represent the changes that occur to the duct geometry from pressure and irradiation effects in the core. This paper describes the methods and results of the industrial, RANS-based CFD code verification, solution verification, and solution validation. The V and V simulations are being performed for the following experiment test geometries: - Isothermal non-deformed bundle; - Heated non-deformed bundle; - Isothermal deformed bundle; - Heated deformed bundle. These flows have been computed at first with the pre-test design specifications, without the aid of any experimental results. Later, post-test calculations will be performed incorporating the actual test conditions and input from some of the experimental results. These RANS simulations of helically wire-wrapped fuel assemblies employ meshes of bare pins without the wire-wrap in the bundle. The effect of wire-wrap on the flow is accounted for by introducing a momentum source into the governing fluid equations. The momentum source is only applied to cells corresponding to the location of the wire-wrap, and its components in each cell are based on the local flow field. This methodology simplifies the geometry and mesh generation and reduces the cell count, while still capturing the effects of the wire-wrap. The momentum source approach is used for all of the assembly simulations performed in this study. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Delayed Gamma-Ray Spectroscopy for Non-Destructive Assay of Nuclear Materials

Modeling capabilities were added to an existing framework and codes were adapted as needed for analyzing experiments and assessing application-specific assay concepts including simulation of measurements over many short irradiation/spectroscopy cycles. The code package was benchmarked against the data collected at the IAC for small targets and assembly-scale data collected at LANL. A study of delayed gamma-ray spectroscopy for nuclear safeguards was performed for a variety of assemblies in the extensive NGSI spent fuel library. The modeling results indicate that delayed gamma-ray responses can be collected from spent fuel assemblies with statistical quality sufficient for analyzing their isotopic composition using a 10 11 n/s neutron generator and COTS detector instrumentation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development and validation of two-phase flow models and critical heat flux prediction for the highly-scalable CFD code NEK-2P

A project is underway to develop, verify and validate an advanced two-phase flow modeling capability for the highly-scalable high-performance CFD code NEK5000[1]. The goal of the project is to develop a new two-phase version of the NEK5000code, named NEK-2P, to simulate the two-phase flow and heat transfer phenomena that occur in a Boiling Water Reactor (BWR) fuel bundle under various operating conditions. The NEK-2P two-phase flow models follow the approach used for the Extended Boiling Framework [2-3] previously developed at Argonne, but include more fundamental physical models of boiling phenomena and advanced numerical algorithms for improved computational accuracy, robustness, and computational speed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Performance Benefits for Wave Rotor-Topped Gas Turbine Engines

The benefits of wave rotor-topping in turboshaft engines, subsonic high-bypass turbofan engines, auxiliary power units, and ground power units are evaluated. The thermodynamic cycle performance is modeled using a one-dimensional steady-state code; wave rotor performance is modeled using one-dimensional design/analysis codes. Design and off-design engine performance is calculated for baseline engines and wave rotor-topped engines, where the wave rotor acts as a high pressure spool. The wave rotor-enhanced engines are shown to have benefits in specific power and specific fuel flow over the baseline engines without increasing turbine inlet temperature. The off-design steady-state behavior of a wave rotor-topped engine is shown to be similar to a conventional engine. Mission studies are performed to quantify aircraft performance benefits for various wave rotor cycle and weight parameters. Gas turbine engine cycles most likely to benefit from wave rotor-topping are identified. Issues of practical integration and the corresponding technical challenges with various engine types are discussed.

Jones, Scott M.↗

Statistical uncertainty analysis applied to the DRAGONv4 code lattice calculations and based on JENDL-4 covariance data

In this paper, multi-group microscopic cross-section uncertainty is propagated through the DRAGON (Version 4) lattice code, in order to perform uncertainty analysis on k{infinity} and 2-group homogenized macroscopic cross-sections predictions. A statistical methodology is employed for such purposes, where cross-sections of certain isotopes of various elements belonging to the 172 groups DRAGLIB library format, are considered as normal random variables. This library is based on JENDL-4 data, because JENDL-4 contains the largest amount of isotopic covariance matrixes among the different major nuclear data libraries. The aim is to propagate multi-group nuclide uncertainty by running the DRAGONv4 code 500 times, and to assess the output uncertainty of a test case corresponding to a 17 x 17 PWR fuel assembly segment without poison. The chosen sampling strategy for the current study is Latin Hypercube Sampling (LHS). The quasi-random LHS allows a much better coverage of the input uncertainties than simple random sampling (SRS) because it densely stratifies across the range of each input probability distribution. Output uncertainty assessment is based on the tolerance limits concept, where the sample formed by the code calculations infers to cover 95% of the output population with at least a 95% of confidence. This analysis is the first attempt to propagate parameter uncertainties of modern multi-group libraries, which are used to feed advanced lattice codes that perform state of the art resonant self-shielding calculations such as DRAGONv4. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A New Safeguards Technique for Spent Fuel Verification in Dry Storage by Using He-4 Detectors - 18262

It is demanded to find key issues and critical paths for verifying and monitoring a spent nuclear fuel cask or storage with a purpose of non-proliferation of nuclear weapon. However, unfortunately there are no reliable methods to verify the containment in a spent nuclear fuel cask or storage. Therefore, research and development (R and D) has been continuously required to develop new and alternative technologies for searching and verifying undeclared nuclear materials, activities, and facilities. Therefore, the present authors have developed a new safeguards technique of passively scanning a spent fuel cask or storage using pressurized He-4 scintillation neutron detectors. A new approach is to quickly verify content in a spent nuclear fuel cask or presence of all special nuclear materials in unknown packages by comparing experimental and computational results. For that, the present authors have developed a new algorithm and user-interface coupled ORIGEN-ARP and MCNP codes, and the prototype safeguards apparatus employing the array of He-4 fast neutron detectors has been designed, built, and tested. The feasibility test of this new apparatus has been performed with Cf-252 neutron sources and the laboratory scale model cask, minimized by a factor of ten compared to the actual dimension of a KN-21 cask. This technique has several advantages over current methods based on neutron and gamma detection; (1) incident neutron energy information is preserved, (2) disturbance of gamma-ray is excellently rejected (discriminated), and (3) much higher neutron counts and operational cost and space-effectiveness are expected since no moderator is required for fast neutrons. The final goal of this study will (1) conduct a prototype instrument and employ it to monitor a nuclear spent fuel cask or storage at a commercial nuclear power reactor facilities, (2) use spectral-specific measurement data to show proof of concept for measurable data successfully converted to quantifiable signatures, such as diversion/removal of fuel, and maintenance of fuel and cask integrity of fuel, and (3) develop a cask imaging capability. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Multi-Objective Optimization of Uranium Target Assembly–3: A Comparison of Genetic and Traditional Methods

Commonly produced as a byproduct of uranium fission, 99 Mo is a key medical isotope that is in high demand in the United States. An international goal is to switch from medical isotope production technologies that require highly enriched uranium to medical isotope production technologies that require only low-enriched uranium. Niowave Inc. is contributing to this goal by developing an accelerator-driven subcritical assembly called the Uranium Target Assembly (UTA). This work compares the performance of Dakota’s Multi-Objective Genetic Algorithm (MOGA) against traditional sensitivity analysis in the neutronic optimization of the UTA-3 system. The design objectives are k-eigenvalue (k eff ) and natural uranium fission power, which are directly correlated with the amount of 99 Mo produced. Dakota:MOGA did not perform as well as human engineering ingenuity in optimization studies with high numbers of input parameters, such as fuel rod type selection and fuel rod placement. However, Dakota:MOGA did outperform traditional sensitivity analysis in optimization studies with fewer than 20 parameters and revealed the degree to which each parameter influences the optimal design space for k eff and natural uranium fission power (to a lesser extent). As the design model became more complex in the final stage of design, the computational resources required to calculate the design objective values in the Monte Carlo N-Particle transport code from selected input parameter combinations limited Dakota:MOGA’s performance, and, unfortunately, human intervention was required to discern the optimal design space. In conclusion, future work will attempt to reduce computational resource constraints by incorporating areduced-order neutronics model into the optimization cycle.

Accelerator-driven systems↗

Characterization of Deficiencies in the Frequency Domain Forced Response Analysis Technique for Supersonic Turbine Bladed Disks

Turbine blades in rocket and jet engine turbomachinery experience enormous harmonic loading conditions. These loads result from the integer number of upstream and downstream stator vanes as well as the other turbine stages. Assessing the blade structural integrity is a complex task requiring an initial characterization of whether resonance is possible and then performing a forced response analysis if that condition is met. The standard technique for forced response analysis in rocket engine turbines is to decompose a computational fluid dynamics (CFD).generated flow field into its harmonic components, and to then perform a frequency response analysis at the problematic natural frequencies using cyclically symmetric structural dynamic models. Recent CFD analysis and water-flow testing at NASA/MSFC, though, indicates that this technique may miss substantial harmonic and non ]harmonic excitation sources that become present in complex flows. This complex content can only be captured by a CFD flow field encompassing at least an entire revolution. A substantial development effort to create a series of software programs to enable application of the 360 degree forcing function in a frequency response analysis on cyclic symmetric models has been completed (to be described in a future paper), but the question still remains whether the frequency response analysis itself is capable of capturing the excitation content sufficiently. Two studies comparing frequency response analysis with transient response analysis, therefore, of bladed-disks undergoing this complex flow environment have been performed. The first is of a bladed disk with each blade modeled by simple beam elements and the disk modeled with plates (using the finite element code MSC/NASTRAN). The focus of this model is to be representative of response of realistic bladed disks, and so the dimensions are roughly equivalent to the new J2X rocket engine 1st stage fuel pump turbine. The simplicity of the model allows the CFD load to be able to be readily applied, along with analytical and experimental variations in both the temporal and spatial fourier components of the excitation. In addition, this model is a first step in identifying response differences between transient and frequency forced response analysis techniques. The second phase assesses this difference for a much more realistic solid model of a bladed-disk in order to evaluate the effect of the spatial variation in loading on blade dominated modes. Neither research on the accuracy of the frequency response method when used in this context or a comprehensive study of the effect of test-observed variation on blade forced response have been found in the literature, so this research is a new contribution to practical structural dynamic analysis of gas turbines. The primary excitation of the upstream nozzles interacts with the blades on fuel pump of the J2X causes the 5th Nodal diameter modes to be excited, as explained by Tyler and Sofrin1, so a modal analysis was first performed on the beam/plate model and the 5ND bladed-disk mode at 40167 hz was identified and chosen to be the one excited at resonance (see figure 1). The first forced response analysis with this model focuses on identifying differences between frequency and transient response analyses. A hypothesis going into the analysis was that perhaps the frequency response was enforcing a temporal periodicity that did not really exist, and so therefore it would overestimate the response. As high dynamic response was a considerable source of stress in the J2X, examining this concept could potentially be beneficial for the program.

Brown, Andrew M.↗