SEARCH · Engineering Papers
Results for “time discretization”
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.
A combined strategy for dynamic probabilistic risk assessment of fission battery designs using EMRALD and DEPM
Not Available
Single Electron Self-coherence and Its Wave/Particle Duality in the Electron Microscope
Intensities in high-resolution phase-contrast images from electron microscopes build up discretely in time by detecting single electrons. A wave description of pulse-like coherent-inelastic interaction of an electron with matter implies a time-dependent coexistence of coherent partial waves. Their superposition forms a wave package by phase decoherence of 0.5 - 1 radian with Heisenbergs energy uncertainty ΔE H = $\hbar$/2 Δt -1 matching the energy loss ΔE of a coherent-inelastic interaction and sets the interaction time Δt. In these circumstances, the product of Planck's constant and the speed of light hc is given by the product of the expression for temporal coherence λ 2 /Δλ and the energy loss ΔE. Experimentally, the self-coherence length was measured by detecting the energy-dependent localization of scattered, plane matter waves in surface proximity exploiting the Goos–Hänchen shift. Chromatic-aberration Cc-corrected electron microscopy on boron nitride (BN) proves that the coherent crystal illumination and phase contrast are lost if the self-coherence length shrinks below the size of the crystal unit cell at ΔE > 200 eV. Finally, in perspective, the interaction time of any matter wave compares with the lifetime of a virtual particle of any elemental interaction, suggesting the present concept of coherent-inelastic interactions of matter waves might be generalizable.
Validation of Numerical Models of the Offshore Wind Turbine From the Alpha Ventus Wind Farm Against Full-Scale Measurements Within OC5 Phase III
The main objective of the Offshore Code Comparison Collaboration Continuation, with Correlation (OC5) project is validation of aero-hydro-servo-elastic simulation tools for offshore wind turbines (OWTs) through comparison of simulated results to the response data of physical systems. Phase III of the OC5 project validates OWT models against the measurements recorded on a Senvion 5M wind turbine supported by the OWEC Quattropod from the alpha ventus offshore wind farm. The following operating conditions of the wind turbine were chosen for the validation: (1) idling below the cut-in wind speed, (2) rotor-nacelle assembly (RNA) rotation maneuver below the cut-in wind speed, (3) power production below and above the rated wind speed, and (4) shutdown. A number of validation load cases were defined based on these operating conditions. The following measurements were used for validation: (1) strains and accelerations recorded on the support structure and (2) pitch, yaw, and azimuth angles, generator speed, and electrical power recorded from the RNA. Strains were not directly available from the majority of the OWT simulation tools; therefore, strains were calculated based on out-of-plane bending moments, axial forces, and cross-sectional properties of the structural members. The simulation results and measurements were compared in terms of time series, discrete Fourier transforms, power spectral densities, and probability density functions of strains and accelerometers. A good match was achieved between the measurements and models setup by OC5 Phase III participants.
Controlling quantum many-body dynamics in driven Rydberg atom arrays
Dynamic stabilization of an array Large-scale systems comprising one-dimensional chains and two-dimensional arrays of excited atoms held in a programmable optical lattice are a powerful platform with which to simulate emergent phenomena. Bluvstein et al. built an array of up to 200 Rydberg atoms and subjected the system to periodic excitation. Under such driven excitation, they found that the array of atoms stabilized, freezing periodically into what looked like time crystals. Understanding and controlling the dynamic interactions in quantum many-body systems lies at the heart of contemporary condensed matter physics and the exotic phenomena that can occur. Science , this issue p. 1355
Lattice Effective Field Theory Simulations of Nuclei
Lattice effective field theory applies the principles of effective field theory in a lattice framework where space and time are discretized. Nucleons are placed on the lattice sites, and the interactions are tuned to replicate the observed features of the nuclear force. Monte Carlo simulations are then employed to predict the properties of nuclear few- and many-body systems. Here, we review the basic methods and several theoretical and algorithmic advances that have been used to further our understanding of atomic nuclei.
Vitrification of Hanford Tank 241-AP-107 with Recycled Condensate
During the vitrification of nuclear waste at the Hanford Waste Treatment and Immobilization Plant (WTP) – the primary mission of the U.S. Department of Energy Office of River Protection – the offgas condensate generated from the waste-to-glass conversion is currently planned to be concentrated by evaporation in the Effluent Management Facility (EMF). This concentrated condensate can then be recycled back to the incoming waste and vitrified. To test the recycle process, a test apparatus was designed to mimic the EMF evaporator and used to concentrate a volume of condensate that had been previously produced during the vitrification of Hanford tank 241-AP-107 (referred to herein as AP-107) waste in a continuous laboratory-scale melter (CLSM). The test apparatus successfully concentrated the AP-107 condensate by a factor of 10 while retaining over 90 % of the technetium-99 ( 99 Tc), Cs, and I inventory. A second portion of AP-107 waste was retrieved by Washington River Protection Solutions, LLC, given to Pacific Northwest National Laboratory, and combined with the AP-107 condensate concentrate after undergoing solids filtration and cesium removal by ion exchange. This combination served to approximate the recycling action to be performed at the WTP. After the addition of glass-forming chemicals (GFCs), the combined AP-107 waste and AP-107 condensate concentrate were processed in the CLSM to produce a glass, called AP-107-1R, that was designed to satisfy the WTP baseline requirements (Kim et al. 2012). During the 8.87 hours of processing, 7.27 kg of AP-107-1R glass were produced for an average glass production rate of 1739 kg m 2 d -1 . Compared to the previous run in the CLSM without recycled condensate, the run with the recycle had a greater average glass production rate, but the rate was within the potential range of variability when processing melter feeds with similar composition in the CLSM. The glass produced from the AP-107 recycle run in the CLSM was within 10 % of the target AP-107-1R glass composition with respect to the primary glass components. Analysis of the minor component impurities revealed that their content in the glass product had approached their nominal target after 2 turnovers of the glass inventory in the CLSM while the activity of the minor radionuclides was retained in the glass product. The 99 Tc and total cesium content in the combined AP-107 waste and recycled condensate were maintained at concentrations expected to be experienced at the WTP. During processing in the CLSM, at discrete sampling time periods, the target 99 Tc/Cs mass ratio in the glass formulation varied from 0.9 to 62.9. Across this range, the Cs retention in the glass ranged from 53 to 60 %, while the retention from the entire runtime totaled 68 %, values which align with Cs retention in other scaled melter systems while processing LAW melter feeds at 99 Tc/Cs mass ratios ranging from 1 to 100. The 99 Tc retention in the glass ranged from 22 to 32 %, primarily due to the cold-cap coverage on the glass melt surface, the area covered by reacting melter feed, varying from ~80 % to ~95 % during processing, demonstrating greater volatility of 99 Tc from the glass while more surface was exposed, as expected based on previous 99 Tc volatility studies.
Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel
Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.
Code for the manuscript titled "Blackout Diffusion: Generative Diffusion Models in Discrete-State Spaces"
We would like to disclose two scripts, written in Jupyter notebook, in which we implement the "Blackout Diffusion Process" described in the manuscript "Blackout Diffusion: Generative Diffusion Models in Discrete-State Spaces" (LA-UR-23-20509), to be submitted to the International Conference of Machine Learning (ICML). The abstract of the manuscript is append below. == Typical generative diffusion models rely on a Gaussian diffusion process for training the backward transformations, which can then be used generate samples from Gaussian noise. However, real world data often takes place in discrete-state spaces, which occur in many scientific applications. Here, we develop a theoretical formulation for arbitrary discrete-state Markov processes in the forward diffusion process. We relate the theory to the existing continuous-state Gaussian diffusion and identify the corresponding reverse-time stochastic process and score function in the continuous-time setting, and the reverse-time mapping in the discrete-time setting. As an example of this framework, we introduce "Blackout Diffusion", which learns to produce samples from an empty image instead of from noise. Numerical experiments on the CIFAR-10 dataset confirm the feasibility of generative diffusion modeling in a discrete space. Generalizing from specific (Gaussian) forward processes to a more general framework also sheds light on how to interpret generative diffusion models and their mathematical structure, which we comment on.
An overview of continuous and discrete phasor analysis of binned or time-gated periodic decays
Time-resolved analysis of periodically excited luminescence decays by the phasor method in the presence of time gating or binning is revisited. Analytical expressions for discrete configurations of square gates are derived and the locus of the phasors of such modified periodic single-exponential decays is compared to the canonical universal semicircle. The effects of IRF offset, decay truncation and gate shape are also discussed. Finally, modified expressions for the phase and modulus lifetimes are provided for some simple cases. A discussion of a modified phasor calibration approach is presented.
Continuous and discrete phasor analysis of binned or time-gated periodic decays
The time-resolved analysis of periodically excited luminescence decays by the phasor method in the presence of time-gating or binning is revisited. Analytical expressions for discrete configurations of square gates are derived, and the locus of the phasors of such modified periodic single-exponential decays is compared to the canonical universal semicircle. The effects of instrument response function offset, decay truncation, and gate shape are also discussed. Finally, modified expressions for the phase and modulus lifetimes are provided for some simple cases. A discussion of a modified phasor calibration approach is presented, and an illustration of the new concepts with examples from the literature concludes this work.
Thermo-hydraulic steam pipe models for district heating simulations: Simplifications to balance accuracy and simulation speed
Steam piping networks are essential for optimizing performance in industrial processes and district heating systems. However, dynamic models that balance thermo-hydraulic accuracy with computational efficiency remain limited. In response, this paper presents a new discretized steam pipe model based on the plug flow approach, capturing key thermo-hydraulic behaviors while simplifying steam phase change processes. Implemented in Modelica, the model accurately calculates temperature and pressure distributions along steam pipelines. To improve computational efficiency for district-scale simulations, five model simplifications are introduced: lumped thermo-hydraulic functions, empirical correlations, fluid state approximations, steady-state dynamics and inclusion of flow derivatives. These simplified models achieve 85%-98% accuracy in predicting pressure drop and condensation losses, including dynamic condensate behavior during pipe warm-up—a factor often overlooked in existing models. The models support diverse network configurations, scaling effectively to systems with multiple distribution pipes and connected building loads. Discrete models provide detailed insights but exhibit a cubic increase in simulation time as the network scales by N connected building O(N 2.42 ). In contrast, lumped models simulate 10–28 times faster than discrete, offering quadratic scaling of simulation time O(N 1.73 ). However, they still require 6 times more computation time than a lossless network, highlighting the inherent computational challenges of modeling compressible fluid flow. In conclusion, the steady-state lumped variant, with its near-linear scalability in computational time O(N 1.01 ), emerges as an efficient solution for preliminary design evaluations and extensive parametric studies.
Simulation-based Performance Evaluation of Model Predictive Control for Building Energy Systems
The performance of model predictive control (MPC) can be significantly affected by different choices of controller parameters such as the time intervals for model discretization and control sampling. Due to the lack of a systematic understanding on how these parameters affect control performance, they are usually selected arbitrarily in practice.In this paper, the combined impacts of selected time intervals for model discretization and control sampling on the performance of MPC are comprehensively investigated for the first time through detailed simulations. Specifically, a typical MPC strategy is first designed to improve building operations based on a reduced-order model of building dynamics. Then, the performance of the designed MPC is evaluated against different choices of time intervals for model discretization and control sampling on a simulated office building. The detailed simulation results reveal that the time interval for model discretization has a much greater influence on the performance of MPC than the time interval for control sampling. Although the time interval for control sampling usually receives more attentions in practice, it turns out that the time interval for model discretization affects the prediction performance, cost saving, and computation time simultaneously and more significantly. Therefore, the simulation-based performance evaluation presented here sheds light on the impacts of different time intervals and facilitates their selection for practical applications of MPC to building operations
Discretization Writeup for Grey Flux-Limited Radiation Diffusion
This report documents the time and space discretizations for grey flux-limited diffusion applied to the thermal radiative transfer (TRT) equations. We begin with a description of the physics being solved before moving into the diffusion approximation. Once we have the TRT system, we show a finite-volume-inspired discretization from Jim Morel (Texas A&M University, NUEN 627 class notes, lecture 8). As systems become hotter they emit more photons in the form of blackbody radiation. Because average photon energy of the blackbody source is proportional to the temperature of the system, we call these thermal photons or thermal radiation. As material temperatures increase, increasing fractions of the total energy in the system go into the radiation field. In addition, radiation can deposit energy and momentum non-locally, making it an important phenomenon for heating and impulse. In order to accurately study systems at high temperatures, we wish to add the physics of thermal radiation to our hydrodynamic system. In practice, coupling radiation and hydrodynamics is often done by operator-splitting each timestep into two consecutive, non-overlapping phases: (1) update the hydrodynamics for a fixed radiation state (2) update the radiation and internal energy for an otherwise fixed hydrodynamic state. Because of this clean separation of physics updates, in this report we show only the latter phase, which involves solely the TRT equations.
Parallel-in-Time Methods for Method-of-Lines Discretizations of Nonlinear Hyperbolic PDEs and Systems (Final Report)
The work for the subcontract is situated in the area of parallel-in-time integration for hyperbolic partial differential equations (PDEs). Parallel-in-time integration is an active area of research due to its ability to enable faster numerical simulations for applications throughout many areas of science. The work in this subcontract builds on a variety of results that were obtained, as part of the work performed for Subcontract No. B648355, for the Multigrid Reduction-in-Time (MGRIT) method from [1] applied to hyperbolic PDEs. This subcontract extends these results further to more efficient methods and to the case of method-of-lines discretizations for nonlinear hyperbolic PDES and systems of PDEs. The following is a summary of the research performed and results achieved during milestone periods 1, 2 and 3 by the PI (Hans De Sterck) and Postdoctoral Research Associate (Oliver Krzysik), for required tasks 1-4 (as listed in the Statement of Work): Research over the previous year has been split into three main projects: (i) solution of acoustic equation system; (ii) solution of nonlinear scalar hyperbolic PDEs; (iii) solution of nonlinear hyperbolic systems of PDEs.
Improving ADAM through an implicit-explicit (IMEX) time-stepping approach
The ADAM optimizer, often used in machine learning for neural network training, corresponds to an underlying ordinary differential equation (ODE) in the limit of very small learning rates. Here, this work shows that the classical ADAM algorithm is a first-order implicit-explicit (IMEX) Euler discretization of the underlying ODE. Employing the time discretization point of view, we propose new extensions of the ADAM scheme obtained by using higher-order IMEX methods to solve the ODE. Based on this approach, we derive a new optimization algorithm for neural network training that performs better than classical ADAM on several regression and classification problems.
Real-Time Analysis of Stochastic Bitstreams from Discrete Tunnel Diodes
Explore the source record for details and available documents.
A Particle-in-cell Method for Plasmas with A Generalized Momentum Formulation, Part III: A family of Gauge Conserving Methods
In this paper, we introduce a new family of spatially co-located field solvers for particle-in-cell applications which evolve the potential formulation of Maxwell’s equations under the Lorenz gauge. Our recent work [2] introduced the concept of time-consistency, which connects charge conservation to the preservation of the gauge at the semi-discrete level. It will be shown that there exists a large family of time discretizations which satisfy this property. Additionally, it will be further shown that for large classes of time marching methods, the satisfaction of the gauge condition automatically implies the satisfaction of Gauss’s law for electricity, with the potential formulation ensuring that that Gauss’s law for magnetism is satisfied by definition. We focus on popular time marching methods including centered differences, backward differences, and diagonally-implicit Runge-Kutta methods, which are coupled to a spectral discretization in space. We demonstrate the theory by testing the methods on a relativistic Weibel instability and a drifting cloud of electrons.