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

Chiral spin liquid and quantum phase transition in the triangular-lattice Hofstadter-Hubbard model

Recent advances in moiré engineering motivate the study of lattice models of strongly correlated electrons subjected to substantial orbital magnetic flux. We analyze the triangular-lattice Hofstadter-Hubbard model at one-quarter flux quantum per plaquette and a density of one electron per site, where a chiral spin liquid phase may exist between weak-coupling integer quantum Hall and strong-coupling 120° antiferromagnetic phases. Here, we use matrix product state methods and analytical arguments to investigate this model compactified to cylinders of finite circumference. We uncover a glide particle-hole symmetry operation which, we argue, is spontaneously broken at the quantum Hall to spin liquid transition on odd-circumference cylinders. We numerically verify the spontaneous symmetry breaking and further demonstrate that this transition is associated with algebraic long-range correlations of various spin-singlet, charge-neutral operators. For even-circumference cylinders, the transition becomes a crossover associated with a large correlation length that grows substantially with circumference. Our findings suggest that in the two-dimensional limit, the transition to a chiral spin liquid phase is continuous and features critical fluctuations of the current.

Divic, Stefan [University of Pennsylvania, Philade↗
NASA NTRS1979

Satellite pattern classification using charge transfer devices

The potential uses of Charge Transfer Devices (CTDs) in pattern classification operations are explored. The needs for a hardware-based pattern classifier are established, and a matrix multiplication subsystem based upon a sum-of-products CTD is presented. Applications of the subsystem to the classification of multi-modal Gaussian distributions in general and to LANDSAT data processing in particular are discussed. Finally, the potential impact of this technology on satellite data processing methodologies is discussed.

Snyder, W. E.↗
NASA NTRS2003

Combustor and Vane Features and Components Tested in a Gas Turbine Environment

The use of ceramic matrix composites (CMCs) as combustor liners and turbine vanes provides the potential of improving next-generation turbine engine performance, through lower emissions and higher cycle efficiency, relative to today s use of superalloy hot-section components. For example, the introduction of film-cooling air in metal combustor liners has led to higher levels of nitrogen oxide (NOx) emissions from the combustion process. An environmental barrier coated (EBC) siliconcarbide- fiber-reinforced silicon carbide matrix (SiC/SiC) composite is a new material system that can operate at higher temperatures, significantly reducing the film-cooling requirements and enabling lower NOx production. Evaluating components and subcomponents fabricated from these advanced CMCs under gas turbine conditions is paramount to demonstrating that the material system can perform as required in the complex thermal stress and environmentally aggressive engine environment. To date, only limited testing has been conducted on CMC combustor and turbine concepts and subelements of this type throughout the industry. As part of the Ultra-Efficient Engine Technology (UEET) Program, the High Pressure Burner Rig (HPBR) at the NASA Glenn Research Center was selected to demonstrate coupon, subcomponent feature, and component testing because it can economically provide the temperatures, pressures, velocities, and combustion gas compositions that closely simulate the engine environments. The results have proven the HPBR to be a highly versatile test rig amenable to multiple test specimen configurations essential to coupon and component testing.

Roinson, R. Craig↗
DOE OSTI2024

Robust Iterative Method for Symmetric Quantum Signal Processing in All Parameter Regimes

Here, this paper addresses the problem of solving nonlinear systems in the context of symmetric quantum signal processing (QSP), a powerful technique for implementing matrix functions on quantum computers. Symmetric QSP focuses on representing target polynomials as products of matrices in SU(2) that possess symmetry properties. We present a novel Newton’s method tailored for efficiently solving the nonlinear system involved in determining the phase factors within the symmetric QSP framework. Our method demonstrates rapid and robust convergence in all parameter regimes, including the challenging scenario with ill-conditioned Jacobian matrices, using standard double precision arithmetic operations. For instance, solving symmetric QSP for a highly oscillatory target function α cos(1000x) (polynomial degree ≈ 1433) takes 6 iterations to converge to machine precision when α = 0.9, and the number of iterations only increases to 18 iterations when α = 1 – 10 -9 with a highly ill-conditioned Jacobian matrix. Leveraging the matrix product state structure of symmetric QSP, the computation of the Jacobian matrix incurs a computational cost comparable to a single function evaluation. Moreover, we introduce a reformulation of symmetric QSP using real-number arithmetics, further enhancing the method’s efficiency. Extensive numerical tests validate the effectiveness and robustness of our approach, which has been implemented in the QSPPACK software package.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Scalable Circuits for Preparing Ground States on Digital Quantum Computers: The Schwinger Model Vacuum on 100 Qubits

The vacuum of the lattice Schwinger model is prepared on up to 100 qubits of IBM’s Eagle-processor quantum computers. A new algorithm to prepare the ground state of a gapped translationally invariant system on a quantum computer is presented, which we call “scalable circuits ADAPT-VQE” (SC-ADAPT-VQE). This algorithm uses the exponential decay of correlations between distant regions of the ground state, together with ADAPT-VQE, to construct quantum circuits for state preparation that can be scaled to arbitrarily large systems. These scalable circuits can be determined with use of classical computers, avoiding the challenging task of optimizing parameterized circuits on a quantum computer. SC-ADAPT-VQE is applied to the Schwinger model, and is shown to be systematically improvable, with an accuracy that converges exponentially with circuit depth. Both the structure of the circuits and the deviations of prepared wave functions are found to become independent of the number of spatial sites, L . This allows a controlled extrapolation of the circuits, determined with use of small or modest-sized systems, to arbitrarily large L . The circuits for the Schwinger model are determined on lattices up to L = 14 (28 qubits) with the Qiskit classical simulator, and are subsequently scaled up to prepare the L = 50 (100 qubits) vacuum on IBM’s 127-superconducting-qubit quantum computers ibm_brisbane and ibm_cusco. After introduction of an improved error-mitigation technique, which we call “operator decoherence renormalization”, the chiral condensate and charge-charge correlators obtained from the quantum computers are found to be in good agreement with classical matrix product state simulations. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1991

Application of Krylov exponential propagation to fluid dynamics equations

This paper presents an application of matrix exponentiation via Krylov subspace projection, to the solution of fluid dynamics problems. The main idea is to approximate the operation exp(A)v by means of a projection-like process onto a Krylov subspace. This results in a computation of an exponential matrix vector product similar to the one above but of a much smaller size. Time integration schemes can then be devised to exploit this basic computational kernel. The motivation of this approach is to provide time-integration schemes that are essentially of an explicit nature but which have good stability properties.

Saad, Y.↗
NASA NTRS1991

Application of Krylov exponential propagation to fluid dynamics equations

An application of matrix exponentiation via Krylov subspace projection to the solution of fluid dynamics problems is presented. The main idea is to approximate the operation exp(A)v by means of a projection-like process onto a krylov subspace. This results in a computation of an exponential matrix vector product similar to the one above but of a much smaller size. Time integration schemes can then be devised to exploit this basic computational kernel. The motivation of this approach is to provide time-integration schemes that are essentially of an explicit nature but which have good stability properties.

Saad, Youcef↗
DOE OSTI2022

Matrix Multiply Performance of GPUs on Exascale-class HPE/Cray Systems

The computation of dense matrix-matrix products (GEMMs) is central to many modeling and simulation workloads as well as AI/ML deep learning campaigns. In fact, millions of dollars are spent annually on computing GEMMs, and large model training demands are increasing exponentially. Specialized processors such as GPUs are designed to perform well for these operations. However, the performance of GEMMs on GPUs can exhibit complex behaviors depending on many factors, making it challenging to optimize the performance of GEMMs on these processors. In this study we undertake an examination of GEMM performance on several leading GPU models taken from product lines of GPUs to be deployed in forthcoming exascale computing systems. We show results to illustrate the many factors that can affect performance of GEMMs on GPUs. We then present data collected from a large number of test runs for an example GEMM operation to show the dependence behaviors of GEMM rate on matrix dimensions. Finally, we show results from machine learning-based performance models using novel feature engineering methods to fit the measured performance, providing a potential basis for GEMM performance tuning and autotuning methods for GPUs. Recommendations are also given for how to achieve high GEMM performance on modern GPUs.

Melesse Vergara, Veronica↗
DOE OSTI2022

Characterization of the Emissions and Crystalline Silica Content of Airborne Dust Generated from Grinding Natural and Engineered Stones

Abstract In this study, we systematically characterized the airborne dust generated from grinding engineered and natural stone products using a laboratory testing system designed and operated to collect representative respirable dust samples. Four stone samples tested included two engineered stones consisting of crystalline silica in a polyester resin matrix (formulations differed with Stones A having up to 90wt% crystalline silica and Stone B up to 50wt% crystalline silica), an engineered stone consisting of recycled glass in a cement matrix (Stone C), and a granite. Aerosol samples were collected by respirable dust samplers, total dust samplers, and a Micro-Orifice Uniform Deposit Impactor. Aerosol samples were analyzed by gravimetric analysis and x-ray diffraction to determine dust generation rates, crystalline silica generation rates, and crystalline silica content. Additionally, bulk dust settled on the floor of the testing system was analyzed for crystalline silica content. Real-time particle size distributions were measured using an Aerodynamic Particle Sizer. All stone types generated similar trimodal lognormal number-weighted particle size distributions during grinding with the most prominent mode at an aerodynamic diameter of about 2.0-2.3 μm, suggesting dust formation from grinding different stones is similar. Bulk dust from Stone C contained no crystalline silica. Bulk dust from Stone A, Stone B, and granite contained 60, 23, and 30wt% crystalline silica, respectively. In Stones A and B, the cristobalite form of crystalline silica was more plentiful than the quartz form. Only the quartz form was detected in granite. The bulk dust, respirable dust, and total dust for each stone had comparable amounts of crystalline silica, suggesting that crystalline silica content in the bulk dust could be representative of that in respirable dust generated during grinding. Granite generated more dust per unit volume of material removed than the engineered stones, which all had similar normalized dust generation rates. Stone A had the highest normalized generation rates of crystalline silica, followed by granite, Stone B, and Stone C (no crystalline silica), which likely leads to the same trend of respirable crystalline silica (RCS) exposure when working with these different stones. Manufacturing and adoption of engineered stone products with formulations such as Stone B or Stone C could potentially lower or eliminate RCS exposure risks. Combining all the effects of dust generation rate, size-dependent silica content, and respirable fraction, the highest normalized generation rate of RCS consistently occurs at 3.2-5.6 µm for all the stones containing crystalline silica. Therefore, removing particles in this size range near the generation sources should be prioritized when developing engineering control measures.

Public, Environmental & Occupational Health↗
DOE OSTI2024

System Modeling and Robust Stability Region Analysis for Multi-inverters Based on VSG

As the number of inverter-interfaced distributed generators (IIDGs) increases, the decrease in the inertia of the microgrid affects the stability of microgrid voltage and current. The virtual synchronous generator (VSG) control method simulates the output characteristics of the synchronous machine to improve the inertia and damping of the microgrid. However, low-frequency oscillations can occur in systems with multi-VSGs connected in parallel, primarily when operating in island mode. Therefore, this article develops the small-signal state-space model for multi-inverters based on VSG and analyses the robust stability region (RSR) of the control parameters. Firstly, an extensible small-signal state-space model is established and simplified in the dq axis, which is a necessary preprocessing for the subsequent stability analysis. Secondly, the guardian map method based on bialternate product operation is constructed, and the RSR is established through translation mapping and rotation mapping. Therefore, the solution of RSR is transformed into the problem of identifying the Hurwitz matrix. Finally, the simulations and experiments verify that the calculated region can ensure stable operation during startup and meet the requirement to suppress low-frequency oscillations. In conclusion, the method of solving RSR proposed in this paper provides a basis for selecting controller parameters in practical engineering.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1995

Numerical mixing calculations of confined reacting jet flows in a cylindrical duct

The results reported in this paper describe some of the main flow characteristics and NOx production results which develop in the mixing process in a constant cross-sectional cylindrical duct. A 3-dimensional numerical model has been used to predict the mixing flow field and NOx characteristics in a mixing section of an RQL combustor. Eighteen configurations have been analyzed in a circular geometry in a fully reacting environment simulating the operating condition of an actual RQL gas turbine combustion liner. The evaluation matrix was constructed by varying three parameter: (1) jet-to-mainstream momentum-flux ration (J), (2) orifice shape or orifice aspect ratio, and (3) slot slant angle. The results indicate that the mixing flow field and NOx production significantly vary with the value of the jet penetration and subsequently, slanting elongated slots generally improve the NOx production at high J conditions. Round orifices produce low NOx at low J due to the strong jet penetration. The NOx production trends do not correlate with the mixing non-uniformity parameters described herein.

Oechsle, Victor L.↗
DOE OSTI2026

Surface-Initiated Atom Transfer Radical Polymerization Using Hydrogel Reactors

Atom transfer radical polymerization (ATRP) is a controlled radical polymerization method that enables the synthesis of tailored polymeric materials with low dispersity, highlighting its immense potential for green fabrication of advanced materials. However, its broader implementation is limited by challenges in product isolation, maintaining catalyst activity, and mitigating atmospheric sensitivity arising from oxygen-sensitive metal catalysts. Here, gelatin hydrogels (GHs) are introduced as a soft “reactor” matrix for interfacial ATRP, operating with minimal metal-catalyst loading while exhibiting possibly an organoreductive behavior. This strategy leverages activator regeneration via electron transfer through a ligand–metal charge-transfer (LMCT) mechanism to reduce oxidized metal catalysts within the GH network. Polymerization is evaluated by growing polymer brushes at an active interface formed between GHs swollen in monomer solution and an initiating surface, and sequential growth experiments confirmed that GH-mediated ATRP preserves living character. Under UV illumination, LMCT is activated, producing polymers both at the desired interface and within the GH bulk. UV–Vis spectroscopy revealed active reduction of Cu(II) to Cu(I) along with concentration-dependent complex formation, indicating dynamic coordination chemistry within the hydrogel. The redox-active arginine- and glutamic acid-rich gelatin backbone coordinates and reduces the metal center, enabling ATRP at ppm-level catalyst concentrations. While polymerization proceeds in GH-Cu(II) reactors, adding external mobile ligands to the GH results in longer polymer brushes. Here, the results reported here are exploratory. More experiments are needed to characterize polymer brush growth in GHs and compare it to conventional surface-initiated polymerization in solution.

Absorption↗
DOE OSTI2023

Influence of extreme fracture flow channels on the thermal performance of open-loop geothermal systems at commercial scale

Abstract Adequate stewardship of geothermal resources requires accurate forecasting of long-term thermal performance. In enhanced geothermal systems and other fracture-dominated reservoirs, predictive models commonly assume constant-aperture fractures, although spatial variations in aperture can greatly affect reservoir permeability, fluid flow distribution, and heat transport. Whereas previous authors have investigated the effects of theoretical random aperture distributions on thermal performance, here we further explore the influence of permeability heterogeneity considering field-constrained aperture distributions from a meso-scale field site in northern New York, USA. Using numerical models of coupled fluid flow and heat transport, we conduct thermal–hydraulic simulations for a hypothetical reservoir consisting of a relatively impervious porous matrix and a single, horizontal fracture. Our results indicate that in highly channelized fields, most well design configurations and operating conditions result in extreme rates of thermal drawdown (e.g., 50% drop in production well temperatures in under 2 years). However, some other scenarios that account for the risks of short-circuiting can potentially enhance heat extraction when mass flow rate is not excessively high, and the direction of geothermal extraction is not aligned with the most permeable features in the reservoir. Through a parametric approach, we illustrate that well separation distance and relative positioning play a major role in the long-term performance of highly channelized fields, and both can be used to help mitigate premature thermal breakthrough.

Rangel-Jurado, Nicolás↗
NASA NTRS2010

NASA Tech Briefs, December 2010

Topics include: Coherent Frequency Reference System for the NASA Deep Space Network; Diamond Heat-Spreader for Submillimeter-Wave Frequency Multipliers; 180-GHz I-Q Second Harmonic Resistive Mixer MMIC; Ultra-Low-Noise W-Band MMIC Detector Modules; 338-GHz Semiconductor Amplifier Module; Power Amplifier Module with 734-mW Continuous Wave Output Power; Multiple Differential-Amplifier MMICs Embedded in Waveguides; Rapid Corner Detection Using FPGAs; Special Component Designs for Differential-Amplifier MMICs; Multi-Stage System for Automatic Target Recognition; Single-Receiver GPS Phase Bias Resolution; Ultra-Wideband Angle-of-Arrival Tracking Systems; Update on Waveguide-Embedded Differential MMIC Amplifiers; Automation Framework for Flight Dynamics Products Generation; Product Operations Status Summary Metrics; Mars Terrain Generation; Application-Controlled Parallel Asynchronous Input/Output Utility; Planetary Image Geometry Library; Propulsion Design With Freeform Fabrication (PDFF); Economical Fabrication of Thick-Section Ceramic Matrix Composites; Process for Making a Noble Metal on Tin Oxide Catalyst; Stacked Corrugated Horn Rings; Refinements in an Mg/MgH2/H2O-Based Hydrogen Generator; Continuous/Batch Mg/MgH2/H2O-Based Hydrogen Generator; Strain System for the Motion Base Shuttle Mission Simulator; Ko Displacement Theory for Structural Shape Predictions; Pyrotechnic Actuator for Retracting Tubes Between MSL Subsystems; Surface-Enhanced X-Ray Fluorescence; Infrared Sensor on Unmanned Aircraft Transmits Time-Critical Wildfire Data; and Slopes To Prevent Trapping of Bubbles in Microfluidic Channels.

Source record↗
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↗