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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Atmospheric model development in support of SEASAT. Volume 3: Prediction models

As part of the SEASAT program of NASA, a set of four hemispheric, atmospheric models are developed. The four models use a polar stereographic grid in the horizontal and a sigma coordinate in the vertical. Conservation forms of the difference equations are integrated using either a fifteen or four minute time step on a 381 km or 127 km grid (at 60 deg N) for the 63 x 63 or 187 x 187 models respectively. A nonlinear pressure smoother is used along with momentum and temperature diffusion to help control computational noise. The horizontal boundary conditions are insulated slippery walls. Centered time differencing with time averaging of the pressure gradient force term in the momentum equations is used. Robert time filtering on of the temperature and moisture solutions is used for computational stability. The moisture and heat source/sink terms are modeled. Terms representing evaporation and large scale condensation, sensible heat exchange, parameterized cumulus convection and precipitation, and solar and terrestrial radiation are included. Initialization of the models is based on a pattern conservation technique to obtain objective analysis of the state parameter structure from the surface to 50 mb.

Wellck, R. E.↗

NASA NTRS1978

A split-recoupled-semidirect computational technique applied to transonic flow over lifting airfoils

A new version of the semidirect iterative method eliminates significant restrictions of previous versions of the method. A semidirect method solves finite-difference equations by a rapid globally implicit iterative process driven by a fast direct elliptic solver. The new approach can treat complex systems of equations in an efficient 'correction form', and allows the use of general, nonorthogonal, boundary-fitted coordinate transformations. These features are expected to lead to significant practical applications with conservation-equation systems in either two or three dimensions. The present application to the full potential equations for steady transonic flow over an airfoil at angle of attack illustrates the utility of the technique.

Martin, E. D.↗

NASA NTRS1978

Steady flow of a non-Newtonian fluid through a contraction

A steady-state analysis is conducted to examine the basic flow structure of a non-Newtonian fluid in a domain including an inflow region, a contraction region, and an outflow region. A Cartesian grid system is used throughout the entire flow domain, including the contraction region, thus creating an irregular grid cell structure adjacent to the curved boundary. At node points adjacent to the curved boundary symmetry conditions are derived for the different flow variables in order to solve the governing difference equations. Attention is given to the motion and non-Newtonian constitutive equations, the boundary conditions, the numerical modeling of the non-Newtonian equations, the stream function contour lines for the non-Newtonian fluid, the vorticity contour lines for the non-Newtonian fluid, the velocity profile across the contraction, and the shear stress contour lines for the non-Newtonian fluid.

Gatski, T. B.↗

NASA NTRS1978

The development of a three-dimensional partially elliptic flow computer program for combustor research

A three dimensional, partially elliptic, computer program was developed. Without requiring three dimensional computer storage locations for all flow variables, the partially elliptic program is capable of predicting three dimensional combustor flow fields with large downstream effects. The program requires only slight increase of computer storage over the parabolic flow program from which it was developed. A finite difference formulation for a three dimensional, fully elliptic, turbulent, reacting, flow field was derived. Because of the negligible diffusion effects in the main flow direction in a supersonic combustor, the set of finite-difference equations can be reduced to a partially elliptic form. Only the pressure field was governed by an elliptic equation and requires three dimensional storage; all other dependent variables are governed by parabolic equations. A numerical procedure which combines a marching integration scheme with an iterative scheme for solving the elliptic pressure was adopted.

Pan, Y. S.↗

NASA NTRS1979

Numerical modeling of two-dimensional confined flows

A numerical model of two-dimensional confined flows is presented. The flow in the duct is partitioned into finite streams. The difference equations are then obtained by applying conservation principles directly to the individual streams. A listing of a computer code based on this approach in FORTRAN 4 language is presented. The code computes two dimensional compressible turbulent flows in ducts when the duct area along the flow is specified and the pressure gradient is unknown.

Greywall, M. S.↗

NASA NTRS1979

Robustness of multiloop linear feedback systems

This paper presents a new approach to the frequency-domain analysis of multiloop linear feedback systems. The properties of the return difference equation are examined using the concepts of singular values, singular vectors and the spectral norm of a matrix. A number of new tools for multiloop systems are developed which are analogous to those for scalar Nyquist and Bode analysis. These provide a generalization of the scalar frequency-domain notions such as gain, bandwidth, stability margins and M-circles, and provide considerable insight into system robustness.

Doyle, J. C.↗

NASA NTRS1979

Time-dependent difference theory for noise propagation in a two-dimensional duct

A time dependent numerical formulation was derived for sound propagation in a two dimensional straight soft-walled duct in the absence of mean flow. The time dependent governing acoustic-difference equations and boundary conditions were developed along with the maximum stable time increment. Example calculations were presented for sound attenuation in hard and soft wall ducts. The time dependent analysis were found to be superior to the conventional steady numerical analysis because of much shorter solution times and the elimination of matrix storage requirements.

Baumeister, K. J.↗

NASA NTRS1979

Accuracy and performance of LACIE crop development models

Of the three principal phenological crop calendar models evaluated for LACIE, Robertson's triquadratic model which predicts the rate of progression of wheat through its biological development, was selected. Daily maximum and minimum temperatures and day length are the input variables, and the principal output is a daily increment of development through six physiological growth stages. Because wheat corresponds differently to the environment during each growth stage, five different equations are required. The estimated and observed crop development data were compared in order to establish a measure of confidence in the model and to identify consistent discrepancies that would adversely affect LACIE operation. Although the model provided reliable estimates for various wheat growing regions of the world, it was found that there are still areas in need of further model improvement or development.

Woolley, S. K.↗

NASA NTRS1978

Multi-Level Adaptive Techniques (MLAT) for singular-perturbation problems

The multilevel (multigrid) adaptive technique, a general strategy of solving continuous problems by cycling between coarser and finer levels of discretization is described. It provides very fast general solvers, together with adaptive, nearly optimal discretization schemes. In the process, boundary layers are automatically either resolved or skipped, depending on a control function which expresses the computational goal. The global error decreases exponentially as a function of the overall computational work, in a uniform rate independent of the magnitude of the singular-perturbation terms. The key is high-order uniformly stable difference equations, and uniformly smoothing relaxation schemes.

Brandt, A.↗

NASA NTRS1980

Time-dependent difference theory for noise propagation in a two-dimensional duct

A time-dependent numerical formulation is derived for sound propagation in a two-dimensional straight soft-walled duct in the absence of mean flow. The time-dependent governing acoustic-difference equations and boundary conditions are developed along with the maximum stable time increment. Example calculations are presented for sound attenuation in hard- and soft-wall ducts. The time-dependent analysis has been found to be superior to the conventional steady numerical analysis because of much shorter solution times and the elimination of matrix storage requirements.

Baumeister, K. J.↗

NASA NTRS1980

Streamwise computation of duct flows

A computational method is presented to calculate momentum and energy transport in two-dimensional viscous compressible duct flows. The flow in the duct is partitioned into finite streams. The difference equations are then obtained by applying momentum and energy conservation principles directly to the individual streams. The method is applicable to laminar and turbulence flows.

Greywall, M. S.↗

NASA NTRS1980

An adaptive computation mesh for the solution of singular perturbation problems

In singular perturbation problems, control of zone size variation can affect the effort required to obtain accurate, numerical solutions of finite difference equations. The mesh is generated by the solution of potential equations. Numerical results for a singular perturbation problem in two dimensions are presented. The mesh was used in calculations of resistive magnetohydrodynamic flow in two dimensions.

Brackbill, J. U.↗

NASA NTRS1981

Influence of exit impedance on finite difference solutions of transient acoustic mode propagation in ducts

The time-dependent governing acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard-wall duct without flow and with spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the steady solutions of the acoustic field. The transient method did not converge to the steady-state solution for cutoff acoustic duct modes. This has implications as to its use in a variable-area duct, where modes may become cutoff in the smal-area portion of the duct. For single cutoff mode propagation the steady-state impedance boundary condition produced acoustic reflections during the initial transient that caused finite instabilities in the numerical calculations. The stability problem is resolved by reformulating the exit boundary condition. Example calculations show good agreement with exact analytical and numerical results for forcing frequencies above, below, and nearly at the cutoff frequency.

Baumeister, K. J.↗

NASA NTRS1981

Basic advances in the finite-volume method for transonic potential flow calculations

The finite-volume method of Jameson and Caughey provides a framework within which it is possible to calculate transonic potential flows about essentially arbitrary geometrical configurations. Improvements designed to increase the accuracy of the basic scheme and its consistency in the far field will be described. These include the incorporation of an artificial viscosity which maintains the formal second-order accuracy of the scheme in supersonic zones, and a modification of the flux balances to allow the free-stream conditions to satisfy the difference equations identically. Results of calculations illustrating the importance of these effects will be presented.

Caughey, D. A.↗

NASA NTRS1981

Galactic civilizations - Population dynamics and interstellar diffusion

A model is developed of the interstellar diffusion of galactic civilizations which takes into account the population dynamics of such civilizations. The problem is formulated in terms of potential theory, with a family of nonlinear partial differential and difference equations specifying population growth and diffusion for an organism with advantageous genes that undergoes random dispersal while increasing in population locally, and a population at zero population growth. In the case of nonlinear diffusion with growth and saturation, it is found that the colonization wavefront from the nearest independently arisen galactic civilization can have reached the earth only if its lifetime exceeds 2.6 million years, or 20 million years if discretization can be neglected. For zero population growth, the corresponding lifetime is 13 billion years. It is concluded that the earth is uncolonized not because interstellar spacefaring civilizations are rare, but because there are too many worlds to be colonized in the plausible colonization lifetime of nearby civilizations, and that there exist no very old galactic civilizations with a consistent policy of the conquest of inhabited worlds.

Newman, W. I.↗

NASA NTRS1982

Optimal Regulator Algorithms for the Control of Linear Systems

Control-theory design package, called Optimal Regulator Algorithms for Control of Linear Systems (ORACLS), aids in design of controllers and optimal filters for systems modeled by linear, time invariant differential and difference equations. ORACLS is particularly attractive rigorous tool for dealing with multi-input and multi-output dynamic systems in both continuous and discrete forms.

Armstrong, E. S.↗

NASA NTRS1981

Electrochemical models for the discharge characteristics of the nickel cadmium cell

The potential time characteristics of a preconditioned fully charged cell discharge at constant current was studied. Electrochemical principles applied to the sealed nickel cadmium cell and its behavior and to predict operating characteristics were described. A thermodynamic approach to arrive at several related but different equations and its discharge are reported.

Spritzer, M. S.↗

NASA NTRS1981

Numerical Modeling of Three-Dimensional Confined Flows

A three dimensional confined flow model is presented. The flow field is computed by calculating velocity and enthalpy along a set of streamlines. The finite difference equations are obtained by applying conservation principles to streamtubes constructed around the chosen streamlines. With appropriate substitutions for the body force terms, the approach computes three dimensional magnetohydrodynamic channel flows. A listing of a computer code, based on this approach is presented in FORTRAN IV language. The code computes three dimensional compressible viscous flow through a rectangular duct, with the duct cross section specified along the axis.

Greywall, M. S.↗