SEARCH · Engineering Papers

Results for “linear systems”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Market optimization and technoeconomic analysis of hydrogen-electricity coproduction systems

Decarbonization efforts across North America, Europe, and beyond rely on variable renewable energy sources such as wind and solar, as well as alternative fuels, such as hydrogen, to support the sustainable energy transition. These advancements have prompted a need for more flexibility in the electric grid to complement non-dispatchable energy sources and increased demand from electrification. Integrated energy systems are well suited to provide this flexibility, but conventional technoeconomic modeling paradigms neglect the time-varying dynamic nature of the grid and thus undervalue resource flexibility. In this work, we develop a computational optimization framework for dynamic market-based technoeconomic comparison of integrated energy systems that coproduce low-carbon electricity and hydrogen (e.g., solid oxide fuel cells, solid oxide electrolysis) against technologies that only produce electricity (e.g., natural gas combined cycle with carbon capture) or only produce hydrogen. Our framework starts with rigorous physics-based process models, built in the open-source Institute for the Design of Advanced Energy Systems (IDAES) modeling and optimization platform, for six energy process concepts. Using these rigorous models and a workflow to optimally design each technology, the framework is shown to be capable of evaluating new and emerging technologies in varying energy markets under a plethora of future scenarios (i.e., renewables penetration, carbon tax, etc.). Ultimately, our framework finds that solid oxide fuel cell-based coproduction systems achieve positive profits for 85% of the analyzed market scenarios. From these market optimization results, we use multivariate linear regression (R 2 values up to 0.99) to determine which electricity price statistics are most significant to predict the optimized annual profit of each system. The proposed framework provides a powerful tool for directly comparing flexible, multi-product energy process concepts to help discern optimal technology and integration options.

08 HYDROGEN↗

DOE OSTI2024

A Hierarchical OPF Algorithm with Improved Gradient Evaluation in Three-Phase Networks

Linear approximation commonly used in solving alternating-current optimal power flow (AC-OPF) simplifies the system models but incurs accumulated voltage errors in large power networks. Such errors will make the primal-dual type gradient algorithms converge to solutions with voltage violation. In this paper, we improve a recent hierarchical OPF algorithm that rested on primal-dual gradients evaluated with a linearized distribution power flow model. Specifically, we propose a more accurate gradient evaluation method based on an unbalanced three-phase nonlinear distribution power flow model to mitigate the errors arising from linearization. The resultant gradients feature a blocked structure that enables our development of an improved hierarchical primal-dual algorithm to solve the OPF problem. Numerical results on the IEEE 123-bus test feeder and a 4,518-node test feeder show that the proposed method can enhance voltage safety at comparable computational efficiency with the linearized algorithm.

approximation algorithms↗

DOE OSTI2024

Stochastic Maxwell-Bloch equations for modeling amplified spontaneous emission

An approach for studying atom-radiation interaction has been developed, associating quantum operators with stochastic variables governed by discrete Heisenberg equations. This framework models general multilevel atomic systems using the paraxial approximation in both single- and multipass configurations. Simulation results for x-ray laser oscillators with parameters of interest are presented. Finally, the stochastic model is efficient to solve, as the required computational resources scale linearly with the number of emitters, and it may be applicable to various other quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Equation-Free Coarse Control of Distributed Parameter Systems via Local Neural Operators

The control of high-dimensional distributed parameter systems (DPS) remains a challenge when explicit coarse-grained equations are unavailable. Classical equation-free (EF) approaches rely on fine-scale simulators treated as black-box timesteppers. However, repeated simulations for steady-state computation, linearization, and control design are often computationally prohibitive, or the microscopic timestepper may not even be available, leaving us with data as the only resource. We propose a data-driven alternative that uses local neural operators, trained on spatiotemporal microscopic/mesoscopic data, to obtain efficient short-time solution operators. These surrogates are employed within Krylov subspace methods to compute coarse steady and unsteady-states, while also providing Jacobian information in a matrix-free manner. Krylov-Arnoldi iterations then approximate the dominant eigenspectrum, yielding reduced models that capture the open-loop slow dynamics without explicit Jacobian assembly. Both discrete-time Linear Quadratic Regulator (dLQR) and pole-placement (PP) controllers are based on this reduced system and lifted back to the full nonlinear dynamics, thereby closing the feedback loop.

93B52, 93C20, 47N70, 65J15, 65M32, 68T07, 68T20, 6↗

DOE OSTI2024

Parametric Instability of Alfvén Waves and Wave Packets in Periodic and Open Systems

The parametric decay instability of Alfvén waves has been widely studied, but few investigations have examined wave packets of finite size and the effect of different boundary conditions on the growth rate. In this paper, we perform a linear analysis of circular and arc-polarized wave trains and wave packets in periodic and open boundary systems in a low- β plasma. We find that both types of wave are 3–5 times more stable in open boundary conditions compared to periodic. Additionally, once the wave packet width ℓ becomes smaller than the system size L , the growth rate decreases nearly with a power law γ ∝ ℓ / L . This study demonstrates that the stability of a pump wave cannot be separated from the laboratory settings, and that the growth rate of daughter waves depends on the conditions downstream and upstream of the pump wave and on the fraction of volume it fills. Our results can explain simulations and experiments of localized Alfvén waves. They also suggest that Alfvénic fluctuations in the solar wind, including sharp impulses known as switchbacks, can be more stable than traditional theory suggests depending on wind conditions.

Alfven waves↗

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 OSTI2021

A strong nonequilibrium bound for sorting of cross-linkers on growing biopolymers

Significance Understanding the connection between microscopic driving forces in biological systems and their emergent structures remains an outstanding challenge. We tackle this problem by deriving a fluctuation–response bound in the context of actin polymerization and bundling and verify its robustness far from equilibrium through simulation. Our results show how actin polymerization, molecular fluxes of cross-linker binding, and the morphology of the growing bundle are associated quantitatively; when the bound is saturated, it becomes a nonequilibrium linear-response relation. Molecular driving forces are difficult to measure directly in biological systems, and this work represents a first step toward estimating them from experimentally accessible quantities such as molecular fluxes and compositions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Using Data-Driven Prediction of Downstream 1D River Flow to Overcome the Challenges of Hydrologic River Modeling

Methods for downstream river flow prediction can be categorized into physics-based and empirical approaches. Although based on well-studied physical relationships, physics-based models rely on numerous hydrologic variables characteristic of the specific river system that can be costly to acquire. Moreover, simulation is often computationally intensive. Conversely, empirical models require less information about the system being modeled and can capture a system’s interactions based on a smaller set of observed data. This article introduces two empirical methods to predict downstream hydraulic variables based on observed stream data: a linear programming (LP) model, and a convolutional neural network (CNN). We apply both empirical models within the Colorado River system to a site located on the Green River, downstream of the Yampa River confluence and Flaming Gorge Dam, and compare it to the physics-based model Streamflow Synthesis and Reservoir Regulation (SSARR) currently used by federal agencies. Results show that both proposed models significantly outperform the SSARR model. Moreover, the CNN model outperforms the LP model for hourly predictions whereas both perform similarly for daily predictions. Although less accurate than the CNN model at finer temporal resolution, the LP model is ideal for linear water scheduling tools.

13 HYDRO ENERGY↗

DOE OSTI2024

Exploring nonlinear Rashba effect and spin Hall conductivity in Janus MXenes W 2 ⁢CO ⁢𝑋 (𝑋=S, Se, Te)

Rashba spin-orbit coupling (RSOC) facilitates spin manipulation without relying on an external magnetic field, opening up exciting possibilities for advanced spintronic devices. In this paper, we examine the effects of crystal momentum (𝑘) nonlinearity and anisotropy on the conventional Rashba effect, with a particular focus on their impact on the spin Hall conductivity (SHC) in a newly predicted family of 2D Janus materials, W 2 ⁢CO⁢𝑋 (𝑋 =S, Se, Te). Using first-principles density functional theory calculations, we confirm the dynamical and mechanical stability of the studied 2D materials. Strikingly, this materials family exhibits pronounced nonlinear Rashba spin splitting at the Γ point of Brillouin zone near the Fermi level, which cannot be adequately described by the linear-𝑘 Rashba model. Therefore, third-order momentum contributions (𝑘 3 ) must be incorporated into the Rashba Hamiltonian. Our analysis reveals that among the studied systems, W 2 ⁢COS exhibits the highest 𝑘 3 contribution of −45.9 eV Å 3 , despite having the lowest linear Rashba constant. Here, a detailed analysis of electronic structure reveals topologically nontrivial behavior in these 2D materials, yielding sizable SHC that is primarily governed by the nonlinear Rashba effect. Notably, these materials also exhibit large spin Hall angle (0.018–2.5 at E 𝐹 ), which is comparable to that of in bulk topological insulators like Bi 2 ⁢Se 3 and Bi 2 ⁢Te 3 , and surpassing those in narrow bandgap bulk semiconductors GeTe and SnTe, as well as heavy metals such as Pt. Sizable SHC, large spin Hall angles, and the ability to tune SHC via electric fields without altering the topological properties, rooted in the crystal field splitting, underscore the potential of these materials for spintronic applications.

Electronic structure↗

DOE OSTI2023

Unified, Geometric Framework for Nonequilibrium Protocol Optimization

Controlling thermodynamic cycles to minimize the dissipated heat is a long-standing goal in thermodynamics, and more recently, a central challenge in stochastic thermodynamics for nanoscale systems. Here, we introduce a theoretical and computational framework for optimizing nonequilibrium control protocols that can transform a system between two distributions in a minimally dissipative fashion. These protocols optimally transport a system along paths through the space of probability distributions that minimize the dissipative cost of a transformation. Furthermore, we show that the thermodynamic metric—determined via a linear response approach—can be directly derived from the same objective function that is optimized in the optimal transport problem, thus providing a unified perspective on thermodynamic geometries. As a result, we investigate this unified geometric framework in two model systems and observe that our procedure for optimizing control protocols is robust beyond linear response.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Synergistic Alignment of Low Aspect‐Ratio π‐Conjugated Molecules Enables Exceptional UV–vis–NIR Polarization Detection

Abstract Polarization detection enhances signal contrast and is widely utilized in diverse advanced applications. An ongoing challenge is the development of high‐performance polarization‐sensitive photodetectors based on optically anisotropic organic semiconductors, particularly in the near‐infrared (NIR) region. While uniaxially aligned π‐conjugated polymers with high aspect ratios exhibit strong linear dichroism and have shown promise, their limited NIR performance and heavy reliance on polymer material now represent critical limitations. Here, a breakthrough is reported in achieving giant linear dichroism and exceptional polarization detection with low aspect‐ratios (AR) non‐fullerene small‐molecule (NFSM) acceptors, extending polarization sensitivity from the UV–vis to the NIR range. An impressive dichroic ratio of 27.1 at 605 nm and 12.0 at 780 nm is demonstrated. The maximum polarization photocurrent ratio is 11.2 at 780 nm under parallel versus perpendicular polarized light. This unprecedented performance originates from synergistic molecular alignment, wherein NFSMs significantly enhance the uniaxial orientation of both the polymer matrix and the NFSMs themselves during self‐assembly and thermal annealing. Besides, such a linear‐polarization‐sensitive photodetectors (LPS‐PDs) are showcased in generating degree‐of‐linear‐polarization imaging. The work establishes NFSMs as a viable material system for next‐generation of organic LPS‐PDs and provides fundamental insights into structural origins of polarization sensitivity in low AR organic semiconductors.

Xue, Yingying↗

DOE OSTI2025

Novel artificial neural network model for instantaneous power losses and operational efficiency mapping of MW-scale vanadium redox flow battery for improved technoeconomic analysis

A novel data-driven, machine-learning-based method for modeling the instantaneous power losses of a distribution-sited 2 MW/8MWh vanadium redox flow battery (VRFB), a grid-scale electrochemical storage technology, is introduced and compared against benchmark empirical modeling approaches, including symmetric and asymmetric models, as well as a recent convex hull modeling approach. The novel loss modeling method introduces several advantages over the benchmark models and over simplistic efficiency estimates, the most significant of which is that the model can accurately reflect the stepwise and non-linear parasitic losses associated with the duty cycles of mechanical auxiliary systems like pump motor drives and blower fans. Residuals of the models are compared; the proposed data driven model features significantly improved accuracy over the benchmark models. The model's coefficient of determination is also improved relative to that of the benchmark models. Furthermore, a novel method for visualization of operational efficiency of the grid-scale storage technology is introduced. To demonstrate the benefits of the novel data-driven method for modeling the VRFB, the benchmark models and the proposed models are embedded into an Open DSS distribution network model to study two applications of the grid-scale electrical storage system: load leveling for grid support and energy arbitrage. This article demonstrates that the accuracy of the instantaneous power loss model significantly impacts the understanding of the state of charge of the VRFB. In turn, the accuracy of the efficiency modeling of the VRFB impacts the understanding of the potential economic value and technical benefits to the distribution network operators. In conclusion, the presented power loss modeling approach is, therefore, highly relevant for utility-stakeholders, battery asset owners, system engineers, system designers, and financial planners interested in evaluating or optimizing the operation of grid-scale VRFBs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗

DOE OSTI2020

Photovoltaic Cleaning Frequency Optimization Under Different Degradation Rate Patterns

Dust accumulation significantly affects the performance of photovoltaic modules and its impact can be mitigated by various cleaning methods. Optimizing the cleaning frequency is essential to minimize the soiling losses and, at the same time, the costs. However, the effectiveness of cleaning lowers with time because of the reduced energy yield due to degradation. Additionally, economic factors such as the escalation in electricity price and inflation can compound or counterbalance the effect of degradation on the soiling mitigation profits. The present study analyzes the impact of degradation, escalation in electricity price and inflation on the revenues and costs of cleanings and proposes a methodology to maximize the profits of soiling mitigation of any system. Here, the energy performance and soiling losses of a 1 MW system installed in southern Spain were analyzed and integrated with theoretical linear and nonlinear degradation rate patterns. The Levelized Cost of Energy and Net Present Value were used as criteria to identify the optimum cleaning strategies. The results showed that the two metrics convey distinct cleaning recommendations, as they are influenced by different factors. For the given site, despite the degradation effects, the optimum cleaning frequency is found to increase with time of operation.

14 SOLAR ENERGY↗

DOE OSTI2020

Observed Warm-Season Characteristics of MCS and Non-MCS Rainfall and Their Recent Changes in the Central United States

Warm-season rainfall characteristics in the central US are investigated as they play important roles in ecohydrology and agricultural productivity. Using rainfall observations, we compare the April-August mesoscale convective systems (MCS) and non-MCS rainfall characteristics and examine their linear trends between 1997 and 2018. We found that MCS rainfall is ~7 times more intense than non-MCS rainfall but it occurs less frequently in time and space. MCS rainfall peaks in nocturnal hours, with synchronized timing of rainfall intensity, area and occurrence, while non-MCS rainfall peaks in late-afternoon hours, mostly attributed to the timing of peak rainfall area. MCS rainfall has increased in the last 22 years due to an increase in frequency and a longer duration per MCS. In contrast, non-MCS rainfall has decreased mainly due to a reduction in rainfall area, leading to fewer total wet days and increased dry intervals between storms

Hu, Huancui↗

DOE OSTI2020

Koopman–von Neumann approach to quantum simulation of nonlinear classical dynamics

Quantum computers can be used to simulate nonlinear non-Hamiltonian classical dynamics on phase space by using the generalized Koopman–von Neumann formulation of classical mechanics. The Koopman–von Neumann formulation implies that the conservation of the probability distribution function on phase space, as expressed by the Liouville equation, can be recast as an equivalent Schrödinger equation on Hilbert space with a Hermitian Hamiltonian operator and a unitary propagator. This Schrödinger equation is linear in the momenta because it derives from a constrained Hamiltonian system with twice the classical phase-space dimension. A quantum computer with finite resources can be used to simulate a finite-dimensional approximation of this unitary evolution operator. Quantum simulation of classical dynamics is exponentially more efficient than a deterministic Eulerian discretization of the Liouville equation if the Koopman–von Neumann Hamiltonian is sparse. Utilizing quantum walk techniques for state preparation and amplitude estimation for the calculation of observables leads to a quadratic improvement over classical probabilistic Monte Carlo algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Nonlinear effects in manybody van der Waals interactions

Van der Waals interactions are ubiquitous and they play an important role for the stability of materials. Current understanding of this type of coupling is based on linear response theory, while optical nonlinearities are rarely considered in this context. Many materials, however, exhibit strong optical nonlinear response, which prompts further evaluation of dispersive forces beyond linear response. Here we present a approach that takes into account linear and nonlinear properties of all dipolar nanoparticles in a given system. This method is based on a Hamiltonian for nonlinear dipoles, which we apply in different systems uncovering a complex interplay of distance, anisotropy, polarizabilities, and hyperpolarizabilities in the vdW energy. This investigation broadens our basic understanding of dispersive interactions, especially in the context of nonlinear materials. Published by the American Physical Society 2024

Le, Dai-Nam (ORCID:0000000307568742)↗

DOE OSTI2025

Probing Curved Spacetime with a Distributed Atomic Processor Clock

Quantum dynamics on curved spacetime has never been directly probed beyond the Newtonian limit. Although we can describe such dynamics theoretically, experiments would provide empirical evidence that quantum theory holds even in this extreme limit. The practical challenge is the minute spacetime curvature difference over the length scale of the typical extent of quantum effects. Here, we propose a quantum network of alkaline earth (like) atomic processors for constructing a distributed quantum state that is sensitive to the differential proper time between its constituent atomic processor nodes, implementing a quantum observable that is affected by post-Newtonian curved spacetime. Conceptually, we propose to delocalize one clock between three locations by encoding the presence or absence of a clock into the state of the local atoms. By separating three atomic nodes over approximately kilometer-scale elevation differences and distributing one clock between them via a 𝑊 state, we demonstrate that the curvature of spacetime is manifest in the interference of the three different proper times that give rise to three distinct beat notes in our nonlocal observable. We further demonstrate that 𝑁-atom entanglement within each node enhances the interrogation bandwidth by a factor of 𝑁. We discuss how our proposed system can probe new facets of fundamental physics, such as the linearity, unitarity, and probabilistic nature of quantum theory on curved spacetime. Our protocol combines several recent advances with neutral atom and trapped ions to realize a novel quantum probe of gravity uniquely enabled by quantum networks.

Covey, Jacob P. [Univ. of Illinois at Urbana-Champ↗