Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Recursive algorithms for vector extrapolation methods

Three classes of recursion relations are devised for implementing some extrapolation methods for vector sequences. One class of recursion relations can be used to implement methods like the modified minimal polynomial extrapolation and the topological epsilon algorithm; another allows implementation of methods like minimal polynomial and reduced rank extrapolation; while the remaining class can be employed in the implementation of the vector E-algorithm. Operation counts and storage requirements for these methods are also discussed, and some related techniques for special applications are also presented. Included are methods for the rapid evaluations of the vector E-algorithm.

Ford, William F.↗
NASA NTRS2014

Validation of Spherically Symmetric Inversion by Use of a Tomographically Reconstructed Three-Dimensional Electron Density of the Solar Corona

Determining the coronal electron density by the inversion of white-light polarized brightness (pB) measurements by coronagraphs is a classic problem in solar physics. An inversion technique based on the spherically symmetric geometry (spherically symmetric inversion, SSI) was developed in the 1950s and has been widely applied to interpret various observations. However, to date there is no study of the uncertainty estimation of this method. We here present the detailed assessment of this method using a three-dimensional (3D) electron density in the corona from 1.5 to 4 solar radius as a model, which is reconstructed by a tomography method from STEREO/COR1 observations during the solar minimum in February 2008 (Carrington Rotation, CR 2066).We first show in theory and observation that the spherically symmetric polynomial approximation (SSPA) method and the Van de Hulst inversion technique are equivalent. Then we assess the SSPA method using synthesized pB images from the 3D density model, and find that the SSPA density values are close to the model inputs for the streamer core near the plane of the sky (POS) with differences generally smaller than about a factor of two; the former has the lower peak but extends more in both longitudinal and latitudinal directions than the latter. We estimate that the SSPA method may resolve the coronal density structure near the POS with angular resolution in longitude of about 50 deg. Our results confirm the suggestion that the SSI method is applicable to the solar minimum streamer (belt), as stated in some previous studies. In addition, we demonstrate that the SSPA method can be used to reconstruct the 3D coronal density, roughly in agreement with the reconstruction by tomography for a period of low solar activity (CR 2066). We suggest that the SSI method is complementary to the 3D tomographic technique in some cases, given that the development of the latter is still an ongoing research effort.

corona↗
NASA NTRS1989

Closure of the cubic tensor polynomial failure surface

An analytical method has been developed to ensure closure of the cubic form of the tensor polynomial strength criterion. The intrinsic complexity of the cubic function is such that special conditions must be met to close the failure surface in three-dimensional stress space. These requirements are derived in terms of non-intersecting conditions for asymptotes and an asymptotic plane. To demonstrate the validity of this approach, closed failure surfaces were derived for two graphite/epoxy material systems (3M SP288-T300 and IM7 8551-7). The agreement of test data with this model clearly shows that it is possible to use a higher order cubic failure theory with confidence.

Jiang, Zhiqing↗
NASA NTRS1994

Separation of lithospheric, external, and core components of the south polar geomagnetic field at satellite altitudes

We present a new approach to producing scalar Magsat magnetic anomaly maps based on correlation coefficient filtering and the use of almost all of the available orbits. Our method differs from earlier techniques with respect to the following: (1) Passes are selected based on their variance properties rather than planetary indices such as Kp. (2) The core field model is least squares fit to individual passes and subsequently removed instead of substracting the model directly. This technique replaces band pass filtering or polynomial trend removal methods. (3) Each selected pass is sorted geographically and by local time, placed into one of four different altitude bands, and correlation coefficient filtered with the two adjacent passes. The filtering is the second step toward isolating the static lithospheric signal from the more dynamic external field signals. (4) Least squares collocation is used to grid the correlated passes; subsequently, the dawn and dusk maps are also correlation filtered providing another step toward removal of external fields. (5) The four resultant total field maps are continued to a common altitude and again correlation filtered for the static lithospheric anomalies. (6) The filtered results are then averaged together to provide a new total field map of the lithosphere south of 40 deg S latitude. Our total field map differs from previous efforts over the crustal blocks of West Antarctica. We obtained a positive anomaly over Edward VII Peninsula, extending into the Byrd subglacial basin and obtained a negative anomaly over the Ellsworth Mountains and parts of the Byrd subglacial basin. Also, a positive anomaly extending from the Ross Sea to offshore Wilkes Land is present in our map; however, this feature is absent in other maps. Positive anomalies marking the Weddell Sea in previous efforts are not present in our map. Prominent external field anomalies in the quadrant offshore of Wilkes Land are present in all previous efforts; however, these signatures are reduced in our total field map. This map significantly advances efforts for the identification and interpretation of magnetic anomalies of the south polar lithosphere.

Alsdorf, Douglas E.↗
NASA NTRS1992

A new model for estimating total body water from bioelectrical resistance

Estimation of total body water (T) from bioelectrical resistance (R) is commonly done by stepwise regression models with height squared over R, H(exp 2)/R, age, sex, and weight (W). Polynomials of H(exp 2)/R have not been included in these models. We examined the validity of a model with third order polynomials and W. Methods: T was measured with oxygen-18 labled water in 27 subjects. R at 50 kHz was obtained from electrodes placed on the hand and foot while subjects were in the supine position. A stepwise regression equation was developed with 13 subjects (age 31.5 plus or minus 6.2 years, T 38.2 plus or minus 6.6 L, W 65.2 plus or minus 12.0 kg). Correlations, standard error of estimates and mean differences were computed between T and estimated T's from the new (N) model and other models. Evaluations were completed with the remaining 14 subjects (age 32.4 plus or minus 6.3 years, T 40.3 plus or minus 8 L, W 70.2 plus or minus 12.3 kg) and two of its subgroups (high and low) Results: A regression equation was developed from the model. The only significant mean difference was between T and one of the earlier models. Conclusion: Third order polynomials in regression models may increase the accuracy of estimating total body water. Evaluating the model with a larger population is needed.

Siconolfi, S. F.↗
NASA NTRS1981

Vibrations of cantilevered shallow cylindrical shells of rectangular planform

A cantilevered, shallow shell of circular cylindrical curvature and rectangular planform exhibits free vibration behavior which differs considerably from that of a cantilevered beam or of a flat plate. Some numerical results can be found for the problem in the previously published literature, mainly obtained by using various finite element methods. The present paper is the first definitive study of the problem, presenting accurate non-dimensional frequency parameters for wide ranges of aspect ratio, shallowness ratio and thickness ratio. The analysis is based upon shallow shell theory. Numerical results are obtained by using the Ritz method, with algebraic polynomial trial functions for the displacements. Convergence is investigated, with attention being given both to the number of terms taken for each co-ordinate direction and for each of the three components of displacement. Accuracy of the results is also established by comparison with finite element results for shallow shells and with other accurate flat plate solutions.

Leissa, A. W.↗
NASA NTRS1983

Vibrations of cantilevered doubly-curved shallow shells

Vibrational characteristics are determined for a previously unsolved class of problems, that of doubly-curved shallow shells having rectangular planforms, clamped along one edge and free on the other three. The solution procedure uses the Ritz method with algebraic polynomial trial functions. Convergence studies are made, and accurate frequencies and contour plots of mode shapes are presented for various curvature ratios, including spherical, circular cylindrical and hyperbolic paraboloidal shells. Particular emphasis is given to the effect of adding spanwise curvature to shells having chordwise curvature; numerous published references already exist for the case of zero spanwise curvature. The effects of changing aspect ratio, thickness ratio and Poisson's ratio are also studied.

Leissa, A. W.↗
NASA NTRS1984

Slave finite elements: The temporal element approach to nonlinear analysis

A formulation method for finite elements in space and time incorporating nonlinear geometric and material behavior is presented. The method uses interpolation polynomials for approximating the behavior of various quantities over the element domain, and only explicit integration over space and time. While applications are general, the plate and shell elements that are currently being programmed are appropriate to model turbine blades, vanes, and combustor liners.

Gellin, S.↗
NASA NTRS1991

Investigation and Control of Rotordynamic Instability in Typical Large Turbogenerators

Described here are the investigation and results of recent studies to solve oil whip in typical large turbogenerators. Included are calculations of the instability speeds and system damping of rotor-bearing systems. The polynomial-transfer matrix method which was developed by the author during the last few years is used in the calculations. Vibration measurements and data indicate the stability of these units. Research indicates that the cause of the instability lies in the three-bole offset bearings. Work was done to solve these problems, and industry tests were performed on one of these abnormal systems.

Lu, Songyuan↗
NASA NTRS2001

Multidisciplinary Design Optimization of a Full Vehicle with High Performance Computing

Multidisciplinary design optimization (MDO) of a full vehicle under the constraints of crashworthiness, NVH (Noise, Vibration and Harshness), durability, and other performance attributes is one of the imperative goals for automotive industry. However, it is often infeasible due to the lack of computational resources, robust simulation capabilities, and efficient optimization methodologies. This paper intends to move closer towards that goal by using parallel computers for the intensive computation and combining different approximations for dissimilar analyses in the MDO process. The MDO process presented in this paper is an extension of the previous work reported by Sobieski et al. In addition to the roof crush, two full vehicle crash modes are added: full frontal impact and 50% frontal offset crash. Instead of using an adaptive polynomial response surface method, this paper employs a DOE/RSM method for exploring the design space and constructing highly nonlinear crash functions. Two NMO strategies are used and results are compared. This paper demonstrates that with high performance computing, a conventionally intractable real world full vehicle multidisciplinary optimization problem considering all performance attributes with large number of design variables become feasible.

Yang, R. J.↗
NASA NTRS1998

Conversion from Engineering Units to Telemetry Counts on Dryden Flight Simulators

Dryden real-time flight simulators encompass the simulation of pulse code modulation (PCM) telemetry signals. This paper presents a new method whereby the calibration polynomial (from first to sixth order), representing the conversion from counts to engineering units (EU), is numerically inverted in real time. The result is less than one-count error for valid EU inputs. The Newton-Raphson method is used to numerically invert the polynomial. A reverse linear interpolation between the EU limits is used to obtain an initial value for the desired telemetry count. The method presented here is not new. What is new is how classical numerical techniques are optimized to take advantage of modem computer power to perform the desired calculations in real time. This technique makes the method simple to understand and implement. There are no interpolation tables to store in memory as in traditional methods. The NASA F-15 simulation converts and transmits over 1000 parameters at 80 times/sec. This paper presents algorithm development, FORTRAN code, and performance results.

Fantini, Jay A.↗
NASA NTRS1992

The algebraic decoding of the (41, 21, 9) quadratic residue code

A new algebraic approach for decoding the quadratic residue (QR) codes, in particular the (41, 21, 9) QR code is presented. The key ideas behind this decoding technique are a systematic application of the Sylvester resultant method to the Newton identities associated with the code syndromes to find the error-locator polynomial, and next a method for determining error locations by solving certain quadratic, cubic and quartic equations over GF(2 exp m) in a new way which uses Zech's logarithms for the arithmetic. The algorithms developed here are suitable for implementation in a programmable microprocessor or special-purpose VLSI chip. It is expected that the algebraic methods developed here can apply generally to other codes such as the BCH and Reed-Solomon codes.

Reed, Irving S.↗
NASA NTRS1990

Implicit application of polynomial filters in a k-step Arnoldi method

The Arnoldi process is a well known technique for approximating a few eigenvalues and corresponding eigenvectors of a general square matrix. Numerical difficulties such as loss of orthogonality and assessment of the numerical quality of the approximations as well as a potential for unbounded growth in storage have limited the applicability of the method. These issues are addressed by fixing the number of steps in the Arnoldi process at a prescribed value k and then treating the residual vector as a function of the initial Arnoldi vector. This starting vector is then updated through an iterative scheme that is designed to force convergence of the residual to zero. The iterative scheme is shown to be a truncation of the standard implicitly shifted QR-iteration for dense problems and it avoids the need to explicitly restart the Arnoldi sequence. The main emphasis of this paper is on the derivation and analysis of this scheme. However, there are obvious ways to exploit parallelism through the matrix-vector operations that comprise the majority of the work in the algorithm. Preliminary computational results are given for a few problems on some parallel and vector computers.

Sorensen, D. C.↗
NASA NTRS1998

On the Numerical Formulation of Parametric Linear Fractional Transformation (LFT) Uncertainty Models for Multivariate Matrix Polynomial Problems

Robust control system analysis and design is based on an uncertainty description, called a linear fractional transformation (LFT), which separates the uncertain (or varying) part of the system from the nominal system. These models are also useful in the design of gain-scheduled control systems based on Linear Parameter Varying (LPV) methods. Low-order LFT models are difficult to form for problems involving nonlinear parameter variations. This paper presents a numerical computational method for constructing and LFT model for a given LPV model. The method is developed for multivariate polynomial problems, and uses simple matrix computations to obtain an exact low-order LFT representation of the given LPV system without the use of model reduction. Although the method is developed for multivariate polynomial problems, multivariate rational problems can also be solved using this method by reformulating the rational problem into a polynomial form.

Belcastro, Christine M.↗
NASA NTRS1986

Full-envelope aerodynamic modeling of the Harrier aircraft

A project to identify a full-envelope model of the YAV-8B Harrier using flight-test and parameter identification techniques is described. As part of the research in advanced control and display concepts for V/STOL aircraft, a full-envelope aerodynamic model of the Harrier is identified, using mathematical model structures and parameter identification methods. A global-polynomial model structure is also used as a basis for the identification of the YAV-8B aerodynamic model. State estimation methods are used to ensure flight data consistency prior to parameter identification.Equation-error methods are used to identify model parameters. A fixed-base simulator is used extensively to develop flight test procedures and to validate parameter identification software. Using simple flight maneuvers, a simulated data set was created covering the YAV-8B flight envelope from about 0.3 to 0.7 Mach and about -5 to 15 deg angle of attack. A singular value decomposition implementation of the equation-error approach produced good parameter estimates based on this simulated data set.

Mcnally, B. David↗
NASA NTRS1974

Method of Characteristics Calculations and Computer Code for Materials with Arbitrary Equations of State and Using Orthogonal Polynomial Least Square Surface Fits

A numerical scheme using the method of characteristics to calculate the flow properties and pressures behind decaying shock waves for materials under hypervelocity impact is developed. Time-consuming double interpolation subroutines are replaced by a technique based on orthogonal polynomial least square surface fits. Typical calculated results are given and compared with the double interpolation results. The complete computer program is included.

Chang, T. S.↗