Search NASA⌕ Search

SEARCH · Search NASA

Results for “LINEARIZATION”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Linear Collider Facility (LCF) at CERN

In this paper we outline a proposal for a Linear Collider Facility as the next flagship project for CERN. It offers the opportunity for a timely, cost-effective and staged construction of a new collider that will be able to comprehensively map the Higgs boson's properties, including the Higgs field potential, thanks to a large span in centre-of-mass energies and polarised beams. A comprehensive programme to study the Higgs boson and its closest relatives with high precision requires data at centre-of-mass energies from the Z pole to at least 1 TeV. It should include measurements of the Higgs boson in both major production mechanisms, ee -> ZH and ee -> vvH, precision measurements of gauge boson interactions as well as of the W boson, Higgs boson and top-quark masses, measurement of the top-quark Yukawa coupling through ee ->ttH, measurement of the Higgs boson self-coupling through HH production, and precision measurements of the electroweak couplings of the top quark. In addition, ee collisions offer discovery potential for new particles complementary to HL-LHC.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

The Linear Point Standard Ruler with DESI DR1 and DR2 Data

The linear point, a purely geometric feature in the monopole of the two-point correlation function, has been proposed as an alternative standard ruler. Compared to the peak in the correlation function, it is more robust to late-time nonlinear effects at the percent level. In light of improved simulations and high quality data, we revisit the robustness of the linear point and use it as an alternative to template-based fitting approaches typically used in BAO analyses. We present the linear point measurements on galaxy samples from the first and second data releases (DR1 and DR2) of the DESI survey. We convert the linear point into a dimensionless parameter $α_{iso,LP}$, defined as the ratio of the linear point in the fiducial cosmology and the observed value, analogous to the isotropic BAO scaling parameter $α_{iso}$ used in previous BAO measurements. Using the 2nd generation of AbacusSummit mock catalogs, we find that linear point measurements are more precise when calculated in the post-reconstruction regime with 15-60% smaller uncertainties than those pre-reconstruction. We find a systematic shift in the linear point measurements compared against the isotropic BAO measurements in mocks; we attribute this to the isotropic damping parameter responsible for smearing the linear point in the nonlinear regime. We propose a sample-dependent correction that mitigates the impact of late-time nonlinear effects. While this introduces a cosmology dependence in an otherwise model-independent measurement, this is necessary given the sub-percent precision dictated by current cosmological surveys. Comparing $α_{iso,LP}$ with isotropic BAO measurements made on the DESI DR1 and DR2 galaxy samples, we find excellent agreement after applying this correction, particularly post-reconstruction. We discuss future scope regarding cosmological inference with linear point measurements.

Uberoi, N. [Yale U.] (ORCID:0000000275179629)↗
DOE OSTI2025

Projection-based multifidelity linear regression for data-scarce applications

Surrogate modeling for systems with high-dimensional quantities of interest remains challenging, particularly when training data are costly to acquire. This work develops multifidelity methods for multiple-input multiple-output linear regression targeting data-limited applications with high-dimensional outputs. Multifidelity methods integrate many inexpensive low-fidelity model evaluations with limited, costly high-fidelity evaluations. We introduce two projection-based multifidelity linear regression approaches with linear and nonlinear features that leverage principal component basis vectors for dimensionality reduction and combine multifidelity data through: (i) a direct data augmentation using low-fidelity data, and (ii) a data augmentation incorporating explicit linear corrections between low-fidelity and high-fidelity data. The data augmentation approaches combine high-fidelity and low-fidelity data into a unified training set and train the linear regression model through weighted least squares with fidelity-specific weights. We introduce a proximity-based weighting scheme with automatic weight selection strategy through cross-validation. Here, the proposed multifidelity linear regression methods are demonstrated on approximating the surface pressure field of a hypersonic vehicle in flight and the temperature field on an aircraft disc braking system. In an ultra low-data regime of no more than twelve high-fidelity samples, multifidelity linear regression achieves approximately 2% – 12% improvement in median accuracy and a higher R 2 score relative to single-fidelity methods at comparable computational cost.

data augmentation↗
DOE OSTI2025

A new tubular ironless linear induction motor using non-interrupted winding method

The ironless linear induction motor device can be used for any application which requires a large motion span or rotation restrictedly by the ironless feature due to the strong magnetic field environment. One of the possible applications is ITER diagnostics shutters. The tokamak reactor systems have a severe radiation issue during the operation. The traditional Kapton-insulated wires have significantly lower radiation resistance than an expected large radiation dose at ITER Motion Stark Effect (MSE) shutter location. Furthermore, the traditional single-sided linear induction motor winds the wire to make the N-turn pitch coils independently and then connect them in series to make a one-phase winding. Therefore, it is a big challenge to make so many radiation-protection connectors in such a limited space as inside the vacuum vessel of fusion reactors. The approved mineral-insulated kinds of cables need a large bend radius for the end turns, so the traditional single-sided linear motor with two-layer three-phase windings need to be upgraded to the tubular linear motors. The new proposed tubular ironless linear induction motor can not only overcome the large bend radius issue with the mineral-insulated cables but also allow the motor coils wound continuously around the tubular-shaped stator, forward and backward N times to make a N-turn winding for each phase, in order to resolve the radiation protection issue with the winding connection. The 3D transient Maxwell electromagnetic model with a tubular three-phase linear induction motor has been analyzed. Furthermore, this paper will present the detailed design and winding method of ironless linear induction motor actuators, and results of the 3D transient numerical simulation of the actuator performance.

Ironless↗
DOE OSTI2026

Multibody for Everybody (M4E) - A Linearization Approach to Enable Frequency Domain Analysis, Time Integration and Control Co-Design

1.1 Background/Objectives: Marine energy represents a promising yet underexploited source of power. To increase the harvested power, significant efforts have been made to improve wave energy converter (WEC) modeling capabilities and optimize power take-off (PTO) performance; however, these efforts have often treated WEC dynamics, PTO design, and controller development sequentially. In contrast, control co-design (CCD) is emerging as a promising strategy to address these issues directly, creating a growing need for fast analysis tools suitable for repeated simulation and parametric studies [1]. To support this need, this work presents the Multibody for Everybody (M4E) [2] linearization module, which employs a symbolic toolbox to provide deeper insight of WEC design parameters. The objective is to demonstrate that a minimal-coordinate linearization of articulated WEC dynamics can provide accurate wave response predictions and substantial computational savings relative to nonlinear time-domain simulation, while preserving compatibility with broader wave-energy analysis workflows, enabling CCD. 1.2 Approach/Activities: The proposed approach linearizes the equations of motion, generated by M4E, in minimal coordinates about a selected operating point and combines the resulting system with frequencydomain hydrodynamic terms to incorporate the reduced mass, damping, stiffness, and forcing operators. The linearized model is used for both impedance-based response amplitude operator (RAO) prediction and rapid regular-wave time integration. The methodology is demonstrated on a single-flap device and a FOSWEC configuration, with linearized M4E responses compared against the corresponding nonlinear M4E simulations and WEC-Sim results. Regular-wave time histories, RAO trends, and runtime differences are assessed. The framework is also compatible with broader wave-energy workflows, including coupling to WecOptTool, although that capability is not the focus of this work [3]. 1.3 Results/Lessons: The linearized M4E model reproduces key regularwave response characteristics such as integration and Response Amplitude over multiple frequencies. This module matches nonlinear M4E and WEC-Sim results while substantially reducing integration cost. Thus, the proposed framework can serve as a rapid analysis layer for articulated WEC design, parameter studies, and controls-oriented workflows. The analysis is most appropriate in the near-equilibrium regime, about the linearization point.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Measuring Topological Constraint Relaxation in Ring-Linear Polymer Blends

Polymers are an effective test bed for studying topological constraints in condensed matter due to a wide array of synthetically available chain topologies. When linear and ring polymers are blended together, emergent rheological properties are observed as the blend can be more viscous than either of the individual components. This emergent behavior arises since ring-linear blends can form long-lived topological constraints as the linear polymers thread the ring polymers. Here, in this work, we demonstrate how the Gauss linking integral can be used to efficiently evaluate the relaxation of topological constraints in ring-linear polymer blends. For majority-linear blends, the relaxation rate of topological constraints depends primarily on reptation of the linear polymers, resulting in the diffusive time τ d,R for rings of length N R blended with linear chains of length N l to scale as τ d,R ~ N$^2_R$N$^{3.4}_L$.

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

Enabling Efficient Sparse Computations using Linear Algebra Aware Compilers

This project developed the LAPIS compiler framework, built on the Multilevel Intermediate Representation (MLIR), to optimize sparse linear algebra operations and support performance portability across diverse architectures. The main innovation of LAPIS is the Kokkos dialect, which allows for lowering codes from a high productivity language to different architectures in an elegant way. The dialect also allows the conversion of lower-level MLIR code to C++ Kokkos code, facilitating the integration of scientific machine learning (SciML) models into applications. To extend LAPIS for distributed memory architectures, a new partition dialect was created to manage the distribution of sparse tensors and express communication patterns for sparse linear algebra operations. This dialect also supports the distributed execution of operators and includes algorithmic optimizations to minimize communication to improve performance. The project also demonstrates that MLIR can enable effective linear algebra-level optimizations, improving performance on different GPUs for both sparse and dense linear algebra kernels. Key applications of LAPIS include sparse linear algebra and graph kernels, TenSQL, a relational database management solution built on GraphBLAS, and the development of subgraph isomorphism and monomorphism kernels, showcasing performance portability. In summary, the LAPIS framework supports productivity, performance, portability, and distributed memory execution, while also enabling linear algebra-level optimizations that are challenging in traditional programming languages, with successful applications ranging from simple sparse linear algebra to complex graph kernels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Hyperplane decision trees as piecewise linear surrogate models for chemical process design

Recent trends in chemical engineering research point towards an increasing reliance on data-driven modeling approaches. Neural networks, for instance, have proven to be accurate when data is plentiful and high-dimensional, but in many cases, they require computationally-intensive training procedures. Here, in this work, we describe hyperplane decision trees (HT) as a highly expressive and low-compute machine learning model architecture. These models are locally linear and have linear decision boundaries, resulting in a piecewise linear model of the data. This property allows them to be converted into mixed-integer linear constraints which can be globally optimized. Our open-source PyTorch implementation of this method is a fast, flexible, and accessible way to build accurate piecewise linear models of data.

Decision trees↗
DOE OSTI2025

Randomized Adiabatic Quantum Linear Solver Algorithm with Optimal Complexity Scaling and Detailed Running Costs

Solving linear systems of equations is a fundamental problem with a wide variety of applications across many fields of science, and there is increasing effort to develop quantum linear solver algorithms. Subaşı et al. [Phys. Rev. Lett. 122, 060504 (2019)] proposed a randomized algorithm inspired by adiabatic quantum computing, based on a sequence of random Hamiltonian simulation steps, with suboptimal scaling in the condition number 𝜅 of the linear system and the target error 𝜖. Here we go beyond these results in several ways. Firstly, using filtering [Lin and Tong, Quantum 4, 361 (2020)] and Poissonization techniques [Cunningham and Roland, ArXiv:2406.03972 (2024)], the algorithm complexity is improved to the optimal scaling 𝑂⁡(𝜅⁢log (1/𝜖))—an exponential improvement in 𝜖, and a shaving of a log 𝜅 scaling factor in 𝜅. Secondly, the algorithm is further modified to achieve constant factor improvements, which are vital as we progress towards hardware implementations on fault-tolerant devices. We introduce a cheaper randomized walk operator method replacing Hamiltonian simulation—which also removes the need for potentially challenging classical precomputations; randomized routines are sampled over optimized random variables; circuit constructions are improved. We obtain a closed formula rigorously upper bounding the expected number of times one needs to apply a block-encoding of the linear system matrix to output a quantum state encoding the solution to the linear system. The upper bound is 837⁢𝜅 at 𝜖 = 10 −10 for Hermitian matrices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Role of Polymer Architecture in CO 2 Capture from Air Using Supported Poly(alkylenimine)s: Linear vs Branched Polymers

Direct air capture (DAC) of CO 2 coupled with geologic storage is a promising climate change mitigation strategy, with some applications employing amines supported on porous solids as CO 2 sorbents. While branched poly(ethylenimine) (PEI) is the standard benchmark amine material, it suffers from limited oxidative stability. Poly(propylenimine) (PPI), as an alternative, has previously demonstrated improved resistance to degradation under harsh oxidative conditions. Linear and branched PEI are commercially available, though at different molecular weights, while PPI is not commercially available. For this reason, a comparative study of all four polymers (linear PEI, branched PEI, linear PPI, branched PPI) has not been reported for DAC. In this study, we synthesize and compare low-molecular-weight (∼800 g/mol) linear (L) and branched (B) PEI and PPI supported on a model support, SBA-15 silica. These materials are evaluated for CO 2 adsorption under dry, DAC-relevant conditions (400 ppm of CO 2 , 30 °C). LPPI exhibited the highest amine efficiency at all loadings, reaching a maximum of 0.14 mmol CO 2 /mmol N, outperforming BPEI, while LPEI consistently showed the lowest uptake capacity. Temperature-programmed desorption reveals that the structure of the amine polymer impacts the CO 2 binding strength, with branched polymers displaying higher desorption energies of 102−111 kJ/mol. In situ infrared spectroscopy experiments show that all sorbents preferentially capture CO 2 as ammonium carbamate. Isobaric CO 2 uptake studies further underscore the influence of polymer mobility and support pore crowding on performance, while demonstrating the sorbents’ performance at elevated temperatures and CO 2 concentrations. All materials demonstrated good stability over 25 adsorption−desorption cycles using thermal regeneration in an inert gas purge, with only BPPI displaying a 10−11% decrease in capacity/amine efficiency during cycling, possibly due to the loss of low molecular weight, oligomeric amines. This is the first side-by-side comparison of the CO 2 sorption properties of linear and branched PEI and PPI with similar molecular weights. These findings highlight the significant role of polymer architecture in CO 2 capture efficiency and inform future designs of durable, high-performance DAC sorbents.

adsorbents↗
DOE OSTI2026

Augmenting subspace optimization methods with linear bandits

In this work, we consider the framework of methods for unconstrained minimization that are, in each iteration, restricted to a model that is only a valid approximation to the objective function on some affine subspace containing an incumbent point. These methods are of practical interest in computational settings where derivative information is either expensive or impossible to obtain. Recent attention has been paid in the literature to employing randomized matrix sketching for generating the affine subspaces within this framework. We consider a relatively straightforward, deterministic augmentation of such a generic subspace optimization method. In particular, we consider a sequential optimization framework where actions consist of one-dimensional linear subspaces and rewards consist of (approximations to) the magnitudes of directional derivatives computed in the direction of the action subspace. Reward maximization in this context is consistent with maximizing lower bounds on descent guaranteed by first-order Taylor models. This sequential optimization problem can be analysed through the lens of dynamic regret. We modify an existing linear upper confidence bound (UCB) bandit method and prove sublinear dynamic regret in the subspace optimization setting. We demonstrate the efficacy of employing this linear UCB method in a setting where forward-mode algorithmic differentiation can provide directional derivatives in arbitrary directions and in a derivative-free setting. For the derivative-free setting, we propose SS-POUNDers, an extension of the derivative-free optimization method POUNDers that employs the linear UCB mechanism to identify promising subspaces. Our numerical experiments suggest a preference, in either computational setting, for employing a linear UCB mechanism within a subspace optimization method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Quantum Algorithm for Linear Non-unitary Dynamics with Near-Optimal Dependence on All Parameters

We introduce a family of identities that express general linear non-unitary evolution operators as a linear combination of unitary evolution operators, each solving a Hamiltonian simulation problem. This formulation can exponentially enhance the accuracy of the recently introduced linear combination of Hamiltonian simulation (LCHS) method [An, Liu, and Lin, Physical Review Letters, 2023]. For the first time, this approach enables quantum algorithms to solve linear differential equations with both optimal state preparation cost and near-optimal scaling in matrix queries on all parameters.

Applied Dynamical Systems↗
DOE OSTI2024

Iterative methods in GPU-resident linear solvers for nonlinear constrained optimization

Linear solvers are major computational bottlenecks in a wide range of decision support and optimization computations. The challenges become even more pronounced on heterogeneous hardware, where traditional sparse numerical linear algebra methods are often inefficient. For example, methods for solving ill-conditioned linear systems have relied on conditional branching, which degrades performance on hardware accelerators such as graphical processing units (GPUs). To improve the efficiency of solving ill-conditioned systems, our computational strategy separates computations that are efficient on GPUs from those that need to run on traditional central processing units (CPUs). Our strategy maximizes the reuse of expensive CPU computations. Iterative methods, which thus far have not been broadly used for ill-conditioned linear systems, play an important role in our approach. In particular, we extend ideas from Arioli et al., (2007) to implement iterative refinement using inexact LU factors and flexible generalized minimal residual (FGMRES), with the aim of efficient performance on GPUs. In conclusion, we focus on solutions that are effective within broader application contexts, and discuss how early performance tests could be improved to be more predictive of the performance in a realistic environment.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Energy recovery proton linear accelerator

High-power proton linear accelerators have important applications in both scientific research and industry. However, the operation of such accelerators with megawatt-level beam power is expensive and limits the broad availability of these facilities. In this Letter, we propose a novel energy recovery proton linear accelerator in which the final GeV-level proton beam is reinjected into the linear accelerator from the accelerator exit and is decelerated in the same accelerator down to about 2 MeV, close to the initial beam energy. This substantially reduces the power consumption of the proton linear accelerator and also avoids the need for a high-power beam dump. We demonstrate this concept through self-consistent simulations and show that the beam-beam effect between the forward-accelerating beam and the backward-decelerating beam would not be a limiting factor for the proposed concept. A potential application of this concept to an electron-ion collider based on energy recovery electron and ion accelerators is discussed.

Energy recovery linacs↗
DOE OSTI2024

Preliminary Design of Ironless Linear Induction Motors for ITER MSE Shutter Actuators

In this article, the shutters need to be closed or opened during the operation and calibration period time for the fusion diagnostic systems, such as the International Thermonuclear Experimental Reactor (ITER) motional stark effect (MSE). Standard electric motors and actuators will not work in a strong magnetic field environment due to the presence of the magnetic field of fusion reactors. The innovative linear induction motor (LIM) with an ironless feature overcomes this kind of challenge and can be used for any application that requires controlled motion with a large stroke in the magnetic field environment. It consists of a high-electrical conductivity plate like copper and three-phase motor windings with nonferrous or stainless-steel stators. The conductor plate attached with the shutter is simply driven by the underneath three-phase linear windings to achieve linear motion. The motor drive can be controlled remotely by a controller using the electrical connection, so no sensitive electronic components are located in the harsh environment where the motor itself is located. The design requirements and test facility have been described. Several 3-D transient Maxwell electromagnetic (EM) models with different locations of three-phase linear motor windings and shutter stators have been analyzed and evolved to meet any applications allowed in the harsh environment inside the vacuum vessel of fusion reactors. The preliminary design results are presented in this article.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗