Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear programming problem”

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 523 records · Page 29

Parameter Design and Optimal Control of an Open Core Flywheel Energy Storage System

In low earth orbit (LEO) satellite applications spacecraft power is provided by photovoltaic cells and batteries. To overcome battery shortcomings the University of Maryland, working in cooperation with NASA/GSFC and NASA/LeRC, has developed a magnetically suspended flywheel for energy storage applications. The system is referred to as an Open Core Composite Flywheel (OCCF) energy storage system. Successful application of flywheel energy storage requires integration of several technologies, viz. bearings, rotor design, motor/generator, power conditioning, and system control. In this paper we present a parameter design method which has been developed for analyzing the linear SISO model of the magnetic bearing controller for the OCCF. The objective of this continued research is to principally analyze the magnetic bearing system for nonlinear effects in order to increase the region of stability, as determined by high speed and large air gap control. This is achieved by four tasks: (1) physical modeling, design, prototyping, and testing of an improved magnetically suspended flywheel energy storage system, (2) identification of problems that limit performance and their corresponding solutions, (3) development of a design methodology for magnetic bearings, and (4) design of an optimal controller for future high speed applications. Both nonlinear SISO and MIMO models of the magnetic system were built to study limit cycle oscillations and power amplifier saturation phenomenon observed in experiments. The nonlinear models include the inductance of EM coils, the power amplifier saturation, and the physical limitation of the flywheel movement as discussed earlier. The control program EASY5 is used to study the nonlinear SISO and MIMO models. Our results have shown that the characteristics and frequency responses of the magnetic bearing system obtained from modeling are comparable to those obtained experimentally. Although magnetic saturation is shown in the bearings, there are good correlations between the theoretical model and experimental data. Both simulation and experiment confirm large variations of the magnetic bearing characteristics due to air gap growth. Therefore, the gap growth effect should be considered in the magnetic bearing system design. Additionally, the magnetic bearing control system will be compared to other design methods using not only parameter design but H-infinity optimal control and mu synthesis.

Pang, D.↗

Riemannian Optimization Applied to AC Optimal Power Flow: Preprint

The nonlinear, nonconvex AC optimal power flow problem is of growing importance as the nature of the power grid evolves. This problem can be difficult to solve for interior point methods. However, the advent of optimization algorithms over smooth Riemannian manifolds presents an alternative approach. The nonlinear, nonconvex constraints in the AC power flow problem form an embedded submanifold of Euclidean space. In this paper, the authors explore the performance of Riemannian optimization algorithms for the ACOPF problem where the optimization is performed directly on the AC power flow manifold. They demonstrate that these are viable computational alternatives to interior point methods. This is done by using Julia and the packages PowerModels.jl and Manopt.jl.

manifold optimization↗

The MODE family of on-orbit experiments: The Middeck Active Control Experiment (MACE)

A flight experiment entitled the Middeck Active Control Experiment (MACE), proposed by the Space Engineering Research Center (SERC) at the Massachusetts Institute of Technology, is described. This is the second in a family of flight experiments being developed at MIT. The first is the Middeck 0-Gravity Dynamics Experiment (MODE) which investigates the nonlinear behavior of contained fluids and truss structures in zero gravity. The objective of the MACE program is to investigate and validate the modeling of the dynamics of an actively controlled flexible, articulating, multibody platform free floating in zero gravity. A rationale and experimental approach for the program are presented. The rationale shows that on-orbit testing, coupled with ground testing and a strong analytical program, is necessary in order to fully understand both how flexibility of the platform affects the pointing problem, as well as how gravity perturbs this structural flexibility causing deviations between 1- and 0-gravity behavior. The experimental approach captures the essential physics of multibody platforms, by identifying the appropriate attributes, tests, and performance metrics of the test article and defines the tests required to successfully validate the analytical framework.

Crawley, Edward F.↗

Linear and nonlinear dynamic analysis of redundant load path bearingless rotor systems

The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.

Murthy, V. R.↗

Problems in large-scale structural optimization

A general design optimization model for large complex systems is defined. Major features of the model that challenge various optimization algorithms are discussed. Requirements of a model optimization algorithm are identified. Objectives of the study of various algorithms are defined and a basis for conducting such a study is developed. Primal as well as transformation methods are analytically studied and a unified viewpoint of various methods is presented. Several numerical examples are solved using different methods to study their performance. Conclusions drawn from the study are presented and discussed. Areas of future research in nonlinear programming as well as structural optimization are identified and discussed.

Arora, J. S.↗

Application of traditional CFD methods to nonlinear computational aeroacoustics problems

This paper describes an implementation of a high order finite difference technique and its application to the category 2 problems of the ICASE/LaRC Workshop on Computational Aeroacoustics (CAA). Essentially, a popular Computational Fluid Dynamics (CFD) approach (central differencing, Runge-Kutta time integration and artificial dissipation) is modified to handle aeroacoustic problems. The changes include increasing the order of the spatial differencing to sixth order and modifying the artificial dissipation so that it does not significantly contaminate the wave solution. All of the results were obtained from the CM5 located at the Numerical Aerodynamic Simulation Laboratory. lt was coded in CMFortran (very similar to HPF), using programming techniques developed for communication intensive large stencils, and ran very efficiently.

Chyczewski, Thomas S.↗

System modeling and optimal design of a Mars-roving vehicle.

The problem of systematically determining the optimal design for an unmanned Mars-roving vehicle is considered. A system model, identifying all feasible designs, is generated by consideration of the physical constraints on the design parameters, and the requirement that the system be deliverable to the Mars surface. An expression which evaluates system performance relative to mission goals is developed. The model and objective function together allow simulation of the effects of design trade-offs upon system performance for all feasible designs. Nonlinear programming techniques are utilized to identify the optimal design.

Smith, E. J.↗

Fabrication and testing of non-graphitic superhybrid composites

A study was conducted to determine the fabrication feasibility and the mechanical properties of adhesively-bonded boron aluminum/titanium and non-graphitic fiber/epoxy resin superhybrid (NGSH) composite laminates for potential aerospace applications. The major driver for this study was the elimination of a potential graphite fiber release problem in the event of a fire. The results of the study show that non-graphitic fibers, such as S-glass and Kevlar 49, may be substituted for the graphite fibers used in superhybrid (SH) composites for some applications. As is to be expected, however, the non-graphitic superhybrids have lower stiffness properties than the graphitic superhybrids. In-plane and flexural moduli of the laminates studied in this program can be predicted reasonably well using linear laminate theory while nonlinear laminate theory is required for strength predictions.

Lark, R. F.↗

Aeroelastic modeling of rotor blades with spanwise variable elastic axis offset: Classic issues revisited and new formulations

In response to a systematic methodology assessment program directed to the aeroelastic stability of hingeless helicopter rotor blades, improved basic aeroelastic reformulations and new formulations relating to structural sweep were achieved. Correlational results are presented showing the substantially improved performance of the G400 aeroelastic analysis incorporating these new formulations. The formulations pertain partly to sundry solutions to classic problem areas, relating to dynamic inflow with vortex-ring state operation and basic blade kinematics, but mostly to improved physical modeling of elastic axis offset (structural sweep) in the presence of nonlinear structural twist. Specific issues addressed are an alternate modeling of the delta EI torsional excitation due to compound bending using a force integration approach, and the detailed kinematic representation of an elastically deflected point mass of a beam with both structural sweep and nonlinear twist.

Bielawa, Richard L.↗

Minimum Weight Design of a Leaf Spring Tapered in Thickness and Width for the Hubble Space Telescope-Space Support Equipment

A linear elastic solution to the problem of minimum weight design of cantilever beams with variable width and depth is presented. The solution shown is for the specific application of the Hubble Space Telescope maintenance mission hardware. During these maintenance missions, delicate instruments must be isolated from the potentially damaging vibration environment of the space shuttle cargo bay during the ascent and descent phases. The leaf springs are designed to maintain the isolation system natural frequency at a level where load transmission to the instruments in a minimum. Nonlinear programming is used for the optimization process. The weight of the beams is the objective function with the deflection and allowable bending stress as the constraint equations. The design variables are the width and depth of the beams at both the free and the fixed ends.

Rodriguez, P. I.↗

A Mixed integer linear programming‐based distributed energy management for networked microgrids considering network operational objectives and constraints

Abstract Mixed integer linear programming (MILP)–based distributed energy management for networked microgrids embedded modern distribution systems is proposed. Considering the diverse ownership of microgrids, distributed energy resources (DERs) that interface directly with utilities and responsive loads, an alternating direction method of multipliers–based distributed framework was formulated for the scheduling of networked microgrids embedded modern distribution systems by adjusting nodal price signals iteratively. In addition, to make the formulated optimization problems resolvable through more accessible and popular MILP solvers, different linearisation techniques were employed to transform the nonlinear terms into linear or mixed integer linear formats. The proposed MILP‐based distributed method preserves all participants' autonomy (e.g., microgrids, DERs that interface directly with utilities and responsive loads), while incentivising them to actively participate in the distribution system operation with price signals. The proposed method is validated with results of numerical simulation using a modern distribution system consisting of multiple networked microgrids, DERs that interface directly with utilities, as well as responsive loads.

24 POWER TRANSMISSION AND DISTRIBUTION↗

New systems in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) serves as a common library of classes between applications developed for advanced reactor analysis, fusion device engineering, spent fuel cask analysis, geochemistry studies, among other fields. These applications drive the development of the framework to meet their needs. Systems in MOOSE group capabilities that share a common purpose and generally common code. They can be leveraged by all downstream applications, providing extensive code re-use and shared maintenance. They facilitate the discovery by new users of the classes meeting at least partially their needs, and offer the same opportunities for customization as other systems. The addition of a new system to MOOSE opens new ways of solving or discretizing nonlinear problems, of performing distributed postprocessing, and a plethora of other needs. While new systems can be introduced in downstream applications rather than at the framework level, the framework team monitors common needs across the community and often triggers their addition. Documentation, training material, development needs can be centralized, limiting duplicated work across the community. The last three years have seen a large expansion in the capabilities of MOOSE. The supporting role of the framework in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has created numerous feature requests to support neutronics, thermal hydraulics, computational fluid dynamics and thermo-mechanics simulations in the Griffin, SAM, Pronghorn and Bison applications respectively. Similarly, laboratory-directed research and development (LDRD) projects in additive manufacturing, high-Reynolds flow simulations, structure optimization also necessitate an expansion of the framework capabilities. This summary reports on the new systems created in MOOSE, their design, their capabilities and some of the relevant interfaces.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Efficient Region of Attraction Characterization for Control and Stabilization of Load Tap Changer Dynamics

In this article, we study the monitoring and control of long-term voltage stability considering load tap changer (LTC) dynamics. We show that under generic conditions, the LTC dynamics admit a unique stable equilibrium. For the stable equilibrium, we characterize an explicit inner approximation of its largest region of attraction (ROA). Compared to existing results, the computational complexity of the ROA characterization is drastically reduced. We propose a quadratically constrained linear program formulation for the ROA characterization problem. In addition, we formulate a second-order cone program for online voltage stability monitoring and control exploiting the proposed ROA characterization. Finally, we demonstrate the efficacy of the proposed formulations on the ROA characterization and stability monitoring and control using a standard IEEE test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Displacement incrementation in non-linear structural analysis by the self-correcting method

The prediction of nonlinear structural behavior by the finite element method wherein buckling does not occur has received considerable attention and, with it, reasonable success has been achieved. However, the post-buckling problem has been less actively pursued probably because of the inherent numerical difficulties encountered. This note reviews very briefly the numerical methods currently being used for pre- and post-buckling analysis and presents a self-correcting approach based on load and displacement incrementation which is shown to be efficient, reliable, and easy to program. Numerical solutions are presented which demonstrate the effectiveness of the method.

Haisler, W. E.↗

An initial investigation into methods of computing transonic aerodynamics sensitivity coefficients

Based upon the program outlined in the modified proposal, the initial effort has concentrated on developing the quasi-analytical approach for determining aerodynamic influence coefficients. In addition, in order to keep these initial efforts computationally efficient and straightforward, it was decided to first consider only two-dimensional flow and to model the problem using the transonic small perturbation equation with surface boundary conditions applied along the slit. The small perturbation equation was selected instead of the full potential equation because it permits the capturing of important nonlinear physical phenomena while being by comparison easy to handle from a coding standpoint. Consequently, the initial objectives were: (1) investigate the feasibility of using the quasi-analytical approach for determining aerodynamic sensitivity coefficients in the transonic regime; (2) determine the accuracy of the quasi-analytical approach by comparing results to those computed using finite difference techniques; and (3) investigate efficient methods for solving the large set of equations associated with the quasi-analytical approach.

Carlson, Leland A.↗

On the Application of a Response Surface Technique to Analyze Roll-over Stability of Capsules with Airbags Using LS-Dyna

As NASA moves towards developing technologies needed to implement its new Exploration program, studies conducted for Apollo in the 1960's to understand the rollover stability of capsules landing are being revisited. Although rigid body kinematics analyses of the roll-over behavior of capsules on impact provided critical insight to the Apollo problem, extensive ground test programs were also used. For the new Orion spacecraft being developed to implement today's Exploration program, new air-bag designs have improved sufficiently for NASA to consider their use to mitigate landing loads to ensure crew safety and to enable re-usability of the capsule. Simple kinematics models provide only limited understanding of the behavior of these air bag systems, and more sophisticated tools must be used. In particular, NASA and its contractors are using the LS-Dyna nonlinear simulation code for impact response predictions of the full Orion vehicle with air bags by leveraging the extensive air bag prediction work previously done by the automotive industry. However, even in today's computational environment, these analyses are still high-dimensional, time consuming, and computationally intensive. To alleviate the computational burden, this paper presents an approach that uses deterministic sampling techniques and an adaptive response surface method to not only use existing LS-Dyna solutions but also to interpolate from LS-Dyna solutions to predict the stability boundaries for a capsule on airbags. Results for the stability boundary in terms of impact velocities, capsule attitude, impact plane orientation, and impact surface friction are discussed.

Horta, Lucas G.↗

Influence of thermal distortion on the anomalous behavior of a gravity gradient satellite

Consideration of the problem of determining the influence of solar heating on the low-frequency dynamic stability of a gravity-gradient satellite. The scope of the effort includes the development of a computer program that uses a quasi-static approach for describing a satellite's unsteady orbital motion. The essential feature of the quasi-static approach is that the satellite's long, slender booms reach their thermal equilibrium position solely as a result of thermal bending without consideration of the dynamic effects of boom vibration. A nonlinear analytical model and a corresponding computer program have been developed to study the effect of absorptivity, sun angle, and hinge stiffness on the stability of a satellite's motion.

Goldman, R. L.↗

Ionosphere/microwave beam interaction study

A solar power satellite microwave power density of 20mw sq cm was confirmed as the level where nonlinear interactions may occur in the ionosphere, particularly at 100 km altitude. Radio wave heating at this altitude, produced at the Arecibo Observatory, yielded negative results for radio wave heating of an underdense ionosphere. Overdense heating produced striations in the ionosphere which may cause severe radio frequency interference problems under certain conditions. The effects of thermal self-focusing are shown to be limited severely geographically. The aspect sensitivity of field-aligned striations makes interference-free regions above magnetic latitude about 60 deg. A test program is proposed to simulate the interaction of the SPS beam with the ionosphere, to measure the effects of the interaction on the ionosphere and on communication and navigation systems, and to interpret the results.

Duncan, L. M.↗