Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Arbitrary functions”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Exact Integrations of Polynomials and Symmetric Quadrature Formulas over Arbitrary Polyhedral Grids

This paper is concerned with two important elements in the high-order accurate spatial discretization of finite volume equations over arbitrary grids. One element is the integration of basis functions over arbitrary domains, which is used in expressing various spatial integrals in terms of discrete unknowns. The other consists of quadrature approximations to those integrals. Only polynomial basis functions applied to polyhedral and polygonal grids are treated here. Non-triangular polygonal faces are subdivided into a union of planar triangular facets, and the resulting triangulated polyhedron is subdivided into a union of tetrahedra. The straight line segment, triangle, and tetrahedron are thus the fundamental shapes that are the building blocks for all integrations and quadrature approximations. Integrals of products up to the fifth order are derived in a unified manner for the three fundamental shapes in terms of the position vectors of vertices. Results are given both in terms of tensor products and products of Cartesian coordinates. The exact polynomial integrals are used to obtain symmetric quadrature approximations of any degree of precision up to five for arbitrary integrals over the three fundamental domains. Using a coordinate-free formulation, simple and rational procedures are developed to derive virtually all quadrature formulas, including some previously unpublished. Four symmetry groups of quadrature points are introduced to derive Gauss formulas, while their limiting forms are used to derive Lobatto formulas. Representative Gauss and Lobatto formulas are tabulated. The relative efficiency of their application to polyhedral and polygonal grids is detailed. The extension to higher degrees of precision is discussed.

Liu, Yen↗
NASA NTRS1981

Sound diffraction around screens and wedges for arbitrary point source locations

A development is presented for the Green's function for a point source in the vicinity of a rigid wedge. The diffraction contributions to the Green's function for arbitrary source and listener location is expressed in a form which can be readily evaluated using the Laguerre technique for numerical integration. The present approach offers the advantages of efficient numerical evaluation and of relatively straightforward reduction to well-known analytical approximations in limiting cases. Comparisons with previously obtained experimental and numerical results obtained by Ambaud and Bergassoli (1972) are presented. The comparison with the experimental results is excellent; the advantages of the present numerical technique, vis a vis that of Ambaud and Bergassoli, are pointed out.

Hadden, W. J., Jr.↗
NASA NTRS1982

Toward autonomous spacecraft

Ways in which autonomous behavior of spacecraft can be extended to treat situations wherein a closed loop control by a human may not be appropriate or even possible are explored. Predictive models that minimize mean least squared error and arbitrary cost functions are discussed. A methodology for extracting cyclic components for an arbitrary environment with respect to usual and arbitrary criteria is developed. An approach to prediction and control based on evolutionary programming is outlined. A computer program capable of predicting time series is presented. A design of a control system for a robotic dense with partially unknown physical properties is presented.

Fogel, L. J.↗
NASA NTRS1972

Response of discrete linear systems to forcing functions with inequality constraints.

An analysis is made of the maximum response of discrete, linear mechanical systems to arbitrary forcing functions which lie within specified bounds. Primary attention is focused on the complete determination of the forcing function which will engender maximum displacement to any particular mass element of a multi-degree-of-freedom system. In general, the desired forcing function is found to be a bang-bang type function, i.e., a function which switches from the maximum to the minimum bound and vice-versa at certain instants of time. Examples of two-degree-of-freedom systems, with and without damping, are presented in detail. Conclusions are drawn concerning the effect of damping on the switching times and the general procedure for finding these times is discussed.

Michalopoulos, C. D.↗
NASA NTRS1984

Bidirectional reflectance spectroscopy. III - Correction for macroscopic roughness

A mathematically rigorous formalism is derived by which an arbitrary photometric function for the bidirectional reflectance of a smooth surface may be corrected to include effects of general macroscopic roughness. The correction involves only one arbitrary parameter, the mean slope angle, and is applicable to surfaces of any albedo. Using physically reasonable assumptions and mathematical approximations, the correction expressions are evaluated analytically to second order in the mean slope angle. The correction is applied to the bidirectional-reflectance function of Hapke (1981). Expressions for both the differential and integral brightnesses are obtained. Photometric profiles on hypothetical smooth and rough planets of low and high albedo are shown to illustrate the effects of macroscopic roughness. The theory is applied to observations of Mercury and predicts the integral phase function, the apparent polar darkening, and the lack of limb brightness surge on the planet. The roughness-corrected bidirectional-reflectance function is sufficiently simple that it can be conveniently evaluated on a programmable hand-held calculator.

Hapke, B.↗
NASA NTRS1995

A 3D High-Order Unstructured Finite-Volume Algorithm for Solving Maxwell's Equations

A three-dimensional finite-volume algorithm based on arbitrary basis functions for time-dependent problems on general unstructured grids is developed. The method is applied to the time-domain Maxwell equations. Discrete unknowns are volume integrals or cell averages of the electric and magnetic field variables. Spatial terms are converted to surface integrals using the Gauss curl theorem. Polynomial basis functions are introduced in constructing local representations of the fields and evaluating the volume and surface integrals. Electric and magnetic fields are approximated by linear combinations of these basis functions. Unlike other unstructured formulations used in Computational Fluid Dynamics, the new formulation actually does not reconstruct the field variables at each time step. Instead, the spatial terms are calculated in terms of unknowns by precomputing weights at the beginning of the computation as functions of cell geometry and basis functions to retain efficiency. Since no assumption is made for cell geometry, this new formulation is suitable for arbitrarily defined grids, either smooth or unsmooth. However, to facilitate the volume and surface integrations, arbitrary polyhedral cells with polygonal faces are used in constructing grids. Both centered and upwind schemes are formulated. It is shown that conventional schemes (second order in Cartesian grids) are equivalent to the new schemes using first degree polynomials as the basis functions and the midpoint quadrature for the integrations. In the new formulation, higher orders of accuracy are achieved by using higher degree polynomial basis functions. Furthermore, all the surface and volume integrations are carried out exactly. Several model electromagnetic scattering problems are calculated and compared with analytical solutions. Examples are given for cases based on 0th to 3rd degree polynomial basis functions. In all calculations, a centered scheme is applied in the interior, while an upwind matching scheme is employed at material interfaces and the Engquist-Majda non-reflecting boundary condition is implemented at the numerical outer boundaries. The staggered leapfrog scheme and the Runge-Kutta methods are utilized for the time integration. Excellent agreements are found between the numerical and analytical solutions.

Liu, Yen↗
NASA NTRS1992

Enhancement of atmospheric radiation by an aerosol layer

The presence of a stratospheric haze layer may produce increases in both the actinic flux and the irradiance below this layer. Such haze layers result from the injection of aerosol-forming material into the stratosphere by volcanic eruptions. Simple heuristic arguments show that the increase in flux below the haze layer, relative to a clear sky case, is a consequence of 'photon trapping'. The magnitude of these flux perturbations, as a function of aerosol properties and illumination conditions, is explored with a new radiative transfer model that can accurately compute fluxes in an inhomogeneous atmosphere with nonconservative scatterers having arbitrary phase function. One calculated consequence of the El Chichon volcanic eruption is an increase in the midday surface actinic flux at 20 deg N latitude, summer, by as much as 45 percent at 2900 A. This increase in flux in the UV-B wavelength range was caused entirely by aerosol scattering, without any reduction in the overhead ozone column.

Michelangeli, Diane V.↗
NASA NTRS2001

Hippocampal place-cell firing during movement in three-dimensional space

"Place" cells of the rat hippocampus are coupled to "head direction" cells of the thalamus and limbic cortex. Head direction cells are sensitive to head direction in the horizontal plane only, which leads to the question of whether place cells similarly encode locations in the horizontal plane only, ignoring the z axis, or whether they encode locations in three dimensions. This question was addressed by recording from ensembles of CA1 pyramidal cells while rats traversed a rectangular track that could be tilted and rotated to different three-dimensional orientations. Cells were analyzed to determine whether their firing was bound to the external, three-dimensional cues of the environment, to the two-dimensional rectangular surface, or to some combination of these cues. Tilting the track 45 degrees generally provoked a partial remapping of the rectangular surface in that some cells maintained their place fields, whereas other cells either gained new place fields, lost existing fields, or changed their firing locations arbitrarily. When the tilted track was rotated relative to the distal landmarks, most place fields remapped, but a number of cells maintained the same place field relative to the x-y coordinate frame of the laboratory, ignoring the z axis. No more cells were bound to the local reference frame of the recording apparatus than would be predicted by chance. The partial remapping demonstrated that the place cell system was sensitive to the three-dimensional manipulations of the recording apparatus. Nonetheless the results were not consistent with an explicit three-dimensional tuning of individual hippocampal neurons nor were they consistent with a model in which different sets of cells are tightly coupled to different sets of environmental cues. The results are most consistent with the statement that hippocampal neurons can change their "tuning functions" in arbitrary ways when features of the sensory input or behavioral context are altered. Understanding the rules that govern the remapping phenomenon holds promise for deciphering the neural circuitry underlying hippocampal function.

Non-NASA Center↗
NASA NTRS2011

Finding Every Root of a Broad Class of Real, Continuous Functions in a Given Interval

One of the most pervasive needs within the Deep Space Network (DSN) Metric Prediction Generator (MPG) view period event generation is that of finding solutions to given occurrence conditions. While the general form of an equation expresses equivalence between its left-hand and right-hand expressions, the traditional treatment of the subject subtracts the two sides, leaving an expression of the form Integral of(x) = 0. Values of the independent variable x satisfying this condition are roots, or solutions. Generally speaking, there may be no solutions, a unique solution, multiple solutions, or a continuum of solutions to a given equation. In particular, all view period events are modeled as zero crossings of various metrics; for example, the time at which the elevation of a spacecraft reaches its maximum value, as viewed from a Deep Space Station (DSS), is found by locating that point at which the derivative of the elevation function becomes zero. Moreover, each event type may have several occurrences within a given time interval of interest. For example, a spacecraft in a low Moon orbit will experience several possible occultations per day, each of which must be located in time. The MPG is charged with finding all specified event occurrences that take place within a given time interval (or pass ), without any special clues from operators as to when they may occur, for the entire spectrum of missions undertaken by the DSN. For each event type, the event metric function is a known form that can be computed for any instant within the interval. A method has been created for a mathematical root finder to be capable of finding all roots of an arbitrary continuous function, within a given interval, to be subject to very lenient, parameterized assumptions. One assumption is that adjacent roots are separated at least by a given amount, xGuard. Any point whose function value is less than ef in magnitude is considered to be a root, and the function values at distances xGuard away from a root are larger than ef, unless there is another root located in this vicinity. A root is considered found if, during iteration, two root candidates differ by less than a pre-specified ex, and the optimum cubic polynomial matching the function at the end and at two interval points (that is within a relative error fraction L at its midpoint) is reliable in indicating whether the function has extrema within the interval. The robustness of this method depends solely on choosing these four parameters that control the search. The roots of discontinuous functions were also found, but at degraded performance.

Tausworthe, Robert C.↗
NASA NTRS1967

Computer program utilizes FORTRAN 4 subroutines for contour plotting

Computer program constructs lists of xy-coordinate pairs that define contour curves for an arbitrary given function of two variables and transmits these lists to plotting equipment to produce contour plots. The principal subroutine, CONTUR, is independent of any specific system of plotting subroutines and equipment.

Block, N.↗
NASA NTRS1973

Response of a rigid aircraft to nonstationary atmospheric turbulence.

The plunging response of an aircraft to a type of nonstationary turbulent excitation is considered. The latter consists of stationary Gaussian noise modulated by a well-defined envelope function. The intent of the investigation is to model the excitation experienced by an airplane flying through turbulence of varying intensity and to examine the influence of intensity variations on exceedance frequencies of the gust velocity and the airplane's plunging velocity and acceleration. One analytical advantage of the proposed model is that the Gaussian assumption for the gust excitation is retained. The analysis described herein is developed in terms of an envelope function of arbitrary form; however, numerical calculations are limited to the case of harmonic modulation.

Verdon, J. M.↗
NASA NTRS1975

JPL transient radiation analysis by computer program (JTRAC)

Digital computer program, JTRAC, simulates time response of electronic circuit to arbitrary forcing functions which may include electrical and/or radiation stimuli. Program is designed to solve linear and nonlinear simultaneous equations which characterize mathematical models used to predict circuit response for electrical and/or radiation input.

Weinstein, S.↗
NASA NTRS1976

Polyatomic molecule vibrations

Polyatomic molecule vibrations are analyzed as harmonic vibrations along normal coordinates. The energy eigenvalues are found for linear and nonlinear symmetric triatomic molecules for valence bond models of the potential function with arbitrary coupling coefficients; such models can usually be fitted to observed energy levels with reasonably good accuracy. Approximate normal coordinates for the H2O molecule are discussed. Degenerate vibrational modes such as occur in CO2 are analyzed and expressions for Fermi resonance between close-lying states of the same symmetry are developed. The bending modes of linear triatomic molecules are expressed in terms of Laguerre polynomials in cylindrical coordinates as well as in terms of Hermite polynomials in Cartesian coordinates. The effects of large-amplitude bending such as occur in the C3 molecule are analyzed, along with anharmonic effects, which split the usually degenerate bending mode energy levels. Finally, the vibrational frequencies, degeneracies, and symmetry properties of XY3, X2Y2, and XY4 type molecules are discussed.

Source record↗
NASA NTRS1977

The general solution to the classical problem of finite Euler Bernoulli beam

An analytical solution is obtained for the problem of free and forced vibrations of a finite Euler Bernoulli beam with arbitrary (partially fixed) boundary conditions. The effects of linear viscous damping, Winkler foundation, constant axial tension, a concentrated mass, and an arbitrary forcing function are included in the analysis. No restriction is placed on the values of the parameters involved, and the solution presented here contains all cited previous solutions as special cases.

Hussaini, M. Y.↗