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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Adaptive methods for nonlinear structural dynamics and crashworthiness analysis

The objective is to describe three research thrusts in crashworthiness analysis: adaptivity; mixed time integration, or subcycling, in which different timesteps are used for different parts of the mesh in explicit methods; and methods for contact-impact which are highly vectorizable. The techniques are being developed to improve the accuracy of calculations, ease-of-use of crashworthiness programs, and the speed of calculations. The latter is still of importance because crashworthiness calculations are often made with models of 20,000 to 50,000 elements using explicit time integration and require on the order of 20 to 100 hours on current supercomputers. The methodologies are briefly reviewed and then some example calculations employing these methods are described. The methods are also of value to other nonlinear transient computations.

Belytschko, Ted↗
NASA NTRS1994

Development of the NASA/FLAGRO computer program for analysis of airframe structures

The NASA/FLAGRO (NASGRO) computer program was developed for fracture control analysis of space hardware and is currently the standard computer code in NASA, the U.S. Air Force, and the European Agency (ESA) for this purpose. The significant attributes of the NASGRO program are the numerous crack case solutions, the large materials file, the improved growth rate equation based on crack closure theory, and the user-friendly promptive input features. In support of the National Aging Aircraft Research Program (NAARP); NASGRO is being further developed to provide advanced state-of-the-art capability for damage tolerance and crack growth analysis of aircraft structural problems, including mechanical systems and engines. The project currently involves a cooperative development effort by NASA, FAA, and ESA. The primary tasks underway are the incorporation of advanced methodology for crack growth rate retardation resulting from spectrum loading and improved analysis for determining crack instability. Also, the current weight function solutions in NASGRO or nonlinear stress gradient problems are being extended to more crack cases, and the 2-d boundary integral routine for stress analysis and stress-intensity factor solutions is being extended to 3-d problems. Lastly, effort is underway to enhance the program to operate on personal computers and work stations in a Windows environment. Because of the increasing and already wide usage of NASGRO, the code offers an excellent mechanism for technology transfer for new fatigue and fracture mechanics capabilities developed within NAARP.

Forman, R. G.↗
NASA NTRS2013

Nonlinear Analysis and Post-Test Correlation for a Curved PRSEUS Panel

The Pultruded Rod Stitched Efficient Unitized Structure (PRSEUS) concept, developed by The Boeing Company, has been extensively studied as part of the National Aeronautics and Space Administration's (NASA s) Environmentally Responsible Aviation (ERA) Program. The PRSEUS concept provides a light-weight alternative to aluminum or traditional composite design concepts and is applicable to traditional-shaped fuselage barrels and wings, as well as advanced configurations such as a hybrid wing body or truss braced wings. Therefore, NASA, the Federal Aviation Administration (FAA) and The Boeing Company partnered in an effort to assess the performance and damage arrestments capabilities of a PRSEUS concept panel using a full-scale curved panel in the FAA Full-Scale Aircraft Structural Test Evaluation and Research (FASTER) facility. Testing was conducted in the FASTER facility by subjecting the panel to axial tension loads applied to the ends of the panel, internal pressure, and combined axial tension and internal pressure loadings. Additionally, reactive hoop loads were applied to the skin and frames of the panel along its edges. The panel successfully supported the required design loads in the pristine condition and with a severed stiffener. The panel also demonstrated that the PRSEUS concept could arrest the progression of damage including crack arrestment and crack turning. This paper presents the nonlinear post-test analysis and correlation with test results for the curved PRSEUS panel. It is shown that nonlinear analysis can accurately calculate the behavior of a PRSEUS panel under tension, pressure and combined loading conditions.

Gould, Kevin↗
NASA NTRS1982

Tenth NASTRAN User's Colloquium

The development of the NASTRAN computer program, a general purpose finite element computer code for structural analysis, was discussed. The application and development of NASTRAN is presented in the following topics: improvements and enhancements; developments of pre and postprocessors; interactive review system; the use of harmonic expansions in magnetic field problems; improving a dynamic model with test data using Linwood; solution of axisymmetric fluid structure interaction problems; large displacements and stability analysis of nonlinear propeller structures; prediction of bead area contact load at the tire wheel interface; elastic plastic analysis of an overloaded breech ring; finite element solution of torsion and other 2-D Poisson equations; new capability for elastic aircraft airloads; usage of substructuring analysis in the get away special program; solving symmetric structures with nonsymmetric loads; evaluation and reduction of errors induced by Guyan transformation.

Source record↗
DOE OSTI2021

Simulating non-native cubic interactions on noisy quantum machines

As a milestone for general-purpose computing machines, we demonstrate that quantum processors can be programed to efficiently simulate dynamics that are not native to the hardware. Moreover, on noisy devices without error correction, we show that simulation results are significantly improved when the quantum program is compiled using modular gates instead of a restricted set of standard gates. We demonstrate the general methodology by solving a cubic interaction problem, which appears in nonlinear optics, gauge theories, as well as plasma and fluid dynamics. To encode the non-native Hamiltonian evolution, we decompose the Hilbert space into a direct sum of invariant subspaces in which the nonlinear problem is mapped to a finite-dimensional Hamiltonian simulation problem. Furthermore, in a three-states example, the resultant unitary evolution is realized by a product of approximately 20 standard gates, using which approximately ten simulation steps can be carried out on state-of-the-art quantum hardware before results are corrupted by decoherence. In comparison, the simulation depth is improved by more than an order of magnitude when the unitary evolution is realized as a single cubic gate, which is compiled directly using optimal control. Alternatively, parametric gates may also be compiled by interpolating control pulses. Modular gates thus obtained provide high-fidelity building blocks for quantum Hamiltonian simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1978

Sonic boom research

It is demonstrated that a supersonic airplane configuration weighing over half a million pounds while creating a maximum sonic boom of less than 1 p.s.f. can be designed. New experimental techniques are developed in the wind tunnel and experiments for the sonic boom measurements were carried out. Theoretical analyses were performed for the effects of sonic boom on structures and pollution problems associated with supersonic flights were investigated. Numerical programs were generated for the sonic boom propagations from the near field of an airplane in supersonic flight at high altitude to the ground, taking into account the nonlinear effects and the asymmetric effects due to lift and the spacewise distributions of lift and volume.

Ting, L.↗
NASA NTRS1981

EET theoretical design techniques

As a part of the EET aerodynamics program an out-of-house program was developed and monitored to provide theoretical procedures useful in the design of transport aircraft. The focus of the effort was to provide tools valid in the nonlinear transonic speed range. The effort was divided into two basic areas, inviscid configuration analysis and design procedures and viscous correction procedures.

Dwoyer, D. L.↗
NASA NTRS1985

Large deflections of a cantilever beam under arbitrarily directed tip load

The nonlinear beam equation was integrated numerically in a direct fashion to obtain results for large deflections of cantilevers under tip loads of arbitrary direction. A short BASIC computer program for performing this integration is presented. Results for selected load cases are presented. The numerical process is performed rapidly on a modern microcomputer, and comparisons with results from closed form solutions show that the process is accurate.

Mccomb, H. E., Jr.↗
NASA NTRS2016

Validation of High-Resolution CFD Method for Slosh Damping Extraction of Baffled Tanks

Determination of slosh damping is a very challenging task as there is no analytical solution. The damping physics involve the vorticity dissipation which requires the full solution of the nonlinear Navier-Stokes equations. As a result, previous investigations and knowledge were mainly carried out by extensive experimental studies. A Volume-Of-Fluid (VOF) based CFD program developed at NASA MSFC was applied to extract slosh damping in a baffled tank from the first principle. First, experimental data using water with subscale smooth wall tank were used as the baseline validation. CFD simulation was demonstrated to be capable of accurately predicting natural frequency and very low damping value from the smooth wall tank at different fill levels. The damping due to a ring baffle at different liquid fill levels from barrel section and into the upper dome was then investigated to understand the slosh damping physics due to the presence of a ring baffle. Based on this study, the Root-Mean-Square error of our CFD simulation in estimating slosh damping was less than 4.8%, and the maximum error was less than 8.5%. Scalability of subscale baffled tank test using water was investigated using the validated CFD tool, and it was found that unlike the smooth wall case, slosh damping with baffle is almost independent of the working fluid and it is reasonable to apply water test data to the full scale LOX tank when the damping from baffle is dominant. On the other hand, for the smooth wall, the damping value must be scaled according to the Reynolds number. Comparison of experimental data, CFD, with the classical and modified Miles equations for upper dome was made, and the limitations of these semi-empirical equations were identified.

Yang, H. Q.↗
NASA NTRS2020

Cloud-top pressure retrieval with DSCOVR EPIC oxygen A- and B-band observations

An analytic transfer inverse model for Earth Polychromatic Imaging Camera (EPIC) observations is proposed to retrieve the cloud-top pressure (CTP) with the consideration of in-cloud photon penetration. In this model, an analytic equation was developed to represent the reflection at the top of the atmosphere from above cloud, in cloud, and below cloud. The coefficients of this analytic equation can be derived from a series of EPIC simulations under different atmospheric conditions using a nonlinear regression algorithm. With estimated cloud pressure thickness, the CTP can be retrieved from EPIC observation data by solving the analytic equation. To simulate the EPIC measurements, a program package using the double-k approach was developed. Compared to line-by-line calculation, this approach can calculate high-accuracy results with a 100-fold computation time reduction. During the retrieval processes, two kinds of retrieval results, i.e., baseline CTP and retrieved CTP, are provided. The baseline CTP is derived without considering in-cloud photon penetration, and the retrieved CTP is derived by solving the analytic equation, taking into consideration in-cloud and below-cloud interactions. The retrieved CTPs for the oxygen A and B bands are smaller than their related baseline CTP. At the same time, both baseline CTP and retrieved CTP at the oxygen B band are larger than those at the oxygen A band. Compared to the difference in baseline CTP between the B band and A band, the difference in retrieved CTP between these two bands is generally reduced. Out of around 10 000 cases, in retrieved CTP between the A and B bands we found an average bias of 93 mb with a standard deviation of 81 mb. The cloud layer top pressure from Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) measurements is used for validation. Under single-layer cloud situations, the retrieved CTPs for the oxygen A band agree well with the CTPs from CALIPSO, the mean difference of which within 5 mb in the case study. Under multiple-layer cloud situations, the CTPs derived from EPIC measurements may be larger than the CTPs of high-level thin clouds due to the effect of photon penetration.

DSCOVR↗
DOE OSTI2018

Control of Three Degrees-of-Freedom Wave Energy Converters Using Pseudo-Spectral Methods

This study presents a solution to the optimal control problem of a three degrees-of-freedom (3DOF) wave energy converter (WEC). The three modes are the heave, pitch, and surge. The dynamic model is characterized by a coupling between the pitch and surge modes, while the heave is decoupled. The heave, however, excites the pitch motion through nonlinear parametric excitation in the pitch mode. This paper uses Fourier series (FS) as basis functions to approximate the states and the control. A simplified model is first used where the parametric excitation term is neglected and a closed-form solution for the optimal control is developed. For the parametrically excited case, a sequential quadratic programming approach is implemented to solve for the optimal control numerically. Numerical results show that the harvested energy from three modes is greater than three times the harvested energy from the heave mode alone. Moreover, the harvested energy using a control that accounts for the parametric excitation is significantly higher than the energy harvested when neglecting this nonlinear parametric excitation term.

16 TIDAL AND WAVE POWER↗
NASA NTRS1994

Techniques to improve maneuver stability characteristics of a nonlinear wide-body transport airplane in cruise flight

The maneuver control stability characteristics of an aircraft are a flying qualities parameter of critical importance, to ensure structural protection as well as adequate predictability to the pilot. Currently, however, maneuver stability characteristics are not uniquely addressed in the Federal Aviation Regulations (FAR) Part 25, for transport aircraft. In past transport category certification programs, the Federal Aviation Administration (FAA) has used a combination of requirements (longitudinal control, vibration and buffeting, high-speed characteristics, and out-of-trim characteristics) to ensure safe and controllable maneuver stability characteristics over a range of flight conditions and airplane configurations. Controversies exist regarding each of these regulations, however, and considerable expenditures in terms of design studies and testing time have resulted from the requirements. It is also recognized that additional engineering guidance is needed for identifying acceptable nonlinear maneuver stability characteristics, particularly as they relate to relaxed stability, highly augmented transport configurations. The current trend in large aircraft design is toward relaxed, or even negative, static margins for improved fuel efficiency. The advanced flight control systems developed for these aircraft, in many instances, have rendered current aforementioned maneuver stability criteria either too stringent or of little practical use. Current design requirements do not account for these advanced designs. The objective was to evaluate a broad spectrum of linear and nonlinear longitudinal stability characteristics to generate data for defining satisfactory and unacceptable maneuver characteristics, as defined by pilot opinion. Primary emphasis was placed on two techniques of varying column force per normal acceleration. This study was a joint venture with four pilots participating; one from NASA, one from the FAA, and two from industry.

Grantham, William D.↗
DOE OSTI2015

Minimax terminal approach problem in two-level hierarchical nonlinear discrete-time dynamical system

We consider a discrete–time dynamical system consisting of three controllable objects. The motions of all objects are given by the corresponding vector nonlinear or linear discrete–time recurrent vector relations, and control system for its has two levels: basic (first or I level) that is dominating and subordinate level (second or II level) and both have different criterions of functioning and united a priori by determined informational and control connections defined in advance. For the dynamical system in question, we propose a mathematical formalization in the form of solving a multistep problem of two-level hierarchical minimax program control over the terminal approach process with incomplete information and give a general scheme for its solving.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Improve earthquake hypocenter using adaptive simulated annealing inversion in regional tectonic, volcano tectonic, and geothermal observation

Observation of earthquakes is routinely used widely in tectonic activity observation, and also in local scale such as volcano tectonic and geothermal activity observation. It is necessary for determining the location of precise hypocenter which the process involves finding a hypocenter location that has minimum error between the observed and the calculated travel times. When solving this nonlinear inverse problem, simulated annealing inversion method can be applied to such global optimization problems, which the convergence of its solution is independent of the initial model. In this study, we developed own program codeby applying adaptive simulated annealing inversion in Matlab environment. We applied this method to determine earthquake hypocenter using several data cases which are regional tectonic, volcano tectonic, and geothermal field. The travel times were calculated using ray tracing shooting method. We then compared its results with the results using Geiger’s method to analyze its reliability. Our results show hypocenter location has smaller RMS error compared to the Geiger’s result that can be statistically associated with better solution. The hypocenter of earthquakes also well correlated with geological structure in the study area. Werecommend using adaptive simulated annealing inversion to relocate hypocenter location in purpose to get precise and accurate earthquake location.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Combined Bernstein Polynomial, Optimal Reciprocal Collision Avoidance, Differential Dynamic Programming for Trajectory Replanning and Collision Avoidance for UAM Vehicles

This paper presents an integration of Differential Dynamic Programming (DDP) with the Optimal Reciprocal Collision Avoidance (ORCA) algorithm as the basis for a new algorithm, titled Combined Bernstein Polynomial Optimal Reciprocal Collision Avoidance DDP (COBRA-DDP), for trajectory replanning and collision avoidance for Urban Air Mobility (UAM) vehicles. State-constrained variants of DDP provide the ability to plan trajectories while avoiding obstacles, but these methods require a large increase in computational time per iteration which hinders the overall speed of the algorithm. ORCA utilizes simplified dynamics to recognize potential collisions along a trajectory and provides an optimal velocity for the avoidance of multiple vehicles. These velocity commands, however, may not result in a dynamically feasible trajectory for DDP to plan around. As such, a Bernstein polynomial curve that considers the general dynamic constraints of the vehicle is generated to approximate a trajectory based on the velocity commands. COBRA-DDP optimizes this suggested trajectory via unconstrained DDP to provide a dynamically feasible trajectory that provides collision avoidance. This new trajectory can be applied to the vehicle or used to warm start the state constrained DDP algorithms to decrease computation time. Its benefits and effectiveness of the algorithm are demonstrated on a UAM Vertical Takeoff and Landing (VTOL) vehicle simulation with highly nonlinear dynamics.

Optimal Reciprocal Collision Avoidance↗
DOE OSTI2015

3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium

Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D aspect, the code was massively parallelized using the single program, multiple data paradigm with the Message Passing Interfaces (MPI) for distributed memory architectures. This allows us to handle problems in the order of a thousand billion mesh points in the four dimensions (3 dimensions of space plus time). The validity of the method has been thoroughly evaluated on many cases with known solutions: linear piston, scattering of plane wave by a heterogeneous sphere, propagation in a waveguide with a shear flow, scattering by a finite amplitude vortex and nonlinear propagation in a thermoviscous medium. This validation process allows for a detailed assessment of the advantages and limitations of the method. Finally, applications to atmospheric propagation of shock waves will be presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

Analytical and finite element solutions of some problems using a viscoplastic model

The feasibility of using a viscoplastic model developed by Robinson to perform a nonlinear structural analysis was investigated. The paper presents analytical solutions for three classical problems: (1) a pressurized thick-walled cylinder, (2) a thin rotating disk, and (3) a pressurized thick-walled sphere. The analytical expressions derived for the stress and the strain rates for these components are general in nature as they consider both the mechanical and thermal loadings to be time-dependent. A computer program VISTAN (VIscoplastic STructural ANalyzer) was developed to obtain the stress and strain distributions. The finite element solutions for these problems are also presented. The numerical results pertaining to isothermal loading conditions are provided. The results obtained demonstrate the feasibility of using the viscoplastic model developed by Robinson to perform nonlinear structural analyses.

Arya, V. K.↗
DOE OSTI2024

HHL algorithm with mapping function and enhanced sampling for model predictive control in microgrids

Here, this paper presents a refined quantum Harrow Hassidim Lloyd (HHL) algorithm for microgrid control. The first novelty of the developed method is that a mapping shift function enables the original HHL algorithm to handle general linear equations with non-singular and indefinite matrix. Second, a method of Matrix Extension for Amplifying Sampling Probabilities of Intended Solution (ME-ASPI) is proposed to design the reformulated linear algebraic equations, allowing for improved sampling efficiency of the quantum tomography in the refined HHL algorithm. Then, we applied the method to solve the model predictive control (MPC) problem in nonlinear dynamical microgrids. Specifically, with the ME-ASPI method, the refined HHL algorithm can effectively obtain the intended partial optimal control inputs for MPC. The optimization of quadratic programming problem in each time step of MPC is transformed into a linear system problem, which is addressed by the proposed quantum solver through using only partial information, with the time complexity improved from $\mathscr{O}(\mathscr{N}^{2.37286})$ classically to $\mathscr{O}(\mathscr{N}^{2} log \mathscr{N}$ x $p$ log $p)$ in quantum. Numerical examples have validated the effectiveness of the refined HHL algorithm with the proposed mapping function and the ME-ASPI method. By leveraging quantum properties, the proposed method provides a hybrid quantum–classical framework for microgrid control. This generic method can also potentially tackle many other challenges in analyzing and controlling general complex engineered systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗