SEARCH · Engineering Papers
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.
This paper gives results of a timing study for computing fast Fourier transforms (FFTs) on LANL based computing platforms. The platforms considered will be a node on Snow and a standard LANL HP laptop. The run-time for 1D FFTs will be measured for different sizes of data and the results are reported in FFTs/s.
Here, we study the problem of designing a distributed observer for an LTI system over a time-varying communication graph. The limited existing work on this topic imposes various restrictions either on the observation model or on the sequence of communication graphs. In contrast, we propose a single-time-scale distributed observer that works under mild assumptions. Specifically, our communication model only requires strong-connectivity to be preserved over non-overlapping, contiguous intervals that are even allowed to grow unbounded over time. We show that under suitable conditions that bound the growth of such intervals, joint observability is sufficient to track the state of any discrete-time LTI system exponentially fast, at any desired rate. We also develop a variant of our algorithm that is provably robust to worst-case adversarial attacks, provided the sequence of graphs is sufficiently connected over time. The key to our approach is the notion of a "freshness-index" that keeps track of the age-of-information being diffused across the network. Such indices enable nodes to reject stale estimates of the state, and, in turn, contribute to stability of the error dynamics.
Explore the source record for details and available documents.
Laser ablation in conjunction with optical emission spectroscopy is a potential non-contact, standoff detection method for all elements in the periodic table and certain isotopes such as radionuclides. Currently, significant development efforts are on-going to use ultrafast laser filaments for remote sensing. The application of filaments is of particular interest in extending the range of stand-off capability associated with elemental and isotopic detection via laser-induced breakdown spectroscopy. In this study, we characterize the expansion dynamics and chemical evolution of filament-produced uranium (U) plasmas. Laser filaments are generated in the laboratory by loosely focusing 35 femtosecond (fs), 6 milli Joule (mJ) pulses in air. Two-dimensional spectrally-integrated and time-resolved imaging was performed to study hydrodynamics and evolution of U atomic and UO molecular emission in filament-produced U plasmas. Our results highlight that filament ablation of U plasmas gives a cylindrical plasma plume morphology with an appearance of plume splitting into slow and fast moving components at later times of its evolution. Emission from the slow-moving region shows no distinct spectral features (i.e. broadband-like) and is contributed in part by nanoparticles generated during ultrafast laser ablation. Additionally, we find U atoms and U oxide molecules (i.e. UO, UxOy) co-exist in the filament produced plasma, which can be attributed to the generation of low-temperature plasma conditions during filament ablation
Three machine learning techniques (multilayer perceptron, random forest, and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modeling and real-time control applications. More accurate simulations with fast electron diffusion are even slower, requiring multiple hours of run time with parallel processing. The machine learning models use a database of 16,000+ GEN-RAY/CQL3D simulations for training, validation, and testing. Latin hypercube sampling methods implemented in πScope ensure that the database covers the range of 9 input parameters (n e0 , T e0 , I p , B t , R 0 , n ∥︀ , Z e f f , V loop , P LHCD ) with sufficient density in all regions of parameter space. The surrogate models reduce the computation time from minutes-hours to ms with high accuracy across the input parameter space. Data-driven surrogate models also allow for solving inverse and “lateral” problems. A surrogate model for the inverse problem maps from a desired current drive or power deposition profile to a set of input parameters that would result in such a profile, while a surrogate model for the lateral problem maps from a measured experimental quantity such as hard x-ray emission to a current drive or power deposition profile. In conclusion, the πScope database creation workflow is flexible and applicable to other RF simulation codes such as TORIC.
The OpenUniverse2024 simulation suite is a cross-collaboration effort to produce matched simulated imaging for multiple surveys as they would observe a common simulated sky. Both the simulated data and associated tools used to produce it are intended to uniquely enable a wide range of studies to maximize the science potential of the next generation of cosmological surveys. We have produced simulated imaging for approximately 70 deg 2 of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Wide-Fast-Deep survey and the Nancy Grace Roman Space Telescope High-Latitude Wide-Area Survey, as well as overlapping versions of the ELAIS-S1 Deep-Drilling Field for LSST and the High-Latitude Time-Domain Survey for Roman. OpenUniverse2024 includes (i) an early version of the updated extragalactic model called Diffsky, which substantially improves the realism of optical and infrared photometry of objects, compared to previous versions of these models; (ii) updated transient models that extend through the wavelength range probed by Roman and Rubin; and (iii) improved survey, telescope, and instrument realism based on up-to-date survey plans and known properties of the instruments. It is built on a new and updated suite of simulation tools that improves the ease of consistently simulating multiple observatories viewing the same sky. The approximately 400 TB of synthetic survey imaging and simulated universe catalogs are publicly available, and we preview some scientific uses of the simulations.
We report on the γ-ray emission above 100 MeV from the GOES M3.3 flare SOL2012-06-03. The hard X-ray (HXR) and microwave emissions have typical time profiles with a fast rise to a well-defined peak followed by a slower decay. The >100 MeV emission during the prompt phase displayed a double-peaked temporal structure with the first peak following the HXR and microwaves, and the second one, about 3 times stronger, occurring 17 ± 2 s later. The time profiles seem to indicate two separate acceleration mechanisms at work, where the second γ-ray peak reveals a potentially pure or at least largely dominant ion acceleration. The Atmospheric Imaging Assembly imaging shows a bright elliptical ribbon and a transient brightening in the northwestern (NW) region. Nonlinear force-free extrapolations at the time of the impulsive peaks show closed field lines connecting the NW region to the southeastern part of the ribbon, and the magnetic topology revealed clusters of nulls. These observations suggest a spine-and-fan geometry, and based on these observations, we interpret the second γ-ray peak as being due to the predominant acceleration of ions in a region with multiple null points. The >100 MeV emission from this flare also exhibits a delayed phase with an exponential decay of roughly 350 s. We find that the delayed emission is consistent with ions being trapped in a closed flux tube with gradual escape via their loss cone to the chromosphere.
Facilitating fast charging in lithium-ion batteries (LiBs) is often linked to Li-plating, which harms performance, longevity, and safety. Early detection of Li-plating is essential for rapid technological development and for preventing performance deterioration and ensuring safety during operation. Fast and real-time detection, using commonly collected measurements like voltage (V), current (I), temperature (T), and pressure (P), is highly desirable. Existing standalone methods relying on electrochemical and mechanical signatures, using half, smaller, or specially designed cells often operated at lower temperatures, fail to account for real-world fast-charging conditions. These signatures may also have inherent unreliability in aged LiBs, a phenomenon currently not-well understood. All these uncertainties have complicated practical implementation of a robust Li-plating detection technique. This study, through multiple case studies involving real-world fast-charging conditions using automotive-grade 11.6 Ah LiBs, shows that many single signal-based diagnostic techniques may be inadequate to detect Li-plating. Among various signatures, end-of-charge rest pressure, differential pressure-sensing, and end-of-charge rest voltage were identified as particularly useful in detecting Li-plating. Furthermore, A multi-signal-based detection technique is shown to be more robust in detecting Li-plating. Using both fresh and aged cells, the results and analysis highlight how detection capabilities are influenced by various factors, including battery design, size, charging speed, operating temperature, and degradation level. Adopting such a multi-signal Li-plating detection approach may be instrumental in the rapid development of battery technology in laboratories, as well as ensuring the enhanced safety of next-generation LiBs in real-world fast-charging applications.
We have used photon pair correlations generated via spontaneous parametric downconversion (SPDC) to measure the fluorescence lifetime of the organic dye rhodamine 6 G, demonstrating that fluorescence lifetime measurements can be achieved using a continuous wave (CW) laser, without pulsed or modulated lasers. Our entangled photon method, quantum fluorescence lifetime (Q-FL) measurements, uses one photon to excite fluorescence and the resulting fluorescence photon is timed and referenced to the arrival time of the other entangled photon. Thus, we can exploit the short timescale of photon pair correlations to conduct experiments that are typically carried out with pulsed lasers and we show that the inherent timing of the photons is fast enough to resolve the nanosecond scale fluorescence lifetime of the sample. This measurement paves the way towards using the time correlations of entangled photons for fluorescence imaging; capitalizing on the presence of fast, sub-100 ps correlations that have not been demonstrated classically.
Chemical looping is a near-zero emission process for generating power from coal. It is based on a multi-phase gas-solid flow and has extremely challenging nonlinear, multi-scale dynamics with jumps, producing large dynamic model uncertainty, which renders traditional robust control techniques, such as linear parameter varying H ∞ design, largely inapplicable. This process complexity is addressed in the present work through the temporal and the spatiotemporal multiresolution modeling along with the corresponding model-based control laws. Namely, the nonlinear autoregressive with exogenous input model structure, nonlinear in the wavelet basis, but linear in parameters, is used to identify the dominant temporal chemical looping process dynamics. The control inputs and the wavelet model parameters are calculated by optimizing a quadratic cost function using a gradient descent method. The respective identification and tracking error convergence of the proposed self-tuning identification and control schemes, the latter using the unconstrained generalized predictive control structure, is separately ascertained through the Lyapunov stability theorem. The rate constraint on the control signal in the temporal control law is then imposed and the control topology is augmented by an additional control loop with self-tuning deadbeat controller which uses the spatiotemporal wavelet riser dynamics representation. The novelty of this work is three-fold: (1) developing the self-tuning controller design methodology that consists in embedding the real-time tunable temporal highly nonlinear, but linearly parametrizable, multiresolution system representations into the classical rate-constrained generalized predictive quadratic optimal control structure, (2) augmenting the temporal multiresolution loop by a more complex spatiotemporal multiresolution self-tuning deadbeat control loop, and (3) demonstrating the effectiveness of the proposed methodology in producing fast recursive real-time algorithms for controlling highly uncertain nonlinear multiscale processes. The latter is shown through the data from the implemented temporal and augmented spatiotemporal solutions of a difficult chemical looping cold flow tracking control problem.
The Imaging Neutral Particle Analyzer (INPA) at DIII-D is a diagnostic system used to accurately resolve the energy and spatial distributions of fast ions in fusion plasmas. A novel artificial intelligence (AI) technique named INPA-net is based on Reservoir Computing Networks and developed here to predict active and passive signals produced by charge-exchange reactions from injected and edge-cold neutrals, respectively, in magnetically confined fusion plasmas. This model is trained using a set of 21 time domain signals between 0 s to 3.35 s that includes injected beam and thermal plasma information, and 6444 real 2D experimental images of the INPA in 12 plasma discharges at DIII-D. The trained neural network is able to forecast experimental images in real-time. The model achieves an R-squared value of 0.91, which is higher than the 0.83 value achieved by a simple linear regression model. This improvement highlights the model's enhanced predictive accuracy for measured images from the validation set. This AI approach is valuable due to its rapid response times and potential for integration into real-time plasma control systems. A version of this model capable of generating syntehic images would be useful for the real-time monitoring of fast-ion transport. A comprehensive sensitivity study reveals that INPA-net maintains high performance even with variations in the input parameters, indicating the model's robustness and reliability. While developed for the INPA, the underlying architecture is adaptable and may be applied to various 2D imaging diagnostics in fusion research.
Model-based optimization offers a systematic approach to advanced scenario planning. In this case, the feedforward-control inputs (actuator trajectories) that are needed to attain and sustain a desired scenario are obtained by solving a nonlinear constrained optimization problem. This class of problems generally minimize a cost function that measures the difference between desired and actual plasma states. Several numerical optimization algorithms, such as sequential quadratic programming, require repeated calculation of the cost function gradients with respect to the input trajectories. Calculating these gradients numerically can be computationally intensive, increasing the time needed to solve the feedforward-control optimization problem. Here, this work introduces a method to analytically calculate these cost function gradients from the current profile evolution model. This can significantly reduce the computational time and allow for fast feedforward-control optimization, which would eventually enable optimal scenario planning between discharges. The performance of the feedforward optimizer with analytical gradients is compared to a traditional optimization algorithm based on numerical gradients for different NSTX-U scenarios. The plasma dynamics in the optimization algorithm are simulated using the Control Oriented Transport SIMulator (COTSIM). Results of the work show that analytical gradients consistently reduce the computation time while achieving trajectories that are comparable to those obtained by traditional optimization algorithms based on numerical gradients.
Abstract Measuring a pulsar’s rotational evolution is crucial to understanding the nature of the pulsar. Here, we provide updated timing models for the rotational evolution of six pulsars, five of which are rotation phase-connected using primarily NICER X-ray data. For the newly discovered fast energetic young pulsar, PSR J0058−7218, we increase the baseline of its timing model from 1.4 days to 8 months and not only measure more precisely its spin-down rate ν ̇ = ( − 6.2324 ± 0.0001 ) × 10 − 11 Hz s − 1 but also for the first time the second time derivative of its spin rate ν ̈ = ( 4.2 ± 0.2 ) × 10 − 21 Hz s − 2 . For the fastest and most energetic young pulsar, PSR J0537−6910 (with a 16 ms spin period), we detect four more glitches, for a total of 15 glitches over 4.5 yr of NICER monitoring, and show that its spin-down behavior continues to set this pulsar apart from all others, including a long-term braking index n = −1.234 ± 0.009 and interglitch braking indices that asymptote to ≲7 for long times after a glitch. For PSR J1101−6101, we measure a much more accurate spin-down rate that agrees with a previous value measured without phase connection. For PSR J1412+7922 (also known as Calvera), we extend the baseline of its timing model from our previous 1 yr model to 4.4 yr, and for PSR J1849−0001, we extend the baseline from 1.5 to 4.7 yr. We also present a long-term timing model of the energetic pulsar PSR J1813−1749, by fitting previous radio and X-ray spin frequencies from 2009–2019 and new ones measured here using 2018 NuSTAR and 2021 Chandra data.
We report the experimental results of controlling the pedestal-top electron density by applying resonant magnetic perturbation (RMP) with in-vessel control coils and the main gas puff in the 2024-2025 KSTAR experimental campaign. The density is reconstructed using a parameterized $ψ_N$ grid and five channels of line-averaged density measured by the two-colored interferometer (TCI). The reconstruction procedure is accelerated by deploying a multi-layer perceptron to run in approximately 120 µ s and is sufficiently fast for real-time control. A proportional-integral controller was adopted, with the controller gains estimated from the system identification procedure. The experimental results demonstrate that the developed controller can follow a dynamic target while exclusively using both actuators. The absolute percentage errors between the electron density at $ψ_N$ = 0.89 and the target were approximately 1.5% median and a 2.5% average, respectively. The developed controller can even lower the density by using the pump-out mechanism under RMP, and it can follow a more dynamic range of density targets than a single actuator controller. The developed controller will enable experimental scenario exploration within a shot by dynamically setting the density target or maintaining a constant electron density within a discharge.
Here, we investigate the body-centered cubic (bcc) to hexagonal close-packed (hcp) phase transition in Fe-10wt. %Ni alloy, combining pressure cycling and fast compression with time-resolved synchrotron x-ray diffraction in a dynamic diamond anvil cell. Three pressure cycles were conducted with compression rates ranging from 0.1 to nearly 10 3 GPa/s. During the first cycle with the slowest compression, the observed orientations in the bcc and hcp phases are consistent with the Burgers mechanism, followed by c-axis rotation of the hcp phase consistent with {$10$ $\bar{1}2$} twinning. During the following cycles with fast compression at 10 2 –10 3 GPa/s, the hcp phase exhibits negligible c-axis rotation with a nearly constant c/a ratio of ∼1.61 up to ∼30 GPa, indicating suppression of plastic deformation (especially, twinning) due to sample confinement from the gasket. Notably, the onset pressure of the transition decreases with additional pressure cycling and faster compression, which normally leads to over-pressurization. This suggests that defects or shear induced from the pressure cycling reduces the transition kinetics even during fast compression. These insights into the deformation and transition behavior in an FeNi alloy under multiple dynamic loading cycles can offer guidance for future design of advanced structural alloys and improve our understanding of planetary core processes.
A CAEN DT5742 digitizer has been investigated to process the fast signals from the photomultiplier tubes of the time-of-flight detectors for fast ion beams. A test setup consisting of two plastic scintillation detectors and a pulsed laser source provided signals which were recorded by the digitizer and systematically analyzed with different algorithms to derive the amplitude, rise time and arrival time of the detection signal. To obtain the best amplitude and time resolutions, various optimization techniques including peak fitting and signal smoothing, interpolation improvement, and time-walk correction by amplitude and rise time have been performed and compared in detail. Here we compared the time resolutions obtained by three digital algorithms of leading-edge, zero-crossing constant-fraction (ZC-CFD) and direct constant-fraction discriminations. Finally, we found that the best time-of-flight resolution between two PMTs can be achieved as 12 ps by the method using sample minimum for peak location and using the time-walk-resistant ZC-CFD with 4-point interpolation for timing.
Oak Ridge National Laboratory (ORNL) collaborated with the Université de Sherbrooke (UdeS) to design a modular and scalable tile of innovative digital photodetectors called Photon-to-Digital Converters (PDC) and the required subsystems to demonstrate the feasibility of large area single-photon detectors to support high-resolution, depth-of-interaction, fast neutron radiography based on scintillation detectors. The goal is to demonstrate the advantages of using PDC over conventional Silicon Photomultiplier (SiPM) for neutron imaging systems. The basic building blocks of analog SiPM and PDC are Single-Photon Avalanche Photodiodes (SPAD). The main difference between conventional SiPMs and PDCs, is in the fact that analog SiPMs sum the charge produced by individual SPAD passively, while in PDCs each SPAD is read out individually by an active electronic circuit. Hence, PDCs provide a direct photon-to-bit conversion where a logic “1” means that there was a detection in a given SPAD. Conversely, an analog SiPM requires a sophisticated preamplifier (current amplifiers or transimpedance amplifier) followed by a shaping amplifier and an analog-to-digital converter. Moreover, since the charge from each SPAD varies slightly, the passive sum of the analog SiPM will have signal fluctuations for the same amount of photon detected. This is one of many issues that is completely eliminated by individual one-to-one SPAD read out in PDCs. The large area required for scintillator readouts for neutron imaging is such that the output capacitance of a large area of SiPM arrays is very high. This creates a burden on the signal-to-noise optimization with respect to the power budget. For PDC, the power consumption is dictated by rate of the incident photon flux, or in photon-starved environments, the dark noise rate of the SPAD array. In other words, the power consumption in the absence of events is extremely low. Hence, the power consumption of a PDC read out system is much lower than its analog counterpart. Further, the large SiPM capacitance may introduce signal distortions on the scintillator fast rise or decay time with undesired effects on timing and pulse shape discrimination. The problem is absent in PDCs, as their readout is independent on device capacitance. UdeS is the world leader in 3D vertical integration of frontside illuminated SPAD arrays over CMOS readouts. These capabilities allow future production of “3D PDC” with maximum photosensitive fill factor and electronic functionalities tailored to the application. UdeS' vast knowledge of SPADs and CMOS design made the team an ideal collaborator in the development of large area high-resolution fast neutron radiography.
During tokamak disruptions, the profile of the net parallel current is observed to flatten on a time scale that is so fast that it must be due to a fast magnetic reconnection. After a fast magnetic reconnection has broken magnetic surfaces, a single magnetic field line covers an entire volume and not just a magnetic surface. The current profile, given by K ≡ μ 0 j | | / B , relaxes to a constant within that volume by Alfven waves propagating along the chaotic magnetic field lines. The time scale for this relaxation determines the commonly observed disruption phenomena of a current spike and a sudden drop in the plasma internal inductance. An efficient method for studying this relaxation is derived, which allows a better understanding of the information encoded in the current spike and the associated sudden drop in the plasma internal inductance. Implications for coronal heating are also discussed.