Search NASA⌕ Search

SEARCH · Search NASA

Results for “DIFFERENCE EQUATION”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

A guidance and navigation system for continuous low-thrust vehicles

A midcourse guidance and navigation system for continuous low thrust vehicles was developed. The equinoctial elements are the state variables. Uncertainties are modelled statistically by random vector and stochastic processes. The motion of the vehicle and the measurements are described by nonlinear stochastic differential and difference equations respectively. A minimum time trajectory is defined; equations of motion and measurements are linearized about this trajectory. An exponential cost criterion is constructed and a linear feedback quidance law is derived. An extended Kalman filter is used for state estimation. A short mission using this system is simulated. It is indicated that this system is efficient for short missions, but longer missions require accurate trajectory and ground based measurements.

Jack-Chingtse, C.↗
NASA NTRS1975

Generalized wave envelope analysis of sound propagation in ducts with stepped noise source profiles and variable axial impedance

A finite difference formulation is presented for sound propagation in a rectangular two-dimensional duct without steady flow. Before the difference equations are formulated, the governing Helmholtz equation is first transformed to a form whose solution tends not to oscillate along the length of the duct. This transformation reduces the required number of grid points by an order of magnitude. Example solutions indicate that stepped noise source profiles have much higher attenuation than plane waves in a uniform impedance liner. Also, multiple stepped impedance liners are shown to have higher attenuation than uniform ducts if the impedances are chosen properly. For optimum noise reduction with axial variations in impedance, the numerical analysis indicates that for a plane wave input the resistance should be near zero at the entrance of a suppressor duct, while the reactance should be near the optimum value associated with the least-attenuated mode in a uniform duct.

Baumeister, K. J.↗
NASA NTRS1974

Subsonic and supersonic indicial aerodynamics and aerodynamic transfer function for complex configurations

A general theory for indicial-potential-compressible aerodynamics around complex configurations is presented. The motion is assumed to consist of constant subsonic or supersonic speed (steady state) and small perturbations around the steady state. Using the finite-element method to discretize the space problem, a set of differential-difference equations in time relating the potential to its normal derivative on the surface of the body was obtained. The aerodynamics transfer function was derived by using standard method of operational calculus.

Morino, L.↗
NASA NTRS1975

Refinement and validation of two digital Microwave Landing System (MLS) theoretical models

Two digital microwave landing system theoretical models are considered which are generic models for the Doppler and scanning-beam frequency reference versions of the MLS. These models represent errors resulting from both system noise and discrete multipath. The data used for the validation effort were obtained from the Texas Instrument conventional scanning beam and the Hazeltine Doppler feasibility hardware versions of the MLS. Topics discussed include tape read software, time history plots, computation of power spectral density, smoothed power spectra, best-fit models, different equations for digital simulation, and discrete multipath errors.

Duff, W. G.↗
NASA NTRS1975

Higher order accurate partial implicitization: An unconditionally stable fourth-order-accurate explicit numerical technique

The previously obtained second-order-accurate partial implicitization numerical technique used in the solution of fluid dynamic problems was modified with little complication to achieve fourth-order accuracy. The Von Neumann stability analysis demonstrated the unconditional linear stability of the technique. The order of the truncation error was deduced from the Taylor series expansions of the linearized difference equations and was verified by numerical solutions to Burger's equation. For comparison, results were also obtained for Burger's equation using a second-order-accurate partial-implicitization scheme, as well as the fourth-order scheme of Kreiss.

Graves, R. A., Jr.↗
NASA NTRS1975

On the convergence of certain finite-difference schemes by an inverse-matrix method

The inverse-matrix method of analyzing the convergence of the solution of a given system of finite-difference equations to the solution of the corresponding system of partial-differential equations is discussed and generalized. The convergence properties of a time- and space-centered differencing of the diffusion equation are analyzed as well as a staggered grid differencing of the Cauchy-Riemann equations. These two schemes are significant since they serve as simplified model algorithms for two recently developed methods used to calculate nonlinear aerodynamic flows.

Steger, J. L.↗
NASA NTRS1975

Relaxation methods in fluid mechanics

The present work considers the iterative solution of a coupled set of difference equations and examines methods that carry successive approximates to a state that is invariant with further iteration and independent of the initial guess. Methods are studied with regard to their efficiency and economy of computer resources. The basic principles of classical relaxation are set forth, with attention confined to linear elliptic equations. This discussion involves the evaluation of the spectral radius that is the magnitude of the eigenvalue with largest modulus. The subject of relaxation is then related to the study of ordinary differential equations and hyperbolic partial differential equations. Problems that occur when linearly dependent eigenvectors appear in the relaxation matrix are discussed, leading to multiply connected eigenvalues in the Jordan canonical form. Finally, a brief survey of relaxation methods used in aerodynamics is given.

Lomax, H.↗
NASA NTRS1975

Analytical method for determining the stability of linear retarded systems with two delays

The stability is considered of the solution differential-difference equations of the retarded type with constant coefficients and two constant time delays. A method that makes use of analytical expressions to determine stability boundaries, and the stability of the system, is derived. The method was applied to a system represented by a second-order differential equation with constant coefficients and time delays in the velocity and displacement terms. The results obtained is in agreement with those obtained by other investigators.

Barker, L. K.↗
NASA NTRS1975

A finite-difference analysis of the nozzle starting process in an expansion tunnel

A suitable finite-difference method for computing the quasi-one-dimensional unsteady flow in an expansion tunnel nozzle was identified. The difference equations are presented along with the appropriate stability limits. A parametric study of the starting process in an expansion tunnel nozzle is made, and acceptable operating conditions were determined.

Weilmuenster, K. J.↗
NASA NTRS1975

Kalman filter estimation of human pilot-model parameters

The parameters of a human pilot-model transfer function are estimated by applying the extended Kalman filter to the corresponding retarded differential-difference equations in the time domain. Use of computer-generated data indicates that most of the parameters, including the implicit time delay, may be reasonably estimated in this way. When applied to two sets of experimental data obtained from a closed-loop tracking task performed by a human, the Kalman filter generated diverging residuals for one of the measurement types, apparently because of model assumption errors. Application of a modified adaptive technique was found to overcome the divergence and to produce reasonable estimates of most of the parameters.

Schiess, J. R.↗
NASA NTRS1975

Cost performance satellite design using queueing theory

A modified Poisson arrival, infinite server queuing model is used to determine the effects of limiting the number of broadcast channels (C) of a direct broadcast satellite used for public service purposes (remote health care, education, etc.). The model is based on the reproductive property of the Poisson distribution. A difference equation has been developed to describe the change in the Poisson parameter. When all initially delayed arrivals reenter the system a (C plus 1) order polynomial must be solved to determine the effective value of the Poisson parameter. When less than 100% of the arrivals reenter the system the effective value must be determined by solving a transcendental equation. The model was used to determine the minimum number of channels required for a disaster warning satellite without degradation in performance. Results predicted by the queuing model were compared with the results of digital simulation.

Hein, G. F.↗
NASA NTRS1976

Transition from incoherence to coherence in the spontaneous emission of extended systems

Using the results of a previous paper, we study the evolution of a completely excited system of N two-level atoms, distributed in a circular cylinder with dimensions much larger than the wavelength, interacting with all modes of radiation field. We emphasize the distinction between r-conserving (rC) and r-nonconserving (rNC) processes (r is the so-called cooperation number). Considering the number of photons emitted as the discrete independent variable, the evolution is described by a partial difference equation. Numerical solutions of this equation show the transition from rNC dominance at the beginning to rC dominance later. This is also a transition from incoherent to coherent emission of radiation.

Lee, C. T.↗
NASA NTRS1977

Fast methods incorporating direct elliptic solvers for nonlinear applications in fluid dynamics

Semidirect methods are discussed, their present role, as well as some developments for their application in computational fluid dynamics. A semidirect method is a computational scheme that uses a fast, direct, elliptic solver as the driving algorithm for the iterative solution of finite difference equations. Specific subtopics include: (1) direct Cauchy Riemann solvers for first order elliptic equations; (2) application of the semidirect method to the mixed elliptic hyperbolic problem of steady, inviscid transonic flow; and (3) the treatment of interior conditions, such as those on an airfoil or wing, in semidirect methods.

Martin, E. D.↗
NASA NTRS1979

Algorithms for linear-systems control

Control-theory design package aids design of linear-quadratic-Gaussian (LOG) controllers and optimal filters. It is applicable to systems that can be modeled by linear time-invariant differential or difference equations.

Armstrong, E. S.↗
NASA NTRS1977

Piping design considerations in a solar-Rankine power plant

Two of the main parameters in sizing the piping of a solar power plant are the working pressure of the vapor leaving the solar collectors, and the type of working fluid used. Numerical examples for each case are given using the graphical Moody friction charts and the analytical Darcy-Weisbach equation. Different working pressures of steam vapor in the solar collector-turbine pipe connection indicate their major role in the design. The size variation was found not to be in linear proportion to vapor density variations. On the other hand, high molecular weight organic fluids such as R-11 and R-113, when compared with water, show insignificant changes in piping sizes.

Lansing, F. L.↗
NASA NTRS1977

Impact on multilayered composite plates

Stress wave propagation in a multilayer composite plate due to impact was examined by means of the anisotropic elasticity theory. The plate was modelled as a number of identical anisotropic layers and the approximate plate theory of Mindlin was then applied to each layer to obtain a set of difference-differential equations of motion. Dispersion relations for harmonic waves and correction factors were found. The governing equations were reduced to difference equations via integral transforms. With given impact boundary conditions these equations were solved for an arbitrary number of layers in the plate and the transient propagation of waves was calculated by means of a Fast Fourier Transform algorithm. The multilayered plate problem was extended to examine the effect of damping layers present between two elastic layers. A reduction of the interlaminar normal stress was significant when the thickness of damping layer was increased but the effect was mostly due to the softness of the damping layer. Finally, the problem of a composite plate with a crack on the interlaminar boundary was formulated.

Kim, B. S.↗
NASA NTRS1978

ORACLS: A system for linear-quadratic-Gaussian control law design

A modern control theory design package (ORACLS) for constructing controllers and optimal filters for systems modeled by linear time-invariant differential or difference equations is described. Numerical linear-algebra procedures are used to implement the linear-quadratic-Gaussian (LQG) methodology of modern control theory. Algorithms are included for computing eigensystems of real matrices, the relative stability of a matrix, factored forms for nonnegative definite matrices, the solutions and least squares approximations to the solutions of certain linear matrix algebraic equations, the controllability properties of a linear time-invariant system, and the steady state covariance matrix of an open-loop stable system forced by white noise. Subroutines are provided for solving both the continuous and discrete optimal linear regulator problems with noise free measurements and the sampled-data optimal linear regulator problem. For measurement noise, duality theory and the optimal regulator algorithms are used to solve the continuous and discrete Kalman-Bucy filter problems. Subroutines are also included which give control laws causing the output of a system to track the output of a prescribed model.

Armstrong, E. S.↗
NASA NTRS1978

On reliable control system designs

A mathematical model for use in the design of reliable multivariable control systems is discussed with special emphasis on actuator failures and necessary actuator redundancy levels. The model consists of a linear time invariant discrete time dynamical system. Configuration changes in the system dynamics are governed by a Markov chain that includes transition probabilities from one configuration state to another. The performance index is a standard quadratic cost functional, over an infinite time interval. The actual system configuration can be deduced with a one step delay. The calculation of the optimal control law requires the solution of a set of highly coupled Riccati-like matrix difference equations. Results can be used for off-line studies relating the open loop dynamics, required performance, actuator mean time to failure, and functional or identical actuator redundancy, with and without feedback gain reconfiguration strategies.

Birdwell, J. D.↗