Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “polynomial method”

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

Inference of response functions with the help of machine-learning algorithms

Response functions are a key quantity to describe the near-equilibrium dynamics of strongly interacting many-body systems. Recent techniques that attempt to overcome the challenges of calculating these ab initio have employed expansions in terms of orthogonal polynomials. We employ a neural network prediction algorithm to reconstruct a response function 𝑆⁡(𝜔) defined over a range in frequencies 𝜔. Here, we represent the calculated response function as a truncated Chebyshev series whose coefficients can be optimized to reduce the representation error. We compare the quality of response functions obtained using coefficients calculated using a neural network (NN) algorithm with those computed using the Gaussian integral transform (GIT) method. In the regime where only a small number of terms in the Chebyshev series are retained, we find that the NN scheme outperforms the GIT method.

Kurkcuoglu, Doga Murat [Fermi National Accelerator↗
DOE OSTI2024

The constructive method for massive particles in QED

The constructive method of determining amplitudes from on-shell pole structure has been shown to be promising for calculating amplitudes in a more efficient way. However, challenges have been encountered when a massless internal photon is involved in the gluing of three-point amplitudes with massive external particles. In this paper, we use the original on-shell method, old-fashioned perturbation theory, to shed light on the constructive method, and show that one can derive the Feynman amplitude by correctly identifying the residue even when there is an internal photon involved.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Polynomial Chaos Surrogate Construction for Random Fields with Parametric Uncertainty

Engineering and applied science rely on computational experiments to rigorously study physical systems. The mathematical models used to probe these systems are highly complex, and sampling-intensive studies often require prohibitively many simulations for acceptable accuracy. Surrogate models provide a means of circumventing the high computational expense of sampling such complex models. In particular, polynomial chaos expansions (PCEs) have been successfully used for uncertainty quantification studies of deterministic models where the dominant source of uncertainty is parametric. We discuss an extension to conventional PCE surrogate modeling to enable surrogate construction for stochastic computational models that have intrinsic noise in addition to parametric uncertainty. We develop a PCE surrogate on a joint space of intrinsic and parametric uncertainty, enabled by Rosenblatt transformations, which are evaluated via kernel density estimation of the associated conditional cumulative distributions. Furthermore, we extend the construction to random field data via the Karhunen–Loève expansion. We then take advantage of closed-form solutions for computing PCE Sobol indices to perform a global sensitivity analysis of the model which quantifies the intrinsic noise contribution to the overall model output variance. Additionally, the resulting joint PCE is generative in the sense that it allows generating random realizations at any input parameter setting that are statistically approximately equivalent to realizations from the underlying stochastic model. The method is demonstrated on a chemical catalysis example model and a synthetic example controlled by a parameter that enables a switch from unimodal to bimodal response distributions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Measurement and calibration of non-linear shear terms in galaxy cluster fields

Introduction Galaxy cluster lensing is a powerful tool for measuring the mass of galaxy clusters, but accurate shear measurement and calibration are critical to obtaining reliable results. This study focuses on the measurement and calibration of weak lensing shears to improve mass estimates in cluster lensing. To deal with the problem, we first developed an image simulation pipeline, jedisim, which utilizes galaxy images extracted from the Hubble Space Telescope (HST) Ultra Deep Field (UDF) and the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Methods The simulations represent realistic galaxy distributions and morphologies as input sources. The foreground halo with a Navarro–Frenk–White (NFW) profile is constructed such that the lensing signals of background galaxies can be measured by the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) Science Pipelines. By comparing the measured reduced shear g meas and the true reduced shear g true , we observe non-linearity up to g ≲ 0.6 . We fit polynomials to the data with quadratic correction adequate to g ≲ 0.4 . Meanwhile, we conduct mass estimates using the pzmassfitter code on four different clusters. Results The mass estimate results are significantly improved after applying the shear calibration derived from the present work—from 4.954 ± 0.504 × 1 0 14 M ⊙ to 10.507 ± 0.498 × 1 0 14 M ⊙ after calibration for a simulated cluster with the mass of 10 × 1 0 14 M ⊙ . In multiple cases of validation, the estimated results are all consistent with true cluster mass. Discussion This study yields the first relationship between reality and shape measurement of the LSST Science Pipelines and serves as the first step toward the overall goal of mass calibration in cluster lensing. By addressing the challenges in shear measurement and calibration, we aim to enhance the accuracy and reliability of mass estimates in galaxy cluster lensing studies.

Astronomy & Astrophysics↗
DOE OSTI2026

Uncertainty quantification for Joule heating processes in fibrous pore-resolved media

Joule heating (JH) is an energy-efficient and sustainable technique for heating materials. Its application for industrial heating, particularly, has been gaining attention due to its potential for increasing the yield of various chemical products. The process involves the use of heating elements (materials that are highly conductive electrically and thermally) to heat up other materials or substances. These conductors, however, can exhbit varying degrees of uncertainty due to non-linearities in their temperature-dependent properties, which could result in variable material behavior. In this work, we carry out uncertainty quantification (UQ) at the pore scale to describe the uncertainty of such materials. In so doing, we applied the non-intrusive polynomial chaos expansion (PCE) technique to quantify the uncertainty within the system. The steady state Joule heating equation was solved numerically at the pore scale mimicking conditions within a heating chamber for propane dehydrogenation, and various electro-thermal profiles were obtained. We also examined the effect of the number of sampling points (20 – 100) and order of the PCE coefficients (2 – 5) on the accuracy of the temperature evaluations. The results were then benchmarked with the standard Monte Carlo (MC) method. The average temperature of the 4th-order global PCE showed good agreement with the MC results (which were positively skewed). Orders greater than 4 gave an underestimation of the temperatures while predictions for the peak temperature improved as the number of sampling points increased.

Fagbemi, Samuel [ORNL] (ORCID:0000000236995025)↗
DOE OSTI2026

Polynomial-time preparation of low-temperature Gibbs states for two-dimensional toric code

In this work, we propose a polynomial-time algorithm for preparing the Gibbs state of the two-dimensional toric code Hamiltonian at any temperature, starting from any initial state, significantly improving upon prior estimates that suggested exponential scaling with inverse temperature. We prove that fast mixing at low temperature for the two-dimensional toric code can be achieved by augmenting local jump operators with simple global jump operators, which enable efficient transitions between logical sectors. To establish tight lower bounds on the spectral gap, we introduce a new reduction method that eventually maps the problem to estimating the spectral gap of a perturbed graph Laplacian on a stair graph. Our proof also shows that the Lindblad dynamics with a digitally implemented low-temperature local Davies generator is able to efficiently drive the quantum state toward the ground state manifold.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Galerkin formulation of path integrals in lattice field theory

We present a mathematical framework for Galerkin formulations of path integrals in lattice field theory. The framework is based on using the degrees of freedom (DOFs) associated to a Galerkin discretization as the fundamental lattice variables. We formulate standard concepts in lattice field theory, such as the partition function and correlation functions, in terms of the DOFs. For example, using continuous finite element spaces, we show that the two-point spatial correlation function can be defined between any two points on the domain (as opposed to at just lattice sites) and furthermore, this two-point function satisfies a weak propagator (or Green’s function) identity, in analogy to the continuum case, as well as a convergence estimate obtained from the standard finite element techniques. Furthermore, this framework leads naturally to higher-order formulations of lattice field theories by considering higher-order finite element spaces for the Galerkin discretization. We consider analytical and numerical examples of scalar field theory to investigate how increasing the order of piecewise polynomial finite element spaces affect the approximation of lattice observables. Finally, we sketch an outline of this Galerkin framework in the context of gauge field theories.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Variance-Reduced Accelerated First-Order Methods: Central Limit Theorems and Confidence Statements

In this paper, we consider a strongly convex stochastic optimization problem and propose three classes of variable sample-size stochastic first-order methods: (i) the standard stochastic gradient descent method, (ii) its accelerated variant, and (iii) the stochastic heavy-ball method. In each scheme, the exact gradients are approximated by averaging across an increasing batch size of sampled gradients. We prove that when the sample size increases at a geometric rate, the generated estimates converge in mean to the optimal solution at an analogous geometric rate for schemes (i)–(iii). Based on this result, we provide central limit statements, whereby it is shown that the rescaled estimation errors converge in distribution to a normal distribution with the associated covariance matrix dependent on the Hessian matrix, the covariance of the gradient noise, and the step length. If the sample size increases at a polynomial rate, we show that the estimation errors decay at a corresponding polynomial rate and establish the associated central limit theorems (CLTs). Under certain conditions, we discuss how both the algorithms and the associated limit theorems may be extended to constrained and nonsmooth regimes. As a result, we provide an avenue to construct confidence regions for the optimal solution based on the established CLTs and test the theoretical findings on a stochastic parameter estimation problem.

Lei, Jinlong↗
DOE OSTI2026

Enabling accurate chemical modeling of shocked energetic materials using a machine learning interatomic potential

Understanding the complex chemistry of organic materials under dynamic compression is important for many applications, but it is challenging due to the large number of reactions occurring at various time scales. Here, in this study, we develop a machine learning potential based on Chebyshev polynomials to study the insensitive energetic material 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) under detonation. We discuss a strategy for constructing diverse training data needed to capture the complex chemistry of TATB. Our potential demonstrates strong transferability across a wide range of thermodynamic conditions and other explosives, enabling accurate and reliable chemical modeling of organic materials under extreme conditions. The efficiency of our approach allows for simulations over several nanoseconds and for large system sizes, providing detailed insights into the chemistry of shocked TATB. The model accurately reproduces experimental Hugoniot equation of state data, and our simulations reveal the rapid formation of nitrogen-rich carbon clusters following shock. The methods and datasets developed here offer a robust framework for accurate chemical modeling of other shocked organic energetic materials.

Chemistry↗
DOE OSTI2024

High-Order Wall-Modeled Large-Eddy Simulation of High-Lift Configuration

This paper presents the assessment of several recent enhancements for a high-order wall-modeled large-eddy simulation (WMLES) approach and demonstrates order independence with a fixed data exchange location in the wall model. The two enhancements include the use of isotropic tetrahedral elements to improve accuracy and an explicit subgrid-scale model, the Vreman model, to improve accuracy and robustness. The [Formula: see text] study focused on the high-lift Common Research Model (HL-CRM) at the angle of attack of 19.57 deg, a benchmark problem from the 4th AIAA High-Lift Prediction Workshop. Solution polynomial orders of [Formula: see text], and 5 were used in the study. The study demonstrated [Formula: see text] independence in integrated forces, pitch moment, velocity profile in the wall-normal direction, and surface flow topology. It also showed that a [Formula: see text] order of at least 3 ([Formula: see text]) was needed to correctly predict the external inviscid flow and the surface flow topology. Thereafter, [Formula: see text] simulations over several other angles of attack demonstrated that the high-order WMLES approach can correctly predict the maximum lift and flow separation regions for HL-CRM with about 40 million degrees of freedom (DOF) compared to at least 250 million DOF required by second-order methods.

Engineering↗
DOE OSTI2025

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗
DOE OSTI2025

Lie-algebraic classical simulations for quantum computing

The classical simulation of quantum dynamics plays an important role in our understanding of quantum complexity and in the development of quantum technologies. Efficient techniques such as those based on the Gottesman-Knill theorem for Clifford circuits, tensor networks for low entanglement-generating circuits, or Wick's theorem for fermionic Gaussian states have become central tools in quantum computing. In this work, we contribute to this body of knowledge by presenting a framework for classical simulations, dubbed “𝔤-sim”, which is based on the underlying Lie algebraic structure of the dynamical process. When the dimension of the algebra grows at most polynomially in the system size, there exist observables for which the simulation is efficient. Indeed, we show that 𝔤-sim enables new regimes for classical simulations, is able to deal with certain forms of noise in the evolution, as well as can be used to tackle several paradigmatic variational and nonvariational quantum computing tasks. For the former, we perform Lie-algebraic simulations to train and optimize parametrized quantum circuits (thus effectively showing that some variational models can be dequantized), design enhanced parameter initialization strategies, solve tasks of quantum circuit synthesis, and train a quantum-phase classifier. For the latter, we report large-scale noiseless and noisy simulations on benchmark problems. By comparing the limitations of 𝔤-sim and certain Wick's theorem-based simulations, we find that the two methods become inefficient for different types of states or observables, hinting at the existence of distinct, nonequivalent resources for classical simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Training quantum neural networks using the quantum information bottleneck method

Abstract We provide in this paper a concrete method for training a quantum neural network to maximize the relevant information about a property that is transmitted through the network. This is significant because it gives an operationally well founded quantity to optimize when training autoencoders for problems where the inputs and outputs are fully quantum. We provide a rigorous algorithm for computing the value of the quantum information bottleneck quantity within error ε that requires O ( log 2 ⁡ ( 1 / ϵ ) + 1 / δ 2 ) queries to a purification of the input density operator if its spectrum is supported on { 0 } ⋃ [ δ , 1 − δ ] for δ > 0 and the kernels of the relevant density matrices are disjoint. We further provide algorithms for estimating the derivatives of the QIB function, showing that quantum neural networks can be trained efficiently using the QIB quantity given that the number of gradient steps required is polynomial.

Çatlı, Ahmet Burak (ORCID:0000000152294141)↗
DOE OSTI2025

Uncertainty propagation in feed-forward neural network models

We develop new uncertainty propagation methods for feed-forward neural network architectures with leaky ReLU activation functions subject to random perturbations in the input vectors. In particular, we derive analytical expressions for the probability density function (PDF) of the neural network output and its statistical moments as a function of the input uncertainty and the parameters of the network, i.e., weights and biases. A key finding is that an appropriate linearization of the leaky ReLU activation function yields accurate statistical results even for large perturbations in the input vectors. This can be attributed to the way information propagates through the network. We also propose new analytically tractable Gaussian copula surrogate models to approximate the full joint PDF of the neural network output. To validate our theoretical results, we conduct Monte Carlo simulations and a thorough error analysis on a multi-layer neural network representing a nonlinear integro-differential operator between two polynomial function spaces. Our findings demonstrate excellent agreement between the theoretical predictions and Monte Carlo simulations.

MLP networks↗
DOE OSTI2025

Exponentially Reduced Circuit Depths Using Trotter Error Mitigation

Product formulas are a popular class of digital quantum simulation algorithms due to their conceptual simplicity, low overhead, and performance, which often exceeds theoretical expectations. Recently, Richardson extrapolation and polynomial interpolation have been proposed to mitigate the Trotter error incurred by the use of these formulas. This work provides a rigorous, general analysis of these techniques for computing time-evolved observables, simplifying the interpolation algorithm in the process, and shows that extrapolation generically improves the performance of product formulas for this task. We demonstrate that, to achieve error 𝜖 in a simulation of time 𝑇 using a 𝑝 ⁢th-order product formula with extrapolation, circuit depths of 𝑂⁡(𝑇 1+1/𝑝 ⁢polylog (1/𝜖)) are sufficient—an exponential improvement in the precision over product formulas alone. Furthermore, we prove that these algorithms achieve commutator scaling, and improve the 𝑇 complexity for the interpolation algorithm. By relaxing the requirement of performing exact Chebyshev interpolation, our simplified algorithm eliminates the need for fractional implementations of Trotter steps, reducing computational overhead. Finally, we show these techniques can be combined with the classical shadows method to estimate many time-evolved local observables. Taken together, our findings provide the strongest evidence yet for the utility of Trotter error-mitigation techniques in algorithmic applications.

quantum algorithms & computation↗
DOE OSTI2024

The 1 → 3 massive splitting functions from QCD factorization and SCET

Splitting functions are universal functions describing the collinear dynamics of gauge theories, and as such are crucial ingredients for a wide variety of calculations in perturbative QCD. We present analytic results for the triple collinear splitting functions in QCD with a single massive parton. We derive the splitting functions using two distinct methods; first by expanding the squared matrix elements in the collinear limit, and secondly by using soft-collinear effective theory with massive quarks. We find agreement between these two approaches, providing a strong check of our results. Additionally, we also check all iterated and soft limits of our results, finding agreement with predictions from factorization. Our results provide an important ingredient for higher order perturbative calculations involving massive partons, and for the description of the collinear dynamics of heavy flavor jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗