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,531 records · Page 85

Topology Optimization for Nonlinear Transient Applications Using a Minimally Invasive Approach (LDRD Final Report)

The purpose of this project was to devise, implement, and demonstrate a method that can use Sandia's existing analysis codes (e.g., Sierra, Alegra, the CTH hydro code) with minimal modification to generate objective function gradients for optimization-based design in transient, non-linear, coupled-physics applications. The approach uses a Moving Least Squares representation of the geometry to substantially reduce the number of geometric degrees of freedom. A Multiple-Program Multiple-Data computing model is then used to compute objective gradients via finite differencing. Details of the formulation and implementation are provided, and example applications are presented that show effectiveness and scalability of the approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multi-level, forming and filament free, bulk switching trilayer RRAM for neuromorphic computing at the edge

CMOS-RRAM integration holds great promise for low energy and high throughput neuromorphic computing. However, most RRAM technologies relying on filamentary switching suffer from variations and noise, leading to computational accuracy loss, increased energy consumption, and overhead by expensive program and verify schemes. We developed a filament-free, bulk switching RRAM technology to address these challenges. We systematically engineered a trilayer metal-oxide stack and investigated the switching characteristics of RRAM with varying thicknesses and oxygen vacancy distributions to achieve reliable bulk switching without any filament formation. We demonstrated bulk switching at megaohm regime with high current nonlinearity, up to 100 levels without compliance current. We developed a neuromorphic compute-in-memory platform and showcased edge computing by implementing a spiking neural network for an autonomous navigation/racing task. Our work addresses challenges posed by existing RRAM technologies and paves the way for neuromorphic computing at the edge under strict size, weight, and power constraints.

97 MATHEMATICS AND COMPUTING↗

Aeroelastic loads prediction for an arrow wing. Task 3: Evaluation of the Boeing three-dimensional leading-edge vortex code

Two separated flow computer programs and a semiempirical method for incorporating the experimentally measured separated flow effects into a linear aeroelastic analysis were evaluated. The three dimensional leading edge vortex (LEV) code is evaluated. This code is an improved panel method for three dimensional inviscid flow over a wing with leading edge vortex separation. The governing equations are the linear flow differential equation with nonlinear boundary conditions. The solution is iterative; the position as well as the strength of the vortex is determined. Cases for both full and partial span vortices were executed. The predicted pressures are good and adequately reflect changes in configuration.

Manro, M. E.↗

Vibration control experiment design for the 15-m hoop/column antenna

A test program is designed for a ground-based vibration control experiment utilizing as the test article the 15-M Hoop/Column Antenna. Overall objectives of the designed ground-based test program include: (1) the validation of large space structure (LSS) control systemm techniques; (2) the validation of LSS parameter identification techniques: (3) the evaluation of actuator of actuator and sensor placement methodology; and (3) the validation of LSS computer models. Critical concerns in LSS Controls and Dynamics are: low frequency vibrational modes, close modal spacing, parameter uncertainties, controller software limitations, nonlinearities and coupling of modes through damping. Analytical results are presented which include compensator designs for varying compensator order.

Ham, F. M.↗

Horizontal vestibuloocular reflex evoked by high-acceleration rotations in the squirrel monkey. III. Responses after labyrinthectomy

The horizontal angular vestibuloocular reflex (VOR) evoked by high-frequency, high-acceleration rotations was studied in four squirrel monkeys after unilateral labyrinthectomy. Spontaneous nystagmus was measured at the beginning and end of each testing session. During the period that animals were kept in darkness (4 days), the nystagmus at each of these times measured approximately 20 degrees /s. Within 18-24 h after return to the light, the nystagmus (measured in darkness) decreased to 2.8 +/- 1.5 degrees /s (mean +/- SD) when recorded at the beginning but was 20.3 +/- 3.9 degrees /s at the end of the testing session. The latency of the VOR measured from responses to steps of acceleration (3,000 degrees /s(2) reaching a velocity of 150 degrees /s) was 8.4 +/- 0.3 ms for responses to ipsilesional rotations and 7.7 +/- 0.4 ms for contralesional rotations. During the period that animals were kept in darkness after the labyrinthectomy, the gain of the VOR measured during the steps of acceleration was 0.67 +/- 0.12 for contralesional rotations and 0.39 +/- 0.04 for ipsilesional rotations. Within 18-24 h after return to light, the VOR gain for contralesional rotations increased to 0.87 +/- 0.08, whereas there was only a slight increase for ipsilesional rotations to 0.41 +/- 0. 06. A symmetrical increase in the gain measured at the plateau of head velocity was noted after the animals were returned to light. The VOR evoked by sinusoidal rotations of 2-15 Hz, +/-20 degrees /s, showed a better recovery of gain at lower (2-4 Hz) than at higher (6-15 Hz) frequencies. At 0.5 Hz, gain decreased symmetrically when the peak amplitude was increased from 20 to 100 degrees /s. At 10 Hz, gain was decreased for ipsilesional half-cycles and increased for contralesional half-cycles when velocity was raised from 20 to 50 degrees /s. A model incorporating linear and nonlinear pathways was used to simulate the data. Selective increases in the gain for the linear pathway accounted for the recovery in VOR gain for responses at the velocity plateau of the steps of acceleration and for the sinusoidal rotations at lower peak velocities. The increase in gain for contralesional responses to steps of acceleration and sinusoidal rotations at higher frequencies and velocities was due to an increase in the contribution of the nonlinear pathway. This pathway was driven into cutoff and therefore did not affect responses for rotations toward the lesioned side.

NASA Discipline Neuroscience↗

Multiplexed Predictive Control of a Large Commercial Turbofan Engine

Model predictive control is a strategy well-suited to handle the highly complex, nonlinear, uncertain, and constrained dynamics involved in aircraft engine control problems. However, it has thus far been infeasible to implement model predictive control in engine control applications, because of the combination of model complexity and the time allotted for the control update calculation. In this paper, a multiplexed implementation is proposed that dramatically reduces the computational burden of the quadratic programming optimization that must be solved online as part of the model-predictive-control algorithm. Actuator updates are calculated sequentially and cyclically in a multiplexed implementation, as opposed to the simultaneous optimization taking place in conventional model predictive control. Theoretical aspects are discussed based on a nominal model, and actual computational savings are demonstrated using a realistic commercial engine model.

Richter, hanz↗

Aeroelasticity Benchmark Assessment: Subsonic Fixed Wing Program

The fundamental technical challenge in computational aeroelasticity is the accurate prediction of unsteady aerodynamic phenomena and the effect on the aeroelastic response of a vehicle. Currently, a benchmarking standard for use in validating the accuracy of computational aeroelasticity codes does not exist. Many aeroelastic data sets have been obtained in wind-tunnel and flight testing throughout the world; however, none have been globally presented or accepted as an ideal data set. There are numerous reasons for this. One reason is that often, such aeroelastic data sets focus on the aeroelastic phenomena alone (flutter, for example) and do not contain associated information such as unsteady pressures and time-correlated structural dynamic deflections. Other available data sets focus solely on the unsteady pressures and do not address the aeroelastic phenomena. Other discrepancies can include omission of relevant data, such as flutter frequency and / or the acquisition of only qualitative deflection data. In addition to these content deficiencies, all of the available data sets present both experimental and computational technical challenges. Experimental issues include facility influences, nonlinearities beyond those being modeled, and data processing. From the computational perspective, technical challenges include modeling geometric complexities, coupling between the flow and the structure, grid issues, and boundary conditions. The Aeroelasticity Benchmark Assessment task seeks to examine the existing potential experimental data sets and ultimately choose the one that is viewed as the most suitable for computational benchmarking. An initial computational evaluation of that configuration will then be performed using the Langley-developed computational fluid dynamics (CFD) software FUN3D1 as part of its code validation process. In addition to the benchmarking activity, this task also includes an examination of future research directions. Researchers within the Aeroelasticity Branch will examine other experimental efforts within the Subsonic Fixed Wing (SFW) program (such as testing of the NASA Common Research Model (CRM)) and other NASA programs and assess aeroelasticity issues and research topics.

Florance, Jennifer P.↗

Extensional, bending and twisting stiffness of Titanium Multiwall Thermal Protection System (TPS)

A test program which determines the extensional, bending and torsional stiffness of various titanium multiwall sandwich configurations is described. It is shown that unlike honeycomb core, the dimpled core is a significant contributor to the stiffness and strength of the sandwich. the extensional stiffness test shows irregularities which are attributed to foil thickness variations and to the difficulty to determine linear values from nonlinear tests.

Meaney, J. E.↗

Stacking sequence optimization of simply supported laminates with stability and strain constraints

An integer programming formulation for the design of symmetric and balanced rectangular composite laminates with simply supported boundary conditions subject to buckling and strain constraints is presented. The design variables that define the stacking sequence of the laminate are ply-identity zero-one integers. The buckling constraint is linear in terms of the ply-identity design variables, but strains are nonlinear functions of these variables. A linear approximation is developed for the strain constraints so that the problem can be solved by sequential linearization using the branch and bound algorithm. Examples of graphite-epoxy plates under biaxial compression are presented. Optimum stacking sequences obtained using the linear approximation are compared with global optimum designs obtained using a genetic search procedure.

Nagendra, S.↗

Modeling and control for closed environment plant production systems

A computer program was developed to study multiple crop production and control in controlled environment plant production systems. The program simulates crop growth and development under nominal and off-nominal environments. Time-series crop models for wheat (Triticum aestivum), soybean (Glycine max), and white potato (Solanum tuberosum) are integrated with a model-based predictive controller. The controller evaluates and compensates for effects of environmental disturbances on crop production scheduling. The crop models consist of a set of nonlinear polynomial equations, six for each crop, developed using multivariate polynomial regression (MPR). Simulated data from DSSAT crop models, previously modified for crop production in controlled environments with hydroponics under elevated atmospheric carbon dioxide concentration, were used for the MPR fitting. The model-based predictive controller adjusts light intensity, air temperature, and carbon dioxide concentration set points in response to environmental perturbations. Control signals are determined from minimization of a cost function, which is based on the weighted control effort and squared-error between the system response and desired reference signal.

NASA Discipline Life Support Systems↗

An Investigation of Nonlinear Controller for Propulsion Controlled Aircraft

Aircraft control systems are usually very reliable because of redundancy and multiple control surfaces. However, there are rare occasions when potentially disastrous flight control system failures do occur. At such times, the use of appropriate modulation of engine thrust to stabilize the aircraft may be the only chance of survival for the people aboard. In several cases where complete loss of control systems has occurred in multi-engine aircraft, pilots used the propulsion system to regain limited control of the aircraft with various degrees of success. In order to evaluate the feasibility of using only engine thrust modulation for emergency backup flight control, the NASA Dryden Flight Research Center has been conducting a series of analytical studies and flight tests on several different types of aircraft in a propulsion controlled aircraft (PCA) program. Simulation studies have included B-720, B-727, MD-11, C-402, C-17, F-18, and F-15, and flight tests have included B-747, B-777, MD-11, T-39, Lear 24, F-18, F-15, T-38, and PA-30. One objective was to determine the degree of control available with manual manipulation (open-loop) of the engine throttles. Flight tests and simulations soon showed that a closed loop controller could improve the chances of making a safe runway landing. The major work to date has concentrated on three aircraft (F-15, F-18, and the MD-11). Successful landings using PCA controllers were performed on the F-15 and MD-11 without the use of control surfaces. During the course of the research, some unique challenges have been identified. Compared to the conventional flight control surfaces, the engines are slow and have limited control effectiveness. Hence the ability of the system to promptly respond to aerodynamic changes is limited. Consequently, many nonlinear effects, which are easily accommodated by a conventional flight control system, become significant issues in the design of an effective controller when the engines are used as the only control effectors. A number of nonlinear behaviors observed during PCA flight tests, which are not meant to be exhaustive, are reviewed in the next section.

Lu, Ping↗

A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments

The computational study of chemical reactions in complex, wet environments is critical for applications in many fields. It is often essential to study chemical reactions in the presence of applied electrochemical potentials, taking into account the non-trivial electrostatic screening coming from the solvent and the electrolytes. As a consequence, the electrostatic potential has to be found by solving the generalized Poisson and the Poisson-Boltzmann equations for neutral and ionic solutions, respectively. In the present work, solvers for both problems have been developed. A preconditioned conjugate gradient method has been implemented for the solution of the generalized Poisson equation and the linear regime of the Poisson-Boltzmann, allowing to solve iteratively the minimization problem with some ten iterations of the ordinary Poisson equation solver. In addition, a self-consistent procedure enables us to solve the non-linear Poisson-Boltzmann problem. Both solvers exhibit very high accuracy and parallel efficiency and allow for the treatment of periodic, free, and slab boundary conditions. The solver has been integrated into the BigDFT and Quantum-ESPRESSO electronic-structure packages and will be released as an independent program, suitable for integration in other codes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Bilevel Voltage Regulation Operation for Distribution Systems With Self-Operated Microgrids

The emerging of microgrids in distribution systems has significantly enhanced the resilience of power grids. However, the operators of a distribution system and microgrids therein can be different and have accessibility to different devices. To model the operation of such a grid, this work proposes a bilevel formulation and probes into the voltage regulation operation, considering the interaction between different systems. The proposed bilevel formulation considers the cooperation of active energy resources (AER), transformer tap-changers, and capacitor banks that are controlled by different operators. To facilitate the solution time of the target bilevel optimization, the lower-level problems with different objectives are modeled using deep neural networks (DNNs) which are then converted into a set of constraints. Hence, the bilevel problem can be reformed to a single-level problem. Lastly, the proposed solution procedures are validated using a customized joint system constructed by the IEEE 123-bus system and a real distribution system in Iowa. According to the numerical validation results, the solution time of the proposed nonlinear activation function based DNN model is 69 times faster than other methods in solving voltage regulation with a bilevel structure.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A nonlinear propulsion system simulation technique for piloted simulators

A real time digital simulation technique providing the capabilities needed to evaluate propulsion system performance and aircraft system interaction on NASA manned flight simulators, is discussed. A parameter correlation technique is used with real and pseudo dynamics in a stable integration convergence loop. The cycle time reported was 2.0 ms on one computer and 5.7 ms on the simulator computer. The model was found to be stable and accurate with time up to 50 ms. It is concluded that the program has generated a valuable simulation technology and flight simulator experience by providing an adequate level of detail to evaluate propulsion systems in a simulated flight environment.

Mihaloew, J. R.↗

Nonlinear potential analysis techniques for supersonic aerodynamic design

A numerical method based on the conservation form of the full potential equation has been applied to the problem of three-dimensional supersonic flows with embedded subsonic regions. The governing equation is cast in a nonorthogonal coordinate system, and the theory of characteristics is used to accurately monitor the type-dependent flow field. A conservative switching scheme is employed to transition from the supersonic marching procedure to a subsonic relaxation algorithm and vice versa. The newly developed computer program can handle arbitrary geometries with fuselage, canard, wing, flow through nacelle, vertical tail and wake components at combined angles of attack and sideslip. Results are obtained for a variety of configurations that include a Langley advanced fighter concept with fuselage centerline nacelle, Rockwell's Advanced Tactical Fighter (ATF) with wing mounted nacelles, and the Shuttle Orbiter configuration. Comparisons with available experiments were good.

Shankar, V.↗

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Extracting the turbulent flow-field from beam emission spectroscopy images using velocimetry

The 2D turbulent E × B flow-field is inferred from density fluctuation images obtained with the beam emission spectroscopy diagnostic on DIII-D using the orthogonal dynamic programming velocimetry algorithm. A synthetic turbulence model is used to test the algorithm and optimize it for measuring zonal flows. Zonal flow measurements are found to require a signal-to-noise ratio above ~10 and a zonal flow wavelength longer than ~2 cm. Comparison between the velocimetry-estimated flow-field and the E × B flow-field using a nonlinear gyrokinetic GENE simulation finds that the flow-fields have identical spatial structure and differ only by the mean turbulence phase velocity, which is spatially uniform in this flux tube simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A numerical calculation of outward propagation of solar disturbances

The responses of the solar atmosphere due to an outward propagation shock are examined by employing the Lax-Wendroff method to solve the set of nonlinear partial differential equations in the model of the solar atmosphere. It is found that this theoretical model can be used to explain the solar phenomena of surge and spray. A criterion to discriminate the surge and spray is established and detailed information concerning the density, velocity, and temperature distribution with respect to the height and time is presented. The complete computer program is also included.

Wu, S. T.↗