Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code Analysis”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI1996

Nuclear Safety [Vol. 37, No. 2, April-June 1996]

Nuclear Safety is a journal that covers significant issues in the field of nuclear safety. Its primary scope is safety in the design, construction, operation, and decommissioning of nuclear power reactors worldwide and the research and analysis activities that promote this goal, but it also encompasses the safety aspects of the entire nuclear fuel cycle, including fuel fabrication, spent-fuel processing and handling, and nuclear waste disposal, the handling of fissionable materials and radioisotopes, and the environmental effects of all these activities. Table of Contents for this issue follows. GENERAL SAFETY CONSIDERATIONS: 97 The Nuclear Community and the Public: Cognitive and Cultural Influences on Thinking About Nuclear Risk, M. A. Meyer; 109 Twenty-Third Water Reactor Safety Information Meeting, D. A. Copinger; ACCIDENT ANALYSIS: 126 Analysis of a PWR LBLOCA Without SCRAM, Trevor N. Tyler, Rafael Macian-Juan and John H. Mahaffy; DESIGN FEATURES: 139 Vulnerability of Multiple-Barrier Systems, N. C. Lind; ENVIRONMENTAL EFFECTS: 149 A Study of Wet Catalytic Oxidation of Radioactive Spent Ion Exchange Resin by Hydrogen Peroxide, Xingchao Jian, Tianbao Wu, and Guichun Yun; 157 A Comparison Study and Resolution of Differences Between Emergency Response and Safety Analysis Codes Used at the Savannah River Site, A. A. Simpkins; OPERATING EXPERIENCES: 164 Reactor Shutdown Experience, Compiled by J. W. Cletcher; RECENT DEVELOPMENTS: 167 Reports, Standards, and Safety Guides, D. S. Queener; 172 Proposed Rule Changes as of Dec. 31,1995; ANNOUNCEMENTS: 178 American Nuclear Society 1997 Annual Meeting; 178 American Nuclear Society Nuclear Criticality and Safety Division Topical Meeting; 176 The Authors.

05 NUCLEAR FUELS↗
DOE OSTI2020

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Pyomo.DOE: An open-source package for model-based design of experiments in Python

Predictive mathematical models are a cornerstone of science and engineering. Yet selecting, calibrating, and validating said science-based models often remains an art in practice. Model-based design of experiments (MBDoE) provides a systematic framework to maximize information gain from experiments while minimizing time and resource costs. But MBDoE remains limited to niche application areas, in part because practitioners must integrate expertise in statistics, computational optimization, and modeling. To help reduce this barrier, we introduce Pyomo.DOE, an open-source package for MBDoE. Pyomo.DOE uses a nonlinear sensitivity analysis code k_aug to quickly approximate the Fisher information matrix and leverages a new stochastic programming abstraction. We demonstrate Pyomo.DOE with the first application of MBDoE to fixed-bed breakthrough experiments, which highlights the power of Pyomo.DOE to quantify the value of experimental modifications a priori for large-scale partial differential-algebraic equation (PDAE) models. Here we also provide a mathematical primer on MBDoE targeted at general chemical engineers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Two-zone stratified wetwell model development and implementation for RELAP-7

The Reactor Core Isolation Cooling (RCIC) system consists of a turbine, pump and wetwell and is used in many US Boiling Water Reactors (BWRs) as an important heat removal system. The wetwell is a major heat sink within containment. Experimental investigations have shown that thermal stratification may occur in the wetwell during RCIC operation. Current systems-level analysis codes lack stratified wetwell models. Here, this paper describes the development of a stratified wetwell model for the RELAP-7 code. This paper uses a simple two-zone model for the wetwell water that can capture thermal stratification. This model assumes that a buoyant plume transports heat and mass from the steam injection site to the upper water zone. The buoyant plume assumption requires that steam is injected at a low enough rate so it adds negligible momentum. Model validation against experimental data shows similar trends of temperature field development in the suppression pool.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

A multi-dimensional dataset for two-phase instability in low pressure natural circulation based on direct transient local measurement

Thermal hydraulic instability can challenge the control and the safety of nuclear reactors. This experimental research captures and investigates the instability in a low-pressure (1–10 atmospheres) natural circulation loop with an adiabatic chimney installed above the heated section. Direct transient local measurement of two-phase flow is performed at multiple axial locations along the flow channel with thermocouples and conductivity probes traversed by customized high-precision mechanisms. Based on the periodicity of oscillations, a post-processing method is proposed to align the data captured non-simultaneously to reconstruct the instability behavior in one representative period. The averaging operation over multiple oscillation cycles is performed to reduce the random errors. Cyclic phenomena are clearly demonstrated by the post-processed data, and rapid void fraction change induced by the flashing is observed to be significant in the driving mechanism for low pressure oscillations. Parametric effects of inlet subcooling, heat flux, inlet flow restriction and system pressure are studied within the experiment range, and qualitative results are explained and compared with existing literature. The ratio between the oscillation period and the travelling time of the fluid particle through the test section is compared with the theoretical value for the density wave oscillation. The new dataset contains directly measured transient multi-dimensional two-phase data along the flow channel for natural circulation under density wave oscillations. The generated dataset is valuable for the validation of system analysis codes and CFD codes in the simulation of natural circulation phenomena and dynamic two-phase flow.

42 ENGINEERING↗
DOE OSTI2023

Hermeticity of SiC/SiC composite and monolithic SiC tubes irradiated under radial high-heat flux

Demonstration of hermetic SiC fiber–reinforced SiC matrix composite cladding under normal operating environments has been identified as one of the most critical feasibility issues for accident-tolerant fuel cladding in light-water reactors. This study provides critical experimental data needed for understanding the effects of irradiation on hermeticity. SiC composite and monolithic tubes were neutron-irradiated to 2 displacements per atom with and without a nominal radial heat flux of 0.6 MW/m2 to produce a simulated in-pile stress state for the normal operation of a light-water reactor. The through-thickness temperature gradient under irradiation results in a gradient in swelling, which causes a significant stress buildup. Such irradiation-induced stress was modeled using a commercial finite element analysis code. The radial heat flux–irradiation synergism was experimentally investigated by constructing a special irradiation capsule and evaluating the helium hermeticity of the specimens. The simulated stress state exhibited a near equi-biaxial tensile axial and hoop stress of ~150 MPa at the inner surface of the SiC composite tube. This stress level is potentially beyond the matrix cracking stress. Degradation of hermeticity of the SiC composite tubes was observed after irradiation, indicating irradiation-induced cracking, whereas the irradiated monolithic SiC tubes remained hermetic. In conclusion, the results indicate that loss of hermeticity caused by radiation-induced microcracking is a potential issue for SiC composite cladding, depending on the magnitude of the temperature gradients. Coating the outer surface of the cladding was identified as a mitigation strategy that might overcome the cracking issue.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modernizing the active neutron collar poison rod correction using simulations to enhance PWR fresh fuel verification

The declared inventory of 235 U in fresh, low enriched uranium, nuclear fuel assemblies is routinely verified by nuclear safeguards inspectorates using the UNCL (Uranium Neutron Collar – Light Water Reactor Fuel) nondestructive assay instrument. The trend in modern fuels is towards higher initial enrichment, which in turn requires a larger number of burnable poison pins with higher Gadolinia concentration to hold down the initial reactivity. UNCL assay error for modern fuel assemblies typically exceeds 10%. The traditional algorithm used to correct the response for burnable poison content needs revision to achieve accuracy comparable to the $\sigma_\text{R}$ of 4.5% expected for poison-free assemblies. Here we review the problem and use a large set of 287 Monte Carlo simulations based on an experimentally benchmarked model corresponding to the 16×16 PWR array used in Brazilian Angra type II and III fuel. Simulated relative responses are used to evaluate the functional form of the poison correction and determine parameters optimal for 16×16 PWR assemblies. This update is conducted using between 4–24 Gadolinia pins ranging 2–11 wt% Gd 2 O 3 with assemblies having mean enrichments ranging 2.5%–5%. Experimental validation of the updated coefficients using seven measurements shows bias is reduced by an order of magnitude and $\sigma_\text{R}$R reduced from 10% to less than 2%. Two updated sets of coefficients are provided, the most accurate of which is directly useable in the existing analysis code (INCC) in use by inspectorate at fuel fabrication facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Development of the reactor core isolation cooling system model for the extended station black-out accident

In order to accurately predict the behavior of the Reactor Core Isolation Cooling System (RCIC) in the context of the FLEX strategies under the extended station black-out accident (SBO), we present a mechanistic model of the RCIC in this work. The control volume formulation of angular momentum conservation is employed to model the Terry turbine in the RCIC system. The two-phase nozzle flow model developed in the previous work is utilized to simulate the convergent-divergent nozzles in the Terry turbine. A nozzle flow table is generated to cover the operation range of the RCIC during beyond design basis accidents (BDBAs). The mechanistic RCIC turbine-pump model and the nozzle flow table are incorporated in a TRACE BWR model. Simulations of short-term and long-term SBO transients are performed to demonstrate dynamic RCIC responses. A system analysis code with the proposed mechanistic model of the RCIC system will have a better prediction of system safety response during SBO accidents. With the addition of such more realistic models, the system code will be capable of providing expanded understanding of safety margin, which can help help to form the technical basis for the successful implementation of the FLEX and SAMG measures under Extended Loss of AC Power (ELAP) conditions for BWR or PWR reactor designs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Core design and performance of the Westinghouse lead fast reactor with UO 2 and MOX configurations

For this work, Westinghouse partnered with Argonne National Laboratory to design, model and optimize UO 2 - and MOX-fueled core designs for a medium size (950 MWt) Lead Fast Reactor that was pursued by Westinghouse. Using Argonne’s suite of reactor analysis codes together with Westinghouse fuel cost economic models, thousands of candidate cores were considered to achieve the economics-optimized cores presented in this paper. This optimization process considered detailed reactor physics, fuel performance, transient performance, and fuel economics models. The reactor performance of the resulting optimized UO 2 - and MOX-fueled core designs are described and compared in this paper. Both cores show fuel performance and transient behavior that is considered acceptable for the optimization presented herein, while further testing campaigns on material performance in high-temperature liquid lead will be required to confirm acceptability at the operating conditions chosen. A multi-batch strategy was selected for the UO 2 core for best fuel utilization with minimum fuel inventory costs. A single-batch fuel management was instead selected for the MOX core to maximize cycle length and minimize the impact of the longer refueling outage resulting from the higher decay heat of the discharged MOX fuel relative to the discharged UO 2 fuel, requiring a longer cooling time before dry-lift of discharged fuel could take place.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Probing the antiknock effect of anisole through an ignition, speciation and modeling study of its blends with isooctane

In order to unravel the reaction pathways relevant to anisole co-oxidation within a fuel blend, a detailed study of isooctane/anisole blends was performed with the ULille RCM. Ignition delays as well as mole fraction profiles were measured during a two-stage ignition delay using sampling and GC techniques. These results are used to validate a kinetic model developed from ab initio calculations for the most relevant rate constants which included H-atom abstraction reactions from anisole, and reactions on the potential energy surfaces of methoxyphenyl + O 2 and anisyl + O 2 . Pressure dependent rate constants were computed for the methoxyphenyl + O 2 and anisyl + O 2 reactive systems using master equation code analysis. The new kinetic model shows good agreement with the experimental data. Dual brute-force sensitivity analysis was performed, on both first- and second-stages of ignition, allowing the identification of the most important reactions in the prediction of both ignition delays. It was observed that while pure anisole does not show NTC behavior, a 60/40 isooctane/anisole blend displays such behavior, as well as two-stage ignition. This suggests anisole addition may not be as beneficial to knock resistance as expected from its high octane number. The kinetic modeling results demonstrate the importance of H-abstraction reactions both from the methoxy group and from the aryl ring in ortho-position and the addition of the resultant radicals to O2, mostly leading to the formation of polar or non-aromatic products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

An automatically curated first-principles database of ferroelectrics

Ferroelectric materials have technological applications in information storage and electronic devices. The ferroelectric polar phase can be controlled with external fields, chemical substitution and size-effects in bulk and ultrathin film form, providing a platform for future technologies and for exploratory research. In this work, we integrate spin-polarized density functional theory (DFT) calculations, crystal structure databases, symmetry tools, workflow software, and a custom analysis toolkit to build a library of known, previously-proposed, and newly-proposed ferroelectric materials. With our automated workflow, we screen over 67,000 candidate materials from the Materials Project database to generate a dataset of 255 ferroelectric candidates, and propose 126 new ferroelectric materials. We benchmark our results against experimental data and previous first-principles results. The data provided includes atomic structures, output files, and DFT values of band gaps, energies, and the spontaneous polarization for each ferroelectric candidate. Finally, we contribute our workflow and analysis code to the open-source python packages atomate and pymatgen so others can conduct analogous symmetry driven searches for ferroelectrics and related phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantifying dislocation-type defects in post irradiation examination via transfer learning

The quantitative analysis of dislocation-type defects in irradiated materials is critical to materials characterization in the nuclear energy industry. The conventional approach of an instrument scientist manually identifying any dislocation defects is both time-consuming and subjective, thereby potentially introducing inconsistencies in the quantification. This work approaches dislocation-type defect identification and segmentation using a standard open-source computer vision model, YOLO11, that leverages transfer learning to create a highly effective dislocation defect quantification tool while using only a minimal number of annotated micrographs for training. This model demonstrates the ability to segment both dislocation lines and loops concurrently in micrographs with high pixel noise levels and on two alloys not represented in the training set. Inference of dislocation defects using transmission electron microscopy on three different irradiated alloys relevant to the nuclear energy industry are examined in this work with widely varying pixel noise levels and with completely unrelated composition and dislocation formations for practical post irradiation examination analysis. Code and models are available at https://github.com/idaholab/PANDA.

36 MATERIALS SCIENCE↗
DOE OSTI2024

ROOT’s RNTuple I/O Subsystem: The Path to Production

The RNTuple I/O subsystem is ROOT’s future event data file format and access API. It is driven by the expected data volume increase at upcoming HEP experiments, e.g. at the HL-LHC, and recent opportunities in the storage hardware and software landscape such as NVMe drives and distributed object stores. RNTuple is a redesign of the TTree binary format and API and has shown to deliver substantially faster data throughput and better data compression both compared to TTree and to industry standard formats. In order to let HENP computing workflows benefit from RNTuple’s superior performance, however, the I/O stack needs to connect efficiently to the rest of the ecosystem, from grid storage to (distributed) analysis frameworks to (multithreaded) experiment frameworks for reconstruction and ntuple derivation. With the RNTuple binary format soon arriving at its first production release, we present RNTuple’s feature set, integration efforts, and its performance impact on the time-to-solution. We show the latest performance figures of RDataFrame analysis code of realistic complexity, comparing RNTuple and TTree as data sources. We discuss RNTuple’s approach to functionality critical to the HENP I/O (such as multithreaded writes, fast data merging, schema evolution) and we provide an outlook on the road to its use in production.

Blomer, Jakob↗
DOE OSTI2021

Initial operation and data processing on a system for real-time evaluation of Thomson scattering signals on the Large Helical Device

A scalable system for real-time analysis of electron temperature and density based on signals from the Thomson scattering diagnostic, initially developed for and installed on the NSTX-U experiment, was recently adapted for the Large Helical Device and operated for the first time during plasma discharges. During its initial operation run, it routinely recorded and processed signals for four spatial points at the laser repetition rate of 30 Hz, well within the system’s rated capability for 60 Hz. We present examples of data collected from this initial run and describe subsequent adaptations to the analysis code to improve the fidelity of the temperature calculations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Thermal energy mitigation and toroidal peaking effects in JET disruptions

Previous investigations on JET suggest half or less of plasma stored thermal energy W th is radiated (f rad,th ≲ 0.5) using either massive gas injection (MGI) or shattered pellet injection (SPI) disruption mitigation. We investigate whether the apparent incomplete f rad,th is explained by radiation peaking near the injection plume. High toroidal peaking throughout the pre-thermal quench is found in argon–deuterium MGI on JET, with typically >3× higher radiation near the injector than toroidally distant. Previously unexplained toroidal bolometry measurements in neon–deuterium SPI are reproduced with similar peaking using the Emis3D radiation analysis code. These observations align with results from Alcator C-Mod and KSTAR. This peaking is not captured by previous JET studies that found poor thermal mitigation. Two sets of neon–deuterium SPI and two sets of argon–deuterium MGI are analyzed using Emis3D. In SPI, f rad,th rises from no-plume estimates of 0.31 and 0.66 to lower bounds of 0.84 and 0.92, respectively, and f rad,th ∼ 1 is possible. In MGI, the toroidal spread of the peaking feature is poorly constrained. f rad,th up to 0.85 and 0.65 are possible using the largest possible spread, increasing from 0.42 and 0.28, although f rad,th ∼ 1 does not appear to be reached. Revised mitigation estimates on JET suggest a lower melt risk to the divertor in mitigated disruptions on ITER and SPARC than previously thought. However, peaking near injectors could increase flash melting risk on nearby plasma facing components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Thermal-Hydraulic Analysis of the AFIP-7 Irradiation Test in the Advanced Test Reactor - Oxide Growth Prediction and Correlation

Knowing the thickness of the oxide layer on the surface of aluminum fuel cladding is vitally important for predicting fuel temperature due to the low thermal conductivity of the oxide layer. Several correlation models for predicting oxide growth can be found in the literature. In previous research, the correlations were combined with heat transfer simulations in Abaqus, a finite element analysis code, to forecast the oxide growth. However, this approach requires heat transfer coefficients for modeling heat exchanges with the external flow field, and such coefficients were obtained through empirical equations. Since different empirical equations yield varying heat transfer coefficients, the cladding temperature and predicted oxide thickness both carry a high degree of uncertainty. Here, this research develops a new approach that integrates the fluid flow, fluid and solid heat transfer, and oxide growth correlation(s) into a single computational fluid dynamics model. We demonstrate this approach’s ability to predict oxide development on the AFIP-7 plates during two Advanced Test Reactor (ATR) irradiation cycles. The projected oxide thickness falls within the experimental measurements taken during post irradiation examination.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Isotopic mass effects of tritium-fueled high-performance TFTR supershots

An increase in total stored energy correlated with the addition of tritium fuel was observed in supershots during the TFTR DT campaign. This supershot regime had strikingly high, centrally peaked ion and electron temperatures, and the largest neutron emission rates observed in TFTR. This paper presents a study of the causes of this increase in stored energy in supershots. Twenty-six supershots have been recently reanalyzed with the TRANSP plasma analysis code. Early TRANSP simulations did not accurately match the measured magnitude and time evolutions of the neutron emission rates. This mismatch is attributed to neglecting apparent increases of trace amounts of heavy impurities during neutral beam injection. The new TRANSP runs were tuned to accommodate this and match the measured global neutron emission rates. These new runs also had improved fidelity in predicting the time histories and radial dependencies of measured DT neutron emission rates. That in turn adds confidence in the simulated thermal deuterium and tritium density profiles that are needed for calculating the average hydrogenic atomic mass profiles. Six subsets of these supershots had well matched toroidal field B tor , plasma current I p , flux geometry, and total injected neutral beam power. The mix of D and T beam ions was varied for different discharges. The magnitude of the increase of the thermal ion energy W i with added tritium was relatively small, and the total W tot increase is dominated by the increase in fast beam ions with T. Analyses at times before the occurrences of deleterious MHD instabilities yielded scaling of W tot with the volume-average isotopic mass consistent with previous publications. The relative fraction of fast energy ions is expected to be small in practical tokamak reactors. Thus the increase in stored energy W tot observed in TFTR supershots does not appear likely to be significantly helpful for producing useful fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗