Search NASA⌕ Search

SEARCH · Search NASA

Results for “Functional 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 109 records · Page 6

A study of the efficiency of various Navier-Stokes solvers

A comparative study of the efficiency of some finite difference methods for the solution of the Navier-Stokes equations was conducted. The study was restricted to the two-dimensional steady, uniform property vorticity-stream function equations. The comparisons were drawn by recording the CPU time required to obtain a solution as well as the accuracy of this solution using five numerical methods: central differences, first order upwind differences, second order upwind differences, exponential differences, and an ADI solution of the central difference equations. Solutions were obtained for two test cases: a recirculating eddy inside a square cavity with a moving top, and an impinging jet flow. The results show that whenever the central difference method is stable it generates results with a given accuracy for less CPU time than any other method.

Atias, M.↗

Marshall system for aerospace system simulation (MARSYAS), user's manual

The capabilities of the Marshall system for aerospace system simulation (MARSYAS) and how to use it are described. MARSYAS is a software system that allows easy setup and control of the simulation of the dynamics of large physical systems on a digital computer. The physical systems are modeled in the form of block diagrams or equations. The blocks can have multiple inputs and multiple outputs, and they can be nested to form hierarchies. The block diagrams can contain transfer functions, nonlinear and logical functions, equations, analog computer elements and FORTRAN programs. The input format of the equations can be combinations of nonlinear, time-varying differential equations and algebraic equations in their original format. MARSYAS could also serve as a storage and retrieval system for models as a basis for a model configuration control system on a central time-shared computer. The outputs of the simulation system can be not only time-responses but also other analysis data such as frequency response, power spectrum and stability parameters. The MARSYAS translator is written in FORTRAN running on the Univac 1108 computer under the EXEC 8 operating system.

Ventre, A.↗

Rapid finite-difference computation of subsonic and transonic aerodynamic flows

Rapid iterative (or semidirect) computation methods are developed for the finite-difference solution of the nonlinear equations of subsonic and transonic aerodynamics. At each iteration, a fast, direct elliptic algorithm solves the entire computation field. In an application to subsonic flow over a lifting airfoil, the full nonlinear stream-function equation is solved. Finally, a direct Cauchy-Riemann solver is used for the nonlinear transonic small-disturbance equations for a biconvex airfoil. At M = 0.7, t/c = 0.1 (subcritical), three iterations on a 39 x 32 mesh (totaling 2.45 sec on an IBM 360/67 computer) obtain convergence within 0.1%. A slightly supercritical case requires seven iterations (6.75 sec) for convergence within 1%.

Martin, E. D.↗

Unsteady Solution of Non-Linear Differential Equations Using Walsh Function Series

Walsh functions form an orthonormal basis set consisting of square waves. The discontinuous nature of square waves make the system well suited for representing functions with discontinuities. The product of any two Walsh functions is another Walsh function - a feature that can radically change an algorithm for solving non-linear partial differential equations (PDEs). The solution algorithm of non-linear differential equations using Walsh function series is unique in that integrals and derivatives may be computed using simple matrix multiplication of series representations of functions. Solutions to PDEs are derived as functions of wave component amplitude. Three sample problems are presented to illustrate the Walsh function series approach to solving unsteady PDEs. These include an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the use of the Walsh function solution algorithms, exploiting Fast Walsh Transforms in multi-dimensions (O(Nlog(N))). Details of a Fast Walsh Reciprocal, defined here for the first time, enable inversion of aWalsh Symmetric Matrix in O(Nlog(N)) operations. Walsh functions have been derived using a fractal recursion algorithm and these fractal patterns are observed in the progression of pairs of wave number amplitudes in the solutions. These patterns are most easily observed in a remapping defined as a fractal fingerprint (FFP). A prolongation of existing solutions to the next highest order exploits these patterns. The algorithms presented here are considered a work in progress that provide new alternatives and new insights into the solution of non-linear PDEs.

Gnoffo, Peter A.↗

Attitude stability of spinning flexible spacecraft

The stability of spinning flexible satellites in a force-free environment was analyzed. The satellite was modeled as a rigid core having attached to it a flexible appendage idealized as a collection of particles (point masses) interconnected by springs. Both Liapunov and Routh-Hurwitz stability procedures are used. In the former, the Hamiltonian of the system, constrained through the angular momentum integral so as to admit complete damping, is used as a testing function. Equations of motion are written using the hybrid coordinate formulation, which readily accepts a modal coordinate transformation ultimately allowing truncation to a level amenable to literal stability analysis. Closed form stability criteria are generated for the first mode of a restricted appendage model lying in a plane containing the system center of mass and orthogonal to the spin axis. The effects of spin on flexible bodies are discussed by considering a very elementary particle model. Control of passively unstable spacecraft is briefly considered.

Likins, P. W.↗

Momentum constraints as integrability conditions for the Hamiltonian constraint in general relativity.

It is shown that if the Hamiltonian constraint of general relativity is imposed as a restriction on the Hamilton principal functional in the classical theory, or on the state functional in the quantum theory, then the momentum constraints are automatically satisfied. This result holds both for closed and open spaces and it means that the full content of the theory is summarized by a single functional equation of the Tomonaga-Schwinger type.

Moncrief, V.↗

Calculations of rotational flows using stream function

The stream function equation is solved for steady two-dimensional (and axisymmetric) rotational flows. Both finite differences and finite volumes discretization techniques are studied, using generalized body fitted coordinates and unstructured staggered grids, respectively. For inviscid transonic flows, a new artificial viscosity scheme which does not produce any artificial vorticity is introduced, for the stability of the mixed flow calculations and for capturing shocks. The solution of Euler equations, in primitive variables, are also considered. The effects of the artificial viscosity and numerical boundary conditions on the total enthalpy and the vorticity distributions are demonstrated.

Hafez, M.↗

Hopf bifurcation in the driven cavity

The algorithm employed in the present incompressible two-dimensional calculations of an impulsively-started lid-driven cavity has its basis in the time-dependent stream-function equation. While a Crank-Nicholson differencing scheme is used for the diffusion terms, the Adams-Bashforth scheme is used for the convection terms. The periodic asymptotic solutions obtained for Reynolds numbers of 5000 and 10,000 are found to be precisely periodic; it is demonstrated that they have reached asymptotic states. The indicators of that achievement are discussed.

Goodrich, John W.↗

A mathematical formulation for the control of the SCOLE configuration

A mathematical formulation for the control of the Spacecraft Control Laboratory Experiment (SCOLE) configuration is given. Two equivalent approaches, one using a functional equation and the other an abstract wave equation, are illustrated. Such a formulation can help in digital computer simulation to evaluate control laws, providing insight, and generating control laws.

Balakrishnan, A. V.↗

Estimating Bias Error Distributions

This paper formulates the general methodology for estimating the bias error distribution of a device in a measuring domain from less accurate measurements when a minimal number of standard values (typically two values) are available. A new perspective is that the bias error distribution can be found as a solution of an intrinsic functional equation in a domain. Based on this theory, the scaling- and translation-based methods for determining the bias error distribution arc developed. These methods are virtually applicable to any device as long as the bias error distribution of the device can be sufficiently described by a power series (a polynomial) or a Fourier series in a domain. These methods have been validated through computational simulations and laboratory calibration experiments for a number of different devices.

Liu, Tian-Shu↗

Observations and Parameterizations of Particle Size Distributions in Deep Tropical Cirrus and Stratiform Precipitation Clouds: Results from In-Situ Observations in TRMM Field Campaigns

In this study, we report on the evolution of particle size distributions (PSDs) and habits as measured during slow, Lagrangian-type spiral descents through deep subtropical and tropical cloud layers in Florida, Brazil, and Kwajalein, Marshall Islands, most of which were precipitating. The objective of the flight patterns was to learn more about how the PSDs evolved in the vertical and to obtain information of the vertical structure of microphysical properties. New instrumentation yielding better information on the concentrations of particles in the size (D) range between 0.2 and 2 cm, as well as improved particle imagery, produced more comprehensive observations for tropical stratiform precipitation regions and anvils than have been available previously. Collocated radar observations provided additional information on the vertical structure of the cloud layers sampled. Most of the spirals began at cloud top, with temperatures (T) as low as -50 C, and ended at cloud base or below the melting layer (ML). The PSDs broadened from cloud top towards cloud base, with the largest particles increasing in size from several millimeters at cloud top to one centimeter or larger towards cloud base. Some continued growth was noted in the upper part of the ML. Concentrations of particles less than 1 mm in size decreased with decreasing height. The result was a consistent change in the PSDs in the vertical. Similarly, systematic changes in the size dependence of the particle cross-sectional area was noted with decreasing height. Aggregation-as ascertained from both the changes in the PSDs and evolution of particle habits as observed in high detail with the cloud particle imager (CPI) probe-was responsible for these trends. The PSDs were generally well-represented by gamma distributions of the form N = N0 gamma D microns e- lambda gamma D that were fitted to the PSDs over 1-km horizontal intervals throughout the spirals. The intercept (N0 gamma), slope (lambda gamma), and dispersion (microns) values were derived for each PSD. Exponential curves (N = N0e- lambdaD; micron = 0) were also fitted to the distributions. The lambda gamma values for given spirals varied systematically with temperature as did the values of lambda (exponential), and the data generally conformed to values found in previous studies involving exponential fits to size distributions in mid-latitude frontal and cirrus layers. Considerable variability often noted in the PSD properties during the loops of individual spirals was manifested primarily in large changes in N0 gamma and N0, but micron, lambda gamma and lambda remained fairly stable. Temperature is not found to be the sole factor controlling lambda gamma or lambda but is a primary one. Direct relationships were found between lambda gamma and N0 gamma or lambda gamma and micron for the gamma distributions and lambda and N0 for the exponential. The latter relationship was not found as distinctly in earlier studies; observed PSDs in this study had better fidelity with less scatter. The micron values changed monotonically with T over the range of temperatures and were directly related to N0 gamma or lambda gamma, thereby reducing the number of variables in the PSD functional equation to two. In the upper part of the ML, N0 and lambda continued to decrease, and in the lower part these values began to increase as the largest particles melted. We developed general expressions relating various bulk microphysical, radar, and radiative transfer-related variables to N0 gamma and lambda gamma, useful for both tropical and mid-latitude clouds. These relationships facilitate the specification of a number of bulk properties in cloud and climate models. The results presented in this paper apply best to temperatures between 0 and -40 C, for which the measured radar reflectivities fall in the range of 0 to 25 dBZe.

Heymsfield, Andrew J.↗

Practical Loop-Shaping Design of Feedback Control Systems

An improved methodology for designing feedback control systems has been developed based on systematically shaping the loop gain of the system to meet performance requirements such as stability margins, disturbance attenuation, and transient response, while taking into account the actuation system limitations such as actuation rates and range. Loop-shaping for controls design is not new, but past techniques do not directly address how to systematically design the controller to maximize its performance. As a result, classical feedback control systems are designed predominantly using ad hoc control design approaches such as proportional integral derivative (PID), normally satisfied when a workable solution is achieved, without a good understanding of how to maximize the effectiveness of the control design in terms of competing performance requirements, in relation to the limitations of the plant design. The conception of this improved methodology was motivated by challenges in designing control systems of the types needed for supersonic propulsion. But the methodology is generally applicable to any classical control-system design where the transfer function of the plant is known or can be evaluated. In the case of a supersonic aerospace vehicle, a major challenge is to design the system to attenuate anticipated external and internal disturbances, using such actuators as fuel injectors and valves, bypass doors, and ramps, all of which are subject to limitations in actuator response, rates, and ranges. Also, for supersonic vehicles, with long slim type of structures, coupling between the engine and the structural dynamics can produce undesirable effects that could adversely affect vehicle stability and ride quality. In order to design distributed controls that can suppress these potential adverse effects, within the full capabilities of the actuation system, it is important to employ a systematic control design methodology such as this that can maximize the effectiveness of the control design in a methodical and quantifiable way. The emphasis is in generating simple but rather powerful design techniques that will allow even designers with a layman s knowledge in controls to develop effective feedback control designs. Unlike conventional ad hoc methodologies of feedback control design, in this approach actuator rates are incorporated into the design right from the start: The relation between actuator speeds and the desired control bandwidth of the system is established explicitly. The technique developed is demonstrated via design examples in a step-by-step tutorial way. Given the actuation system rates and range limits together with design specifications in terms of stability margins, disturbance rejection, and transient response, the procedure involves designing the feedback loop gain to meet the requirements and maximizing the control system effectiveness, without exceeding the actuation system limits and saturating the controller. Then knowing the plant transfer function, the procedure involves designing the controller so that the controller transfer function together with the plant transfer function equate to the designed loop gain. The technique also shows what the limitations of the controller design are and how to trade competing design requirements such as stability margins and disturbance rejection. Finally, the technique is contrasted against other more familiar control design techniques, like PID control, to show its advantages.

Kopasakis, George↗

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.↗

On identifying transfer functions and state equations for linear systems.

Two methods are established for identifying constant-coefficient, C to the 2n power type of noise-free linear systems if the time response data of the input-output or of all states are known. 2n response data are required to identify an nth-order transfer function or state equation for an unknown linear system. The order of the unknown system can be identified by checking a sequence of determinants. The Z transform and its inversion are mainly used.

Shieh, L. S.↗

A new stream function formulation for the Euler equations

A new stream function formulation is developed for the solution of Euler's equations in the transonic flow region. The stream function and the density are the dependent variables in this method, while the governing equations for adiabatic flow are the momentum equations which are solved in the strong conservation law form. The application of this method does not require a knowledge of the vorticity. The algorithm is combined with the automatic grid solver (GRAPE) of Steger and Sorenson (1979) in order to study arbitrary geometries. Results of the application of this method are presented for the NACA 0012 airfoil at various Mach numbers and angles of attack, and cylinders. In addition, detailed comparisons are made with other solutions of the Euler equations.

Atkins, H. L.↗

Correction function in the Lidar equation and the solution techniques for CO2 Lidar date reduction

For lidar systems with long laser pulses the unusual behavior of the near-range signals causes serious difficulties and large errors in reduction. The commonly used lidar equation is no longer applicable since the convolution of the laser pulse with the atmospheric parameter distributions should be taken into account. It is important to give more insight into this problem and find the solution techniques. Starting from the original equation, a general form is suggested for the single scattering lidar equation where a correction function Cr is introduced. The correction Function Cr(R) derived from the original equation indicates the departure from the normal lidar equation. Examples of Cr(R) for a coaxial CO2 lidar system are presented. The Differential Absorption Lidar (DIAL) errors caused by the differences of Cr(R) for H2O measurements are plotted against height.

Zhao, Y.↗