Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “non-linear models”

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 109 records · Page 6

Filtering of non-linear instabilities

For Courant numbers larger than one and cell Reynolds numbers larger than two, oscillations and in some cases instabilities are typically found with implicit numerical solutions of the fluid dynamics equations. This behavior has sometimes been associated with the loss of diagonal dominance of the coefficient matrix. It is shown that these problems can be related to the choice of the spatial differences, with the resulting instability related to aliasing or nonlinear interaction. Appropriate filtering can reduce the intensity of these oscillations and possibly eliminate the instability. These filtering procedures are equivalent to a weighted average of conservation and nonconservation differencing. The entire spectrum of filtered equations retains a three point character as well as second order spatial accuracy. Burgers equation was considered as a model.

Khosla, P. K.↗

Multiple Equilibria of the Barotropic Vorticity Equation on a Sphere

Multiple states of the barotropic vorticity equation in which the balance is between the first and third terms on the r.h.s. of (1) are given. Solutions of this type were also considered by Charney and Devore (1979) in their discussion of thermally, rather than topographically, forced waves. Whereas in the case of topographic forcing the multiplicity arises from what they called form-drag instability in the wave-zonal flow interaction, in the thermal forcing case the associated instability appears to be the Rossby-wave instability discussed by Lorenz (1972), and the multiple states of the highly truncated model proved to be unstable when more degrees of freedom were added. The multiple statistically steady solutions described are thus novel in that they do not involve form-drag instability, they have stable statistics in calculations with a large number of degrees of freedoms, and they occur on the sphere, with no artificial confinement in a resonant cavity. Furthermore, the solutions are obtained with no external forcing of the zonal flow. The full non-linear equations for two-dimensional non-divergent motion between smooth, rigid boundaries on a sphere were used.

Suarez, M. J.↗

A finite difference solution for the propagation of sound in near sonic flows

An explicit time/space finite difference procedure is used to model the propagation of sound in a quasi one-dimensional duct containing high Mach number subsonic flow. Nonlinear acoustic equations are derived by perturbing the time-dependent Euler equations about a steady, compressible mean flow. The governing difference relations are based on a fourth-order, two-step (predictor-corrector) MacCormack scheme. The solution algorithm functions by switching on a time harmonic source and allowing the difference equations to iterate to a steady state. The principal effect of the non-linearities was to shift acoustical energy to higher harmonics. With increased source strengths, wave steepening was observed. This phenomenon suggests that the acoustical response may approach a shock behavior at at higher sound pressure level as the throat Mach number aproaches unity. On a peak level basis, good agreement between the nonlinear finite difference and linear finite element solutions was observed, even through a peak sound pressure level of about 150 dB occurred in the throat region. Nonlinear steady state waveform solutions are shown to be in excellent agreement with a nonlinear asymptotic theory.

Hariharan, S. I.↗

A finite difference solution for the propagation of sound in near sonic flows

An explicit time/space finite difference procedure is used to model the propagation of sound in a quasi one-dimensional duct containing high Mach number subsonic flow. Nonlinear acoustic equations are derived by perturbing the time-dependent Euler equations about a steady, compressible mean flow. The governing difference relations are based on a fourth-order, two-step (predictor-corrector) MacCormack scheme. The solution algorithm functions by switching on a time harmonic source and allowing the difference equations to iterate to a steady state. The principal effect of the non-linearities was to shift aocustical energy to higher harmonics. With increased source strength, wave steepening was observed. This phenomenon suggests that the acoustical response may approach a shock behavior at higher sound pressure level as the throat Mach number approaches unity. On a peak level basis, good agreement between the nonlinear finite difference and linear finite element solutions was observed, even through a peak sound pressure level of about 150 dB occurred in the throat region. Nonlinear steady state waveform solutions are shown to be in excellent agreement with a nonlinear asymptotic theory. Previously announced in STAR as N83-30167

Hariharan, S. I.↗

A numerical solution of the Navier-Stokes equations for supercritical fluid thermodynamic analysis

An explicit numerical solution of the compressible Navier-Stokes equations is applied to the thermodynamic analysis of supercritical oxygen in the Apollo cryogenic storage system. The wave character is retained in the conservation equations which are written in the basic fluid variables for a two-dimensional Cartesian coordinate system. Control-volume cells are employed to simplify imposition of boundary conditions and to ensure strict observance of local and global conservation principles. Non-linear real-gas thermodynamic properties responsible for the pressure collapse phenomonon in supercritical fluids are represented by tabular and empirical functions relating pressure and temperature to density and internal energy. Wall boundary conditions are adjusted at one cell face to emit a prescribed mass flowrate. Scaling principles are invoked to achieve acceptable computer execution times for very low Mach number convection problems. Detailed simulations of thermal stratification and fluid mixing occurring under low acceleration in the Apollo 12 supercritical oxygen tank are presented which model the pressure decay associated with de-stratification induced by an ordinary vehicle maneuver and heater cycle operation.

Heinmiller, P. J.↗

A new approach to resonance line scattering in planetary atmospheres

The invariant imbedding technique is applied to the problem of resonance line scattering with frequency redistribution. By applying the principles of invariance to the resonance line scattering problem a set of non-linear, integrodifferential equations for the scattering and transmission functions are derived. Thus, an alternate approach to the resonance line scattering problem is presented. Rather than solving for the specific intensity directly the invariant imbedding equations may be solved for the scattering and transmission functions from which the specific intensity may be calculated. The technique is very general and in principle, can accommodate scattering in inhomogeneous atmospheres with arbitrary frequency redistribution and scattering phase functions. Solutions to the invariant imbedding equations in some simple cases are presented. If the scattering is isotropic and completely noncoherent and the atmosphere may be modeled as homogeneous, the scattering and transmission functions may be expressed in terms of generalized H, X and Y functions. This technique is extended to the problems of coupled spectral lines and scattering of radiation from internal sources.

Yelle, Roger V.↗

Lateral-directional stability investigation of the space shuttle orbiter at Mach 6

Lateral-directional aerodynamic data and oil-flow visualization results are pesented from four tests (conducted in two facilities using three models) which were designed to verify the hypersonic aerodynamics of the Space Shuttle orbiter at M=6. Comparisons of measured lateral-directional stability data and oil-flow results between the tests show excellent agreement, especially considering the nonuniform, unpredictable flow which occurs in the vicinity of the vertical tail. Results ere shown to be sensitive to Reynolds number with the higher Reynolds number cases producing more stable values and also showing good agreement with flight values. The results also show that the effects of Reynolds number, angle of attack, and angle of sideslip on the lateral-directional stability of future entry configuration should be carefully assessed for non-linearities.

Calloway, R. L.↗

A digital strategy for manometer dynamic enhancement

Application of digital signal processing techniques to improve the non-linear dynamic characteristics of a sonar-type mercury manometer is described. The dynamic enhancement strategy quasi-linearizes the manometer characteristics and improves the effective bandwidth in the context of a wind-tunnel pressure regulation system. Model identification data and real-time hybrid simulation data demonstrate feasibility of approach.

Stoughton, J. W.↗

Filtering of non-linear instabilities

For Courant numbers larger than one and cell Reynolds numbers larger than two, oscillations and in some cases instabilities are typically found with implicit numerical solutions of the fluid dynamics equations. This behavior has sometimes been associated with the loss of diagonal dominance of the coefficient matrix. It is shown here that these problems can in fact be related to the choice of the spatial differences, with the resulting instability related to aliasing or nonlinear interaction. Appropriate 'filtering' can reduce the intensity of these oscillations and in some cases possibly eliminate the instability. These filtering procedures are equivalent to a weighted average of conservation and non-conservation differencing. The entire spectrum of filtered equations retains a three-point character as well as second-order spatial accuracy. Burgers equation has been considered as a model. Several filters are examined in detail, and smooth solutions have been obtained for extremely large cell Reynolds numbers.

Khosla, P. K.↗

Advanced turbo-prop airplane interior noise reduction-source definition

Acoustic pressure amplitudes and phases were measured in model scale on the surface of a rigid semicylinder mounted in an acoustically treated wind tunnel near a prop-fan (an advanced turboprop with many swept blades) model. Operating conditions during the test simulated those of a prop-fan at 0.8 Mach number cruise. Acoustic pressure amplitude and phase contours were defined on the semicylinder surface. Measurements obtained without the semi-cylinder in place were used to establish the magnitude of pressure doubling for an aircraft fuselage located near a prop-fan. Pressure doubling effects were found to be 6dB at 90 deg incidence decreasing to no effect at grazing incidence. Comparisons of measurements with predictions made using a recently developed prop-fan noise prediction theory which includes linear and non-linear source terms showed good agreement in phase and in peak noise amplitude. Predictions of noise amplitude and phase contours, including pressure doubling effects derived from test, are included for a full scale prop-fan installation.

Magliozzi, B.↗

Identification and dual adaptive control of a turbojet engine

The objective of this paper is to utilize the design methods of modern control theory to realize a 'dual-adaptive' feedback control unit for a highly non-linear single spool airbreathing turbojet engine. Using a very detailed and accurate simulation of the non-linear engine as the data source, linear operating point models of unspecified dimension are identified. Feedback control laws are designed at each operating point for a prespecified set of sampling rates using sampled-data output regulator theory. The control system sampling rate is determined by an adaptive sampling algorithm in correspondence with turbojet engine performance. The result is a 'dual-adpative' control law that is functionally dependent upon the sampling rate selected and environmental operating conditions. Simulation transients demonstrate the utility of the dual-adaptive design to improve on-board computer utilization while maintaining acceptable levels of engine performance.

Merrill, W.↗

Improving the performance of soft decision Viterbi decoding in a non-Gaussian environment through non-linear quantization

The performance of Viterbi decoding in a non-Gaussian environment is investigated using a nonlinear quantization strategy. The channel model consists of a convolutionally encoded BPSK signal transmitted to a satellite where it is corrupted with additive white Gaussian noise and pulsed radio frequency interference (RFI). The resultant signal is then passed through a satellite nonlinearity and transmitted to a ground station where it is coherently detected. Interleaving is assumed in order to make the channel memoryless. The presence of RFI makes the channel statistics non-Gaussian, leading to a nonlinear log-likelihood function. A near optimum quantization scheme is found by maximizing a channel parameter, or by matching the quantizer to the log-likelihood function in a mean square error sense. Bit error rate performance improvement is achieved by using such nonlinear quantization.

Mcgregor, D. N.↗

The Construction of Curves and Surfaces Using Numerical Optimization Techniques

Numerical optimization techniques are playing an increasing role in curve and surface construction. Often difficult problems in curve and surface construction, especially when some aspect of shape control is involved, can be phrased as a constrained optimization problem. Four such classes of problems are explored: parametric curve fitting with non-linear shape constraints; explicit surface fitting with linear shape constraints; surface fitting to scattered data giving rise to ill-posed problems; finally, variable knot problems. In each of these problems there is a nonlinear aspect: either the shape of the curve or surface is important for manufacturing or engineering reasons or the shape affects the convergence of numerical algorithms which use the curve or surface or the placement of knots affects the accuracy of the fits. In all cases the class of functions used is that of parametric spline curves and tensor or direct product spline surfaces. The reason for choosing this class is that splines provide flexible models that are easily evaluated and stored. Furthermore, the B-spline representation of splines leads to convenient expressions for shape control over regions.

Ferguson, D. R.↗

Position Control for Non-linear, Multiple-Link Robots

Several approaches to complicated control problem examined. Report surveys methods for controlling motion of robot manipulators. Applies to coupled, highly nonlinear multiple-link robots. Presents adaptive control technique based on discrete linear model of robot arm, obtained by linearizing nonlinear dynamical equations about nominal trajectory. Parameters of model calculated with adaptive parameter-identification algorithm. Based on system model, optimal control input computed using position and velocity feedback.

Seraji, H.↗

Near UV atmospheric absorption measurements from the DC-8 aircraft during the 1987 airborne Antarctic ozone experiment

During the Airborne Antarctic Ozone Experiment from 28 August to 30 September 1987 near UV zenith scattered sky measurements were made over Antarctic from the NASA DC-8 aircraft using a one third m spectrograph equipped with a diode-array detector. Scattered sky light data in the wavelength range 348 nm to 388 nm was spectrally analyzed for O3, NO2, OClO, and BrO column abundances. Slant column abudances of O3, NO2, OClO and BrO were determined, using a computer algorithm of non-linear and linear least square correlation of Antarctic scattered sky spectra to laboratory absorption cross section data. Using measured vertical electrochemical sonde ozone profiles from Palmer, Halley Bay, and the South Pole Stations the slant columns of O3 were converted into vertical column abundances. The vertical column amounts of NO2, OClO, and BrO were derived using vertical profiles calculated by a chemical model appropriate for Antarctica. NO2 vertical column abundances show steep latitudinal decrease with increasing latitude for all 13 flights carried out during the mission. In the regions where NO2 abudances are low, OClO and BrO were observed. The spatial and temporal vertical column abundances of these species are discussed in the context of the chemistry and dynamics in the antarctic polar vortex during the austral spring.

Wahner, A.↗

Effect of wave localization on plasma instabilities

The Anderson model of wave localization in random media is involved to study the effect of solar wind density turbulence on plasma processes associated with the solar type III radio burst. ISEE-3 satellite data indicate that a possible model for the type III process is the parametric decay of Langmuir waves excited by solar flare electron streams into daughter electromagnetic and ion acoustic waves. The threshold for this instability, however, is much higher than observed Langmuir wave levels because of rapid wave convection of the transverse electromagnetic daughter wave in the case where the solar wind is assumed homogeneous. Langmuir and transverse waves near critical density satisfy the Ioffe-Reigel criteria for wave localization in the solar wind with observed density fluctuations -1 percent. Numerical simulations of wave propagation in random media confirm the localization length predictions of Escande and Souillard for stationary density fluctations. For mobile density fluctuations localized wave packets spread at the propagation velocity of the density fluctuations rather than the group velocity of the waves. Computer simulations using a linearized hybrid code show that an electron beam will excite localized Langmuir waves in a plasma with density turbulence. An action principle approach is used to develop a theory of non-linear wave processes when waves are localized. A theory of resonant particles diffusion by localized waves is developed to explain the saturation of the beam-plasma instability. It is argued that localization of electromagnetic waves will allow the instability threshold to be exceeded for the parametric decay discussed above.

Levedahl, William Kirk↗

Report on the multicolor pyrometry session at the Second NASA Non-Contact Temperature-Measurement Workshop

The session met jointly with the other special sessions. Due to time constraints, the concerns of multicolor pyrometry were not addressed in depth. Multicolor pyrometry attempts to determine an object's temperature and emittance as a function of wavelength from simultaneous measurement of the spectral radiance of the object at multiple wavelengths. The main issues are those of the most appropriate model for the wavelength dependence of the emittance, the number of undetermined parameters, and the error in the measured temperature. It has long been known that it is inappropriate to attempt to model the emittance with a functional form that contains as many unknowns as wavelengths. Such attempts will lead to gross errors in the measured temperature. However, more recent work by Hunter and Khan has indicated that emittance can be modeled by a function with many fewer degrees of freedom than sampled wavelengths (e.g., it may be appropriate to use a second order polynomial to fit data from a forty-color pyrometer). In a platform lecture, Khan described algorithms for determination of the parameters by minimization of the errors between the model and the data in a least-squares sense. The first was a linear method and the second was a non-linear method where the emittance could assume any functional form. Both methods were shown to be very sensitive to random noise and non-linearity in the measurement system.

Spjut, R. Erik↗

Experimental Investigation of a Hypersonic Glider Configuration at a Mach Number of 6 and at Full-Scale Reynolds Numbers

The aerodynamic characteristics of a hypersonic glider configuration, consisting of a slender ogive cylinder with three highly swept wings, spaced 120 apart, with the wing chord equal to the body length, were investigated experimentally at a Mach number of 6 and at Reynolds numbers from 6 to 16 million. The objectives were to evaluate the theoretical procedures which had been used to estimate the performance of the glider, and also to evaluate the characteristics of the glider itself. A principal question concerned the viscous drag at full-scale Reynolds number, there being a large difference between the total drags for laminar and turbulent boundary layers. It was found that the procedures which had been applied for estimating minimum drag, drag due to lift, lift curve slope, and center of pressure were generally accurate within 10 percent. An important exception was the non-linear contribution to the lift coefficient which had been represented by a Newtonian term. Experimentally, the lift curve was nearly linear within the angle-of-attack range up to 10 deg. This error affected the estimated lift-drag ratio. The minimum drag measurements indicated that substantial amounts of turbulent boundary layer were present on all models tested, over a range of surface roughness from 5 microinches maximum to 200 microinches maximum. In fact, the minimum drag coefficients were nearly independent of the surface smoothness and fell between the estimated values for turbulent and laminar boundary layers, but closer to the turbulent value. At the highest test Reynolds numbers and at large angles of attack, there was some indication that the skin friction of the rough models was being increased by the surface roughness. At full-scale Reynolds number, the maximum lift-drag ratio with a leading edge of practical diameter (from the standpoint of leading-edge heating) was 4.0. The configuration was statically and dynamically stable in pitch and yaw, and the center of pressure was less than 2-percent length ahead of the centroid of plan-form area.

Seiff, Alvin↗