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 325 records · Page 18

Identification of Linear and Nonlinear Aerodynamic Impulse Responses Using Digital Filter Techniques

This paper discusses the mathematical existence and the numerically-correct identification of linear and nonlinear aerodynamic impulse response functions. Differences between continuous-time and discrete-time system theories, which permit the identification and efficient use of these functions, will be detailed. Important input/output definitions and the concept of linear and nonlinear systems with memory will also be discussed. It will be shown that indicial (step or steady) responses (such as Wagner's function), forced harmonic responses (such as Theodorsen's function or those from doublet lattice theory), and responses to random inputs (such as gusts) can all be obtained from an aerodynamic impulse response function. This paper establishes the aerodynamic impulse response function as the most fundamental, and, therefore, the most computationally efficient, aerodynamic function that can be extracted from any given discrete-time, aerodynamic system. The results presented in this paper help to unify the understanding of classical two-dimensional continuous-time theories with modern three-dimensional, discrete-time theories. First, the method is applied to the nonlinear viscous Burger's equation as an example. Next the method is applied to a three-dimensional aeroelastic model using the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) code and then to a two-dimensional model using the CFL3D Navier-Stokes code. Comparisons of accuracy and computational cost savings are presented. Because of its mathematical generality, an important attribute of this methodology is that it is applicable to a wide range of nonlinear, discrete-time problems.

Silva, Walter A.↗

Research study on stabilization and control: Modern sampled-data control theory. Design of the large space telescope system

Conditions of self-sustained oscillations in a two-axis model of the nonlinear LST system are studied. The describing function of the CMG frictional nonlinearity of the LST system is used for the analysis, as well as continuous-data and discrete-data models of the simplified LST control system. A numerical-iterative method is described for the analysis of the two-axis system. Approximation methods and the direct plotting of the stability equation are implemented in the study. It is shown that although the dynamics of the two axes are identical, the amplitudes of self-sustained oscillations in the two axes may in principle be different. Analysis shows that the LST systems are of equal amplitudes but with 180-degree phase shift.

Kuo, B. C.↗

SHARP simulation of discontinuities in highly convective steady flow

For steady multidimesional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

Simple high-accuracy resolution program for convective modelling of discontinuities

For steady multidimensional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

Survey and development of finite elements for nonlinear structural analysis. Volume 2: Nonlinear shell finite elements

The development of two new shell finite elements for applications to large deflection problems is considered. The elements in question are doubly curved and of triangular and quadrilateral planform. They are restricted to small strains of elastic materials, and can accommodate large rotations. The elements described, which are based on relatively simple linear elements, make use of a new displacement function approach specifically designed for strongly nonlinear problems. The displacement function development for nonlinear applications is based on certain beam element formulations, and the strain-displacement equations are of a shallow shell type. Additional terms were included in these equations in an attempt to avoid the large errors characteristic of shallow shell elements in certain types of problems. An incremental nonlinear solution procedure specifically adopted to the element formulation was developed. The solution procedure is of combined incremental and total Lagrangian type, and uses a new updating scheme. A computer program was written to evaluate the developed formulations. This program can accommodate small element groups in arbitrary arrangements. Two simple programs were successfully solved. The results indicate that this new type of element has definite promise and should be a fruitful area for further research.

Source record↗

Stability analysis of the discrete-data Large-Space-Telescope system

This paper considers the application of the discrete describing function to the stability analysis of the Large-Space Telescope (LST) system. An analytical model of the CMG gimbal friction is derived, which is then used to arrive at a closed-form analytical expression for the discrete describing function of the nonlinearity. The analysis is used for the study of the fine-pointing stability of the LST vehicle. Simulation results corroborate the conclusions from the analytical analysis.

Kuo, B. C.↗

Satellite observations of suprathermal electron bursts

Data from the ISIS 2 satellite are used to study the spatial distribution of the suprathermal electron bursts discovered by high-altitude sounding rockets. It is found that the bursts form an integral part of the boundary plasma sheet and are linked with the formation of discrete auroral arcs. The maximum energy of the bursts is always below the peak energy in neighboring discrete forms, and it appears that if the bursts fail to reach the minimum energy found to be necessary to sustain a discrete arc they collapse within a few seconds. It is suggested that in a pass across the auroral oval the overall envelope of the maximum energy in bursts and the peak energy in discrete auroral forms reflect the variation of the available acceleration (perhaps a parallel electric potential) across the auroral oval while the energy structure in individual forms is a function of the nonlinear processes which sustain them. This explanation appears to resolve the apparent discrepancy in scale thickness between rocket and satellite inverted V's.

Johnstone, A. D.↗

New calculation of the properties of the positronium ion

The positronium negative ion Ps(-), the system composed of two electrons and a positron, has been reinvestigated theoretically. Using a Hylleraas wave function with two nonlinear parameters and more than 200 linear terms, excellent values of binding energy and annihilation lifetime of the particle-stable 1Se ground state have been obtained. In addition, the question of stability of the 3Pe state discussed by Mills has been examined and it is agreed with his conclusion that the state is probably not stable against breakup into Ps(n = 2) + e(-). Improved limits on the critical 'positron' mass for binding the 3Pe state have also been obtained.

Bhatia, A. K.↗

On Fofonoff's Mode

The classical barotropic ocean model of Veronis (1966) is reexamined to verify the possible existence of Fofonoff's (1954) steady free inertial mode by considering the balance where J is the Jacobian, eta is the streamfunction of the horizontal velocity, and R is a parameter of nonlinearity. A linear functional relationship between the potential vorticity and the streamfunction and inertial boundary layers along the ocean basin was suggested by Fofonoff. However, the numerical integration showed that the system is barotropically unstable. The typical instantaneous field is given. The statistical equilibrium state eta is also given. Although it resembles Fofonoff's solution, it is maintained by the Reynolds stress field induced by the eddies in addition to the wind stress and dissipation. The plot of eta versus R delta sq +y is given. It is shown that the Fofonoff's linear relation between the potential vorticity and the streamfunction does not hold uniformly over the entire basin. The above results show that the Fofonoff's solution can not be realized physically.

Merkine, L. O.↗

An accurate and efficient satellite long-term orbit predictor employing 'fictitious' mean orbital elements

By using Von Zeipel's generating function procedure the perturbing earth gravitational potential is averaged with respect to the fast variable (mean anomaly) and a set of 'fictitous' mean orbital elements which can be used as a long-term satellite orbit predictor is obtained. The set of elements is shown to be a function of the nonlinear square of the second zonal harmonic coefficient. It is found that the long-term orbit prediction using the 'fictitous' mean elements is as accurate as that using the osculating elements, but has a computing speed about two orders of magnitude faster. For short-term orbit predictions, the osculating elements approach must be used.

Tang, Charles C. H.↗

The role of shock induced trailing-edge separation in limit cycle oscillations

The potential role of shock induced trailing edge separation (SITES) in limit cycle oscillations (LCO) was established. It was shown that the flip-flop characteristics of transition to and from SITES as well as its hysteresis could couple with wing modes with torsional motion and low damping. This connection led to the formulation of a very simple nonlinear math model using the linear equations of motion with a nonlinear step forcing function with hysteresis. A finite difference solution with time was developed and calculations were made for the F-111 TACT were used to determine the step forcing function due to SITES transition. Since no data were available for the hysteresis, a parameter study was conducted allowing the hysteresis effect to vary. Very small hysteresis effects, which were within expected bounds, were required to obtain reasonable response levels that essentially agreed with flight test results. Also in agreement with wind tunnel tests, LCO calculations for the 1/6 scale F-111 model showed that the model should have not experienced LCO.

Cunningham, Atlee M., Jr.↗

Transfer Functions Via Laplace- And Fourier-Borel Transforms

Approach to solution of nonlinear ordinary differential equations involves transfer functions based on recently-introduced Laplace-Borel and Fourier-Borel transforms. Main theorem gives transform of response of nonlinear system as Cauchy product of transfer function and transform of input function of system, together with memory effects. Used to determine responses of electrical circuits containing variable inductances or resistances. Also possibility of doing all noncommutative algebra on computers in such symbolic programming languages as Macsyma, Reduce, PL1, or Lisp. Process of solution organized and possibly simplified by algebraic manipulations reducing integrals in solutions to known or tabulated forms.

Can, Sumer↗

Roll plus maneuver load alleviation control system designs for the active flexible wing wind-tunnel model

Three designs for controlling loads while rolling for the Active Flexible Wing (AFW) are discussed. The goal is to provide good roll control while simultaneously limiting the torsion and bending loads experienced by the wing. The first design uses Linear Quadratic Gaussian/Loop Transfer Recovery (LQG/LTR) modern control methods to control roll rate and torsional loads at four different wing locations. The second design uses a nonlinear surface command function to produce surface position commands as a function of current roll rate and commanded roll rate. The final design is a flutter suppression control system. This system stabilizes both symmetric and axisymmetric flutter modes of the AFW.

Moore, Douglas B.↗

Linear dynamic analysis of the Space Station Freedom thermal/environmental control system

A methodology for dynamic analysis of the Space Station Freedom thermal/environmental control system is described which uses linear system analysis techniques to determine the ability of the system to accommodate thermal load transients. A nonlinear integrated model of the system as it exists at Stage 6, Man Tended Capability (MTC), is described which encompasses the cabin air cooling system, avionics air cooling system, internal and external thermal control systems. The nonlinear model is linearized at an assumed design operating point, and the linear model is evaluated by comparison with the nonlinear model. Transfer functions are derived to determine the response of the heat rejection system to time-varying thermal loads applied at various points in the heat acquisition system. The frequency response calculated from the transfer functions is used to determine allowable operating envelopes of thermal load magnitude and frequency for each point of application of thermal load. The dynamic operating envelopes provide a quantitative indication of the degree and duration of peak loads which the thermal/environmental system is able to transport and reject.

Tandler, John↗

Computation of a controlled store separation from a cavity

Coupling of the Reynolds-averaged Navier-Stokes equations, rigid-body dynamics, and a pitch attitude control law is demonstrated in two- and three-dimensions. The application problem was the separation of a canard-controlled store from an open-flow rectangular cavity bay at a freestream Mach number of 1.2. The transient flowfield was computed using a diagonal scheme in an overset mesh framework, with the resultant aerodynamic loads used as the forcing functions in the nonlinear dynamics equations. The proportional and rate gyro sensitivities were computed a priori using pole placement techniques for the linearized dynamical equations. These fixed gain values were used in the controller for the nonlinear simulation. Reasonable comparison between the full and linearized equations for a perturbed two-dimensional missile was found. Also in two-dimensions, a controlled store was found to possess improved separation characteristics over a canard-fixed store. In three-dimensions, trajectory comparisons with wind-tunnel data for the canard-fixed case will be made. In addition, it will be determined if a canard-controlled store is an effective means of improving cavity store separation characteristics.

Atwood, Christopher A.↗

Contrast Gain Control Model Fits Masking Data

We studied the fit of a contrast gain control model to data of Foley (JOSA 1994), consisting of thresholds for a Gabor patch masked by gratings of various orientations, or by compounds of two orientations. Our general model includes models of Foley and Teo & Heeger (IEEE 1994). Our specific model used a bank of Gabor filters with octave bandwidths at 8 orientations. Excitatory and inhibitory nonlinearities were power functions with exponents of 2.4 and 2. Inhibitory pooling was broad in orientation, but narrow in spatial frequency and space. Minkowski pooling used an exponent of 4. All of the data for observer KMF were well fit by the model. We have developed a contrast gain control model that fits masking data. Unlike Foley's, our model accepts images as inputs. Unlike Teo & Heeger's, our model did not require multiple channels for different dynamic ranges.

Watson, Andrew B.↗