Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear 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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

A Constitutive Equation for Stratospheric Balloon Materials

The selection of a suitable material for use as a reliable stratospheric balloon gas barrier and structural component is based on a variety of properties. Due to a more desirable combination of properties, the low density polyethylene that has been used for the last half century has been replaced during the last decade by linear low density polyethylene (LLDPE). This paper describes the effort to characterize the time dependent properties of a 38 micron coextrusion of LLDPE. The nonlinear viscoelastic constitutive equation presented may be used to accurately describe the creep and/or relaxation of this film when subjected to a biaxial state of stress, such as might be required for an extended balloon flight. Recent laboratory data have been used to mod@ an existing model of LLDPE to account for differences caused by the coextrusion process. The new model will facilitate structural design optimization and reliability assessment, and may be further utilized as a predictive tool to benefit in-flight operations. Current structural analysis tech&ques based on linear elastic properties have predicted stresses in excess of those which would actually exist.

Rand, J. L.↗
NASA NTRS1976

Nonlinear behavior of acoustic waves in combustion chambers. I, II

The general problem of the nonlinear growth and limiting amplitude of acoustic waves in a combustion chamber is treated in three parts: (1) the general conservation equations are expanded in two small parameters, and then combined to yield a nonlinear inhomogeneous wave equation, (2) the unsteady pressure and velocity fields are expressed as a synthesis of the normal modes of the chamber, but with unknown time-varying amplitudes, and (3) the system of nonlinear equations is treated by the method of averaging to produce a set of coupled nonlinear first order differential equations for the amplitudes and phases of the modes. This approximate analysis is applied to the investigation of the unstable motions in a solid propellant rocket engine and in a T burner.

Culick, F. E. C.↗
NASA NTRS1986

Locally bound constrained Newton-Raphson solution algorithms

This paper develops strategies which enable the automatic adjustment of the constraint surfaces recently used to extend the range and numerical stability/efficiency of nonlinear finite-element equation solvers. In addition to handling kinematic and material induced nonlinearity, both pre- and postbuckling behavior can be treated. The scheme developed employs localized bounds on various hierarchical partitions of the field variables. These are used to resize, shape, and orient the global constraint surface, thereby enabling essentially automatic load/deflection incrementation. Due to the generality of the approach taken, it can be implemented in conjunction with constraints of arbitrary functional type. To benchmark the method, several numerical experiments are presented. These include problems involving kinematic and material nonlinearity, as well as, pre- and postbuckling characteristics.

Padovan, J.↗
NASA NTRS1989

An investigation of a mathematical model of an optically pumped Ti(3+):Al2O3 laser system

During the last several years, solid state lasers were developed that have the potential for meeting rigorous performance requirements for space-based remote sensing of the atmosphere. In order to design a stable and efficient laser and to understand the effect on laser output of changes in the physical and design parameters, an understanding of the development of the dynamical processes of the laser is necessary. Typically, the dynamical processes in a laser system are investigated via rate equations describing the evolution of the occupancy in the electronic levels and of the photon density in the laser cavity. There are two approaches to this type of study. Most often, for the sake of simplicity, the spatial variations of the dynamic variables in the laser system are disregarded and the mathematical model consists of a system of first order nonlinear ordinary differential equations (ODE). The second approach is to take into account both spatial and temporal variations in the dynamic variables in the laser cavity. The resulting model consists of a first order semilinear system of partial differential equations (PDE). The model which was studied was studied was generic in the sense that it was a four-level laser system, but the parameters used in the numerical study were specific to Titanium-doped sapphire. For simplicity, a constant, spatially uniform pumping scheme was considered. In addition, a simplification of the model was made so that it treats a single lasing wavelength with a narrow bandwidth. The purpose was to investigate both versions of the mathematical model and to determine whether the numerical solutions are similar both qualitatively and quantitatively. The systems of ordinary differential equations were solved numerically using a Runge-Kutta-Fehlberg algorithm which was very efficient for typical values of the physical parameters. A numerical scheme, based on the Modified Euler method, for computing solutions to the system of partial differential equations was developed and implemented. The PDE model was solved numerically at the expense of greatly increased computer time.

Roberts, Lila F.↗
NASA NTRS1993

Hypersonic flutter of a curved shallow panel with aerodynamic heating

The general equations describing the nonlinear fluttering oscillations of shallow, curved, heated orthotropic panels have been derived. The formulation takes into account the location of the panel on the surface of a generic hypersonic vehicle, when calculating the aerodynamic loads. It is also shown that third order piston theory produces unsteady aerodynamic loading which is in close agreement with that based upon direct solution of the Euler equations. Results, for simply supported panels, are obtained using Galerkin's method combined with direct numerical integration in time to compute stable limit cycle amplitudes. These results illustrate the sensitivity of the aeroelastic behavior to the unsteady aerodynamic assumptions, temperature, orthotropicity and flow orientation.

Bein, T.↗
NASA NTRS1988

Modified radiative transfer theory for a two-layer anisotropic random medium

The modified radiative transfer (MRT) equations which describe propagation and scattering of the electromagnetic field intensity in a layered anisotropic random medium are derived from the Bethe-Salpeter equation with the ladder approximation and the Dyson equation with the nonlinear approximation. The extinction matrices and scattering phase matrices for both mean and incoherent field intensities are obtained. Backscattering enhancement is observed due to the spatial coherent effect between upward and downward propagating waves. It also occurs for the half-space case, because of coupling between ordinary and extraordinary waves. The depolarization effect is predicted in the first-order renormalization approximation to the MRT equation.

Lee, J. K.↗
NASA NTRS1989

Nonlinear interactions in mixing layers and compressible heated round jets

The nonlinear interactions between a fundamental instability mode and both its harmonics and the changing mean flow are studied using the weakly nonlinear stability theory of Stuart and Watson, and numerical solutions of coupled nonlinear partial differential equations. The first part focuses on incompressible cold (or isothermal; constant temperature throughout) mixing layers, and for these, the first and second Landau constants are calculated as functions of wavenumber and Reynolds number. It is found that the dominant contribution to the Landau constants arises from the mean flow changes and not from the higher harmonics. In order to establish the range of validity of the weakly nonlinear theory, the weakly nonlinear and numerical solutions are compared and the limitation of each is discussed. At small amplitudes and at low-to-moderate Reynolds numbers, the two results compare well in describing the saturation of the fundamental, the distortion of the mean flow, and the initial stages of vorticity roll-up. At larger amplitudes, the interaction between the fundamental, second harmonic, and the mean flow is strongly nonlinear and the numerical solution predicts flow oscillations, whereas the weakly nonlinear theory yields saturation. In the second part, the weakly nonlinear theory is extended to heated (or nonisothermal; mean temperature distribution) subsonic round jets where quadratic and cubic nonlinear interactions are present, and the Landau constants also depend on jet temperature ratio, Mach number and azimuthal mode number. Under exponential growth and nonlinear saturation, it is found that heating and compressibility suppress the growth of instability waves, that the first azimuthal mode is the dominant instability mode, and that the weakly nonlinear solution describes the early stages of the roll-up of an axisymmetric shear layer. The receptivity of a typical jet flow to pulse type input disturbance is also studied by solving the initial value problem and then examining the behavior of the long-time solution.

Jarrah, Yousef Mohd↗
NASA NTRS1994

Crystal growth and fluid mechanics problems in directional solidification

Broadly speaking, our efforts have been concentrated in two aspects of directional solidification: (A) a more complete theoretical understanding of convection effects in a Bridgman apparatus; and (B) a clear understanding of scalings of various features of dendritic crystal growth in the sensitive limit of small capillary effects. For studies that fall within class A, the principal objectives are as follows: (A1) Derive analytical formulas for segregation, interfacial shape and fluid velocities in mathematically amenable asymptotic limits. (A2) Numerically verify and extend asymptotic results to other ranges of parameter space with a view to a broader physical understanding of the general trends. With respect to studies that fall within class B, the principal objectives include answering the following questions about dendritic crystal growth: (B1) Are there unsteady dendrite solutions in 2-D to the completely nonlinear time evolving equations in the small surface tension limit with only a locally steady tip region with well defined tip radius and velocity? Is anisotropy in surface tension necessary for the existence of such solutions as it is for a true steady state needle crystal? How does the size of such a local region depend on capillary effects, anisotropy and undercooling? (B2) How do the different control parameters affect the nonlinear amplification of tip noise and dendritic side branch coarsening?

Tanveer, Saleh↗
NASA NTRS1990

Implications of chaos, scale-invariance, and fractal statistics in geology

A set of three nonlinear total differential equations (Lorenz equations) exhibiting deterministic chaos is considered, and it is shown that these equations demonstrate that deterministic equations with deterministic initial conditions can yield stocastic solutions with fractal statistics. The logistic map, fractal distributions, and fragmentation are discussed. It is pointed out that well-defined fractal distributions of earthquakes are found both regionally and globally, and that the general applicability of the fractal relation for seismicity can provide the basis for a quantitative seismic hazard assessment. It is suggested that the governing physics of erosional topography is nonlinear and may be related to a fractal distribution of storms and floods that generate and renew erosional feature such as gullies and drainage systems.

Turcotte, D. L.↗
NASA NTRS1973

Vertical distribution of vibrational energy of molecular nitrogen in a stable auroral red arc and its effect on ionospheric electron densities

Previous solutions of the problem of the distribution of vibrationally excited molecular nitrogen in the thermosphere have either assumed a Boltzmann distribution and considered diffusion as one of the loss processes or solved for the energy level populations and neglected diffusion. Both of the previous approaches are combined by solving the time dependent continuity equations, including the diffusion process, for the first six energy levels of molecular nitrogen for conditions in the thermosphere corresponding to a stable auroral red arc. The primary source of molecular nitrogen excitation was subexcitation, and inelastic collisions between thermal electrons and molecular nitrogen. The reaction rates for this process were calculated from published cross section calculations. The loss processes for vibrational energy were electron and atomic oxygen quenching and vibrational energy exchange. The coupled sets of nonlinear, partial differential equations were solved numerically by employing finite difference equations.

Newton, G. P.↗
NASA NTRS1996

Low-Dimensional Dynamical Models of Thermal Convection

A low-dimensional dynamic model for transitional buoyancy-driven flow in a differentially heated tall enclosure is presented. The full governing partial differential equations with the associated boundary conditions are solved by a spectral element method for a cavity of aspect ratio A=20. Proper orthogonal decomposition is applied to the oscillatory solution at Prandtl number Pr=P tau (omega) = 0.71 and Grashof number G tau (omega) = 3.2 x 10 (exp 4) to construct empirical eigenfunctions. Using the four most energetic empirical eigenfunctions for the velocity and temperature as basis functions and applying Galerkin's method, a reduced model consisting of eight nonlinear ordinary differential equations is obtained. Close to the 'design' conditions (P tau(omega) G tau(omega)), the low-order model (LOM) predictions are in excellent agreement with the predictions of the full model. In particular, the critical Grashof number at the onset of the first temporal flow instability (Hopf bifurcation) was well as the frequency and amplitude of oscillations at supercritical conditions are in excellent agreement with the predictions of the full model. Far from the 'design' conditions, the LOM predicts the existence of multiple stable steady solutions at large values of G tau, and a unique stable steady solution at small values of G tau, and exhibits hysteretic behavior that is qualitatively similar to that observed in direct numerical simulations based on the full model.

Liakopoulos, Anthony↗
NASA NTRS1975

Stability of localized plasma model in two and three dimensions

Stability against catastrophic radial collapse is demonstrated for initially well-behaved, localized solutions in a class of two- and three-dimensional wave equations characterized by saturating nonlinearities, which model various electrostatic and electromagnetic 'caviton' structures in an asymptotically uniform plasma. This result contrasts sharply with mathematical predictions of equations with low-order nonlinearities. Discrete classes of possible quasistationary modes exhibit features comparable to atomic wave functions and differ qualitatively from one-dimensional structures.

Wilcox, J. Z.↗
NASA NTRS1972

Lagrangian methods in the analysis of nonlinear wave interactions in plasma

An averaged-Lagrangian method is developed for obtaining the equations which describe the nonlinear interactions of the wave (oscillatory) and background (nonoscillatory) components which comprise a continuous medium. The method applies to monochromatic waves in any continuous medium that can be described by a Lagrangian density, but is demonstrated in the context of plasma physics. The theory is presented in a more general and unified form by way of a new averaged-Lagrangian formalism which simplifies the perturbation ordering procedure. Earlier theory is extended to deal with a medium distributed in velocity space and to account for the interaction of the background with the waves. The analytic steps are systematized, so as to maximize calculational efficiency. An assessment of the applicability and limitations of the method shows that it has some definite advantages over other approaches in efficiency and versatility.

Galloway, J. J.↗
NASA NTRS1998

A Discontinuous Galerkin Finite Element Method for Hamilton-Jacobi Equations

In this paper, we present a discontinuous Galerkin finite element method for solving the nonlinear Hamilton-Jacobi equations. This method is based on the Runge-Kutta discontinuous Galerkin finite element method for solving conservation laws. The method has the flexibility of treating complicated geometry by using arbitrary triangulation, can achieve high order accuracy with a local, compact stencil, and are suited for efficient parallel implementation. One and two dimensional numerical examples are given to illustrate the capability of the method.

Hu, Changqing↗
NASA NTRS1974

Cosmic ray scintillations. II - General theory of interplanetary scintillations

The motion of charged particles in a stochastic magnetic field with nonzero mean is considered via a generalized quasi-linear expansion of Liouville's equation. The general result is an equation relating cosmic ray scintillations to magnetic fluctuations and to cosmic ray gradients. The resonant interaction between particles and the random magnetic field is considered in detail, and the effect of nonlinear terms in the equations is considered. The nonlinear terms are important in damping out initial conditions and in determining conditions near cyclotron resonances. The application of the theory to the propagation of cosmic rays during quiet times in interplanetary space is considered. It is concluded that cosmic ray scintillations in interplanetary space may provide useful information about interplanetary particles and fields and also about nonlinear plasma interactions.

Owens, A. J.↗
NASA NTRS1975

Gravitational collapse of a turbulent vortex with application to star formation

The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitation and rotation terms in the equations are of first order in the space variables, the pressure gradient terms are of second order, and the turbulent viscosity term is of third order. The presence of a turbulent viscosity insures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the initial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial inflow plays an important role and allows collapse to occur even when the rotation is large. An approximate solution of the governing partial differential equations is also given; the equations are used to study the spacial distributions of the density and velocity.

Deissler, R. G.↗
NASA NTRS1975

Gravitational collapse of a turbulent vortex with application to star formation

The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitation and rotation terms in the equations are of first order in the space variables, the pressure gradient terms are of second order, and the turbulent viscosity term is of third order. The presence of a turbulent viscosity insures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the intial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial inflow plays an important role and allows collapse to occur even when the rotation is large. An approximate solution of the governing partial differential equations is also given, in order to study the spacial distributions of the density and velocity.

Deissler, R. G.↗
NASA NTRS1976

Gravitational collapse of a turbulent vortex with application to star formation

The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitational and rotational terms in the equations are of first order in the space variables, the pressure-gradient terms are of second order, and the turbulent-viscosity term is of third order. The presence of turbulent viscosity ensures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the initial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial flow is always inward and allows collapse to occur (axially) even when the rotation is large. An approximate solution of the governing partial differential equations is also given in order to study the spatial distributions of the density and velocity.

Deissler, R. G.↗