Search NASA⌕ Search

SEARCH · Search NASA

Results for “polynomial method”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Normalization of Gravitational Acceleration Models

Unlike the uniform density spherical shell approximations of Newton, the con- sequence of spaceflight in the real universe is that gravitational fields are sensitive to the nonsphericity of their generating central bodies. The gravitational potential of a nonspherical central body is typically resolved using spherical harmonic approximations. However, attempting to directly calculate the spherical harmonic approximations results in at least two singularities which must be removed in order to generalize the method and solve for any possible orbit, including polar orbits. Three unique algorithms have been developed to eliminate these singularities by Samuel Pines [1], Bill Lear [2], and Robert Gottlieb [3]. This paper documents the methodical normalization of two1 of the three known formulations for singularity-free gravitational acceleration (namely, the Lear [2] and Gottlieb [3] algorithms) and formulates a general method for defining normalization parameters used to generate normalized Legendre Polynomials and ALFs for any algorithm. A treatment of the conventional formulation of the gravitational potential and acceleration is also provided, in addition to a brief overview of the philosophical differences between the three known singularity-free algorithms.

Eckman, Randy A.↗
NASA NTRS2016

Normalization and Implementation of Three Gravitational Acceleration Models

Unlike the uniform density spherical shell approximations of Newton, the consequence of spaceflight in the real universe is that gravitational fields are sensitive to the asphericity of their generating central bodies. The gravitational potential of an aspherical central body is typically resolved using spherical harmonic approximations. However, attempting to directly calculate the spherical harmonic approximations results in at least two singularities that must be removed to generalize the method and solve for any possible orbit, including polar orbits. Samuel Pines, Bill Lear, and Robert Gottlieb developed three unique algorithms to eliminate these singularities. This paper documents the methodical normalization of two of the three known formulations for singularity-free gravitational acceleration (namely, the Lear and Gottlieb algorithms) and formulates a general method for defining normalization parameters used to generate normalized Legendre polynomials and Associated Legendre Functions (ALFs) for any algorithm. A treatment of the conventional formulation of the gravitational potential and acceleration is also provided, in addition to a brief overview of the philosophical differences between the three known singularity-free algorithms.

Eckman, Randy A.↗
NASA NTRS1986

A comparative analysis of static and dynamic load balancing strategies

The problem of uniformly distributing the load of a parallel program over a multiprocessor system was considered. A program was analyzed whose structure permits the computation of the optimal static solution. Then four strategies for load balancing were described and their performance compared. The strategies are: (1) the optimal static assignment algorithm which is guaranteed to yield the best static solution, (2) the static binary dissection method which is very fast but suboptimal, (3) the greedy algorithm, a static fully polynomial time approximation scheme, which estimates the optimal solution to arbitrary accuracy, and (4) the predictive dynamic load balancing heuristic which uses information on the precedence relationships within the program and outperforms any of the static methods. It is also shown that the overhead incurred by the dynamic heuristic is reduced considerably if it is started off with a static assignment provided by either of the three strategies.

Iqbal, M. Ashraf↗
NASA NTRS2002

Three-Dimensional High-Order Spectral Finite Volume Method for Unstructured Grids

Many areas require a very high-order accurate numerical solution of conservation laws for complex shapes. This paper deals with the extension to three dimensions of the Spectral Finite Volume (SV) method for unstructured grids, which was developed to solve such problems. We first summarize the limitations of traditional methods such as finite-difference, and finite-volume for both structured and unstructured grids. We then describe the basic formulation of the spectral finite volume method. What distinguishes the SV method from conventional high-order finite-volume methods for unstructured triangular or tetrahedral grids is the data reconstruction. Instead of using a large stencil of neighboring cells to perform a high-order reconstruction, the stencil is constructed by partitioning each grid cell, called a spectral volume (SV), into 'structured' sub-cells, called control volumes (CVs). One can show that if all the SV cells are partitioned into polygonal or polyhedral CV sub-cells in a geometrically similar manner, the reconstructions for all the SVs become universal, irrespective of their shapes, sizes, orientations, or locations. It follows that the reconstruction is reduced to a weighted sum of unknowns involving just a few simple adds and multiplies, and those weights are universal and can be pre-determined once for all. The method is thus very efficient, accurate, and yet geometrically flexible. The most critical part of the SV method is the partitioning of the SV into CVs. In this paper we present the partitioning of a tetrahedral SV into polyhedral CVs with one free parameter for polynomial reconstructions up to degree of precision five. (Note that the order of accuracy of the method is one order higher than the reconstruction degree of precision.) The free parameter will be determined by minimizing the Lebesgue constant of the reconstruction matrix or similar criteria to obtain optimized partitions. The details of an efficient, parallelizable code to solve three-dimensional problems for any order of accuracy are then presented. Important aspects of the data structure are discussed. Comparisons with the Discontinuous Galerkin (DG) method are made. Numerical examples for wave propagation problems are presented.

Liu, Yen↗
NASA NTRS1989

Composite structure dynamics and space applicability

Guided wave modes in cladded or uncladded fiber-reinforced composite plates and tubes have been analyzed using a stiffness method in which the displacement variation through the thickness is approximated by polynomial interpolation functions. This allows for an arbitrary number of laminations and fiber orientations different from lamina to lamina. It is shown that dispersive behavior of guided modes depends significantly on the cladding, number of laminae, and interfaces between the adjacent laminae. A hybrid modeling technique is described in which an inner region containing cracks (or other defects) is discretized by finite elements, and the field in the exterior region is represented in terms of modes that are found using the stiffness method described above. It is found that the reflected and transmitted amplitudes of modes vary significantly with the size of a transverse or longitudinal (delamination) crack and frequency. We have also studied the impact response of a unidirectional fiber-reinforced plate. Received signals at the epicentral and other locations are shown. Strong longitudinal anisotropy of the graphite/epoxy plate causes the signal to be considerably different from that in an isotropic plate.

Datta, Subhendu K.↗
NASA NTRS1995

Statistically generated weighted curve fit of residual functions for modal analysis of structures

A statistically generated weighting function for a second-order polynomial curve fit of residual functions has been developed. The residual flexibility test method, from which a residual function is generated, is a procedure for modal testing large structures in an external constraint-free environment to measure the effects of higher order modes and interface stiffness. This test method is applicable to structures with distinct degree-of-freedom interfaces to other system components. A theoretical residual function in the displacement/force domain has the characteristics of a relatively flat line in the lower frequencies and a slight upward curvature in the higher frequency range. In the test residual function, the above-mentioned characteristics can be seen in the data, but due to the present limitations in the modal parameter evaluation (natural frequencies and mode shapes) of test data, the residual function has regions of ragged data. A second order polynomial curve fit is required to obtain the residual flexibility term. A weighting function of the data is generated by examining the variances between neighboring data points. From a weighted second-order polynomial curve fit, an accurate residual flexibility value can be obtained. The residual flexibility value and free-free modes from testing are used to improve a mathematical model of the structure. The residual flexibility modal test method is applied to a straight beam with a trunnion appendage and a space shuttle payload pallet simulator.

Bookout, P. S.↗
NASA NTRS2023

Towards Soil Moisture Profile Estimation in the Root Zone Using L- and P-Band Radiometer Observations: A Coherent Modelling Approach

Precision irrigation management and crop water stress assessment rely on accurate estimation of root zone soil moisture. However, only the top 5cm soil moisture can be estimated using the two current passive microwave satellite missions, Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP), which operate at L-band (wavelength of ~21cm). Since the contributing depth of the soil to brightness temperature increases with observation wavelength, it is expected that a P-band (wavelength of ~40cm) radiometer could potentially provide soil moisture information from deeper layers of the soil profile. Moreover, by combining both L- and P- bands, it is hypothesized that the soil moisture profile can be estimated even beyond their individual observation depths. The aim of this study was to demonstrate the potential of combined L-band and P-band radiometer observations to estimate the soil moisture profile under flat bare soil using a stratified coherent forward model. Brightness temperature observations at L-band and P-band from a tower based experimental site across a dry (April 2019) and a wet (March 2020) period, covering different soil moisture profile shapes, were used in this study. Results from an initial synthetic study showed that the performance of a combined L-band and P-band approach was better than the performance of using either band individually, with an average depth over which reliable soil moisture profile information could be estimated (i.e. with a target root mean square error (RMSE) of less than 0.04 m3/m3) being 20cm for linear and 15cm for second-order polynomial functions. Other functions were also tested but found to have a poorer performance. Applying the method to the tower-based brightness temperature achieved an average estimation depth of 28cm (20cm) and 5cm (5cm) during the dry and wet periods respectively when using a second-order polynomial (linear) function. These findings highlight the opportunity of a satellite mission with L-band and P-band observations to accurately estimate the soil moisture profile to as deep as 30cm globally.

Foad Brakhasi↗
NASA NTRS1993

Microwave Scattering Model for Grass Blade Structures

In this paper, the electromagnetic scattering solution for a grass blade with complex cross-section geometry is considered. It is assumed that the blade cross section is electrically small, but its length is large compared to the incident wavelength. In a recent study it has been shown that the scattering solution for such problems, in the form of a polarizability tensor, can be obtained using the low-frequency approximation in conjunction with the method of moments. In addition, the study shows that the relationship between the polarizability tensor of a dielectric cylinder and its dielectric constant can be approximated by a simple algebraic expression. The results of this study are used to show that this algebraic approximation is valid also for cylinders with cross sections the shape of grass blades, providing that proper values am selected for each of three constants appearing in the expression. These constants are dependent on cylinder shape, and if the relationship between the constants and the three parameters describing a grass blade shape can be determined, an algebraic approximation relating polarizability tensor to blade shape, as well as dielectric constant, can be formed. Since the elements of the polarizability tensor are dependent on only these parameters, this algebraic approximation can replace the cumbersome method of moments model. A conjugate gradient method is then implemented to correctly determine the three constants of the algebraic approximation for each blade shape. A third-order polynomial fit to the data is then determined for each constant, thus providing a complete analytic replacement to the numerical (moment method) scattering model. Comparisons of this approximation to the numerical model show an average error of less than 3%.

Stiles, James M.↗
NASA NTRS2003

Analytic Method for Computing Instrument Pointing Jitter

A new method of calculating the root-mean-square (rms) pointing jitter of a scientific instrument (e.g., a camera, radar antenna, or telescope) is introduced based on a state-space concept. In comparison with the prior method of calculating the rms pointing jitter, the present method involves significantly less computation. The rms pointing jitter of an instrument (the square root of the jitter variance shown in the figure) is an important physical quantity which impacts the design of the instrument, its actuators, controls, sensory components, and sensor- output-sampling circuitry. Using the Sirlin, San Martin, and Lucke definition of pointing jitter, the prior method of computing the rms pointing jitter involves a frequency-domain integral of a rational polynomial multiplied by a transcendental weighting function, necessitating the use of numerical-integration techniques. In practice, numerical integration complicates the problem of calculating the rms pointing error. In contrast, the state-space method provides exact analytic expressions that can be evaluated without numerical integration.

Bayard, David↗
NASA NTRS2007

Analytic Guidance for the First Entry in a Skip Atmospheric Entry

This paper presents an analytic method to generate a reference drag trajectory for the first entry portion of a skip atmospheric entry. The drag reference, expressed as a polynomial function of the velocity, will meet the conditions necessary to fit the requirements of the complete entry phase. The generic method proposed to generate the drag reference profile is further simplified by thinking of the drag and the velocity as density and cumulative distribution functions respectively. With this notion it will be shown that the reference drag profile can be obtained by solving a linear algebraic system of equations. The resulting drag profile is flown using the feedback linearization method of differential geometric control as guidance law with the error dynamics of a second order homogeneous equation in the form of a damped oscillator. This approach was first proposed as a revisited version of the Space Shuttle Orbiter entry guidance. However, this paper will show that it can be used to fly the first entry in a skip entry trajectory. In doing so, the gains in the error dynamics will be changed at a certain point along the trajectory to improve the tracking performance.

Garcia-Llama, Eduardo↗
NASA NTRS1973

Geometric aspects in digital analysis of Multi-Spectral Scanner (MSS) data

Present automated systems of interpretation which apply pattern recognition techniques on MSS data do not fully consider the geometry of the acquisition system. In an effort to improve the usefulness of the MSS data when digitally treated, geometric aspects are analyzed and discussed. Attempts to correct for scanner instabilities in position and orientation by affine and polynomial transformations, as well as by modified collinearity equations are described. Methods of accounting for panoramic and relief effects are also discussed. It is anticipated that reliable area as well as position determinations can be accomplished during the process of automatic interpretation. A concept for a unified approach to the treatment of remote sensing data, both metric and nonmetric is presented.

Mikhail, E. M.↗
NASA NTRS1976

A method for the selection of a functional form for a thermodynamic equation of state using weighted linear least squares stepwise regression

A method was developed for establishing a rational choice of the terms to be included in an equation of state with a large number of adjustable coefficients. The methods presented were developed for use in the determination of an equation of state for oxygen and nitrogen. However, a general application of the methods is possible in studies involving the determination of an optimum polynomial equation for fitting a large number of data points. The data considered in the least squares problem are experimental thermodynamic pressure-density-temperature data. Attention is given to a description of stepwise multiple regression and the use of stepwise regression in the determination of an equation of state for oxygen and nitrogen.

Jacobsen, R. T.↗
NASA NTRS1983

Convergence and stability properties of minimal polynomial and reduced rank extrapolation algorithms

The minimal polynomial and reduced rank extrapolation algorithms are two acceleration of convergence methods for sequences of vectors. In a recent survey these methods were tested and compared with the scalar, vector, topological epsilon algorithms, and were observed to be more efficient than the latter. It was also observed that the two methods have similar convergence properties. The convergence and stability properties of these methods are analyzed and the performance of the acceleration methods when applied to a class of vector sequences that includes those sequences obtained from systems of linear equations by using matrix iterative methods is discussed.

Sidi, A.↗
NASA NTRS1982

On the geodetic applications of simultaneous range-differencing to LAGEOS

The possibility of improving the accuracy of geodetic results by use of simultaneously observed ranges to Lageos, in a differencing mode, from pairs of stations was studied. Simulation tests show that model errors can be effectively minimized by simultaneous range differencing (SRD) for a rather broad class of network satellite pass configurations. The methods of least squares approximation are compared with monomials and Chebyshev polynomials and the cubic spline interpolation. Analysis of three types of orbital biases (radial, along- and across track) shows that radial biases are the ones most efficiently minimized in the SRC mode. The degree to which the other two can be minimized depends on the type of parameters under estimation and the geometry of the problem. Sensitivity analyses of the SRD observation show that for baseline length estimations the most useful data are those collected in a direction parallel to the baseline and at a low elevation. Estimating individual baseline lengths with respect to an assumed but fixed orbit not only decreases the cost, but it further reduces the effects of model biases on the results as opposed to a network solution. Analogous results and conclusions are obtained for the estimates of the coordinates of the pole.

Pablis, E. C.↗
NASA NTRS1983

On the three-dimensional vibrations of the cantilevered rectangular parallelepiped

A solution is presented for the three-dimensional problem of determining the free vibration frequencies and mode shapes for a rectangular parallelepiped which is completely fixed on one face and free on the other five faces. This problem apparently is previously unsolved in the published literature. The Ritz method is used, with displacements assumed in the form of algebraic polynomials. Convergence is studied. Numerical results are given for the first five frequencies of each of the four symmetry classes of vibration, for five thick parallelepiped configurations, including the cube. Contour plots are exhibited for the modal displacements of the cube. The effects of varying Poisson's ratio are also observed.

Leissa, A.↗
NASA NTRS1990

Interpolation schemes in the control of systems with unknown dynamics

Time-delay control for systems with unknown dynamics involves estimations. In its present form, these estimations include not only the interpolation of a time-varying function but also its derivative between data points. Presently, the control law is reformulated, taking into account the computation delay. A window-shifting scheme is then devised to view n data points as they are sampled in order to perform the function estimation. At each new sampling time, the window shifts to include the new data point as well as the last n-1 points. Several interpolation methods are considered. These methods use data points from a sampled function to construct a polynomial estimate of the function. The control system performance was experimentally tested using a servosystem. The results show that a Newtonian interpolation provided best results when the computation time was 0.4 times the sampling period.

Youcef-Toumi, K.↗
NASA NTRS1995

Improved L-C resonant decay technique for Q measurement of quasilinear power inductors: New results for MPP and ferrite powdered cores

The L-C resonant decay technique for measuring circuit Q or losses is improved by eliminating the switch from the inductor-capacitor loop. A MOSFET switch is used instead to momentarily connect the resonant circuit to an existing voltage source, which itself is gated off during the decay transient. Very reproducible, low duty cycle data could be taken this way over a dynamic voltage range of at least 10:1. Circuit Q is computed from a polynomial fit to the sequence of the decaying voltage maxima. This method was applied to measure the losses at 60 kHz in inductors having loose powder cores of moly permalloy and an Mn-Zn power ferrite. After the copper and capacitor losses are separated out, the resulting specific core loss is shown to be roughly as expected for the MPP powder, but anomalously high for the ferrite powder. Possible causes are mentioned.

Niedra, Janis M.↗