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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS2007

On Flow Stagnation in a Tube Radiator

An analysis of the physical process for occurrence of flow stagnation in a space tube-radiator is performed and the mechanism and mathematic description for the flow stagnation are presented. Two causes for pressure drop unbalance between tubes of the radiator are identified: non-uniform cooling environment and different local flow resistances between the tubes. This analysis provides a theoretical basis for experimental simulations of the flow stagnation in a ground-based lab as well as two suggested methods to experimentally simulate flow stagnation. Criteria for the flow stagnation, depending on the viscosity data regressive polynomial, are derived from the extreme condition of the pressure drop in colder tubes. A preliminary numerical calculation is conducted for a space tube-radiator model which confirms the physical and mathematical analyses. The prediction by the criteria for flow stagnation in the tube-radiator model coincides with the numerical calculation result.

Motil, Brian↗

NASA NTRS2010

Model Robust Calibration: Method and Application to Electronically-Scanned Pressure Transducers

This article presents the application of a recently developed statistical regression method to the controlled instrument calibration problem. The statistical method of Model Robust Regression (MRR), developed by Mays, Birch, and Starnes, is shown to improve instrument calibration by reducing the reliance of the calibration on a predetermined parametric (e.g. polynomial, exponential, logarithmic) model. This is accomplished by allowing fits from the predetermined parametric model to be augmented by a certain portion of a fit to the residuals from the initial regression using a nonparametric (locally parametric) regression technique. The method is demonstrated for the absolute scale calibration of silicon-based pressure transducers.

Walker, Eric L.↗

NASA NTRS2012

Method for Pre-Conditioning a Measured Surface Height Map for Model Validation

This software allows one to up-sample or down-sample a measured surface map for model validation, not only without introducing any re-sampling errors, but also eliminating the existing measurement noise and measurement errors. Because the re-sampling of a surface map is accomplished based on the analytical expressions of Zernike-polynomials and a power spectral density model, such re-sampling does not introduce any aliasing and interpolation errors as is done by the conventional interpolation and FFT-based (fast-Fourier-transform-based) spatial-filtering method. Also, this new method automatically eliminates the measurement noise and other measurement errors such as artificial discontinuity. The developmental cycle of an optical system, such as a space telescope, includes, but is not limited to, the following two steps: (1) deriving requirements or specs on the optical quality of individual optics before they are fabricated through optical modeling and simulations, and (2) validating the optical model using the measured surface height maps after all optics are fabricated. There are a number of computational issues related to model validation, one of which is the "pre-conditioning" or pre-processing of the measured surface maps before using them in a model validation software tool. This software addresses the following issues: (1) up- or down-sampling a measured surface map to match it with the gridded data format of a model validation tool, and (2) eliminating the surface measurement noise or measurement errors such that the resulted surface height map is continuous or smoothly-varying. So far, the preferred method used for re-sampling a surface map is two-dimensional interpolation. The main problem of this method is that the same pixel can take different values when the method of interpolation is changed among the different methods such as the "nearest," "linear," "cubic," and "spline" fitting in Matlab. The conventional, FFT-based spatial filtering method used to eliminate the surface measurement noise or measurement errors can also suffer from aliasing effects. During re-sampling of a surface map, this software preserves the low spatial-frequency characteristic of a given surface map through the use of Zernike-polynomial fit coefficients, and maintains mid- and high-spatial-frequency characteristics of the given surface map by the use of a PSD model derived from the two-dimensional PSD data of the mid- and high-spatial-frequency components of the original surface map. Because this new method creates the new surface map in the desired sampling format from analytical expressions only, it does not encounter any aliasing effects and does not cause any discontinuity in the resultant surface map.

Sidick, Erkin↗

NASA NTRS1987

Lewis inverse design code (LINDES): Users manual

The method of complex characteristics and hodograph transformation for the design of shockless airfoils was introduced by Bauer, Garabedian, and Korn and has been extended by the author to design subcritical and supercritical cascades with high solidities and large inlet angles. This new capability was achieved by introducing a new conformal mapping of the hodograph domain onto an ellipse and expanding the solution in terms of Chebyshev polynomials. A new computer code, the NASA Lewis inverse design code, was developed based on this idea. This new design code is an efficient method for the design of airfoils in cascade. In particular, the design of subcritical cascades of airfoils is a very fast, robust, and versatile process. The inverse design code can be made to interact with a turbulent boundary layer calculation to obtain airfoils with no separated flows at the design condition. This report is intended to serve as a users manual for this design code. Material previously reported by the author is included here for completeness and quick access to the user. The manual contains a description of the method followed by a discussion of the design procedure and examples. The input parameters necessary to run the code are then described and their default values given. Output listings corresponding to six different blade shapes designed with the code are given, as well as the necessary input data to reproduce the computer runs. The examples have been chosen to show that a wide range of applications can be covered with the code, ranging from supercritical propeller sections to wind tunnel turning vanes that can operate with a large inlet flow angle range.

Sanz, Jose M.↗

NASA NTRS1991

A Lagrange multiplier based divide and conquer finite element algorithm

A novel domain decomposition method based on a hybrid variational principle is presented. Prior to any computation, a given finite element mesh is torn into a set of totally disconnected submeshes. First, an incomplete solution is computed in each subdomain. Next, the compatibility of the displacement field at the interface nodes is enforced via discrete, polynomial and/or piecewise polynomial Lagrange multipliers. In the static case, each floating subdomain induces a local singularity that is resolved very efficiently. The interface problem associated with this domain decomposition method is, in general, indefinite and of variable size. A dedicated conjugate projected gradient algorithm is developed for solving the latter problem when it is not feasible to explicitly assemble the interface operator. When implemented on local memory multiprocessors, the proposed methodology requires less interprocessor communication than the classical method of substructuring. It is also suitable for parallel/vector computers with shared memory and compares favorably with factorization based parallel direct methods.

Farhat, C.↗

NASA NTRS1996

Element Library for Three-Dimensional Stress Analysis by the Integrated Force Method

The Integrated Force Method, a recently developed method for analyzing structures, is extended in this paper to three-dimensional structural analysis. First, a general formulation is developed to generate the stress interpolation matrix in terms of complete polynomials of the required order. The formulation is based on definitions of the stress tensor components in term of stress functions. The stress functions are written as complete polynomials and substituted into expressions for stress components. Then elimination of the dependent coefficients leaves the stress components expressed as complete polynomials whose coefficients are defined as generalized independent forces. Such derived components of the stress tensor identically satisfy homogenous Navier equations of equilibrium. The resulting element matrices are invariant with respect to coordinate transformation and are free of spurious zero-energy modes. The formulation provides a rational way to calculate the exact number of independent forces necessary to arrive at an approximation of the required order for complete polynomials. The influence of reducing the number of independent forces on the accuracy of the response is also analyzed. The stress fields derived are used to develop a comprehensive finite element library for three-dimensional structural analysis by the Integrated Force Method. Both tetrahedral- and hexahedral-shaped elements capable of modeling arbitrary geometric configurations are developed. A number of examples with known analytical solutions are solved by using the developments presented herein. The results are in good agreement with the analytical solutions. The responses obtained with the Integrated Force Method are also compared with those generated by the standard displacement method. In most cases, the performance of the Integrated Force Method is better overall.

Kaljevic, Igor↗

DOE OSTI2025

Customizable wave tailoring nonlinear materials enabled by bilevel inverse design

Abstract Passive wave transformation via nonlinearity is ubiquitous in settings from acoustics to optics and electromagnetics. It is well known that different nonlinearities yield different effects on propagating signals, which raises the question of “what precise nonlinearity is the best for a given wave tailoring application?” In this work, considering a one-dimensional spring-mass chain connected by polynomial springs (a variant of the Fermi-Pasta-Ulam-Tsingou system), we introduce a bilevel inverse design method which couples the shape optimization of structures for tailored constitutive responses with reduced-order nonlinear dynamical inverse design. We apply it to two qualitatively distinct problems—minimization of peak transmitted kinetic energy from impact, and pulse shape transformation—demonstrating our method’s breadth of applicability. For the impact problem, we obtain two fundamental insights. First, small differences in nonlinearity can drastically change the dynamic response of the system, from severely under- to outperforming a comparative linear system. Second, the oft-used strategy of impact mitigation via “energy locking” bistability can be significantly outperformed by our optimal nonlinearity. We validate this case with impact experiments and find excellent agreement. This study establishes a framework for broader passive nonlinear mechanical wave tailoring material design, with applications to computing, signal processing, shock mitigation, and autonomous materials.

Science & Technology - Other Topics↗

DOE OSTI2025

Reassessment of Boiling Water Reactor Rod Cusping Capability of VERA-MPACT

This technical report assesses the current modeling capabilities of the Virtual Environment for Reactor Applications (VERA) for minimizing control rod cusping in boiling water reactors. Assessment is performed on a 10 × 10 General Electric (GE)-14 assembly as well as a simplified version of the same assembly. Runs were performed using VERA for three available cusping treatments in the Michigan PArallel Characteristics Transport (MPACT) code: no treatment, subplane treatment, and fitted polynomial correction treatment. Subplane treatment is discussed, but specific results are not provided because of failures in the method of implementation for boiling water reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Reassessment of Boiling Water Reactor Rod Cusping Capability of VERA-MPACT

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1982

Tolerance analysis of optical telescopes using coherent addition of wavefront errors

A near diffraction-limited telescope requires that tolerance analysis be done on the basis of system wavefront error. One method of analyzing the wavefront error is to represent the wavefront error function in terms of its Zernike polynomial expansion. A Ramsey-Korsch ray trace package, a computer program that simulates the tracing of rays through an optical telescope system, was expanded to include the Zernike polynomial expansion up through the fifth-order spherical term. An option to determine a 3 dimensional plot of the wavefront error function was also included in the Ramsey-Korsch package. Several assimulation runs were analyzed to determine the particular set of coefficients in the Zernike expansion that are effected by various errors such as tilt, decenter and despace. A 3 dimensional plot of each error up through the fifth-order spherical term was also included in the study. Tolerance analysis data are presented.

Davenport, J. W.↗

NASA NTRS1992

Eleven-Species Thermochemical Model Of Air

SPECIES program computes thermodynamic and transport properties and equilibrium constants. Uses least-squares polynomial curve fits of most accurate data available to provide requested values more quickly than table-lookup methods. Subroutines computing transport coefficients and collision cross sections include additional code to correct for any electron pressure when working with ionic species. Species represented: diatomic nitrogen, diatomic oxygen, atomic nitrogen, atomic oxygen, nitric oxide, ionized nitric oxide, free electron, ionized atomic nitrogen, ionized atomic oxygen, ionized diatomic nitrogen, and ionized diatomic oxygen. Written in FORTRAN 77.

Thompson, Richard A.↗

NASA NTRS2001

Distribution functions of probabilistic automata

Each probabilistic automaton M over an alphabet A defines a probability measure Prob sub(M) on the set of all finite and infinite words over A. We can identify a k letter alphabet A with the set {0, 1,..., k-1}, and, hence, we can consider every finite or infinite word w over A as a radix k expansion of a real number X(w) in the interval [0, 1]. This makes X(w) a random variable and the distribution function of M is defined as usual: F(x) := Prob sub(M) { w: X(w) < x }. Utilizing the fixed-point semantics (denotational semantics), extended to probabilistic computations, we investigate the distribution functions of probabilistic automata in detail. Automata with continuous distribution functions are characterized. By a new, and much more easier method, it is shown that the distribution function F(x) is an analytic function if it is a polynomial. Finally, answering a question posed by D. Knuth and A. Yao, we show that a polynomial distribution function F(x) on [0, 1] can be generated by a prob abilistic automaton iff all the roots of F'(x) = 0 in this interval, if any, are rational numbers. For this, we define two dynamical systems on the set of polynomial distributions and study attracting fixed points of random composition of these two systems.

Denotational semantics↗

NASA NTRS2013

Bayesian Revision of Residual Detection Power

This paper addresses some issues with quality assessment and quality assurance in response surface modeling experiments executed in wind tunnels. The role of data volume on quality assurance for response surface models is reviewed. Specific wind tunnel response surface modeling experiments are considered for which apparent discrepancies exist between fit quality expectations based on implemented quality assurance tactics, and the actual fit quality achieved in those experiments. These discrepancies are resolved by using Bayesian inference to account for certain imperfections in the assessment methodology. Estimates of the fraction of out-of-tolerance model predictions based on traditional frequentist methods are revised to account for uncertainty in the residual assessment process. The number of sites in the design space for which residuals are out of tolerance is seen to exceed the number of sites where the model actually fails to fit the data. A method is presented to estimate how much of the design space in inadequately modeled by low-order polynomial approximations to the true but unknown underlying response function.

DeLoach, Richard↗

NASA NTRS2022

P- and L-Band Retrieval of Subsurface Soil Moisture and Temperature Profiles as First-Order Polynomial Function

This paper demonstrates the potential use of P and L band passive measurements to determine root zone soil moisture (SM) and soil temperature(ST). SM and ST data have been taken as a function of depth during the NASA GSFC PLEX19 experiment in the summer of 2019 at Beltsville, MD, USA. Using these data, a coherent model has been used to compute H and V brightness temperatures at frequencies of 0.8 and 1.4 GHz with an observation angle of 35 degrees. These synthetic brightness data are then used to estimate the SM and ST profiles which are represented by linear polynomials. The inversion problem is formulated as a least square problem that is solved by a global optimization method known as the Adaptive Simulated Annealing(ASA) method. Four inversion examples having different SM and ST profiles are presented. Selected results show that the standard deviation between the retrieved and measured data is less than 0.077 cm3/cm3 for SM, and 2.245 °C for ST.

Ming Li↗

NASA NTRS1989

Theory of wide-angle photometry from standard stars

Wide angle celestial structures, such as bright comet tails and nearby galaxies and clusters of galaxies, rely on photographic methods for quantified morphology and photometry, primarily because electronic devices with comparable resolution and sky coverage are beyond current technological capability. The problem of the photometry of extended structures and of how this problem may be overcome through calibration by photometric standard stars is examined. The perfect properties of the ideal field of view are stated in the guise of a radiometric paraxial approximation, in the hope that fields of view of actual telescopes will conform. Fundamental radiometric concepts are worked through before the issue of atmospheric attenuation is addressed. The independence of observed atmospheric extinction and surface brightness leads off the quest for formal solutions to the problem of surface photometry. Methods and problems of solution are discussed. The spectre is confronted in the spirit of standard stars and shown to be chimerical in that light, provided certain rituals are adopted. After a brief discussion of Baker-Sampson polynomials and the vexing issue of saturation, a pursuit is made of actual numbers to be expected in real cases. While the numbers crunched are gathered ex nihilo, they demonstrate the feasibility of Newton's method in the solution of this overdetermined, nonlinear, least square, multiparametric, photometric problem.

Usher, Peter D.↗

NASA NTRS2013

An Extension of the Time-Spectral Method to Overset Solvers

Relative motion in the Cartesian or overset framework causes certain spatial nodes to move in and out of the physical domain as they are dynamically blanked by moving solid bodies. This poses a problem for the conventional Time-Spectral approach, which expands the solution at every spatial node into a Fourier series spanning the period of motion. The proposed extension to the Time-Spectral method treats unblanked nodes in the conventional manner but expands the solution at dynamically blanked nodes in a basis of barycentric rational polynomials spanning partitions of contiguously defined temporal intervals. Rational polynomials avoid Runge's phenomenon on the equidistant time samples of these sub-periodic intervals. Fourier- and rational polynomial-based differentiation operators are used in tandem to provide a consistent hybrid Time-Spectral overset scheme capable of handling relative motion. The hybrid scheme is tested with a linear model problem and implemented within NASA's OVERFLOW Reynolds-averaged Navier- Stokes (RANS) solver. The hybrid Time-Spectral solver is then applied to inviscid and turbulent RANS cases of plunging and pitching airfoils and compared to time-accurate and experimental data. A limiter was applied in the turbulent case to avoid undershoots in the undamped turbulent eddy viscosity while maintaining accuracy. The hybrid scheme matches the performance of the conventional Time-Spectral method and converges to the time-accurate results with increased temporal resolution.

Leffell, Joshua Isaac↗

NASA NTRS2017

On High-Order Upwind Methods for Advection

In the fourth installment of the celebrated series of five papers entitled "Towards the ultimate conservative difference scheme", Van Leer (1977) introduced five schemes for advection, the first three are piecewise linear, and the last two, piecewise parabolic. Among the five, scheme I, which is the least accurate, extends with relative ease to systems of equations in multiple dimensions. As a result, it became the most popular and is widely known as the MUSCL scheme (monotone upstream-centered schemes for conservation laws). Schemes III and V have the same accuracy, are the most accurate, and are closely related to current high-order methods. Scheme III uses a piecewise linear approximation that is discontinuous across cells, and can be considered as a precursor of the discontinuous Galerkin methods. Scheme V employs a piecewise quadratic approximation that is, as opposed to the case of scheme III, continuous across cells. This method is the basis for the on-going "active flux scheme" developed by Roe and collaborators. Here, schemes III and V are shown to be equivalent in the sense that they yield identical (reconstructed) solutions, provided the initial condition for scheme III is defined from that of scheme V in a manner dependent on the CFL number. This equivalence is counter intuitive since it is generally believed that piecewise linear and piecewise parabolic methods cannot produce the same solutions due to their different degrees of approximation. The finding also shows a key connection between the approaches of discontinuous and continuous polynomial approximations. In addition to the discussed equivalence, a framework using both projection and interpolation that extends schemes III and V into a single family of high-order schemes is introduced. For these high-order extensions, it is demonstrated via Fourier analysis that schemes with the same number of degrees of freedom 𝐾 per cell, in spite of the different piecewise polynomial degrees, share the same sets of eigenvalues and thus, have the same stability and accuracy. Moreover, these schemes are accurate to order 2𝐾−1, which is higher than the expected order of 𝐾.

high-order methods↗

NASA NTRSDate not supplied

Tailoring the Stiffness of Large Mirrors to Reduce Specific Optical Aberrations

Presented is a method that tailors the geometry of a mirror’s core structure as well as the mounting interface such that a desired set of aberrations is reduced. It will be shown that targeting specific aberrations as opposed to reducing overall surface figure error (SFE) is a valid method to reduce weight or improve performance. Analyzed first is the self-weight deflection of a thin circular plate of uniform thickness continuously supported along its outer edge. The Kirchhoff-Love plate bending equations are related to the radially symmetric Zernike polynomials of piston, defocus, and primary spherical aberration. Next, the mirror’s radial thickness is optimized to eliminate primary and secondary spherical aberrations when the outer edge is simply supported. As an example, the method is demonstrated using a 2-meter diameter aluminum mirror where surface deflections can have an appreciable effect at mid-to-near infrared wavelengths

mirror stiffness optimization↗