Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex systems theory”

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

Dynamic Modeling, Trajectory Optimization, and Linear Control of Cable-Driven Parallel Robots for Automated Panelized Building Retrofits

The construction industry faces a growing need for automation to reduce costs, improve accuracy and productivity, and address labor shortages. One area that stands to benefit significantly from automation is panelized prefabricated building envelope retrofits, which can improve a building’s energy efficiency in heating and cooling interior spaces. In this paper, we propose using cable-driven parallel robots (CDPRs), which can effectively lift and handle large objects, to install these panels. However, implementing CDPRs presents significant challenges because of their nonlinear dynamics, complex trajectory planning, and precise control requirements. To tackle these challenges, this work focuses on a new application of established control and trajectory optimization theories in a CDPR simulation of a building envelope retrofit under real-world conditions. We first model the dynamics of CDPRs, highlighting the critical role of damping in system behavior. Building on this dynamic model, we formulate a trajectory optimization problem to generate feasible and efficient motion plans for the robot under operational and environmental constraints. Given the high precision required in the construction industry, accurately tracking the optimized trajectory is essential. However, challenges such as partial observability and external vibrations complicate this task. To address these issues, a Linear Quadratic Gaussian control framework is applied, enabling the robot to track the optimized trajectories with precision. Simulation results show that the proposed controller enables precise end effector positioning with errors under 4 mm, even in the presence of external wind disturbances. Through comprehensive simulations, our approach allows for an in-depth exploration of the system’s nonlinear dynamics, trajectory optimization, and control strategies under controlled yet highly realistic conditions. The results demonstrate the feasibility of CDPRs for automating panel installation and provide insights into their practical deployment.

CDPR↗
NASA NTRS1972

General chemical kinetics computer program for static and flow reactions, with application to combustion and shock-tube kinetics

A general chemical kinetics program is described for complex, homogeneous ideal-gas reactions in any chemical system. Its main features are flexibility and convenience in treating many different reaction conditions. The program solves numerically the differential equations describing complex reaction in either a static system or one-dimensional inviscid flow. Applications include ignition and combustion, shock wave reactions, and general reactions in a flowing or static system. An implicit numerical solution method is used which works efficiently for the extreme conditions of a very slow or a very fast reaction. The theory is described, and the computer program and users' manual are included.

Bittker, D. A.↗
NASA NTRS1989

Minimal complexity control law synthesis

A paradigm for control law design for modern engineering systems is proposed: Minimize control law complexity subject to the achievement of a specified accuracy in the face of a specified level of uncertainty. Correspondingly, the overall goal is to make progress towards the development of a control law design methodology which supports this paradigm. Researchers achieve this goal by developing a general theory of optimal constrained-structure dynamic output feedback compensation, where here constrained-structure means that the dynamic-structure (e.g., dynamic order, pole locations, zero locations, etc.) of the output feedback compensation is constrained in some way. By applying this theory in an innovative fashion, where here the indicated iteration occurs over the choice of the compensator dynamic-structure, the paradigm stated above can, in principle, be realized. The optimal constrained-structure dynamic output feedback problem is formulated in general terms. An elegant method for reducing optimal constrained-structure dynamic output feedback problems to optimal static output feedback problems is then developed. This reduction procedure makes use of star products, linear fractional transformations, and linear fractional decompositions, and yields as a byproduct a complete characterization of the class of optimal constrained-structure dynamic output feedback problems which can be reduced to optimal static output feedback problems. Issues such as operational/physical constraints, operating-point variations, and processor throughput/memory limitations are considered, and it is shown how anti-windup/bumpless transfer, gain-scheduling, and digital processor implementation can be facilitated by constraining the controller dynamic-structure in an appropriate fashion.

Bernstein, Dennis S.↗
NASA NTRS1997

Mixed H2/H-Infinity Control of a Flexible Space Structure

As theory progresses for design and analysis of robust multivariable control laws, synthesis procedures and to a larger extent, experimental verification generally lags behind. Recent developments in robust control theory have extended the H-infinity and mu-synthesis methods to incorporate H2 properties in the control synthesis. A major difficulty in implementing robust controllers is the order of the compensator and associated complexity of the computation required for synthesis, especially when order constraints are imposed. This paper presents results of system identification and robust control design for the Controls/Structures Interaction Ground Test Facility at NASA/Marshall Space Flight Center.

Whorton, Mark↗
DOE OSTI2023

The theory of kinetic effects on resistive wall mode stability in tokamaks

Tokamak fusion plasmas benefit from high pressures but are then susceptible to modes of instability. These magnetohydrodynamic (MHD) modes are macroscopic distortions of the plasma, but certain collective motions of individual particles can provide stabilizing effects opposing them. The presence of a resistive wall slows the mode growth, converting a kink to a resistive wall mode (RWM). A kinetic MHD model includes Maxwell's equations, ideal MHD constraints, and kinetic effects included through the pressure tensor, calculated with the perturbed drift-kinetic distribution function of the particles. The kinetic stabilizing effects on the RWM arise through resonances between the plasma rotation and particle drift motions: precession, bounce, and transit. A match between particle motions and the mode allows efficient transfer of energy that would otherwise drive the growth of the mode, thus damping the growth. The first approach to calculating RWM stability is to write a set of equations for the complex mode frequency in terms of known quantities and then to solve the system. The “energy principle” approach, which has the advantage of clarity in distinguishing the various stabilizing and destabilizing effects, is to change the force balance equation into an equation in terms of changes of kinetic and potential energies, and then to write a dispersion relation for the mode frequency in terms of those quantities. These methods have been used in various benchmarked codes to calculate kinetic effects on RWM stability. Importantly, the theory has illuminated the important roles of plasma rotation, energetic particles, and collisions in RWM stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗
NASA NTRS1993

Flight control application of new stability robustness bounds for linear uncertain systems

This paper addresses the issue of obtaining bounds on the real parameter perturbations of a linear state-space model for robust stability. Based on Kronecker algebra, new, easily computable sufficient bounds are derived that are much less conservative than the existing bounds since the technique is meant for only real parameter perturbations (in contrast to specializing complex variation case to real parameter case). The proposed theory is illustrated with application to several flight control examples.

Yedavalli, Rama K.↗
DOE OSTI2020

New H2O–H2O collisional rate coefficients for cometary applications

ABSTRACT We re-introduce a semiclassical methodology based on theories developed for the determination of broadening coefficients. We show that this simple and extremely fast methodology provides results that are in good agreement with results obtained using the more sophisticate MQCT approach. This semiclassical methodology could be an alternative approach which allows to provide large sets of collisional data for very complex molecular systems. It saves time both on the determination of potential energy surfaces and on the collisional dynamical calculations. In addition, this paper provides more complete sets of rotational de-excitation cross-sections and rate coefficients of H2O perturbed by a thermal average of water molecules. Those data can be used in the radiative transfer modelling of cometary atmospheres.

Boursier, C.↗
DOE OSTI2022

Strong suppression of heat conduction in a laboratory replica of galaxy-cluster turbulent plasmas

In conventional gases and plasmas, it is known that heat fluxes are proportional to temperature gradients, with collisions between particles mediating energy flow from hotter to colder regions and the coefficient of thermal conduction given by Spitzer’s theory. However, this theory breaks down in magnetized, turbulent, weakly collisional plasmas, although modifications are difficult to predict from first principles due to the complex, multiscale nature of the problem. Understanding heat transport is important in astrophysical plasmas such as those in galaxy clusters, where observed temperature profiles are explicable only in the presence of a strong suppression of heat conduction compared to Spitzer’s theory. To address this problem, we have created a replica of such a system in a laser laboratory experiment. Our data show a reduction of heat transport by two orders of magnitude or more, leading to large temperature variations on small spatial scales (as is seen in cluster plasmas).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Hydrological connectivity: a review and emerging strategies for integrating measurement, modeling, and management

This review synthesizes methods for measuring, modeling, and managing hydrologic connectivity, offering pathways to improve practices and address environmental challenges (e.g., climate change) and sustainability. As a key driver of water movement and nutrient cycling, hydrologic connectivity influences flood mitigation, water quality regulation, and biodiversity conservation. However, traditional field-based methods (e.g., dye tracing), indirect measurements (e.g., runoff analysis), and remote sensing techniques (e.g., InSAR) often struggle to capture the complexity of catchment-scale interactions. Similarly, modeling approaches—including process-based and percolation theory-based models, graph theory, and entropy-based metrics—face limitations in fully representing these interconnected processes. Both modeling and measurement techniques are constrained by inadequate spatial and temporal coverage, high data demands, computational complexity, and difficulties in representing subsurface connectivity. Subsequently, we critique current management practices that prioritize isolated variables (e.g., streamflow, sediment transport) over system-wide strategies and emphasize the need for adaptive, connectivity-based approaches in water resource planning and restoration. Moving forward, we highlight the importance of interdisciplinary collaboration, technological innovations (e.g., AI-driven modeling, real-time monitoring), and integrated frameworks to improve connectivity measurement, modeling, and adaptive management to restore fragmented hydrologic networks. This integrated approach sets the stage for transformative water resource management, fostering proactive policy development and stakeholder engagement.

Dwivedi, Dipankar↗
NASA NTRS1989

Parametric Cost Analysis: A Design Function

Parametric cost analysis uses equations to map measurable system attributes into cost. The measures of the system attributes are called metrics. The equations are called cost estimating relationships (CER's), and are obtained by the analysis of cost and technical metric data of products analogous to those to be estimated. Examples of system metrics include mass, power, failure_rate, mean_time_to_repair, energy _consumed, payload_to_orbit, pointing_accuracy, manufacturing_complexity, number_of_fasteners, and percent_of_electronics_weight. The basic assumption is that a measurable relationship exists between system attributes and the cost of the system. If a function exists, the attributes are cost drivers. Candidates for metrics include system requirement metrics and engineering process metrics. Requirements are constraints on the engineering process. From optimization theory we know that any active constraint generates cost by not permitting full optimization of the objective. Thus, requirements are cost drivers. Engineering processes reflect a projection of the requirements onto the corporate culture, engineering technology, and system technology. Engineering processes are an indirect measure of the requirements and, hence, are cost drivers.

Dean, Edwin B.↗
DOE OSTI2021

Topological Machine Learning Methods for Power System Responses to Contingencies: Preprint

While deep learning tools, coupled with the emerging machinery of topological data analysis, are proven to deliver various performance gains in a broad range of applications, from image classification to biosurveillance to blockchain fraud detection, their utility in areas of high societal importance such as power system modeling and, particularly, resilience quantification in the energy sector yet remain untapped. To provide fast acting synthetic regulation and contingency reserve services to the grid while having minimal disruptions on customer quality of service, we propose a new topologybased system that depends on neural network architecture for impact metrics classification and prediction in power systems. This novel topology-based system allows one to evaluate the impact of three power system contingency types, namely, in conjunction with transmission lines, transformers, and transmission lines combined with transformers. We show that the proposed new neural network architecture equipped with local topological measures facilitates both more accurate classification of unserved load as well as the amount of unserved load. In addition, we are able to learn complex relationships between electrical properties and local topological measurements on the simulated response to contingencies for NREL-SIIP power system.

graph theory↗
DOE OSTI2025

Multi-modal characterization of the B2 phase in the Ta-Re binary system

The energy and transportation industries demand materials that retain their mechanical property at high temperatures. Refractory complex concentrated alloys (RCCAs) with a BCC + B2 microstructure offer a potential solution, where maintaining the high temperature mechanical properties can be achieved by precipitation strengthening. This depends on the B2 phase in RCCAs being thermodynamically stable with a high solvus temperature. Recently, we predicted the high temperature stability of the B2 structure in the Ta-Re binary system, using density functional theory. Here, in this study, we provide experimental evidence for the existence of this phase for the first time, using a Ta 65 Re 35 alloy. Despite Ta-Re binary phase diagrams predicting a single-phase BCC microstructure for Ta 65 Re 35 , we show that a high Z nanoscale secondary phase appears after heat treatment at 1550 °C and 1100 °C. Scanning transmission electron microscopy (STEM) revealed that this phase has a cubic structure and is equiatomic TaRe though B2 superlattice reflections were absent in fast Fourier transforms (FFT) and diffraction patterns (DPs). DP simulations indicate that the B2 TaRe superlattice reflections are up to two orders of magnitude weaker than their fundamental reflections, making their detection challenging via electron microscopy. Neutron diffraction confirmed the second phase had a B2 structure. This study identified a previously unobserved high temperature stable B2 phase in the Ta-Re system, enabling the development of new high temperature BCC + B2 RCCAs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sterile neutrinos, black hole vacuum and holographic principle

We construct an effective field theory (EFT) model that describes matter field interactions with Schwarzschild mini-black-holes (SBH’s), treated as a scalar field, B 0 (x). Fermion interactions with SBH’s require a complex spurion field, θ ij , which we interpret as the EFT description of “holographic information,” which is correlated with the SBH as a composite system. We consider Hawking’s virtual black hole vacuum (VBH) as a Higgs phase, $\langle$B 0 $\rangle$=V. Integrating sterile neutrino loops, the information field θij is promoted to a dynamical field, necessarily developing a tachyonic instability and acquiring a VEV of order the Planck scale. For N sterile neutrinos this breaks the vacuum to SU(N) × U(1)/SO(N) with N degenerate Majorana masses, and 12N(N+1) Nambu-Goldstone neutrino-Majorons. The model suggests many scalars fields, corresponding to all fermion bilinears, may exist bound nonperturbatively by gravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Quantum Choreography of the Nucleus: Rotations, Vibrations, and Emergent Structure

Abstract Nuclei are complex many-body quantum systems where interactions of the neutrons and protons via the strong, the weak, and the electromagnetic forces lead to the emergence of simple patterns of energy states that have been described by various theoretical approaches. One of the goals of all the theoretical models is the development of a universal theory that can be applied across the entire chart of nuclides. Significant progress has been made by experiments as well as the increasing sophistication of models, but a universal theory has yet to be established. A recent review of nuclei in the Z = 50–82 region of the chart of nuclides has analyzed all the available compiled data from several decades of studies towards a clarification of the low-lying structure of nuclei. Other reviews have reported and explained the emergence of multiple different shapes in nuclei at somewhat higher excitation energies than the ground state. Some have challenged the interpretation of the first excited K π = 0 + band as a vibration of ground state. This work attempts to provide a guide to determining the nature of the first excited K π = 0 + band in nuclei by the combined use of nuclear lifetimes, energy level evolutions, dynamic moments of inertia, and intrinsic quadrupole moments extracted from transition probabilities. The result is that for a subset of the nuclei in this region, the K π = 0 + band is consistent with the traditional 𝛽-vibration description of an oscillation built on the ground state.

CHFB+5DCH↗
NASA NTRS1992

Calculation of Mg(+)-ligand relative binding energies

The calculated relative binding energies of 16 organic molecules to Mg(+) are compared with experimental results where available. The geometries of the ligands and the Mg(+)-ligand complexes arc optimized at the self-consistent field level using a 6-31G* basis set. The Mg(+) binding energies are evaluated using second-order perturbation theory and basis sets of triple-sigma quality augmented with two sets of polarization functions. This level of theory is calibrated against higher levels of theory for selected systems. The computed binding energies are accurate to about 2 kcal/mol.

Partridge, Harry↗
NASA NTRS2008

Control Architecture for Robotic Agent Command and Sensing

Control Architecture for Robotic Agent Command and Sensing (CARACaS) is a recent product of a continuing effort to develop architectures for controlling either a single autonomous robotic vehicle or multiple cooperating but otherwise autonomous robotic vehicles. CARACaS is potentially applicable to diverse robotic systems that could include aircraft, spacecraft, ground vehicles, surface water vessels, and/or underwater vessels. CARACaS incudes an integral combination of three coupled agents: a dynamic planning engine, a behavior engine, and a perception engine. The perception and dynamic planning en - gines are also coupled with a memory in the form of a world model. CARACaS is intended to satisfy the need for two major capabilities essential for proper functioning of an autonomous robotic system: a capability for deterministic reaction to unanticipated occurrences and a capability for re-planning in the face of changing goals, conditions, or resources. The behavior engine incorporates the multi-agent control architecture, called CAMPOUT, described in An Architecture for Controlling Multiple Robots (NPO-30345), NASA Tech Briefs, Vol. 28, No. 11 (November 2004), page 65. CAMPOUT is used to develop behavior-composition and -coordination mechanisms. Real-time process algebra operators are used to compose a behavior network for any given mission scenario. These operators afford a capability for producing a formally correct kernel of behaviors that guarantee predictable performance. By use of a method based on multi-objective decision theory (MODT), recommendations from multiple behaviors are combined to form a set of control actions that represents their consensus. In this approach, all behaviors contribute simultaneously to the control of the robotic system in a cooperative rather than a competitive manner. This approach guarantees a solution that is good enough with respect to resolution of complex, possibly conflicting goals within the constraints of the mission to be accomplished by the vehicle(s).

Huntsberger, Terrance↗
DOE OSTI2026

Structural and compositional complexities of hierarchical self-assembly: A hypergraph approach

Programmable self-assembly enables the construction of complex molecular, supramolecular, and crystalline architectures from well-designed building blocks. In this work, we introduce a hypergraph-based formalism, Blocks & Bonds (B&B), which generalizes classical chemical graph theory by incorporating directed and multicolored interactions, internal symmetries, and hierarchical organization. Within this framework, we develop the Structure Code (SC), a compact and versatile language for describing self-assembled architectures. We define a Kolmogorov-style structural complexity as the total information content of SC, obtained through its tokenization and Shannon information assignment. Complementing this encoding-based measure, we introduce a much simpler quantity, the compositional complexity, which depends only on the number and cumulative usage of block and bond types in the construction set. A central result of this work is a strong empirical correlation between the token-based structural complexity and the compositional complexity across all examined systems. Owing to this agreement, the compositional complexity emerges as the most practical and broadly applicable measure: it is easy to compute, requires no explicit encoding, and yet closely tracks the actual information content of structurally diverse architectures. Applications to molecular systems (ethylene glycol and glucose), DNA-origami lattices, and crystalline assemblies show that B&B hypergraphs provide a unified, scalable, and information-efficient representation of structural organization, naturally capturing symmetry, modularity, and stereochemistry. This framework establishes a quantitative foundation for complexity-aware classification and inverse design of programmable matter.

36 MATERIALS SCIENCE↗