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 559 records · Page 31

Development of a Monte Carlo code for the data analysis of the {sup 18}F(p,α){sup 15}O reaction at astrophysical energies

Novae are astrophysical events (violent explosion) occurring in close binary systems consisting of a white dwarf and a main-sequence star or a star in a more advanced stage of evolution. They are called 'narrow systems' because the two components interact with each other: there is a process of mass exchange with resulting in the transfer of matter from the companion star to the white dwarf, leading to the formation of this last of the so-called accretion disk, rich mainly of hydrogen. Over time, more and more material accumulates until the pressure and the temperature reached are sufficient to trigger nuclear fusion reactions, rapidly converting a large part of the hydrogen into heavier elements. The products of 'hot hydrogen burning' are then placed in the interstellar medium as a result of violent explosions. Studies on the element abundances observed in these events can provide important information about the stages of evolution stellar. During the outbursts of novae some radioactive isotopes are synthesized: in particular, the decay of short-lived nuclei such as {sup 13}N and {sup 18}F with subsequent emission of gamma radiation energy below 511 keV. The gamma rays from products electron-positron annihilation of positrons emitted in the decay of {sup 18}F are the most abundant and the first observable as soon as the atmosphere of the nova starts to become transparent to gamma radiation. Hence the importance of the study of nuclear reactions that lead both to the formation and to the destruction of {sup 18}F. Among these, the {sup 18}F(p,α){sup 15}O reaction is one of the main channels of destruction. This reaction was then studied at energies of astrophysical interest. The experiment done at Riken, Japan, has as its objective the study of the {sup 18}F(p,α){sup 15}O reaction, using a beam of {sup 18}F produced at CRIB, to derive important information about the phenomenon of novae. In this paper we present the experimental technique and the Monte Carlo code developed to be used in the data analysis process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

PDC: Plant Dynamics Code for Design and Transient Analysis of Supercritical Brayton Cycles

The current report documents the Plant Dynamics Code developed at the Argonne National Laboratory for steady-state and transient analysis of supercritical Brayton cycles. The report consists of three major parts. Code Overview provides the general description of the code, including the purpose of the code development and the code’s main features. User’s Guide explains how to work with the code, including its user interface. Code Manual provides the detailed description of the code, including modeling, assumptions, and the solution schemes.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

$Z$-Pinch Interferometry Analysis With the Fourier-Based TNT Code

Here, we present the analysis of interferometry diagnostics with the user-friendly Talbot Numerical Tool (TNT), a Fourier-based postprocessing code that enables real-time assessment of plasma systems. TNT performance was explored with visible and infrared interferometry in pulsed-power-driven Z -pinch configurations to expand its capabilities beyond Talbot X-ray interferometry in the high-intensity laser environment. TNT enabled accurate electron density characterization of magnetically driven plasma flows and shocks through phase-retrieval methods that did not require data modification or masking. TNT demonstrated enhanced resolution, detecting below 4 % fringe shift, which corresponds to 8.7 × 10 15 cm –2 within 28 μ m, approaching the laser probing system limit. TNT was tested against a well-known interferometry analysis software, delivering an average resolving power nearly ten times better (~28 μ m versus ~ 210 μ m) when resolving plasma ablation features. TNT demonstrated higher sensitivity when probing sharp electron density gradients in supersonic shocks. A maximum electron areal density of 4.1 × 10 17 cm –2 was measured in the shocked plasma region, and a minimum electron density detection of ~ 1.0 × 10 15 cm –2 was achieved. When probing colliding plasma flows, the calculations of the effective adiabatic index and the associated errors were improved from γ* = 2.6 ± 1.6 – 1.4 ± 0.2 with TNT postprocessing, contributing valuable data for the interpretation of radiative transport. Additional applications of TNT in the characterization of pulsed-power plasmas and beyond are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development code for sensitivity and uncertainty analysis of input on the MCNPX for neutronic calculation in PWR core

This research was carried out on the development of code for uncertainty analysis is based on a statistical approach for assessing the uncertainty input parameters. In the butn-up calculation of fuel, uncertainty analysis performed for input parameters fuel density, coolant density and fuel temperature. This calculation is performed during irradiation using Monte Carlo N-Particle Transport. The Uncertainty method based on the probabilities density function. Development code is made in python script to do coupling with MCNPX for criticality and burn-up calculations. Simulation is done by modeling the geometry of PWR terrace, with MCNPX on the power 54 MW with fuel type UO2 pellets. The calculation is done by using the data library continuous energy cross-sections ENDF / B-VI. MCNPX requires nuclear data in ACE format. Development of interfaces for obtaining nuclear data in the form of ACE format of ENDF through special process NJOY calculation to temperature changes in a certain range.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Application of Performance Analysis Tools on SNL ASC Codes

This milestone 1) exercised a broad set of performance profiling and analysis tools, including tools whose development has been promoted by the ASC program; 2) exercised the tools on two different SNL ASC codes, one Sierra code (Sierra/Aria, a C++ codebase) and one RAMSES code (ITS, a Fortran codebase); and 3) exercised the tools on multiple platforms, including the CTS-1 (e.g., Serrano) and ATS-1 Trinity (e.g., Mutrino) platforms. The milestone generated a plethora of strong and weak scaling, trend and profile data for multiple versions and problem cases for each of the two codes. A wealth of experience was gained with the various tools that included identification of problems, an improved understanding of feature sets, enhanced usage documentation, and insights for future tool-development. Results are provided from a large number and variety of performance analysis runs with the target codes, together with instructions for how to make use of the tools with the codes.

97 MATHEMATICS AND COMPUTING↗

Cross-Code verification and sensitivity analysis to effectively model the electrothermal instability

This manuscript presents verification cases that are developed to study the electrothermal instability (ETI). Specific verification cases are included to ensure that the unit physics components necessary to model the ETI are accurate, providing a path for fluid-based codes to effectively simulate ETI in the linear and nonlinear growth regimes. Two software frameworks with different algorithmic approaches are compared for accuracy in their ability to simulate diffusion of a magnetic field, linear growth of the ETI, and a fully nonlinear ETI evolution. The nonlinear ETI simulations show early time agreement, with some differences emerging, as noted in the wavenumber spectrum, late into the nonlinear development of ETI. Here, a sensitivity study explores the role of equation-of-state (EOS), vacuum density, and vacuum resistivity. EOS and vacuum resistivity are found to be the most critical factors in the modeling of nonlinear ETI development.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Investigations into alternative radiation transport codes for ITER neutronics analysis

As radiation transport models for ITER become increasingly complex, the use of traditional constructive solid geometry (CSG) modelling approaches such as MCWNP presents significant disadvantages during model creation and integration. Studies were performed to assess the capabilities of alternative mesh-based neutronics workflows offered by the MCNP, DAG-MCNP and Serpent-2 codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of SAM Code Capabilities for Safety Analysis of GCR Air-ingress Events

Air-ingress following a depressurized loss-of-forced-cooling (DLOFC) event is a challenging, multiphysics safety scenario for High-Temperature Gas-Cooled Reactors (HTGRs), involving coupled gas composition transport, buoyancy-driven flow redistribution, graphite oxidation, and structural heat-up. Despite its importance — air ingress is a key scenario identified in the PIRT process for the HTGRs — existing system-level safety codes have lacked the integrated capability to simulate the complete event sequence with high confidence. This report documents the development, validation, and demonstration of three new capabilities in the SAM code to address this gap: (1) a multi-component gas mixture flow model with binary diffusion to track the helium-air composition and its effect on system density and flow; (2) a 0-D graphite oxidation model based on the Roes correlation, including oxygen consumption and exothermic heat release; and (3) an isentropic critical flow model for accurate representation of primary system depressurization through a break. These capabilities are validated against two benchmark experiments. The NSTF heavy-gas ingress experiment validates the multi-component flow model: SAM correctly reproduces the rapid buoyancydriven flow stagnation and subsequent natural circulation recovery driven by composition-dependent density changes. The NACOK graphite oxidation experiment validates the oxidation model: SAM predicts a bottom-level graphite weight loss of 25%, in close agreement with the measured 24%, and reproduces the strong axial nonuniformity and block-geometry dependence of oxidation, at a level comparable to the SPECTRA and TINTE codes. The validated capabilities are then exercised together in an integrated, reactor-scale simulation of a DLOFC air-ingress transient in a simplified HTR-PM pebble-bed reactor. In a single calculation spanning approximately 8 days, SAM reproduces the complete accident sequence: rapid depressurization, densityand diffusion-driven air ingress over ˜15 hours, onset of buoyancy-driven natural circulation, exothermic graphite oxidation with a peak fuel temperature at ˜62 hours, and eventual passive cooldown. These results demonstrate that SAM now provides the nuclear community with a preliminarily validated, modern systemlevel tool for HTGR air-ingress safety analysis, filling a recognized capability gap. Future extensions to broaden species tracking, improve oxidation chemistry, and refine the reactor model are discussed.

Yang, Gang↗

Testing the Jump Finding Code for NEXUS Qubit Analysis

At Fermilab, the NEXUS dilution refrigerator is repeating the same experiment from Wilen et al. 1 underground to further investigate the effect of ionizing radiation and cosmic ray muons on qubit decoherence. This refrigerator is 100m underground in the MINOS cavern to significantly lower the incident rate of cosmic ray muons, meaning the majority of incidents is due to ionizing radiation. To gain a better understanding of how the rate of ionizing radiation affects the qubits, we use 4 radiation configurations to collect data from the chip. In this work, we examine the efficiency of our jump detection code to find jumps of various sizes. We do so using simulated data with a known number of injected jumps of known sizes. These jumps are charge offsets of the qubit, thought to be due to quasiparticle poisoning. To test this, we run the jump finding code on an array of simulated data in which two main parameters, the smoothing factor and the threshold at which a charge offset is considered a jump, are changed to test an array of these parameters. We then perform a χ2 analysis to find the best combination of these two parameters.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Genetic algorithms applied to reconstructing coded imaging of neutrons and analysis of residual watermark

Monte-Carlo simulation of neutron coded imaging based on encoding aperture for Z-pinch of large field-of-view with 5 mm radius has been investigated, and then the coded image has been obtained. Reconstruction method of source image based on genetic algorithms (GA) has been established. 'Residual watermark,' which emerges unavoidably in reconstructed image, while the peak normalization is employed in GA fitness calculation because of its statistical fluctuation amplification, has been discovered and studied. Residual watermark is primarily related to the shape and other parameters of the encoding aperture cross section. The properties and essential causes of the residual watermark were analyzed, while the identification on equivalent radius of aperture was provided. By using the equivalent radius, the reconstruction can also be accomplished without knowing the point spread function (PSF) of actual aperture. The reconstruction result is close to that by using PSF of the actual aperture.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Validation and Code-to-Code Verification of FAST for a Megawatt-Scale Wind Turbine with Aeroelastically Tailored Blades

This paper presents validation and code-to-code verification of the latest version of the U.S. Department of Energy, National Renewable Energy Laboratory wind turbine aeroelastic engineering simulation tool, FAST v8. A set of 1,141 test cases, for which experimental data from a Siemens 2.3 MW machine have been made available and were in accordance with the International Electrotechnical Commission 61400-13 guidelines, were identified. These conditions were simulated using FAST as well as the Siemens in-house aeroelastic code, BHawC. This paper presents a detailed analysis comparing results from FAST with those from BHawC as well as experimental measurements, using statistics including the means and the standard deviations along with the power spectral densities of select turbine parameters and loads. Results indicate a good agreement among the predictions using FAST, BHawC, and experimental measurements. Here, these agreements are discussed in detail in this paper, along with some comments regarding the differences seen in these comparisons relative to the inherent uncertainties in such a model-based analysis.

17 WIND ENERGY↗

RAVEN, a New Software for Dynamic Risk Analysis

RAVEN is a generic software driver to perform parametric and probabilistic analysis of code simulating complex systems. Initially developed to provide dynamic risk analysis capabilities to the RELAP-7 code [1] is currently being generalized with the addition of Application Programming Interfaces (APIs). These interfaces are used to extend RAVEN capabilities to any software as long as all the parameters that need to be perturbed are accessible by inputs files or directly via python interfaces. RAVEN is capable to investigate the system response probing the input space using Monte Carlo, grid strategies, or Latin Hyper Cube schemes, but its strength is its focus toward system feature discovery like limit surfaces separating regions of the input space leading to system failure using dynamic supervised learning techniques. The paper will present an overview of the software capabilities and their implementation schemes followed by same application examples.

97 MATHEMATICS AND COMPUTING↗

In-Depth Analysis of Simulation Engine Codes for Comparison with DOE's Roof Savings Calculator and Measured Data

The Roof Savings Calculator (RSC) was developed through collaborations among Oak Ridge National Laboratory (ORNL), White Box Technologies, Lawrence Berkeley National Laboratory (LBNL), and the Environmental Protection Agency in the context of a California Energy Commission Public Interest Energy Research project to make cool-color roofing materials a market reality. The RSC website and a simulation engine validated against demonstration homes were developed to replace the liberal DOE Cool Roof Calculator and the conservative EPA Energy Star Roofing Calculator, which reported different roof savings estimates. A preliminary analysis arrived at a tentative explanation for why RSC results differed from previous LBNL studies and provided guidance for future analysis in the comparison of four simulation programs (doe2attic, DOE-2.1E, EnergyPlus, and MicroPas), including heat exchange between the attic surfaces (principally the roof and ceiling) and the resulting heat flows through the ceiling to the building below. The results were consolidated in an ORNL technical report, ORNL/TM-2013/501. This report is an in-depth inter-comparison of four programs with detailed measured data from an experimental facility operated by ORNL in South Carolina in which different segments of the attic had different roof and attic systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Design Analysis of SNS Target StationBiological Shielding Monoligh with Proton Power Uprate

This report documents the analysis of the dose rate in the experiment area outside the Spallation Neutron Source (SNS) target station shielding monolith with proton beam energy of 1.3 GeV. The analysis implemented a coupled three dimensional (3D)/two dimensional (2D) approach that used both the Monte Carlo N-Particle Extended (MCNPX) 3D Monte Carlo code and the Discrete Ordinates Transport (DORT) two dimensional deterministic code. The analysis with proton beam energy of 1.3 GeV showed that the dose rate in continuously occupied areas on the lateral surface outside the SNS target station shielding monolith is less than 0.25 mrem/h, which complies with the SNS facility design objective. However, the methods and codes used in this analysis are out of date and unsupported, and the 2D approximation of the target shielding monolith does not accurately represent the geometry. We recommend that this analysis is updated with modern codes and libraries such as ADVANTG or SHIFT. These codes have demonstrated very high efficiency in performing full 3D radiation shielding analyses of similar and even more difficult problems.

42 ENGINEERING↗

Sensitivity Studies of the Inter-subchannel Mixing in Subchannel Analysis

In the PWR reactor core, the power profile is nonuniform and the inter-subchannel mixing is of great importance for flatting the temperature distribution of the coolant and for reactor core safety analysis. Subchannel codes are widely used in the reactor core safety analysis to assess the thermal hydraulic behavior of fuel assemblies. One of basic assumptions in subchannel analysis is that axial flow predominates in the channel flow. Under this assumption, the lateral momentum equation in subchannel codes can be simplified without bringing in significant errors. The inter-subchannel mixing leading to the exchange of mass, energy and momentum are modeled by adding source terms to the hydrodynamic equations or by separate models, and most of source terms and models are empirical correlations. With the advancement of mixing vane grid design, the lateral mixing has becoming a critical components in the evaluation of the inter-subchannel mixing in flow field. The accuracy of the subchannel analysis is strongly dependent on the modeling of inter-subchannel interactions between adjacent subchannels. So the inter-subchannel mixing models should be well considered and developed in subchannel analysis to better reflect their influence on flow field and temperature field in rod bundle. Generally speaking, the inter-subchannel mixing in subchannel analysis contains turbulent mixing and crossflow mixing in single phase flow. When it comes to two phase flow, the void drift will also contribute to inter-subchannel mixing [2]. Both the mechanisms and the generation conditions of each type of inter-subchannel mixing are different from each other. Turbulent mixing is a result of random turbulent flow and pressure fluctuations. Crossflow is caused by a net transverse pressure difference between adjacent subchannels. Lateral pressure differences can result from different subchannel hydraulic diameters, mass flux, heat flux distributions and changes in flow areas caused, for example, by wire wraps and spacer grids, respectively. Two-phase void drift is due to the tendency of the vapor phase to redistribute itself to a preferred equilibrium void distribution. This paper focuses on the sensitivity study of the inter-subchannel mixing with mixing vane grids for single phase flow in subchannel analysis. Both the turbulent mixing and crossflow mixing of the 5 x 5 rod bundle with mixing vane grids are predicted, and the results are compared and analyzed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Binary Analysis with Architecture and Code Section Detection using Supervised Machine Learning

When presented with an unknown binary, which may or may not be complete, having the ability to determine information about it is critical to future reverse engineering, particularly in discovering the binary’s intended use and potentially malicious nature. This paper details techniques to both identify the machine architecture of the binary, as well as to locate the important code segments within the file. This identification of unknown binaries makes use of a technique called byte histogram in addition to various machine learning (ML) techniques, which we call “What is it Binary” or WiiBin. Benefits of byte histograms reflect the simplicity of calculation and do not rely on file headers or metadata, allowing for acceptable results when only a small portion of the original file is available. Utilizing WiiBin, we were able to accurately (>80%) determine the architecture of test binaries with as little as a 20% contagious portion of the file present. We were also able to determine the location of code sections within a binary by utilizing the WiiBin framework. Ultimately, the more information that can be gleaned from a binary file, the easier it is to successfully reverse engineer.

99 GENERAL AND MISCELLANEOUS↗