Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,639 records · Page 91

Advances in Multi-Physics Modeling for Parallel Finite Element Code Suite ACE3P

SLAC developed ACE3P code suite is a comprehensive set of parallel finite-element codes for multi-physics modeling of particle accelerators. Running on massively parallel computer platforms for high fidelity and high accuracy simulation, ACE3P enables rapid virtual prototyping of accelerator and RF component design, optimization and analysis. Furthermore, recent advances of ACE3P have been focused on the implementation of advanced numerical algorithms, enhancement of multi-physics capabilities, and improvement of code performance on state-of-the-art high performance computing (HPC) platforms for large-scale RF modeling in accelerator applications and beyond. In this paper, we will present a nonlinear electromagnetic (EM) eigensolver for damping calculation of resonant modes for cavities coupled with waveguides, the thermal solver for accurate evaluation of heat loads in superconducting RF (SRF) cavity power coupler, the elastic solver for investigation of SRF cavity cryomodule (CM) bowing due to thermal effect, the frequency domain mechanical solver for studying SRF cavity mechanical oscillation modes, the particle tracking code for dark current simulation in SRF CM using a hybrid MPI+OpenMP parallel programing model, the time domain solver for high fidelity modeling of EM wave propagation for medical device in human body.

97 MATHEMATICS AND COMPUTING↗

Simulation studies using multibody dynamics code DART

DART is a multibody dynamics code developed by Photon Research Associates for the Air Force Astronautics Laboratory (AFAL). The code is intended primarily to simulate the dynamics of large space structures, particularly during the deployment phase of their missions. DART integrates nonlinear equations of motion numerically. The number of bodies in the system being simulated is arbitrary. The bodies' interconnection joints can have an arbitrary number of degrees of freedom between 0 and 6. Motions across the joints can be large. Provision for simulating on-board control systems is provided. Conservation of energy and momentum, when applicable, are used to evaluate DART's performance. After a brief description of DART, studies made to test the program prior to its delivery to AFAL are described. The first is a large angle reorientating of a flexible spacecraft consisting of a rigid central hub and four flexible booms. Reorientation was accomplished by a single-cycle sine wave shape torque input. In the second study, an appendage, mounted on a spacecraft, was slewed through a large angle. Four closed-loop control systems provided control of this appendage and of the spacecraft's attitude. The third study simulated the deployment of the rim of a bicycle wheel configuration large space structure. This system contained 18 bodies. An interesting and unexpected feature of the dynamics was a pulsing phenomena experienced by the stays whole playout was used to control the deployment. A short description of the current status of DART is given.

Keat, James E.↗

Loads Model Development and Analysis for the F/A-18 Active Aeroelastic Wing Airplane

The Active Aeroelastic Wing airplane was successfully flight-tested in March 2005. During phase 1 of the two-phase program, an onboard excitation system provided independent control surface movements that were used to develop a loads model for the wing structure and wing control surfaces. The resulting loads model, which was used to develop the control laws for phase 2, is described. The loads model was developed from flight data through the use of a multiple linear regression technique. The loads model input consisted of aircraft states and control surface positions, in addition to nonlinear inputs that were calculated from flight-measured parameters. The loads model output for each wing consisted of wing-root bending moment and torque, wing-fold bending moment and torque, inboard and outboard leading-edge flap hinge moment, trailing-edge flap hinge moment, and aileron hinge moment. The development of the Active Aeroelastic Wing loads model is described, and the ability of the model to predict loads during phase 2 research maneuvers is demonstrated. Results show a good match to phase 2 flight data for all loads except inboard and outboard leading-edge flap hinge moments at certain flight conditions. The average load prediction errors for all loads at all flight conditions are 9.1 percent for maximum stick-deflection rolls, 4.4 percent for 5-g windup turns, and 7.7 percent for 4-g rolling pullouts.

Allen, Michael J.↗

Loads Model Development and Analysis for the F/A-18 Active Aeroelastic Wing Airplane

The Active Aeroelastic Wing airplane was successfully flight-tested in March 2005. During phase 1 of the two-phase program, an onboard excitation system provided independent control surface movements that were used to develop a loads model for the wing structure and wing control surfaces. The resulting loads model, which was used to develop the control laws for phase 2, is described. The loads model was developed from flight data through the use of a multiple linear regression technique. The loads model input consisted of aircraft states and control surface positions, in addition to nonlinear inputs that were calculated from flight-measured parameters. The loads model output for each wing consisted of wing-root bending moment and torque, wing-fold bending moment and torque, inboard and outboard leading-edge flap hinge moment, trailing-edge flap hinge moment, and aileron hinge moment. The development of the Active Aeroelastic Wing loads model is described, and the ability of the model to predict loads during phase 2 research maneuvers is demonstrated. Results show a good match to phase 2 flight data for all loads except inboard and outboard leading-edge flap hinge moments at certain flight conditions. The average load prediction errors for all loads at all flight conditions are 9.1 percent for maximum stick-deflection rolls, 4.4 percent for 5-g windup turns, and 7.7 percent for 4-g rolling pullouts.

Allen, Michael J.↗

Polarized Bidirectional Reflectance of Optically Thick Sparse Particulate Layers: an Efficient Numerically Exact Radiative-Transfer Solution

We describe a simple yet efficient numerical algorithm for computing polarized bidirectional reflectance of an optically thick (semi-infinite), macroscopically flat layer composed of statistically isotropic and mirror symmetric random particles. The spatial distribution of the particles is assumed to be sparse, random, and statistically uniform. The 44 Stokes reflection matrix is calculated by iterating the Ambartsumian's vector nonlinear integral equation. The result is a numerically exact solution of the vector radiative transfer equation and as such fully satisfies the energy conservation law and the fundamental reciprocity relation. Since this technique bypasses the computation of the internal radiation field, it is very fast and highly accurate. The FORTRAN implementation of the technique is publicly available on the World Wide Web at http://www.giss.nasa.gov/staff/ mmishchenko/brf. It can be combined with several existing computer programs providing the requisite single-scattering properties of spherical or morphologically complex particles and applied to a wide range of optical characterization problems. Benchmark results obtained with this program can be used for testing alternative solvers of the vector radiative transfer equation.

radiative transfer↗

Flight Performance and Stability of Space Launch System Core Stage Thrust Vector Control

The Space Launch System (SLS) Core Stage (CS) Thrust Vector Control (TVC) system is comprised of eight mechanical feedback Shuttle heritage Type III TVC actuators and four RS-25 engines, each attached to a Shuttle heritage gimbal block/bearing. The Core Stage TVC shares vehicle control authority with the SLS 5-segment Solid Rocket Boosters (SRBs) during boost phase flight, and is the sole means of vehicle flight control during in exoatmospheric flight following SRB separation. TVC responses during Green Run Hot Fire (GRHF) testing revealed that the TVC did not meet its performance specifications. Step and frequency responses exhibited unexpected departures from prior laboratory data and modeled behavior. Post-test analysis determined that the characteristics of the structure and gimbal friction are significantly influenced by the thrust-loaded conditions, and the command avionics exhibited a small but important gain nonlinearity. Using the available test data, the design team augmented the flight control TVC models to bound the observed results and include the additional fidelity needed for vehicle flight control analysis so as to build sufficient rationale for flight certification. Prior to the Green Run tests, “simplex” linear models typically used for flight control analysis did not include gimbal friction and other nonlinearities owing to long-standing assumptions that these effects were negligible in the Shuttle Orbiter TVC system. Following the Green Run findings, simulation analysis of the flight dynamics in the time and frequency domain revealed the propensity for a flight control limit cycle oscillation (LCO) if friction and structural compliances fell near the edges of test-predicted bounds. While the “most probable” models did not predict an in-flight LCO, the SLS Program conservatively proceeded with a system-wide evaluation and ultimate acceptance of the possibility for a small amplitude, low-frequency TVC LCO in flight. A final validation of the extensive test and modeling effort occurred when the first flight of SLS successfully demonstrated the fully integrated performance of the vehicle’s TVC system This paper is the final installment in a seven-paper series surveying the design, engineering, test validation, and flight performance of the Core Stage Thrust Vector Control system. In this paper, the development of flight rationale in light of the TVC responses observed in Green Run is discussed, along with a review of the flight telemetry illustrating the correlation of the preflight predictions with the observed performance.

John H. Wall↗

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗

SCI Identification (SCIDNT) program user's guide

The computer program Linear SCIDNT which evaluates rotorcraft stability and control coefficients from flight or wind tunnel test data is described. It implements the maximum likelihood method to maximize the likelihood function of the parameters based on measured input/output time histories. Linear SCIDNT may be applied to systems modeled by linear constant-coefficient differential equations. This restriction in scope allows the application of several analytical results which simplify the computation and improve its efficiency over the general nonlinear case.

Source record↗

Derivatives of buckling loads and vibration frequencies with respect to stiffness and initial strain parameters

A uniform variational approach to sensitivity analysis of vibration frequencies and bifurcation loads of nonlinear structures is developed. Two methods of calculating the sensitivities of bifurcation buckling loads and vibration frequencies of nonlinear structures, with respect to stiffness and initial strain parameters, are presented. A direct method requires calculation of derivatives of the prebuckling state with respect to these parameters. An adjoint method bypasses the need for these derivatives by using instead the strain field associated with the second-order postbuckling state. An operator notation is used and the derivation is based on the principle of virtual work. The derivative computations are easily implemented in structural analysis programs. This is demonstrated by examples using a general purpose, finite element program and a shell-of-revolution program.

Haftka, Raphael T.↗

Resource Letter GrW-1: Gravitational Waves

The phenomenon of gravitational radiation was one of the first predictions of Einstein's general theory of relativity. Progress in understanding this radiation theoretically was slow at first, owing to the difficulty of the nonlinear field equations and the subtleties of their physical effects. The experimental side of this subject also has taken a long time to develop, with efforts at detection severely challenged by the extreme weakness of the waves impinging on the Earth. However, as the 21st century begins, observations of the gravitational waves from astrophysical sources such as black holes, neutron stars, and stellar collapse are expected to open a new window on the universe. Vigorous experimental programs centered on ground-based detectors are being carried out worldwide, and a space-based detector is in the planning stages. On the theoretical side, much effort is being expended to produce robust models of the astrophysical sources and accurate calculations of the waveforms they produce. In this Resource Letter, a set of basic references will be presented first, to provide a general introduction to and overview of the literature in this field. The focus then will shift to highlighting key resources in more specialized areas at the forefront of current research.

White, Nicholas E.↗

Pressure Testing of a Minimum Gauge PRSEUS Panel

Advanced aircraft configurations that have been developed to increase fuel efficiency require advanced, novel structural concepts capable of handling the unique load conditions that arise. One such concept is the Pultruded Rod Stitched Efficient Unitized Structure (PRSEUS) developed by the Boeing Company. The PRSEUS concept is being investigated by NASA s Environmentally Responsible Aviation (ERA) Program for use in a hybrid-wing body (HWB) aircraft. This paper summarizes the analysis and test of a PRSEUS panel subjected to internal pressure, the first such pressure test for this structural concept. The pressure panel used minimum gauge skin, with stringer and frame configurations consistent with previous PRSEUS tests. Analysis indicated that for the minimum gauge skin panel, the stringer locations exhibit fairly linear response, but the skin bays between the stringers exhibit nonlinear response. Excellent agreement was seen between nonlinear analysis and test results in the critical portion at the center of the panel. The pristine panel was capable of withstanding the required 18.4 psi pressure load condition without exhibiting any damage. The impacted panel was capable of withstanding a pressure load in excess of 28 psi before initial failure occurred at the center stringer, and the panel was capable of sustaining increased pressure load after the initial failure. This successful PRSEUS panel pressure panel test was a critical step in the building block approach for enabling the use of this advanced structural concept on future aircraft, such as the HWB.

Lovejoy, Andrew J.↗

Solving the Dynamics-Aware Economic Dispatch Problem with the Koopman Operator

The dynamics-aware economic dispatch (DED) problem embeds low-level generator dynamics and operational constraints to enable near real-time scheduling of generation units in a power network. DED produces a more dynamic supervisory control policy than traditional economic dispatch (T-ED) that reduces overall generation costs. However, in contrast to T-ED, DED is a nonlinear, non-convex optimization problem that is computationally prohibitive to solve. We introduce a machine learning-based operator-theoretic approach for solving the DED problem efficiently. Specifically, we develop a novel discrete-time Koopman Operator (KO) formulation that embeds domain information into the structure of the KO to learn high-fidelity approximations of the generator dynamics. Using the KO approximation, the DED problem can be reformulated as a computationally tractable linear program (abbreviated DED-KO). We demonstrate the high solution quality and computational-time savings of the DED-KO model over the original DED formulation on a 9-bus test system.

King, Ethan↗

Economic Optimal Control for Minimizing Fuel Consumption of Heavy-Duty Trucks in a Highway Environment

This paper provides a comparative assessment of three economic optimal control strategies, aimed at minimizing the fuel consumption of heavy-duty trucks in a highway environment, under a representative lead vehicle model informed by traffic data. These strategies fuse a global, offline dynamic programming (DP) optimization with online model predictive control (MPC). We then show how two of the three strategies can be adapted to accommodate the presence of traffic and optimally navigate signalized intersections using infrastructure-to-vehicular (I2V) communication. The MPC optimization, which is local in nature, makes refinements to a coarsely (but globally, subject to grid resolution) optimized target velocity profile from the DP optimization. The three candidate economic MPC formulations that are evaluated include: a nonlinear time-based formulation that directly penalizes predicted fuel consumption, a nonlinear time-based formulation that penalizes braking effort as a surrogate for fuel consumption, and a linear distance-based convex formulation that maintains a tradeoff between energy expenditure and tracking of the coarsely optimized velocity profile obtained from DP. Using a medium-fidelity Simulink model, based on a Volvo truck's longitudinal and engine dynamics, we analyze the optimization's performance on four highway routes under various traffic scenarios. Results demonstrate 3.7-8.3% fuel economy improvement on highway routes without traffic and 6.5-10% on the same routes with traffic included. Furthermore, we present a detailed analysis of energy usage by "type" (aerodynamic losses, braking losses, and comparison of brake-specific fuel consumption), under each candidate control strategy.

33 ADVANCED PROPULSION SYSTEMS↗

One- and Two-Qubit Gates: Rabi Technique and Single Unipolar Pulses

Methods for controlling states of interacting superconducting flux qubits using power-efficient devices of fast single-quantum logic (Josephson nonlinearity cavities) are studied. One- and two-qubit quantum logical operations performed within the conventional control technique using Rabi pulses and using picosecond single unipolar magnetic field pulses are comparatively analyzed. It is shown that all main operations can be implemented with an accuracy of better than 97% due to optimization of the shape and parameters of unipolar control pulses (associated with, e.g., propagation of fluxons in transmission lines). The efficiency of the developed technique for programming a two-qubit quantum processor implementing the simplest Deutsch–Jozsa algorithm is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Seismic Retrofit for Electric Power Systems

Our paper develops a two-stage stochastic program and solution procedure to optimize the selection of seismic retrofit strategies to increase the resilience of electric power systems against earthquake hazards. The model explicitly considers the range of earthquake events that are possible and, for each, an approximation of the distribution of damage experienced. Furthermore, this is important because electric power systems are spatially distributed and so their performance is driven by the distribution of component damage. We also test this solution procedure against the nonlinear integer solver in LINGO 13 and apply the formulation and solution strategy to the Eastern Interconnection, where seismic hazard stems from the New Madrid seismic zone.

58 GEOSCIENCES↗

Optimal feedback control for nonlinear control moment gyro systems

A study was made, leading to development and analysis of a technique allowing the control moment gyro (CMG) momentum vectors to avoid adverse orientations and thus to provide optimal control for the Apollo Telescope Mount. The technique was formulated in terms of a feedback scheme using a check factor, K, formed as the product of the dot products of the three CMG momentum unit vectors. The parameter K approaches unity for all adverse orientations and was used to initiate a single correction maneuver to rotate the CMG momenta away from a pending adverse orientation and to reduce K to an acceptable value. A digital computer program for ATM simulations was modified to incorporate this correction maneuver and performance was evaluated in simulated cases. Antiparallelism was used as the primary example for the evaluation of the feedback optimization technique.

Kyle, R. G.↗

Computer studies of baroclinic flow

Programs necessary for computing the transition curve on the regime diagram for the atmospheric general circulation experiment (AGOE) were completed and used to determine the regime diagram for the rotating annulus and some axisymmetric flows for one possible AGOE configuration. The effect of geometrical constraints on the size of eddies developing from a basic state is being examined. In AGOE, the geometric constraint should be the width of the shear zone or the baroclinic zone. Linear and nonlinear models are to be used to examine both barotropic and baroclinic flows. The results should help explain the scale selection mechanism of baroclinic eddies in the atmosphere experimental models such as AGOE, and the multiple vortex phenomenon in tornadoes.

Gall, R.↗

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