Search NASA⌕ Search

SEARCH · Search NASA

Results for “tracking code”

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 235 records · Page 13

Simulation scaling studies of reactor core two-phase flow using direct numerical simulation

Tremendous growth in supercomputing power in recent years has resulted in the emergence of high-resolution flow analysis methods as an advanced research tool to evaluate single and two-phase flow behavior. In particular, unstructured mesh-based methods have been applied to analyze flows in complex reactor core geometries, including those of light water reactors (LWR). The finite-element based code, PHASTA, is utilized to perform large-scale simulations of two-phase bubbly flows in LWR geometries. Given the large computational cost of direct numerical simulation (DNS) coupled with interface tracking methods (ITM), typical domains encompass a portion of a single subchannel. In the presented research, the state-of-the-art analysis of turbulent two-phase flows in complex LWR subchannel geometries are demonstrated at both prototypical reactor parameters as well as scaled low pressure conditions. Three different cases are studied, a high-pressure simulation in prototypical reactor subchannel geometry, a low-pressure case in prototypical geometry and a final low-pressure case in a geometry scaled up to conserve the ratio between the bubble size and the domain pitch. Utilizing advanced statistical processing tools, these simulation conditions are compared to shed light on the relevancy of two-phase flow characteristics given the significant differences between LWR and low-pressure conditions. These findings can lead to the generation of useful guiding principles when researchers need to scale the two-phase flow behavior captured at low pressure and temperature conditions to those at reactor operating conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Impact of Nuclear Physics Uncertainties on Galactic Chemical Evolution Predictions

Modeling the evolution of the elements in the Milky Way is a multidisciplinary and challenging task. In addition to simulating the ~13 billion years evolution of our Galaxy, chemical evolution simulations must keep track of the elements synthesized and ejected from every astrophysical site of interest (e.g., supernova, compact binary merger). The elemental abundances of such ejecta, which are a fundamental input for chemical evolution codes, are usually taken from theoretical nucleosynthesis calculations performed by the nuclear astrophysics community. Therefore, almost all chemical evolution predictions rely on the nuclear physics behind those calculations. In this proceedings article, we highlight the impact of nuclear physics uncertainties on galactic chemical evolution predictions. We demonstrate that nuclear physics and galactic evolution uncertainties both have a significant impact on interpreting the origin of neutron-capture elements in our Solar System. Those results serve as a motivation to create and maintain collaborations between the fields of nuclear astrophysics and galaxy evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Fission fragment angular anisotropy in neutron-induced fission of U 235 measured with a time projection chamber

Fission fragment angular distributions can provide an important constraint on fission theory, improving predictive fission codes, and are a prerequisite for a precise ratio cross section measurement. Available anisotropy data are sparse, especially at neutron energies above 5 MeV. For the first time, a three-dimensional tracking detector is employed to study fragment emission angles and provide a direct measurement of angular anisotropy. The Neutron Induced Fission Fragment Tracking Experiment (NIFFTE) collaboration has deployed the fission time projection chamber to measure nuclear data with unprecedented precision. The fission fragment anisotropy of 235 U has been measured over a wide range of incident neutron energies from 180 keV to 200 MeV; a careful study of the systematic uncertainties complement the data.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

A miniaturized long-life acoustic transmitter for animal tracking in marine environments

Marine and hydrokinetic (MHK) energy devices are early-stage technologies and their impacts to ocean life and environment are largely unknown, which consequently limits the progress of deploying these devices at commercial scale. Optical and acoustic imaging have been the primary technologies for monitoring marine animals’ behavior at existing or potential MHK sites. However, these technologies lack the ability to identify individual animals, which makes them less useful in identifying protected species as well as studying the animals’ long-term behavior. Implantable acoustic transmitters with unique identification codes could be used as a complementary technology. To minimize the tag burden on the tagged animal (and hence the resulting bias) and to satisfy the long-term monitoring needs, it is desired that the implanted transmitter be as small and light as possible and can last for as long as possible. To date, the smallest acoustic transmitters suitable for marine animal tracking have relatively short service life and detection range (up to 200 m). In this report, we introduce a new miniaturized acoustic transmitter which operates at 208.3 kHz. It offers a smaller dry weight (0.45 g), a significantly longer service life (140 days at a 5-second ping rate) as well as a greater detection range than the existing commercial counterparts, which will help obtain more comprehensive movement and behavioral information of animals around MHK energy devices.

Li, Huidong↗

Machine learning-driven descriptions of protein dynamics at solid-liquid interfaces

This chapter has described how ML has enabled quantitative analysis of HS-AFM data to discover the physical phenomena governing protein dynamics and ordering at solid-liquid interfaces. The research detailed in this chapter modeled the rotation models of protein nanorods, the discovery of which would otherwise not be possible. By tracking the trajectories of individual protein rods from frame to frame, it was possible to model Brownian type motion and behaviors and Levy-flight dynamics that had not previously been shown. We also described the application of the Python package AtomAI, which has been developed specifically to analyze and extract physical phenomena, providing exemplar code for training an ensemble of deep neural networks to produce the semantic segmentation of AFM data and functions for encoding and decoding local environments. We last described a combinatorial approach to analyze very noisy data with a densely covered substrate where the emergence of order for the protein liquid crystals could be elucidated. By combining the methods from Case 1 and 2, it was possible to obtain the center of mass and angle for each rod in the images and track the assembly of the rods over time into a 2D liquid crystal array on the surface of mica.

protein dynamics, solid-liquid interfaces, atomic ↗

bifacialVF - Bifacial PV View Factor model 2022 [SWR 17-58]

bifacialvf is a self-contained view factor (or configuration factor) model which replicates a 5-row PV system of infinite extent perpendicular to the module rows. Single-axis tracking is supported, and hourly output files based on TMY inputs are saved. Spatial nonuniformity is reported, with multiple rear-facing irradiances collected on the back of each module row. Bilinear interpolation code add-on to bifacialvf (description below) to pre-generate IV arrays and bifacial coefficients, and to examine the energy production with back side irradiance mismatch for either a portrait or landscape module. Included are IV curves and bifacial info for a Yingli (standard) module. Based on the publication: "A Practical Irradiance Model for Bifacial PV Modules" B. Marion, S. MacAlpine, C. Deline, A. Asgharzadeh, F. Toor, D. Riley, J. Stein, C. Hansen 2017 IEEE Photovoltaic Specialists Conference, Washington DC, 2017 URL: https://www.nrel.gov/docs/fy17osti/67847.pdf

Deline, Chris↗

Implementation of the Windowed Multipole Method in Shift

The windowed multipole (WMP) method has been implemented in the Shift Monte Carlo (MC) radiation transport code with support for both CPU and GPU execution. With this method, small WMP data libraries (~100 MB) can be used to accurately Doppler broaden cross sections to arbitrary temperatures “on the fly” during an MC simulation. This approach yields significant memory savings relative to traditional methods, making it ideal for high-fidelity analysis such as coupled multiphysics simulations. This document provides the exact forms of the WMP equations used by Shift, as well as a detailed description of the structure of WMP HDF5 data files provided by the Massachusetts Institute of Technology (MIT). The Shift implementation has been validated against the OpenMC radiation transport code, with excellent agreement demonstrated for 70 nuclides across an operative range of temperatures. CPU and GPU performance testing using a small module reactor (SMR) problem demonstrated that this method decreases the neutron tracking rate by a factor of ~2 on the Summitdev machine. A new set of WMP data being developed in-house will employ novel methods to improve tracking rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Battery Lifecycle Framework

The Battery Lifecycle (BLC) Framework is an open-source platform that provides tools to visualize and share battery data from material characterization, cell testing, manufacturing, and field testing through the technology development cycle. BLC has three components: data importers, a front-end for querying the data and creating visualizations, and an application programming interface to provide access to the data from Python. BLC has been deployed for tracking the development of a battery from the lab to a manufacturing line and systems installed in the field and for comparing studies of multiple cells of the same battery chemistry and configuration. The code was developed around Redash, a robust open-source extract-transform-load engine. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-4546 O

De Angelis, Valerio↗

Generator Scorecard

Code built on Archive Walker. The Generator Scorecard analyzes power grid measurements to automate the process of evaluating the performance of generators. The tool includes three aspects: frequency response, voltage response, and voltage schedule tracking.

Follum, Jim↗

Numerical study on aerosol sampling in a nuclear facility duct with a 90-degree elbow

Due to the challenging design requirements, elbows are often unavoidable in duct configuration, and these 90-degree bends introduce swirl, velocity variations, recirculation, and secondary flow. These disturbances make it difficult for nuclear facilities to meet particle sampling standards. A series of numerical analyses are conducted to track aerosols in a nuclear facility duct having a 90-degree elbow with the assistance of computational fluid dynamics (CFD). A turbulence model, a continuity, and a momentum, a discrete phase model, and species transport equations are solved simultaneously to track aerosols in the duct. The effect of turbulence models, turbulent dispersion models, droplet drag model, aerosol amount, aerosol spray configuration, guide vanes, and mixers are investigated. Simulation results are analyzed per relevant testing codes such as DOE-HDBK-1169, ASME AG1, ISO 14644-3, ACGIH, and ANSI/HPS N13.1.

Han, Kai [Savannah River Nuclear Solutions (SRNS),↗

Opening the reheating box in multifield inflation

The robustness of multi-field inflation to the physics of reheating is investigated. In order to carry out this study, reheating is described in detail by means of a formalism which tracks the evolution of scalar fields and perfect fluids in interaction (the inflatons and their decay products). This framework is then used to establish the general equations of motion of the background and perturbative quantities controlling the evolution of the system during reheating. Next, these equations are solved exactly by means of a new numerical code. Moreover, new analytical techniques, allowing us to interpret and approximate these solutions, are developed. As an illustration of a physical prediction that could be affected by the micro-physics of reheating, the amplitude of non-adiabatic perturbations in double inflation is considered. It is found that ignoring the fine-structure of reheating, as usually done in the standard approach, can lead to differences as big as ~ 50%, while our semi-analytic estimates can reduce this error to ~ 10%. Finally, we conclude that, in multi-field inflation, tracking the perturbations through the details of the reheating process is important and, to achieve good precision, requires the use of numerical calculations.

79 ASTRONOMY AND ASTROPHYSICS↗

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Simulation of electron Bernstein waves using FullWave with a 2D non-local hot plasma model

Hot plasma wave simulation capability is expanded in the FullWave code by updating the hybrid iterative solver in the code with a semi-implicit time stepping method. The new approach is used to simulate Electron Bernstein Wave (EBW) heating in over-dense spherical tokamak plasmas. The code’s hybrid iterative solver circumvents the prohibitive memory cost of direct methods by combining a time evolution of Maxwell’s equations with frequency-domain relaxation, while the conductivity kernel, calculated via 3D particle tracking, captures the essential non-local wave–particle interactions. One-dimensional EBW simulations verify the algorithm’s accuracy by demonstrating mode conversion from X-mode wave to EBW at the upper hybrid resonance and a strong cyclotron damping near the plasma core. Two-dimensional simulation reproduces the predicted short EBW wavelength and quantitatively matches the hot-plasma dispersion relation. This study demonstrates the fidelity of the hybrid solver for the electron cyclotron frequency range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of the RF system with reversed phasing

In this paper we report on the ELEGANT simulation for the Electron Storage Ring (ESR) RF system with reversed phasing using the lattice version of 5.3 at 5 GeV energy. The full beam dynamics results, including energy spread, bunch length and centroid offset as a function of bunch number in the train are presented. Some of the RF system related parameters, required as the input for particle tracking simulations, have been calculated analytically using the RF system related equations. The obtained results has been compared with the results simulated by Tianmu Xin, Michael Blaskiewicz and Gabriele Bassi using C++ and the SPACE codes.

43 PARTICLE ACCELERATORS↗

Overview of System Integration Analysis Activities for Integrated Waste Management

Spent nuclear fuel (SNF) generated by the current fleet of commercial nuclear reactors is being stored at the reactor sites in spent fuel pools (SFPs) and in dry independent spent fuel storage installations (ISFSIs). The US Department of Energy Office of Nuclear Energy (DOE-NE) is developing an Integrated Waste Management Program (IWMP) comprising a suite of options and supporting analyses to enable future informed choices. The IWMP is organized into the following five major areas: 1) Consent-Based Siting, 2) IWM Facilities and Equipment Concepts and Development, 3) Transportation Capability Analysis and Support, 4) Information Technology Solutions and Support, and 5) System Integration Analysis and Support. This paper discusses the activities ongoing in the IWMP System Integration Analysis and Support area. Two main areas of research in system integration are data and tools development, as well as system analysis assessments. One of the tools being developed in the system integration area is the Used Nuclear Fuel-Storage, Transportation & Disposal Analysis Resource and Data System (UNF-ST&DARDS) tool. It is being developed as a foundational resource for DOE-NE to manage SNF data, along with several compatible analysis tools for time-dependent characterization of SNF and related systems. UNF-ST&DARDS has the unparalleled ability to track SNF through the entire back end of the fuel cycle—from the time the fuel is discharged from a reactor through its disposal in a geological repository. UNF ST&DARDS interfaces with the SCALE code system for nuclear analysis and COBRA-SFS for thermal analysis. Another main tool being developed is the Next Generation System Analysis Model (NGSAM). NGSAM is an agent-based simulation software tool expressly designed to be capable of modeling the waste management system. NGSAM has been developed to enable informed decision-making by providing the capability of analyzing various potential system options for the management of SNF and HLW. Using NGSAM, system architecture analyses are being conducted to support the future deployment of a comprehensive nuclear waste management system that considers all major back-end aspects of the nuclear fuel cycle (i.e., transportation, storage, and disposal). System analysis assessments may investigate the implications of various strategies such as different acceptance rates, acceptance queues, facility capacities and options, standardized canisters, and different assumed system operation start dates. Recently, some system analysis effort has begun to look at how the waste management system might operate for advanced reactor fuel cycles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Extension of OpenMC for Fixed Source Transmutation Calculations

This report documents work performed under a Strategic Partnership Project between Argonne National Laboratory (ANL) and the United Kingdom Atomic Energy Authority (UKAEA). The overall goal of this project is to extend OpenMC [1], a community-developed Monte Carlo particle transport code, to be able to perform fixed-source transmutation calculations. In fission and fusion reactors, the high flux of energetic neutrons causes materials within the reactor to “transmute,” or undergo a nuclear reaction that results in the addition/removal of neutrons and protons from the nucleus of an atom. If subjected to these reactions for long enough, the overall composition and physical properties of the material itself begin to change as a result of transmutation. Such a feature is vital for predicting the decrease in tritium production rate within a breeder blanket during the lifetime of a fusion reactor. OpenMC is capable of simulating neutron transport in fission/fusion systems, thereby allowing it to estimate the flux that causes transmutation. It is also capable of solving the transmutation equations, which determine how the composition of a material changes over time due to neutron irradiation and radioactive decay. However, solving the transmutation equations was previously only possible when the source of neutrons came from a fission system. In a fusion system, the source of neutrons is typically determined by a separate code and then given as an input to the particle transport simulation. This is known as a fixed source calculation. Through this project, we have extended OpenMC to solve the transmutation equations for a fixed source calculation. Evaluating the change in material compositions due to transmutation and its effect on physical properties is of key importance to a range of engineering analyses for fission and fusion systems. For example, in a fusion reactor, estimating the dose rate at different physical locations resulting from irradiated materials in the reactor allows designers to ensure that workers are not exposed to doses beyond applicable regulations. In order to properly dispose of irradiated materials, designers also need to estimate the radiotoxicity, which again relies on knowledge of the material composition at some future time. The specific tasks for this project that were agreed to between ANL and UKAEA were as follows: 1. Make changes and additions in the openmc.deplete and related modules in OpenMC to support transmutation calculations following a fixed source transport simulation. 2. Make necessary changes to OpenMC to model transmutation due to an arbitrary set of reactions needed for fusion applications. Use this new capability to generate a depletion chain file based on the TENDL nuclear data library. 3. Improve the openmc.deplete module in OpenMC to keep track of gases produced as a by-product of nuclear reactions during transmutation calculations. 4. Validate the new capabilities by carrying out fixed-source transmutation calculations on a suitable benchmark problem using OpenMC and a comparable Monte Carlo neutron transport code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dust Formation and Destruction in Core Collapse Supernovae Explosions [Dissertation]

Dust is a crucial component of the Interstellar Medium (ISM) and an important driver in many astrophysical processes. Dust formed from evolved stars increases the metallicity of the ISM, enriching it with heavier elements. Dust makes up the source material of collapsing clouds in the ISM leading to stellar and galactic formation. The stabilization of these systems are aided by dust through absorption and scattering of excess energy. As stellar systems evolve and undergo supernova and mass loss events, enriched dust is produced, further enhancing the metallicity of the ISM. However, despite the importance of dust, the processes that produce dust are not well understood. Core-Collapse Supernovae are one source of dust. Dying massive stars explode, sending their constituent material into the Interstellar Medium, which as the ejecta expands and cools, material condenses and nucleates forming dust grains. I investigate the formation, survivability, and properties of dust. In Chapter 2, dust grain nucleation in CCSNe is discussed. A revised formulation of Kinetic Nucleation Theory (KNT) is used to track the nucleation and growth of dust grains in the ejecta. The affects of the progenitor system on the amount and type of dust is explored, showing a strong dependence on explosion energy. I discuss the code and methods used to model dust nucleation and destruction in Chapter 3. The sputtering and erosion of dust grains in CCSNe shocks are explored in Chapter 4. The erosion and survival of this nucleated dust is modeled through thermal and non-thermal sputtering by the surrounding gas and shocks. The optical properties and scattering of light by surviving dust grains is modeled by Mie scattering in Chapter 5. A discussion of observations of dust in variable star systems and supernovae is outlined in Chapter 6.

79 ASTRONOMY AND ASTROPHYSICS↗