Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗
DOE OSTI2020

Self-consistent GW method: O(N) algorithm for polarizability and self energy

In this work, an efficient implementation of the self-consistent GW method in the FlapwMBPT code [1] is presented. It features the evaluation of polarizability and self-energy which scales linearly with respect to the system size. Altogether the computational time scaling was measured to be be-tween linear and quadratic in the applications to silicon supercells with up to 72 atoms. Application to such materials as paracostibite CoSbS, supercells of La 2 CuO 4 (up to 56 atoms) and SmB 6 , illustrate the potential of the approach in computational material science.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A Contextually-Aware Sensitivity Analysis to Guide the Design of Randomized Least Squares Solvers in Applications

Our work on the DOE-sponsored project “A Contextually-Aware Sensitivity Analysis to Guide the Design of Randomized Least Squares Solvers in Applications,” was an effort to address critical challenges in nu merical computing and its applications to optimization. The increasing demand for robust and scalable solutions to large-scale linear algebra problems has highlighted the limitations of traditional approaches, particularly in heterogeneous and extreme-scale computing environments. Randomized Numerical Linear Algebra (RandNLA) offers a promising framework to address these challenges, and this proposal builds on this foundation by introducing innovations in sensitivity analysis and computational adaptability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Quantifying Seasonal and Diurnal Cycles of Solar‐Induced Fluorescence With a Novel Hyperspectral Imager

Abstract Solar‐induced fluorescence (SIF) is a proxy of ecosystem photosynthesis that often scales linearly with gross primary productivity (GPP) at the canopy scale. However, the mechanistic relationship between GPP and SIF is still uncertain, especially at smaller temporal and spatial scales. We deployed a ultra‐hyperspectral imager over two grassland sites in California throughout a soil moisture dry down. The imager has high spatial resolution that limits mixed pixels, enabling differentiation between plants and leaves within one scene. We find that imager SIF correlates well with diurnal changes in leaf‐level physiology and gross primary productivity under well‐watered conditions. These relationships deteriorate throughout the dry down event. Our results demonstrate an advancement in SIF imaging with new possibilities in remotely sensing plant canopies from the leaf to the ecosystem. These data can be used to resolve outstanding questions regarding SIF's meaning and usefulness in terrestrial ecosystem monitoring.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Solving nuclear structure problems with the adaptive variational quantum algorithm

Here we use the Lipkin-Meshkov-Glick (LMG) model and the valence-space nuclear shell model to examine the likely performance of variational quantum eigensolvers in nuclear-structure theory. The LMG model exhibits both a phase transition and spontaneous symmetry breaking at the mean-field level in one of the phases, features that characterize collective dynamics in medium-mass and heavy nuclei. We show that with appropriate modifications, the ADAPT-VQE algorithm, a particularly flexible and accurate variational approach, is not troubled by these complications. We treat up to 12 particles and show that the number of quantum operations needed to approach the ground-state energy scales linearly with the number of qubits. We find similar scaling when the algorithm is applied to the nuclear shell model with realistic interactions in the sd and pf shells. Although most of these simulations contain no noise, we use a noise model from real IBM hardware to show that for the LMG model with four particles, weak noise has no effect on the efficiency of the algorithm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A Preliminary Feasibility Study on Multi-Cavity Cryomodule Integration for the Electron Ion Collider Energy Recover Linac Cooler

The Electron-Ion Collider (EIC) is a state-of-the-art accelerator designed for collisions between highly polarized electrons and ions. To achieve optimal luminosity, the ion beam is cooled using an electron beam sourced from an energy recovery linac (ERL). In the current ERL baseline design, one BNL type RF cavity is used per cryomodule, leading to spatial and cost challenges. This study examines the feasibility of using more compact PERLE (Powerful Energy Recovery Linac for Experiments) type cavities, which can house multiple cavities within a single cryomodule, by evaluating their Beam Breakup (BBU) instability performance. Higher Order Modes (HOM) parameters were obtained through frequency scaling, assuming constant quality factor Q_L, R/Q, and linearly scaling HOM frequencies. To predict the BBU threshold current and ensure accuracy, two different BBU tracking simulations are used for cross-verification. Although the reduced footprint of PERLE-type cavities is advantageous, maintaining sufficient damping of HOMs remains crucial. Finally, we compare the HOM damping efficiency of both cavity types and suggest a pathway forward.

Deitrick, K.↗
DOE OSTI2022

Nature and Scalings of Density Fluctuations of Compressible Magnetohydrodynamic Turbulence with Applications to the Solar Wind

Abstract The solar wind is a magnetized and turbulent plasma. Its turbulence is often dominated by Alfvénic fluctuations and often deemed as nearly incompressible far away from the Sun, as shown by in situ measurements near 1 au. However, for solar wind closer to the Sun, the plasma β decreases (often lower than unity) while the turbulent Mach number M t increases (can approach unity, e.g., transonic fluctuations). These conditions could produce significantly more compressible effects, characterized by enhanced density fluctuations, as seen by several space missions. In this paper, a series of 3D MHD simulations of turbulence are carried out to understand the properties of compressible turbulence, particularly the generation of density fluctuations. We find that, over a broad range of parameter space in plasma β , cross helicity, and polytropic index, the turbulent density fluctuations scale linearly as a function of M t , with the scaling coefficients showing weak dependence on parameters. Furthermore, through detailed spatiotemporal analysis, we show that the density fluctuations are dominated by low-frequency nonlinear structures, rather than compressible MHD eigenwaves. These results could be important for understanding how compressible turbulence contributes to solar wind heating near the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

The (RH+t) aging correlation. Electrical resistivity of PVB at various temperatures and relative humidities

Electrical products having organic materials functioning as pottants, encapsulants, and insulation coatings are commonly exposed to elevated conditions of temperature and humidity. In order to assess service life potential from this method of accelerated aging, it was empirically observed that service life seems proportional to an aging correlation which is the sum of temperature in degrees Celsius (t), and the relative humidity (RH) expressed in percent. Specifically, the correlation involves a plot of time-to-failure on a log scale versus the variable RH + T plotted on a linear scale. A theoretical foundation is provided for this empirically observed correlation by pointing out that the correlation actually involves a relationship between the electrical resistivity (or conductivity) of the organic material, and the variable RH + t. If time-to-failure is a result of total number of coulombs conducted through the organic material, then the correlation of resistivity versus RH + t is synonymous with the empirical correlation of time-to-failure versus RH + t.

Cuddihy, E. F.↗
DOE OSTI2021

Prediction of Transition-State Scaling Relationships and Universal Transition-State Vibrational and Entropic Correlations for Dehydrogenations

Linear scaling relationships (LSRs) and Brønsted–Evans–Polanyi (BEP) or transition-state scaling (TSS) relations aid with the prediction of electronic energies. However, temperature effects and pre-exponentials are often taken as constants across metal surfaces or a homologous series. Vibrational scaling relationships (VSRs) offer a way to determine such parameters. Transition-state VSRs (TSVSRs) between local minima and transition states of AH X (A = C, N, O) surface diffusions correlate with BEP relations and broaden to thermochemical property scaling. Using density functional theory, we extend TSVSRs to AHX dehydrogenation reactions on transition-metal surfaces, relating vibrational modes of local minima to transition states. We first predict the slopes of the TSS relations by incorporating bond angles using the Slater–Koster structure factors and hybridization through crystal orbital overlap population analysis and energy overlap integrals between adsorbates and metal surfaces. Additionally, we uncover universal thermochemical property scaling, enabling the estimation of entropies and temperature corrections to enthalpies across a homologous series. Here, we demonstrate both significant vibrational corrections in reactions with low intrinsic electronic barriers and considerable variation in the pre-exponential of a simple dehydrogenation reaction across metals and AH X adsorbates.

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

Frontal Slice Approaches for Tensor Linear Systems

Inspired by the row and column action methods for solving large-scale linear systems, in this work, we explore the use of frontal slices for solving tensor linear systems. In particular, this paper presents a novel approach for using frontal slices of a tensor $\mathcal{A}$ to solve tensor linear systems $\mathcal{A} ∗\mathcal{X} = \mathcal{B}$ where ∗ denotes the $t$-product. In addition, we consider variations of this method, including cyclic, block, and randomized approaches, each designed to optimize performance in different operational contexts. Our primary contribution lies in the development and convergence analysis of these methods. Experimental results on synthetically generated and real-world data, including applications such as image and video deblurring, demonstrate the efficacy of our proposed approaches and validate our theoretical findings.

Luo, Hengrui↗
DOE OSTI2022

Evidence for Nonlinear Gluon Effects in QCD and Their Mass Number Dependence at STAR

The STAR Collaboration reports measurements of back-to-back azimuthal correlations of di-π 0 s produced at forward pseudorapidities (2.6<η<4.0) in p+p, p+Al, and p+Au collisions at a center-of-mass energy of 200 GeV. We observe a clear suppression of the correlated yields of back-to-back π 0 pairs in p+Al and p+Au collisions compared to the p+p data. The observed suppression of back-to-back pairs as a function of transverse momentum suggests nonlinear gluon dynamics arising at high parton densities. Furthermore, the larger suppression found in p+Au relative to p+Al collisions exhibits a dependence of the saturation scale $Q^{2}_{s}$ on the mass number A . A linear scaling of the suppression with A 1/3 is observed with a slope of -0.09±0.01.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Sparsity of the electron repulsion integral tensor using different localized virtual orbital representations in local second-order Møller–Plesset theory

Utilizing localized orbitals, local correlation theory can reduce the unphysically high system-size scaling of post-Hartree–Fock (post-HF) methods to linear scaling in insulating molecules. The sparsity of the four-index electron repulsion integral (ERI) tensor is central to achieving this reduction. For second-order Møller–Plesset theory (MP2), one of the simplest post-HF methods, only the (ia|jb) ERIs are needed, coupling occupied orbitals i, j and virtuals a, b. In this paper, we compare the numerical sparsity (called the “ragged list”) and two other approaches revealing the low-rank sparsity of the ERI. The ragged list requires only one set of (localized) virtual orbitals, and we find that the orthogonal valence virtual-hard virtual set of virtuals originally proposed by Subotnik et al. gives the sparsest ERI tensor. To further compress the ERI tensor, the pair natural orbital (PNO) type representation uses different sets of virtual orbitals for different occupied orbital pairs, while the occupied-specific virtual (OSV) approach uses different virtuals for each occupied orbital. Here, our results indicate that while the low-rank PNO representation achieves significant rank reduction, it also requires more memory than the ragged list. The OSV approach requires similar memory to that of the ragged list, but it involves greater algorithmic complexity. An approximation (called the “fixed sparsity pattern”) for solving the local MP2 equations using the numerically sparse ERI tensor is proposed and tested to be sufficiently accurate and to have highly controllable error. A low-scaling local MP2 algorithm based on the ragged list and the fixed sparsity pattern is therefore promising.

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

Computational Discovery of Active and Selective Metal‐Nitrogen‐Graphene Catalysts for Electrooxidation of Water to H 2 O 2

A direct electrosynthesis of H 2 O 2 from either O 2 or H 2 O is an attractive strategy to replace the energy-intensive industrial anthraquinone process. Two-electron water oxidation reaction (2e-WOR) offers several advantages over the oxygen reduction reaction such as better mass transfer due to the absence of gas-phase reactants. However, 2e-WOR is a more challenging and less studied process with only a handful of metal oxides exhibiting reasonable activity/selectivity properties. Herein, we employ density-functional-theory calculations to screen a variety of metal-nitrogen-graphene structures for 2e-WOR. As a consequence of scaling between the adsorption energies of reaction intermediates, we determine a linear relation between selectivities for the first and second reaction steps of 2e-WOR, viz. that if selectivity toward adsorbed OH is improved, then selectivity toward H 2 O 2 at the subsequent step is decreased. We also find that selectivity and activity are linearly scaled in such a way that a higher activity (i.e., a lower overpotential) leads to a lower selectivity for the H 2 O 2 formation step. Based on the obtained results several chemistries, e.g., containing NiN x -C moieties, are predicted to rival the best-performing metal oxides such as ZnO and CaSnO 3 in terms of combination of their activity/selectivity characteristics for 2e-WOR.

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

Block encoding of the three-dimensional heterogeneous Poisson equation with application to fracture flow

Quantum linear system (QLS) algorithms offer the potential to solve large-scale linear systems exponentially faster than classical methods. However, applying QLS algorithms to real-world problems remains challenging due to issues such as state preparation, data loading, and efficient information extraction. In this work, we study the feasibility of applying QLS algorithms to solve discretized three-dimensional (3D) heterogeneous Poisson equations, with specific examples relating to groundwater flow through geologic fracture networks. We explicitly construct a block encoding for the 3D heterogeneous Poisson matrix by leveraging the sparse local structure of the discretized operator. While classical solvers benefit from preconditioning, we show that block encoding the system matrix and preconditioner separately does not improve the effective condition number that dominates the QLS run-time. This differs from classical approaches where the preconditioner and the system matrix can often be implemented independently. Nevertheless, due to the structure of the problem in three dimensions, the quantum algorithm achieves a run-time of 𝑂⁡(𝑁 2/3 polylog 𝑁 ⋅log (1/𝜖)), outperforming the best classical methods (with run times of 𝑂⁡(𝑁⁢log 𝑁 ⋅log (1/𝜖))) and offering exponential memory savings. These results highlight both the promise and limitations of QLS algorithms for practical scientific computing, and point to effective condition-number reduction as a key barrier in achieving quantum advantages.

58 GEOSCIENCES↗
DOE OSTI2024

Gradient Coding With Iterative Block Leverage Score Sampling

Gradient coding is a method for mitigating straggling servers in a centralized computing network that uses erasure-coding techniques to distributively carry out first-order optimization methods. Randomized numerical linear algebra uses randomization to develop improved algorithms for large-scale linear algebra computations. In this study, we propose a method for distributed optimization that combines gradient coding and randomized numerical linear algebra. The proposed method uses a randomized ℓ 2 -subspace embedding and a gradient coding technique to distribute blocks of data to the computational nodes of a centralized network, and at each iteration the central server only requires a small number of computations to obtain the steepest descent update. The novelty of our approach is that the data is replicated according to importance scores, called block leverage scores, in contrast to most gradient coding approaches that uniformly replicate the data blocks. Furthermore, we do not require a decoding step at each iteration, avoiding a bottleneck in previous gradient coding schemes. We show that our approach results in a valid ℓ 2 -subspace embedding, and that our resulting approximation converges to the optimal solution.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

The aging correlation (RH + t): Relative humidity (%) + temperature (deg C)

An aging correlation between corrosion lifetime, and relative humidity RH (%) and temperature t (C) has been reported in the literature. This aging correlation is a semi-log plot of corrosion lifetime on the log scale versus the interesting summation term RH(%) + t(C) on the linear scale. This empirical correlation was derived from observation of experimental data trends and has been referred to as an experimental law. Using electrical resistivity data of polyvinyl butyral (PVB) measured as a function of relative humidity and temperature, it was found that the electrical resistivity could be expressed as a function of the term RH(%) t(C). Thus, if corrosion is related to leakage current through an organic insulator, which, in turn, is a function of RH and t, then some partial theoretical validity for the correlation is indicated. This article describes the derivation of the term RH(%) t(C) from PVB electrical resistivity data.

Cuddihy, E. F.↗
DOE OSTI2026

Quantum Ising model on (2+1)-dimensional anti–de Sitter space using tensor networks

We study the quantum Ising model on (2+1)-dimensional anti-de Sitter space using matrix product states (MPS) and matrix product operators (MPOs). We explore the bulk phase diagram of the theory on regular tessellations of hyperbolic space with coordination number seven and find disordered and ordered phases separated by a phase transition. We find that the boundary-boundary spin correlation function exhibits power law scaling deep in the disordered phase of the Ising model consistent with holography. At the critical point, we find the boundary entanglement entropy scales logarithmically with subsystem size but away from this, we see a linear scaling. In comparison, the full system exhibits a volume law scaling, which is expected in chaotic and/or highly connected systems. We also measure out of time ordered correlators (OTOCs) to explore the scrambling behavior of the theory.

Quantum spin models↗