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,099 records · Page 61

An exploratory study of finite difference grids for transonic unsteady aerodynamics

Unsteady aerodynamic forces are calculated by the XTRAN2L finite difference program which solves the complete two dimensional unsteady transonic small perturbation equation. The unsteady forces are obtained using a pulse transfer function technique which assumes the flow field behaves in a locally linear fashion about a mean condition. Forces are calculated for a linear flat plate using the default grids from the LTRAN2-NLR, LTRAN2-HI, and XTRAN3S programs. The forces are compared to the exact theoretical values for flat plate, and grid generated boundary and internal numerical reflections are observed to cause significant errors in the unsteady airloads. Grids are presented that alleviate the reflections while reducing computational time up to fifty-three percent and program size up to twenty-eight percent. Forces are presented for a six percent thick parabolic arc airfoil which demonstrate that the transform technique may be successfully applied to nonlinear transonic flows.

Seidel, D. A.↗

A predictive model for failure properties of thermoset resins

A predictive model for the three-dimensional failure behavior of engineering polymers has been developed in a recent NASA-sponsored research program. This model acknowledges the underlying molecular deformation mechanisms and thus accounts for the effects of different chemical compositions, crosslink density, functionality of the curing agent, etc., on the complete nonlinear stress-strain response including yield. The material parameters required by the model can be determined from test-tube quantities of a new resin in only a few days. Thus, we can obtain a first-order prediction of the applicability of a new resin for an advanced aerospace application without synthesizing the large quantities of material needed for failure testing. This technology will effect order-of-magnitude reductions in the time and expense required to develop new engineering polymers.

Caruthers, James M.↗

A vectorized algorithm for 3D dynamics of a tethered satellite

Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.

Wilson, Howard B.↗

Structural reliability methods: Code development status

The Probabilistic Structures Analysis Method (PSAM) program integrates state of the art probabilistic algorithms with structural analysis methods in order to quantify the behavior of Space Shuttle Main Engine structures subject to uncertain loadings, boundary conditions, material parameters, and geometric conditions. An advanced, efficient probabilistic structural analysis software program, NESSUS (Numerical Evaluation of Stochastic Structures Under Stress) was developed as a deliverable. NESSUS contains a number of integrated software components to perform probabilistic analysis of complex structures. A nonlinear finite element module NESSUS/FEM is used to model the structure and obtain structural sensitivities. Some of the capabilities of NESSUS/FEM are shown. A Fast Probability Integration module NESSUS/FPI estimates the probability given the structural sensitivities. A driver module, PFEM, couples the FEM and FPI. NESSUS, version 5.0, addresses component reliability, resistance, and risk.

Millwater, Harry R.↗

Semi-analytical static nonlinear structural sensitivity analysis

The paper discusses two uses of the semianalytical methods for nonlinear sensitivity analysis. First, the application to noncritical response is analyzed. It is shown that the semianalytical method is equivalent to a particular form of the overall finite difference approach. It is concluded that the difference between overall finite differences and the semianalytical approach is blurred for the case of nonlinear static response. Next, the application of the semianalytical method to calculating sensitivities of limit loads is discussed. A formulation that is easy to implement in general purpose finite element programs is derived. Three examples are used to demonstrate the application of the formulation and to explore its accuracy.

Haftka, Raphael T.↗

Optimal System Realization in Frequency Domain

Several approaches are presented to identify an experimental system model directly from frequency response data. The formulation begins with a matrix-fraction description as the model structure. Frequency weighting such as exponential weighting is introduced to solve a weighted least-squares problem to obtain the coefficient matrices for the matrix-fraction description. A multi-variable state-space model can then be formed using the coefficient matrices of the matrix-fraction description. An approach is introduced to fine-tune the model using non-linear programming methods to minimize the desired cost function. The method deals with the model in the real Schur or modal form and reassigns a subset of system poles using a nonlinear optimizer. At every optimization step, the input and output influence matrices are refined through least-squares procedures. The proposed approaches are used to identify an analytical model for a NASA testbed from experimental data.

Juang, Jer-Nan↗

Near-Simultaneous Spectroscopic and Broadband Polarimetric Observations of Be Stars

Near simultaneous optical spectroscopic (on four nights) and broadband linear continuum (B, V, R, and I bands) polarimetric (on seven nights) observations of 29 Be stars were carried out during 1993 November-December. The program Be stars displayed wavelength dependence of intrinsic polarizations with no frequency dependence of polarimetric position angles. Some of the Be stars displayed long-term polarization variability. The Be and Be-shell stars could not be distinguished from one another solely on the basis of their polarization values. Full widths at half-maximum of the H.alpha profiles and the intrinsic linear continuum polarizations are closely correlated with the projected rotational velocities of the program stars. Photospheric-absorption-corrected equivalent widths of H.alpha profiles [W(alpha)] and the radii of H.alpha-emitting or -absorbing envelopes (R(sub e) or R(sub a)) are nonlinearly correlated with the intrinsic continuum polarizations of these stars. However, W(alpha) and R(sub e) are linearly correlated. With large uncertainties, there is a trend of spectral dependence of polarization. Detailed discussion of these results is presented in this paper.

Ghosh, K.↗

Dynamics Impact Tolerance of Shuttle RCC Leading Edge Panels Using LS-DYNA

This paper describes a research program conducted to enable accurate prediction of the impact tolerance of the shuttle Orbiter leading-edge wing panels using physics-based codes such as LS-DYNA, a nonlinear, explicit transient dynamic finite element code. The shuttle leading-edge panels are constructed of Reinforced-Carbon-Carbon (RCC) composite material, which is used because of its thermal properties to protect the shuttle during reentry into the Earth's atmosphere. Accurate predictions of impact damage from insulating foam and other debris strikes that occur during launch required materials characterization of expected debris, including strain-rate effects. First, analytical models of individual foam and RCC materials were validated. Next, analytical models of foam cylinders impacting 6- in. x 6-in. RCC flat plates were developed and validated. LS-DYNA pre-test models of the RCC flat plate specimens established the impact velocity of the test for three damage levels: no-detectable damage, non-destructive evaluation (NDE) detectable damage, or visible damage such as a through crack or hole. Finally, the threshold of impact damage for RCC on representative Orbiter wing panels was predicted for both a small through crack and for NDE-detectable damage.

Fasanella, Edwin L.↗

Dynamic Impact Tolerance of Shuttle RCC Leading Edge Panels using LS-DYNA

This paper describes a research program conducted to enable accurate prediction of the impact tolerance of the shuttle Orbiter leading-edge wing panels using 'physics-based- codes such as LS-DYNA, a nonlinear, explicit transient dynamic finite element code. The shuttle leading-edge panels are constructed of Reinforced-Carbon-Carbon (RCC) composite material, which issued because of its thermal properties to protect the shuttle during re-entry into the Earth's atmosphere. Accurate predictions of impact damage from insulating foam and other debris strikes that occur during launch required materials characterization of expected debris, including strain-rate effects. First, analytical models of individual foam and RCC materials were validated. Next, analytical models of individual foam cylinders impacting 6-in. x 6-in. RCC flat plates were developed and validated. LS-DYNA pre-test models of the RCC flat plate specimens established the impact velocity of the test for three damage levels: no-detectable damage, non-destructive evaluation (NDE) detectable damage, or visible damage such as a through crack or hole. Finally, the threshold of impact damage for RCC on representative Orbiter wing panels was predicted for both a small through crack and for NDE-detectable damage.

Fasanella, Edwin↗

Applying cumulative effects to strategically advance large-scale ecosystem restoration

Programs restoring large terrestrial-aquatic ecosystems are expected to increase in area and number over the next decade. Factors ensuring long-term effectiveness at the program scale differ from those at the individual restoration project scale, necessitating a new synthesis. Here we contextualize hydrological, geomorphological, and biological effects of large-scale ecosystem restoration activities within a cumulative-effects paradigm. Evidence of eight modes of cumulative effects benefiting species and ecosystems was elucidated by reviewing restoration efforts in the greater Florida Everglades, Gulf Coast, lower Columbia River and estuary, Puget Sound estuary, Missouri River, northeastern coastal states, San Francisco Bay and Sacramento Delta. Additive and nonlinear cumulative effects of interacting restoration projects propagated across landscape and regional geographic scales, and throughout trophic levels, affecting population structures. Applying a cumulative-effects paradigm to planning, implementing, and evaluating large-scale restoration will avoid countervailing effects and capitalize on positive feedback, advancing evidence-based programs to recover priority species and ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Fast model-based scenario optimization in NSTX-U enabled by analytic gradient computation

Model-based optimization offers a systematic approach to advanced scenario planning. In this case, the feedforward-control inputs (actuator trajectories) that are needed to attain and sustain a desired scenario are obtained by solving a nonlinear constrained optimization problem. This class of problems generally minimize a cost function that measures the difference between desired and actual plasma states. Several numerical optimization algorithms, such as sequential quadratic programming, require repeated calculation of the cost function gradients with respect to the input trajectories. Calculating these gradients numerically can be computationally intensive, increasing the time needed to solve the feedforward-control optimization problem. Here, this work introduces a method to analytically calculate these cost function gradients from the current profile evolution model. This can significantly reduce the computational time and allow for fast feedforward-control optimization, which would eventually enable optimal scenario planning between discharges. The performance of the feedforward optimizer with analytical gradients is compared to a traditional optimization algorithm based on numerical gradients for different NSTX-U scenarios. The plasma dynamics in the optimization algorithm are simulated using the Control Oriented Transport SIMulator (COTSIM). Results of the work show that analytical gradients consistently reduce the computation time while achieving trajectories that are comparable to those obtained by traditional optimization algorithms based on numerical gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toward performance-portable PETSc for GPU-based exascale systems

The Portable Extensible Toolkit for Scientific computation (PETSc) library delivers scalable solvers for nonlinear time-dependent differential and algebraic equations and for numerical optimization. The PETSc design for performance portability addresses fundamental GPU accelerator challenges and stresses flexibility and extensibility by separating the programming model used by the application from that used by the library, and it enables application developers to use their preferred programming model, such as Kokkos, RAJA, SYCL, HIP, CUDA, or OpenCL, on upcoming exascale systems. Furthermore, a blueprint for using GPUs from PETSc-based codes is provided, and case studies emphasize the flexibility and high performance achieved on current GPU-based systems.

97 MATHEMATICS AND COMPUTING↗

Shock and vibration technology with applications to electrical systems

A survey is presented of shock and vibration technology for electrical systems developed by the aerospace programs. The shock environment is surveyed along with new techniques for modeling, computer simulation, damping, and response analysis. Design techniques based on the use of analog computers, shock spectra, optimization, and nonlinear isolation are discussed. Shock mounting of rotors for performance and survival, and vibration isolation techniques are reviewed.

Eshleman, R. L.↗

Leading-edge force features of the aerodynamic finite element method.

Description of a practical procedure for computing the wing leading-edge thrust distribution by the finite element method. When incorporated into a wing-body aerodynamic computer program, the technique is capable of predicting (at subsonic and supersonic speeds) the leading-edge thrust distribution (and therefore, the lateral-directional stability derivatives due to roll) and the nonlinear aerodynamic characteristics of low aspect-ratio wings with leading edge separation due to the application of suction technology.

Lan, C.-T.↗

Buckling behavior and structural efficiency of open-section stiffened composite compression panels

Several experiments with J- and blade-stiffened graphite/epoxy panels were conducted to obtain insight into how well experimental data could be correlated with analysis for the buckling behavior of open-section stiffened composite compression panels. Although some nonlinear behavior was observed during the experiments, adequate correlation with analysis was obtained to justify the use of linear, thin-plate buckling analysis in a minimum-weight design synthesis program for J- and blade-configurations. Results from two design studies using this program are presented. In the first study the minimum weights of Jand blade-configurations for two different material systems (graphite/epoxy and aluminum) are determined subject to buckling and strength constraints for a wide range of the compressive load index. In the second study the minimum weights required for graphite/epoxy blade-stiffened panels to satisfy additional stiffness constraints typical of medium-size commercial aircraft wing structures are determined. Both minimum-weight studies indicate that graphite/epoxy open-section stiffened panels can be designed so that weight savings of 30 to 50% are possible compared with the most efficient aluminum designs.

Williams, J. G.↗

A study of coupled rotor-fuselage vibration with higher harmonic control using a symbolic computing facility

A fundamental study of vibration prediction and vibration reduction in helicopters using active controls was performed. The nonlinear equations of motion for a coupled rotor/flexible fuselage system have been derived using computer algebra on a special purpose symbolic computing facility. The details of the derivation using the MACSYMA program are described. The trim state and vibratory response of the helicopter are obtained in a single pass by applying the harmonic balance technique and simultaneously satisfying the trim and the vibratory response of the helicopter for all rotor and fuselage degrees of freedom. The influence of the fuselage flexibility on the vibratory response is studied. It is shown that the conventional single frequency higher harmonic control (HHC) capable of reducing either the hub loads or only the fuselage vibrations but not both simultaneously. It is demonstrated that for simultaneous reduction of hub shears and fuselage vibrations a new scheme called multiple higher harmonic control (MHHC) is required. The fundamental aspects of this scheme and its uniqueness are described in detail, providing new insight on vibration reduction in helicopters using HHC.

Papavassiliou, I.↗

Implementation and Testing of Turbulence Models for the F18-HARV Simulation

This report presents three methods of implementing the Dryden power spectral density model for atmospheric turbulence. Included are the equations which define the three methods and computer source code written in Advanced Continuous Simulation Language to implement the equations. Time-history plots and sample statistics of simulated turbulence results from executing the code in a test program are also presented. Power spectral densities were computed for sample sequences of turbulence and are plotted for comparison with the Dryden spectra. The three model implementations were installed in a nonlinear six-degree-of-freedom simulation of the High Alpha Research Vehicle airplane. Aircraft simulation responses to turbulence generated with the three implementations are presented as plots.

Yeager, Jessie C.↗

Application of a Split-Fiber Probe to Velocity Measurement in the NASA Research Compressor

A split-fiber probe was used to acquire unsteady data in a research compressor. The probe has two thin films deposited on a quartz cylinder 200 microns in diameter. A split-fiber probe allows simultaneous measurement of velocity magnitude and direction in a plane that is perpendicular to the sensing cylinder, because it has its circumference divided into two independent parts. Local heat transfer considerations indicated that the probe direction characteristic is linear in the range of flow incidence angles of +/- 35. Calibration tests confirmed this assumption. Of course, the velocity characteristic is nonlinear as is typical in thermal anemometry. The probe was used extensively in the NASA Glenn Research Center (GRC) low-speed, multistage axial compressor, and worked reliably during a test program of several months duration. The velocity and direction characteristics of the probe showed only minute changes during the entire test program. An algorithm was developed to decompose the probe signals into velocity magnitude and velocity direction. The averaged unsteady data were compared with data acquired by pneumatic probes. An overall excellent agreement between the averaged data acquired by a split-fiber probe and a pneumatic probe boosts confidence in the reliability of the unsteady content of the split-fiber probe data. To investigate the features of unsteady data, two methods were used: ensemble averaging and frequency analysis. The velocity distribution in a rotor blade passage was retrieved using the ensemble averaging method. Frequencies of excitation forces that may contribute to high cycle fatigue problems were identified by applying a fast Fourier transform to the absolute velocity data.

Lepicovsky, Jan↗