Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Parallels between control PDE's and systems of ODE's

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 differential 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 paralled 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 NTRS1992

Asymptotic integration algorithms for first-order ODEs with application to viscoplasticity

When constructing an algorithm for the numerical integration of a differential equation, one must first convert the known ordinary differential equation (ODE), which is defined at a point, into an ordinary difference equation (O(delta)E), which is defined over an interval. Asymptotic, generalized, midpoint, and trapezoidal, O(delta)E algorithms are derived for a nonlinear first order ODE written in the form of a linear ODE. The asymptotic forward (typically underdamped) and backward (typically overdamped) integrators bound these midpoint and trapezoidal integrators, which tend to cancel out unwanted numerical damping by averaging, in some sense, the forward and backward integrations. Viscoplasticity presents itself as a system of nonlinear, coupled first-ordered ODE's that are mathematically stiff, and therefore, difficult to numerically integrate. They are an excellent application for the asymptotic integrators. Considering a general viscoplastic structure, it is demonstrated that one can either integrate the viscoplastic stresses or their associated eigenstrains.

Freed, Alan D.↗
NASA NTRS1992

Modeling of the WSTF frictional heating apparatus in high pressure systems

In order to develop a computer program able to model the frictional heating of metals in high pressure oxygen or nitrogen a number of additions have been made to the frictional heating model originally developed for tests in low pressure helium. These additions include: (1) a physical property package for the gases to account for departures from the ideal gas state; (2) two methods for spatial discretization (finite differences with quadratic interpolation or orthogonal collocation on finite elements) which substantially reduce the computer time required to solve the transient heat balance; (3) more efficient programs for the integration of the ordinary differential equations resulting from the discretization of the partial differential equations; and (4) two methods for determining the best-fit parameters via minimization of the mean square error (either a direct search multivariable simplex method or a modified Levenburg-Marquardt algorithm). The resulting computer program has been shown to be accurate, efficient and robust for determining the heat flux or friction coefficient vs. time at the interface of the stationary and rotating samples.

Skowlund, Christopher T.↗
NASA NTRS1995

The Effect of Nonlinear Critical Layers on Boundary Layer Transition

Asymptotic methods are used to describe the nonlinear self-interaction between pairs of oblique instability modes that eventually develops when initially linear and spatially growing instability waves evolve downstream in nominally two-dimensional and spanwise periodic laminar boundary layers. The first nonlinear reaction takes place locally within a so-called 'critical layer' with the flow outside this layer consisting of a locally parallel mean flow plus an appropriate superposition of linear instability waves. The amplitudes of these waves are determined by either a single integro-differential equation or by a pair of integro-differential equations with quadratic to quartic-type nonlinearities.

Goldstein, Marvin E.↗
NASA NTRS1986

Nonlinear filtering for spacecraft attitude estimation

Nonlinear filtering techniques are applied to spacecraft attitude estimation using quaternion parameterization for the attitude kinematics. By replacing the angular velocity vector by the gyro output vector, a state dependent noise vector is introduced in the seven-dimensional system equations. The resulting conditional probability density function from the Ito differential rule is governed by the Fokker Planck partial differential equation which is approximated by the second order mean and covariance differential equations. In order to minimize computer loading, the covariance propagation is carried out in six-dimensional state space using a matrix transformation. The star tracker data is used to update the covariance matrix in the seven-dimensional space. The algorithm is simulated for an earth pointing spacecraft mission, using Monte Carlo samples of gyro and star measurements. The performance of the second order filter is compared with the extended Kalman Filter through several simulation runs and drift rates have been identified.

Vathsal, S.↗
NASA NTRS2008

Simulation of Stochastic Processes by Coupled ODE-PDE

A document discusses the emergence of randomness in solutions of coupled, fully deterministic ODE-PDE (ordinary differential equations-partial differential equations) due to failure of the Lipschitz condition as a new phenomenon. It is possible to exploit the special properties of ordinary differential equations (represented by an arbitrarily chosen, dynamical system) coupled with the corresponding Liouville equations (used to describe the evolution of initial uncertainties in terms of joint probability distribution) in order to simulate stochastic processes with the proscribed probability distributions. The important advantage of the proposed approach is that the simulation does not require a random-number generator.

Zak, Michail↗
NASA NTRS1973

Coincidence degree and periodic solutions of neutral equations

The problem of existence of periodic solutions for some nonautonomous neutral functional differential equations is examined. It is an application of a basic theorem on the Fredholm alternative for periodic solutions of some linear neutral equations and of a generalized Leray-Schauder theory. Although proofs are simple, the results are nontrivial extensions to the neutral case of existence theorems for periodic solutions of functional differential equations.

Hale, J. K.↗
NASA NTRS1973

Analysis and synthesis of distributed-lumped-active networks by digital computer

The use of digital computational techniques in the analysis and synthesis of DLA (distributed lumped active) networks is considered. This class of networks consists of three distinct types of elements, namely, distributed elements (modeled by partial differential equations), lumped elements (modeled by algebraic relations and ordinary differential equations), and active elements (modeled by algebraic relations). Such a characterization is applicable to a broad class of circuits, especially including those usually referred to as linear integrated circuits, since the fabrication techniques for such circuits readily produce elements which may be modeled as distributed, as well as the more conventional lumped and active ones.

Source record↗
NASA NTRS1976

Mathematical model of one-man air revitalization system

A mathematical model was developed for simulating the steady state performance in electrochemical CO2 concentrators which utilize (NMe4)2 CO3 (aq.) electrolyte. This electrolyte, which accommodates a wide range of air relative humidity, is most suitable for one-man air revitalization systems. The model is based on the solution of coupled nonlinear ordinary differential equations derived from mass transport and rate equations for the processes which take place in the cell. The boundary conditions are obtained by solving the mass and energy transport equations. A shooting method is used to solve the differential equations.

Source record↗
NASA NTRS1983

Parallel architectures for iterative methods on adaptive, block structured grids

A parallel computer architecture well suited to the solution of partial differential equations in complicated geometries is proposed. Algorithms for partial differential equations contain a great deal of parallelism. But this parallelism can be difficult to exploit, particularly on complex problems. One approach to extraction of this parallelism is the use of special purpose architectures tuned to a given problem class. The architecture proposed here is tuned to boundary value problems on complex domains. An adaptive elliptic algorithm which maps effectively onto the proposed architecture is considered in detail. Two levels of parallelism are exploited by the proposed architecture. First, by making use of the freedom one has in grid generation, one can construct grids which are locally regular, permitting a one to one mapping of grids to systolic style processor arrays, at least over small regions. All local parallelism can be extracted by this approach. Second, though there may be a regular global structure to the grids constructed, there will be parallelism at this level. One approach to finding and exploiting this parallelism is to use an architecture having a number of processor clusters connected by a switching network. The use of such a network creates a highly flexible architecture which automatically configures to the problem being solved.

Gannon, D.↗
NASA NTRS1984

Non-Gaussian approach for parametric random vibration of non-linear structures

The dynamic response of a nonlinear, single degree of freedom structural system subjected to a physically white noise parametric excitation is investigated. The Ito stochastic calculus is employed to derive a general differential equation for the moments of the response coordinates. The differential equations of moments of any order are found to be coupled with higher order moments. A non-Gaussian closure scheme is developed to truncate the moment equations up to fourth order. The statistical of the stationary response are computed numerically and compared with analytical solutions predicted by a Gaussian closure scheme and the stochastic averaging method. It is found that the computed results exhibit the jump phenomenon which is typical of the characteristics of deterministic nonlinear systems. In addition, the numerical algorithm leads to multiple solutions all of which give positive mean squares. However, two of these solutions are found to violate the properties of high order moments. One solution preserves the moments properties and demonstrates that the system achieves a stationary response.

Ibrahim, R. A.↗
NASA NTRS1996

An aeroelastic analysis of helicopter rotor blades incorporating piezoelectric fiber composite twist actuation

A simple aeroelastic analysis of a helicopter rotor blade incorporating embedded piezoelectric fiber composite, interdigitated electrode blade twist actuators is described. The analysis consist of a linear torsion and flapwise bending model coupled with a nonlinear ONERA based unsteady aerodynamics model. A modified Galerkin procedure is performed upon the rotor blade partial differential equations of motion to develop a system of ordinary differential equations suitable for numerical integration. The twist actuation responses for three conceptual full-scale blade designs with realistic constraints on blade mass are numerically evaluated using the analysis. Numerical results indicate that useful amplitudes of nonresonant elastic twist, on the order of one to two degrees, are achievable under one-g hovering flight conditions for interdigitated electrode poling configurations. Twist actuation for the interdigitated electrode blades is also compared with the twist actuation of a conventionally poled piezoelectric fiber composite blade. Elastic twist produced using the interdigitated electrode actuators was found to be four to five times larger than that obtained with the conventionally poled actuators.

Wilkie, W. Keats↗
NASA NTRS1996

Aeroelastic Analysis of Helicopter Rotor Blades Incorporating Anisotropic Piezoelectric Twist Actuation

A simple aeroelastic analysis of a helicopter rotor blade incorporating embedded piezoelectric fiber composite, interdigitated electrode blade twist actuators is described. The analysis consists of a linear torsion and flapwise bending model coupled with a nonlinear ONERA based unsteady aerodynamics model. A modified Galerkin procedure is performed upon the rotor blade partial differential equations of motion to develop a system of ordinary differential equations suitable for dynamics simulation using numerical integration. The twist actuation responses for three conceptual fullscale blade designs with realistic constraints on blade mass are numerically evaluated using the analysis. Numerical results indicate that useful amplitudes of nonresonant elastic twist, on the order of one to two degrees, are achievable under one-g hovering flight conditions for interdigitated electrode poling configurations. Twist actuation for the interdigitated electrode blades is also compared with the twist actuation of a conventionally poled piezoelectric fiber composite blade. Elastic twist produced using the interdigitated electrode actuators was found to be four to five times larger than that obtained with the conventionally poled actuators.

Wilkie, W. Keats↗
NASA NTRS1972

Smoothing properties of neutral equations

A neutral functional differential equation is defined as d/dt(Dx sub t) = f(x sub t), with D linear, continuous, atomic at zero. The solution generally is no smoother than the initial data after any finite number of steps. A more restrictive class of D-operators for which some smoothing takes place after an infinite number of steps is defined. This result indicates that a solution can be in an omega-limit set only if it corresponds to initial data which are smooth. A space which can be considered as a Banach space with the topology of uniform convergence is also defined.

Hale, J. K.↗
NASA NTRS1985

Chandrasekhar equations for infinite dimensional systems

Chandrasekhar equations are derived for linear time invariant systems defined on Hilbert spaces using a functional analytic technique. An important consequence of this is that the solution to the evolutional Riccati equation is strongly differentiable in time and one can define a strong solution of the Riccati differential equation. A detailed discussion on the linear quadratic optimal control problem for hereditary differential systems is also included.

Ito, K.↗
NASA NTRS1976

Stability of elastic bending and torsion of uniform cantilever rotor blades in hover with variable structural coupling

The stability of elastic flap bending, lead-lag bending, and torsion of uniform, untwisted, cantilever rotor blades without chordwise offsets between the elastic, mass, tension, and areodynamic center axes is investigated for the hovering flight condition. The equations of motion are obtained by simplifying the general, nonlinear, partial differential equations of motion of an elastic rotating cantilever blade. The equations are adapted for a linearized stability analysis in the hovering flight condition by prescribing aerodynamic forces, applying Galerkin's method, and linearizing the resulting ordinary differential equations about the equilibrium operating condition. The aerodynamic forces are obtained from strip theory based on a quasi-steady approximation of two-dimensional unsteady airfoil theory. Six coupled mode shapes, calculated from free vibration about the equilibrium operating condition, are used in the linearized stability analysis. The study emphasizes the effects of two types of structural coupling that strongly influence the stability of hingeless rotor blades. The first structural coupling is the linear coupling between flap and lead-lag bending of the rotor blade. The second structural coupling is a nonlinear coupling between flap bending, lead-lag bending, and torsion deflections. Results are obtained for a wide variety of hingeless rotor configurations and operating conditions in order to provide a reasonably complete picture of hingeless rotor blade stability characteristics.

Hodges, D. H., Roberta.↗
NASA NTRS1972

Implications of a quadratic stream definition in radiative transfer theory.

An explicit definition of the radiation-stream concept is stated and applied to approximate the integro-differential equation of radiative transfer with a set of twelve coupled differential equations. Computational efficiency is enhanced by distributing the corresponding streams in three-dimensional space in a totally symmetric way. Polarization is then incorporated in this model. A computer program based on the model is briefly compared with a Monte Carlo program for simulation of horizon scans of the earth's atmosphere. It is found to be considerably faster.

Whitney, C.↗