Engineering Papers⌕ Search

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.

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

At least 163 records · Page 9

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

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

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

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

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

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

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

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

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

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

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

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↗

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

Rotordynamics analysis of a Jeffcott model with deadband

A method is developed for determining the stability margins of a simple Jeffcott model with deadband via analysis of the discrete Fourier transform of the system response. The model in question is of a uniform, unbalanced, flexible shaft that is supported by a bearing as it rotates about its x axis. This model is represented by a system of coupled nonlinear differential equations.

Zalik, R. A.↗

A multiple scales approach to sound generation by vibrating bodies

The problem of determining the acoustic field in an inviscid, isentropic fluid generated by a solid body whose surface executes prescribed vibrations is formulated and solved as a multiple scales perturbation problem, using the Mach number M based on the maximum surface velocity as the perturbation parameter. Following the idea of multiple scales, new 'slow' spacial scales are introduced, which are defined as the usual physical spacial scale multiplied by powers of M. The governing nonlinear differential equations lead to a sequence of linear problems for the perturbation coefficient functions. However, it is shown that the higher order perturbation functions obtained in this manner will dominate the lower order solutions unless their dependence on the slow spacial scales is chosen in a certain manner. In particular, it is shown that the perturbation functions must satisfy an equation similar to Burgers' equation, with a slow spacial scale playing the role of the time-like variable. The method is illustrated by a simple one-dimenstional example, as well as by three different cases of a vibrating sphere. The results are compared with solutions obtained by purely numerical methods and some insights provided by the perturbation approach are discussed.

Geer, James F.↗

Dynamics of aging magnetic clouds

The dynamics of radially expanding magnetic clouds is rigorously analyzed within the framework of ideal MHD. The cloud is modelled as a cylindrically symmetric magnetic flux rope. In the force balance we include the gas pressure gradient and the Lorentz force. Interaction with the ambient solar wind due to expansion of the magnetic cloud is represented by a drag force proportional to the bulk velocity. We consider the self-similar expansion of a polytrope, and reduce the problem to an ordinary nonlinear differential equation for the evolution function. Analyzing the asymptotic behavior of the evolution function, we formulate theoretical expectations for the long-term behavior of cloud parameters. We focus on the temporal evolution of (1) the magnetic field strength; (2) the twist of the field lines; (3) the asymmetry of the total field profile; and (4) the bulk flow speed. We present data from two magnetic clouds observed at 1 AU and 2 AU, respectively, and find good agreement with theoretical expectations. For a peak magnetic field strength at 1 AU of 25 nT and a polytropic index of 0.5, we find that a magnetic cloud can be distinguished from the background interplanetary field up to a distance of about 5 AU. Taking larger magnetic fields and bigger polytropic indices this distance can double.

Osherovich, V. A.↗

Two-dimensional compressible flow in turbomachines with conic flow surfaces

A general method of analysis is developed for two-dimensional, steady, compressible flow in stators or rotors of radial and mixed flow turbomachines with conic flow surfaces (surfaces of right circular cones generated by center line of flow passage in the axial-radial plane). The variables taken into account are: (1) tip speed of the rotor, (2) flow rate, (3) blade shape, (4) variation in passage height with radius, (5) number of blades, and (6) cone angle of the flow surface. Relaxation methods are used to solve the nonlinear differential equation for the stream function. Two numerical examples are presented; one for compressible and the other for incompressible flow in a centrifugal compressor with thin, straight blades. The results of these examples are given by plots of the streamlines, constant velocity-ratio lines, and constant pressure-ratio line.

Stanitz, John D↗

On transonic flow past a wave-shaped wall

This report is an extension of a previous investigation (described in NACA rep. 1069) concerned with the solution of the nonlinear differential equation for transonic flow past a wavy wall. In the present work several new notions are introduced which permit the solution of the recursion formulas arising from the method of integration in series. In addition, a novel numerical tests of convergence, applied to the power series (in transonic similarity parameter) representing the local Mach number distribution at the boundary, indicates that smooth symmetrical potential flow past the wavy wall is no longer possible once the critical value of the stream Mach number has been exceeded.

Kaplan, Carl↗