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

Adaptive Linear State Estimation for Unbalanced Distribution System

The inclusion of PMU functionality in distribution relays enables the implementation of a linear state estimator (LSE) in Distribution Systems (DS). However, the unbalanced topology and phase coupling in distribution lines necessitate modifications to the LSE formulation. Additionally, the higher fault frequency in distribution systems requires a state estimation approach that is resilient to contingencies. This work proposes an adaptive linear state estimation algorithm tailored for unbalanced distribution systems with single-phase and two-phase laterals. Furthermore, a modified Optimal PMU Placement (OPP) strategy is introduced to ensure full observability in distribution systems with single-phase and two-phase buses. To maintain adaptability to topology changes, the state estimator incorporates circuit breaker status data provided by PMUs, ensuring robust performance during topology changes triggered by faults. The performance of the algorithm is verified on the IEEE 13-bus, 34-bus, and 123-bus systems.

PMUs↗
DOE OSTI2025

Linear Solver for Electromagnetic Simulation of General Distribution Feeders

High-fidelity electromagnetic transient (EMT) modeling is required for accurate simulation and analysis of power system dynamics in modern distribution feeders. However, the high-fidelity of EMT models often leads to significant computational challenges, particularly in terms of computational resources and simulation time. This paper investigates the development and application of a detailed EMT model for general distribution feeders, with a focus on improving computational efficiency. A direct linear solver is proposed for a bordered block diagonal (BBD) matrix structure commonly encountered in a EMT model of distribution feeders. The solver integrates the Schur complement method with the block tridiagonal matrix algorithm to enhance the computational performance. The proposed solver is validated using the primary feeder of the IEEE 342-node test system, demonstrating its accuracy and efficiency in EMT simulations. Furthermore, the solver’s performance is benchmarked against MATLAB’s built-in linear solvers, showing significant improvements in computation time while maintaining high fidelity and accuracy in simulation results.

Choi, Jongchan [ORNL] (ORCID:000000025952455X)↗
DOE OSTI2024

Losing Control of Your Linear Network? Try Resilience Theory

Resilience of cyber-physical networks to unexpected failures is a critical need widely recognized across domains. For instance, power grids, telecommunication networks, transportation infrastructures, and water treatment systems have all been subject to disruptive malfunctions and catastrophic cyberattacks. Following such adverse events, we investigate scenarios where a node of a linear network suffers a loss of control authority over some of its actuators. These actuators are not following the controller's commands and are instead producing undesirable outputs. The repercussions of such a loss of control can propagate and destabilize the whole network despite the malfunction occurring at a single node. To assess system vulnerability, we establish resilience conditions for networks with a subsystem enduring a loss of control authority over some of its actuators. Furthermore, we quantify the destabilizing impact on the overall network when such a malfunction perturbs a nonresilient subsystem. We illustrate our resilience conditions on two academic examples, on an islanded microgrid and on the linearized IEEE 39-bus system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Conditional Sixth-Harmonic Injection to Improve the Linear Modulation Range of Three-Phase Voltage-Source Converters

Overmodulation techniques are used to improve dc bus voltage utilization for three-phase voltage-source converters (VSCs), resulting in significant reduction in fundamental voltage gain. Conventional pulse width modulation (PWM) techniques for three-phase VSCs cannot enhance the linear range of dc bus utilization by more than 15.5%. Here, in this letter, a conditional sixth-harmonic injection technique is introduced, which can theoretically improve the linear modulation range by up to 19.2% for three-phase ac systems without additional hardware components. Experimental validation of the proposed scheme is presented, with its performance compared against the conventional third-harmonic PWM technique.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Research and Development to Reduce Impurity Production and Transport of the Impurities to the Target in Linear Plasma Devices Using Helicon Plasma Sources

Linear plasma devices used to test plasma facing materials (PFMs) and components for fusion reactors are often suffering under the production of intrinsic impurities from the plasma source system. Most linear plasma devices use internal electrodes (hollow cathodes, reflex arc, or cascading arc), which typically are the source of the main impurities. The next-generation plasma generators use radio frequency (RF) plasma sources like helicons to avoid internal electrodes. However, high power operation of helicons has proven to result in impurity production due to the high rectified sheath voltages created. Depending on the plasma parameters and magnetic configuration, these impurities can be transported to the target and deposited there to unacceptable high levels. Here, in this contribution, the experimental results from Proto-MPEX are summarized, and the conclusions of the impurity source physics are given. Methods to reduce the impurity production, the impurity transport, and the net deposition on the target are presented. These methods to reduce the impurity production include Faraday screens to reduce the sheath voltage drop, high-Z refractory coatings to reduce the erosion yield, and wall conditioning methods. Methods to reduce the impurity transport include changes in the magnetic configuration as well as electron heating to change axial and radial temperature profiles. Preliminary results on the effectiveness of some of these methods are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Novel Approach for Computing Rigid Body Motion Using Linear Accelerations

Here, a novel approach is presented for computing general rigid body motion based on a few known linear accelerations. This method utilizes linear acceleration data obtained from three distinct points on the body, all within a body-fixed reference frame. The only requirement is that the three chosen points must not be collinear. A system of differential-algebraic equations is derived, combining principles of rigid body kinematics with theory of the rotation group SO(3). These equations provide a framework for numerically computing various motion parameters, including angular velocity, angular acceleration, body orientation, velocity field, acceleration field, and displacement field. By numerically solving this system of equations, we can fully characterize rigid body motion in three-dimensional space. A numerical example is provided to demonstrate the practical implementation and efficacy of the proposed technique, illustrating its potential for accurate motion computation in various applications.

42 ENGINEERING↗
DOE OSTI2025

Fast and Scalable FFT-Based GPU-Accelerated Algorithms for Block-Triangular Toeplitz Matrices with Application to Linear Inverse Problems Governed by Autonomous Dynamical Systems

In this work, we present an efficient and scalable algorithm for performing matrix-vector multiplications (matvecs) for block Toeplitz matrices. Such matrices, which are shift-invariant with respect to their blocks, arise in the context of solving inverse problems governed by autonomous systems, and time-invariant systems in particular. In this article, we consider inverse problems that infer unknown parameters from observational data of a linear time-invariant dynamical system given in the form of partial differential equations (PDEs). Matrix-free Newton-conjugate-gradient methods are often the gold standard for solving these inverse problems, but they require numerous actions of the Hessian on a vector. Matrix-free adjoint-based Hessian matvecs require solution of a pair of linearized forward/adjoint PDE solves per Hessian action, which may be prohibitive for large-scale inverse problems. Time invariance of the forward PDE problem leads to a block Toeplitz structure of the discretized parameter-to-observable (p2o) map defining the mapping from inputs (parameters) to outputs (observables) of the PDEs. This block Toeplitz structure enables us to exploit two key properties: (1) compact storage of the p2o map and its adjoint, and (2) efficient fast Fourier transform–based Hessian matvecs. The proposed algorithm is mapped onto large multi-GPU clusters and achieves more than 80% of peak bandwidth on NVIDIA A100 GPUs. Excellent weak scaling is shown for up to 48 A100 GPUs. For the targeted problems, the implementation executes Hessian matvecs within fractions of a second, which is orders of magnitude faster than can be achieved by conventional matrix-free Hessian matvecs via forward/adjoint PDE solves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Ipopt Interface to Re::Solve Linear Solver

The software provides Ipopt optimization package an interface to the Re::Solve linear solver library. Re::Solve features GPU-resident direct and iterative linear solvers that could be used to accelerate optimization computations.

Alam, Maksudul [Oak Ridge National Laboratory (ORN↗
DOE OSTI2024

LAPIS: Linear Algebra Performance for Intermediate Subprograms

SAND2025-11594O LAPIS (Linear Algebra Performance for Intermediate Subprograms) is a compiler infrastructure for linear algebra that targets both high productivity and performance portability. It is based on the open-source MLIR package from the LLVM project. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Kelley, Brian↗
DOE OSTI2025

Bayesian Linear Regression for Hugoniot Data

This repository provides the code and datasets used in the paper Bayesian Analysis of Linear Shock Compression Data. This paper analyzes publicly available shock compression datasets on copper, argon, and nickel from Marsh (1980) using Bayesian linear regression, and compares the results with those obtained using bootstrapping methods. References: - Marsh, S. P. (1980). LASL shock Hugoniot data (Vol. 5). Univ of California Press.

Bernstein, JasonA [Lawrence Livermore National Lab↗
DOE OSTI2024

Consistent linearization of micromorphic continuum theories

We outline the procedure of consistent linearization and apply it to the micromorphic, microstretch, and micropolar theories of continua. This yields tractable linear theories for nonlinearly elastic microstructured materials undergoing finite deformations. The results may be readily utilized in computational mechanics, stability and bifurcation analyses, and small-deformation problems in the context of these types of continua. Our results generalize those existing in current literature and facilitate their recovery upon incorporating appropriate kinematic and constitutive assumptions.

McAvoy, Ryan C.↗
DOE OSTI2024

Evolution of the SLATE linear algebra library

SLATE (Software for Linear Algebra Targeting Exascale) is a distributed, dense linear algebra library targeting both CPU-only and GPU-accelerated systems, developed over the course of the Exascale Computing Project (ECP). While it began with several documents setting out its initial design, significant design changes occurred throughout its development. In some cases, these were anticipated: an early version used a simple consistency flag that was later replaced with a full-featured consistency protocol. In other cases, performance limitations and software and hardware changes prompted a redesign. Sequential communication tasks were parallelized; host-to-host MPI calls were replaced with GPU device-to-device MPI calls; more advanced algorithms such as Communication Avoiding LU and the Random Butterfly Transform (RBT) were introduced. Early choices that turned out to be cumbersome, error prone, or inflexible have been replaced with simpler, more intuitive, or more flexible designs. Applications have been a driving force, prompting a lighter weight queue class, nonuniform tile sizes, and more flexible MPI process grids. Of paramount importance has been building a portable library that works across several different GPU architectures – AMD, Intel, and NVIDIA – while keeping a clean and maintainable codebase. Here we explore the evolving design choices and their effects, both in terms of performance and software sustainability.

Gates, Mark↗
DOE OSTI2026

SEGUID v2: Extending SEGUID checksums for circular, linear, single- and double-stranded biological sequences

Background Synthetic biology involves combining different DNA fragments, each containing functional biological parts, to address specific problems. Fundamental gene-function research often requires cloning and propagating DNA fragments, such as those from the iGEM Parts Registry or Addgene, typically distributed as circular plasmids. Addgene’s repository alone offers around 150,000 plasmids. To ensure data integrity, cryptographic checksums can be calculated for the sequences. Each sequence has a unique checksum, making checksums useful for validation and quick lookups of associated annotations. For example, the SEGUID checksum uniquely identifies protein sequences with a 27-character string. Objectives The original SEGUID, while effective for protein sequences and single-stranded DNA (ssDNA), is not suitable for circular DNA since there is no natural starting position nor for double-stranded DNA (dsDNA) since two separate sequences are present. Challenges include how to uniquely represent linear dsDNA, circular ssDNA, and circular dsDNA. To meet these needs, we propose SEGUID v2, which extends the original SEGUID to handle additional types of sequences. Conclusions SEGUID v2 produces orientation and rotation invariant checksums for single-stranded, double-stranded, possibly staggered, linear, and circular DNA and RNA sequences. Customizable alphabets allow for other types of sequences. In contrast to the original SEGUID, which uses Base64, SEGUID v2 uses Base64url to encode the SHA-1 hash. This ensures SEGUID v2 checksums can be used as-is in filenames, regardless of platform, and in URLs, with minimal friction. Availability SEGUID v2 is readily available for major programming languages, distributed under the MIT license. JavaScript package seguid is available on npm, Python package seguid on PyPi, R package seguid on CRAN, and a Tcl script on GitHub. These tools, along with documentation, examples, and an online SEGUID Calculator , can be found at https://www.seguid.org .

Pereira, Humberto↗
DOE OSTI2024

Bringing randomized algorithms to mainstream numerical linear algebra

Numerical linear algebra (NLA) underpins huge swaths of computational science and engineering. For scientists and engineers to make the most of the DOE’s computing resources, it is essential that they have access to high-performance implementations of algorithms with best-in-class scalability and reliability. Despite this, prevailing NLA libraries have little to no support for breakthrough algorithms from the field of randomized numerical linear algebra (RandNLA) that have been developed over the past twenty years. The goal of this LDRD was to break a log-jam that had prevented broad adoption of RandNLA. Our work had two thrusts. The first was to develop RandBLAS: a trustworthy and high-performance C++ library for randomized dimension reduction (an operation widely known as sketching). The second was the development of a novel randomized algorithm for computing a challenging type of matrix decomposition known as Householder QR with column pivoting (Householder QRCP). In this one-year late-start LDRD we successfully delivered RandBLAS 1.0 and new CPU and GPU codes for Householder QRCP. RandBLAS has extensive documentation at https://randblas.readthedocs.io/en/stable/. Papers on RandBLAS and and our high-performance QRCP codes are forthcoming.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

High Linearity RF Circuits in CMOS Platforms (CRADA Final Report)

The goal of this CRADA was to understand high linearity silicon-on-insulator (SOI) radio frequency switches. The first direction investigated was high linearity RF mixers as these components are composed of a network of switches. During the course of the investigation, it was determined that while the SOI RF switches were effective in the application in mixers and produced state of the art performance, the use of SOI for mixers is not economically feasible because the technology only worked with advanced CMOS nodes <65nm which can have an NRE cost (from lithography masks) in excess of $\$$500,000.

42 ENGINEERING↗
DOE OSTI2025

Applying Linear Optics from Closed Orbits Modulation for finding beam-based alignment of harmonic sextupoles

A fast and accurate beam-based alignment (BBA) method for harmonic sextupoles has been devel oped at NSLS-II using Linear Optics from Closed Orbit Modulation (LOCOM). The approach excites the beam with simultaneous sine-wave signals at two fast correctors, chosen with an appropriate phase advance to span the full betatron phase space. This strategy suppresses systematic errors from hysteresis, while additional errors from orbit drift and power-supply calibration are minimized by the short measurement time (a few minutes) and reliance solely on beam-based current-to-field conver sion of the sextupoles, with hysteresis explicitly included. Simulations indicate that Linear Optics from Closed Orbits (LOCO) combined with 0.5 mm local orbit bumps can resolve relative sextupole field offsets (∆k₂) with precision better than 10% of k₂, reflecting to beam-based alignment (BBA) accuracy finer than 50 µm. Moreover, employing machine-learning-optimized local orbit bumps en hances sextupole-induced quadrupole signals in a deterministic manner and maintains orbit stability under large sextupole strength variations (±40%). The proposed method is experimentally validated through proof-of-principle measurements at NSLS-II, demonstrating its potential as a fast, precise, and robust tool for harmonic sextupole alignment.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Performance Testing of Enhanced Linear Variable Integrated Sensor (ELVIS III) for LVDTs

The Enhanced Linear Variable Intrinsic Sensor (ELVIS) was a 2024 breakthrough discovery by Idaho National Laboratory (INL) in terms of further innovating the usage of all linear variable differential transformers (LVDTs) and enhancing their performance in harsh environments. In irradiation tests conducted within material test reactors (MTRs), LVDTs with an internal temperature sensing capability can address several critical challenges. This report focuses on performance testing of the ELVIS III device with new hybrid LVDTs provided by the Institute for Energy Technology (IFE), which is the world’s sole supplier of nuclear-grade LVDTs. The ELVIS III device was evaluated in terms of temperature showcasing promising results comparable to a typical type-K thermocouple (TC).

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2024

Initial PIP-II Beam Current Monitor Fault Case Analyses & Beam Position Monitor Linearity Studies in CST Studio Suite

The use of non-invasive sensors & systems to measure particle beam characteristics is a crucial part of modern accelerator control systems due to their ability to return real time beam data while minimizing negative effects on beam quality. To ensure that one can be reasonably confident these sensors will behave as desired upon be-ing implemented within the beamline, simulations pre-dicting the performance of these sensors under beamline conditions can be used as a valuable tool for checking sensor functionality without a physical test bench. This paper details the design, testing, and results of two sensor models developed using CST Studio Suite soft-ware designed to mimic two sensors to be implemented within the PIP-III beamline: an elliptical, large-aperture beam position monitor (BPM) for which vertical & hori-zontal position signal linearity was analyzed, and an AC current transformer (ACCT) beam current monitor (BCM) used to search for potential fault cases within the BCM and beam pipe flange gaps. Special focus is given to the discovery of linearity variations within the BPM and the use of frequency domain techniques in the BCM fault case analyses.

Rouzky, A. R.↗