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 145 records · Page 8

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.↗

Thermal mechanical analysis of sprag clutches

Work done at Case Western Reserve University on the Thermal Mechanical analysis of sprag helicopter clutches is reported. The report is presented in two parts. The first part is a description of a test rig for the measurement of the heat generated by high speed sprag clutch assemblies during cyclic torsional loading. The second part describes a finite element modeling procedure for sliding contact. The test rig provides a cyclic torsional load of 756 inch-pounds at 5000 rpm using a four-square arrangement. The sprag clutch test unit was placed between the high speed pinions of the circulating power loop. The test unit was designed to have replaceable inner ad outer races, which contain the instrumentation to monitor the sprag clutch. The torque loading device was chosen to be a water cooled magnetic clutch, which is controlled either manually or through a computer. In the second part, a Generalized Eulerian-Lagrangian formulation for non-linear dynamic problems is developed for solid materials. This formulation is derived from the basic laws and axioms of continuum mechanics. The novel aspect of this method is that we are able to investigate the physics in the spatial region of interest as material flows through it without having to follow material points. A finite element approximation to the governing equations is developed. Iterative Methods for the solution of the discrete finite element equations are explored. A FORTRAN program to implement this formulation is developed and a number of solutions to problems of sliding contact are presented.

Mullen, Robert L.↗

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.↗

Solder creep-fatigue interactions with flexible leaded parts

With flexible leaded parts, the solder-joint failure process involves a complex interplay of creep and fatigue mechanisms. To better understand the role of creep in typical multi-hour cyclic loading conditions, a specialized non-linear finite-element creep simulation computer program has been formulated. The numerical algorithm includes the complete part-lead-solder-PWB system, accounting for strain-rate dependence of creep on applied stress and temperature, and the role of the part-lead dimensions and flexibility that determine the total creep deflection (solder strain range) during stress relaxation. The computer program has been used to explore the effects of various solder creep-fatigue parameters such as lead height and stiffness, thermal-cycle test profile, and part/board differential thermal expansion properties. One of the most interesting findings is the strong presence of unidirectional creep-ratcheting that occurs during thermal cycling due to temperature dominated strain-rate effects. To corroborate the solder fatigue model predictions, a number of carefully controlled thermal-cycle tests have been conducted using special bimetallic test boards.

Ross, R. G., Jr.↗

Combustion noise and combustion instabilities in propulsion systems

This paper is concerned with some aspects of non-linear behavior of unsteady motions in combustion chambers. The emphasis is on conditions under which organized oscillations having discrete frequencies may exist in the presence of random motions. In order to treat the two types of motions together, and particularly to investigate coupling between noise and combustion instabilities, the unsteady field is represented as a synthesis of acoustic modes having time-varying amplitudes. Each of the amplitudes are written as the sum of two parts, one associated with the random field and the remainder representing the organized oscillations. After spatial averaging, the general problem is reduced to solution of a set of second-order ordinary differential equations whose structure depends on the sorts of nonlinear processes accounted for. This formulation accommodates any physical process; in particular, terms are included to represent noise sources, although only limited modeling is discussed. Our results suggest that random sources of noise have only small effects on combustion instabilities and seem not to be a cause of unstable motions. However, the coupling between the two sorts of unsteady motions may be important as an essential process in a proposed scheme for noise control. It is now a familiar observation that many nonlinear deterministic systems are capable of exhibiting apparently random motions called 'chaos.' This is a particularly interesting possibility for systems which also executed non-deterministic random motions. In combustion chambers, a nonlinear deterministic system (acoustical motions) exists in the presence of noise produced by flow separation, turbulent motions, and energy released by combustion processes. The last part of the paper is directed to the matter of discovering whether or not chaotic motions exist in combustion systems. Analysis has not progressed sufficiently far to answer the question. We report here recent results of processing data taken in one combustor to determine the dimensions of any attractors in the motions. No evidence has been found for chaos in the strict sense, but the method seems to be an important means of investigating the nonlinear behavior of combustion systems.

F. E. C. Culick↗

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↗

Bypass transition in compressible boundary layers

Transition to turbulence in aerospace applications usually occurs in a strongly disturbed environment. For instance, the effects of free-stream turbulence, roughness and obstacles in the boundary layer strongly influence transition. Proper understanding of the mechanisms leading to transition is crucial in the design of aircraft wings and gas turbine blades, because lift, drag and heat transfer strongly depend on the state of the boundary layer, laminar or turbulent. Unfortunately, most of the transition research, both theoretical and experimental, has focused on natural transition. Many practical flows, however, defy any theoretical analysis and are extremely difficult to measure. Morkovin introduced in his review paper the concept of bypass transition as those forms of transition which bypass the known mechanisms of linear and non-linear transition theories and are currently not understood by experiments. In an effort to better understand the mechanisms leading to transition in a disturbed environment, experiments are conducted studying simpler cases, viz. the effects of free stream turbulence on transition on a flat plate. It turns out that these experiments are very difficult to conduct, because generation of free stream turbulence with sufficiently high fluctuation levels and reasonable homogeneity is non trivial. For a discussion see Morkovin. Serious problems also appear due to the fact that at high Reynolds numbers the boundary layers are very thin, especially in the nose region of the plate where the transition occurs, which makes the use of very small probes necessary. The effects of free-stream turbulence on transition are the subject of this research and are especially important in a gas turbine environment, where turbulence intensities are measured between 5 and 20 percent, Wang et al. Due to the fact that the Reynolds number for turbine blades is considerably lower than for aircraft wings, generally a larger portion of the blade will be in a laminar transitional state. The effects of large free stream turbulence in compressible boundary layers at Mach numbers are examined both in the subsonic and transonic regime using direct numerical simulations. The flow is computed over a flat plate and curved surface. while many applications operate in the transonic regime. Due the nature of their numerical scheme, a non-conservation formulation of the Navier-Stokes equations, it is a non-trivial extension to compute flow fields in the transonic regime. This project aims at better understanding the effects of large free-stream turbulence in compressible boundary layers at mach number both in the subsonic and transonic regime using direct numerical simulations. The present project aims at computing the flow over a flat plate and curved surface. This research will provide data which can be used to clarify mechanisms leading to transition in an environment with high free stream turbulence. This information is useful for the development of turbulence models, which are of great importance for CFD applications, and are currently unreliable for more complex flows, such as transitional flows.

Vandervegt, J. J.↗

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↗