Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Performance evaluations of signed and unsigned noisy approximate quantum Fourier arithmetic

The Quantum Fourier Transform (QFT) grants competitive advantages, especially in resource usage and circuit approximation, for performing arithmetic operations on quantum computers, and offers a potential route toward a numerical quantum-computational paradigm. In this paper, we utilize efficient techniques to implement QFT-based integer addition and multiplications. These operations are fundamental to various quantum applications including Shor’s algorithm, weighted-sum optimization problems in data processing and machine learning, and quantum algorithms requiring inner products. We carry out performance evaluations of these implementations based on IBM’s superconducting-qubit architecture using different compatible noise models. We isolate the sensitivity of the component quantum circuits on both one-/two-qubit gate error rates, and the number of the arithmetic operands’ superposed integer states. We analyze performance and identify the most effective approximation depths for unsigned quantum addition and quantum multiplication within the given context. We then perform a similar analysis of signed addition and compare to the unsigned results. We observe significant dependency of the optimal approximation depth on the degree of machine noise and the number of superposed states in certain performance regimes. Finally, we elaborate on the algorithmic challenges—relevant to signed, unsigned, modular and non-modular versions—that could also be applied to current implementations of QFT-based subtraction, division, exponentiation, and their potential tensor extensions. Here, we analyze the performance trends in our results and speculate on possible future developments within this computational paradigm.

Computational models↗
DOE OSTI2023

An assessment of the utility of multirate time integration for the modeling of laser powder bed fusion

Finite element simulation of the laser powder bed fusion process is made challenging by the inherently multiscale nature of the process. When using the typical global time stepping techniques, slowly-evolving regions of the domain receive the same numerical treatment as the regions with the highest temperature rates. The current work details the implementation and evaluation of an implicit multirate method which is able to advance different regions of the domain with distinct time steps, depending on their current solution rate. Previous work indicates that this representation of the temporal scales of the problem can lead to significant reductions in wall clock run time, and it is shown herein that multirate time integration, when used with a uniform process-scale mesh, can result in speedups of approximately 19, 42, and 87 times, for domains with edge lengths of approximately 1 mm, 2 mm, and 5 mm, respectively. When used in conjunction with $\ h$-refinement (limited to two levels of refinement), the resulting speedups (taken relative to the uniform mesh) are around 28, 85, and 665 times. The method is demonstrated to converge as indicated by the literature, and its use with an AM-Bench domain is demonstrated.

42 ENGINEERING↗
DOE OSTI2020

Time-series machine-learning error models for approximate solutions to parameterized dynamical systems

This work proposes a machine-learning framework for modeling the error incurred by approximate solutions to parameterized dynamical systems. In particular, we extend the machine-learning error models (MLEM) framework proposed in Ref. Freno and Carlberg (2019) to dynamical systems. The proposed Time-Series Machine-Learning Error Modeling (T-MLEM) method constructs a regression model that maps features – which comprise error indicators that are derived from standard a posteriori error-quantification techniques – to a random variable for the approximate-solution error at each time instance. The proposed framework considers a wide range of candidate features, regression methods, and additive noise models. We consider primarily recursive regression techniques developed for time-series modeling, including both classical time-series models (e.g., autoregressive models) and recurrent neural networks (RNNs), but also analyze standard non-recursive regression techniques (e.g., feed-forward neural networks) for comparative purposes. Finally, numerical experiments conducted on multiple benchmark problems illustrate that the long short-term memory (LSTM) neural network, which is a type of RNN, outperforms other methods and yields substantial improvements in error predictions over traditional approaches.

42 ENGINEERING↗
DOE OSTI2024

a priori uncertainty quantification of reacting turbulence closure models using Bayesian neural networks

While many physics-based closure model forms have been posited for the sub-filter scale (SFS) in large eddy simulation (LES), vast amounts of data available from direct numerical simulations (DNS) create opportunities to leverage data-driven modeling techniques. Albeit flexible, data-driven models still depend on the dataset and the functional form of the model chosen. Increased adoption of such models requires reliable uncertainty estimates both in the data-informed and out-of-distribution regimes. Here, in this work, we employ Bayesian neural networks (BNNs) to capture both epistemic and aleatoric uncertainties in a reacting flow model. In particular, we model the filtered progress variable scalar dissipation rate which plays a key role in the dynamics of turbulent premixed flames. We demonstrate that BNN models can provide unique insights about the structure of uncertainty of the data-driven closure models. We also propose a method for the incorporation of out-of-distribution information in a BNN, which can be used for out-of-distribution query detection. The efficacy of the model is demonstrated by a priori evaluation on a dataset consisting of a variety of flame conditions and fuels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Rheinheimera sp . T2C2 Bacterial Biofilm for Bioremediation of Cobalt(II)

Toxic metals, including cobalt, are often the cause of the contamination of rivers and lakes in mining regions. Heavy metal water pollution has been linked to numerous human health problems, prompting the need for environmental remediation. Existing techniques for removing heavy metals from water, such as chemical precipitation and filtration, produce toxic waste, are costly, or require high power consumption for pumping. Biosorption is a potential alternative strategy that is cost-effective and uses readily available and naturally produced biomass and living material to absorb pollutants. Engineering living materials, such as biofilms, which consist of living cells and a secreted polymer matrix, offer the potential to integrate toxin sensing, sequestration, and metabolism capabilities of cells to improve pollution remediation strategies. Alternative biofilm producing candidates need to be explored to implement these material capabilities. Previous biosorption studies have primarily used bacterial biofilms from known pathogens and/or generated toxic waste in the form of the absorbent material combined with the heavy metal. Here, we describe a recently isolated bacterium called Rheinheimera sp. T2C2 that forms biofilms with promising biosorption characteristics. T2C2 is an aquatic bacterium with low nutrient requirements and high biofilm production that is not known to be pathogenic. We demonstrate (1) the efficacy of Rheinheimera sp. T2C2 as a biosorbent for cobalt bioremediation; (2) how biosorption is altered by water conditions to establish the efficacy of this strategy in different environments; and (3) how the metal can be released from the biofilm for metal recycling. Our findings will provide a living materials strategy that overcomes the existing barriers for bioremediation and improves the health of ecosystems and humans through heavy metal removal and recycling.

Rheinheimera↗
DOE OSTI2022

Yield Estimation of the August 2020 Beirut Explosion by Using Physics‐Based Propagation Simulations of Regional Infrasound

Abstract The August 2020 Beirut explosion is the largest single‐fired ammonium nitrate explosion documented in history. The massive explosion excited loud infrasound in the atmosphere, and clear waveforms were recorded by a regional infrasound array at an epicentral distance of 100 km, allowing for accurate measurements of explosion energy. We estimate the explosion size based on the infrasound waveform inversion. Unlike conventional inversions using empirical models, we perform full 3‐D finite‐difference simulations to obtain a physics‐based propagation model for the inversion. Accurate numerical modeling of infrasound is challenging as the propagation is substantially affected by the turbulent atmosphere. Instead of a single deterministic prediction, we provide a range of waveform predictions by running multiple simulations with stochastic weather forecast models, which allows for comprehensive uncertainty analysis of numerical modeling and estimated yields. Finally, we expand the yield estimation technique for seismoacoustic analysis and demonstrate the substantial advantage of the joint approach.

58 GEOSCIENCES↗
DOE OSTI2023

Quantum error mitigation by Pauli check sandwiching

Abstract We describe and analyze an error mitigation technique that uses multiple pairs of parity checks to detect the presence of errors. Each pair of checks uses one ancilla qubit to detect a component of the error operator and represents one layer of the technique. We build on the results on extended flag gadgets and put it on a firm theoretical foundation. We prove that this technique can recover the noiseless state under the assumption of noise not affecting the checks. The method does not incur any encoding overhead and instead chooses the checks based on the input circuit. We provide an algorithm for obtaining such checks for an arbitrary target circuit. Since the method applies to any circuit and input state, it can be easily combined with other error mitigation techniques. We evaluate the performance of the proposed methods using extensive numerical simulations on 1850 random input circuits composed of Clifford gates and non-Clifford single-qubit rotations, a class of circuits encompassing most commonly considered variational algorithm circuits. We observe average improvements in fidelity of 34 percentage points with six layers of checks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Performance Evaluations of Noisy Approximate Quantum Fourier Arithmetic

The Quantum Fourier Transform (QFT) grants competitive advantages, especially in resource usage and circuit approximation, for performing arithmetic operations on quantum computers, and offers a potential route towards a numerical quantum-computational paradigm. In this paper, we utilize efficient techniques to implement QFT-based integer addition and multiplications. These operations are fundamental to various quantum applications including Shor’s algorithm, weighted sum optimization problems in data processing and machine learning and quantum algorithms requiring inner products. We carry out performance evaluations of these implementations based on IBM’s superconducting qubit architecture using different compatible noise models. We isolate the sensitivity of the component quantum circuits on both one-/two-qubit gate error rates, and the number of the arithmetic operands’ superposed integer states. We analyze performance, and identify the most effective approximation depths for quantum add and quantum multiply within the given context. We observe significant dependency of the optimal approximation depth on the degree of machine noise and the number of superposed states in certain performance regimes. Finally, we elaborate on the algorithmic challenges - relevant to signed, unsigned, modular and non-modular versions - that could also be applied to current implementations of QFT-based subtraction, division, exponentiation, and their potential tensor extensions. Here, we analyze performance trends in our results and speculate on possible future development within this computational paradigm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Two-Dimensional Non-Linear Magnetic Equivalent Circuit Model to Facilitate the Preliminary Design of a Normal Conducting Quadrupole

Normal conducting quadrupoles have been used to focus charged particle beams in synchrotrons, beam transfer lines, medical linacs, etc. for a long time. Optimization techniques based on analytical expressions combined with the use of a numerical field analysis tool exist in the literature for the design of an optimal pole tip shape. However, the initial shape and dimensions of the remaining yoke (including the pole itself and yoke base that act as the return path) are usually less well defined. This article discusses a design methodology for a normal conducting quadrupole based on a two-dimensional magnetic equivalent circuit. Here, this approach considers the geometry of the entire magnet and the non-linear behavior of the yoke material, thereby eliminating the initial iterations. The design outcome of this exercise serves as a competent starting point and can then be used to refine the pole tip, pole taper, yoke size, and add other geometrical features to achieve the required field quality, gradient, coil considerations, etc. by employing a finite element analysis tool. An example design study to demonstrate the proposed methodology is presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Investigation of Automated Corrective Actions for Special Protection Schemes

The constantly evolving nature of the grid is compelling the design process of Remedial Action Schemes (RAS) to keep up with the changes. This document proposes a flexible and computationally efficient approach to automatically determine RAS corrective actions that alleviate line overloading violations. Statistical and functional characteristics summarized from RAS implemented in real power systems are used to guide the design parameters. This report also leverages sensitivity-based techniques to determine corrective actions for specific contingencies quickly without repeated numerical simulations. Finally, future directions for implementing this approach for a fully automated, online RAS are discussed.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Development and Evaluation of a General Drag Model for Gas-Solid Flows via Deep Learning

This project presents the development and evaluation of a general drag model for gas–solid multiphase flows using deep learning techniques. A comprehensive database of more than 4,000 experimental and numerical data points for spherical and non spherical particles was compiled, incorporating geometric features such as sphericity, aspect ratio, and orientation. Several predictive approaches—including traditional em pirical correlations, machine learning, and deep neural networks—were benchmarked, with the proposed Drag Coefficient Correlation-aided Deep Neural Network (DCC DNN) demonstrating superior accuracy. To account for particle–particle interactions, additional drag data were generated using CFD-based simulations of packed and flu idized beds, leading to the development of a retrained model capable of incorporat ing volume fraction effects. Integration of the trained model with the MFiX CFD solver was achieved using FTorch, enabling drag predictions during discrete element method (DEM) simulations. Validation against experimental data for single particles and fluidized beds confirmed the model’s improved predictive ability, particularly for non-spherical geometries. While the model performed strongly under fluidized con ditions, limitations remained in unfluidized regimes, suggesting a need for expanded datasets. Overall, this study demonstrates the feasibility of combining deep learning with physics-informed CFD to improve drag modeling for gas–solid flows, with promis ing implications for scaling multiphase simulations in industrial applications.

42 ENGINEERING↗
DOE OSTI2020

Optimization and Prediction of Spectral Response of Metasurfaces Using Artificial Intelligence

Hot-electron generation has been a topic of intense research for decades for numerous applications ranging from photodetection and photochemistry to biosensing. Recently, the technique of hot-electron generation using non-radiative decay of surface plasmons excited by metallic nanoantennas, or meta-atoms, in a metasurface has attracted attention. These metasurfaces can be designed with thicknesses on the order of the hot-electron diffusion length. The plasmonic resonances of these ultrathin metasurfaces can be tailored by changing the shape and size of the meta-atoms. One of the fundamental mechanisms leading to generation of hot-electrons in such systems is optical absorption, therefore, optimization of absorption is a key step in enhancing the performance of any metasurface based hot-electron device. Here we utilized an artificial intelligence-based approach, the genetic algorithm, to optimize absorption spectra of plasmonic metasurfaces. Using genetic algorithm optimization strategies, we designed a polarization insensitive plasmonic metasurface with 90% absorption at 1550 nm that does not require an optically thick ground plane. We fabricated and optically characterized the metasurface and our experimental results agree with simulations. Finally, we present a convolutional neural network that can predict the absorption spectra of metasurfaces never seen by the network, thereby eliminating the need for computationally expensive simulations. Our results suggest a new direction for optimizing hot-electron based photodetectors and sensors.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Automatic Differentiation of C++ Codes on Emerging Manycore Architectures with Sacado

Automatic differentiation (AD) is a well-known technique for evaluating analytic derivatives of calculations implemented on a computer, with numerous software tools available for incorporating AD technology into complex applications. However, a growing challenge for AD is the efficient differentiation of parallel computations implemented on emerging manycore computing architectures such as multicore CPUs, GPUs, and accelerators as these devices become more pervasive. In this work, we explore forward mode, operator overloading-based differentiation of C++ codes on these architectures using the widely available Sacado AD software package. In particular, we leverage Kokkos, a C++ tool providing APIs for implementing parallel computations that is portable to a wide variety of emerging architectures. Here we describe the challenges that arise when differentiating code for these architectures using Kokkos, and two approaches for overcoming them that ensure optimal memory access patterns as well as expose additional dimensions of fine-grained parallelism in the derivative calculation. We describe the results of several computational experiments that demonstrate the performance of the approach on a few contemporary CPU and GPU architectures. We then conclude with applications of these techniques to the simulation of discretized systems of partial differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Analysis of the Feasibility of UAS-Based EMI Sensing for Underground Utilities Detection and Mapping

This paper investigates the feasibility of using a linear current sensing (LCS) technique integrated on an unmanned aerial system (UAS) for detecting and mapping underground infrastructure rapidly and cost-effectively. The LCS technique is based on data from a wide band of electromagnetic induction frequencies (50 kHz to 2 MHz) using a vector magnetic field gradiometer. This technique takes advantage of a slowly decaying secondary magnetic field in order to achieve greater standoff detection distance ($\frac{1}{R^2}$ vs. $\frac{1}{R^6}$ for compact metallic targets during EMI sensing, where R is the distance from a target to the sensor). These secondary magnetic fields are produced by the excite current on long conductors, allowing detection at a distance of 10 meters or more. The system operates between tens of kHz to a few MHz and uses either an active EMI source or existing EM fields to excite this linear current on a long metallic subsurface target. Once excited, these linear currents produce a secondary magnetic field that is detected with an above ground triaxial magnetic field gradiometer. By moving and tracking its geolocation, the system outputs rich datasets sufficient to support high-fidelity forward and inverse EMI models for estimating the depth and orientation of deep underground long linear metallic infrastructure. The system’s hardware and its integration to a UAS system are outlined, along with the formulation of LCS theory, and numerical and experimental data are presented. The results illustrate that the LCS technique offers large standoff detection, is adaptable to UAS, and could be used effectively for detecting deep underground infrastructure such as wires and pipes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Revealing hidden defects through stored energy measurements of radiation damage

With full knowledge of a material’s atomistic structure, it is possible to predict any macroscopic property of interest. In practice, this is hindered by limitations of the chosen characterization techniques. For example, electron microscopy is unable to detect the smallest and most numerous defects in irradiated materials. Instead of spatial characterization, we propose to detect and quantify defects through their excess energy. Differential scanning calorimetry of irradiated Ti measures defect densities five times greater than those determined using transmission electron microscopy. Our experiments also reveal two energetically distinct processes where the established annealing model predicts one. Molecular dynamics simulations discover the defects responsible and inform a new mechanism for the recovery of irradiation-induced defects. The combination of annealing experiments and simulations can reveal defects hidden to other characterization techniques and has the potential to uncover new mechanisms behind the evolution of defects in materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Cardinal: Seismic and Geoacoustic Array Processing

Data collected via seismic and infrasound array deployments are leveraged in the geosciences to detect and characterize a myriad of natural and anthropogenic sources. These deployments consist of numerous sensors placed in a predetermined configuration to amplify signal strength and improve the efficacy of array processing techniques used to measure signal directionality and waveform coherence. High‐fidelity feature extraction is often predicated on interstation distance as well as the frequency content and wavelength of an incident signal. Numerous array processing softwares analyze data in sequential frequency bands to obtain a more detailed characterization of a signal. However, current algorithms are limited in their ability to determine optimal array configuration for each band. We introduce an open‐source Python code, called Cardinal, to process seismic and infrasound array data in discretized time–frequency space with the option of applying an adaptive array design to determine optimal subarray configuration for each frequency band. To reduce computational time, the array processing step can be run in parallel using multithreading. Furthermore, the software has the capability to aggregate array processing results from different time–frequency pixels to produce separate sets of detections, or families, with added utility via the application of an adaptive semblance threshold, which aids in isolating signals‐of‐interest from coherent background noise. Upon appropriate configuration, Cardinal exhibits the potential to combine distinct seismic and infrasound phases into separate families.

Adaptive Array↗
DOE OSTI2026

An analytic result for the 0 → ggHHH amplitude

We present a fully analytic calculation of the leading-order one-loop amplitude for triple Higgs production via gluon fusion, gg → HHH, retaining full dependence on the mass of the heavy quark circulating in the loop. This amplitude provides a direct probe of the triple and quartic Higgs self-couplings, the measurement of which is a central goal of current and future colliders. The amplitude can be presented in compact form thanks to the use of analytic reconstruction techniques, based on finite-field and p-adic evaluations, multivariate partial fraction decompositions, and primary decompositions to identify common numerator factors. Although full analytic results are given in the text and in the supplementary material, the main thrust of this paper is to further test and illustrate these analytic reconstruction techniques in a concrete physical example. Our results provide a compact and efficient representation of the matrix element for this process, enabling evaluations that are more than an order of magnitude faster than existing numerical alternatives. Full analytic control of the leading-order, loop-induced amplitude is an important step towards handling more complex 2-loop or real-radiation corrections to this and related processes.

Campbell, John M. [Fermilab] (ORCID:00000002136193↗