Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained 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.

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

At least 37 records · Page 2

Adaptive Multi-Agent Systems for Constrained Optimization

Product Distribution (PD) theory is a new framework for analyzing and controlling distributed systems. Here we demonstrate its use for distributed stochastic optimization. First we review one motivation of PD theory, as the information-theoretic extension of conventional full-rationality game theory to the case of bounded rational agents. In this extension the equilibrium of the game is the optimizer of a Lagrangian of the (probability distribution of) the joint state of the agents. When the game in question is a team game with constraints, that equilibrium optimizes the expected value of the team game utility, subject to those constraints. The updating of the Lagrange parameters in the Lagrangian can be viewed as a form of automated annealing, that focuses the MAS more and more on the optimal pure strategy. This provides a simple way to map the solution of any constrained optimization problem onto the equilibrium of a Multi-Agent System (MAS). We present computer experiments involving both the Queen s problem and K-SAT validating the predictions of PD theory and its use for off-the-shelf distributed adaptive optimization.

Macready, William↗

The expanded Lagrangian system for constrained optimization problems

Smooth penalty functions can be combined with numerical continuation/bifurcation techniques to produce a class of robust and fast algorithms for constrained optimization problems. The key to the development of these algorithms is the Expanded Lagrangian System which is derived and analyzed in this work. This parameterized system of nonlinear equations contains the penalty path as a solution, provides a smooth homotopy into the first-order necessary conditions, and yields a global optimization technique. Furthermore, the inevitable ill-conditioning present in a sequential optimization algorithm is removed for three penalty methods: the quadratic penalty function for equality constraints, and the logarithmic barrier function (an interior method) and the quadratic loss function (an interior method) for inequality constraints. Although these techniques apply to optimization in general and to linear and nonlinear programming, calculus of variations, optimal control and parameter identification in particular, the development is primarily within the context of nonlinear programming.

Poore, A. B.↗

The expanded LaGrangian system for constrained optimization problems

Smooth penalty functions can be combined with numerical continuation/bifurcation techniques to produce a class of robust and fast algorithms for constrainted optimization problems. The key to the development of these algorithms is the Expanded Lagrangian System which is derived and analyzed in this work. This parameterized system of nonlinear equations contains the penalty path as a solution, provides a smooth homotopy into the first-order necessary conditions, and yields a global optimization technique. Furthermore, the inevitable ill-conditioning present in a sequential optimization algorithm is removed for three penalty methods: the quadratic penalty function for equality constraints, and the logarithmic barrier function (an interior method) and the quadratic loss function (an interior method) for inequality constraints. Although these techniques apply to optimization in general and to linear and nonlinear programming, calculus of variations, optimal control and parameter identification in particular, the development is primarily within the context of nonlinear programming.

Poore, A. B.↗

Using Observations of Deep Convective Systems to Constrain Atmospheric Column Absorption of Solar Radiation in the Optically Thick Limit

Atmospheric column absorption of solar radiation A(sub col) is a fundamental part of the Earth's energy cycle but is an extremely difficult quantity to measure directly. To investigate A(sub col), we have collocated satellite-surface observations for the optically thick Deep Convective Systems (DCS) at the Department of Energy Atmosphere Radiation Measurement (ARM) Tropical Western Pacific (TWP) and Southern Great Plains (SGP) sites during the period of March 2000 December 2004. The surface data were averaged over a 2-h interval centered at the time of the satellite overpass, and the satellite data were averaged within a 1 deg X 1 deg area centered on the ARM sites. In the DCS, cloud particle size is important for top-of-atmosphere (TOA) albedo and A(sub col) although the surface absorption is independent of cloud particle size. In this study, we find that the A(sub col) in the tropics is approximately 0.011 more than that in the middle latitudes. This difference, however, disappears, i.e., the A(sub col) values at both regions converge to the same value (approximately 0.27 of the total incoming solar radiation) in the optically thick limit (tau greater than 80). Comparing the observations with the NASA Langley modified Fu_Liou 2-stream radiative transfer model for optically thick cases, the difference between observed and model-calculated surface absorption, on average, is less than 0.01, but the model-calculated TOA albedo and A(sub col) differ by 0.01 to 0.04, depending primarily on the cloud particle size observation used. The model versus observation discrepancies found are smaller than many previous studies and are just within the estimated error bounds. We did not find evidence for a large cloud absorption anomaly for the optically thick limit of extensive ice cloud layers. A more modest cloud absorption difference of 0.01 to 0.04 cannot yet be ruled out. The remaining uncertainty could be reduced with additional cases, and by reducing the current uncertainty in cloud particle size.

Dong, Xiquan↗

Constraining Solar System Bombardment Using In Situ Radiometric Dating

The leading, but contentious, model for lunar impact history includes a pronounced increase in impact events at around 3.9 Ga. This late heavy bombardment would have scarred Mars and the terrestrial planets, influenced the course of biologic evolution on the early Earth, and rearranged the very architecture of our Solar System. But what if it's not true? In the last decade, new observations and sample analyses have reinterpreted basin ages and "pulled the pin" on the cataclysm - we may only have the age of one large basin (Imbrium). The Curie mission would constrain the onset of the cataclysm by determining the age of a major pre-Imbrium lunar basin (Nectaris or Crisium), characterize new lunar lithologies far from the Apollo and Luna landing sites, including the basalts in the basin-filling maria and olivine-rich lithologies in the basin margins, and provide a unique vantage point to assess volatiles in the lunar regolith from dawn to dusk.

Cohen, Barbara A.↗

Application Deconfliction Characterization and Alternatives Analysis

This report provides an overview of the domain space and solution techniques that could be used to create a robust, flexible app deconfliction service. Three approaches are reviewed with summaries of the characteristics, elements, and results from preliminary demonstrations of solution techniques based on each approach: 1) rules and heuristics, 2) cooperation, and 3) optimization. The strengths and weaknesses of each solution technique were explored through a set of numerical demonstrations on modified IEEE 123 node and 9500 node test feeders. An alternatives analysis of individual deconfliction elements was performed with each solution technique element evaluated against criteria reflecting the dynamic app environment, need to balance app objectives, and scalability issues versus the number of applications, setpoints, and distributed control areas. It is anticipated that a combined solution for a GridAPPS-D Deconfliction Service can be formulated using a combination of elements from each solution technique. The combined solution would combine 1) device control budgets to reduce the size of the solution space by constraining system setpoints to those will not result in accelerated degradation of physical assets, 2) system operations rules to constrain the solution space by eliminating setpoints that result in violations of system limits or operational best practices, 3) contextual status signals shared with or among apps such that they could update their desired setpoints based on the evolving context, 4) a mediator that incentivizes apps to come to a cooperative solution, and 5) Setpoint-informed optimization as a fallback mechanism if a cooperative solution cannot be agreed upon by applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Stochastic Constrained Extended System Dynamics for Solving Charge Equilibration Models

In this work, we present a new stochastic extended Lagrangian molecular dynamics solution to charge equilibration that eliminates self-consistent field (SCF) calculations, thus eliminating the computational bottleneck in solving the charge distribution with standard SCF solvers. By formulating both charges and chemical potential as latent variables and introducing a holonomic constraint that satisfies charge conservation, the SC-XLMD method accurately reproduces thermodynamic, dynamic, and structural properties within the framework of ReaxFF for a bulk water system and highly reactive RDX molecules simulated at high temperature. The SC-XLMD method shows excellent computational performance and is available in the publicly available LAMMPS package.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TRIM: AI Guided Random Number Generation for Resource-Constrained IoT Systems

Random numbers often serve as the backbone for many security solutions in diverse domains such as cryptography, side channel leakage prevention, and moving target defense. However, generating true random numbers requires a physical source of entropy (e.g. hardware, quantum, environmental phenomenon) making it difficult to realize at a large scale and at a low cost. On the flip side, pseudorandom number generators (easy to implement) following a specific distribution (e.g. Gaussian) can be easily compromised given a sufficient amount of traces. In this work, we have developed a machine learning-guided generative approach that can be used to create portable, resource-efficient, and cost-effective random number generators with high throughput and true randomness characteristics. We implement the proposed approach as a highly parameterized framework and perform extensive evaluation for different settings. The framework was able to learn from true random sources such as irrational numbers and environmental audio noise and imitate those sources towards generating new good quality random numbers on demand. We have generated more than 1 billion bits and observed robust performance in terms of true randomness metrics obtained from NIST SP 800-22 and FIPS 140-1 randomness test suites achieving a throughput of up to 142.85 Mbps. Compared to the state-of-the-art (SOTA) technique, the iso-cost setup of our framework can achieve more than 500 Mbps in a distributed setting. We have evaluated the efficacy of running the true randomness imitation AI models on target edge devices such as Raspberry Pi 4 (Model B), Nvidia Jetson Nano, Nvidia Jetson Orin Nano and Nvidia Jetson Xavier. We have also looked at the security of the TRIM framework itself against different adversarial threat models.

Cybersecurity↗

Revisiting the Effects of Photovoltaic Module Efficiency on Installed System Cost and Market Adoption

We conduct bottom-up analysis of how PV module efficiency could influence residential rooftop installed system costs, using models with improved representation of how soft costs and power electronics are affected by efficiency. We show that the amount of system installed cost savings resulting from high efficiency modules depends on whether or not the system is area-constrained, with much larger savings for the area-constrained case. We then estimate the fraction of physically area-constrained rooftops in the United States to be 61%, with variation between 20% to over 80% of rooftops depending on the state. We find that fraction of area-constrained systems could grow to over 88% by 2050 if building and transportation loads are electrified.

efficiency↗

Stabilization of computational procedures for constrained dynamical systems

A new stabilization method of treating constraints in multibody dynamical systems is presented. By tailoring a penalty form of the constraint equations, the method achieves stabilization without artificial damping and yields a companion matrix differential equation for the constraint forces; hence, the constraint forces are obtained by integrating the companion differential equation for the constraint forces in time. A principal feature of the method is that the errors committed in each constraint condition decay with its corresponding characteristic time scale associated with its constraint force. Numerical experiments indicate that the method yields a marked improvement over existing techniques.

Park, K. C.↗

Time variant analysis of large scale constrained rotorcraft systems dynamics - An exploitation of IBM-3090 vector-processor's pipe-lining feature

A generalized algorithmic procedure is presented for handling the constraints in transmissions, which are treated as a multibody system of interconnected rigid/flexible bodies. The type of constraints are classified based on the interconnection of the bodies, assuming one or more points of contact to exist between them. The method is explained through flow charts and configuration/interaction tables. A significant increase in speed of execution is achieved by vectorizing the developed code in computationally intensive areas. The study of an example consisting of two meshing disks rotating at high angular velocity is carried out. The dynamic behavior of the constraint forces associated with the generalized coordinates of the system are plotted by selecting various modes. Applications are intended for the study of dynamic and subsequent prediction of constraint forces at the gear teeth contacting points in helicopter transmissions with the aim of improving performance dependability.

Amirouche, F. M. L.↗

Introducing Kynema, an Open-Source Performance-Portable Flexible-Multibody-Dynamics Solver

In this talk we introduce Kynema, an open-source general flexible-multibody-dynamics solver that is well suited for simulating wind turbine structural dynamics. Kynema uses a Lie-group time integrator for constrained systems and runs on both CPUs and GPUs. Timing results for simulations are presented for the IEA 15-MW turbine with and without aerodynamic forces.

17 WIND ENERGY↗