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 1,081 records · Page 60

Electro-optical flight deflection measurement system

This paper describes an electrooptical flight-deflection measurement system (FDMS) developed for use on highly-maneuverable-aircraft-technology (HiMAT) RPV flight research program. The FDMS provides in-flight measurements of the aircraft structural deflections and magnetic-tape recordings for automated data processing. The capabilities and limitations, requirements for installation on an aircraft, analytical considerations, and typical flight data acquired from the HiMAT research program are examined. The flight data indicate that the background light is the major obstacle to acquiring high-quality data and that the relationship between the target displacement and the output of the FDMS is nonlinear; however, the nonlinear effects can be minimized with judicious planning of the installation of the FDMS on the aircraft. Excellent flight deflection data were obtained from both the HiMAT and X-29 A flight research program with very little data lost as a result of encounters with severe background light.

Deangelis, V. Michael↗

Implicit shock tracking for unsteady flows by the method of lines

Here, a recently developed high-order implicit shock tracking (HOIST) framework for resolving discontinuous solutions of inviscid, steady conservation laws is extended to the unsteady case. Central to the framework is an optimization problem which simultaneously computes a discontinuity-aligned mesh and the corresponding high-order approximation to the flow, which provides nonlinear stabilization and a high-order approximation to the solution. This work extends the implicit shock tracking framework to the case of unsteady conservation laws using a method of lines discretization via a diagonally implicit Runge-Kutta method by “solving a steady problem at each timestep”. We formulate and solve an optimization problem that produces a feature-aligned mesh and solution at each Runge-Kutta stage of each timestep, and advance this solution in time by standard Runge-Kutta update formulas. A Rankine-Hugoniot based prediction of the shock location together with a high-order, untangling mesh smoothing procedure provides a high-quality initial guess for the optimization problem at each time, which results in rapid convergence of the sequential quadratic programing (SQP) optimization solver. This method is shown to deliver highly accurate solutions on coarse, high-order discretizations without nonlinear stabilization and recover the design accuracy of the Runge-Kutta scheme. We demonstrate this framework on a series of inviscid, unsteady conservation laws in both one- and two- dimensions. We also verify that our method is able to recover the design order of accuracy of our time integrator in the presence of a strong discontinuity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Identification of Linear and Nonlinear Aerodynamic Impulse Responses Using Digital Filter Techniques

This paper discusses the mathematical existence and the numerically-correct identification of linear and nonlinear aerodynamic impulse response functions. Differences between continuous-time and discrete-time system theories, which permit the identification and efficient use of these functions, will be detailed. Important input/output definitions and the concept of linear and nonlinear systems with memory will also be discussed. It will be shown that indicial (step or steady) responses (such as Wagner's function), forced harmonic responses (such as Tbeodorsen's function or those from doublet lattice theory), and responses to random inputs (such as gusts) can all be obtained from an aerodynamic impulse response function. This paper establishes the aerodynamic impulse response function as the most fundamental, and, therefore, the most computationally efficient, aerodynamic function that can be extracted from any given discrete-time, aerodynamic system. The results presented in this paper help to unify the understanding of classical two-dimensional continuous-time theories with modem three-dimensional, discrete-time theories. First, the method is applied to the nonlinear viscous Burger's equation as an example. Next the method is applied to a three-dimensional aeroelastic model using the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) code and then to a two-dimensional model using the CFL3D Navier-Stokes code. Comparisons of accuracy and computational cost savings are presented. Because of its mathematical generality, an important attribute of this methodology is that it is applicable to a wide range of nonlinear, discrete-time problems.

Silva, Walter A.↗

Identification of Linear and Nonlinear Aerodynamic Impulse Responses Using Digital Filter Techniques

This paper discusses the mathematical existence and the numerically-correct identification of linear and nonlinear aerodynamic impulse response functions. Differences between continuous-time and discrete-time system theories, which permit the identification and efficient use of these functions, will be detailed. Important input/output definitions and the concept of linear and nonlinear systems with memory will also be discussed. It will be shown that indicial (step or steady) responses (such as Wagner's function), forced harmonic responses (such as Theodorsen's function or those from doublet lattice theory), and responses to random inputs (such as gusts) can all be obtained from an aerodynamic impulse response function. This paper establishes the aerodynamic impulse response function as the most fundamental, and, therefore, the most computationally efficient, aerodynamic function that can be extracted from any given discrete-time, aerodynamic system. The results presented in this paper help to unify the understanding of classical two-dimensional continuous-time theories with modern three-dimensional, discrete-time theories. First, the method is applied to the nonlinear viscous Burger's equation as an example. Next the method is applied to a three-dimensional aeroelastic model using the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) code and then to a two-dimensional model using the CFL3D Navier-Stokes code. Comparisons of accuracy and computational cost savings are presented. Because of its mathematical generality, an important attribute of this methodology is that it is applicable to a wide range of nonlinear, discrete-time problems.

Silva, Walter A.↗

Large displacements and stability analysis of nonlinear propeller structures

The use of linear rigid formats in COSMIC NASTRAN without DMAP procedures for the analysis of nonlinear propeller structures is described. Approaches for updating geometry and applying follower forces for incremental loading are demonstrated. The COSMIC NASTRAN rigid formats and other independent finite element programs are compared. The comparisons include results from the four approaches for updating the geometry using RIGID FORMAT 1, RIGID FORMATS 4 and 13, MARC and MSC/NASTRAN. It is shown that user friendly updating approaches can be used to predict the large displacements and instability of these nonlinear structures. The approaches are easily implemented by the user and predict conservative results.

Aiello, R. A.↗

ALESQP: An Augmented Lagrangian Equality-Constrained SQP Method for Optimization with General Constraints

Here we present a new algorithm for infinite-dimensional optimization with general constraints, called ALESQP. In short, ALESQP is an augmented Lagrangian method that penalizes inequality constraints and solves equality-constrained nonlinear optimization subproblems at every iteration. The subproblems are solved using a matrix-free trust-region sequential quadratic programming (SQP) method that takes advantage of iterative, i.e., inexact linear solvers, and is suitable for large-scale applications. A key feature of ALESQP is a constraint decomposition strategy that allows it to exploit problem-specific variable scalings and inner products. We analyze convergence of ALESQP under different assumptions. We show that strong accumulation points are stationary. Consequently, in finite dimensions ALESQP converges to a stationary point. In infinite dimensions we establish that weak accumulation points are feasible in many practical situations. Under additional assumptions we show that weak accumulation points are stationary. We present several infinite-dimensional examples where ALESQP shows remarkable discretization-independent performance in all of its iterative components, requiring a modest number of iterations to meet constraint tolerances at the level of machine precision. Also, we demonstrate a fully matrix-free solution of an infinite-dimensional problem with nonlinear inequality constraints.

97 MATHEMATICS AND COMPUTING↗

NASA Tech Briefs, February 2004

Topics include: Simulation Testing of Embedded Flight Software; Improved Indentation Test for Measuring Nonlinear Elasticity; Ultraviolet-Absorption Spectroscopic Biofilm Monitor; Electronic Tongue for Quantitation of Contaminants in Water; Radar for Measuring Soil Moisture Under Vegetation; Modular Wireless Data-Acquisition and Control System; Microwave System for Detecting Ice on Aircraft; Routing Algorithm Exploits Spatial Relations; Two-Finger EKG Method of Detecting Evasive Responses; Updated System-Availability and Resource-Allocation Program; Routines for Computing Pressure Drops in Venturis; Software for Fault-Tolerant Matrix Multiplication; Reproducible Growth of High-Quality Cubic-SiC Layers; Nonlinear Thermoelastic Model for SMAs and SMA Hybrid Composites; Liquid-Crystal Thermosets, a New Generation of High-Performance Liquid-Crystal Polymers; Formulations for Stronger Solid Oxide Fuel-Cell Electrolytes; Simulation of Hazards and Poses for a Rocker-Bogie Rover; Autonomous Formation Flight; Expandable Purge Chambers Would Protect Cryogenic Fittings; Wavy-Planform Helicopter Blades Make Less Noise; Miniature Robotic Spacecraft for Inspecting Other Spacecraft; Miniature Ring-Shaped Peristaltic Pump; Compact Plasma Accelerator; Improved Electrohydraulic Linear Actuators; A Software Architecture for Semiautonomous Robot Control; Fabrication of Channels for Nanobiotechnological Devices; Improved Thin, Flexible Heat Pipes; Miniature Radioisotope Thermoelectric Power Cubes; Permanent Sequestration of Emitted Gases in the Form of Clathrate Hydrates; Electrochemical, H2O2-Boosted Catalytic Oxidation System; Electrokinetic In Situ Treatment of Metal-Contaminated Soil; Pumping Liquid Oxygen by Use of Pulsed Magnetic Fields; Magnetocaloric Pumping of Liquid Oxygen; Tailoring Ion-Thruster Grid Apertures for Greater Efficiency; and Lidar for Guidance of a Spacecraft or Exploratory Robot.

Source record↗

Variational Trajectory Optimization Tool Set: Technical description and user's manual

The algorithms that comprise the Variational Trajectory Optimization Tool Set (VTOTS) package are briefly described. The VTOTS is a software package for solving nonlinear constrained optimal control problems from a wide range of engineering and scientific disciplines. The VTOTS package was specifically designed to minimize the amount of user programming; in fact, for problems that may be expressed in terms of analytical functions, the user needs only to define the problem in terms of symbolic variables. This version of the VTOTS does not support tabular data; thus, problems must be expressed in terms of analytical functions. The VTOTS package consists of two methods for solving nonlinear optimal control problems: a time-domain finite-element algorithm and a multiple shooting algorithm. These two algorithms, under the VTOTS package, may be run independently or jointly. The finite-element algorithm generates approximate solutions, whereas the shooting algorithm provides a more accurate solution to the optimization problem. A user's manual, some examples with results, and a brief description of the individual subroutines are included.

Bless, Robert R.↗

Measured and predicted structural behavior of the HiMAT tailored composite wing

A series of load tests was conducted on the HiMAT tailored composite wing. Coupon tests were also run on a series of unbalanced laminates, including the ply configuration of the wing, the purpose of which was to compare the measured and predicted behavior of unbalanced laminates, including - in the case of the wing - a comparison between the behavior of the full scale structure and coupon tests. Both linear and nonlinear finite element (NASTRAN) analyses were carried out on the wing. Both linear and nonlinear point-stress analyses were performed on the coupons. All test articles were instrumented with strain gages, and wing deflections measured. The leading and trailing edges were found to have no effect on the response of the wing to applied loads. A decrease in the stiffness of the wing box was evident over the 27-test program. The measured load-strain behavior of the wing was found to be linear, in contrast to coupon tests of the same laminate, which were nonlinear. A linear NASTRAN analysis of the wing generally correlated more favorably with measurements than did a nonlinear analysis. An examination of the predicted deflections in the wing root region revealed an anomalous behavior of the structural model that cannot be explained. Both hysteresis and creep appear to be less significant in the wing tests than in the corresponding laminate coupon tests.

Nelson, Lawrence H.↗

Pyomo.DOE: An open-source package for model-based design of experiments in Python

Predictive mathematical models are a cornerstone of science and engineering. Yet selecting, calibrating, and validating said science-based models often remains an art in practice. Model-based design of experiments (MBDoE) provides a systematic framework to maximize information gain from experiments while minimizing time and resource costs. But MBDoE remains limited to niche application areas, in part because practitioners must integrate expertise in statistics, computational optimization, and modeling. To help reduce this barrier, we introduce Pyomo.DOE, an open-source package for MBDoE. Pyomo.DOE uses a nonlinear sensitivity analysis code k_aug to quickly approximate the Fisher information matrix and leverages a new stochastic programming abstraction. We demonstrate Pyomo.DOE with the first application of MBDoE to fixed-bed breakthrough experiments, which highlights the power of Pyomo.DOE to quantify the value of experimental modifications a priori for large-scale partial differential-algebraic equation (PDAE) models. Here we also provide a mathematical primer on MBDoE targeted at general chemical engineers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimal design of microtube recuperators for an indirect supercritical carbon dioxide recompression closed Brayton cycle

This paper presents a baseline design and optimization approach developed in Aspen Custom Modeler (ACM) for microtube shell-and-tube exchangers (MSTEs) used for high- and low-temperature recuperation in a 10 MWe indirect supercritical carbon dioxide (sCO 2 ) recompression closed Brayton cycle (RCBC). The MSTE-type recuperators are designed using one-dimensional models with thermal-hydraulic correlations appropriate for sCO 2 and properties models that capture considerable nonlinear changes in CO 2 properties near the critical and pseudo-critical points. Using the successive quadratic programming (SQP) algorithm in ACM, optimal recuperator designs are obtained for either custom or industry-standard microtubes considering constraints based on current advanced manufacturing techniques. The three decision variables are the number of tubes, tube pitch-to-diameter ratio, and tube diameter. Five different objective functions based on different key design measures are considered: minimization of total heat transfer area, heat exchanger volume, metal weight, thermal residence time, and maximization of compactness. Sensitivities studies indicate the constraint on the maximum number of tubes per shell does affect the number of parallel heat exchanger trains but not the tube selection, total number of tubes, tube length and other key design measures in the final optimal design when considering industry-standard tubes. In this study, the optimally designed high- and low-temperature recuperators have 47,000 3/32 inch tubes and 63,000 1/16 inch tubes, respectively. In addition, sensitivities to the design temperature approach and maximum allowable pressure drop are studied, since these specifications significantly impact the optimal design of the recuperators as well as the thermal efficiency and the economic performance of the entire sCO 2 Brayton cycle.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Distributed Load Restoration in Unbalanced Active Distribution Systems

We present that the resilience of modern power systems has been constantly threatened by various natural disasters and manmade attacks. Under those threats, fast and automatic restoration actions are critical for restoring impacted systems to normal operating conditions and preventing from additional disastrous consequences. Smart grid technologies, such as distributed energy resources (DERs) and microgrids, provide both opportunities and challenges for distribution system restoration. In this paper, distributed load restoration (DLR), in unbalanced active distribution networks is developed using the alternating direction method of multipliers (ADMM) algorithm. The step-by-step restoration actions are provided to schedule the output of DERs and distribution system control devices, as well as the power consumption from transmission system restoration. The nonlinearities from three-phase unbalanced power flow and distribution components modeling are relaxed into a convex quadratic programming model. Then, the problem is decomposed into subproblems for each node by applying ADMM-DLR; and solves through each agent by exchanging limited information with neighboring nodes in an iterative procedure. The developed models and algorithms are validated and demonstrated through testing of the IEEE 13-node and 123-node test feeders. Lastly, simulation results demonstrate the effectiveness of restoring unbalanced distribution system with the hybrid bottom-up and top-down restoration strategies.

14 SOLAR ENERGY↗

On visual-vestibular interaction

Experimental evidence is presented regarding visual vestibular interaction, and the results of three studies on the subject are briefly noted. An attempt to put together some of these observations with elementary notions of a visual vestibular interaction program is shown in the form of a flow chart representation of a possible model. This is a nonlinear model in which visual and vestibular influences are linearly weighted when they are in relative agreement but switch to the more believable one when they are in disagreement. A solution to the human space orientation problem is depicted by a schema for optimal subjective orientation based on several sensory modalities.

Young, L. R.↗

An improved method for the prediction of completely three-dimensional aerodynamic load distributions of configurations with leading edge vortex separation

The application of a higher-order subsonic potential flow panel method to the solution of three-dimensional flow about wing and wing-body combinations with leading-edge vortex separation is presented. The governing equations are the linear flow differential equation and nonlinear boundary conditions which require that the flow be parallel to the wing and body surfaces and that the free vortex sheet, springing from the leading and trailing edges, be aligned with the local flow and support no pressure jump. The vortex core is modeled as a simple line vortex which receives vorticity from the free sheet through a connecting sheet. The Kutta condition is imposed on all appropriate edges of the wing. This set of nonlinear equations is solved by an iterative procedure. The Goethert rule accounts for compressibility. The method has been programmed for the CDC 6600. Delta wings, gothic wings, arrow wings, cambered wings, and wing with body have been analyzed. Initial studies involving variations of panel density, vortex sheet sizing, Jacobian update, and initial geometry demonstrate that the present method generally exhibits good convergence characteristics.

Rubbert, p. E.↗

Estimate of Orbiter static aeroelasticity properties via FLEXSTAB

The theory underlying the use of rigid-model experimental data in FLEXSTAB aeroelasticity calculations is summarized. A description is given of the aerodynamic and structural models used in the calculation of Space Shuttle Orbiter aeroelasticity properties during the re-entry phase. Experience gained with the use of pressure, nonlinear force and moment, and area-ratio correction methods in the Orbiter application is described and typical program output data are included. The inability of the program to correctly estimate differential control surface pressures and to accept external corrections for these terms is shown to be a major limitation in the Orbiter application.

Cavin, R. K., III↗

Hydrodynamic stability problem formulated by numerical solutions of the Navier-Stokes equations

The problem of hydrodynamic stability and the transition from laminar to turbulent flows are reformulated by seeking numerical solutions of the full, unapproximated Navier-Stokes equations. This method differs significantly from the well known Orr-Sommerfeld equation approach. The oncoming laminar flow is disturbed by forced, time dependent perturbations. The magnitudes of these perturbations are arbitrary. Then, the ensuing spatial and temporal development of the imposed perturbations on the basic flow is calculated by direct numerical solutions of the time dependent Navier-Stokes equations. One of the main advantages of this method is its ability to simulate nonlinear processes. As a specific application of this technique to SSME (Space Shuttle Main Engine) flow configurations, computer programs have been written for the two dimensional flow over a backward facing step. This numerical code will be tested for operational use as part of a continued research collaboration effort with the NASA/MSFC counterparts.

Hyun, J. M.↗

Paneling methods with vorticity effects and corrections for nonlinear compressibility

The prediction of missile aerodynamic characteristics is presently undertaken through the application of supersonic paneling methods and nonlinear corrections to the prediction of missile aerodynamic characteristics. Attention is given to supersonic panel methods and line-singularity methods for the modeling of axisymmetric bodies, in combination with corrections for nonlinear flow phenomena, which are applied to complete missile, inlets, and wing-body combinations. The LRCDM2 computer program is used as an example of the methods presented.

Dillenius, Marnix F. E.↗

Turbomachinery aeroelasticity at NASA Lewis Research Center

The turbomachinery aeroelastic effort is focused on unstalled and stalled flutter, forced response, and whirl flutter of both single rotation and counter rotation propfans. It also includes forced response of the Space Shuttle Main Engine (SSME) turbopump blades. Because of certain unique features of propfans and the SSME turbopump blades, it is not possible to directly use the existing aeroelastic technology of conventional propellers, turbofans or helicopters. Therefore, reliable aeroelastic stability and response analysis methods for these propulsion systems must be developed. The development of these methods for propfans requires specific basic technology disciplines, such as 2-D and 3-D steady and unsteady aerodynamic theories in subsonic, transonic and supersonic flow regimes; modeling of composite blades; geometric nonlinear effects; and passive and active control of flutter and response. These methods are incorporated in a computer program, ASTROP. The program has flexibility such that new and future models in basic disciplines can be easily implemented.

Kaza, Krishna Rao V.↗