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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Computational modeling of femtosecond optical solitons from Maxwell's equations

An algorithm is developed that permits the direct time integration of full-vector nonlinear Maxwell's equations. This capability permits the modeling of both linear and nonlinear instantaneous and dispersive effects in the electric polarization in material media. The modeling of the optical carrier is retained. The fundamental innovation is to notice that it is possible to treat the linear and nonlinear convolution integrals, which describe the dispersion, as new dependent variables. A coupled system of nonlinear second-order ordinary differential equations can then be derived for the linear and nonlinear convolution integrals, by differentiating them in the time domain. These equations, together with Maxwell's equations, are solved to determine the electromagnetic fields in nonlinear dispersive media. Results are presented of calculations in one dimension of the propagation and collision of femtosecond electromagnetic solitons that retain the optical carrier, taking into account as the Kerr and Raman interactions.

Goorjian, Peter M.↗
NASA NTRS1988

A multivariate variational objective analysis-assimilation method. Part 1: Development of the basic model

The variational method of undetermined multipliers is used to derive a multivariate model for objective analysis. The model is intended for the assimilation of 3-D fields of rawinsonde height, temperature and wind, and mean level temperature observed by satellite into a dynamically consistent data set. Relative measurement errors are taken into account. The dynamic equations are the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation. The model Euler-Lagrange equations are eleven linear and/or nonlinear partial differential and/or algebraic equations. A cyclical solution sequence is described. Other model features include a nonlinear terrain-following vertical coordinate that eliminates truncation error in the pressure gradient terms of the horizontal momentum equations and easily accommodates satellite observed mean layer temperatures in the middle and upper troposphere. A projection of the pressure gradient onto equivalent pressure surfaces removes most of the adverse impacts of the lower coordinate surface on the variational adjustment.

Achtemeier, Gary L.↗
NASA NTRS1972

A numerical investigation of film cooling as an alternate thermal protection system for spacecraft during atmospheric reentry

Research was undertaken to determine the mass required to protect a surface during atmospheric reentry by injecting a cold gas tangentially along the surface. A consideration of this problem revealed that the boundary layer equations could be employed to approximate the most important flow configurations which occur in film cooling. However, even with this simplification, the equations to be solved are complex, nonlinear, coupled partial differential equations. Due to the relative weakness of classical methods of boundary layer analysis, finite difference methods were chosen as the best approach to the formidable problem under consideration. Instead of blindly adopting a proposed technique and trying to force it to the problem, a more systematic study was deemed to be not only desirable but necessary. A relatively basic study of the application of finite difference methods was undertaken with the objective of developing an algorithm capable of obtaining accurate solutions to film cooling problems.

Clausing, A. M.↗
NASA NTRS1973

Nonlinear response of plates subjected to inplane and lateral pressure pulses.

The nonlinear response of a rectangular plate exposed to a far-field sonic boom disturbance is studied. The plate is subjected to both lateral and in-plane disturbances. The lateral disturbance is in the form of an N-shaped pressure pulse, and the in-plane disturbance is represented by a sinusoidal pulse. The equations of motion are reduced to a set of nonlinear coupled ordinary differential equations using Galerkin's method. These equations are solved numerically using Hamming's (1959) modified predictor-corrector integration method. The effects of in-plane boundary conditions and in-plane inertia are investigated. The nonlinear results, when compared to the linear theory, serve to delineate the realm of validity of the linear theory.

Knapp, L. J.↗
NASA NTRS1979

Stability analysis of a liquid fuel annular combustion chamber

The problems of combustion instability in an annular combustion chamber are investigated. A modified Galerkin method was used to produce a set of modal amplitude equations from the general nonlinear partial differential acoustic wave equation. From these modal amplitude equations, the two variable perturbation method was used to develop a set of approximate equations of a given order of magnitude. These equations were modeled to show the effects of velocity sensitive combustion instabilities by evaluating the effects of certain parameters in the given set of equations. By evaluating these effects, parameters which cause instabilities to occur in the combustion chamber can be ascertained. It is assumed that in the annular combustion chamber, the liquid propellants are injected uniformly across the injector face, the combustion processes are distributed throughout the combustion chamber, and that no time delay occurs in the combustion processes.

Mcdonald, G. H.↗
NASA NTRS1978

Stability analysis of a liquid fuel annular combustion chamber

High frequency combustion instability problems in a liquid fuel annular combustion chamber are examined. A modified Galerkin method was used to produce a set of modal amplitude equations from the general nonlinear partial differential acoustic wave equation in order to analyze the problem of instability. From these modal amplitude equations, the two variable perturbation method was used to develop a set of approximate equations of a given order of magnitude. These equations were modeled to show the effects of velocity sensitive combustion instabilities by evaluating the effects of certain parameters in the given set of equations.

Mcdonald, G. H.↗
NASA NTRS1992

Comparison of level-1 and level-2 buckling and postbuckling solutions

The accuracy of the solutions obtained with the Level-1 computational modules of DISDECO (Delft Interactive Shell Design Code) are investigated by comparing them with solutions obtained with ANILISA, a Level-2 computational module. Whereas the Level-1 solutions are based on the use of trigonometric functions to represent the out-of-plane displacement W and involve Galerkin type approximations to reduce the governing partial differential equations to simple algebraic expressions, the Level-2 solutions employ a truncated Fourier expansion in the circumferential direction and solve the resulting ordinary differential equations forming a nonlinear eigenvalue problem by a Stodola like technique. This involves the numerical integration of the set of ordinary differential equations, which in turn makes it possible to satisfy the specified boundary conditions exactly. It is shown that in many applications the Level-1 predictions concerning the critical buckling load of the perfect structure and its imperfection sensitivity are quite accurate and reliable. For those cases where this is not true, it is shown how the pre and postprocessing facilities of DISDECO can be used efficiently to arrive at a reliable prediction. This involves always the use of the Level-2 computational module ANILISA with its rigorous prebuckling solution and its feature which makes it possible to satisfy the specified boundary conditions exactly.

Arbocz, J.↗
NASA NTRS1993

The rollup of a vortex layer near a wall

The behavior of an inviscid vortex layer of non-zero thickness near a wall is studied, both through direct numerical simulation of the two-dimensional vorticity equation at high Reynolds numbers, and using an approximate ordinary nonlinear integro-differential equation which is satisfied in the limit of a thin layer under long-wavelength perturbations. For appropriate initial conditions the layer rolls up and breaks into compact vortices which move along the wall at constant speed. Because of the effect of the wall, they correspond to equilibrium counter-rotating vortex dipoles. This breakup can be related to the disintegration of the initial conditions of the approximate nonlinear dispersive equation into solitary waves. The study is motivated by the formation of longitudinal vortices from vortex sheets in the wall region of a turbulent channel.

Jimenez, Javier↗
NASA NTRS1988

Estimation of time- and state-dependent delays and other parameters in functional differential equations

A parameter estimation algorithm is developed which can be used to estimate unknown time- or state-dependent delays and other parameters (e.g., initial condition) appearing within a nonlinear nonautonomous functional differential equation. The original infinite dimensional differential equation is approximated using linear splines, which are allowed to move with the variable delay. The variable delays are approximated using linear splines as well. The approximation scheme produces a system of ordinary differential equations with nice computational properties. The unknown parameters are estimated within the approximating systems by minimizing a least-squares fit-to-data criterion. Convergence theorems are proved for time-dependent delays and state-dependent delays within two classes, which say essentially that fitting the data by using approximations will, in the limit, provide a fit to the data using the original system. Numerical test examples are presented which illustrate the method for all types of delay.

Murphy, K. A.↗
NASA NTRS1990

Estimation of time- and state-dependent delays and other parameters in functional differential equations

A parameter estimation algorithm is developed which can be used to estimate unknown time- or state-dependent delays and other parameters (e.g., initial condition) appearing within a nonlinear nonautonomous functional differential equation. The original infinite dimensional differential equation is approximated using linear splines, which are allowed to move with the variable delay. The variable delays are approximated using linear splines as well. The approximation scheme produces a system of ordinary differential equations with nice computational properties. The unknown parameters are estimated within the approximating systems by minimizing a least-squares fit-to-data criterion. Convergence theorems are proved for time-dependent delays and state-dependent delays within two classes, which say essentially that fitting the data by using approximations will, in the limit, provide a fit to the data using the original system. Numerical test examples are presented which illustrate the method for all types of delay.

Murphy, K. A.↗
NASA NTRS1993

On the dynamics of approximating schemes for dissipative nonlinear equations

Since one can rarely write down the analytical solutions to nonlinear dissipative partial differential equations (PDE's), it is important to understand whether, and in what sense, the behavior of approximating schemes to these equations reflects the true dynamics of the original equations. Further, because standard error estimates between approximations of the true solutions coming from spectral methods - finite difference or finite element schemes, for example - and the exact solutions grow exponentially in time, this analysis provides little value in understanding the infinite time behavior of a given approximating scheme. The notion of the global attractor has been useful in quantifying the infinite time behavior of dissipative PDEs, such as the Navier-Stokes equations. Loosely speaking, the global attractor is all that remains of a sufficiently large bounded set in phase space mapped infinitely forward in time under the evolution of the PDE. Though the attractor has been shown to have some nice properties - it is compact, connected, and finite dimensional, for example - it is in general quite complicated. Nevertheless, the global attractor gives a way to understand how the infinite time behavior of approximating schemes such as the ones coming from a finite difference, finite element, or spectral method relates to that of the original PDE. Indeed, one can often show that such approximations also have a global attractor. We therefore only need to understand how the structure of the attractor for the PDE behaves under approximation. This is by no means a trivial task. Several interesting results have been obtained in this direction. However, we will not go into the details. We mention here that approximations generally lose information about the system no matter how accurate they are. There are examples that show certain parts of the attractor may be lost by arbitrary small perturbations of the original equations.

Jones, Donald A.↗
NASA NTRS1993

The nonlinear evolution of modes on unstable stratified shear layers

The nonlinear development of disturbances in stratified shear flows (having a local Richardson number of value less than one quarter) is considered. Such modes are initially fast growing but, like related studies, we assume that the viscous, non-parallel spreading of the shear layer results in them evolving in a linear fashion until they reach a position where their amplitudes are large enough and their growth rates have diminished sufficiently so that amplitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer theories. Four different basic integro-differential amplitude equations are possible, including one due to a novel mechanism; the relevant choice of amplitude equation, at a particular instance, being dependent on the relative sizes of the disturbance amplitude, the growth rate of the disturbance, its wavenumber, and the viscosity of the fluid. This richness of choice of possible nonlinearities arises mathematically from the indicial Frobenius roots of the governing linear inviscid equation (the Taylor-Goldstein equation) not, in general, differing by an integer. The initial nonlinear evolution of a mode will be governed by an integro-differential amplitude equations with a cubic nonlinearity but the resulting significant increase in the size of the disturbance's amplitude leads on to the next stage of the evolution process where the evolution of the mode is governed by an integro-differential amplitude equations with a quintic nonlinearity. Continued growth of the disturbance amplitude is expected during this stage, resulting in the effects of nonlinearity spreading to outside the critical level, by which time the flow has become fully nonlinear.

Blackaby, Nicholas↗
NASA NTRS1992

Discrete time learning control in nonlinear systems

In this paper digital learning control methods are developed primarily for use in single-input, single-output nonlinear dynamic systems. Conditions for convergence of the basic form of learning control based on integral control concepts are given, and shown to be satisfied by a large class of nonlinear problems. It is shown that it is not the gross nonlinearities of the differential equations that matter in the convergence, but rather the much smaller nonlinearities that can manifest themselves during the short time interval of one sample time. New algorithms are developed that eliminate restrictions on the size of the learning gain, and on knowledge of the appropriate sign of the learning gain, for convergence to zero error in tracking a feasible desired output trajectory. It is shown that one of the new algorithms can give guaranteed convergence in the presence of actuator saturation constraints, and indicate when the requested trajectory is beyond the actuator capabilities.

Longman, Richard W.↗
NASA NTRS1987

Parallels between control PDE's (Partial Differential Equations) and systems of ODE's (Ordinary Differential Equations)

System theorists understand that the same mathematical objects which determine controllability for nonlinear control systems of ordinary differential equations (ODEs) also determine hypoellipticity for linear partial differentail equations (PDEs). Moreover, almost any study of ODE systems begins with linear systems. It is remarkable that Hormander's paper on hypoellipticity of second order linear p.d.e.'s starts with equations due to Kolmogorov, which are shown to be analogous to the linear PDEs. Eigenvalue placement by state feedback for a controllable linear system can be paralleled for a Kolmogorov equation if an appropriate type of feedback is introduced. Results concerning transformations of nonlinear systems to linear systems are similar to results for transforming a linear PDE to a Kolmogorov equation.

Hunt, L. R.↗
NASA NTRS1983

Postbuckling of long orthotropic plates in combined shear and compression

The nonlinear large-deflection partial differential equations of von karman for orthotropic plates loaded in combined shear and compression are converted into a set of first-order nonlinear ordinary differential equations by assuming trigonometric functions in one direction. These equations are solved numerically using a two point boundary problem solver which makes use of Newton's method. Results are obtained which determine the postbuckling behavior of rectangular plates with loading up to about three times the buckling load. Both isotropic and orthotropic composite plates are considered. Results show that orthotropic plates may behave quite differently than isotropic plates and that in-plane boundary conditions are important for plates loaded in shear.

Stein, M.↗