Search NASA⌕ Search

SEARCH · Search NASA

Results for “Polynomials”

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 415 records · Page 23

Positivity-preserving numerical schemes for multidimensional advection

This report describes the construction of an explicit, single time-step, conservative, finite-volume method for multidimensional advective flow, based on a uniformly third-order polynomial interpolation algorithm (UTOPIA). Particular attention is paid to the problem of flow-to-grid angle-dependent, anisotropic distortion typical of one-dimensional schemes used component-wise. The third-order multidimensional scheme automatically includes certain cross-difference terms that guarantee good isotropy (and stability). However, above first-order, polynomial-based advection schemes do not preserve positivity (the multidimensional analogue of monotonicity). For this reason, a multidimensional generalization of the first author's universal flux-limiter is sought. This is a very challenging problem. A simple flux-limiter can be found; but this introduces strong anisotropic distortion. A more sophisticated technique, limiting part of the flux and then restoring the isotropy-maintaining cross-terms afterwards, gives more satisfactory results. Test cases are confined to two dimensions; three-dimensional extensions are briefly discussed.

Leonard, B. P.↗

Design of an essentially non-oscillatory reconstruction procedure in finite-element type meshes

An essentially non oscillatory reconstruction for functions defined on finite element type meshes is designed. Two related problems are studied: the interpolation of possibly unsmooth multivariate functions on arbitary meshes and the reconstruction of a function from its averages in the control volumes surrounding the nodes of the mesh. Concerning the first problem, the behavior of the highest coefficients of two polynomial interpolations of a function that may admit discontinuities of locally regular curves is studied: the Lagrange interpolation and an approximation such that the mean of the polynomial on any control volume is equal to that of the function to be approximated. This enables the best stencil for the approximation to be chosen. The choice of the smallest possible number of stencils is addressed. Concerning the reconstruction problem, two methods were studied: one based on an adaptation of the so called reconstruction via deconvolution method to irregular meshes and one that lies on the approximation on the mean as defined above. The first method is conservative up to a quadrature formula and the second one is exactly conservative. The two methods have the expected order of accuracy, but the second one is much less expensive than the first one. Some numerical examples are given which demonstrate the efficiency of the reconstruction.

Abgrall, Remi↗

A fast algorithm for image-based ranging

Image-based ranging has emerged as a critical issue in the low altitude operation of flight vehicles such as rotorcraft and planetary landers. These flight regimes require ranging systems for recovering the geometry of the terrain and obstacles for use with guidance algorithms. The development of a ranging equation combining image irradiance together with various order spatial partial derivatives and the vehicle motion parameters is discussed. The ranging equation is in the form of a polynomial in scene depth. Two-dimensional linear filters are then used to compute the coefficients of this polynomial to result in a fast image-based ranging algorithm. Performance of the algorithm is demonstrated using laboratory images.

Menon, P. K. A.↗

Random interactions in higher order neural networks

Recurrent networks of polynomial threshold elements with random symmetric interactions are studied. Precise asymptotic estimates are derived for the expected number of fixed points as a function of the margin of stability. In particular, it is shown that there is a critical range of margins of stability (depending on the degree of polynomial interaction) such that the expected number of fixed points with margins below the critical range grows exponentially with the number of nodes in the network, while the expected number of fixed points with margins above the critical range decreases exponentially with the number of nodes in the network. The random energy model is also briefly examined and links with higher order neural networks and higher order spin glass models made explicit.

Baldi, Pierre↗

Symbolic discrete event system specification

Extending discrete event modeling formalisms to facilitate greater symbol manipulation capabilities is important to further their use in intelligent control and design of high autonomy systems. An extension to the DEVS formalism that facilitates symbolic expression of event times by extending the time base from the real numbers to the field of linear polynomials over the reals is defined. A simulation algorithm is developed to generate the branching trajectories resulting from the underlying nondeterminism. To efficiently manage symbolic constraints, a consistency checking algorithm for linear polynomial constraints based on feasibility checking algorithms borrowed from linear programming has been developed. The extended formalism offers a convenient means to conduct multiple, simultaneous explorations of model behaviors. Examples of application are given with concentration on fault model analysis.

Zeigler, Bernard P.↗

Nonlinear aerodynamic modeling using multivariate orthogonal functions

A technique was developed for global modeling of nonlinear aerodynamic coefficients using multivariate orthogonal functions based on the data. Each orthogonal function retained in the model was decomposed into an expansion of ordinary polynomials in the independent variables, so that the final model could be interpreted as selectively retained terms from a multivariable power series expansion. A predicted squared-error metric was used to determine the orthogonal functions to be retained in the model; analytical derivatives were easily computed. The approach was demonstrated on the Z-body axis aerodynamic force coefficient (Cz) wind tunnel data for an F-18 research vehicle which came from a tabular wind tunnel and covered the entire subsonic flight envelope. For a realistic case, the analytical model predicted experimental values of Cz very well. The modeling technique is shown to be capable of generating a compact, global analytical representation of nonlinear aerodynamics. The polynomial model has good predictive capability, global validity, and analytical differentiability.

Morelli, Eugene A.↗

Importance of Geosat orbit and tidal errors in the estimation of large-scale Indian Ocean variations

To improve the estimate accuracy of large-scale meridional sea-level variations, Geosat ERM data on the Indian Ocean for a 26-month period were processed using two different techniques of orbit error reduction. The first technique removes an along-track polynomial of degree 1 over about 5000 km and the second technique removes an along-track once-per-revolution sine wave about 40,000 km. Results obtained show that the polynomial technique produces stronger attenuation of both the tidal error and the large-scale oceanic signal. After filtering, the residual difference between the two methods represents 44 percent of the total variance and 23 percent of the annual variance. The sine-wave method yields a larger estimate of annual and interannual meridional variations.

Perigaud, Claire↗

Generalized indical forces on deforming rectangular wings in supersonic flight

A method is presented for determining the time-dependent flow over a rectangular wing moving with a supersonic forward speed and undergoing small vertical distortions expressible as polynomials involving spanwise and chordwise distances. The solution for the velocity potential is presented in a form analogous to that for steady supersonic flow having the familiar "reflected area" concept discovered by Evvard. Particular attention is paid to indicial-type motions and results are expressed in terms of generalized indicial forces. Numerical results for Mach numbers equal to 1.1 and 1.2 are given for polynomials of the first and fifth degree in the chordwise and spanwise directions, respectively, on a wing having an aspect ratio of 4.

Lomax, Harvard↗

On well-partial-order theory and its application to combinatorial problems of VLSI design

We nonconstructively prove the existence of decision algorithms with low-degree polynomial running times for a number of well-studied graph layout, placement, and routing problems. Some were not previously known to be in p at all; others were only known to be in p by way of brute force or dynamic programming formulations with unboundedly high-degree polynomial running times. Our methods include the application of the recent Robertson-Seymour theorems on the well-partial-ordering of graphs under both the minor and immersion orders. We also briefly address the complexity of search versions of these problems.

Fellows, M.↗

Application specific serial arithmetic arrays

High performance systolic arrays of serial-parallel multiplier elements may be rapidly constructed for specific applications by applying hardware description language techniques to a library of full-custom CMOS building blocks. Single clock pre-charged circuits have been implemented for these arrays at clock rates in excess of 100 Mhz using economical 2-micron (minimum feature size) CMOS processes, which may be quickly configured for a variety of applications. A number of application-specific arrays are presented, including a 2-D convolver for image processing, an integer polynomial solver, and a finite-field polynomial solver.

Winters, K.↗

Flutter analysis using transversality theory

A new method of calculating flutter boundaries of undamped aeronautical structures is presented. The method is an application of the weak transversality theorem used in catastrophe theory. In the first instance, the flutter problem is cast in matrix form using a frequency domain method, leading to an eigenvalue matrix. The characteristic polynomial resulting from this matrix usually has a smooth dependence on the system's parameters. As these parameters change with operating conditions, certain critical values are reached at which flutter sets in. Our approach is to use the transversality theorem in locating such flutter boundaries using this criterion: at a flutter boundary, the characteristic polynomial does not intersect the axis of the abscissa transversally. Formulas for computing the flutter boundaries and flutter frequencies of structures with two degrees of freedom are presented, and extension to multi-degree of freedom systems is indicated. The formulas have obvious applications in, for instance, problems of panel flutter at supersonic Mach numbers.

Afolabi, D.↗

Gravity investigation of the Manson impact structure, Iowa

The Manson crater, of probable Cretaceous/Tertiary age, is located in northwestern Iowa (center at 42 deg. 34.44 min N; 94 deg. 33.60 min W). A seismic reflection profile along an east west line across the crater and drill hole data indicate a crater about 35 km in diameter having the classic form for an impact crater, an uplifted central peak composed of uplifted Proterozoic crystalline bedrock, surrounded by a 'moat' filled with impact produced breccia and a ring graben zone composed of tilted fault blocks of the Proterozoic and Paleozoic country rocks. The structure has been significantly eroded. This geologic structure would be expected to produce a significant gravity signature and study of that signature would shed additional light on the details of the crater structure. A gravity study was undertaken to better resolve the crustal structure. The regional Bouguer gravity field is characterized by a southeastward decreasing field. To first order, the Bouguer gravity field can be understood in the context of the geology of the Precambrian basement. The high gravity at the southeast corner is associated with the mid-continent gravity high; the adjacent low to the northwest results from a basin containing low-density clastic sediments shed from the basement high. Modeling of a simple basin and adjacent high predicts much of the observed Bouguer gravity signature. A gravity signature due to structure associated with the Manson impact is not apparent in the Bouguer data. To resolve the gravity signature of the impact, a series of polynomial surfaces were fit to the Bouguer gravity field to isolate the small wavelength residual anomalies. The residual gravity obtained after subtracting a 5th- or 6th-order polynomial seems to remove most of the regional effects and isolate local anomalies. The pattern resolved in the residual gravity is one of a gravity high surrounded by gravity lows and in turn surrounded by isolated gravity highs. The central portion of the crater is characterized by two positive anomalies having amplitudes of about plus 4 mGal separated by a gentle saddle located approximately at the crater center.

Plescia, J. B.↗

On the Gibbs phenomenon 3: Recovering exponential accuracy in a sub-interval from a spectral partial sum of a piecewise analytic function

The investigation of overcoming Gibbs phenomenon was continued, i.e., obtaining exponential accuracy at all points including at the discontinuities themselves, from the knowledge of a spectral partial sum of a discontinuous but piecewise analytic function. It was shown that if we are given the first N expansion coefficients of an L(sub 2) function f(x) in terms of either the trigonometrical polynomials or the Chebyshev or Legendre polynomials, an exponentially convergent approximation to the point values of f(x) in any sub-interval in which it is analytic can be constructed.

Gottlieb, David↗

In-flight scale/distortion calibration of the Hubble Space Telescope fixed-head star trackers

This paper describes an in-flight scale and distortion calibration procedure that has been developed for the Ball Aerospace Systems Division Fixed-Head Star Trackers (FHST's) used on the Hubble Space Telescope (HST). The FHST is a magnetically focused and deflected imaging sensor that is designed to track stars as faint as m(sub v) = 5.7 over an 8 degree by 8 degree field of view. Raw FHST position measurements are accurate to approximately 200 arcseconds, but this can be improved to 10-15 arcseconds by processing the raw measurements through calibration polynomials that correct for flat field, temperature intensity, and magnetic field effects. The coefficients for these polynomials were initially determined using ground test data. On HST the use of three FHST's is an integral part of the preliminary attitude update procedures required before the acquisition of guide stars for science observations. To this end, FHST-based attitude determination having single-axis errors no worse than 22 arcseconds (1 sigma) is required. In early 1991 it became evident that one of the HST FHST's was experiencing a significant change in its optical scale. By mid-1993 the size of this error had grown to a point that, if not corrected, it would correspond to a maximum position error on the order of 100 arcseconds. Subsequent investigations demonstrated that substantial, uncompensated cubic distortion effects had also developed, the maximum contribution to position errors from the cubic terms being on the order of 30 arcseconds. To ensure accurate FHST-based attitude updates, procedures have been developed to redetermine the FHST scale and distortion calibration coefficients based on in-flight data gathered during normal HST operations. These scale and distortion calibrations have proven very effective operationally, and procedures are in place to monitor FHST calibration changes on a continuing basis.

Mccutcheon, R.↗

On the Gibbs phenomenon 5: Recovering exponential accuracy from collocation point values of a piecewise analytic function

The paper presents a method to recover exponential accuracy at all points (including at the discontinuities themselves), from the knowledge of an approximation to the interpolation polynomial (or trigonometrical polynomial). We show that if we are given the collocation point values (or a highly accurate approximation) at the Gauss or Gauss-Lobatto points, we can reconstruct a uniform exponentially convergent approximation to the function f(x) in any sub-interval of analyticity. The proof covers the cases of Fourier, Chebyshev, Legendre, and more general Gegenbauer collocation methods.

Gottlieb, David↗

Satellite passive-microwave observations and the distribution of net surface accumulation in the interior of Greenland

The rate of net accumulation at the surface and its distribution in ice sheet areas bounded by the dry snow line may be approximated by a determination of either firn external emissivity or various other emissivity terms. A compilation of mean annual brightness temperature for the Greenland ice sheet based on Nimbus-5 Electrically Scanning Microwave Radiometer (19 GHz) data is used to approximate the distribution of the accumulation rate in the area bounded by the dry snow line, applying a simple conversion involving a polynomial function. This is compared with a compilation produced elsewhere using Nimbus-6 SCAMS (31 GHz) data, applying a complex conversion involving the temperature dependence of the absorption coefficient, and with the latest compilation based on surface measurements of the accumulation rate. The findings indicate that in Greenland, the use of polynomial functions to convert for net surface accumulation is a reliable and simple method for estimating both the rate of accumulation and its distribution in areas where there is no melting at the surface.

Giovinetto, M. B.↗

A uniform geometrical optics and an extended uniform geometrical theory of diffraction for evaluating high frequency EM fields near smooth caustics and composite shadow boundaries

A uniform geometrical optics (UGO) and an extended uniform geometrical theory of diffraction (EUTD) are developed for evaluating high frequency electromagnetic (EM) fields within transition regions associated with a two and three dimensional smooth caustic of reflected rays and a composite shadow boundary formed by the caustic termination or the confluence of the caustic with the reflection shadow boundary (RSB). The UGO is a uniform version of the classic geometrical optics (GO). It retains the simple ray optical expressions of classic GO and employs a new set of uniform reflection coefficients. The UGO also includes a uniform version of the complex GO ray field that exists on the dark side of the smooth caustic. The EUTD is an extension of the classic uniform geometrical theory of diffraction (UTD) and accounts for the non-ray optical behavior of the UGO reflected field near caustics by using a two-variable transition function in the expressions for the edge diffraction coefficients. It also uniformly recovers the classic UTD behavior of the edge diffracted field outside the composite shadow boundary transition region. The approach employed for constructing the UGO/EUTD solution is based on a spatial domain physical optics (PO) radiation integral representation for the fields which is then reduced using uniform asymptotic procedures. The UGO/EUTD analysis is also employed to investigate the far-zone RCS problem of plane wave scattering from two and three dimensional polynomial defined surfaces, and uniform reflection, zero-curvature, and edge diffraction coefficients are derived. Numerical results for the scattering and diffraction from cubic and fourth order polynomial strips are also shown and the UGO/EUTD solution is validated by comparison to an independent moment method (MM) solution. The UGO/EUTD solution is also compared with the classic GO/UTD solution. The failure of the classic techniques near caustics and composite shadow boundaries is clearly demonstrated and it is shown that the UGO/EUTD results remain valid and uniformly reduce to the classic results away from the transition regions. Mathematical details on the asymptotic properties and efficient numerical evaluation of the canonical functions involved in the UGO/EUTD expressions are also provided.

Constantinides, E. D.↗

Splines and control theory

In this work, the relationship between splines and the control theory has been analyzed. We show that spline functions can be constructed naturally from the control theory. By establishing a framework based on control theory, we provide a simple and systematic way to construct splines. We have constructed the traditional spline functions including the polynomial splines and the classical exponential spline. We have also discovered some new spline functions such as trigonometric splines and the combination of polynomial, exponential and trigonometric splines. The method proposed in this paper is easy to implement. Some numerical experiments are performed to investigate properties of different spline approximations.

Zhang, Zhimin↗