Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Real numbers”

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 325 records · Page 18

Real Time Observation of Lithium Insertion into Pre-Cycled Conversion-Type Materials

Conversion-type electrode materials for lithium-ion batteries experience significant structural changes during the first discharge–charge cycle, where a single particle is taken apart into a number of nanoparticles. This structural evolution may affect the following lithium insertion reactions; however, how lithiation occurs in pre-cycled electrode materials is elusive. In this work, in situ transmission electron microscopy was employed to see the lithium-induced structural and chemical evolutions in pre-cycled nickel oxide as a model system. The introduction of lithium ions induced the evolution of metallic nickel, with volume expansion as a result of a conversion reaction. After pre-cycling, the phase evolutions occurred in two separate areas almost at the same time. This is different from the first lithiation, where the phase change takes place successively, with a boundary dividing the reacted and unreacted areas. Structural changes were restricted at the areas having large amount of fluorine, implying the residuals from the decomposition of electrolytes may have hindered the electrochemical reactions. This work provides insights into phase and chemical evolutions in pre-cycled conversion-type materials, which govern electrochemical properties during operation.

25 ENERGY STORAGE↗

Scaling of two-phase water-steam relative permeability and thermal fluxes in porous media

Two-phase water-steam flow conditions are frequently encountered in many engineering applications, including geothermal reservoirs. Although routine calculations are based on the multiphase Darcy’s law, the role of the topology of the flowing phases at the pore-scale is usually neglected in the estimation of relative permeabilities. Instead, the latter are frequently computed using empirical models like the Corey correlation. In this work, we first apply the model for relative permeabilities based on pore-scale flow regimes developed by Picchi and Battiato (2019), Relative permeability scaling from pore-scale flow regimes, Water Resour. Res. 55, 3215–3233, to scenarios typical of geothermal reservoirs and then extend it by deriving the scaling laws for the transmissibilities and the thermal properties as a function of temperature. First, we discuss the scaling behavior of normalized relative permeabilities in terms of viscosity ratio and capillary number of water-steam systems and, then, we provide a validation of the model against experimental data available in the literature. The model captures well the data trends collected in real 3D porous media. These results suggest that water-steam relative permeabilities follow the same scaling behavior of gas-liquid systems where the non-wetting phase is much less viscous than the wetting phase. Finally, we investigate the impact that relative permeabilities have on heat transfer rates at two-phase flow conditions and the scaling of mass and energy transmissibility and thermal properties of the mixture. An estimation of the exergy carried by the two-phase water-steam mixture is also included.

42 ENGINEERING↗

Ultrafast Spin Crossover in a Room-Temperature Ferrimagnet: Element-Specific Spin Dynamics in Photoexcited Cobalt Ferrite

Transition metal complexes capable of photoinduced spin crossover have been greatly investigated because of their potential to enable ultrafast optical control of information processing. Yet, any real application of photoswitchable molecules requires that spin crossover be paired with additional functionality such as long-range magnetic order. Important advances combining these functions are notably reported for a number of bimetallic Prussian Blue analogues; however, to date, PBA-based magnetic photoswitches can only operate below 150 K due to loss of magnetic order. In contrast, cobalt ferrite is a ferrimagnetic semiconductor with a Curie temperature of 790 K and extremely favorable magnetic properties by comparison to state-of-the-art PBAs. The mixed valence electronic structure of cobalt ferrite is reminiscent of cobalt–iron PBA, which is a well-known photoswitch. To investigate the potential for photoswitching in this material, we employ transient XUV spectroscopy to probe charge and spin dynamics with element-specific resolution on the femtosecond time scale. Results show that 400 nm light excites a metal-to-metal charge transfer transition, which drives the crossover of high-spin Co 2+ to low-spin Co 3+ with a time constant of 405 ± 29 fs and an internal quantum efficiency of unity. This result establishes the existence of efficient photoswitching in a new class of robust ferrimagnetic spinel ferrites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Edge Reconstruction in a Quantum Spin Hall Insulator

We study interaction-driven edge reconstruction in a quantum spin Hall insulator described by the Bernevig-Hughes-Zhang model with Kanamori-Hubbard interactions using the real-space density matrix renormalization group method in both the grand-canonical and canonical ensembles. For a two-dimensional cylinder with a smooth edge, we identify discrete particle-number transitions that lead to a spin-polarized edge state stabilized by an emergent ferromagnetic exchange interaction. The reconstruction is orbital-selective, occurring predominantly in the 𝑠-orbital channel. Our results reveal a microscopic mechanism for emergent fluctuating moments at the edge that could compromise the topological protection of helical edge states by time-reversal symmetry.

Soni, Rahul [ORNL] (ORCID:0000000317714299)↗

Quasiprobabilistic Readout Correction of Midcircuit Measurements for Adaptive Feedback via Measurement Randomized Compiling

Quantum measurements are a fundamental component of quantum computing. However, on present-day quantum computers, measurements can be more error prone than quantum gates and are susceptible to nonunital errors as well as nonlocal correlations due to measurement crosstalk. While readout errors can be mitigated in postprocessing, this is inefficient in the number of qubits due to a combinatorially large number of possible states that need to be characterized. In this work, we show that measurement errors can be tailored into a simple stochastic error model using randomized compiling, enabling the efficient mitigation of readout errors via quasiprobability distributions reconstructed from the measurement of a single preparation state in an exponentially large confusion matrix. We demonstrate the scalability and power of this approach by correcting readout errors without matrix inversion on a large number of different preparation states applied to a register of eight superconducting transmon qubits. Moreover, we show that this method can be extended to midcircuit measurements used for active feedback via quasiprobabilistic error cancellation, and we demonstrate the correction of measurement errors on an ancilla qubit used to detect and actively correct bit-flip errors on an entangled memory qubit. Our approach enables the correction of readout errors on large numbers of qubits and offers a strategy for correcting readout errors in adaptive circuits in which the results of midcircuit measurements are used to perform conditional operations on nonlocal qubits in real time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reoptimization of Quantum Circuits via Hierarchical Synthesis

The current phase of quantum computing is in the Noisy Intermediate-Scale Quantum (NISQ) era. On NISQ devices, two-qubit gates such as CNOTs are much noisier than single-qubit gates, so it is essential to minimize their count. Quantum circuit synthesis is a process of decomposing an arbitrary unitary into a sequence of quantum gates, and can be used as an optimization tool to produce shorter circuits to improve overall circuit fidelity. However, the time-to-solution of synthesis grows exponentially with the number of qubits. As a result, synthesis is intractable for circuits on a large qubit scale. In this paper, we propose a hierarchical, block-by-block opti-mization framework, QGo, for quantum circuit optimization. Our approach allows an exponential cost optimization to scale to large circuits. QGo uses a combination of partitioning and synthesis: 1) partition the circuit into a sequence of independent circuit blocks; 2) re-generate and optimize each block using quantum synthesis; and 3) re-compose the final circuit by stitching all the blocks together. We perform our analysis and show the fidelity improvements in three different regimes: small-size circuits on real devices, medium-size circuits on noisy simulations, and large-size circuits on analytical models. Our technique can be applied after existing optimizations to achieve higher circuit fidelity. Further, using a set of NISQ benchmarks, we show that QGo can reduce the number of CNOT gates by 29.9% on average and up to 50% when compared with industrial compiler optimizations such as t|ket). When executed on the IBM Athens system, shorter depth leads to higher circuit fidelity. We also demonstrate the scalability of our QGo technique to optimize circuits of 60+ qubits, Our technique is the first demonstration of successfully employing and scaling synthesis in the compilation tool chain for large circuits. Overall, our approach is robust for direct incorporation in production compiler toolchains to further improve the circuit fidelity.

97 MATHEMATICS AND COMPUTING↗

Physics-informed latent neural operator for real-time predictions of time-dependent parametric PDEs

Deep operator network (DeepONet) has shown significant promise as surrogate models for systems governed by partial differential equations (PDEs), enabling accurate mappings between infinite-dimensional function spaces. However, when applied to systems with high-dimensional input-output mappings arising from large numbers of spatial and temporal collocation points, these models often require heavily overparameterized networks, leading to long training times. Latent DeepONet addresses some of these challenges by introducing a two-step approach: first learning a reduced latent space using a separate model, followed by operator learning within this latent space. While efficient, this method is inherently data-driven and lacks mechanisms for incorporating physical laws, limiting its robustness and generalizability in data-scarce settings. Here, in this work, we propose PI-Latent-NO, a physics-informed latent neural operator framework that integrates governing physics directly into the learning process. Our architecture features two coupled DeepONets trained end-to-end: a Latent-DeepONet that learns a low-dimensional representation of the solution, and a Reconstruction-DeepONet that maps this latent representation back to the physical space. By embedding PDE constraints into the training via automatic differentiation, our method eliminates the need for labeled training data and ensures physics-consistent predictions. The proposed framework is both memory and compute-efficient, exhibiting near-constant scaling with problem size and demonstrating significant speedups over traditional physics-informed operator models. We validate our approach on a range of parametric PDEs, showcasing its accuracy, scalability, and suitability for real-time prediction in complex physical systems.

Latent representations↗

Stochastic Real-Time Second-Order Green’s Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green’s function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(N e 3 ), where N e is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

In situ thermal conductivity measurement revealing kinetics of thermochemical reactions

Utilizing thermochemical reactions for thermal energy storage and solar fuel production has been an emerging research topic. Thermal transport properties of the materials are an important parameter that can determine the kinetics and efficiency of thermochemical reactions. With the increasing number of new thermochemical materials (TCMs); however, there is a lack of reliable techniques to monitor the thermal transport property of the materials and their changes as a function of reactions in real time. In this work, we report the in situ monitoring of thermochemical reactions using modulated photothermal radiometry (MPR). The thermal conductivities of two TCMs, namely, calcium hydroxide (Ca(OH) 2 ) and Ba 0.15 Sr 0.85 FeO 3–δ (BSF1585), were measured as a function of temperature and time using the MPR technique. The measured thermal conductivities were correlated to the reaction. The work has two significant contributions to the research communities. First, it provides a non-invasive diagnostic tool for monitoring the thermal transport properties of TCMs that can potentially be a high-throughput measurement technique conducive to optimizing TCMs, reactors, and related thermal systems. Second, for TCMs that show observable changes in thermal transport properties, a correlation between the measured thermal conductivity and the conversion fraction of the reaction can be established for monitoring the reaction kinetics based on thermal characterization.

14 SOLAR ENERGY↗

Core-excited states of SF 6 probed with soft-x-ray femtosecond transient absorption of vibrational wave packets

A vibrational wavepacket in SF 6 is created by impulsive stimulated Raman scattering with a few-cycle infrared pulse and mapped simultaneously onto five sulfur core-excited states using table-top soft x-ray transient absorption spectroscopy between 170 to 200 eV. The femtosecond vibrations induce real-time energy shifts of the x-ray absorption, whose amplitude depend strongly on the nature of the core-excited state. The pump laser intensity is used to control the number of vibrational states in the superposition, thereby accessing core-excited levels for various extensions of the S-F stretching motion. This enables the determination of the relative core-level potential energy gradients for the symmetric stretching mode, in good agreement with TDDFT calculations. This experiment demonstrates a new means of characterizing core-excited potential energy curves.

74 ATOMIC AND MOLECULAR PHYSICS↗

Meta-optic accelerators for object classifiers

Rapid advances in deep learning have led to paradigm shifts in a number of fields, from medical image analysis to autonomous systems. These advances, however, have resulted in digital neural networks with large computational requirements, resulting in high energy consumption and limitations in real-time decision-making when computation resources are limited. Here, we demonstrate a meta-optic–based neural network accelerator that can off-load computationally expensive convolution operations into high-speed and low-power optics. In this architecture, metasurfaces enable both spatial multiplexing and additional information channels, such as polarization, in object classification. End-to-end design is used to co-optimize the optical and digital systems, resulting in a robust classifier that achieves 93.1% accurate classification of handwriting digits and 93.8% accuracy in classifying both the digit and its polarization state. This approach could enable compact, high-speed, and low-power image and information processing systems for a wide range of applications in machine vision and artificial intelligence.

42 ENGINEERING↗

Systems, methods and computer program products for self-tuning sensor data processing

Systems and methods are disclosed that include tools that utilize Dynamic Detector Tuning (DDT) software that identifies near-optimal parameter settings for each sensor using a neuro-dynamic programming (reinforcement learning) paradigm. DDT adapts parameter values to the current state of the environment by leveraging cooperation within a neighborhood of sensors. The key metric that guides the dynamic tuning is consistency of each sensor with its nearest neighbors: parameters are automatically adjusted on a per station basis to be more or less sensitive to produce consistent agreement of detections in its neighborhood. The DDT algorithm adapts in near real-time to changing conditions in an attempt to automatically self-tune a signal detector to identify (detect) only signals from events of interest. The disclosed systems and methods reduce the number of missed legitimate detections and the number of false detections, resulting in improved event detection.

Draelos, Timothy J.↗

Using Best Basis Inventory Data to Direct Strategies for Real Time Monitoring of Hanford High Level Waste – 26226

The potential to accelerate the processing of low- and high-level tank waste by applying real-time monitoring (RTM) of chemical and physical properties has prompted research into the suitability of multiple analytical methods for that purpose. The broad variety of waste stream properties and the large number of analytes of interest (as evidenced by Waste Acceptance Criteria (WAC) and Process Control Limit (PCL) lists) lead to an overwhelming set of possible analytical scenarios. This report describes the use of Best Basis Inventory (BBI) data to find the most relevant analytical targets for the specific case of monitoring the blending of High Level Waste from multiple tanks prior to introduction into a vitrification facility. Campaigns for blending this waste to minimize the risk of exceeding WACs and PCLs have been proposed. However, the predicted compositions of the blended materials do not incorporate any uncertainties that may be associated with the representativeness of the waste layer samples or the laboratory analyses that generated the BBI data. Also, they do not include any uncertainty associated with the precision of collecting highly specific fractions of the layers during a blending campaign or any inhomogeneities that may exist in those layers. Monte Carlo methods are used to apply uncertainties to the compositions of the individual layers specified in the campaign recipes. The resulting variations in the compositions of the blended materials allow estimation of the risks of exceeding WACs and PCLs for each campaign. A critical subset of WACs/PCLs – NOx, NaK, AlFeZr, and S – are especially at risk of being exceeded in multiple campaigns. These analytes should be the focus of instrument development. We also have extracted the expected solid/supernate distribution for these analytes, which establishes important performance criteria for individual analytical methods. The BBI data also permits an understanding of the different chemical forms in which the analytes appear. Thus, the need to establish instrumental sensitivity to these forms can be gainfully addressed. Although concentrating on one specific application – the blending of tank waste - this approach should be generalizable for the analysis of other possible RTM applications for waste processing.

Lascola, Robert [Savannah River National Laborator↗

Neural network-based classification and regression of magnetohydrodynamic modes in tokamaks

We present a machine learning-based magnetohydrodynamic (MHD) classifier and regressor that utilizes real or complex-valued 3D magnetic sensor array data to determine neoclassical tearing mode (NTM) onset times in tokamaks with millisecond accuracy. The input dataset consists of poloidal profiles of complex Fourier amplitudes with an n = 1 toroidal mode number from 144 human-labeled ITER Baseline Scenario discharges in the DIII-D tokamak, spanning both tearing-dominated and sawtooth-dominated regimes. Since m, n = 2,1 NTMs frequently emerge alongside sawteeth at the same frequency in this scenario, the focus is on isolating the m = 1 and m = 2 components of the n = 1 MHD mode near the tearing onset. To improve model regularization and prediction stability, singular value decomposition was applied to balance the sawtooth and tearing datasets. The enriched datasets facilitated training neural networks that learn the key distinguishing features of sawtooth and tearing modes in the poloidal profiles of their magnetic amplitude and phase. When the modes occur independently, the networks achieve perfect classification due to the modes’ distinct characteristics and low measurement noise. In the more experimentally relevant case where both modes coexist, the networks maintain exceptional performance across key metrics. Tests on synthetic data with known ground truth demonstrate the superior accuracy of the neural network trained on complex-valued input compared to models using real amplitude, phase, or pseudo-complex data, achieving both a mean time delay and standard deviation below 1 ms. Notably, standard linear regression methods fitting the dominant singular modes to the data closely match the neural network’s performance. Applying these methods across a broad range of H-mode scenarios will enable future studies to systematically identify dominant NTM triggers as scenario-specific variables, paving the way for more effective tearing mode avoidance strategies in future fusion reactor designs.

machine learning↗

Parametric Study of SANSMIC Input and Resulting Impact on Predicted Cavern Geometry and Leaching Efficiency

The Sandia Solution Mining Code (SANSMIC) has been used for many years to examine the development of salt cavern geometry, both in a confirmatory manner with comparisons made to real-world sonar data and in a predictive manner when updated sonar data are not available. SANSMIC models require some modeling choices in order to incorporate real-world data. Key modeling choices include the vertical resolution of cavern geometry to implement, as well as how to incorporate daily raw water injection data into the SANSMIC model. This report documents five studies that address the impact of the modeling choices on the predicted cavern geometries and calculated leaching efficiencies. In most cases, hypothetical cylindrical initial cavern geometries are used to provide a common baseline against which to test the systematic variation of input variables including cavern radius, oil-brine-interface (OBI) depth, vertical cell size, raw water injection rate, raw water injection duration, workover time, and number of leaching stages. The use of smaller cell sizes is recommended moving forward to provide a better one-to-one relationship between sonar data and the modeled cavern. A new methodology for incorporating raw water injection data is also recommended, in order to more closely model real-world injection and workover times. Overall, the systematic studies performed here have increased our confidence in previous SANSMIC model results, as well future use of the code for predicting leaching effects on cavern geometries. Some minor changes to modeling choices are recommended, which can easily be applied with the version of SANSMIC currently under development.

42 ENGINEERING↗

Real-time time-dependent density functional theory in Quantum Espresso and related codes

A set of subroutines that interact with an existing computer program, Quantum Espresso (QE), that takes a set of electronic orbitals from QE and propagates them forward in time in response to a number of perturbations. These perturbations include the motion of ions and a wide spectrum of applied electromagnetic fields (e.g., static magnetic fields, slowly-varying electric fields, x-rays, etc.). From the time evolution of the electronic orbitals comes time-varying values of physical observables like the electronic charge and current densities and atomic forces. These are further post-processed to yield properties of interest like stopping powers, conductivities, and the dynamic structure factor. The electronic orbitals are described in a plane wave basis and the equations of motion are carried out using a unitary or approximately unitary time propagation algorithm. Related work will include subroutines for interfacing this capability, as well as others that already exist within QE, with a fork of QE that is more extensible and modular.SAND2019-12611 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Baczewski, Andrew↗

A General Framework for Progressive Data Compression and Retrieval

In scientific simulations, observations, and experiments, the transfer of data to and from disk and across networks has become a major bottleneck for data analysis and visualization. Compression techniques have been employed to tackle this challenge, but traditional lossy methods often demand conservative error tolerances to meet the numerical accuracy requirements of both anticipated and unknown data analysis tasks. Progressive data compression and retrieval has emerged as a promising solution, where each analysis task dictates its own accuracy needs. However, few analysis algorithms inherently support progressive data processing, and adapting compression techniques, file formats, client/server frameworks, and APIs to support progressivity can be challenging. Here, this paper presents a framework that enables progressive-precision data queries for any data compressor or numerical representation. Our strategy hinges on a multi-component representation that successively reduces the error between the original and compressed field, allowing each field in the progressive sequence to be expressed as a partial sum of components. We have implemented this approach with four established scientific data compressors and assessed its effectiveness using real-world data sets from the SDRBench collection. The results show that our framework competes in accuracy with the standalone compressors it is based upon. Additionally, (de)compression time is proportional to the number of components requested by the user. Finally, our framework allows for fully lossless compression using lossy compressors when a sufficient number of components are employed.

97 MATHEMATICS AND COMPUTING↗