Search NASA⌕ Search

SEARCH · Search NASA

Results for “RATIONAL FUNCTION”

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 145 records · Page 8

GENIE(++): A Multi-Block Structured Grid System

The computer code GENIE++ is a continuously evolving grid system containing a multitude of proven geometry/grid techniques. The generation process in GENIE++ is based on an earlier version. The process uses several techniques either separately or in combination to quickly and economically generate sculptured geometry descriptions and grids for arbitrary geometries. The computational mesh is formed by using an appropriate algebraic method. Grid clustering is accomplished with either exponential or hyperbolic tangent routines which allow the user to specify a desired point distribution. Grid smoothing can be accomplished by using an elliptic solver with proper forcing functions. B-spline and Non-Uniform Rational B-splines (NURBS) algorithms are used for surface definition and redistribution. The built in sculptured geometry definition with desired distribution of points, automatic Bezier curve/surface generation for interior boundaries/surfaces, and surface redistribution is based on NURBS. Weighted Lagrance/Hermite transfinite interpolation methods, interactive geometry/grid manipulation modules, and on-line graphical visualization of the generation process are salient features of this system which result in a significant time savings for a given geometry/grid application.

Williams, Tonya↗

Partial Reflection D-region Electron Densities

The differential absorption technique of measuring electron densities as a function of height in the D region is discussed. In the basic experiment, pulses of medium or high frequency, usually at a fixed frequency (2 to MHz), are radiated upwards with known wave polarizatin (usually linear or circular) from a transmitter at ground level. Partial reflections, from ionospheric scatterers at heights below the E region, are received at the ground, and are resolved into two characteristic components, the ordinary (0) and extraordinary (E) modes whose amplitude ration A(x)/A(o) is then measured as a function of height, h. The heights of these are determined by delay times, the group retardation being minimal in the undisturbed D-region. The electronic system can be very simple. Power splitters and quadrature networks to separate the A(x) and A(o) components are commercially available at low prices and an A-D converter, height-gate system, and microcomputer allows the real-time calculation of mean amplitudes. The ratio of the coefficients of reflection of the two modes, as they originate at each reflection height is then calculable.

Manson, A. H.↗

Round-off error in long-term orbital integrations using multistep methods

Techniques for reducing roundoff error are compared by testing them on high-order Stormer and summetric multistep methods. The best technique for most applications is to write the equation in summed, function-evaluation form and to store the coefficients as rational numbers. A larger error reduction can be achieved by writing the equation in backward-difference form and performing some of the additions in extended precision, but this entails a larger central processing unit (cpu) cost.

Quinlan, Gerald D.↗

A Rayleigh-Ritz approach to the synthesis of large structures with rotating flexible components

The equations of motion for large structures with rotating flexible components are derived by regarding the structure as an assemblage of substructures. Based on a stationarity principle for rotating structures, it is shown that each continuous or discrete substructure can be simulated by a suitable set of admissible functions or admissible vectors. This substructure synthesis approach provides a rational basis for truncating the number of degrees of freedom both of each substructure and of the assembled structure.

Meirovitch, L.↗

A Rayleigh-Ritz approach to the synthesis of large gyroscopic structures

The equations of motion for large structures with rotating flexible components are derived by regarding the structure as an assemblage of substructures. Based on a stationarity principle for rotating structures, it is shown that each continuous or discrete substructure can be simulated by a suitable set of admissible functions or admissible vectors. This substructure synthesis approach provides a rational basis for truncating the number of degrees of freedom both of each substructure and of the assembled structure.

Meirovitch, L.↗

Accurate numerical evaluation of modified Struve functions occurring in unsteady aerodynamics

A method is developed for computing the modified Struve functions that occur in unsteady aerodynamics. The method uses a rational approximation supplemented by an asymptotic series for large argument. Simple recursive formulas for generating the coefficients are derived. The method is capable of generating results of arbitrary accuracy. It can also be used for complex argument and order. For greater computing speed, a method is presented that uses the rational and asymptotic approximations to generate Chebyshev coefficients.

Desmarais, R. N.↗

Computational Models of the Cardiovascular System and Its Response to Microgravity

Computational models of the cardiovascular system are powerful adjuncts to ground-based and in-flight experiments. We will provide NSBRI with a model capable of simulating the short-term effects of gravity on cardiovascular function. The model from this project will: (1) provide a rational framework which quantitatively defines interactions among complex cardiovascular parameters and which supports the critical interpretation of experimental results and testing of hypotheses. (2) permit predictions of the impact of specific countermeasures in the context of various hypothetical cardiovascular abnormalities induced by microgravity. Major progress has been made during the first 18 months of the program: (1) We have developed an operational first-order computer model capable of simulating the cardiovascular response to orthostatic stress. The model consists of a lumped parameter hemodynamic model and a complete reflex control system. The latter includes cardiopulmonary and carotid sinus reflex limbs and interactions between the two. (2) We have modeled the physiologic stress of tilt table experiments and lower body negative pressure procedures (LBNP). We have verified our model's predictions by comparing them with experimental findings from the literature. (3) We have established collaborative efforts with leading investigators interested in experimental studies of orthostatic intolerance, cardiovascular control, and physiologic responses to space flight. (4) We have established a standardized method of transferring data to our laboratory from the ongoing NSBRI bedrest studies. We use this data to estimate input parameters to our model and compare our model predictions to actual data to further verify our model. (5) We are in the process of systematically simulating current hypotheses concerning the mechanism underlying orthostatic intolerance by matching our simulations to stand test data from astronauts pre- and post-flight. (6) We are in the process of developing a JAVA version of the simulator which will be distributed amongst the cardiovascular team members. Future work on this project involves modifications of the model to represent a rodent (rat) model, further evaluation of the bedrest astronaut and animal data, and systematic investigation of specific countermeasures.

Kamm, Roger D.↗

Application of vector-valued rational approximations to the matrix eigenvalue problem and connections with Krylov subspace methods

Let F(z) be a vectored-valued function F: C approaches C sup N, which is analytic at z=0 and meromorphic in a neighborhood of z=0, and let its Maclaurin series be given. We use vector-valued rational approximation procedures for F(z) that are based on its Maclaurin series in conjunction with power iterations to develop bona fide generalizations of the power method for an arbitrary N X N matrix that may be diagonalizable or not. These generalizations can be used to obtain simultaneously several of the largest distinct eigenvalues and the corresponding invariant subspaces, and present a detailed convergence theory for them. In addition, it is shown that the generalized power methods of this work are equivalent to some Krylov subspace methods, among them the methods of Arnoldi and Lanczos. Thus, the theory provides a set of completely new results and constructions for these Krylov subspace methods. This theory suggests at the same time a new mode of usage for these Krylov subspace methods that were observed to possess computational advantages over their common mode of usage.

Sidi, Avram↗

Macular Bioaccelerometers on Earth and in Space

Space flight offers the opportunity to study linear bioaccelerometers (vestibular maculas) in the virtual absence of a primary stimulus, gravitational acceleration. Macular research in space is particularly important to NASA because the bioaccelerometers are proving to be weighted neural networks in which information is distributed for parallel processing. Neural networks are plastic and highly adaptive to new environments. Combined morphological-physiological studies of maculas fixed in space and following flight should reveal macular adaptive responses to microgravity, and their time-course. Ground-based research, already begun, using computer-assisted, 3-dimensional reconstruction of macular terminal fields will lead to development of computer models of functioning maculas. This research should continue in conjunction with physiological studies, including work with multichannel electrodes. The results of such a combined effort could usher in a new era in understanding vestibular function on Earth and in space. They can also provide a rational basis for counter-measures to space motion sickness, which may prove troublesome as space voyager encounter new gravitational fields on planets, or must re-adapt to 1 g upon return to earth.

Ross, M. D.↗

Modeling Humans as Reinforcement Learners: How to Predict Human Behavior in Multi-Stage Games

This paper introduces a novel framework for modeling interacting humans in a multi-stage game environment by combining concepts from game theory and reinforcement learning. The proposed model has the following desirable characteristics: (1) Bounded rational players, (2) strategic (i.e., players account for one anothers reward functions), and (3) is computationally feasible even on moderately large real-world systems. To do this we extend level-K reasoning to policy space to, for the first time, be able to handle multiple time steps. This allows us to decompose the problem into a series of smaller ones where we can apply standard reinforcement learning algorithms. We investigate these ideas in a cyber-battle scenario over a smart power grid and discuss the relationship between the behavior predicted by our model and what one might expect of real human defenders and attackers.

Lee, Ritchie↗

A highly accurate Voigt function algorithm

A complex Voigt lineshape algorithm is presented whose maximum relative error over the complex plane is less than 1 x 10(exp -8). The algorithm consists of series, rational approximations and Gauss-Hermite integrations which makes it suitable as a general purpose software module for a wide variety of uses, including a Voigt function standard.

Shippony, Z.↗

A Highly Accurate Voight Function Algorithm

A complex Voight lineshape algorithm is presented whose maximum relative error over the complex plane is less than 1 x 10 super minus 8. The algorithm consists of series, rational approximations and Gauss-Hermite integrations which makes it suitable as a general purpose software module for a wide variety of uses, including a Voight function standard.

Voight lineshape Doppler velocity distribution col↗

l1-optimal control of multivariable systems with output norm constraints

This paper considers the l1-optimal control problem for general rational plants. It is shown that, for plants with no poles or zeros on the unit circle, an optimal compensator exists and that the resulting closed loop transfer function is polynomial whenever there are at least as many controls as regulated outputs and at least as many measurements as exogeneous inputs. Exactly or approximately, optimal rational compensators can be obtained by solving a sequence of finite linear programs for the coefficients of a polynomial closed-loop transfer function. No assumptions on plant poles or zeros are required to obtain at least approximately optimal compensators. It is shown that constrained problems in which a set of outputs is regulated subject to l(infinity)-norm constraints on another set of outputs can be solved using a slight modification of the same algorithm.

Mcdonald, J. S.↗

Algebraic grid adaptation method using non-uniform rational B-spline surface modeling

An algebraic adaptive grid system based on equidistribution law and utilized by the Non-Uniform Rational B-Spline (NURBS) surface for redistribution is presented. A weight function, utilizing a properly weighted boolean sum of various flow field characteristics is developed. Computational examples are presented to demonstrate the success of this technique.

Yang, Jiann-Cherng↗

Water vapor variance measurements using a Raman lidar

Because of the importance of atmospheric water vapor variance, we have analyzed data from the NASA/Goddard Raman lidar to obtain temporal scales of water vapor mixing ratio as a function of altitude over observation periods extending to 12 hours. The ground-based lidar measures water vapor mixing ration from near the earth's surface to an altitude of 9-10 km. Moisture profiles are acquired once every minute with 75 m vertical resolution. Data at each 75 meter altitude level can be displayed as a function of time from the beginning to the end of an observation period. These time sequences have been spectrally analyzed using a fast Fourier transform technique. An example of such a temporal spectrum obtained between 00:22 and 10:29 UT on December 6, 1991 is shown in the figure. The curve shown on the figure represents the spectral average of data from 11 height levels centered on an altitude of 1 km (1 plus or minus .375 km). The spectra shows a decrease in energy density with frequency which generally follows a -5/3 power law over the spectral interval 3x10 (exp -5) to 4x10 (exp -3) Hz. The flattening of the spectrum for frequencies greater than 6x10 (exp -3) Hz is most likely a measure of instrumental noise. Spectra like that shown in the figure are calculated for other altitudes and show changes in spectral features with height. Spectral analysis versus height have been performed for several observation periods which demonstrate changes in water vapor mixing ratio spectral character from one observation period to the next. The combination of these temporal spectra with independent measurements of winds aloft provide an opportunity to infer spatial scales of moisture variance.

Evans, K.↗

Module theoretic zero structures for system matrices

The coordinate-free module-theoretic treatment of transmission zeros for MIMO transfer functions developed by Wyman and Sain (1981) is generalized to include noncontrollable and nonobservable linear dynamical systems. Rational, finitely-generated-modular, and torsion-divisible interpretations of the Rosenbrock system matrix are presented; Gamma-zero and Omega-zero modules are defined and shown to contain the output-decoupling and input-decoupling zero modules, respectively, as submodules; and the cases of left and right invertible transfer functions are considered.

Wyman, Bostwick F.↗

System evaluation from impulse response fragments

By use of rational approximation, a procedure was found which allows calculation of the z-transfer function from a restricted number of time sequence values. The order of the system need not be known. Transformation to the s-domain is easily provided.

Baumeister, J. K.↗