SEARCH · Engineering Papers

Results for “nonlinear differential equations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Prediction of a propellant tank pressure history using state space methods

An analysis of the time response of a propellant supply system operating in the blowdown mode is presented. The supply system is part of a pump-fed propulsion system intended for use on interplanetary spacecraft. As such, the supply system must provide the pump with propellant at sufficient pressure to avoid pump cavitation. The system, consisting of the tank, the liquid propellant, the pressurant gas and propellant vapor mixture, and a film layer separating the liquid and vapor phases, is analyzed using the principles of mass and energy conservation. The resulting set of ordinary, coupled, nonlinear differential equations for the thermodynamic state variables is integrated as an initial value problem. The resulting histories of total pressure, propellant vapor pressure, propellant liquid temperature, film layer temperature, propellant vapor/pressurant gas temperature, propellant vapor mass, and propellant liquid mass enable the calculation of the net positive suction head available at the pump which determines the viability of the pump-fed system concept when operated in the blowdown mode.

Estey, P. N.↗

NASA NTRS1986

Predicting the Dynamic Behavior of High-Speed Roller Bearings

Computer program developed for solving nonlinear differential equations of motion of high-speed, lightly loaded, cylindrical roller bearing. Program organized into three major parts: Main program, controls input and output; differential equation solver; and subroutines, used to evaluate vector function. Program has capability of performing either twodimensional or three-dimensional simulation.

Conry, T. F.↗

NASA NTRS1984

Chemoviscosity modeling for thermosetting resins - I

A new analytical model for chemoviscosity variation during cure of thermosetting resins was developed. This model is derived by modifying the widely used WLF (Williams-Landel-Ferry) Theory in polymer rheology. Major assumptions involved are that the rate of reaction is diffusion controlled and is linearly inversely proportional to the viscosity of the medium over the entire cure cycle. The resultant first order nonlinear differential equation is solved numerically, and the model predictions compare favorably with experimental data of EPON 828/Agent U obtained on a Rheometrics System 4 Rheometer. The model describes chemoviscosity up to a range of six orders of magnitude under isothermal curing conditions. The extremely non-linear chemoviscosity profile for a dynamic heating cure cycle is predicted as well. The model is also shown to predict changes of glass transition temperature for the thermosetting resin during cure. The physical significance of this prediction is unclear at the present time, however, and further research is required. From the chemoviscosity simulation point of view, the technique of establishing an analytical model as described here is easily applied to any thermosetting resin. The model thus obtained is used in real-time process controls for fabricating composite materials.

Hou, T. H.↗

NASA NTRS1985

Contributions on the Subject of Longitudinal Movements of Aircraft in Wind Shears

The effect of downburst-type wind shears on the longitudinal dynamic behavior of an unguided aircraft is simulated numerically on the basis of published meteorological data and the flight characteristics of an A300-B passenger jet. The nonlinear differential equations of the aircraft motion are linearized by conventional methods, and the wind effects are introduced via the linear derivatives of the wind components referred to the wind gradients to obtain simplified technical models of the longitudinal response to all possible types of constant-gradient wind shears during the first 20-60 sec. Graphs, maps, and diagrams are provided, and a number of accidents presumed to have involved wind shears are analyzed in detail.

Krauspe, P.↗

NASA NTRS1985

Viscoplasticity based on overstress with a differential growth law for the equilibrium stress

Two coupled, nonlinear differential equations are proposed for the modeling of the elastic and rate (time) dependent inelastic behavior of structural metals in the absence of recovery and aging. The structure of the model is close to the unified theories but contains essential differences. It is shown that the model reproduces almost elastic regions upon initial loading and in the unloading regions of the hysteresis loop. Under loading, unloading and reloading in strain control the model simulated the experimentally observed sharp transition from nearly elastic to inelastic behavior. When a formulation akin to existing unified theories is adopted the almost elastic regions reduce the points and the transition upon reloading is very gradual. For different formulations the behavior under sudden in(de)creases of the strain rate by two orders of magnitude is simulated by numerical experiments and differences are noted. The model represents cyclically neutral behavior and contains three constants and two positive, decreasing functions. The determination of constants and functions from monotonic loading with strain rate changes and relaxation periods is described.

Krempl, E.↗

NASA NTRS1985

Relativistic electron in curved magnetic fields

Making use of the perturbation method based on the nonlinear differential equation theory, the author investigates the classical motion of a relativistic electron in a class of curved magnetic fields which may be written as B=B(O,B sub phi, O) in cylindrical coordinates (R. phi, Z). Under general astrophysical conditions the author derives the analytical expressions of the motion orbit, pitch angle, etc., of the electron in their dependence upon parameters characterizing the magnetic field and electron. The effects of non-zero curvature of magnetic field lines on the motion of electrons and applicabilities of these results to astrophysics are also discussed.

An, S.↗

NASA NTRS1984

Thermosolutal convection during directional solidification

During solidification of a binary alloy at constant velocity vertically upward, thermosolutal convection can occur if the solute rejected at the crystal-melt interface decreases the density of the melt. We assume that the crystal-melt interface remains planar and that the flow field is periodic in the horizontal direction. The time-dependent nonlinear differential equations for fluid flow, concentration, and temperature are solved numerically in two spatial dimensions for small Prandtl numbers and moderately large Schmidt numbers. For slow solidification velocities, the thermal field has an important stabilizing influence: near the onset of instability the flow is confined to the vicinity of the crystal-melt interface. Further, for slow velocities, as the concentration increases, the horizontal wavelength of the flow decreases rapidly - a phenomenon also indicated by linear stability analysis. The lateral inhomogeneity in solute concentration due to convection is obtained from the calculations. For a narrow range of solutal Rayleigh numbers and wavelengths, the flow is periodic in time.

Mcfadden, G. B.↗

NASA NTRS1983

Robustness studies in adaptive control

It has been shown that several similar Model Reference Adaptive Controllers (MRAC's) are globally stable under certain restrictive assumptions, including the assumption that the order and relative degree of the plant are exactly known, and that no disturbances are present. Under certain 'sufficiently rich' excitation conditions, the origin is globally asymptotically stable for the adaptive controller parameter error. The time-evolution of the parameter error vector represents the key factor to the stability of the adaptive controller. This time-evolution is described by a complicated set of time-varying nonlinear differential equations. The present investigation is concerned with an approximate analysis technique for studying the long term trends of the parameter error vector trajectory for an adaptive controller under periodic excitation.

Krause, J.↗

NASA NTRS1986

Flap-lag-torsional dynamics of helicopter rotor blades in forward flight

A perturbation/numerical methodology to analyze the flap-lead/lag motion of a centrally hinged spring restrained rotor blade that is valid for both hover and for forward flight was developed. The derivation of the nonlinear differential equations of motion and the analysis of the stability of the steady state response of the blade were conducted entirely in a Symbolics 3670 Machine using MACSYMA to perform all the lengthy symbolic manipulations. It also includes generation of the fortran codes and plots of the results. The Floquet theory was also applied to the differential equations of motion in order to compare results with those obtained from the perturbation analysis. The results obtained from the perturbation methodology and from Floquet theory were found to be very close to each other, which demonstrates the usefullness of the perturbation methodology. Another problem under study consisted in the analysis of the influence of higher order terms in the response and stability of a flexible rotor blade in forward flight using Computerized Symbolic Manipulation and a perturbation technique to bypass the Floquet theory. The derivation of the partial differential equations of motion is presented.

Crespodasilva, M. R. M.↗

NASA NTRS1986

Viscoplasticity based on overstress with a differential growth law for the equilibrium stress

Krempl, E.↗

NASA NTRS1986

Analysis of the harmonic gyrotron traveling wave amplifier

A single nonlinear differential equation is derived for describing the spatial evolution of the wave field in a gyrotron amplifier on the cylindrical guide geometry. The equation is then used to determine the efficiencies and the optimized interaction lengths for gyrotron operated with TE sub N1 modes at the N-th cyclotron harmonic. The equation offers a significant advantage in simplifying numerical work in comparison with the brute force particle simulation approach. The efficiencies of the harmonic operations by using a hollow beam decrease quickly with increasing harmonic number, but are enhanced by using an axis-encircling beam instead.

Kuo, S. P.↗

NASA NTRS1987

Automating the parallel processing of fluid and structural dynamics calculations

The NASA Lewis Research Center is actively involved in the development of expert system technology to assist users in applying parallel processing to computational fluid and structural dynamic analysis. The goal of this effort is to eliminate the necessity for the physical scientist to become a computer scientist in order to effectively use the computer as a research tool. Programming and operating software utilities have previously been developed to solve systems of ordinary nonlinear differential equations on parallel scalar processors. Current efforts are aimed at extending these capabilities to systems of partial differential equations, that describe the complex behavior of fluids and structures within aerospace propulsion systems. This paper presents some important considerations in the redesign, in particular, the need for algorithms and software utilities that can automatically identify data flow patterns in the application program and partition and allocate calculations to the parallel processors. A library-oriented multiprocessing concept for integrating the hardware and software functions is described.

Arpasi, Dale J.↗

NASA NTRS1986

Toward a systematic design theory for silicon solar cells using optimization techniques

This work is a first detailed attempt to systematize the design of silicon solar cells. Design principles follow from three theorems. Although the results hold only under low injection conditions in base and emitter regions, they hold for arbitrary doping profiles and include the effects of drift fields, high/low junctions and heavy doping concentrations of donor or acceptor atoms. Several optimal designs are derived from the theorems, one of which involves a three-dimensional morphology in the emitter region. The theorems are derived from a nonlinear differential equation of the Riccati form, the dependent variable of which is a normalized recombination particle current.

Misiakos, K.↗

NASA NTRS1988

On the generation of mean flows by the interaction of Goertler vortices and Tollmien-Schlichting waves in curved channel flows

There are many fluid flows where the onset of transition can be caused by different instability mechanisms which compete in the nonlinear regime. The interaction of a centrifugal instability mechanism with the viscous mechanism which causes Tollmein-Schlichting waves is discussed. The interaction between these modes can be strong enough to drive the mean state; the interaction is investigated in the context of curved channel flows so as to avoid difficulties associated with boundary layer growth. Essentially it was found that the mean state adjusts itself so that any modes present are neutrally stable even at finite amplitude. In the first instance, the mean state driven by a vortex of short wavelength in the absence of Tollmein-Schlichting waves is considered. It is shown that for a given channel curvature and vortex wavelength, there is an upper limit to the mass flow rate which the channel can support as the pressure gradient is increased. When the Tollmein-Schlichting waves are present, then the nonlinear differential equation to determine the mean state is modified.

Bassom, Andrew P.↗

NASA NTRS1987

Automating the parallel processing of fluid and structural dynamics calculations

The NASA Lewis Research Center is actively involved in the development of expert system technology to assist users in applying parallel processing to computational fluid and structural dynamic analysis. The goal of this effort is to eliminate the necessity for the physical scientist to become a computer scientist in order to effectively use the computer as a research tool. Programming and operating software utilities have previously been developed to solve systems of ordinary nonlinear differential equations on parallel scalar processors. Current efforts are aimed at extending these capabilties to systems of partial differential equations, that describe the complex behavior of fluids and structures within aerospace propulsion systems. This paper presents some important considerations in the redesign, in particular, the need for algorithms and software utilities that can automatically identify data flow patterns in the application program and partition and allocate calculations to the parallel processors. A library-oriented multiprocessing concept for integrating the hardware and software functions is described.

Arpasi, Dale J.↗

NASA NTRS1988

Rotordynamic analysis of a bearing tester

The properties of the solutions of a system of four coupled nonlinear differential equations that model the behavior of the rotating shaft of a bearing tester are studied. In particular, it is shown how the bounds for the rotations of these equations can be obtained from bounds for the solutions of the linearized equations. By studying the behavior of the Fourier transforms of the solution, the approach to the stability boundary can also be predicted. These conclusions are verified by means of numerical solutions of the equations, and of power spectrum density (PSD) plots.

Zalik, Richard A.↗

NASA NTRS1990

Plasma density enhancements created by the ionization of the Earth's upper atmosphere by artificial electron beams

Analytical calculations and experimental observations relating to the interaction with the Earth's upper atmosphere of electron beams emitted from low altitude spacecraft are presented. The problem is described by two coupled nonlinear differential equations in the up-going (along a magnetic field line) and down-going differential energy flux. The equations are solved numerically, using the MSIS atmospheric model and the IRI ionospheric model. The results form the model compare well with recent observations from the CHARGE 2 sounding rocket experiment. Two aspects of the beam-neutral atmosphere interaction are discussed. First, the limits on the electron beam current that can be emitted from a spacecraft without substantial spacecraft charging are investigated. This is important because the charging of the spacecraft to positive potentials limits the current and the escape energy of the beam electrons and thereby limits the ionization of the neutral atmosphere. As an example, we find from CHARGE 2 observations and from the model calculations that below about 180 km, secondary electrons generated through the ionization of the neutral atmosphere by 1 to 10 keV electron beams from sounding rockets, completely balance the beam current, thereby allowing the emission of very high beam currents. Second, the amount of plasma production in the beam-streak is discussed. Results are shown for selected values of the beam energy, spacecraft velocity, and spacecraft altitude.

Neubert, Torsten↗

NASA NTRS1989

Nonlinear mechanical behavior of thermoplastic matrix materials for advanced composites

Two recent theories of nonlinear mechanical response are quantitatively compared and related to experimental data. Computer techniques are formulated to handle the numerical integration and iterative procedures needed to solve the associated sets of coupled nonlinear differential equations. Problems encountered during these formulations are discussed and some open questions described. Bearing in mind these cautions, the consequences of changing parameters that appear in the formulations on the resulting engineering properties are discussed. Hence, engineering approaches to the analysis of thermoplastic matrix material can be suggested.

Arenz, R. J.↗