Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix product operators”

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 OSTI2020

A matrix-free hyperviscosity formulation for high-order ALE hydrodynamics

The numerical approximation of compressible hydrodynamics is at the core of high-energy density (HED) multiphysics simulations as shocks are the driving force in experiments like inertial confinement fusion (ICF). In this work, we describe our extension of the hyperviscosity technique, originally developed for shock treatment in finite difference simulations, for use in arbitrarily high-order finite element methods for Lagrangian hydrodynamics. Hyperviscosity enables shock capturing while preserving the high-order properties of the underlying discretization away from the shock region. Specifically, we compute a high-order term based on a product of the mesh length scale to a high power scaled by a hyper-Laplacian operator applied to a scalar field. We then form the total artificial viscosity by taking a non-linear blend of this term and a traditional artificial viscosity term. We also present a matrix-free formulation for computing the finite element based hyper-Laplacian operator. Such matrix-free methods have superior performance characteristics compared to traditional full matrix assembly approaches and offer advantages for GPU based HPC hardware. We demonstrate the numerical convergence of our method and its application to complex, multi-material ALE simulations on high-order (curved) meshes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Numerically exact configuration interaction at quadrillion-determinant scale

The combinatorial growth of configuration interaction (CI) has long limited this formally exact quantum chemistry method to only the smallest molecules. Here, we report a numerically exact CI calculation exceeding one quadrillion (10 15 ) determinants, made possible by a lossless categorical compression strategy within the small-tensor-product distributed active space (STP-DAS) framework. This approach overcomes the traditional memory bottlenecks of CI by a numerically exact compression of the wavefunction representation and reformulating the most computationally demanding matrix–vector operations. Using this method, we performed a fully relativistic CI calculation of the ground state of HBrTe with over 10 15 complex-valued determinants in just 34.5 h on 1000 computing nodes—the largest CI calculation ever reported. We further achieved fast computation for systems with hundreds of billions of determinants on only a few compute nodes. Extensive benchmarks confirm that the method retains full numerical exactness while cutting memory and computational cost by orders of magnitude. Compared to previous state-of-the-art CI calculations, this work achieves a 1000 times increase in CI space, a 10 6 -fold increase in floating-point operations performed, and a 10 6 -fold improvement in computational speed.

Computational chemistry↗
DOE OSTI2020

Online Data-Enabled Predictive Control

We develop an online data-enabled predictive (ODeePC) control method for trajectory tracking of unknown systems, building upon the recently proposed DeePC. Our proposed ODeePC method leverages a primal-dual algorithm with real-time measurement feedback to iteratively compute the corresponding real-time optimal control policy as system conditions change. Specifically, our developed ODeePC: a) records data from the unknown system and updates the underlying primal-dual algorithm dynamically, b) can track changes in the system's operating point and adjust the control inputs, and c) is computationally efficient as it deploys a Fast Fourier Transform-based algorithm enabling the fast computation of the product of a non-square Hankel matrix with a vector. We provide theoretical guarantees regarding the asymptotic behavior of ODeePC and demonstrate its performance through a power system application.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

Denoising of imaginary time response functions with Hankel projections

Imaginary-time response functions of finite-temperature quantum systems are often obtained with methods that exhibit stochastic or systematic errors. Reducing these errors comes at a large computational cost—in quantum Monte Carlo simulations, the reduction of noise by a factor of two incurs a simulation cost of a factor of four. In this paper, we relate certain imaginary-time response functions to an inner product on the space of linear operators on Fock space. We then show that data with noise typically does not respect the positive definiteness of its associated Gramian. The Gramian has the structure of a Hankel matrix. As a method for denoising noisy data, we introduce an alternating projection algorithm that finds the closest positive definite Hankel matrix consistent with noisy data. We test our methodology at the example of fermion Green's functions for continuous-time quantum Monte Carlo data and show remarkable improvements of the error, reducing noise by a factor of up to 20 in practical examples. We argue that Hankel projections should be used whenever finite-temperature imaginary-time data of response functions with errors is analyzed, be it in the context of quantum Monte Carlo, quantum computing, or in approximate semianalytic methodologies. Published by the American Physical Society 2024

Yu, Yang (ORCID:0000000186178878)↗
DOE OSTI2023

Conclusions from 3 years of continuous capture plant operation without exchange of the AMP/PZ-based solvent at Niederaussem – insights into solvent degradation management

A many times heard mantra of solvent degradation management in amine-based post combustion capture is “keep the solvent clean” to minimize solvent consumption. It is assumed that the amine losses would decrease by the removal of metals, degradation products, and reactive trace components which are captured from the flue gas, like NO 2 (as potentially driving components of the amine degradation besides dissolved O 2 ). However, this theoretical hypothesis – based on results from laboratory experiments typically generated with fresh amines – disregards the complexity of the solvent matrix, interaction of potential metal catalysts with degradation products and oxidizing agents, and specific chemical requirements which must be fulfilled before a degradation mechanism can proceed. Degradation of the solvent CESAR1 (aqueous solution of 3.0 molar 2-amino-2-methylpropan-1-ol (AMP) and 1.5 molar piperazine (PZ)) is investigated in a unique long-time test campaign (testing time up to now 40 months; 24/7 operation) without replacement of the solvent inventory at the capture pilot plant at the lignite-fired power plant in Niederaussem. Three solvent management strategies with different effect mechanisms are investigated and evaluated: (a) removal of only anionic compounds and trace elements (within 75 days solvent inventory treated two times) and anionic as well as cationic compounds and trace elements (114 days, inventory treated four times) from the solvent by ion exchange, (b) adsorptive removal of trace elements from the solvent by active carbon in 35% of the operating time, and (c) removal of >80% NO 2 by flue gas pretreatment with thiosulfate/sulfite solution (dosing for 2,000 h). The results of the testing program clearly show that “solvent cleanliness” is not a well-defined parameter and that results from laboratory tests, tests without fully representative industrial flue gasses, and short-term testing of monoethanolamine cannot be generalized for other solvents and industrial application. Furthermore, these results showcase that specific degradation management considering solvent, capture plant and flue gas quality is reasonable. Overshooting efforts for solvent management are contra-productive and produce unnecessary waste streams, efficiency losses and costs.

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

Matrix-free preconditioning for high-order H (curl) discretizations

The greater arithmetic intensity of high-order finite element discretizations makes them attractive for implementation on next-generation hardware, but assembly of high-order finite element operators as matrices is prohibitively expensive. As a result, the development of general algebraic solvers for such operators has been an open research challenge. Fast matrix-free application of high-order operators has received significant attention in the literature in the context of Poisson-type problems, but preconditioners and solvers for inverting more general operators are not very well-developed. In this paper, we consider the problem of preconditioning a definite Maxwell operator at high polynomial order without assembling a matrix. We show that given efficient preconditioners for high-order H 1 finite element problems on the same mesh, efficient H(curl) preconditioners can be constructed in an auxiliary space framework. We demonstrate the resulting preconditioners in a practical setting with tensor-product basis functions on an unstructured mesh of quadrilaterals. Overall, our approach uses a sparsified H 1 solver constructed on a low-order mesh of the nodal points of the underlying high-order space, and we show that the resulting H(curl) preconditioner is effective at very high polynomial orders for two-dimensional model problems with complicated geometry, varying piecewise constant coefficients, and curved elements. The resulting preconditioner scales with nearly optimal O(p d+1 ) floating point operation count and optimal O(p d ) memory transfer requirements, outperforming existing Maxwell preconditioners in the high-order regime.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

Quantum Search in Hilbert Space

A proposed quantum-computing algorithm would perform a search for an item of information in a database stored in a Hilbert-space memory structure. The algorithm is intended to make it possible to search relatively quickly through a large database under conditions in which available computing resources would otherwise be considered inadequate to perform such a task. The algorithm would apply, more specifically, to a relational database in which information would be stored in a set of N complex orthonormal vectors, each of N dimensions (where N can be exponentially large). Each vector would constitute one row of a unitary matrix, from which one would derive the Hamiltonian operator (and hence the evolutionary operator) of a quantum system. In other words, all the stored information would be mapped onto a unitary operator acting on a quantum state that would represent the item of information to be retrieved. Then one could exploit quantum parallelism: one could pose all search queries simultaneously by performing a quantum measurement on the system. In so doing, one would effectively solve the search problem in one computational step. One could exploit the direct- and inner-product decomposability of the unitary matrix to make the dimensionality of the memory space exponentially large by use of only linear resources. However, inasmuch as the necessary preprocessing (the mapping of the stored information into a Hilbert space) could be exponentially expensive, the proposed algorithm would likely be most beneficial in applications in which the resources available for preprocessing were much greater than those available for searching.

Zak, Michail↗
NASA NTRS1973

Multivariable system synthesis with step disturbance rejection

The primary objective of this paper is to present a constructive procedure for the synthesis of linear multivariable systems whose entire internal state and output are corrupted by unknown step disturbances. It is assumed that the dynamical behavior of the system is expressed in any one of three equivalent ways - i.e., the state space representation, the controllable and observable differential operator representation, and the transfer matrix representation. In terms of the transfer matrix representation, T(s), the synthesis procedure is shown to be capable of producing any stable, desired closed loop transfer matrix, Td(s), which can be expressed as the product of T(s) and any proper rational matrix, Tc(s), while simultaneously eliminating the steady-state effect of step disturbances at the output of the system. Furthermore, the synthesis scheme outlined employes only the known, directly measurable input and output signals.

Wolovich, W. A.↗
DOE OSTI2026

Unveiling the sea: universality of the transverse momentum dependent quark distributions at small x

Within the Colour Glass Condensate effective theory, we demonstrate that back-to-back dijet correlations in dilute-dense collisions involving a small-quark from the nuclear target can be factorised in terms of universal transverse momentum dependent distributions (TMDs) for the sea quarks. Two building blocks are needed to construct all these TMDs at the operator level: the sea quark TMD operator which appears in semi-inclusive Deep-Inelastic Scattering (SIDIS) or in the Drell-Yan process and the elastic-matrix for a quark-antiquark dipole. Compared to SIDIS, the saturation effects are stronger for dijet production in forward proton-nucleus collisions, due to additional scattering in the initial and final state, effectively resulting in a larger value for the nuclear saturation momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Thermalization from quantum entanglement: Jet simulations in the massive Schwinger model

We investigate the emergence of thermalization in a quantum-field-theoretic model mimicking the production of jets in QCD: the massive lattice Schwinger model coupled to external sources. Specifically, we compute the expectation values of local operators as functions of time and compare them to their thermal counterparts, quantify the overlap between the evolving density matrix and the thermal one, and compare the dynamics of the energy-momentum tensor to predictions from relativistic hydrodynamics. Through these studies, we find that the system approaches thermalization at late times and elucidate the mechanisms by which quantum entanglement drives thermalization in closed field-theoretic systems. Our results show how thermodynamic behavior emerges in real time from unitary quantum dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Mass Spectrometer Transient Analysis

This software implements a complete preprocessing pipeline for transient mass spectrometry (MS) data collected during TAP (Temporal Analysis of Products) experiments. It is designed to extract chemically meaningful fluxes from overlapping ion signals by applying a calibrated defragmentation matrix and solving the resulting linear system using non-negative least squares (NNLS) regression. The core script, preprocess_mass_spec.py, performs the following operations: Gain correction: Applies amplifier gain scalars derived from inert-packed calibration pulses to normalize signal intensities across AMUs and acquisition settings. Background subtraction: Removes experiment baselines using user-defined time windows, ensuring compatibility with slow-diffusing species and preventing negative values that would interfere with NNLS. Options to subtract before and after defragmentation. Defragmentation: Constructs a fragmentation matrix A from zeroth moments of calibration pulses (equal molar gas:inert mixtures) and solves Ax=b at each time point, where b is the raw MS signal and x is the estimated species flux. The matrix is normalized to inert signals and accounts for instrument-specific fragmentation behavior. Pulse-mode handling: Supports both averaged and individual pulse modes, enabling statistical treatment of fluxes and calculation of standard deviations. Integration and output: Computes zeroth moments (integrated fluxes) and exports time-resolved and integrated data in CSV format, suitable for downstream kinetic modeling. The software is validated using both virtual TAP simulations (VTAP) and experimental data from propane dehydrogenation (PDH) on CrOx/Al2O3 catalysts. It preserves temporal resolution by applying NNLS point-by-point across the pulse duration (typically 6,000+ time slices per pulse), leveraging the linear superposition principle to reconstruct full flux profiles. The defragmented outputs are compatible with kinetic extraction methods such as the G and Y procedures, which are used to derive rate–concentration relationships from TAP data. The details of these validations are discussed in detail in the supporting manuscript and supporting information. Example data and output files are also included. The methodology is robust to experimental noise and drift, with calibration protocols that account for pulse size effects, MS aging, and inert gas normalization. The software is modular, reproducible, and tailored for high-throughput TAP-MS workflows in catalysis research.

Kristy, Stephen [Idaho National Laboratory (INL), ↗
DOE OSTI2020

Sparse Matrix-Based HPC Tomography

Tomographic imaging has benefited from advances in X-ray sources, detectors and optics to enable novel observations in science, engineering and medicine. These advances have come with a dramatic increase of input data in the form of faster frame rates, larger fields of view or higher resolution, so high performance solutions are currently widely used for analysis. Tomographic instruments can vary significantly from one to another, including the hardware employed for reconstruction: from single CPU workstations to large scale hybrid CPU/GPU supercomputers. Furthermore, flexibility on the software interfaces and reconstruction engines are also highly valued to allow for easy development and prototyping. This paper presents a novel software framework for tomographic analysis that tackles all aforementioned requirements. The proposed solution capitalizes on the increased performance of sparse matrix-vector multiplication and exploits multi-CPU and GPU reconstruction over MPI. Furthermore, the solution is implemented in Python and relies on CuPy for fast GPU operators and CUDA kernel integration, and on SciPy for CPU sparse matrix computation. As opposed to previous tomography solutions that are tailor-made for specific use cases or hardware, the proposed software is designed to provide flexible, portable and high-performance operators that can be used for continuous integration at different production environments, but also for prototyping new experimental settings or for algorithmic development. The experimental results demonstrate how our implementation can even outperform state-of-the-art software packages used at advanced X-ray sources worldwide.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Developing an experimental-computational framework to investigate the deformation mechanisms and mechanical properties of Al-8Ce-10Mg alloys at micro and macroscales

There is a promising future for the use of aluminum-cerium-magnesium alloys in a broad range of applications, including devices that operate at high temperatures. With cerium currently considered a waste product of rare earth mining validation studies of possible applications are essential to reduce the environmental waste. Here, a computational-experimental framework is developed to investigate the role both the intermetallic and matrix have on the mechanical properties of these alloys. A set of experiments, including SEM/EBSD imaging, nanoindentation, in-situ SEM tensile testing, and in-situ SEM-DIC tests are performed to characterize the microstructure and mechanical properties of these alloys. Furthermore, the elastic, plastic, and failure deformation mechanisms of the microstructure, and their correlation with the bulk scale mechanical properties are investigated. Experimental results are also used to calibrate parameters for a crystal plasticity finite element model, by performing a computational framework that minimizes the error between the computational and experimental results. This model is then utilized to investigate how the area percentage of intermetallics and the crystallographic texture control the mechanical properties of the alloy. Simulation results show that an increase in the percentage of intermetallics increase the strength but decrease the ductility of the alloy. Also, a change in material texture improves strength and reduces damage that leads to material failure. The development of the crystal plasticity model, as discussed in this work, opens opportunities for future investigations of similar aluminum-cerium-magnesium alloys.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Eigenvalue properties of structural mean-axis systems

The authors review some of the properties of the pseudoinverse and oblique pseudoinverse of a linear transformation T from one finite-dimensional inner-product space into another, and then use these properties and a theorem of Milne (1968), which states that the oblique pseudoinverse can be expressed in terms of a weak generalized inverse and two projection operators, in order to compute a mean-axis influence coefficient matrix for the dynamic analysis of an elastic body. Some eigenvalue invariance properties of the mean-axis structural dynamics equations are demonstrated, and on the simple example of a uniform beam, it is shown that the finite frequencies and mode shapes of the mean axis structural system are identical to the nonzero frequencies and mode shapes of the free structure.

Cavin, R. K., III↗
DOE OSTI2026

Adaptive Stiffness Power Take-Off System for a Hydraulic Hose-Pump Point Absorber Wave Energy Converter

Two crucial challenges faced in the design of ocean wave energy converters (WECs) are maintaining high capacity factors and remaining robust in everchanging, harsh ocean environments. This study introduces a novel mechanism for achieving reactive WEC control with a soft hydraulic power take-off, thereby addressing both challenges. The device of interest involves a hydraulic pump, constructed from a fluidic flexible matrix composite (F2MC), which is cyclically stretched by driving waves and pressurizes working fluid through a turbine-generator. This class of fiber-reinforced stretch hoses is established for ocean use as robust moorings and shock absorbers. The dynamic properties of the F2MC pump depend on internal pressure, which can be controlled by a continuously variable transmission between the turbine and generator, allowing the system to be adapted for better performance over a wide range of wave conditions. This work presents a dynamic model for operation and adaptive behavior of this power take-off system, experimentally validates components of the model, and performs dynamic simulations of the system that show increases in power production and capacity factor compared to non-adaptive operation for a variety of driving wave conditions. This work develops the basis for future large-scale experimental testing and controller design of this system.

Wave energy conversion, fluidic flexible matrix co↗
NASA NTRS2001

Microgravity Emissions Laboratory Developed

The Microgravity Emissions Laboratory (MEL) was developed for the support, simulation, and verification of the International Space Station microgravity environment. The MEL utilizes an inertial measurement system using acceleration emissions generated by various operating components of the space station. These emissions, if too large, could hinder the science performed on the space station by disturbing the microgravity environment. Typical test components are disk drives, pumps, motors, solenoids, fans, and cameras. These components will produce inertial forces, which disturb the microgravity on-orbit station environment. These components, usually housed within a station rack, must meet acceleration limits imposed at the rack interface for minimizing the onboard station-operating environment. The NASA Glenn Research Center developed this one-of-a-kind laboratory for testing components and, eventually, rack-level configurations. The MEL approach is to measure the component's generated inertial forces. This force is a product of the full diagonal mass matrix including the test setup (the center of gravity, mass moment of inertia, and weight) and the resolved diagonal rigid-body acceleration determined from measurements using the 10 apparatus accelerometers. The mass matrix can be test derived. The bifilar torsional pendulum method is used to measure the moment of inertia for the test component.

Goodnight, Thomas W.↗
NASA NTRS2024

Northern Rockies Ecological Conservation: Leveraging Earth Observations to Monitor and Predict Populations of Federally Threatened Whitebark Pine (Pinus albicaulis) across the Intermountain West

Whitebark pine (WBP; Pinus albicaulis) is an ecologically important species in North America. As a federally listed threatened species, an understanding of WBP habitat, distribution, and health is important for the natural resource managers of the National Park Service, United States Forest Service, Bureau of Land Management, Fish and Wildlife Service, and non-profit organizations such as the Whitebark Pine Ecosystem Foundation. Previous attempts to develop models of WBP habitat suitability and distribution lack confidence in their validity and integrity for these organizations. The updated models of habitat suitability and distribution developed by this study would provide managers with a capability to be employed in the conservation and future research direction for WBP. Thus, we developed a habitat suitability model of WBP at a high spatial resolution (Landsat 9 Operational Land Image-2, National Land Cover Database, NASA Shuttle Radar Topography Mission; 30m pixels) using a generalized logistic regression with an area under the curve value of 0.754. We extracted spectral reflectance signatures from overlapped ground sample points and Sentinel-2 Multispectral Instrument. The spectral signature analysis indicates WBP is separable from other tree species. We also utilized a visual validation approach and random forest (RF) modeling to separate WBP from limber pine. Through visual validation the RF classifier successfully identified 8out of 10 WBP trees gathered through ground truth points. Additionally, we achieved an overall accuracy of 91%in our confusion matrix for the distribution model using a dependent validation approach. The derived products from this study allow project partners to assess current suitable habitat and apparent health status in areas of identified WBP occurrence, providing data to aid future research regarding WBP health.

Sentinel-2↗