Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear approximation”

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 OSTI2025

A 291-day Evaluation of the Performance of a Consumer-grade Temporal Radon Detector

Affordable, accurate, and robust temporal measurement devices are desirable for screening and assessment of radon levels in private homes and workplaces. This research expands upon prior research, using the RadonFTlab RadonEye device through a comparison of multiple samples of this instrument with a laboratory-grade instrument, the Saphymo AlphaGUARD, over a more extensive period than reported previously. Data were collected over 291 d in a poorly ventilated basement space in an occupied building. Environmental conditions varied naturally, changing both the radon source term and radon entry into the space approximating typically deployed conditions. The R-squared linear regression correlation coefficient and relative sensitivities of each RadonEye with the AlphaGUARD were computed. Altogether temporal and diurnal variations were also studied. The sensitivities of all RadonEyes and the AlphaGUARD agreed to within 22% throughout the entire deployment period.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Stochastic Microgrid Scheduling With Chance‐Constrained Resilience Consideration

Traditionally, it is assumed that microgrids transition seamlessly from grid‐connected operation to islanded mode in the event of sudden main grid outages. In reality, the islanding process, especially unintentional islanding, is rarely seamless. Instead, it is subject to voltage and frequency fluctuations caused by the instantaneous disconnection of the point of common coupling (PCC) switch, variations in loads and renewable generation output and even the protection tripping of distributed energy resources (DERs). To mitigate these fluctuations and facilitate a smooth islanding process, we propose a stochastic microgrid scheduling model that incorporates chance‐constrained resilience measures. Specifically, the resilience measure is defined as the probability of successful islanding (PSI), that is, the probability that a microgrid can mitigate the generation‐demand imbalance caused by the disconnection of the PCC switch, variations in load and renewable generation and DER tripping. This measure is modelled using chance constraints. Unlike existing reliability and resilience indices, which typically neglect the possibility of microgrid/DER failure under extreme events and assume their survival while primarily focussing on reducing impact duration or magnitude, the proposed PSI‐based framework explicitly addresses microgrid and DER survival during the islanding transition. The formulated nonlinear chance constraints are approximated using a multiinterval approach and equivalently represented as a mixed‐integer linear programming (MILP) formulation. Case study results validate the proposed method, showing that the PSI estimation error is reduced to less than 8%, compared to approximately 28% with existing methods. Various sensitivity analyses on the DER tripping rate and PSI settings were performed to validate the robustness of the proposed method. In particular, the necessity of accounting for DER tripping in the PSI calculation was demonstrated.

chance constrained optimization↗
DOE OSTI2024

Optimizing spectral phase transfer in four-wave mixing with gas-filled capillaries

Four-wave mixing (FWM) in gas-filled hollow-core capillaries, a nonlinear optical process that mixes signal and pump photon frequencies to generate idler frequency photons, offers a method for precise spectral phase transfer from signal to idler at ultrashort timescales and extreme powers. However, this regime is challenged by competing linear and nonlinear dynamics, leading to significant trade-offs between spectral phase transfer and conversion efficiency. Our computational investigation focuses on the upconversion of femtosecond pulses from the infrared (IR) to the ultraviolet (UV), a range notoriously difficult to manipulate. We explore an intermediate energy regime that strikes an optimal balance between FWM-mediated phase-transfer fidelity and nonlinear conversion efficiency. By adjusting the energy ratios and spectral phase profiles of the input signal, we achieve conversion efficiencies of approximately 5-15% while maintaining an effective quasi-linear spectral phase transfer. These findings will contribute to establishing first-principles and scaling laws essential for applications such as high-precision imaging, spectroscopy, quantum transduction, and distributed entangled interconnects, facilitating advanced control of ultrafast photonic and electronic wavepackets in quantum materials with unprecedented spatial and temporal precision.

Zhang, Hao↗
DOE OSTI2025

A neural-network-enhanced parameter-varying framework for multi-objective model predictive control applied to buildings

Management of the electrical grid is becoming more complex due to the increased penetration of alternative energy generation technologies and a broadening diversity of electric loads. This complexity creates challenges in balancing demand and generation that can increase the potential for grid instabilities. One effective way to address this issue is to leverage previously unexploited demand flexibility through advanced control strategies. In this work, we propose an advanced control method, called adaptive neural parameter-varying model predictive control (ANPV-MPC), to control the temperature and energy consumption of a building via its Heating, Ventilation, and Air Conditioning system. ANPV-MPC combines key ideas in parameter-varying control, adaptive control, and online learning strategies to bridge the gap between computationally efficient linear model predictive control and more accurate nonlinear model predictive control. The novelty in ANPV-MPC is the use of a physics-inspired Bayesian neural network to estimate the coefficients of the parameter-varying linear control model. The Bayesian neural network additionally provides uncertainty estimates, triggering online training to capture evolving building system conditions. We show that ANPV-MPC can approximate the building system dynamics with a 28.39% higher accuracy than traditional linear model predictive control, resulting in 36.23% better control performance without increasing complexity of the optimal control problem. ANPV-MPC also adapts in real time to previously unseen conditions using online learning, further improving its performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Heterogeneous Mixtures of Dictionary Functions to Approximate Subspace Invariance in Koopman Operators: Why Deep Koopman Operators Work

Abstract Koopman operators model nonlinear dynamics as a linear dynamic system acting on a nonlinear function as the state. This nonstandard state is often called a Koopman observable and is usually approximated numerically by a superposition of functions drawn from a dictionary . In a widely used algorithm, extended dynamic mode decomposition (EDMD), the dictionary functions are drawn from a fixed class of functions. Deep learning combined with EDMD has been used to learn novel dictionary functions in an algorithm called deep dynamic mode decomposition (deepDMD). The learned representation both (1) accurately models and (2) scales well with the dimension of the original nonlinear system. In this paper, we analyze the learned dictionaries from deepDMD and explore the theoretical basis for their strong performance. We explore State-Inclusive Logistic Lifting (SILL) dictionary functions to approximate Koopman observables. Error analysis of these dictionary functions show they satisfy a property of subspace approximation, which we define as uniform finite approximate closure. Typically, a Koopman dictionary’s nonlinear functions are homogeneous. In this paper, we discover that structured mixing of heterogeneous dictionary functions drawn from different classes of nonlinear functions achieve the same accuracy and dimensional scaling as the deep-learning-based deepDMD algorithm Yeung et al. ( In: 2019 American Control Conference (ACC), 2019). We specifically show this by building a heterogeneous dictionary comprised of SILL functions and conjunctive radial basis functions (RBFs). This mixed dictionary achieves similar accuracy and dimensional scaling to deepDMD with an order of magnitude reduction in parameters, while maintaining geometric interpretability. These results strengthen the viability of dictionary-based Koopman models to solving high-dimensional nonlinear learning problems.

Johnson, Charles A.↗
DOE OSTI2025

Temperature Evolution of the Activation Barriers Leads to Meyer−Neldel Rules for Structural Relaxation and Transport in Polymers

Understanding activation barriers controlling structural relaxation in glass-forming liquids, molecular transport, and ionic conductivity in amorphous polymers is a grand challenge of fundamental scientific and materials engineering interest across disciplines. Over decades, intriguing but puzzling empirical correlations between the elementary time scale of activated barrier crossing and the apparent Arrhenius activation energy, the so-called Meyer−Neldel (MN) rules, have been discovered in diverse liquids and glasses. Here, in this study, we formulate and successfully apply a new experimental analysis and an explicitly dynamical theoretical framework which provides an understanding of the origin, validity, and failure of such correlations, that bridge and unify the three fields of structural relaxation, molecular transport, and ionic conductivity in liquids and quenched glasses. Distinct quasi-universal laws are predicted in equilibrated liquids and nonequilibrium glasses, consistent with experiments. Our analysis reveals that even if the relaxation appears Arrhenius over a limited temperature range, the physical activation barrier is generally temperature-dependent in polymeric systems even below glass transition temperature. In addition, we show that the approximate validity of classical MN rules hinges on a linear temperature dependence of this barrier and the temperature range probed in experiments. Our findings are relevant for controlling the activation barrier in functional soft polymeric materials relevant to molecular separations, barrier coatings, and charge transport, and also provide new constraints on the theoretical understanding of the mechanism underlying slow activated dynamics in glass-forming condensed matter.

Meyer-Neldel rules↗
DOE OSTI2025

Disoriented isospin condensates in heavy-ion collisions

Anomalous neutral to charged kaon correlations measured by the ALICE collaboration have defied usual explanations. We propose that the large fluctuations could arise because of a disoriented isospin condensate where there is an imbalance between up and down condensates at the time kaons hadronize. This could happen in heavy-ion collisions when the quark condensate re-forms as the system cools and the approximate chiral symmetry of QCD is broken. Within the linear sigma model, we show that the energy cost of forming a disoriented isospin condensate is small making it very plausible.

Kapusta, Joseph [University of Minnesota, Minneapo↗
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

Parameterized anomalous transport model for current-carrying collisionless plasmas in pulsed power inertial confinement fusion

Current delivery in pulsed power inertial confinement fusion is influenced by collisionless current-carrying microturbulent plasmas, which are sourced from electrode surfaces. In this setting, the lower hybrid drift instability—triggered by plasma acceleration—is a leading candidate driver of difficult-to-predict momentum and energy transport. To characterize the nonlinear state of the microturbulent plasma, a parameterized anomalous transport model is developed for the instability, with analytic formulas for anomalous collision frequency, resistivity, and species heating rates. The formulas are expressed in terms of linear-theory variables and four dimensionless parameters that characterize the macroscopic plasma state. The model is built on linear theory analysis, power law analysis, and quasilinear theory analysis, and is validated using a series of nonlinear continuum kinetic Vlasov–Poisson simulations. The theoretical and computational investigation demonstrates that the anomalous collision frequency associated with the instability can be reliably approximated, within about a factor of five or better, by the unscaled linear theory growth rate of the fastest-growing wavenumber mode. This finding enables efficient calculation of anomalous resistivity and species heating rates over a wide range of plasma conditions, resulting in improved predictive capabilities.

Complex functions↗
DOE OSTI2025

Thermodynamics and collisionality in firehose-susceptible high- β plasmas

We study the evolution of collisionless plasmas that, due to their macroscopic evolution, are susceptible to the firehose instability, using both analytic theory and hybrid-kinetic particle-in-cell simulations. We establish that, depending on the relative magnitude of the plasma β, the characteristic time scale of macroscopic evolution and the ion-Larmor frequency, the saturation of the firehose instability in high-β plasmas can result in three qualitatively distinct thermodynamic (and electromagnetic) states. By contrast with the previously identified ‘ultra-high-beta’ and ‘Alfvén-inhibiting’ states, the newly identified ‘Alfvén-enabling’ state, which is realised when the macroscopic evolution time τ exceeds the ion-Larmor frequency by a β-dependent critical parameter, can support linear Alfvén waves and Alfvénic turbulence because the magnetic tension associated with the plasma’s macroscopic magnetic field is never completely negated by anisotropic pressure forces. We characterise these states in detail, including their saturated magnetic-energy spectra. The effective collision operator associated with the firehose fluctuations is also described; we find it to be well approximated in the Alfvén-enabling state by a simple quasi-linear pitch-angle scattering operator. The box-averaged collision frequency is ν eff ∼ β/τ, in agreement with previous results, but certain subpopulations of particles scatter at a much larger (or smaller) rate depending on their velocity in the direction parallel to the magnetic field. Our findings are essential for understanding low-collisionality astrophysical plasmas including the solar wind, the intracluster medium of galaxy clusters and black hole accretion flows. We show that all three of these plasmas are in the Alfvén-enabling regime of firehose saturation and discuss the implications of this result.

astrophysical plasmas↗
DOE OSTI2025

Microstructurally Strained Pyrochlore–Perovskite Biphasic Electrocatalysts for the Oxygen Evolution Reaction

Efficiency of water splitting for hydrogen production is often limited by the sluggish kinetics of multiple electronic transfers required in the heterogeneous oxygen evolution reaction (OER). Catalyst design for reducing the high OER overpotential remains a major scientific challenge. Lattice-strain engineering, a method for tuning the electronic structure and surface geometric configuration of active sites, may greatly affect the interaction between adsorbates and catalytic surfaces for high activity and stability. Here, in this study, we present the synthesis of biphasic oxides of YPrSrRuMnO x , which consists of distinct phases of Y 2 Ru 2 O 7 pyrochlore and (Pr 0.7 Sr 0.3 )MnO 3 perovskite, and the development of a suitable analytical approach to study the strain–catalytic property relationship. Linear sweep voltammetry results reveal that the biphasic oxide exhibits approximately 3.1 times greater mass activity and 2.4 times larger turnover frequency (TOF) than single-phase Y 2 Ru 2 O 7 in the 0.1 M HClO 4 electrolyte. The biphasic catalyst is also about 3 times more stable than the single-phase oxide under acidic conditions. X-ray photoelectron spectroscopy, nitrogen isotherm, and electrochemical surface area analyses indicate that the oxidation state, specific surface area, and electrochemical surface area do not cause enough difference in the observed enhancement of OER performance. We examined the effects of microstrain on electrocatalysis, originating from lattice mismatch between different phases, using three different structural models. Specifically, we compared the Williamson–Hall method, standard stress–strain analysis, and Rietveld refinement in analyzing the structure–property relationship. Strain mapping using geometric phase analysis (GPA) further revealed significant microstrain and lattice dislocations localized near phase boundaries in the biphasic oxide, in contrast to the uniform strain in single-phase materials. The results reveal that the increased microstrain correlates well with the improved OER performance, as the biphasic oxide catalyst exhibits 2–3 times greater microstrain than Y 2 Ru 2 O 7 pyrochlore.

electrocatalysts↗
DOE OSTI2026

Expansion-driven Weibel instability in magnetohydrodynamics: Linear theory for static systems

The Weibel dispersion relation is obtained in the magnetohydrodynamic (MHD) approximation by including a tensor expression for temperature. MHD gives an upper cutoff wavenumber for Weibel growth identical to kinetic theory if the electron drift velocity is included in the evolution of the electron temperature. MHD overestimates growth rates compared to kinetic theory and gives maximum growth at larger wavenumbers, but it only leads to divergent results as temperature anisotropy tends to infinity. Thermal conduction in MHD lowers the growth rates and shifts maximum growth to smaller wavenumbers, but flux-limited or nonlocal thermal conduction is found to significantly limit these effects. The results are used as a first step in evaluating the potential of MHD to simulate expansion-driven Weibel instability, which occurs due to the temperature tensor only cooling in the direction of expansion. Finally, the small temperature anisotropy found in expansion-driven Weibel instability means that MHD could be an adequate model because growth rates are much less than the electron plasma frequency, and wavelengths are much greater than the Debye length.

Davies, J. R. [Univ. of Rochester, NY (United Stat↗
DOE OSTI2024

A Performance and Energy Study of GPU-Resident Preconditioners for Conjugate Gradient Solvers: In the Context of Existing and Novel Approaches

Optimizing a particular subprogram out of the set of Basic (sparse) Linear Algebra Subprograms (BLAS) for a given architecture is a common topic of research. In applications, however, these BLAS functions rarely appear in isolation; usually, many of them are used together, in various combinations and with varying inputs. As the need to solve a large, sparse linear system is ubiquitous throughout HPC applications, linear solvers constitute a realistic, sufficiently complex and well-defined representative use case for composite BLAS routines. To this end, based on a representative set of matrices drawn from a diverse set of fields, we present a framework to study, from the performance and energy perspective, the efficacy of GPU- resident parallel Conjugate Gradient (CG) linear solver with different preconditioner options, including Gauss-Seidel, Jacobi, and incomplete Cholesky. We also propose a novel GPU-based preconditioner, in which the triangular solves are approximated by an iterative process. The development of this preconditioner was motivated by solving large graph Laplacian linear systems, for which the existing preconditioners either perform slow on GPU-based platforms or are not applicable. We compare the performance of these preconditioners on different hardware accelerator architectures, i.e., AMD MI250X, MI100, Nvidia A100, V100, and Jetson. Our experiments reveal performance trade-offs and provide information on how to select the best strategy for the given linear system, dictated by its properties, and the platform of interest. We demonstrate the application of our novel preconditioner for solving CG and graph Laplacian systems. Overall, the framework can be utilized as a benchmark to guide informed decisions in choosing a specific preconditioner, i.e., whether it is better to rely on the performance of a triangular solver or on the performance of sparse matrix-vector product. Finally, by considering power consumption to solve the linear systems, we report the energy footprint for the solvers.

Preconditioned Conjugate Gradient, GPUs, iterative↗
DOE OSTI2025

Quantitative Modeling of High-Energy Electron Scattering in Thick Samples Using Monte Carlo Techniques

Cryo-electron microscopy (cryo-EM) is a powerful tool for imaging biological samples but is typically limited by sample thickness, which is restricted to a few hundred nanometers depending on the electron energy. However, there is a growing need for imaging techniques capable of studying biological samples up to 10 µm in thickness while maintaining nanoscale resolution. This need motivates the use of mega-electron-volt scanning transmission electron microscopy (MeV-STEM), which leverages the high penetration power of MeV electrons to generate high-resolution images of thicker samples. In this study, we employ Monte Carlo simulations to model electron–sample interactions and explore the signal decay of imaging electrons through thick specimens. By incorporating material properties, interaction cross-sections for energy loss, and experimental parameters, we investigate the relationship between the incident and transmitted beam intensities. Key factors such as detector collection angle, convergence semi-angle, and the material properties of samples were analyzed. Our results demonstrate that the relationship between incident and transmitted beam intensities follows the Beer–Lambert law over thicknesses ranging from a few microns to several tens of microns, depending on material composition, electron energy, and collection angles. The linear depth of silicon dioxide reaches 3.9 µm at 3 MeV, about 6 times higher than that at 300 keV. Meanwhile, the linear depth of amorphous ice reaches 17.9 µm at 3 MeV, approximately 11.5 times higher than that at 300 keV. These findings are crucial for advancing the study of thick biological and semiconductor samples using MeV-STEM.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reduced basis approximations of parameterized dynamical partial differential equations via neural networks

Projection-based reduced order models are effective at approximating parameter-dependent differential equations that are parametrically separable. When parametric separability is not satisfied, which occurs in both linear and nonlinear problems, projection-based methods fail to adequately reduce the computational complexity. Devising alternative reduced order models is crucial for obtaining efficient and accurate approximations to expensive high-fidelity models. In this work, we develop a timestepping procedure for dynamical parameter-dependent problems, in which a neural-network is trained to propagate the coefficients of a reduced basis expansion. This results in an online stage with a computational cost independent of the size of the underlying problem. Here, we demonstrate our method on several parabolic partial differential equations, including a problem that is not parametrically separable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Upcycling Linear Low-Density Polyethylene Waste into Graphene for High Mass Loading Supercapacitors

Polyethylene is notoriously difficult to upcycle because it decomposes into light gases at approximately 350-400 °C which prevents processing it at higher temperatures to convert it into electrode materials. We address this challenge by using an air-based, thermos-oxidative process, which heats linear low-density polyethylene (LLDPE) just below the decomposition point to initiate oxidation and cross-linking of LLDPE alkyl chains. These molecular transformations allow LLDPE to be graphenized at 950 °C without decomposing. The LLDPE-derived graphene (LLDPE-G) has a BET specific surface area up to 1,800 m2/g and Raman ID/IG ratio of 0.85. When used as electrode material for a symmetric supercapacitor with the 1 M H2SO4 electrolyte, LLDPE-G possesses an specific capacitance up to 175 F/g at mass loading of 20 mg/cm2, yielding an areal capacitance of 3.5 F/cm2. The cycling stability with capacitance demonstrates a retention of 95.8% after 100,000 cycles at high current density of 4.0 A/g.

Gao, Yuan↗
DOE OSTI2024

Constrained Local Approximate Ideal Restriction for Advection-Diffusion Problems

Herein this paper focuses on developing a reduction-based algebraic multigrid (AMG) method that is suitable for solving general (non)symmetric linear systems and is naturally robust from pure advection to pure diffusion. Initial motivation comes from a new reduction-based AMG approach, $\ell \text{AIR}$ (local approximate ideal restriction), that was developed for solving advection-dominated problems. Though this new solver is very effective in the advection-dominated regime, its performance degrades in cases where diffusion becomes dominant. This is consistent with the fact that in general, reduction-based AMG methods tend to suffer from growth in complexity and/or convergence rates as the problem size is increased, especially for diffusion-dominated problems in two or three dimensions. Motivated by the success of $\ell \text{AIR}$ in the advective regime, our aim in this paper is to generalize the AIR framework with the goal of improving the performance of the solver in diffusion-dominated regimes. To do so, we propose a novel way to combine mode constraints as used commonly in energy-minimization AMG methods with the local approximation of ideal operators used in $\ell \text{AIR}$. The resulting constrained $\ell \text{AIR}$ algorithm is able to achieve fast scalable convergence on advective and diffusive problems. In addition, it is able to achieve standard low complexity hierarchies in the diffusive regime through aggressive coarsening, something that was previously difficult for reduction-based methods.

97 MATHEMATICS AND COMPUTING↗