Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “temporal computing”

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

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗

Final report- UFL - RAPIDS2: A SciDAC Institute for Computer Science, Data, and Artificial Intelligence

The research initiatives supported by the U.S. Department of Energy (DOE) Grant DE-SC0022265 are fundamentally aimed at pioneering advanced machine learning (ML) techniques for scientific data compression within high-performance computing (HPC) environments. This comprehensive body of work addresses the critical challenge posed by the exponential growth of data generated by scientific simulations in domains such as fusion energy, climate modeling, and computational fluid dynamics (CFD). A core objective is to develop compression algorithms that achieve substantial data reduction—often by orders of magnitude—while rigorously ensuring the fidelity of both the primary data (PD) and scientifically crucial derived quantities of interest (QoI). The methodologies deployed under this grant integrate sophisticated deep learning architectures, prominently featuring autoencoders, advanced generative models like conditional diffusion, and hybrid learning techniques. Key innovations include the development of Guaranteed Autoencoders (GAE) and the Guaranteed Conditional Diffusion with Tensor Correction (GCDTC) framework, which provide explicit, instance-level error bounds on reconstructed data. Furthermore, specialized strategies such as nonlinear constraint satisfaction are employed to preserve the integrity of QoI, a vital requirement for the trustworthiness of downstream scientific analyses. This research also focuses on the design and implementation of scalable, GPU-accelerated software pipelines that seamlessly integrate into existing HPC workflows, ensuring both computational efficiency and practical applicability. The CAESAR framework, for example, unifies foundation and generative models to create an adaptive and efficient compression solution for spatio-temporal scientific data. Collectively, these efforts represent a significant advancement in mitigating the scientific data deluge, enabling more effective data management, accelerated scientific discovery, and optimized utilization of HPC resources.

97 MATHEMATICS AND COMPUTING↗

HPC Molecular Simulation Tries Out a New GPU: Experiences on Early AMD Test Systems for the Frontier Supercomputer

Molecular simulation is an important tool for nu- merous efforts in physics, chemistry, and the biological sciences. Simulating molecular dynamics requires extremely rapid cal- culations to enable sufficient sampling of simulated temporal molecular processes. The Hewlett Packard Enterprise (HPE) Cray EX Frontier supercomputer installed at the Oak Ridge Leadership Computing Facility (OLCF) will provide an exascale resource for open science, and will feature graphics processing units (GPUs) from Advanced Micro Devices (AMD). The future LUMI supercomputer in Finland will be based on an HPE Cray EX platform as well. Here we test the ports of several widely used molecular dynamics packages that have each made substantial use of acceleration with NVIDIA GPUs, on Spock, the early Cray pre-Frontier testbed system at the OLCF which employs AMD GPUs. These programs are used extensively in industry for pharmaceutical and materials research, as well as academia, and are also frequently deployed on high-performance computing (HPC) systems, including national leadership HPC resources. We find that in general, performance is competitive and installation is straightforward, even at these early stages in a new GPU ecosystem. Our experiences point to an expanding arena for GPU vendors in HPC for molecular simulation.

Sedova, Ada↗

Advancement of Functional Expansion Tallies Capabilities in Serpent 2

The need for accident tolerant fuels (ATFs) became self-evident after the series of unfortunate nuclear events at Fukushima Daiichi in 2011. Restarting the Transient Reactor Test Facility (TREAT) at INL is one component in this concerted effort to establish the ATF development and testing program under the United States Department of Energy. Prior to the closure of TREAT in 1994, each experiment required multiple low-power calibration tests before and after each run. Since then, computing advancements have set computer-based modeling and simulation (M and S) as the preferred calibration method. The anticipated savings are both temporal and monetary. MAMMOTH, a collection of M and S tools built from the finite element-based MOOSE Framework, has been established to support the designed experiments intended for TREAT. This collection of tools includes: 1) Rattlesnake, a striking diffusion-based neutronics code, 2) BISON, an advanced materials and fuels performance code, and 3) RELAP-7, a powerful thermal-hydraulics solver. Each is capable of independently performing its function accurately and effectively; unified, they create a powerful multiphysics reactor core M and S application. One unfortunate weakness in this configuration lies in the unique geometry and design of TREAT's core. Although simple to model, this poses a challenge to the diffusion solver in Rattlesnake. One alternative is to replace Rattlesnake with Serpent, a fully-fledged Monte Carlo (MC) neutronics code. Coincidentally, Rattlesnake calculations are frequently validated against Serpent simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

gdess: A framework for evaluating simulated atmospheric CO 2 in Earth System Models

Atmospheric carbon dioxide (CO 2 ) plays a key role in the global carbon cycle and global warming. Climate-carbon feedbacks are often studied and estimated using Earth System Models (ESMs), which couple together multiple model components—including the atmosphere, ocean, terrestrial biosphere, and cryosphere—to jointly simulate mass and energy exchanges within and between these components. Despite tremendous advances, model intercomparisons and benchmarking are aspects of ESMs that warrant further improvement (Fer et al., 2021; Smith et al., 2014). Such benchmarking is critical because comparing the value of state variables in these simulations against observed values provides evidence for appropriately refining model components; moreover, researchers can learn much about Earth system dynamics in the process (Randall et al., 2019). We introduce `gdess` (a.k.a., Greenhouse gas Diagnostics for Earth System Simulations), which parses observational datasets and ESM simulation output, combines them to be in a consistent structure, computes statistical metrics, and generates diagnostic visualizations. In its current incarnation, `gdess` facilitates evaluating a model's ability to reproduce observed temporal and spatial variations of atmospheric CO 2 . The diagnostics implemented modularly in `gdess` support more rapid assessment and improvement of model-simulated global CO 2 sources and sinks associated with land and ocean ecosystem processes. We intend for this set of automated diagnostics to form an extensible, open source framework for future comparisons of simulated and observed concentrations of various greenhouse gases across Earth system models.

97 MATHEMATICS AND COMPUTING↗

An Optimal Pufferfish Privacy Mechanism for Temporally Correlated Trajectories

Temporally correlated trajectories are ubiquitous, and it has been a challenging problem to protect the temporal correlation from being used against users’ privacy. In this work, we propose an optimal Pufferfish privacy mechanism to achieve better data utility while providing guaranteed privacy of temporally correlated daily trajectories. First, a Laplace noise mechanism is realized through geometric sum of noisy Fourier coefficients of temporally correlated daily trajectories. Then, we prove that the proposed noisy Fourier coefficients’ geometric sum satisfies Pufferfish privacy, i.e. , the so-called FGS-Pufferfish privacy mechanism. Furthermore, we achieve better data utility for a given privacy budget by solving a constrained optimization problem of the noisy Fourier coefficients via the Lagrange multiplier method. What is more, a rigorous mathematical formula has been obtained for the Fourier coefficients’ Laplace noise scale parameters. At last, we evaluate our FGS-Pufferfish privacy mechanism on both simulated and real-life data and find that our proposed mechanism achieves better data utility and privacy compared with the other state-of-the-art existing approach.

97 MATHEMATICS AND COMPUTING↗

Sharing leaky-integrate-and-fire neurons for memory-efficient spiking neural networks

Spiking Neural Networks (SNNs) have gained increasing attention as energy-efficient neural networks owing to their binary and asynchronous computation. However, their non-linear activation, that is Leaky-Integrate-and-Fire (LIF) neuron, requires additional memory to store a membrane voltage to capture the temporal dynamics of spikes. Although the required memory cost for LIF neurons significantly increases as the input dimension goes larger, a technique to reduce memory for LIF neurons has not been explored so far. To address this, we propose a simple and effective solution, EfficientLIF-Net, which shares the LIF neurons across different layers and channels. Our EfficientLIF-Net achieves comparable accuracy with the standard SNNs while bringing up to ~4.3× forward memory efficiency and ~21.9× backward memory efficiency for LIF neurons. We conduct experiments on various datasets including CIFAR10, CIFAR100, TinyImageNet, ImageNet-100, and N-Caltech101. Furthermore, we show that our approach also offers advantages on Human Activity Recognition (HAR) datasets, which heavily rely on temporal information. The code has been released at https://github.com/Intelligent-Computing-Lab-Yale/EfficientLIF-Net.

60 APPLIED LIFE SCIENCES↗

Fast-ion physics in SPARC

Potential loss of energetic ions including alphas and radio-frequency tail ions due to classical orbit effects and magnetohydrodynamic instabilities (MHD) are central physics issues in the design and experimental physics programme of the SPARC tokamak. The expected loss of fusion alpha power due to ripple-induced transport is computed for the SPARC tokamak design by the ASCOT and SPIRAL orbit-simulation codes, to assess the expected surface heating of plasma-facing components. We find good agreement between the ASCOT and SPIRAL simulation results not only in integrated quantities (fraction of alpha power loss) but also in the spatial, temporal and pitch-angle dependence of the losses. If the toroidal field (TF) coils are well-aligned, the SPARC edge ripple is small (0.15–0.30 %), the computed ripple-induced alpha power loss is small ( ~0.25% ) and the corresponding peak surface power density is acceptable ( 244 kW m -2 ). However, the ripple and ripple-induced losses increase strongly if the TF coils are assumed to suffer increasing magnitudes of misalignment. Surface heat loads may become problematic if the TF coil misalignment approaches the centimetre level. Ripple-induced losses of the energetic ion tail driven by ion cyclotron range of frequency (ICRF) heating are not expected to generate significant wall or limiter heating in the nominal SPARC plasma scenario. Because the expected classical fast-ion losses are small, SPARC will be able to observe and study fast-ion redistribution due to MHD including sawteeth and Alfvén eigenmodes (AEs). SPARC's parameter space for AE physics even at moderate Q is shown to reasonably overlap that of the demonstration power plant ARC (Sorbom et al., Fusion Engng Des., vol. 100, 2015, p. 378), and thus measurements of AE mode amplitude, spectrum and associated fast-ion transport in SPARC would provide relevant guidance about AE behaviour expected in ARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Segmented slant hole collimator for stationary cardiac SPECT: Monte Carlo simulations

Purpose: This work is a preliminary study of a stationary cardiac SPECT system. The goal of this research is to propose a stationary cardiac SPECT system using segmented slant-hole collimators and to perform computer simulations to test the feasibility. Compared to the rotational SPECT, a stationary system has a benefit of acquiring temporally consistent projections. The most challenging issue in building a stationary system is to provide sufficient projection view-angles. Methods: A GATE (GEANT4 application for tomographic emission) Monte Carlo model was developed to simulate a two-detector stationary cardiac SPECT that uses segmented slant-hole collimators. Each detector contains seven segmented slant-hole sections that slant to a common volume at the rotation center. Consequently, 14 view-angles over 180° were acquired without any gantry rotation. The NCAT phantom was used for data generation and a tailored maximum-likelihood expectation-maximization algorithm was used for image reconstruction. Effects of limited number of view-angles and data truncation were carefully evaluated in the paper. Results: Simulation results indicated that the proposed segmented slant-hole stationary cardiac SPECT system is able to acquire sufficient data for cardiac imaging without a loss of image quality, even when the uptakes in the liver and kidneys are high. Seven views are acquired simultaneously at each detector, leading to 5-fold sensitivity gain over the conventional dual-head system at the same total acquisition time, which in turn increases the signal-to-noise ratio by 19%. The segmented slant-hole SPECT system also showed a good performance in lesion detection. In our prototype system, a short hole-length was used to reduce the dead zone between neighboring collimator segments. The measured sensitivity gain is about 17-fold over the conventional dual-head system. Conclusions: The GATE Monte Carlo simulations confirm the feasibility of the proposed stationary cardiac SPECT system with segmented slant-hole collimators. The proposed collimator consists of combined parallel and slant holes, and the image on the detector is not reduced in size.

60 APPLIED LIFE SCIENCES↗

GCAM Regional Tuning: A framework to tune GCAM parameters

GCAM assumptions typically generate scenarios that are designed to be internally consistent and globally coherent. The gcamdata tool which facilitates the compilation of data sets and user assumptions is not well suited to tailoring to specific country or regional realities, sponsor requirements, or perform harmonization for model intercomparison needs. As described in this report, the GCAM Regional Tuning project develops a computational framework that enables users to adjust GCAM parameters, so model outputs match targeted outcomes at user-defined spatial, temporal, and sectoral resolutions. The framework integrates GCAM, gcamdata, and gcamwrapper with a set of flexible “tuning directives” and an iterative numerical solver. Users can define targets (e.g., technology shares in power generation, BEV uptake, sectoral service demands), select tuners that manipulate relevant GCAM parameters (e.g., share weights, cost adders, elasticities), and export tuned parameters as reusable GCAM XML inputs for future runs. We demonstrate the approach and document usage, diagnostics, and known limitations, and we outline potential future directions.

97 MATHEMATICS AND COMPUTING↗

Parallel Time Integration for Constrained Optimization

The number of transistors in an average processor continues to increase, but individual clock speeds have plateaued. Those transistors are instead going into additional cores, increasing the number of different things that a processor can do at once and placing an emphasis on parallel computation. Many problems in scientific computing follow a time-evolution model, and it can be difficult to solve such problems in parallel across the temporal domain. The Multi-Grid Reduction In Time (MGRIT) algorithm, developed at Lawrence Livermore National Laboratory (LLNL), solves differential equations with a method designed specifically to take advantage of extreme numbers of processors by parallelizing across time. The Tri-diagonal MGRIT (TriMGRIT) algorithm, also developed at LLNL, is a generalization of MGRIT which enables parallel-in-time solving of a greater number of problems. Constrained optimization problems, in particular, may be solved in parallel using TriMGRIT. These consist of choosing a control function such that an objective functional is minimized, constrained by a differential-equation. We consider two such problems: applying torque to a pendulum to bring it to a gentle stop and moving a crowd of people from one distribution into another. We also perform some miscellaneous theoretical and practical research, including investigating the use of a line-search subroutine to refine intermediate TriMGRIT results and preliminary work on strategies for choosing operators for TriMGRIT to use.

97 MATHEMATICS AND COMPUTING↗

Automated anomalous behaviour detection in soil radon gas prior to earthquakes using computational intelligence techniques

Here in this article, three computational intelligence (CI) models were developed to automatically detect anomalous behaviour in soil radon gas ( 222 Rn) time series data. Data were obtained at a fault line and analysed using three machine learning techniques with the aim at identifying anomalies in temporal radon data prompted by seismic events. Radon concentrations were modelled with corresponding meteorological and statistical parameters. This leads to the estimation of soil radon gas without and with meteorological parameters. The comparison between computed radon concentration and actual radon concentrations was used in finding radon anomaly based upon automated system. The anomaly in radon time series data could be considered due to noise or seismic activity. Findings of study show that under meticulously characterized environments, on exclusion of noise contribution, seismic activity is responsible for anomalous behaviour seen in radon time series data.

58 GEOSCIENCES↗

An interregional optimization approach for time series aggregation in continent-scale electricity system models

Modeling electric power systems with high shares of weather-dependent resources requires tradeoffs between temporal, spatial, and operational resolution. Many studies perform time series aggregation using clustering algorithms to reduce the temporal dimension, but when modeling continent-scale electricity systems that are large enough to contain multiple independent weather systems, this approach requires large numbers of representative periods to minimize errors in regional wind and solar capacity factors. Here, a new optimization-based approach for representative period selection and weighting is introduced that minimizes regional errors in average renewable capacity factors and electricity demand. The method delivers higher regional fidelity with fewer representative periods than alternative clustering methods when applied to wind, solar, and demand profiles for the contiguous United States. When representative periods are selected from multiple weather years, the optimized method reproduces regional averages with lower error than a complete 365-day time series from any single weather year. The method identifies only representative (as opposed to outlying) periods but can be combined with an iterative "stress period" identification approach to guide efficient decision-making considering both average and high-risk weather conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Computational Study of Variable Fuel-Air Ratio and Hydrogen Doping in a Rotating Detonation Engine

Rotating detonation engines (RDEs) hold great promise as a technology for powering aviation engines. The pressure gain that results from continuously traversing detonation waves enables high theoretical performance compared to deflagration-based combustion. However, the harsh operating conditions inside an RDE restrict the availability of suitable experimental diagnostics to study these systems in detail and the range of spatial and temporal scales render numerical simulations of these systems difficult. In this study, we have combined adaptive mesh refinement and a robust shock-capturing approach to study RDE systems computationally using the compressible, reactive solver PeleC. Multi-species transport along with compressible Navier-Stokes equations are solved in the model along with finite-rate chemistry. Complex geometries are represented using an embedded boundary method with second-order spatial accuracy and the system is evolved in time using a second-order Runge-Kutta method. We present studies of a methane-air RDE system that is doped with varying levels of hydrogen. Our analysis focusses on how fuel mixture and doping levels effect the detonation flame structure and the generation of multiple wave modes.

ADVANCED PROPULSION SYSTEMS↗

Use of volumetric features for temporal comparison of mass lesions in full field digital mammograms

Purpose: Temporal comparison of lesions might improve classification between benign and malignant lesions in full-field digital mammograms (FFDM). The authors compare the use of volumetric features for lesion classification, which are computed from dense tissue thickness maps, to the use of mammographic lesion area. Use of dense tissue thickness maps for lesion characterization is advantageous, since it results in lesion features that are invariant to acquisition parameters. Methods: The dataset used in the analysis consisted of 60 temporal mammogram pairs comprising 120 mediolateral oblique or craniocaudal views with a total of 65 lesions, of which 41 were benign and 24 malignant. The authors analyzed the performance of four volumetric features, area, and four other commonly used features obtained from temporal mammogram pairs, current mammograms, and prior mammograms. The authors evaluated the individual performance of all features and of different feature sets. The authors used linear discriminant analysis with leave-one-out cross validation to classify different feature sets. Results: Volumetric features from temporal mammogram pairs achieved the best individual performance, as measured by the area under the receiver operating characteristic curve (A{sub z} value). Volume change (A{sub z} = 0.88) achieved higher A{sub z} value than projected lesion area change (A{sub z} = 0.78) in the temporal comparison of lesions. Best performance was achieved with a set that consisted of a set of features extracted from the current exam combined with four volumetric features representing changes with respect to the prior mammogram (A{sub z} = 0.90). This was significantly better (p = 0.005) than the performance obtained using features from the current exam only (A{sub z} = 0.77). Conclusions: Volumetric features from temporal mammogram pairs combined with features from the single exam significantly improve discrimination of benign and malignant lesions in FFDM mammograms compared to using only single exam features. In the comparison with prior mammograms, use of volumetric change may lead to better performance than use of lesion area change.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

In-situ TEM investigation of void swelling in nickel under irradiation with analysis aided by computer vision

Understanding the stability of irradiation-induced voids in materials is important for engineering material's swelling behavior under irradiation. In-situ TEM offers a spatial and temporal resolution that is suitable for investigating the evolution of voids under irradiation. However, the in-situ videos have often been too large to be analyzed manually, leaving the valuable data underutilized. For this work, we developed a deep learning-based semantic segmentation model to consistently study the growth and shrinkage of voids in nickel under 1 MeV krypton ion irradiation at various temperatures from 525 °C to 650 °C. With a foil thickness near 100 nm and ion flux of 6.3 × 10 11 ions∙ cm –2 ∙s –1 , the pre-existing voids, which were created beforehand by irradiation at 600 °C to 0.5 dpa, shrank at low temperatures and grew at high temperatures under irradiation, where the transition occurred at 575 °C (~0.5 T M ). The observed stability transition provided new insight for the shrinkage mechanism of voids under irradiation. In addition, an annealing experiment on nickel, previously irradiated at 600 °C to 3 dpa, was performed sequentially at 650 °C to 720 °C to reveal the shrinkage rate of void as a function of temperature and void size. The advantage of combining computer vision and in-situ TEM to obtain comprehensive void evolution was demonstrated.

36 MATERIALS SCIENCE↗

High quantum efficiency photocathode simulation for the investigation of novel structured designs

A computer model in CST Studio Suite has been developed to evaluate several novel geometrically enhanced photocathode designs. This work was aimed at identifying a structure that would increase the total electron yield by a factor of two or greater in the 1–30 keV range. The modeling software was used to simulate the electric field and generate particle tracking for several potential structures. The final photocathode structure has been tailored to meet a set of detector performance requirements, namely, a spatial resolution of <40 μm and a temporal spread of 1–10 ps. We present the details of the geometrically enhanced photocathode model and resulting static field and electron emission characteristics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗