SEARCH · Engineering Papers

Results for “nonlinear constraints”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Resolution enhancement of multichannel microwave imagery from the Nimbus-7 SMMR for maritime rainfall analysis

A restoration of the 37, 21, 18, 10.7, and 6.6 GHz satellite imagery from the scanning multichannel microwave radiometer (SMMR) aboard Nimbus-7 to 22.2 km resolution is attempted using a deconvolution method based upon nonlinear programming. The images are deconvolved with and without the aid of prescribed constraints, which force the processed image to abide by partial a priori knowledge of the high-resolution result. The restored microwave imagery may be utilized to examine the distribution of precipitating liquid water in marine rain systems.

Olson, William S.↗

NASA NTRS1990

The dynamics of orbital maneuvering: Design and evaluation of a visual display aid for human controllers

An interactive proximity operations planning system, which allows on-site planning of fuel-efficient, multi-burn maneuvers in a potential multi-spacecraft environment was developed and tested. Though this display system most directly assists planning by providing visual feedback to aid visualization of the trajectories and constraints, its most significant features include an inverse dynamics algorithm that removes control nonlinearities facing the operator and a trajectory planning technique that reduces the order of control and creates, through a geometric spread-sheet the illusion of an inertially stable environment. This synthetic environment provides the user with control of relevant static and dynamic properties of way-points during small orbital changes allowing independent solutions to the normally coupled problems of orbital maneuvering. An experiment was carried out in which experienced operators were required to plan a trajectory to retrieve an object accidently separated from a dual-keel space station. The time required to plan these maneuvers was found to be predicted by the direction of the insertion thrust and did not depend on the point of separation from the space station.

Ellis, Stephen R.↗

NASA NTRS1997

Enhanced Attitude Control Experiment for SSTI Lewis Spacecraft

The enhanced attitude control system experiment is a technology demonstration experiment on the NASA's small spacecraft technology initiative program's Lewis spacecraft to evaluate advanced attitude control strategies. The purpose of the enhanced attitude control system experiment is to evaluate the feasibility of designing and implementing robust multi-input/multi-output attitude control strategies for enhanced pointing performance of spacecraft to improve the quality of the measurements of the science instruments. Different control design strategies based on modern and robust control theories are being considered for the enhanced attitude control system experiment. This paper describes the experiment as well as the design and synthesis of a mixed H(sub 2)/H(sub infinity) controller for attitude control. The control synthesis uses a nonlinear programming technique to tune the controller parameters and impose robustness and performance constraints. Simulations are carried out to demonstrate the feasibility of the proposed attitude control design strategy. Introduction

Maghami, Peoman G.↗

DOE OSTI2024

Advanced Features and Recent Developments in the WEC-Sim Open-Source Design Tool: Preprint

WEC-Sim is an open-source tool developed in Matlab/Simulink to simulate wave energy converters. WEC-Sim has been developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (SNL) since its initial release in 2014. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). Within this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. WEC-Sim functionalities are developed using Simulink blocks and Matlab functions. WEC-Sim is capabilities include simulation of different wave conditions, constraints and PTOs, body dynamics, mooring systems, among other advanced features such as nonlinear hydrodynamics, nonlinear PTOs, body-to-body interactions, etc.

hydrodynamics modeling↗

DOE OSTI2024

Cool dark sector, concordance, and a low σ 8

We investigate a cosmological model in which a fraction of the dark matter is atomic dark matter (ADM). This ADM consists of dark versions of the electron and of the proton, interacting with each other and with dark photons just as their light sector versions do, but interacting with everything else only gravitationally. We find constraints given current cosmic microwave background (CMB) and baryon acoustic oscillation (BAO) data, with and without an H 0 prior, and with and without enforcing a big bang nucleosynthesis consistent helium abundance. We find that, at low dark photon temperature, one can have consistency with BAO and CMB data, with a fraction of dark matter that is ADM (f adm ) as large as ~0.1 . Such a large f adm leads to a suppression of density fluctuations today on scales below about 60 Mpc that may be of relevance to the σ 8 tension. Our work motivates calculation of nonlinear corrections to matter power spectrum predictions in the ADM model. Here, we forecast parameter constraints to come from future ground-based CMB surveys, and find that if ADM is indeed the cause of the σ 8 tension, the influence of the ADM, primarily on CMB lensing, will likely be detectable at high significance.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

AToM: Advanced Tokamak Modeling Environment

Stability constraints play an important role in integrated modeling. The global plasma stability and macro-instabilities such as internal kink modes, neoclassical tearing modes, and edge localized perturbations limit the plasma performance and can result in large-scale transient events and plasma disruptions. These macro-instabilities can be also beneficial to the plasma performance. For example, the peeling models often lead to the edge localized modes (ELMs). However, the formation of stationary edge harmonic oscillations (EHOs) due to nonlinear interaction of peeling modes can result in a transition to the Quiescent H-mode (QH-mode) without ELMs. The profiles predicted with transport models need to be refined using the MHD stability constraints. This is important especially for transient and nonlinear stages of discharges such as ramp-up access to hybrid and steady state operation, and L- to H-mode transition as well as transition to the QH-mode. Several stability and MHD codes are already included in the OMFIT framework. These codes include BALOO, ELITE, GATO, M3DC1, MARS, and NIMROD. However, the verifications of stability constraints are currently mostly excluded from transport modeling workflows. The only exception is the EPED module which includes the ELITE predictions to limit the pedestal height. Here, we improved the OMFIT workflow to include stability calculations. Including the stability conditions to the integrated modeling workflow improved the robustness of the predictive modeling discharges. Being implemented in the workflow, the stability conditions can be also used for the stability analysis of experimental data. This improved the physics understanding of various discharge scenarios and can be used in experiment planning.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1986

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.↗

DOE OSTI2022

A Julia Framework for Graph-Structured Nonlinear Optimization

Graph theory provides a convenient framework for modeling and solving structured optimization problems. Under this framework, the modeler can arrange/assemble the components of an optimization model (variables, constraints, objective functions, and data) within nodes and edges of a graph, and this representation can be used to visualize, manipulate, and solve the problem. In this work, we present a Julia framework for modeling and solving graph-structured nonlinear optimization problems. Our framework integrates the modeling package Plasmo.jl (which facilitates the construction and manipulation of graph models) and the nonlinear optimization solver MadNLP.jl (which provides capabilities for exploiting graph structures to accelerate solution). We illustrate with a simple example how model construction and manipulation can be performed in an intuitive manner using Plasmo.jl and how the model structure can be exploited by MadNLP.jl. We also demonstrate the scalability of the framework by targeting a large-scale, stochastic gas network problem that contains over 1.7 million variables.

Cole, David↗

NASA NTRS1990

Numerical modeling of the atmosphere with an isentropic vertical coordinate

A theta-coordinate model simulating the nonlinear evolution of a baroclinic wave is presented. In the model, vertical discretization maintains important integral constraints such as conservation of the angular momentum and total energy. A massless-layer approach is used in the treatment of the intersections of coordinate surfaces with the lower boundary. This formally eliminates the intersection problem, but raises other computational problems. Horizontal discretization of the continuity and momentum equations in the model are designed to overcome these problems. Selected results from a 10-day integration with the 25-layer, beta-plane version of the model are presented. It is concluded that the model can simulate the nonlinear evolution of a baroclinic wave and associated dynamical processes without major computational difficulties.

Hsu, Yueh-Jiuan G.↗

DOE OSTI2025

Computational Algorithms for Unit Commitment with AC Power Flows (Final Report)

Security-constrained unit commitment (SCUC) is a key component in power system operations. When AC power flow constraints are considered in the SCUC model (AC-SCUC), the problem becomes extremely difficult due to its discrete and non-convex nature, as described in “Grid Optimization Competition Challenge 3 Problem Formulation (GOCC)”. There are four main challenges: (i) Discrete decisions regarding unit online/offline status and start-up/shut-down procedures for every single unit. The number of discrete decision variables increases considerably when a system integrates multiple generators; (ii) Configuration-based combined-cycle formulations, and multi-commodity models that include ramping products, spin/non-spin products, and regulation up/down products. The combined-cycle units introduce additional discrete decision variables and auxiliary service products further complicate the model by connecting multi-commodity products’ continuous and discrete variables; (iii) SCUC models with AC power flow constraints are far more complex due to massive bilinear terms in the large-scale nonlinear power balance equations. The nonlinear power balance equations are further complicated by the discrete step control variables of shunts; (iv) N − 1 contingency analysis. The size of the model increases linearly with the number of contingencies considered, greatly increasing the size of the optimization model. Accordingly, there is an emergent need to develop a robust algorithm capable of deriving a high-quality solution in a short time and passing through contingency tests simultaneously. In this project, we explore innovative techniques to address this challenging problem by integrating advanced polyhedral theory, approximation methods, relaxation strategies, decomposition techniques, and parallel computing. Each technique approaches the problem from a different perspective, leveraging its specific strengths to tackle distinct challenges. Each individual method has demonstrated its effectiveness in the PI’s previous research. Their integration is expected to significantly reduce the computational time required to solve the proposed complex problem. Successful completion of this project has the potential to transform the industry by enhancing optimization solvers capable of handling large-scale day-ahead energy market clearing models within strict time constraints, while incorporating AC power flow constraints. This advancement will lead to reduced overall generation costs and, consequently, increased social welfare.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

NASA NTRS1974

Computer-aided design of control systems to meet many requirements

A method is described for using nonlinear programing in the computer-aided design of airplane control systems. It is assumed that the quality of such systems depends on many criteria. These criteria are included in the constraints vector (instead of attempting to combine them into a single scalar criterion, as is usually done), and the design proceeds through a sequence of nonlinear programing solutions in which the designer varies the specification of sets of requirements levels. The method is applied to design of a lateral stability augmentation system (SAS) for a fighter airplane, in which the requirements vector is chosen from the official handling qualities specifications. Results are shown for several simple SAS configurations designed to obtain desirable handling qualities over all design flight conditions with minimum feedback gains. The choice of the final design for each case is not unique but depends on the designer's decision as to which achievable set of requirements levels represents the best for that system. Results indicate that it may be possible to design constant parameter SAS which can satisfy the most stringent handling qualities requirements for fighter airplanes in all flight conditions. The role of the designer as a decision maker, interacting with the computer program, is discussed. Advantages of this type of designer-computer interaction are emphasized. Desirable extensions of the method are indicated.

Schy, A. A.↗

DOE OSTI2024

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗

NASA NTRS1975

Passive three-axis stabilization of the Long Duration Exposure Facility

This paper presents an analysis of the attitude dynamics of the Long Duration Exposure Facility (LDEF). LDEF is a large cylindrical gravity gradient stabilized earth satellite which is planned to be delivered to a 270-n mi circular orbit by the space shuttle. The fundamental linear stability, capture requirements, and pitch bias constraints generated by the Garber instability are discussed. Numerical simulations, based on the full nonlinear equations for the coupled orbital and attitude motion of the vehicle and the viscous magnetic damper, show stable behavior of the spacecraft and a damping time constant of 30 to 70 orbits.

Huckins, E. K., III↗

NASA NTRS1981

A four dimensional variational analysis experiment

A demonstration that the four dimensional variational analysis method using the governing equations as exact constraints can be successfully employed for a perfect model and for a simple, but nonlinear, system is presented. The method is stable in an assimilation cycle. The method reconstructs the unobservable variables; in the case demonstrated, no velocity data was observed. The analysis errors are much smaller than the observing system errors.

Hoffman, R.↗

NASA NTRS1983

A new algorithm for constrained nonlinear least-squares problems, part 1

A Gauss-Newton algorithm is presented for solving nonlinear least squares problems. The problem statement may include simple bounds or more general constraints on the unknowns. The algorithm uses a trust region that allows the objective function to increase with logic for retreating to best values. The computations for the linear problem are done using a least squares system solver that allows for simple bounds and linear constraints. The trust region limits are defined by a box around the current point. In its current form the algorithm is effective only for problems with small residuals, linear constraints and dense Jacobian matrices. Results on a set of test problems are encouraging.

Hanson, R. J.↗

NASA NTRS1991

A synthetic environment for visualization and planning of orbital maneuvers

An interactive proximity operations planning system, which allows on-site planning of fuel-efficient, multi-burn maneuvers in a potential multi-space-craft environment has been developed and tested. This display system most directly assists planning by providing visual feedback in a synthetic virtual space that aids visualization of trajectories and their constraints. Its most significant features include (1) an 'inverse dynamics' algorithm that removes control nonlinearities facing the operator and (2) a stack-oriented action-editor that reduces the order of control and creates, through a 'geometric spreadsheet,' the illusion of an inertially stable environment. This synthetic environment provides the user with control of relevant static and dynamic properties of way-points during small orbital changes allowing independent solutions to otherwise coupled problems of orbital maneuvering.

Ellis, Stephen R.↗

NASA NTRS1999

Design Optimization of Stiffened Panels with Postbuckling Constraints

The funding provided by the grant is used to complete the final stages of the development of a geometrically nonlinear analysis and design capability for the static response of compressively loaded prismatic plate structures. The analysis is based on the nonlinear finite strip method and is applicable for structures, such as stiffened panels or box columns, that can be modeled as assemblages of finite length plate strips. In an effort to reduce the computational cost of the nonlinear finite strip method, thus making it suitable for use in the design optimization environment, reduced basis techniques as described in various references by Noor are used in conjunction with the finite strip method. In addition, an efficient scheme for tracing the nonlinear equilibrium paths through highly nonlinear response curves was implemented. The new scheme, which is referred to as the normal flow algorithm, is based on homotopy methods, and is capable of negotiating highly nonlinear limit point instabilities with a reduced computational cost compared to the popular Ricks/Wempner and Chrisfield algorithms.

Guerdal, Zafer↗

DOE OSTI2022

Modeling design and control problems involving neural network surrogates

Here, we consider nonlinear optimization problems that involve surrogate models represented by neural networks. We demonstrate first how to directly embed neural network evaluation into optimization models, highlight a difficulty with this approach that can prevent convergence, and then characterize stationarity of such models. We then present two alternative formulations of these problems in the specific case of feedforward neural networks with ReLU activation: as a mixed-integer optimization problem and as a mathematical program with complementarity constraints. For the latter formulation we prove that stationarity at a point for this problem corresponds to stationarity of the embedded formulation. Each of these formulations may be solved with state-of-the-art optimization methods, and we show how to obtain good initial feasible solutions for these methods. We compare our formulations on three practical applications arising in the design and control of combustion engines, in the generation of adversarial attacks on classifier networks, and in the determination of optimal flows in an oil well network.

97 MATHEMATICS AND COMPUTING↗