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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

A hybrid perturbation-Galerkin technique that combines multiple expansions

A two-step hybrid perturbation-Galerkin method for the solution of a variety of differential equations type problems is found to give better results when multiple perturbation expansions are employed. The method assumes that there is parameter in the problem formulation and that a perturbation method can be used to construct one or more expansions in this perturbation coefficient functions multiplied by computed amplitudes. In step one, regular and/or singular perturbation methods are used to determine the perturbation coefficient functions. The results of step one are in the form of one or more expansions each expressed as a sum of perturbation coefficient functions multiplied by a priori known gauge functions. In step two the classical Bubnov-Galerkin method uses the perturbation coefficient functions computed in step one to determine a set of amplitudes which replace and improve upon the gauge functions. The hybrid method has the potential of overcoming some of the drawbacks of the perturbation and Galerkin methods as applied separately, while combining some of their better features. The proposed method is applied, with two perturbation expansions in each case, to a variety of model ordinary differential equations problems including: a family of linear two-boundary-value problems, a nonlinear two-point boundary-value problem, a quantum mechanical eigenvalue problem and a nonlinear free oscillation problem. The results obtained from the hybrid methods are compared with approximate solutions obtained by other methods, and the applicability of the hybrid method to broader problem areas is discussed.

Geer, James F.↗
NASA NTRS2002

The Statistical Mechanics of Ideal Homogeneous Turbulence

Plasmas, such as those found in the space environment or in plasma confinement devices, are often modeled as electrically conducting fluids. When fluids and plasmas are energetically stirred, regions of highly nonlinear, chaotic behavior known as turbulence arise. Understanding the fundamental nature of turbulence is a long-standing theoretical challenge. The present work describes a statistical theory concerning a certain class of nonlinear, finite dimensional, dynamical models of turbulence. These models arise when the partial differential equations describing incompressible, ideal (i.e., nondissipative) homogeneous fluid and magnetofluid (i.e., plasma) turbulence are Fourier transformed into a very large set of ordinary differential equations. These equations define a divergenceless flow in a high-dimensional phase space, which allows for the existence of a Liouville theorem, guaranteeing a distribution function based on constants of the motion (integral invariants). The novelty of these particular dynamical systems is that there are integral invariants other than the energy, and that some of these invariants behave like pseudoscalars under two of the discrete symmetry transformations of physics, parity, and charge conjugation. In this work the 'rugged invariants' of ideal homogeneous turbulence are shown to be the only significant scalar and pseudoscalar invariants. The discovery that pseudoscalar invariants cause symmetries of the original equations to be dynamically broken and induce a nonergodic structure on the associated phase space is the primary result presented here. Applicability of this result to dissipative turbulence is also discussed.

Shebalin, John V.↗
NASA NTRS2018

Curvilinear Displacement Transfer Functions for Deformed Shape Predictions of Curved Structures Using Distributed Surface Strains

Curvilinear Displacement Transfer Functions were formulated for deformed shape predictions of different curved structures using surface strains. The embedded curved beam (depth-wise cross section of a curved structure along a surface strain-sensing line) was discretized into multiple small domains, with domain junctures matching the strain-sensing stations. Thus, the surface strain distribution can be described with a piecewise linear or a piecewise nonlinear function. The discrete approach enabled piecewise integrations of a curvature-strain differential equation for the embedded curved beam to yield closed-form Curvilinear Displacement Transfer Functions, which are written in terms of embedded curved-beam geometrical parameters and surface strains. By inputting the surface strain data, the Curvilinear Displacement Transfer Functions can transform surface strains into deflections along each embedded curved beam for mapping out the overall structural deformed shapes. The finite-element method was used to analytically generate the surface strains of the curved beams. The deformed shape prediction accuracies were then determined by comparing the theoretical deflections with the finite-element-generated deflections, which were used as yardsticks. By introducing the correction factors in simple mathematical forms, the Curvilinear Displacement Transfer Functions can be quite accurate for shape predictions of different curved-beam structures ranging from limit case of straight beam up to semicircular curved beam.

Ko, William L.↗
NASA NTRS1973

Program for the solution of multipoint boundary value problems of quasilinear differential equations

Linear equations are solved by a method of superposition of solutions of a sequence of initial value problems. For nonlinear equations and/or boundary conditions, the solution is iterative and in each iteration a problem like the linear case is solved. A simple Taylor series expansion is used for the linearization of both nonlinear equations and nonlinear boundary conditions. The perturbation method of solution is used in preference to quasilinearization because of programming ease, and smaller storage requirements; and experiments indicate that the desired convergence properties exist although no proof or convergence is given.

Source record↗
NASA NTRS1995

On the Nonlinear Evolution of a Stationary Cross-Flow Vortex in a Fully Three-Dimensional Boundary Layer Flow

We consider the nonlinear stability of a fully three-dimensional boundary layer flow in an incompressible fluid and derive an equation governing the nonlinear development of a stationary cross-flow vortex. The amplitude equation is a novel integro-differential equation which has spatial derivatives of the amplitude occurring in the kernal function. It is shown that the evolution of the cross-flow vortex is strongly coupled to the properties of an unsteady wall layer which is in fact driven by an unknown slip velocity, proportional to the amplitude of the cross-flow vortex. The work is extended to obtain the corresponding equation for rotating disk flow. A number of special cases are examined and the numerical solution for one of cases, and further analysis, demonstrates the existence of finite-distance as well as focussing type singularities. The numerical solutions also indicate the presence of a new type of nonlinear wave solution for a certain set of parameter values.

Gajjar, J. S. B.↗
NASA NTRS1971

Generalized invariance principles and the theory of stability.

Description of some recent extensions of the invariance principle to more generalized dynamical systems where the state space is not locally compact and the flow is unique only in the forward direction of time. A sufficient condition for asymptotic stability of an invariant set is obtained which does not require that the Liapunov function be positive-definite. A recently developed generalized invariance principle is described which is applicable to functional differential equations, partial differential equations, and, in particular, to certain stability problems arising in thermoelasticity, viscoelasticity, and distributed nonlinear networks.

Lasalle, J. P.↗
NASA NTRS1984

Nonlinear control of aircraft

Transformations of nonlinear systems were used to design automatic flight controllers for vertical and short takeoff aircraft. Under the assumption that a nonlinear system can be mapped to a controllable linear system, a method using partial differential equations was constructed to approximate transformations in cases where exact ones cannot be found. An application of the design theory to a rotorcraft, the UH-1H helicopter, was presented.

Hunt, L. R.↗
NASA NTRS1976

A strictly Markovian expansion for plasma turbulence theory

The collision operator that appears in the equation of motion for a particle distribution function that was averaged over an ensemble of random Hamiltonians is non-Markovian. It is non-Markovian in that it involves a propagated integral over the past history of the ensemble averaged distribution function. All formal expansions of this nonlinear collision operator to date preserve this non-Markovian character term by term yielding an integro-differential equation that must be converted to a diffusion equation by an additional approximation. An expansion is derived for the collision operator that is strictly Markovian to any finite order and yields a diffusion equation as the lowest nontrivial order. The validity of this expansion is seen to be the same as that of the standard quasilinear expansion.

Jones, F. C.↗
NASA NTRS1978

A strictly Markovian expansion for plasma turbulence theory

The collision operator that appears in the equation of motion for a particle distribution function that has been averaged over an ensemble of random Hamiltonians is non-Markovian. It is non-Markovian in that it involves a propagated integral over the past history of the ensemble averaged distribution function. All formal expansions of this nonlinear collision operator to date preserve this non-Markovian character term by term yielding an integro-differential equation that must be converted to a diffusion equation by an additional approximation. In this note we derive an expansion of the collision operator that is strictly Markovian to any finite order and yields a diffusion equation as the lowest non-trivial order. The validity of this expansion is seen to be the same as that of the standard quasi-linear expansion.

Jones, F. C.↗
NASA NTRS1974

Inelastic, nonlinear analysis of stiffened shells of revolution by numerical integration

This paper describes the latest addition to the STARS system of computer programs, STARS-2P, for the plastic, large deflection analysis of axisymmetrically loaded shells of revolution. The STARS system uses a numerical integration scheme to solve the governing differential equations. Several unique features for shell of revolution programs that are included in the STARS-2P program are described. These include orthotropic nonlinear kinematic hardening theory, a variety of shell wall cross sections and discrete ring stiffeners, cyclic and nonproportional mechanical and thermal loading capability, the coupled axisymmetric large deflection elasto-plastic torsion problem, an extensive restart option, arbitrary branching capability, and the provision for the inelastic treatment of smeared stiffeners, isogrid, and waffle wall constructions. To affirm the validity of the results, comparisons with available theoretical and experimental data are presented.

Levine, H. S.↗
NASA NTRS1973

A square-root data array solution of the continuous-discrete filtering problem

The Dyer-McReynolds (1969) discrete square-root filtering algorithm is extended to accommodate continuous dynamics. Differential equations are given to represent the time evolution of the filter data array. These equations are nonlinear, but it is shown that the nonlinearities act to enhance the stability of the solution.

Bierman, G. J.↗
NASA NTRS1978

The practical application of a finite difference method for analyzing transonic flow over oscillating airfoils and wings

Separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady equations for small disturbances was performed. The steady velocity potential was obtained first from the well known nonlinear equation for steady transonic flow. The unsteady velocity potential was then obtained from a linear differential equation in complex form with spatially varying coefficients. Since sinusoidal motion is assumed, the unsteady equation is independent of time. The results of an investigation into the relaxation-solution-instability problem was discussed. Concepts examined include variations in outer boundary conditions, a coordinate transformation so that the boundary condition at infinity may be applied to the outer boundaries of the finite difference region, and overlapping subregions. The general conclusion was that only a full direct solution in which all unknowns are obtained at the same time will avoid the solution instabilities of relaxation. An analysis of the one-dimensional form of the unsteady transonic equation was studied to evaluate errors between exact and finite difference solutions. Pressure distributions were presented for a low-aspect-ratio clipped delta wing at Mach number of 0.9 and for a moderate-aspect-ratio rectangular wing at a Mach number of 0.875.

Weatherill, W. H.↗
NASA NTRS1980

Further investigation of a finite difference procedure for analyzing the transonic flow about harmonically oscillating airfoils and wings

Analytical and empirical studies of a finite difference method for the solution of the transonic flow about harmonically oscillating wings and airfoils are presented. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady equations for small disturbances. The steady velocity potential is obtained first from the well-known nonlinear equation for steady transonic flow. The unsteady velocity potential is then obtained from a linear differential equation in complex form with spatially varying coefficients. Since sinusoidal motion is assumed, the unsteady equation is independent of time. An out-of-core direct solution procedure was developed and applied to two-dimensional sections. Results are presented for a section of vanishing thickness in subsonic flow and an NACA 64A006 airfoil in supersonic flow. Good correlation is obtained in the first case at values of Mach number and reduced frequency of direct interest in flutter analyses. Reasonable results are obtained in the second case. Comparisons of two-dimensional finite difference solutions with exact analytic solutions indicate that the accuracy of the difference solution is dependent on the boundary conditions used on the outer boundaries. Homogeneous boundary conditions on the mesh edges that yield complex eigenvalues give the most accurate finite difference solutions. The plane outgoing wave boundary conditions meet these requirements.

Weatherill, W. H.↗
NASA NTRS1982

Improved Numerical Differencing Analyzer

SINDA, Systems Improved Numerical Differencing Analyzer, solves differential and algebric equations representing physical systems. SINDA solves numerically almost any set of ordinary differential equations that represent transient behavior of a lumped-parameter system or any set of nonlinear algebraic equations that represents the steady state conditions of a physical system.

Skladany, J. T.↗
NASA NTRS1988

Continuous Optimization on Constraint Manifolds

This paper demonstrates continuous optimization on the differentiable manifold formed by continuous constraint functions. The first order tensor geodesic differential equation is solved on the manifold in both numerical and closed analytic form for simple nonlinear programs. Advantages and disadvantages with respect to conventional optimization techniques are discussed.

Dean, Edwin B.↗