Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear function”

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,243 records · Page 69

Computation of optimal output-feedback compensators for linear time-invariant systems

The control of linear time-invariant systems with respect to a quadratic performance criterion was considered, subject to the constraint that the control vector be a constant linear transformation of the output vector. The optimal feedback matrix, f*, was selected to optimize the expected performance, given the covariance of the initial state. It is first shown that the expected performance criterion can be expressed as the ratio of two multinomials in the element of f. This expression provides the basis for a feasible method of determining f* in the case of single-input single-output systems. A number of iterative algorithms are then proposed for the calculation of f* for multiple input-output systems. For two of these, monotone convergence is proved, but they involve the solution of nonlinear matrix equations at each iteration. Another is proposed involving the solution of Lyapunov equations at each iteration, and the gradual increase of the magnitude of a penalty function. Experience with this algorithm will be needed to determine whether or not it does, indeed, possess desirable convergence properties, and whether it can be used to determine the globally optimal f*.

Platzman, L. K.↗

Application of the MNA design method to a nonlinear turbofan engine

Using nonlinear digital simulation as a representative model of the dynamic operation of the QCSEE turbofan engine, a feedback control system is designed by variable frequency design techniques. Transfer functions are generated for each of five power level settings covering the range of operation from approach power to full throttle (62.5% to 100% full power). These transfer functions are then used by an interactive control system design synthesis program to provide a closed loop feedback control using the multivariable Nyquist array and extensions to multivariable Bode diagrams and Nichols charts.

Leininger, G. G.↗

A nonparametric identification method for a class of non-linear systems

A relatively simple nonparametric method for the identification of a class of close-coupled nonlinear multi-degree-of-freedom systems has been developed. The identification of arbitrary memoryless nonlinearities is possible through knowledge of the accelerations, velocities, and displacements of the various masses. These quantities are used to obtain the surfaces of the restoring forces as functions of the intermass displacements and velocities. The method was demonstrated by application to a four-degree-of-freedom system to identify the restoring forces. It is found that the identification results are relatively insensitive to measurement noise.

Udwadia, F. E.↗

A baroclinic quasigeostrophic open ocean model

A baroclinic quasigeostrophic open ocean model is presented, calibrated by a series of test problems, and demonstrated to be feasible and efficient for application to realistic mid-oceanic mesoscale eddy flow regimes. Two methods of treating the depth dependence of the flow, a finite difference method and a collocation method, are tested and intercompared. Sample Rossby wave calculations with and without advection are performed with constant stratification and two levels of nonlinearity, one weaker than and one typical of real ocean flows. Using exact analytical solutions for comparison, the accuracy and efficiency of the model is tabulated as a function of the computational parameters and stability limits set; typically, errors were controlled between 1 percent and 10 percent RMS after two wave periods. Further Rossby wave tests with realistic stratification and wave parameters chosen to mimic real ocean conditions were performed to determine computational parameters for use with real and simulated data. Finally, a prototype calculation with quasiturbulent simulated data was performed successfully, which demonstrates the practicality of the model for scientific use.

Miller, R. N.↗

A new algorithm for constrained nonlinear least-squares problems, part 1

A Gauss-Newton algorithm is presented for solving nonlinear least squares problems. The problem statement may include simple bounds or more general constraints on the unknowns. The algorithm uses a trust region that allows the objective function to increase with logic for retreating to best values. The computations for the linear problem are done using a least squares system solver that allows for simple bounds and linear constraints. The trust region limits are defined by a box around the current point. In its current form the algorithm is effective only for problems with small residuals, linear constraints and dense Jacobian matrices. Results on a set of test problems are encouraging.

Hanson, R. J.↗

Frequency domain laser velocimeter signal processor

The features and performance parameters of the frequency domain laser velocimeter signal processor (FDSP) for examining flow phenomena are described. A candidate for replacing the standard signal processor, the FDSP can function receiving bursts of as few as 150 photons in an instantaneous mode or down to a photon resolved regime for an averaging mode. the FDSP has a nonlinear A/D converter and data storage shift registers (high speed recorder group), an elliptic digital filter bank, a controlling microprocessor and frequency and signal gain control circuits (low speed data processing group). The high speed group amplifies, digitizes, identifies and stores burst waveforms and the low speed group identifies the signal frequency within the burst after the signal passes through digital filters. The FDSP can operate in a near-autonomous mode using photon correlation techniques and requires no human intervention for calibration.

Meyers, J. F.↗

Statistics and dynamics of persistent anomalies

Persistent anomalies with recurrent spatial patterns play an important role in the atmosphere's low-frequency variability. A connection between statistical and dynamical methods of description and prediction of persistent anomalies is established by computing and analyzing the empirical orthogonal functions (EOFs) in a simple deterministic model, on the one hand, and in Southern Hemisphere geopotential heights, on the other. The dynamical model is governed by the fully nonlinear, equivalent-barotropic vorticity equation on the sphere, with simplified forcing, dissipation and topography. The Southern Hemisphere data consist in gridded daily maps of 500 mb heights from June 1972 to July 1983. Two types of persistent anomalies appear in this time series, both having a strong wavenumber-three component; they differ by the value of the constant phase of this wave and by the strength of the wavenumber-one component. The first two EOFs bear a striking resemblance to these two patterns. It is concluded that the dynamical interpretation of EOFs is their pointing from the time mean to the most populated regions of the system's phase space. Pursuing this interpretation, a Markov-chain formulation of transitions from one persistent anomaly regime to another is introduced, and the implications for long-range forecasting are discussed.

Mo, Kingtse C.↗

Spillover, nonlinearity, and flexible structures

Many systems whose evolution in time is governed by Partial Differential Equations (PDEs) are linearized around a known equilibrium before Computer Aided Control Engineering (CACE) is considered. In this case, there are infinitely many independent vibrational modes, and it is intuitively evident on physical grounds that infinitely many actuators would be needed in order to control all modes. A more precise, general formulation of this grave difficulty (spillover problem) is due to A.V. Balakrishnan. A possible route to circumvention of this difficulty lies in leaving the PDE in its original nonlinear form, and adding the essentially finite dimensional control action prior to linearization. One possibly applicable technique is the Liapunov Schmidt rigorous reduction of singular infinite dimensional implicit function problems to finite dimensional implicit function problems. Omitting details of Banach space rigor, the formalities of this approach are given.

Bass, Robert W.↗

Identification of aerodynamic models for maneuvering aircraft

The method based on Fourier functional analysis and indicial formulation for aerodynamic modeling as proposed by Chin and Lan is extensively examined and improved for the purpose of general applications to realistic airplane configurations. Improvement is made to automate the calculation of model coefficients, and to evaluate more accurately the indicial integral. Test data of large angle-of-attack ranges for two different models, a 70 deg. delta wing and an F-18 model, are used to further verify the applicability of Fourier functional analysis and validate the indicial formulation. The results show that the general expression for harmonic motions throughout a range of k is capable of accurately modeling the nonlinear responses with large phase lag except in the region where an inconsistent hysteresis behavior from one frequency to the other occurs. The results by the indicial formulation indicate that more accurate results can be obtained when the motion starts from a low angle of attack where hysteresis effect is not important.

Lan, C. Edward↗

Two-dimensional viscous accretion disk models. I - On meridional circulations in radiative regions

It has been suggested that thin Keplerian disks, like rotating stars, cannot attain hydrostatic and thermal equilibrium simultaneously if the angular velocity is a function of radius only. To address the problem of vertical structure and stability of thin accretion disks, we performed 2D axisymmetric calculations by solving the set of fully nonlinear hydrodynamic equations, including radiation transport. The results indicate that the variation of the rotational velocity with height is too small to sustain hydrostatic equilibrium. Instead, the variation of the radial (inflow) velocity with height suffices to establish a stationary state. In particular, we find that for low values of the viscosity parameter alpha, there is mass outflow in the central parts of the disk close to the equatorial plane, and inflow only near the surface. Only for modestly high values of alpha, the flow throughout the disk is directed inward. Thus, the flow within the disk is the result of viscous radial inflow plus a circulatory flow directed outward in the midplane and inward near the surface.

Kley, W.↗

Experimental Investigation of Temperature Feedback Control Systems Applicable to Turbojet-engine Control

Temperature - fuel-flow and temperature-area feedback control systems were investigated as means of controlling tailpipe gas temperature of a turbojet engine during transient operation in the high-speed region. Proportional-plus-integral control was used in both systems, but in the temperature-area control system it was necessary to add nonlinear components to the basic proportional-plus-integral control to provide satisfactory transient response to a desired step increase in temperature. Time integral of temperature-error functions were used as criteria for determining optimum transient response. A description of engine dynamics was obtained from frequency-response data.

ENGINES, TURBOJET↗

Dynamic tests on the NASA Langley CSI evolutionary model

A modal analysis study, representing one of the anticipated 'Cooperative Spacecraft Structural Dynamics Experiments on the NASA Langley CSI Evolutionary Model', was carried out as a sub-task under the NASA/DLR collaboration in dynamics and control of large space systems. The CSI evolutionary testbed (CEM) is designed for the development of Controls-Structures Interaction (CSI) technology to improve space science platform pointing. For orbiting space structures like large flexible trusses, new identification challenges arise due to their specific dynamic characteristics (low frequencies and high modal density) on the one hand, and the limited possibilities of exciting such structures and measuring their responses on orbit on the other. The main objective was to investigate the modal identification potential of several different types of forcing functions that could possibly be realized with on-board excitation equipment using a minimum number of exciter locations as well as response locations. These locations were defined in an analytical test prediction process used to study the implications of measuring and analyzing the responses thus produced. It turned out that broadband excitation is needed for a general modal survey, but if only certain modes are of particular interest, combinations of exponentially decaying sine functions provide favorable excitation conditions as they allow to concentrate the available energy on the modes being of special interest. From a practical point-of-view structural nonlinearities as well as noisy measurements make the analysis more difficult, especially in the low frequency range and when the modes are closely spaced.

Troidl, H.↗

Unsteady three-dimensional thermal field prediction in turbine blades using nonlinear BEM

A time-and-space accurate and computationally efficient fully three dimensional unsteady temperature field analysis computer code has been developed for truly arbitrary configurations. It uses boundary element method (BEM) formulation based on an unsteady Green's function approach, multi-point Gaussian quadrature spatial integration on each panel, and a highly clustered time-step integration. The code accepts either temperatures or heat fluxes as boundary conditions that can vary in time on a point-by-point basis. Comparisons of the BEM numerical results and known analytical unsteady results for simple shapes demonstrate very high accuracy and reliability of the algorithm. An example of computed three dimensional temperature and heat flux fields in a realistically shaped internally cooled turbine blade is also discussed.

Martin, Thomas J.↗

Compression failure of angle-ply laminates

Test results from the compression loading of (+ or - Theta/ - or + Theta)(sub 6s) angle-ply IM7-8551-7a specimens, 0 less than or = Theta less than or = 90 degs, are presented. The observed failure strengths and modes are discussed, and typical stress-strain relations shown. Using classical lamination theory and the maximum stress criterion, an attempt is made to predict failure stress as a function of Theta. This attempt results in poor correlation with test results and thus a more advanced model is used. The model, which is based on a geometrically nonlinear theory, and which was taken from previous work, includes the influence of observed layer waviness. The waviness is described by the wave length and the wave amplitude. The theory is briefly described and results from the theory are correlated with test results. It is shown that by using levels of waviness observed in the specimens, the correlation between predictions and observations is good.

Peel, L. D.↗

Modeling of delamination damage evolution in laminated composites subjected to low velocity impact

This study examines the delamination evolution, under quasi-static conditions, of laminated polymeric composites with mechanically nonlinear resin rich interfaces. The constitutive behavior of the interface is represented by two models developed by Needleman and Tvegaard. These models assumed that the interfacial tractions, a function of only the interfacial displacement, will behave similarly to the interatomic forces generated during the interatomic seperation. The interface material's parameters control the load at which the delamination growth initiates and the final delamination size. A wide range of damage accumulation responses have been obtained by varying the model parameters. These results show that Tvergaard's model has been found to be better suited of the two models in predicting damage evolution for the configurations examined.

Lo, David C.↗

PLATSIM: A Simulation and Analysis Package for Large-Order Flexible Systems

The software package PLATSIM provides efficient time and frequency domain analysis of large-order generic space platforms. PLATSIM can perform open-loop analysis or closed-loop analysis with linear or nonlinear control system models. PLATSIM exploits the particular form of sparsity of the plant matrices for very efficient linear and nonlinear time domain analysis, as well as frequency domain analysis. A new, original algorithm for the efficient computation of open-loop and closed-loop frequency response functions for large-order systems has been developed and is implemented within the package. Furthermore, a novel and efficient jitter analysis routine which determines jitter and stability values from time simulations in a very efficient manner has been developed and is incorporated in the PLATSIM package. In the time domain analysis, PLATSIM simulates the response of the space platform to disturbances and calculates the jitter and stability values from the response time histories. In the frequency domain analysis, PLATSIM calculates frequency response function matrices and provides the corresponding Bode plots. The PLATSIM software package is written in MATLAB script language. A graphical user interface is developed in the package to provide convenient access to its various features.

Maghami, Peiman G.↗

Computation of Sensitivity Derivatives of Navier-Stokes Equations using Complex Variables

Accurate computation of sensitivity derivatives is becoming an important item in Computational Fluid Dynamics (CFD) because of recent emphasis on using nonlinear CFD methods in aerodynamic design, optimization, stability and control related problems. Several techniques are available to compute gradients or sensitivity derivatives of desired flow quantities or cost functions with respect to selected independent (design) variables. Perhaps the most common and oldest method is to use straightforward finite-differences for the evaluation of sensitivity derivatives. Although very simple, this method is prone to errors associated with choice of step sizes and can be cumbersome for geometric variables. The cost per design variable for computing sensitivity derivatives with central differencing is at least equal to the cost of three full analyses, but is usually much larger in practice due to difficulty in choosing step sizes. Another approach gaining popularity is the use of Automatic Differentiation software (such as ADIFOR) to process the source code, which in turn can be used to evaluate the sensitivity derivatives of preselected functions with respect to chosen design variables. In principle, this approach is also very straightforward and quite promising. The main drawback is the large memory requirement because memory use increases linearly with the number of design variables. ADIFOR software can also be cumber-some for large CFD codes and has not yet reached a full maturity level for production codes, especially in parallel computing environments.

Vatsa, Veer N.↗

Tensile Stress Acoustic Constants of Unidirectional Graphite/Epoxy Composites

Previously, the stress acoustic constants (SAC's) of unidirectional graphite/epoxy composites were measured to determine the nonlinear moduli of this material. These measurements were made under compressive loading in order to obtain the sufficient number of values needed to calculate these moduli. However, because their strength in tension along fiber directions can be several times greater, most composites are used under tensile loading. Thus, it is important to characterize the nonlinear properties of these materials in tension as well. The SAC's which are defined as the slope of the normalized change in ultrasonic "natural" velocity as a function of stress were measured in a unidirectional laminate of T300/5208 graphite/epoxy. Tensile load was applied along the fiber axis with the ultrasonic waves propagating perpendicular to the fiber direction. Changes in velocity were measured using a pulsed phase locked loop ultrasonic interferometer with the nominal frequency of the ultrasonic waves being 2.25 MHz.

Prosser, W. H.↗