Search NASA⌕ Search

SEARCH · Search NASA

Results for “time coding”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Report on technical development in ANL for innovative creep testing and measurement

This report summarizes a research effort undertaken by Argonne National Laboratory to accelerate the development of creep data package for ASME code qualification. Progress on the fabrication of a 3-sample load train and the application of a direct current potential drop (DCPD) technique to measure large creep strains are reported. Specifically: (1) The 3-sample load train, with three sets of Linear Variable Differential Transformers (LVDT), has been fabricated and assembled. (2) A DCPD method for creep-strain measurement up to 50% has been demonstrated. The DCPD data matches very well to LVDT results, suggesting its potential application in future creep testing. The technical developments contribute towards shortening the time needed to generate a creep data package for an ASME code case of the laser powder bed fusion 316H stainless steel, an effort under the US Department of Energy’s Advanced Materials and Manufacturing Technologies (AMMT) program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

FlowPM: Distributed TensorFlow implementation of the FastPM cosmological N-body solver

Here, we present FlowPM, a Particle-Mesh (PM) cosmological N-body code implemented in Mesh-TensorFlow for GPU-accelerated, distributed, and differentiable simulations. We implement and validate the accuracy of a novel multi-grid scheme based on multiresolution pyramids to compute large-scale forces efficiently on distributed platforms. We explore the scaling of the simulation on large-scale supercomputers and compare it with corresponding Python based PM code, finding on an average 10x speed-up in terms of wallclock time. We also demonstrate how this novel tool can be used for efficiently solving large scale cosmological inference problems, in particular reconstruction of cosmological fields in a forward model Bayesian framework with hybrid PM and neural network forward model. We provide skeleton code for these examples and the entire code is publicly available at https://github.com/modichirag/flowpm[Formula presented].

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗
DOE OSTI2023

Impact of hohlraum cooling on ignition metrics for inertial fusion implosions

This paper extends the evaluation of ignition metrics to include the impact of hohlraum cooling before peak implosion velocity in radiation driven implosions. First, we provide an extension of the results for the key hot spot stagnation quantities from the 2018 paper. The modified analytic expressions presented here match the Hydra results for these National Ignition Facility scale implosions both with and without hohlraum cooling before peak velocity if the effective ablation pressure P abl (effective) = P abl (t pv - 0.5 ns) is used in the analytic formulas, where t pv is the time of peak implosion velocity. Second, we provide an analysis that enables a comparison of the Hydra radiation hydrodynamics code calculations utilized here with the predictions of the analytic piston model of an ICF implosion, which focused on sensitivity to time duration of the hohlraum cooling phase before peak velocity (often called the “coast time”) and the shell radius at peak velocity R pv . Third, we provide a set of ignition metrics that are valid across a wide range of capsule designs valid for implosions both with and without hohlraum cooling before peak implosion velocity is reached.

42 ENGINEERING↗
DOE OSTI2023

OpenACC Unified Programming Environment for Multi-hybrid Acceleration with GPU and FPGA

Accelerated computing in HPC such as with GPU, plays a central role in HPC nowadays. However, in some complicated applications with partially different performance behavior is hard to solve with a single type of accelerator where GPU is not the perfect solution in these cases. We are developing a framework and transpiler allowing the users to program the codes with a single notation of OpenACC to be compiled for multi-hybrid accelerators, named MHOAT (Multi-Hybrid OpenACC Translator) for HPC applications. MHOAT parses the original code with directives to identify the target accelerating devices, currently supporting NVIDIA GPU and Intel FPGA, dispatching these specific partial codes to background compilers such as NVIDIA HPC SDK for GPU and OpenARC research compiler for FPGA, then assembles binaries for the final object with FPGA bitstream file. In this paper, we present the concept, design, implementation, and performance evaluation of a practical astrophysics simulation code where we successfully enhanced the performance up to 10 times faster than the GPU-only solution.

Boku, Taisuke↗
DOE OSTI2021

Community-scale interaction of energy efficiency and demand flexibility in residential buildings

Demand-side management (DSM) strategies, including energy efficiency (EE) and demand flexibility (DF), contribute to cost-effective operation of the electricity grid. From a system-level perspective, such programs reduce costs, enhance reliability, and reduce network issues. Similarly, DSM programs help participating customers reduce utility bills while maintaining occupant comfort. Understanding the relationship between EE and DF is key to realizing the full potential of DSM programs. In this study, we modeled an all-electric residential community based on a 498-home community that is planned for construction in Fort Collins, Colorado in the United States. We used this community model to study the relationship between different EE measures, including building envelope upgrades and smart appliances, and DF enabled by a home energy management system (HEMS) responding to a time-varying tariff. Various EE levels in the homes – code-minimum, zero energy ready, and even higher levels of envelope efficiency – were simulated. DF is enabled by the HEMS, which coordinates behind-the-meter resources, including flexible building loads, PV, and home battery systems, to minimize utility bills while maintaining occupant comfort. When comparing to the code-minimum homes, EE upgrades alone reduce HVAC energy use during peak hours by up to 50% and the HVAC utility bill by up to $312/year. With the addition of HEMS, the average daily peak demand can be reduced by up to 0.58 MW or 1.2 kW/home in the higher envelope efficiency homes. The combination of EE upgrades, HEMS, and home battery systems is expected to save homeowners up to $590/year while increasing community load flexibility. However, HEMS and home battery systems are less effective in increasing the DF in the more efficient homes due to the lower load.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
DOE OSTI2025

GPR_calculator: An on-the-fly surrogate model to accelerate massive nudged elastic band calculations

We present GPR_calculator, a package based on Python and C++ programming languages to build an on-the-fly surrogate model using Gaussian Process Regression (GPR) to approximate computationally expensive electronic structure calculations. The key idea is to dynamically train a GPR model during the simulation that can accurately predict energies and forces with uncertainty quantification. When the uncertainty is high, the costly electronic structure calculation is performed to obtain the ground truth data, which is then used to update the GPR model. To illustrate the effectiveness of GPR_calculator, we demonstrate its application in Nudged Elastic Band (NEB) simulations of surface diffusion and reactions, achieving 3-10 times acceleration compared to pure ab initio calculations. The source code is available at https://github.com/MaterSim/GPR_calculator.

Gaussian process regression↗
DOE OSTI2025

Self-consistent equilibrium and transport simulations for NSTX-U plasmas enhanced via machine learning surrogate models

The Control-Oriented Transport SIMulator (COTSIM) is an advanced equilibrium and transport code designed for simulating tokamak discharges at computational speeds suitable for control applications. COTSIM’s modular framework enables users to select models that balance accuracy with speed according to specific needs, allowing the code to operate from fast to faster-than-real-time performance levels. This work presents recent enhancements to COTSIM’s predictive accuracy for NSTX-U scenarios, achieved by integrating neural-network-based surrogate models and self-consistent equilibrium calculations. To improve source deposition predictions, a surrogate model for NUBEAM has been incorporated. Additionally, a surrogate model for the Multi-Mode Module (MMM) now supports predictions of anomalous thermal, momentum, and particle diffusivities—key factors for modeling the evolution of temperature and rotation. Each surrogate model was specifically trained for the NSTX-U operational regime to enhance COTSIM’s accuracy while maintaining computational efficiency. Moreover, COTSIM now couples fixed-boundary equilibrium solvers with its transport solvers, enabling self-consistent predictions of plasma profiles and equilibrium evolution over the discharge. Simulation results demonstrate strong agreement between COTSIM and TRANSP predictions for NSTX-U discharges. These substantial advancements expand COTSIM’s utility in model-based control applications for NSTX-U. Potential applications include simultaneous optimization of equilibrium and transport scenarios, integration into digital twins, real-time profile estimation (e.g., temperature and rotation) from limited or noisy measurements, and advanced feedback-based scenario control.

Equilibrium and transport modeling↗
DOE OSTI2021

Numerical simulation study of impurity B transport during real-time B powder injection in EAST

To understand the beneficial effects of impurity boron (B) on the plasma, the transport of B species in both atomic and ionic states in the edge plasma and the ensuing evolution of the background plasma are simulated during real-time B powder injection by using two-dimensional fluid code SOLPS-ITER. The initial background plasma profiles for the simulation are reconstructed based on an EAST upper-single null discharge (shot No. 70601). With the injection rate of B atoms, 4.0 × 10 21 atoms/s, the simulation shows that effective charge number is less than 1.1 at the inner boundary R - R sep = -5.2 cm (background plasma density, 5.1 × 10 19 m -3 ) at the outer mid-plane indicating that the B impurity does not have much influence on the energy confinement. The present work finds that the impurity B exhibits very specific features in the edge plasma, mostly present in the narrow region just outside the separatrix. The particular distribution pattern can reduce the heat flux greatly, leading to the divertor detachment. The divertor detachment can be achieved easily under an injection rate of B powder in a reasonable range. This study suggests that the impurity B has a promising potential application in future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Speeding Up Hartree–Fock in JuliaChem with Density Fitting

In this work, the density fitting (DF) approximation is added to the restricted Hartree–Fock (RHF) implementation in the JuliaChem computational chemistry code. Utilizing a DF algorithm that uses symmetry and integral screening, a significant reduction in time to compute the Fock matrix is achieved. The symmetry and screening DF-RHF techniques were adapted to be performed on graphics processing units (GPUs), which are well suited to perform the matrix multiplications that comprise the bulk of the Fock build time in DF-RHF. The JuliaChem DF-RHF GPU algorithm employs a novel approach that automatically switches between two DF-RHF algorithms depending on the number of basis functions in the calculation. The JuliaChem GPU DF-RHF implementation demonstrates up to 2× speedup for Fock build times compared to the existing best-in-class GPU DF-RHF implementation by operating directly on screened intermediate matrices. Due to the high portability of the Julia language code, the JuliaChem CPU and GPU DF-RHF implementations could be benchmarked on a variety of CPU and GPU architectures from multiple hardware vendors.

Hayes, John J. [Ames Laboratory, and Iowa State Un↗
DOE OSTI2021

Tokamak ITG-KBM transition benchmarking with the mixed variables/pullback transformation electromagnetic gyrokinetic scheme

Electromagnetic gyrokinetic simulation of high temperature plasma is required to predict confinement in magnetic fusion devices and has posed challenges for existing codes. In this paper, we demonstrate successful global gyrokinetic simulation of the ion temperature gradient-driven mode-kinetic ballooning mode transition in a toroidal fusion plasma test case using the mixed variables/pullback transformation (MV/PT) scheme with the particle-in-cell codes XGC and ORB5, and compare to results from a conventional continuum code from the literature. Furthermore, the MV/PT scheme combines explicit time integration with mitigation of the well-known electromagnetic gyrokinetic “cancelation problem.” We calculate eigenmodes in the electrostatic and parallel vector potentials, and find good agreement in growth rate, real frequency, and the normalized plasma pressure of mode transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

High-Throughput Computing: Case Study of Medical Image Processing Applications

HPC is designed for large-scale simulations using monolithic codes of tightly coupled processes highly optimized to deliver decreased time to solution. Medical image processing is not a traditional field of HPC. Similar to AI applications, medical image processing parses large datasets, typically multiple times, to support a variety of studies for classification, diagnosis or monitoring purposes. The convergence of AI, HPC and Big Data encouraged more fields using image processing to transition to HPC. However, not all applications benefit from the same optimizations. In this paper we focus on high throughput medical image processing applications that analyze a huge dataset of small MRI images and that require HPC systems to decrease the time of parsing the entire dataset and not individual MRIs. We show in this research the performance of running SLANT, an image processing application for a whole brain segmentation, on large-scale systems and highlight performance limitations. We present optimizations prioritizing throughput that exhibit a 3.5x speed-up on the Summit Supercomputer that can be used as a baseline for building a high-throughput execution framework for other HPC systems.

Predescu, Maria↗
DOE OSTI2024

Annual IC Progress Report

The goal of this project is to computationally investigate the progenitors of astrophysical gamma-ray bursts (GRB), the most extreme explosions in the universe across length and time scales. We have developed a cutting-edge, code-bridging approach to exploring the physics behind these luminous events to test the hypothesis that GRBs with electromagnetic radio emission originate from massive star (MS) and black hole (BH) binary systems. Our multi-scale, multi-physics theoretical framework for GRBs is based on solving well-defined and testable hydrodynamics problems associated with the dynamical evolution of the system. We solve stages of the evolution in a tractable way with a suite of well-developed state-of-the-art computational tools that include the stellar evolution code MESA, the general relativistic magnetohydrodynamics (GRMHD) code Athena++, and the binary black hole (BBH) population synthesis code COSMIC, as well as a number of self-written python post-processing tools.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Computing the Instantaneous Collision Probability between Satellites using Characteristic Function Inversion

The probability that two satellites overlap in space at a specified instant of time is called their instantaneous collision probability. Assuming Gaussian uncertainties and spherical satellites, this probability is the integral of a Gaussian distribution over a sphere. This paper shows how to compute the probability using an established numerical procedure called characteristic function inversion. The collision probability in the short-term encounter scenario is also evaluated with this approach, where the instant at which the probability is computed is the time of closest approach between the objects. Python and R code is provided to evaluate the probability in practice. Overall, the approach has been established for over fifty years, is implemented in existing software, does not rely on analytical approximations, and can be used to evaluate two and three dimensional collision probabilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Demonstration of Monte Carlo Reference Equilibrium Core Solutions for Pebble Bed Reactors Using Shift in kugelpy

This report describes how to produce reference solutions for the equilibrium core search of pebble bed reactors using the pyrates Python library. The pyrates library uses the kugelpy methodology to perform running-in calculations to reach an equilibrium core. This approach relies on full-core Monte Carlo neutron transport calculations and the shifting of fuel pebble inventory through the reactor to simulate pebble movement in the core from one step to the next of a running-in scenario. In addition to the use of Serpent for the Monte Carlo neutron transport calculations, as part of this work, updates for the use of the Shift Monte Carlo code were integrated into the official pyrates GitHub repository. Comparisons of running-in simulations for a generic pebble bed reactor produced by using pyrates with Serpent and Shift are shown. Although consistent trends are shown between the use of the two Monte Carlo codes, the limiting factor in either calculation is the computation time due to the need to perform several hundred full-core neutron transport calculations before an equilibrium core is reached.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Demonstration of Monte Carlo Reference Equilibrium Core Solutions for Pebble Bed Reactors Using Shift in kugelpy

This report describes how to produce reference solutions for the equilibrium core search of pebble bed reactors using the pyrates Python library. The pyrates library uses the kugelpy methodology to perform running-in calculations to reach an equilibrium core. This approach relies on full-core Monte Carlo neutron transport calculations and the shifting of fuel pebble inventory through the reactor to simulate pebble movement in the core from one step to the next of a running-in scenario. In addition to the use of Serpent for the Monte Carlo neutron transport calculations, as part of this work, updates for the use of the Shift Monte Carlo code were integrated into the official pyrates GitHub repository. Comparisons of running-in simulations for a generic pebble bed reactor produced by using pyrates with Serpent and Shift are shown. Although consistent trends are shown between the use of the two Monte Carlo codes, the limiting factor in either calculation is the computation time due to the need to perform several hundred full-core neutron transport calculations before an equilibrium core is reached.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗