Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear”

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 73 records · Page 4

Learning linear optical circuits with coherent states

We analyze the energy and training data requirements for supervised learning of an M-mode linear optical circuit by minimizing an empirical risk defined solely from the action of the circuit on coherent states. When the linear optical circuit acts non-trivially only on k < M unknown modes (i.e. a linear optical k-junta), we provide an energy-efficient, adaptive algorithm that identifies the junta set and learns the circuit. We compare two schemes for allocating a total energy, E, to the learning algorithm. In the first scheme, each of the T random training coherent states has energy E/T. In the second scheme, a single random MT-mode coherent state with energy E is partitioned into T training coherent states. The latter scheme exhibits a polynomial advantage in training data size sufficient for convergence of the empirical risk to the full risk due to concentration of measure on the $(2MT-1)$-sphere. Specifically, generalization bounds for both schemes are proven, which indicate that for ε-approximation of the full risk by the empirical risk with high probability, $O(E^{2/3}M^{2/3}/\epsilon^{2/3})$ training states are sufficient for the first scheme and $O(E^{1/3}M^{1/3}/\epsilon^{2/3})$ training states are sufficient for the second scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sparse Linear Solvers for Large-scale Electromagnetic Transient Simulations

Linear solvers form the basis for electromagnetic transient (EMT) simulations. There is a need to speed up EMT simulations as larger regions are analyzed using EMT simulations. For the same, the performance of linear solvers plays an important role. Exploiting the sparsity of the matrices generated in EMT simulations could assist with speed-up. Scalability is also crucial as power grids expand, demanding solutions capable of accommodating the increasing system size. Recent studies from the North American Electric Reliability Corporation (NERC) increasingly emphasize that EMT simulation models of the power grid will grow larger with the inclusion of power electronics components. Parallelisms in sparsity patterns exploit modern central processing units (CPUs), multi-core CPUs, and graphics processing units (GPUs) architectures in sparse solver designs. Therefore, this paper explores publicly available existing linear solvers and investigates their efficiency in large-scale power grid simulations. A large-scale power grid is developed by increasing the size of the IEEE 39 bus test system to up to 39000 bus systems.

Hsu, Kuan-Chieh↗

A Linear Programming Approach to Backtracking for Single-Axis Trackers on Rolling Terrain

In this article, we present a computationally efficient method for determining optimal backtracking rotations for single-axis solar trackers on nonuniform terrain. The method allows for ganged tracking, mechanical rotation constraints, uneven row spacing, and arbitrary maximum allowable shaded fractions (to enable “fractional backtracking”). As with previous 2-D approaches, the method is suitable for terrain that varies in the transverse direction with respect to the rotation axis of the trackers. The novelty of the method lies in formulating the problem of shade avoidance as a linear problem, which is achieved by using the row interception width as the optimization variable instead of rotation angles. Formulating backtracking as a linear problem enables the use of extremely efficient linear programming algorithms, making the method highly scalable, requiring less than 1 min to compute optimal rotation schedules for hundreds of trackers. It also produces more effective backtracking rotations, reducing the frequency of shading by 4× and improving system energy output by 1%–2%.

Optimization↗

Early Exploration of a Flexible Framework for Efficient Quantum Linear Solvers in Power Systems

The rapid integration of renewable energy resources presents formidable challenges in managing power grids. While advanced computing and machine learning techniques offer some solutions for accelerating grid modeling and simulation, there remain complex problems that classical computers cannot effectively address. Quantum computing, a promising technology, has the potential to fundamentally transform how we manage power systems, especially in scenarios with a higher proportion of renewable energy sources. One critical aspect is solving linear systems of equations, crucial for power system applications like power flow analysis, for which the Harrow-Hassidim-Lloyd (HHL) algorithm is a well-known quantum solution. However, HHL quantum circuits often exhibit excessive depth, making them impractical for current Noisy-Intermediate-Scale-Quantum (NISQ) devices. In this paper, we introduce a versatile framework, powered by NWQSim, that bridges the gap between power system applications and quantum linear solvers available in Qiskit. This framework empowers researchers to efficiently explore power system applications using quantum linear solvers. Through innovative gate fusion strategies, reduced circuit depth, and GPU acceleration, our simulator significantly enhances resource efficiency. Power flow case studies have demonstrated up to a eight-fold speedup compared to Qiskit Aer, all while maintaining comparable levels of accuracy.

quantum computing, Harrow-Hassidim-Lloyd, high-per↗

Limitations of Fault-Tolerant Quantum Linear System Solvers for Quantum Power Flow

Quantum computers hold promise for solving problems intractable for classical computers, especially those with high time or space complexity. Practical quantum advantage can be said to exist for such problems when the end-to-end time for solving such a problem using a classical algorithm exceeds that required by a quantum algorithm. Reducing the power flow (PF) problem into a linear system of equations allows for the formulation of quantum PF (QPF) algorithms, which are based on solving methods for quantum linear systems such as the Harrow-Hassidim-Lloyd (HHL) algorithm. Speedup from using QPF algorithms is often claimed to be exponential when compared to classical PF solved by state-of-the-art algorithms. Here, we investigate the potential for practical quantum advantage in solving QPF compared to classical methods on gate-based quantum computers. Notably, this paper does not present a new QPF solving algorithm but scrutinizes the end-to-end complexity of the QPF approach, providing a nuanced evaluation of the purported quantum speedup in this problem. Our analysis establishes a best-case bound for the HHL-based quantum power flow complexity, conclusively demonstrating that the HHL-based method has higher runtime complexity compared to the classical algorithm for solving the direct current power flow (DCPF) and fast decoupled load flow (FDLF) problem. Notably, our analysis and conclusions can be extended to any quantum linear system solver with rigorous performance guarantees, based on the known complexity lower bounds for this problem. Additionally, we establish that for potential practical quantum advantage (PQA) to exist it is necessary to consider DCPF-type problems with a very narrow range of condition number values and readout requirements.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Linear viscoelastic measurements of commerical therapy putties for design and intuition

This dataset complements the publication Mapping Linear Viscoelasticity for Design and Tactile Intuition (Corman & Ewoldt, 2019). It contains rigorous shear rheometry measurements for six commercial silicone putties, including creep recovery, SAOS, LAOS, and stress relaxation. Additional, unpublished data was generated as part of a graduate rheology course taught by Randy Ewoldt at University of Illinois Urbana-Champaign. Together, these datasets provide a comprehensive example of linear viscoelasticity of a commercial material. The collection is intended both as a pedagogical resource for developing physical intuition in linear viscoelasticity and as a reference dataset for data interpretation, model fitting, and comparative rheological analysis.

C. Marsh, Maxwell↗

Randomized Algorithms for Linear Solvers

Recently, randomized algorithms in numerical linear algebra, specifically those centered around random sketching, have gained traction in primarily theoretical research due to their potential to significantly reduce problem dimensionality at the cost of an O(1) multiplicative distortion factor. It has been assumed that this sketching can be done efficiently, but thorough investigation into how precisely to do it has been neglected. Moreover, the theory-based community has argued for sketching’s ability to reduce computational cost via complexity analysis, but has not researched how it affects the stability of the algorithms. At Sandia, efficient linear solvers that scale well on modern HPC architectures while maintaining stability are imperative for practical applications. In this LDRD, we developed a random sketching strategy that is substantially faster than existing ones, and demonstrate its superior performance in practice on a NVIDIA H100 GPU. Moreover, we show how this can be used to significantly outperform existing linear least squares solvers while improving the solver’s stability as well. Additionally, we demonstrate how this sketching strategy can be used to make a fast, stable QR factorization that can subsequently be used in s-step and block Krylov solvers. Finally, we incorporate a sketching-based block orthogonalization scheme into s-step GMRES, which is stable and faster than existing approaches on the Perlmutter supercomputer.

97 MATHEMATICS AND COMPUTING↗

Compensating for Amplifier Non-Linearity in a SEL Controller

The SEL architecture used commonly in CW superconducting linacs is a positive feedback system that requires determining the correct limits for the operating conditions to ensure stable operation. These limits are determined from the results of the amplifier calibration which characterizes the dac drive(forward power) against the cavity field. When the amplifier non-linearity precludes a simple linear fit for this characteristic, its modification to account for the amplifiers non-linearity is essential to operate the system in the SEL architecture. This process is described with test data from a 32kW SSA for the PIP-II linac.

Sankar Raman, Shrividhyaa [Fermilab]↗

Geometric Interpretation of a Non-Linear Extension of Quantum Mechanics

We recently introduced a particular non-linear generalization of quantum mechanics that has the property that it is exactly solvable in terms of the eigenvalues and eigenfunctions of the Hamiltonian of the usual linear quantum mechanics problem. In this paper, we suggest that the two components of the wave function represent the system described by the Hamiltonian H in two different asymptotic regions of spacetime and we show that the non-linear terms can be viewed as giving rise to gravitational effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

ESPPU INPUT: C$^3$ within the "Linear Collider Vision"

The Linear Collider Vision calls for a Linear Collider Facility with a physics reach from a Higgs Factory to the TeV-scale with $e^+e^{-}$ collisions. One of the technologies under consideration for the accelerator is a cold-copper distributed-coupling linac capable of achieving high gradient. This technology is being pursued by the C$^3$ collaboration to understand its applicability to future colliders and broader scientific applications. In this input we share the baseline parameters for a C$^3$ Higgs-factory and the energy reach of up to 3 TeV in the 33 km tunnel foreseen under the Linear Collider Vision. Recent results, near-term plans and future R&D needs are highlighted.

43 PARTICLE ACCELERATORS↗

Status of the International Linear Collider

This paper is not a proposal for a CERN future project but provides information on the International Linear Collider (ILC) considered for Japan in order to facilitate the European Strategy discussion in a global context. It describes progress to date, ongoing engineering studies, updated cost estimate for the machine at $\sqrt{s}=250~\rm GeV$ and the situation in Japan. The physics of the ILC is not presented here, but jointly for all Linear Collider projects in a separate document ``A Linear Collider Vision for the Future of Particle Physics'' submitted for the forthcoming European Strategy deliberations.

Abe, Y. [KEK, Tsukuba]↗

Compensating for Amplifier Non-Linearity in a SEL Controller

The Self-Excited Loop (SEL) architecture, used in some continuous-wave (CW) superconducting linacs, relies on a positive feedback mechanism that requires carefully defined operating limits to ensure stable operation. These limits are typically derived from amplifier calibration, which characterizes the relationship between forward power and DAC drive. However, amplifier non-linearity often prevents a simple linear fit of this characteristic, introducing errors that can compromise stability. To address this, we present a modified calibration procedure that incorporates amplifier non-linearity into the SEL framework. The approach is validated with test data from a 32 kW solid-state amplifier (SSA) and a cavity emulator developed for the Fermilab PIP-II linac.

Raman, S. [Fermilab] (ORCID:0009000735937206)↗

Rational Design of Methylated Triazine‐Based Linear Conjugated Polymers for Efficient CO 2 Photoreduction with Water

The development of semiconducting conjugated polymers for photoredox catalysis holds great promise for sustainable utilization of solar energy. In this work, a new family of porous methylated triazine‐based linear conjugated polymers is reported that enable efficient photoreduction of carbon dioxide (CO 2 ) with water (H 2 O) vapor, in the absence of any additional photosensitizer, sacrificial agents or cocatalysts. It is demonstrated that the key lies in the generation of methylated triazine linkages through a facile condensation reaction between benzamidine and acetic anhydride, which impedes the formation of conventional triazine‐based frameworks. It is also shown that regulating conjugated linear backbones with different lengths of electron‐donated benzyl units provides a facile means to modulate their optical properties and the exciton dissociation, thereby affording more long‐lived photogenerated charge carriers and boosting charge separation and transfer. A high‐performance carbon monoxide (CO) production rate of 218.9 µmol g −1 h −1 is achieved with ≈ 100% CO selectivity, which is accompanied by exceptional H 2 O oxidation to oxygen (O 2 ). It anticipates this new study will advance synthetic approaches toward polymeric semiconductors and facilitate new possibilities for triazine‐based conjugated polymers with promising potential in artificial photocatalysis.

36 MATERIALS SCIENCE↗

Solution of the linear wave-particle kinetic equation for global modes of arbitrary frequency in a tokamak

The linear response of a plasma to perturbations of arbitrary frequency and wavelength is derived for any axisymmetric magnetized toroidal plasma. An explicit transformation to action-angle coordinates is achieved using orthogonal magnetic coordinates and the Littlejohn Lagrangian, establishing the validity of this result to arbitrary order in normalized Larmor radius. The global resonance condition for compressional modes is clarified in more detail than in previous works, confirming that the poloidal orbit-average of the cyclotron frequency gives the desired result at lowest order in Larmor radius. The global plasma response to the perturbation at each resonance is captured by a poloidal and gyroaverage of the perturbing potential. A “global gyroaveraging” of the potential is a natural by-product of this analysis which takes into account the changing of the magnetic field over an orbit. The resonance condition depends on two arbitrary integers which completely separately capture the effects poloidal non-uniformity and finite Larmor radius in generating sidebands. We learn that poloidal sidebands generated for compressional modes are dominated by the change in gyrofrequency over the orbit, which is very different to shear modes where the gyrofrequency only contributes via a finite Larmor radius effect. This increases the number of bounce harmonics required to compute the linear drive, giving a more complicated resonance map. An example calculation is given comparing resonance of shear and compressional modes in a published DIII-D case.

Compressional↗

RandONets: Shallow networks with random projections for learning linear and nonlinear operators

Deep neural networks have been extensively used for the solution of both the forward and the inverse problem for dynamical systems. However, their implementation necessitates optimizing a high-dimensional space of parameters and hyperparameters. This fact, along with the requirement of substantial computational resources, pose a barrier to achieving high numerical accuracy, but also interpretability. Here, to address the above challenges, we present Random Projection-based Operator Networks (RandONets): shallow networks with random projections and tailor-made numerical analysis methods that learn accurately and fast linear and nonlinear operators. Building on previous works, we prove that RandOnets are universal approximators of linear and nonlinear operators. Due to their simplicity, RandONets provide a one-step transformation of the input space, facilitating interpretability. For the evaluation of their performance, we focus on operators of PDEs. We show, that RandONets outperform by several orders of magnitude, both in terms of numerical approximation accuracy and computational cost, the “vanilla” DeepONets. Hence, we believe that our method will trigger further developments in the field of scientific machine learning, for the development of new ‘’light”schemes that will provide high accuracy while reducing dramatically the computational cost. A MATLAB toolbox for RandONets, including demos, is available on GitHub at https://github.com/GianlucaFabiani/RandONets.

Interpretable machine learning↗

Summary of Working Group 7: Linear colliders

The contributions presented during the Working Group 7: Linear Colliders (WG7) sessions at the 2024 Advanced Accelerator Concepts Workshop are briefly summarized, as well as the discussions regarding a design initiative for a linear collider based on wakefield accelerator technology.

Advanced Accelerator Concepts↗

γIn-beam angular distribution and linear polarization measurements with GRETINA using a simple energy-ordering approach

Angular distribution and linear polarization measurements are powerful tools for inferring the spins and parities of nuclear levels. In this work, the performance of the Gamma-Ray Energy Tracking In-beam Nuclear Array (GRETINA) as a Compton polarimeter was characterized in a fusion-evaporation reaction experiment using a simple energy-ordering approach for the interaction points assigned in the signal decomposition process. A variety of multipolarities and characters for γ-ray transitions in the reaction products 25 Mg, 25 Na, and 22 Ne, formed from fusion-evaporation of an 18 O beam on a 9 Be target, were examined. The experimental angular distributions and linear polarization asymmetries were consistent with predictions using the theoretical formalism accounting for the Lorentz boost.

Angular distribution↗