Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Robust stability”

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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Controls of interaction dynamics of orbital assembly

Building structures and spacecraft in orbit will require technologies for positioning, docking/berthing, and joining orbital structures. A fundamental problem underlying the operation of docking and berthing is that of controlling the contact dynamics of mechanical structures actuated by active mechanisms such as robotic devices. Control systems must be designed to control these active mechanisms so that both the free space motions and contact motions are stable and satisfy specifications on position accuracy and bounds on contact forces. For the large orbital structures of the future, the problem of interactive dynamics and control is fundamentally different in several ways than it was for spacecraft docking in the past. First, future space structures must be treated as flexible structures - the operations of docking, berthing, and assembly will need to respect the vibrations of the structures. Second, the assembly of these structures will require multiple-point contact, rather than the essentially single-point positioning of conventional spacecraft docking. Third, some assembly operations require the subassemblies to be brought and held in contact so that successful joining can be accomplished. A preliminary study of contact stability and compliance control design has resulted in the development of an analytical method and a design method to analyze stability. The analytical method analyzes the problem of stability when an actively-controlled structure contacts a passive structure. This method makes it possible to accurately estimate the stiffness of the passive structures with which the contact motion will become unstable. The analytic results suggest that passivity is neither achievable in practice, nor necessary as a design concept. A contact control system need only be passive up to a certain frequency; beyond that frequency the system can be stabilized with sufficiently small gains. With this concept the Center developed a design methodology for achieving desired compliant contact motions. This design method is based on H-infinity norm optimization, which makes it possible to consider both driving point mechanical impedance and systems robustness to modeling uncertainty. A laboratory facility was set up to verify experimentally the analytical and design theory.

Su, Renjeng↗
NASA NTRS1994

Flight Dynamics and Control of Elastic Hypersonic Vehicles Uncertainty Modeling

It has been shown previously that hypersonic air-breathing aircraft exhibit strong aeroelastic/aeropropulsive dynamic interactions. To investigate these, especially from the perspective of the vehicle dynamics and control, analytical expressions for key stability derivatives were derived, and an analysis of the dynamics was performed. In this paper, the important issue of model uncertainty, and the appropriate forms for representing this uncertainty, is addressed. It is shown that the methods suggested in the literature for analyzing the robustness of multivariable feedback systems, which as a prerequisite to their application assume particular forms of model uncertainty, can be difficult to apply on real atmospheric flight vehicles. Also, the extent to which available methods are conservative is demonstrated for this class of vehicle dynamics.

Chavez, Frank R.↗
DOE OSTI2025

Safe Physics-Informed Machine Learning for Dynamics and Control

This tutorial paper focuses on safe physics-informed machine learning in the context of dynamics and control, providing a comprehensive overview of how to integrate physical models and safety guarantees. As machine learning techniques enhance the modeling and control of complex dynamical systems, ensuring safety and stability remains a critical challenge, especially in safety-critical applications like autonomous vehicles, robotics, medical decision-making, and energy systems. We explore various approaches for embedding and ensuring safety constraints, including structural priors, Lyapunov and Control Barrier Functions, predictive control, projections, and robust optimization techniques. Additionally, we delve into methods for uncertainty quantification and safety verification, including reachability analysis and neural network verification tools, which help validate that control policies remain within safe operating bounds even in uncertain environments. The paper includes illustrative examples demonstrating the implementation aspects of safe learning frameworks that combine the strengths of data-driven approaches with the rigor of physical principles, offering a path toward the safe control of complex dynamical systems.

Drgona, Jan↗
NASA NTRS1992

On slowly varying systems - l-infinity to l-infinity performance and implications to robust adaptive control

The authors present a result on the l-infinity to l-infinity performance of slowly time varying systems. In particular, it is shown that the performance of such systems cannot be much worse than that of the frozen-time systems which are time invariant. This result is used to characterize a class of indirect adaptive controllers that can stabilize a time-invariant system subjected to both parametric and l-infinity to l-infinity bounded unstructured uncertainty. Pertaining to this class of controllers, a particular indirect adaptive scheme is proposed that provides the greatest upper bound on the size of the unstructured uncertainty for which stability is ensured.

Voulgaris, Petros G.↗
DOE OSTI2021

Stress- and interface-compatible red phosphorus anode for high-energy and durable sodium-ion batteries

Sodium-ion battery is a promising candidate for energy storage application, but the absence of high-capacity and low-cost anode materials significantly limit their practical specific energy and cost. Red phosphorus (RP) possesses a high theoretical specific capacity but suffers from large volume change, low electronic conductivity and unstable solid-electrolyte interphase (SEI). Herein, a hierarchical micro/nanostructured antimony-doped RP/carbon anode was developed, which demonstrate extraordinary electrochemical performance with high initial Coulombic efficiency of ~90%, high areal (~1.7 mAh cm-2) capacity together with good cycle stability and rate capability. Combined experimental and computational studies consistently revealed that such a unique structural design can dramatically accommodate the mechanical stress and moreover effectively restrain the undesired decomposition of electrolyte solvents regardless of electrolytes formulation, resulting in superior structural integrity and thin/robust SEI formation during cycling. The present finding has offered an alternative strategy for stress management and interface engineering on high-capacity alloying-based anode materials.

Liu, Xiang↗
DOE OSTI2021

Bilayer Wigner crystals in a transition metal dichalcogenide heterostructure

One of the first theoretically predicted manifestations of strong interactions in many-electron systems was the Wigner crystal, in which electrons crystallize into a regular lattice. The crystal can melt via either thermal or quantum fluctuations. Quantum melting of the Wigner crystal is predicted to produce exotic intermediate phases and quantum magnetism because of the intricate interplay of Coulomb interactions and kinetic energy. However, studying two-dimensional Wigner crystals in the quantum regime has often required a strong magnetic field or a moiré superlattice potential, thus limiting access to the full phase diagram of the interacting electron liquid. Here we report the observation of bilayer Wigner crystals without magnetic fields or moiré potentials in an atomically thin transition metal dichalcogenide heterostructure, which consists of two MoSe 2 monolayers separated by hexagonal boron nitride. We observe optical signatures of robust correlated insulating states at symmetric (1:1) and asymmetric (3:1, 4:1 and 7:1) electron doping of the two MoSe 2 layers at cryogenic temperatures. Here, we attribute these features to bilayer Wigner crystals composed of two interlocked commensurate triangular electron lattices, stabilized by inter-layer interaction. The Wigner crystal phases are remarkably stable, and undergo quantum and thermal melting transitions at electron densities of up to 6 × 10 12 per square centimetre and at temperatures of up to about 40 kelvin. Our results demonstrate that an atomically thin heterostructure is a highly tunable platform for realizing many-body electronic states and probing their liquid–solid and magnetic quantum phase transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling of thermal-ion-driven internal kink in DIII-D high- T i plasmas

Toroidal modeling utilizing the non-perturbative magnetohydrodynamic-kinetic hybrid code MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503), to a DIII-D high-performance discharge with the core thermal ion temperature of T i ~ 13 keV, finds multiple branches of unstable internal kink (IK) modes triggered by the bounce resonance of trapped thermal ions. The range of the computed mode frequency aligns well with the bounce frequency of trapped thermal ions, confirming the nature of resonance drive for the instability. The computed mode frequency is close to that observed in experiments. Plasma toroidal flow, in particular the fast flow ($ > $10% of Alfvén speed) as in the considered DIII-D discharge, strongly stabilizes the fluid IK branch but not all of the thermal-ion-driven internal kink (TIK) branches, leaving one dominant unstable TIK branch under the experimental conditions. Plasma resistivity also significantly affects stability of the TIKs (but not much on the IK which is ideally strongly unstable) in combination with the drift kinetic drive by thermal ions. High resistivity can drive additional unstable TIK branches. Scanning the on-axis safety factor over a large range confirms robustness of the TIK triggering by thermal ions. The computed TIKs contain both the IK and the tearing components as the dominant contributions to the eigenmode structure. Furthermore, these findings, together with the experimental results reported in reference (Du X.D. et al 2021 Phys. Rev. Lett. 127 025001), strongly suggest the possibility of thermal-ion driven instabilities in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Efficient and Robust Dynamic Simulation of Power Systems With Holomorphic Embedding

Dynamic simulation is vitally important in power system analysis, but traditional approaches based on numerical integration over small time steps are time-consuming. Also, the Newton-Raphson method suffers from difficulty in convergence when solving nonlinear algebraic equations. In this paper, we propose a novel dynamic simulation approach based on holomorphic embedding. By obtaining a high-order approximation of system dynamics, it achieves a much larger time step and thus enhances the computational efficiency significantly. In addition, the new approach avoids non-convergence issues in solving algebraic equations, which improves robustness. The approach includes flexible modeling of synchronous generators and controllers, and we propose a method for modeling generator coordinate transformations. The approach is tested on the IEEE 39-bus, 10-generator system and a Polish 2383-bus, 327-generator system. The results demonstrate promising computational efficiency and satisfactory numerical robustness for the analysis of large-scale power systems.

42 ENGINEERING↗
DOE OSTI2023

A geospatial risk analysis graphical user interface for identifying hazardous chemical emission sources

Background: Performing back trajectory and forward trajectory using the Hybrid Single-Particle Lagrangian Integrated Trajectory Model (HYSPLIT) is a reliable approach for assessing particle transport after release among mid-field atmospheric models. HYSPLIT has an externally facing online interface that allows non-expert users to run the model trajectories without requiring extensive training or programming. However, the existing HYSPLIT interface is limited if simulations have a large amount of meteorological data and timesteps that are not coincident. The objective of this study is to design and develop a more robust tool to rapidly evaluate hazard transport conditions and to perform risk analysis, while still maintaining an intuitive and user-friendly interface. Methods: HYSPLIT calculates forward and backward trajectories of particles based on wind speed, wind direction, and the corresponding location, timestamp, and Pasquill stability classes of the regions of the atmosphere in terms of the wind speed, the amount of solar radiation, and the fractional cloud cover. The computed particle transport trajectories, combined with the online Proton Transfer Reaction-Mass Spectrometry (PTR-MS) data (https://figshare.com/articles/dataset/ARL_Data_from_PROS_station_at_Hanford_site/19993964), can be used to identify and quantify the sources and affected area of the hazardous chemicals’ emission using the potential source distribution function (PSDF). PSDF is an improved statistical function based on the well-known potential source contribution function (PSCF) in establishing the air pollutant source and receptor relationship. Performing this analysis requires a range of meteorological and pollutant concentration measurements to be statistically meaningful. The existing HYSPLIT graphical user interface (GUI) does not easily permit computations of trajectories of a dataset of meteorological data in high temporal frequency. To improve the performance of HYSPLIT computations from a large dataset and enhance risk analysis of the accidental release of material at risk, a geospatial risk analysis tool (GRAT-GUI) is created to allow large data sets to be processed instantaneously and to provide ease of visualization. Results: The GRAT-GUI is a native desktop-based application and can be run in any Windows 10 system without any internet access requirements, thus providing a secure way to process large meteorological datasets even on a standalone computer. GRAT-GUI has features to import, integrate, and convert meteorological data with various formats for hazardous chemical emission source identification and risk analysis as a self-explanatory user interface. The tool is available at https://figshare.com/articles/software/GRAT/19426742.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Active flutter suppression using invariant zeros/eigensystem assignment

An active control law is synthesized to raise the flutter dynamic pressure boundary of an active flexible wing wind tunnel model. The multi-input/multi-output controller has a two-degree-of-freedom structure consisting of: (1) an output dynamic compensator selected to assign system invariant zeros to shape the multivariable root loci, and (2) an output gain feedback controller, based on eigensystem assignment theory, using the compensator augmented system to stabilize the flutter mode and optimize a specified performance index. Evaluation of a digital implementation of a constant gain controller in a general batch simulation of the system reveals that the feedback system is stable over a dynamic pressure range of 50-400 psf, indicating robustness of design to dynamic pressure variations.

Srinathkumar, S.↗
NASA NTRS2003

High-Time-Resolution Study of Magnetic Holes in the Solar Wind

The objectives of this investigation are to determine the internal plasma structure of kinetic-scale and larger scale magnetic holes, and to determine their stability, their source mechanism(s), and their spatial extent. It is also of importance to determine the relationship between kinetic-scale holes and long-duration holes. As smaller and smaller magnetic depressions are investigated in order to make this a complete study, a robust criterion is necessary for distinguishing magnetic holes from random or unresolvable fluctuations in the interplanetary magnetic field. In order to resolve this ambiguity, we obtained from the MFI experiments magnetic field measurements from the WIND spacecraft at a time resolution of 46 to 184 ms over certain periods. We have also devised a measure of certainty for magnetic hole detections. The certainty factor, q, is defined as the difference between the mean magnetic field in the hole and the local magnetic field, in units of the local standard deviation of the field strength. For fullest generality, it is necessary to calculate this q over the range of available scales of interest, from 60 ms up to 300 s. This technique establishes a two dimensional matrix of relative probabilities that a hole of some duration (d) might exist in the data set at a given time (t). In identifying q-peaks in time and duration, we also come upon a natural method for distinguishing holes with internal structure from multiple holes in close proximity or holes nested inside of others. If two q-peaks are more than a half-width apart, they are simply said to be separate events.

Lazarus, Alan↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2025

Unified Universal Control and Coordination of Inverter-Based Resources, and Validation for a PV + Battery Hybrid Plant

As renewable energy deployment grows, hybrid power plants (HPPs) combining photovoltaic (PV) and battery systems must evolve to offer both energy and grid stability services. These systems typically include a mix of grid-following (GFL) and grid-forming (GFM) inverters, presenting unique coordination and control challenges. This Department of Energy–funded project developed and validated a Unified Universal Control and Coordination (UUCC) framework for such PV + battery hybrid plants, enabling seamless and stable operation, including ultrafast black start, autonomous synchronization, and robust frequency and voltage regulation, under different grid conditions. The project significantly advanced the understanding of inverter-based resource (IBR) control by developing and validating three complementary system-level approaches for hybrid GFL/GFM operation: 1. A combined Virtual Resistance (VR)-based GFL and Virtual Oscillator Control (VOC)-based GFM method, where each inverter type is governed by a specialized control strategy. Together, these achieve stable, fast-response coordination, eliminating inrush current and enabling smooth black start and grid synchronization across a wide range of grid strengths. 2. A Deadbeat-based UUCC strategy, which uses discrete-time, switching-cycle-level control for both GFL and GFM inverters. This approach replaces traditional PI/PLL control with a control parameter-free, high-bandwidth framework that supports stable LVRT and instantaneous synchronization under all conditions. 3. A benchmark comparison with Siemens’ commercial GFM microgrid controller, which provided a fast baseline platform. The commercial approach decoupled v & f control was implemented on a commercial microgrid controller.The baseline commercial benchmark helped highlight superior transient response and black start performance offered by the deadbeat and VOC approaches. These technical contributions offer substantial improvements over conventional inverter control schemes, which often rely on slow phase-locked loop (PLL)-based synchronization, require careful control parameters tuning, and prone to unstable in weak grids with GFL inverters and in stiff grid with GFM inverters therefore challenging for hybrid GFL+GFM under all grid conditions. The deadbeat-based UUCC framework enables simpler, faster, and more robust operation of hybrid IBR systems using wide-bandgap (WBG) devices such as SiC power semiconductors. The rapid expansion of hybrid distributed energy resources (DERs), including residential and commercial PV-BESS installations such as Tesla Powerwall, PV with vehicle-to-grid (V2G) capability, and other integrated configurations, presents complex operational challenges for medium-voltage radial distribution feeders. These networks are subject to frequent disturbances such as faults, switching operations, rapid reclosing sequences, and feeder reconfigurations, all of which introduce dynamic stress on IBRs. In addition, planned feeder segmentation and deliberate islanding for resilience will require DERs that can autonomously perform blackstart, establish voltage and frequency references, and resynchronize with the main grid. The advanced deadbeat-based UUCC control and blackstart functionalities developed in this project directly address these requirements, enabling decentralized and autonomous operation of inverter-dominated DERs in distribution systems under a wide range of fault and reconfiguration scenarios. From a public benefit perspective, these innovations enable more reliable and cost-effective integration of renewable energy into distribution networks. The ability to autonomously black start and stabilize grids under varying grid conditions support accelerates recovery from outages and support decentralized resilient energy systems. By reducing system complexity and improving performance, this project lays critical groundwork for future inverter-dominated power grids that are clean, reliable, and accessible to all.

14 SOLAR ENERGY↗
DOE OSTI2023

Solvent-Mediated Formation of Quasi-2D Dion-Jacobson Phases on 3D Perovskites for Inverted Solar Cells Over 23% Efficiency

2D-on-3D (2D/3D) perovskite heterostructures present a promising strategy to realize efficient and stable photovoltaics. However, their applicability in inverted solar cells is limited due to the quantum confinement of the 2D-layer and solvent incompatibilities that disrupt the underlying 3D layer, hampering electron transport at the 2D/3D interface. Herein, solvent-dependent formation dynamics and structural evolution of 2D/3D heterostructures are investigated via in situ X-ray scattering. In this study, it is revealed that solvent interaction with the 3D surface determines the formation sequence and spatial distribution of quasi-2D phases with n = 2–4. Isopropanol (IPA) reconstructs the perovskite into a PbI 2 -rich surface, forming a strata with smaller n first, followed by a thinner substratum of larger n. In contrast, 2,2,2-Trifluoroethanol (TFE) preserves the 3D surface, promoting the formation of uniformly distributed larger n domains first, and smaller n last. Leveraging these insights, Dion–Jacobson perovskites are used with superior charge transport properties and structural robustness to fabricate 2D/3D heterostructures dominated by n ≥ 3 and engineer a favorable energy landscape for electron tunneling. Inverted solar cells based on 3-Aminomethylpyridine and TFE achieve a champion efficiency of 23.60%, with V oc and FF of 1.19 V and 84.5%, respectively, and superior stabilities with t 94 of 960 h under thermal stress.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

How non-ohmic contact-layer diodes in perovskite pinholes affect abrupt low-voltage reverse-bias breakdown and destruction of solar cells

Perovskite solar cells (PSCs) rapidly degrade under reverse bias, a condition that may occur during partial shading. Here, in this study, we use electrical measurements, electron microscopy, and optical and thermal imaging to investigate abrupt breakdown and hotspotting under low reverse potentials (<|-2| V). We show that microscopic pinholes in the perovskite layer cause rapid, destructive breakdown under reverse bias despite minimally reducing power conversion efficiencies. Measurements on miniature (200-micrometer diameter) PSCs and perovskite-free transport-layer diodes indicate that abrupt, low-voltage breakdown occurs in nanoscale to micrometer-scale defects and that metal migration and filamentation are unlikely causes. Reverse-bias stability substantially improves when pinholes in the perovskite and transport layers are eliminated. Atomic layer deposition of tin oxide prevents abrupt breakdown by ensuring physical separation between electrodes-not by blocking metal ion migration. Perovskite researchers should adopt cleaner, more uniform deposition techniques to enable robust PSCs for further research and commercial applications.

14 SOLAR ENERGY↗
DOE OSTI2022

Salton Sea Geothermal Development: Nontechnical Barriers to Entry – Analysis and Perspectives

Geothermal energy offers an opportunity to generate baseload, renewable energy that can help support the transition to an energy economy with reduced impacts on climate change and replace older, more expensive, nonrenewable, and more resource-impacting energy-generation facilities. The United States has the largest known geothermal resource in the world, with over 31 GW of conventional geothermal potential. However, due to market conditions, an inability to properly quantify both electrical grid benefits and resource stability, and the difficulty of exploring and developing the geothermal resource, few new geothermal projects have come online over the past three decades. The Salton Sea, in Imperial County, California, provides a prime location and opportunity to develop new geothermal resources. The Salton Sea contains a robust, well-mapped, geothermal resource, with opportunities for concurrent development of lithium and other mineral resources. This report describes the history of geothermal development at the Salton Sea and compares geothermal to other renewable energy sources in the area. The report then uses a techno-economic analysis (TEA) model to analyze the relative benefits and costs of various challenges and opportunities and provides recommendations for streamlining geothermal development at the Salton Sea and elsewhere. The challenges and opportunities analyzed in the TEA model were informed by stakeholder interviews and literature reviews. Based upon the identified challenges and opportunities and the results of the TEA model, primary findings are that certain nontechnical barriers such as permitting costs play only a minor role in determining the viability of development of the geothermal resource at the Salton Sea. Other barriers such as permitting timelines, government/agency coordination, and the potential co-location of lithium extraction with a geothermal plant may result in much larger impacts on project viability.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Accelerated Selectrion of Optimal Perovskite Alloys for Solar PV using a Combined Quantum and Machine Learning Hierachiral Approach

The project aims to: (i) accelerate the discovery of “Missing HP alloys” by combining quantum mechanics and artificial intelligence machine learning approaches, and (ii) analyze the stabilities of candidate alloys, including those that do not pass selection filters (and are thus expected to degrade over time) to decipher the nature of the instabilities to guide the development of durable solar cell materials. Successful candidates will be subjected to validation experiments at NREL's state-of-the-art facilities. The discoveries this effort will provide will be directly testable and implementable and will greatly impact U.S. progress in HP PV as they will provide clear direction and motivation for experimental studies including specific material synthetic targets, device optimization, and device stability protocols. A key advantage of this effort is the feedback and guidance provided by the Industry Collaborative Work Group that we established to coordinate academic and national lab research with industry needs. The proposed work will provide a basis for and direct the development of robust and reliable HP PV. It will also provide a timely, valuable and extensive roadmap to the experimental HP PV community to enable it to focus its efforts on improving and fine-tuning promising HP compositions that this effort predicts will likely be the best performers rather than wandering in the vast chemical space for decades spending enormous resources mostly evaluating unpromising candidate materials.

14 SOLAR ENERGY↗
NASA NTRS1989

Algorithms for computing the multivariable stability margin

Stability margin for multiloop flight control systems has become a critical issue, especially in highly maneuverable aircraft designs where there are inherent strong cross-couplings between the various feedback control loops. To cope with this issue, we have developed computer algorithms based on non-differentiable optimization theory. These algorithms have been developed for computing the Multivariable Stability Margin (MSM). The MSM of a dynamical system is the size of the smallest structured perturbation in component dynamics that will destabilize the system. These algorithms have been coded and appear to be reliable. As illustrated by examples, they provide the basis for evaluating the robustness and performance of flight control systems.

Tekawy, Jonathan A.↗