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,477 records · Page 82

Development of the reentry flight dynamics simulator for evaluation of space shuttle orbiter entry systems

A nonlinear, six degree of freedom, digital computer simulation of a vehicle which has constant mass properties and whose attitudes are controlled by both aerodynamic surfaces and reaction control system thrusters was developed. A rotating, oblate Earth model was used to describe the gravitational forces which affect long duration Earth entry trajectories. The program is executed in a nonreal time mode or connected to a simulation cockpit to conduct piloted and autopilot studies. The program guidance and control software used by the space shuttle orbiter for its descent from approximately 121.9 km to touchdown on the runway.

Rowell, L. F.↗

The application of TSIM software to act design and analysis on flexible aircraft

The TSIM software is described. This is a package which uses an interactive FORTRAN-like simulation language for the simulation on nonlinear dynamic systems and offers facilities which include: mixed continuous and discrete time systems, time response calculations, numerical optimization, automatic trimming of nonlinear aircraft systems, and linearization of nonlinear equations for eigenvalues, frequency responses and power spectral response evaluation. Details are given of the application of TSIM to the analysis of aeroelastic systems under the RAE Farborough extension FLEX-SIM. The aerodynamic and structural data for the equations of motion of a flexible aircraft are prepared by a preprocessor program for incorporation in TSIM simulations. Within the simulation, the flexible aircraft model may then be selected interactively for different flight conditions and modal reduction techniques applied. The use of FLEX-SIM is demonstrated by an example of the flutter prediction for a simple aeroelastic model. By utilizing the numerical optimization facility of TSIM, it is possible to undertake identification of required parameters in the TSIM model within the simulation.

Kaynes, Ian W.↗

A Comparison of Trajectory Optimization Methods for the Impulsive Minimum Fuel Rendezvous Problem

In this paper we present a comparison of optimization approaches to the minimum fuel rendezvous problem. Both indirect and direct methods are compared for a variety of test cases. The indirect approach is based on primer vector theory. The direct approaches are implemented numerically and include Sequential Quadratic Programming (SQP), Quasi-Newton, Simplex, Genetic Algorithms, and Simulated Annealing. Each method is applied to a variety of test cases including, circular to circular coplanar orbits, LEO to GEO, and orbit phasing in highly elliptic orbits. We also compare different constrained optimization routines on complex orbit rendezvous problems with complicated, highly nonlinear constraints.

Hughes, Steven P.↗

Phonon-informed Neural Thermal Scattering (NeTS) Optimization for Crystalline Graphite and Beryllium Metal

Fast neutrons born from fission lose energy through scattering interactions in the process of slowing-down. As neutrons thermalize to the order of $k$ $b$ $T$ (where $k$ $b$ is the Boltzmann constant, and $T$ is the temperature of the medium), their de Broglie wavelength and energy approaches the order of inter-atomic spacing and quantized lattice vibrations, i.e., phonons. At thermal energies, the thermal scattering law (TSL), i.e., $S$($α, β$), captures crystal binding contributions to the total reaction rate, or cross section. This dimensionless material property describes the energy ($β$) and momentum ($α$) exchanges available in a medium. Currently, $S$($α, β$) is evaluated in the Full Law Analysis Scattering System Hub (FLASSH) code for discrete inputs and stored as ENDF/B File 7 for 0-phonon elastic (MT 2) and n-phonon inelastic (MT 4) processes. Further processing recasts $S$($α, β$) into cumulative distribution functions for sampling post-collision scattering kinematics. In practice, interpolation schemes are employed to access data between tabulated values. An improvement to this juncture of the nuclear data pipeline is supplying cross sections on-the-fly (OTF), as has been developed for the un-resolved resonance region to minimize non-physical interpolation errors. This capability may improve simulation accuracy for accident and transient analyses, where rapidly varying changes in temperature and pressure are difficult to predict beforehand. To do so, deep artificial neural networks (ANNs) can be employed which collapse non-linear, complex data into a lightweight dictionary of neural weights and biases. This has been successfully demonstrated for the hydrogen in light water $S$($α, β$) dataset in the form of a Neural Thermal Scattering (NeTS) module. In this work, the NeTS framework is extended to consider the impact of material-dependent dynamical features on optimal neural pre-processing and architecture design decisions, such as number of neurons per hidden layer, residual skip connections and neural depth. New NeTS modules for crystalline graphite and beryllium metal illuminate a novel correlation between dynamical nonlinearity and optimal neural parametrization when deploying $S$($α, β$) on-the-fly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Overview of the New Beam Physics Research at the IOTA/FAST Facility

The Fermilab Accelerator Science and Technology (FAST) facility is dedicated to the exploration of novel concepts in accelerator and beam physics, and the development of a robust workforce, in order to enable and enhance next-generation particle accelerators. FAST comprises a high-brightness superconducting electron linac, and a storage ring, the Integrable Optics Test Accelerator (IOTA). Experiments in the most recent operational run include studies of nonlinear integrable lattices; tracking of single electrons; precise characterization of undulator radiation; studies with low-momentum-compaction lattices; and ultra-wide range beam diagnostics based on Photomultiplier tubes. In the linac, experiments on noise in intense electron bunches were conducted. The IOTA proton injector, currently being commissioned, will enable a diverse program on space-charge-dominated beams. Research areas include non-invasive beam profile monitoring for proton beams; beam dynamics with electron lenses; halo suppression, feedback systems, and electron cooling. In this presentation, we provide an overview of the recent results and highlight future plans together with opportunities for collaboration.

43 PARTICLE ACCELERATORS↗

Sonic-boom-induced building structure responses including damage.

Concepts of sonic-boom pressure loading of building structures and the associated responses are reviewed, and results of pertinent theoretical and experimental research programs are summarized. The significance of sonic-boom load time histories, including waveshape effects, are illustrated with the aid of simple structural elements such as beams and plates. Also included are discussions of the significance of such other phenomena as three-dimensional loading effects, air cavity coupling, multimodal responses, and structural nonlinearities. Measured deflection, acceleration, and strain data from laboratory models and full-scale building tests are summarized, and these data are compared, where possible, with predicted values. Damage complaint and claim experience due both to controlled and uncontrolled supersonic flights over communities are summarized with particular reference to residential, commercial, and historic buildings. Sonic-boom-induced building responses are compared with those from other impulsive loadings due to natural and cultural events and from laboratory simulation tests.

Clarkson, B. L.↗

A numerical procedure for predicting creep and delayed failures in laminated composites

A numerical procedure is described for predicting the viscoelastic response of general laminates. A nonlinear compliance model is used to predict the creep response of the individual laminae. A biaxial delayed failure model predicts ply failure. The numerical procedure, based on lamination theory, increases by increments through time to predict creep compliance and delayed failures in laminates. Numerical stability problems and experimental verification are discussed. Although the program has been quite successful in predicting creep of general laminates, the assumptions associated with lamination theory have resulted in erroneous bounds on the predicted material response. Delayed failure predictions have been conservative. Several improvements are suggested to increase the accuracy of the procedure.

Dillard, D. A.↗

System IDentification Programs for AirCraft (SIDPAC)

A collection of computer programs for aircraft system identification is described and demonstrated. The programs, collectively called System IDentification Programs for AirCraft, or SIDPAC, were developed in MATLAB as m-file functions. SIDPAC has been used successfully at NASA Langley Research Center with data from many different flight test programs and wind tunnel experiments. SIDPAC includes routines for experiment design, data conditioning, data compatibility analysis, model structure determination, equation-error and output-error parameter estimation in both the time and frequency domains, real-time and recursive parameter estimation, low order equivalent system identification, estimated parameter error calculation, linear and nonlinear simulation, plotting, and 3-D visualization. An overview of SIDPAC capabilities is provided, along with a demonstration of the use of SIDPAC with real flight test data from the NASA Glenn Twin Otter aircraft. The SIDPAC software is available without charge to U.S. citizens by request to the author, contingent on the requestor completing a NASA software usage agreement.

Morelli, Eugene A.↗

JuTrack: A Julia package for auto-differentiable accelerator modeling and particle tracking

Efficient accelerator modeling and particle tracking are key for the design and configuration of modern particle accelerators. In this work, we present JuTrack, a nested accelerator modeling package developed in the Julia programming language and enhanced with compiler-level automatic differentiation (AD). With the aid of AD, JuTrack enables rapid derivative calculations in accelerator modeling, facilitating sensitivity analyses and optimization tasks. Here we demonstrate the effectiveness of AD-derived derivatives through several practical applications, including sensitivity analysis of space-charge-induced emittance growth, nonlinear beam dynamics analysis for a synchrotron light source, and lattice parameter tuning of the future Electron-Ion Collider (EIC). Through the incorporation of automatic differentiation, this package opens up new possibilities for accelerator physicists in beam physics studies and accelerator design optimization.

43 PARTICLE ACCELERATORS↗

Numerical and Experimental Investigation of Turbulent Transport Control via Shaping of Radial Plasma Flow Profiles (Final Technical Report)

The strong coupling between the different physical variables involved in the plasma transport phenomenon and the high complexity of its dynamics call for a model-based, multivariable approach to profile control where those predictive models could be exploited. The overall objective of this project has been to extend the existing body of work by investigating numerically and experimentally active control of unstable fluctuations, including fully developed turbulence and the associated cross-field particle transport, via manipulation of flow profiles in a magnetized laboratory plasma device. Fluctuations and particle transport can be monitored by an array of electrostatic probes, and E x B flow profiles can be controlled via a set of biased concentric ring electrodes that terminate the plasma column. The goals of the proposed research have been threefold: i- to develop a predictive code to simulate plasma transport in the linear HELCAT (HELicon-CAThode) plasma device at the University of New Mexico (UNM), where the experimental component of the proposed research has been carried out; ii- to establish the feasibility of using advanced model-based control algorithms to control cross-field turbulence-driven particle transport through appropriate manipulation of radial plasma flow profiles, iii- to investigate the fundamental nonlinear dynamics of turbulence and transport physics. Lehigh University (LU), including Prof. Eugenio Schuster and one full-time graduate student, has been primarily responsible for control-oriented modeling and model-based control design. Undergraduate students have also participated in this project through the National Science Foundation Research Experience for Undergraduate (REU) program. The main goal of the LU Plasma Control Group has been to study the feasibility of controlling turbulence-driven transport by shaping the radial poloidal flow profile (i.e., by controlling flow shear) via biased concentric ring electrodes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computer program for analysis of imperfection sensitivity of ring stiffened shells of revolution

A FORTRAN 4 digital computer program is presented for the initial postbuckling and imperfection sensitivity analysis of bifurcation buckling modes for ring-stiffened orthotropic multilayered shells of revolution. The boundary value problem for the second-order contribution to the buckled state was solved by the forward integration technique using the Runge-Kutta method. The effects of nonlinear prebuckling states and live pressure loadings are included.

Cohen, G. A.↗

Computation of acoustic waves in a jet

A numerical treatment of acoustic waves in a jet is described. The full time dependent Euler equations are used in both linear and nonlinear formulations. The computational region of integration is artificially bounded and boundary conditions are developed to simulate outgoing waves and to enable the computational domain to be substantially restricted. Higher order methods and coordinate transformations are introduced to further reduce the number of grid points as well as to increase the efficiency of the program. Numerical results are presented for time harmonic sources as well as for sources with more complicated time dependence.

Bayliss, A.↗

Residual Strength Analysis Methodology: Laboratory Coupons to Structural Components

The NASA Aircraft Structural Integrity (NASIP) and Airframe Airworthiness Assurance/Aging Aircraft (AAA/AA) Programs have developed a residual strength prediction methodology for aircraft fuselage structures. This methodology has been experimentally verified for structures ranging from laboratory coupons up to full-scale structural components. The methodology uses the critical crack tip opening angle (CTOA) fracture criterion to characterize the fracture behavior and a material and a geometric nonlinear finite element shell analysis code to perform the structural analyses. The present paper presents the results of a study to evaluate the fracture behavior of 2024-T3 aluminum alloys with thickness of 0.04 inches to 0.09 inches. The critical CTOA and the corresponding plane strain core height necessary to simulate through-the-thickness effects at the crack tip in an otherwise plane stress analysis, were determined from small laboratory specimens. Using these parameters, the CTOA fracture criterion was used to predict the behavior of middle crack tension specimens that were up to 40 inches wide, flat panels with riveted stiffeners and multiple-site damage cracks, 18-inch diameter pressurized cylinders, and full scale curved stiffened panels subjected to internal pressure and mechanical loads.

Dawicke, D. S.↗

Orion Parachute Riser Cutter Development

This paper presents the tests and analytical approach used on the development of a steel riser cutter for the CEV Parachute Assembly System (CPAS) used on the Orion crew module. Figure 1 shows the riser cutter and the steel riser bundle which consists of six individual cables. Due to the highly compressed schedule, initial unavailability of the riser material and the Orion Forward Bay mechanical constraints, JSC primarily relied on a combination of internal ballistics analysis and LS-DYNA simulation for this project. Various one dimensional internal ballistics codes that use standard equation of state and conservation of energy have commonly used in the development of CAD devices for initial first order estimates and as an enhancement to the test program. While these codes are very accurate for propellant performance prediction, they usually lack a fully defined kinematic model for dynamic predictions. A simple piston device can easily and accurately be modeled using an equation of motion. However, the accuracy of analytical models is greatly reduced on more complicated devices with complex external loads, nonlinear trajectories or unique unlocking features. A 3D finite element model of CAD device with all critical features included can vastly improve the analytical ballistic predictions when it is used as a supplement to the ballistic code. During this project, LS-DYNA structural 3D model was used to predict the riser resisting load that was needed for the ballistic code. A Lagrangian model with eroding elements shown in Figure 2 was used for the blade, steel riser and the anvil. The riser material failure strain was fine tuned by matching the dent depth on the anvil with the actual test data. LS-DYNA model was also utilized to optimize the blade tip design for the most efficient cut. In parallel, the propellant type and the amount were determined by using CADPROG internal ballistics code. Initial test results showed a good match with LS-DYNA and CADPROG simulations. Final paper will present a detailed roadmap from initial ballistic modeling and LS-DYNA simulation to the performance testing. Blade shape optimization study will also be presented.

Oguz, Sirri↗

Programmable soft valves for digital and analog control

In soft devices, complex actuation sequences and precise force control typically require hard electronic valves and microcontrollers. Existing designs for entirely soft pneumatic control systems are capable of either digital or analog operation, but not both, and are limited by speed of actuation, range of pressure, time required for fabrication, or loss of power through pull-down resistors. Using the nonlinear mechanics intrinsic to structures composed of soft materials—in this case, by leveraging membrane inversion and tube kinking—two modular soft components are developed: a piston actuator and a bistable pneumatic switch. These two components combine to create valves capable of analog pressure regulation, simplified digital logic, controlled oscillation, nonvolatile memory storage, linear actuation, and interfacing with human users in both digital and analog formats. Three demonstrations showcase the capabilities of systems constructed from these valves: 1) a wearable glove capable of analog control of a soft artificial robotic hand based on input from a human user’s fingers, 2) a human-controlled cushion matrix designed for use in medical care, and 3) an untethered robot which travels a distance dynamically programmed at the time of operation to retrieve an object. This work illustrates pathways for complementary digital and analog control of soft robots using a unified valve design.

42 ENGINEERING↗

Studies on L-citrulline doped potassium dihydrogen phosphate- A non linear crystal with significant nonlinear properties

Potassium Dihydrogen Phosphate (KDP) single crystal is considered as one of the best representative of nonlinear optical crystals. Recently, amino acids having excellent nonlinear optical characteristics are being investigated as prospective dopants to improve the non linear optical characteristics of KDP. The present work is an attempt in this direction and L citrulline, one of the non essential amino acids showing good non linear optical characteristics is used as the dopant for KDP. Good quality crystals of L-citrulline doped KDP crystals were grown by slow evaporation technique. From the powder X-ray diffraction studies of doped KDP crystal, the structure of the doped crystals was determined by direct method and refined by Pawley method employing Topaz version program using the single crystal X-ray data for pure KDP. The lattice parameters for L citrulline doped KDP are a=7.467A{sup 0}, b=7.467 A{sup 0}, c=6.977 A{sup 0}. The crystal falls into the tetragonal crystal system with space group I42 d. The presence of carbon and oxygen, which are primary components of amino acids, in the EDAX spectrum confirms the effectiveness of doping. The absorption spectra of the doped samples show that the crystals are transparent in the entire visible region. The second harmonic generation efficiency of the doped samples was determined by Kurtz powder technique using the Q-switched Nd:YAG laser beam and is found to be 2.2 times that of KDP. The nonlinear optical properties can be well studied by the open aperture Z scan technique. The open aperture curve exhibits a normalized transmittance valley. The nonlinear absorption coefficient β is obtained by theoretical fitting for two photon absorption. It is inferred that doping KDP with L citrulline has enhanced the nonlinearity considerably. This obviously suggests the potentiality of the crystal as an optical power limiter and also for various optical device applications.

36 MATERIALS SCIENCE↗

Estimating epistemic uncertainty in soil parameters for nonlinear site response analyses: Introducing the Latin Hypercube Sampling technique

This study quantifies the effects of epistemic uncertainty in soil parameters on nonlinear (NL) site response analysis (SRA) results, validated against the data recorded at a well-instrumented geotechnical downhole array located in Japan. To this end, a one-dimensional soil column model of the Service Hall Array (SHA) near the Kashiwazaki-Kariwa Nuclear Power Plant (KKNPP) is developed using the finite element (FE) program LS-DYNA. The dynamic stress–strain relationship is characterized by a modified two-stage hyperbolic (MTH) NL backbone curve formulation capable of capturing soil behavior at both small- and large-shear strains. The model is then validated against the ground motion recordings to capture the model bias. The uncertainties associated with the shear-wave velocity profile (a small-strain soil property) and soil shear strength (a large-strain soil property) are incorporated in NL SRA to quantify their separate and joint randomization effects on the results. This study proposes using the Latin Hypercube Sampling (LHS) method as an efficient alternative to commonly used methods, such as Standard Monte Carlo (SMC), to account for uncertainty propagation in such reliability analysis. Both low-intensity and design-level records from the recordings at the SHA are applied to study the contribution of the small- and large-strain NL dynamic soil properties. Results from 46,200 NL FE analyses (23,100 per input motion) are presented. Measured and predicted site response, using recorded ground motions at this downhole array, is compared to assess the significance of soil parameter uncertainty on the observed ground motion dispersions. It is demonstrated that increasing the number of soft realizations and implementing higher level earthquake intensity lead to higher ground motion dispersion. Unlike past studies in randomization of Vs profiles with the SMC method, the LHS method is shown to have no significant effect on the predicted median surface response spectra and amplification factors (AFs) for this case study.

Engineering↗

Unsteady flow through compressor stages

The application of the nonsteady Lax-Wendroff technique to problems with asymptotically periodic solution which offers a potentially powerful method for the investigation of the interaction of rotating and stationary blade rows in turbomachinery is reported. A technique for specifying boundary conditions with phase lag was developed to accomplish this. A complete nonlinear analysis is carried out numerically to determine the entire flow field without recourse to the assumption of small disturbances of linear equations which underlie the previous acoustic theories. The result, obtained for the case of equal number of rotor and stator blades shows that transonic flow can be handled without difficulty. In addition, the program is not limited with regard to blade thickness, camber or loading. Extension of this method to incorporate viscous wakes and to analysis of fully three dimensional configuration is feasible, and would greatly expand its utility in practical applications.

Alzner, E.↗