Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear programming.”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Transonic Flutter Suppression Control Law Design, Analysis and Wind-Tunnel Results

The benchmark active controls technology and wind tunnel test program at NASA Langley Research Center was started with the objective to investigate the nonlinear, unsteady aerodynamics and active flutter suppression of wings in transonic flow. The paper will present the flutter suppression control law design process, numerical nonlinear simulation and wind tunnel test results for the NACA 0012 benchmark active control wing model. The flutter suppression control law design processes using classical, and minimax techniques are described. A unified general formulation and solution for the minimax approach, based on the steady state differential game theory is presented. Design considerations for improving the control law robustness and digital implementation are outlined. It was shown that simple control laws when properly designed based on physical principles, can suppress flutter with limited control power even in the presence of transonic shocks and flow separation. In wind tunnel tests in air and heavy gas medium, the closed-loop flutter dynamic pressure was increased to the tunnel upper limit of 200 psf. The control law robustness and performance predictions were verified in highly nonlinear flow conditions, gain and phase perturbations, and spoiler deployment. A non-design plunge instability condition was also successfully suppressed.

Mukhopadhyay, Vivek↗
NASA NTRS1999

Transonic Flutter Suppression Control Law Design Using Classical and Optimal Techniques with Wind-Tunnel Results

The benchmark active controls technology and wind tunnel test program at NASA Langley Research Center was started with the objective to investigate the nonlinear, unsteady aerodynamics and active flutter suppression of wings in transonic flow. The paper will present the flutter suppression control law design process, numerical nonlinear simulation and wind tunnel test results for the NACA 0012 benchmark active control wing model. The flutter suppression control law design processes using (1) classical, (2) linear quadratic Gaussian (LQG), and (3) minimax techniques are described. A unified general formulation and solution for the LQG and minimax approaches, based on the steady state differential game theory is presented. Design considerations for improving the control law robustness and digital implementation are outlined. It was shown that simple control laws when properly designed based on physical principles, can suppress flutter with limited control power even in the presence of transonic shocks and flow separation. In wind tunnel tests in air and heavy gas medium, the closed-loop flutter dynamic pressure was increased to the tunnel upper limit of 200 psf. The control law robustness and performance predictions were verified in highly nonlinear flow conditions, gain and phase perturbations, and spoiler deployment. A non-design plunge instability condition was also successfully suppressed.

Mukhopadhyay, Vivek↗
NASA NTRS1982

A doublet lattice method for the determination of rotor induced empennage vibration airloads. Analysis description and program documentation

An efficient state-of-the-art method was developed to determine the unteady vibratory airloads produced by the interaction of the main rotor wake with a helicopter empennage. This method was incorporated into a computer program, Rotor Induced Empennage Vibration Analysis (RIEVA). The program requires the main rotor wake position and the strength of the vortices locates near the empennage surfaces. A nonlinear lifting surface analysis is utilized to predict the aerodynamic loads on the empennage surfaces in the presence of these concentrated vortices. The analysis was formulated to include all pertinent effects such as suction of the interacting vortices and the shed vorticity behind the empennage surfaces. The analysis employs a time domain solution. The output of the program consists of chordwise and spanwise airload distributions on the empennage surfaces. The airload distributions are harmonically analyzed and formulated for input into the Coupled Rotor/Airframe Vibration Analysis.

Gangwani, S. T.↗
DOE OSTI2021

Iteration-based Linearized Distribution-level Locational Marginal Price for Three-phase Unbalanced Distribution Systems

Distributed energy resources (DERs) are rocking the utilities’ business landscape. It calls for competitive market environments that incentivize DERs to form maximum operating efficiency. Among proposed pricing schemes, distribution-level locational marginal price (DLMP) is effective in signaling the marginal generation cost differences driven by energy losses and network constraints. It can be derived from a distribution-level optimal power flow (OPF) framework, as it essentially presents the sensitivity of optimized generation cost towards incremental loads. However, due to the high resistance-to-inductance ratio and unbalanced characteristics of distribution networks, computational affordable DLMPs are highly challenged. This article provides a linear-approximated DLMP that can be solved efficiently and generalized to account for reactive power flow, three-phase unbalanced loads and meshed network structure. The successive linear programming technique is introduced to enhance the model accuracy. Case studies on an IEEE 123-Bus system validate its accuracy against a nonlinear benchmark and capability in offering proper incentives.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1991

Supercomputer optimizations for stochastic optimal control applications

Supercomputer optimizations for a computational method of solving stochastic, multibody, dynamic programming problems are presented. The computational method is valid for a general class of optimal control problems that are nonlinear, multibody dynamical systems, perturbed by general Markov noise in continuous time, i.e., nonsmooth Gaussian as well as jump Poisson random white noise. Optimization techniques for vector multiprocessors or vectorizing supercomputers include advanced data structures, loop restructuring, loop collapsing, blocking, and compiler directives. These advanced computing techniques and superconducting hardware help alleviate Bellman's curse of dimensionality in dynamic programming computations, by permitting the solution of large multibody problems. Possible applications include lumped flight dynamics models for uncertain environments, such as large scale and background random aerospace fluctuations.

Chung, Siu-Leung↗
NASA NTRS1988

Nonlinear dynamics in cardiac conduction

Electrical conduction in the heart shows many phenomena familiar from nonlinear dynamics. Among these phenomena are multiple basins of attraction, phase locking, and perhaps period-doubling bifurcations and chaos. We describe a simple cellular-automation model of electrical conduction which simulates normal conduction patterns in the heart as well as a wide range of disturbances of heart rhythm. In addition, we review the application of percolation theory to the analysis of the development of complex, self-sustaining conduction patterns.

NASA Discipline Cardiopulmonary↗
NASA NTRS1981

The Numerical Studies Program for the Atmospheric General Circulation Experiment (AGCE) for Spacelab Flights

The atmospheric general circulation experiment (AGCE) numerical design for Spacelab flights was studied. A spherical baroclinic flow experiment which models the large scale circulations of the Earth's atmosphere was proposed. Gravity is simulated by a radial dielectric body force. The major objective of the AGCE is to study nonlinear baroclinic wave flows in spherical geometry. Numerical models must be developed which accurately predict the basic axisymmetric states and the stability of nonlinear baroclinic wave flows. A three dimensional, fully nonlinear, numerical model and the AGCE based on the complete set of equations is required. Progress in the AGCE numerical design studies program is reported.

Fowlis, W. W.↗
NASA NTRS1992

Focal plane arrays for submillimeter waves using two-dimensional electron gas elements: A grant under the Innovative Research Program

This final report describes a three-year research effort, aimed at developing new types of THz low noise receivers, based on bulk effect ('hot electron') nonlinearities in the Two-Dimensional Electron Gas (2DEG) Medium, and the inclusion of such receivers in focal plane arrays. 2DEG hot electron mixers have been demonstrated at 35 and 94 GHz with three orders of magnitude wider bandwidth than previous hot electron mixers, which use bulk InSb. The 2DEG mixers employ a new mode of operation, which was invented during this program. Only moderate cooling is required for this mode, to temperatures in the range 20-77 K. Based on the results of this research, it is now possible to design a hot electron mixer focal plane array for the THz range, which is anticipated to have a DSB receiver noise temperature of 500-1000K. In our work on this grant, we have found similar results the the Cronin group (resident at the University of Bath, UK). Neither group has so far demonstrated heterodyne detection in this mode, however. We discovered and explored some new effects in the magnetic field mode, and these are described in the report. In particular, detection of 94 GHz and 238 GHz, respectively, by a new effect, 'Shubnikov de Haas detection', was found to be considerably stronger in our materials than the cyclotron resonance detection. All experiments utilized devices with an active 2DEG region of size of the order of 10-40 micrometers long, and 20-200 micrometers wide, formed at the heterojunction between AlGaAs and GaAs. All device fabrication was performed in-house. The materials for the devices were also grown in-house, utilizing OMCVD (Organo Metallic Chemical Vapor Deposition). In the course of this grant, we developed new techniques for growing AlGaAs/GaAs with mobilities equalling the highest values published by any laboratory. We believe that the field of hot electron mixers and detectors will grow substantially in importance in the next few years, partly as a result of the opportunity given us through this grant, which represents the major effort in the US so far. We note, however, that parallel research on hot electron mixers in thin film superconductors in Russia, and recently in Sweden, have demonstrated mixing up to 1 THz, with the potential for low-noise receivers for frequencies up to many THz. The three groups recently assessed the relative adtantages of 2DEG and superconducting film mixers in a joint paper (Kollberg et al., 1992; see Appendix II).

Yngvesson, K. Sigfrid↗
DOE OSTI2016

Solution of monotone complementarity and general convex programming problems using a modified potential reduction interior point method

We present a homogeneous algorithm equipped with a modified potential function for the monotone complementarity problem. We show that this potential function is reduced by at least a constant amount if a scaled Lipschitz condition (SLC) is satisfied. A practical algorithm based on this potential function is implemented in a software package named iOptimize. The implementation in iOptimize maintains global linear and polynomial time convergence properties, while achieving practical performance. It either successfully solves the problem, or concludes that the SLC is not satisfied. When compared with the mature software package MOSEK (barrier solver version 6.0.0.106), iOptimize solves convex quadratic programming problems, convex quadratically constrained quadratic programming problems, and general convex programming problems in fewer iterations. Moreover, several problems for which MOSEK fails are solved to optimality. In addition, we also find that iOptimize detects infeasibility more reliably than the general nonlinear solvers Ipopt (version 3.9.2) and Knitro (version 8.0).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Vibration computer programs E13101, E13102, E13104, and E13112 and application to the NERVA program. Project 187: Methodology documentation

Programs for the analyses of the free or forced, undamped vibrations of one or two elastically-coupled lumped parameter teams are presented. Bearing nonlinearities, casing and rotor distributed mass and elasticity, rotor imbalance, forcing functions, gyroscopic moments, rotary inertia, and shear and flexural deformations are all included in the system dynamics analysis. All bearings have nonlinear load displacement characteristics, the solution is achieved by iteration. Rotor imbalances allowed by such considerations as pilot tolerances and runouts as well as bearing clearances (allowing concail or cylindrical whirl) determine the forcing function magnitudes. The computer programs first obtain a solution wherein the bearings are treated as linear springs of given spring rates. Then, based upon the computed bearing reactions, new spring rates are predicted and another solution of the modified system is made. The iteration is continued until the changes to bearing spring rates and bearing reactions become negligibly small.

Mironenko, G.↗
NASA NTRS2020

A Deformed Geometry Synthesis Technique for Determining Stacking and Cryogenically Induced Preloads for the Space Launch System

The Space Launch System (SLS) stacking and Core Stage (CS) fueling induce significant preloads that contribute to the liftoff pad separation “twang”. To accurately capture this, an approach is required that can replicate the physics of all SLS physical stacking steps, CS cryogenic shrinkage, associated geometric nonlinearities, and the transient behavior and decay of the preloads with changing boundary conditions as the vehicle separates from the pad. The Deformed Geometry Synthesis (DGS) approach presented here satisfies the above requirements. DGS determines induced preloads by modeling components in their deformed geometry states and then enforcing compatibility by closing the resulting “deadbands”. DGS seamlessly integrates into the multibody modal synthesis framework and does not require the use of artificial external loads to enforce preloads or post-processing steps to remove their influence. Since DGS iterates to solve for the deformed state inclusive of geometric nonlinearities, running linearized parametrics to exercise different potential orientations of ball jointed struts that connect the CS to Boosters for cryo-shrinkage analyses is entirely avoided. Relative to the transient behavior and decay of stacking and cryo-induced preloads with SLS liftoff pad separation, this is an area of considerable interest to the SLS program. To capture this in the most accurate way possible, DGS algorithms are designed to work with Henkel-Mar nonlinear pad separation algorithms which operate on the separating longitudinal and lateral degrees of freedom (DoFs) between the vehicle and the pad. As the separating DoFs release, in whatever manner as dictated by the interface geometries, interface loads and interface flexibilities as well as the external loading on the vehicle, the subject preloads generate a complex twang/decay time-trace as dictated by the physics of the problem. This paper presents DGS numerical verification against the closed-form solution for Timoshenko’s 3 ball-jointed strut preload problem. This problem is then extended by the authors to the geometric nonlinear case where DGS is compared to the Newton-Raphson solution of the nonlinear equations. Next, DGS is utilized to solve the SLS stacking and cryogenic shrinkage coupled loads analyses. Finally, Henkel-Mar pad separation simulations are executed that isolate the impact of the induced preloads’ twang and decay characteristics.

Space Launch System↗
DOE OSTI2018

Model Identification of Inverter Nonlinear Control Dynamics

To comprehensively predict in silico potential disruptive dynamic performance of converter-based generation and develop solutions, the application of electromagnetic transient programs is desirable. However, these accurate simulations highly depend on inverter models of actual field implementations, which are often unavailable because they are proprietary. To address this problem, we propose a time-domain model identification method that can be used to determine the control models of inverters in continuous time from multiple experimental studies. The focus is on identifying the dynamic representation of controls because they dominate the dynamic response of renewable assets. To illustrate the applicability of the technique, we consider a subclass of inverter controllers for autonomous grid applications that rely on nonlinear oscillators.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1989

The effects of nonlinear loading upon the Space Station Freedom 20 kHz power system

The Space Station Freedom power distribution system, which consists of dual redundant 20-kHz, 440-V RMS, single-phase power systems, is discussed. The effect of a typical space station nonlinear load on the measurement of RMS current and voltage at various points in the space station power system has been investigated using the Electromagnetic Transients Program (EMTP). The load current distortion at the user interface, its effect on the distribution system, and its relationship to power factor have been studied. Modeling results are compared to test data. The differences under nonlinear loading are evaluated and presented as a measure of distribution voltage distortion and current measurement accuracy.

Leskovich, R. Thomas↗
DOE OSTI2017

Physics-based control-oriented modeling of the current density profile evolution in NSTX-Upgrade

Active control of the toroidal current density profile is among those plasma control milestones that the National Spherical Torus eXperiment-Upgrade (NSTX-U) program must achieve to realize its next-step operational goals. Motivated by the coupled, nonlinear, multivariable, distributed-parameter plasma dynamics, the first step towards control design is the development of a physics-based, control-oriented model for the current profile evolution in response to non-inductive current drives and heating systems. The evolution of the toroidal current density profile is closely related to the evolution of the poloidal magnetic flux profile, whose dynamics is modeled by a nonlinear partial differential equation (PDE) referred to as the magnetic-flux diffusion equation (MDE). The proposed control-oriented model predicts the spatial-temporal evolution of the current density profile by combining the nonlinear MDE with physics-based correlations obtained at NSTX-U for the electron density, electron temperature, and non-inductive current drives (neutral beams). The resulting first-principles-driven, control-oriented model is tailored for NSTX-U based on predictions by the time-dependent transport code TRANSP. Furthermore, main objectives and possible challenges associated with the use of the developed model for the design of both feedforward and feedback controllers are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Investigation of a New Monte Carlo Method for the Transitional Gas Flow

The Direct Simulation Monte Carlo method (DSMC) is well developed for rarefied gas flow in transition flow regime when 0.01<Kn<1. However, such a simulation for a complex 3D vacuum system is still a challenging task because of the huge demand on the memory and long computational time. On the other hand, if Kn>10, the gas flow is free molecular and can be simulated by the Test Particle Monte Carlo method (TPMC) without any problem even for a complex 3D vacuum system. In this paper we will investigate the approach to extend the TPMC to the transition flow regime by considering the collision between gas molecules as an interaction between a probe molecule and the gas background. Recently this collision mechanism has been implemented into ProVac3D, a new TPMC simulation program developed by Karlsruhe Institute of Technology (KIT). The preliminary simulation result shows a correct nonlinear increasing of the gas flow. However, there is still a quantitative discrepancy with the experimental data, which means further improvement is needed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Adaptive PID Gain Scheduling Control for Hydropower Turbine Using Neural CDE and Stochastic Distribution Shaping

This paper introduces a gain-scheduling PID controller design strategy for hydroturbine frequency control mode. This scheme first uses real data to learn the nonlinear dynamics of the hydroturbine using neural controlled differential equations and then perturbs the obtained nonlinear system at different equilibrium points, based on which a static output feedback adaptive dynamic programming algorithm is then used to optimize the PID gains for each equilibrium point. Moreover, a continuous-time version of stochastic distribution control is proposed to further fine-tune the optimized PID gains. Finally, the controller is obtained by implementing linear interpolation between the optimized PID control gains. The simulation results show that the proposed gain-scheduling PID controller can control a larger range of operation points compared with the given fixed PID controller and the baseline method. Compared with the given fixed PID controller, the proposed gain-scheduling PID controller can regulate hydroturbine frequency against disturbances induced by power-load variation with over 50% less overshoot for some operation points.

13 HYDRO ENERGY↗
DOE OSTI2021

Convex Optimization of Integrated Power-Gas Energy Flow Model With Applications to Probabilistic Energy Flow

Energy flow calculation is a fundamental problem of the integrated power and gas system (IPGS) operation and planning. However, the nonlinear gas flow model introduces major challenges to the energy flow calculation. In this paper, we propose a tractably convex optimization model to solve the energy flow problem in IPGSs. It is demonstrated that the proposed optimization model has the same optimal solution as the original nonlinear steady energy flow model. Also, piecewise linearization is adopted to tightly linearize the nonlinear objective function of the model, which transforms the formulated convex optimization into a linear program one. Thus, the computation complexity of the proposed energy flow model is significantly reduced as compared with the existing methods. In addition, the proposed model can be extended to probabilistic energy flow estimation. Extensive case studies are conducted to validate the effectiveness of the proposed energy flow model using three IPGSs.

42 ENGINEERING↗
DOE OSTI2022

Mitigating the Impacts of Uncertain Geomagnetic Disturbances on Electric Grids: A Distributionally Robust Optimization Approach

Severe geomagnetic disturbances (GMDs) increase the magnitude of the electric field on the Earth’s surface (E-field) and drive geomagnetically-induced currents (GICs) along the transmission lines in electric grids. These additional currents can pose severe risks, such as current distortions, transformer saturation and increased reactive power losses, each of which can lead to system unreliability. Today several mitigation actions (e.g., changing grid topology) exist that can reduce the harmful GIC effects on the grids. Making such decisions can be challenging, however, because the magnitude and direction of the E-field are uncertain and non-stationary. In this paper, we model uncertain E-fields using the distributionally robust optimization (DRO) approach that determines optimal transmission grid operations such that the worst-case expectation of the system cost is minimized. We also capture the effect of GICs on the nonlinear AC power flow equations. For solution approaches, we develop an accelerated column-and-constraint generation (CCG) algorithm by exploiting a special structure of the support set of uncertain parameters representing the E-field. Extensive numerical experiments based on “epri-21” and “uiuc-150” systems, designed for GMD studies, demonstrate (i) the computational performance of the accelerated CCG algorithm, (ii) the superior performance of distributionally robust grid operations that satisfy nonlinear, nonconvex AC power flow equations and GIC constraints, in comparison with standard stochastic programming-based methods during the out-of-sample testing.

42 ENGINEERING↗