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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Quantification of Error Associated with Using Misaligned Meshes in Continuum Damage Mechanics Material Models for Matrix Crack Growth Predictions in Composites

The ability of a material model to capture in-plane matrix mode I and mode II crack growth is an essential component for modeling ply level damage evolution in composite structures. Previous studies using a continuum damage mechanics (CDM) approach have shown success in satisfying benchmark solutions for mode I and II crack growth. However, success was shown using a fiber-aligned meshing strategy, which encourages matrix cracks to propagate in a single band of elements, along the fiber direction. Generating a fiber-aligned mesh becomes a highly involved process for laminates including off-axis (non 0° or 90°) plies. The objective of this study is to quantify the effect of non-fiber aligned mesh discretization on predictions of in-plane matrix crack propagation. The approach taken incrementally varies the mesh orientation angle relative to the fiber orientation; more specifically, misaligned meshes are used to quantify the effect of element angle orientation relative to the initial crack orientation on the energy released during matrix crack propagation simulations using a CDM method. CDM solutions obtained with the misaligned meshes are evaluated against known benchmarks for mode I and II matrix crack growth. The CDM solutions reveal a near-polynomial trend of increased predicted failure stress with increased mesh misalignment angle; hence implying a potential relationship between element orientation angle and apparent fracture toughness.

Justusson, Brian↗
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 NTRS1992

Faster Computation Of Far Fields Of Dish Antennas

Modified Zernike-polynomial approach reduces number of terms in numerical integration. Yields greater or equal accuracy with fewer terms in series expansion and fewer sampling points across aperture and, less computation. Method based on Jacobi-Bessel-expansion concept.

Jamnejad, Vahraz↗
NASA NTRS2017

Fuzzy Logic Controller Stability Analysis Using a Satisfiability Modulo Theories Approach

While many widely accepted methods and techniques exist for validation and verification of traditional controllers, at this time no solutions have been accepted for Fuzzy Logic Controllers (FLCs). Due to the highly nonlinear nature of such systems, and the fact that developing a valid FLC does not require a mathematical model of the system, it is quite difficult to use conventional techniques to prove controller stability. Since safety-critical systems must be tested and verified to work as expected for all possible circumstances, the fact that FLC controllers cannot be tested to achieve such requirements poses limitations on the applications for such technology. Therefore, alternative methods for verification and validation of FLCs needs to be explored. In this study, a novel approach using formal verification methods to ensure the stability of a FLC is proposed. Main research challenges include specification of requirements for a complex system, conversion of a traditional FLC to a piecewise polynomial representation, and using a formal verification tool in a nonlinear solution space. Using the proposed architecture, the Fuzzy Logic Controller was found to always generate negative feedback, but inconclusive for Lyapunov stability.

Fuzzy Logic Controller↗
NASA NTRS2017

Fuzzy Logic Controller Stability Analysis Using a Satisfiability Modulo Theories Approach

While many widely accepted methods and techniques exist for validation and verification of traditional controllers, at this time no solutions have been accepted for Fuzzy Logic Controllers (FLCs). Due to the highly nonlinear nature of such systems, and the fact that developing a valid FLC does not require a mathematical model of the system, it is quite difficult to use conventional techniques to prove controller stability. Since safety-critical systems must be tested and verified to work as expected for all possible circumstances, the fact that FLC controllers cannot be tested to achieve such requirements poses limitations on the applications for such technology. Therefore, alternative methods for verification and validation of FLCs needs to be explored. In this study, a novel approach using formal verification methods to ensure the stability of a FLC is proposed. Main research challenges include specification of requirements for a complex system, conversion of a traditional FLC to a piecewise polynomial representation, and using a formal verification tool in a nonlinear solution space. Using the proposed architecture, the Fuzzy Logic Controller was found to always generate negative feedback, but inconclusive for Lyapunov stability.

Fuzzy Logic Controller↗
NASA NTRS1974

A quasi-vortex-lattice method in thin wing theory

A quasi-continuous method is developed for solving thin-wing problems. For the purpose of satisfying the wing boundary conditions, the spanwise vortex distribution is assumed to be stepwise-constant, while the chordwise vortex integral is reduced to a finite sum through a modified trapezoidal rule and the theory of Chebyshev polynomials. Wing-edge and Cauchy singularities are acounted for. The total aerodynamic characteristics are obtained by an appropriate quadrature integration. The two-dimensional results for airfoils without flap deflection reproduce the exact solutions in lift and pitching moment coefficients, the leading edge suction, and the pressure difference at a finite number of points. For a flapped airfoil, the present results are more accurate than those given by the vortex-lattice method. The three-dimensional results also show an improvement over the results of the vortex-lattice method. Extension to nonplanar applications is discussed.

Lan, C. E.↗
NASA NTRS1976

Numerical Airfoil Optimization Using a Reduced Number of Design Coordinates

A method is presented for numerical airfoil optimization whereby a reduced number of design coordinates are used to define the airfoil shape. The approach is to define the airfoil as a linear combination of shapes. These basic shapes may be analytically or numerically defined, allowing the designer to use his insight to propose candidate designs. The design problem becomes one of determining the participation of each such function in defining the optimum airfoil. Examples are presented for two-dimensional airfoil design and are compared with previous results based on a polynomial representation of the airfoil shape. Four existing NACA airfoils are used as basic shapes. Solutions equivalent to previous results are achieved with a factor of more than 3 improvements in efficiency, while superior designs are demonstrated with an efficiency greater than 2 over previous methods. With this shape definition, the optimization process is shown to exploit the simplifying assumptions in the inviscid aerodynamic analysis used here, thus demonstrating the need to use more advanced aerodynamics for airfoil optimization.

Vanderplaats, G. N.↗
NASA NTRS1976

Prediction of jump phenomena in roll-coupled maneuvers of airplanes

An easily computerized analytical method is developed for identifying critical airplane maneuvers in which nonlinear rotational coupling effects may cause sudden jumps in the response to pilot's control inputs. Fifth and ninth degree polynomials for predicting multiple pseudo-steady states of roll-coupled maneuvers are derived. The program calculates the pseudo-steady solutions and their stability. The occurrence of jump-like responses for several airplanes and a variety of maneuvers is shown to correlate well with the appearance of multiple stable solutions for critical control combinations. The analysis is extended to include aerodynamics nonlinear in angle of attack.

Schy, A. A.↗
NASA NTRS1977

Expressions for the precession quantities based upon the IAU /1976/ system of astronomical constants

The structure of the expressions usually employed in calculating the effects of precession is examined, and a method is outlined for revising the expressions to account for changes in the fundamental astronomical constants. It is shown that the basic set of parameters, upon which depend the lengthy polynomials for computing the mean obliquity of data and the elements of the precession matrix, consists of the mean obliquity, the speed of general precession in longitude at a fixed epoch, and the system of planetary masses. Special attention is given to the motion of the ecliptic pole, formulations for a basic epoch as well as an arbitrary epoch, and ecliptic motion relative to the basic epoch. Numerical precession quantities at epoch J2000.0 (JED 2451545.0) are presented which result from the revision of astronomical constants adopted at the XVI General Assembly of the IAU.

Lieske, J. H.↗
NASA NTRS1999

Development of Response Surface Models for Rapid Analysis & Multidisciplinary Optimization of Launch Vehicle Design Concepts

Multdisciplinary design optimization (MDO) is an important step in the design and evaluation of launch vehicles, since it has a significant impact on performance and lifecycle cost. The objective in MDO is to search the design space to determine the values of design parameters that optimize the performance characteristics subject to system constraints. Vehicle Analysis Branch (VAB) at NASA Langley Research Center has computerized analysis tools in many of the disciplines required for the design and analysis of launch vehicles. Vehicle performance characteristics can be determined by the use of these computerized analysis tools. The next step is to optimize the system performance characteristics subject to multidisciplinary constraints. However, most of the complex sizing and performance evaluation codes used for launch vehicle design are stand-alone tools, operated by disciplinary experts. They are, in general, difficult to integrate and use directly for MDO. An alternative has been to utilize response surface methodology (RSM) to obtain polynomial models that approximate the functional relationships between performance characteristics and design variables. These approximation models, called response surface models, are then used to integrate the disciplines using mathematical programming methods for efficient system level design analysis, MDO and fast sensitivity simulations. A second-order response surface model of the form given has been commonly used in RSM since in many cases it can provide an adequate approximation especially if the region of interest is sufficiently limited.

Unal, Resit↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1973

Experimental and data analysis techniques for deducing collision-induced forces from photographic histories of engine rotor fragment impact/interaction with a containment ring

An analysis method termed TEJ-JET is described whereby measured transient elastic and inelastic deformations of an engine-rotor fragment-impacted structural ring are analyzed to deduce the transient external forces experienced by that ring as a result of fragment impact and interaction with the ring. Although the theoretical feasibility of the TEJ-JET concept was established, its practical feasibility when utilizing experimental measurements of limited precision and accuracy remains to be established. The experimental equipment and the techniques (high-speed motion photography) employed to measure the transient deformations of fragment-impacted rings are described. Sources of error and data uncertainties are identified. Techniques employed to reduce data reading uncertainties and to correct the data for optical-distortion effects are discussed. These procedures, including spatial smoothing of the deformed ring shape by Fourier series and timewise smoothing by Gram polynomials, are applied illustratively to recent measurements involving the impact of a single T58 turbine rotor blade against an aluminum containment ring. Plausible predictions of the fragment-ring impact/interaction forces are obtained by one branch of this TEJ-JET method; however, a second branch of this method, which provides an independent estimate of these forces, remains to be evaluated.

Yeghiayan, R. P.↗
NASA NTRS1975

Geometric analysis and restitution of digital multispectral scanner data arrays

An investigation was conducted to define causes of geometric defects within digital multispectral scanner (MSS) data arrays, to analyze the resulting geometric errors, and to investigate restitution methods to correct or reduce these errors. Geometric transformation relationships for scanned data, from which collinearity equations may be derived, served as the basis of parametric methods of analysis and restitution of MSS digital data arrays. The linearization of these collinearity equations is presented. Algorithms considered for use in analysis and restitution included the MSS collinearity equations, piecewise polynomials based on linearized collinearity equations, and nonparametric algorithms. A proposed system for geometric analysis and restitution of MSS digital data arrays was used to evaluate these algorithms, utilizing actual MSS data arrays. It was shown that collinearity equations and nonparametric algorithms both yield acceptable results, but nonparametric algorithms possess definite advantages in computational efficiency. Piecewise polynomials were found to yield inferior results.

Baker, J. R.↗
NASA NTRS1986

Total ionospheric electron content calibration using SERIES GPS satellite data

The current status of the Deep Space Network advanced systems research into ionospheric calibration techniques, based on Global Positioning System (GPS) data is described. A GPS-based calibration system is planned to replace the currently used Faraday rotation method by 1989. The SERIES receiver system determines the differential group delay of signals transmitted at two different carrier frequencies. This differential delay includes an ionospheric component and a GPS transmitter offset. The transmitter offsets are different for each GPS satellite. Tests were conducted to assess the effect of the offsets on the ionospheric calibration accuracy. From the obtained data, the total electron content and GPS transmitter offsets are calculated by a least squares estimation method employing a local model of total ionospheric electron content. The end product is an estimation of the total ionospheric content for an arbitrary line-of-sight direction. For the presented polynomial fitting technique, the systematic error due to mismodeling is estimated to be approximately 6 x 10 to the 16th power el/sq m, while the formal error is approximately 2 x 10 to the 16th power el/sq m. The final goal is an error of 3 x 10 to the 16th power el/sq m (approximately 0.7 ns at 2.3 GHz).

Lanyi, G.↗
NASA NTRS1990

The CFL condition for spectral approximations to hyperbolic initial-boundary value problems

The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.

Gottlieb, David↗
NASA NTRS1991

The CFL condition for spectral approximations to hyperbolic initial-boundary value problems

The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.

Gottlieb, David↗
NASA NTRS2021

Linear and Nonlinear Trend Analyzes in Global Satellite‐Based Precipitation, 1998–2017

Precipitation varies spatio‐temporally in amount, intensity, and frequency. Although, much research has been conducted on analyzing precipitation patterns and variation at the global scale, trend types have still not received much attention. This study developed a new polynomial‐based model for detecting nonlinear and linear trends in a satellite precipitation product (TRMM 3B43) for the 1998–2017 period at a near‐global scale. We used an automated trend classification method that detects significant trends and classifies them into linear and nonlinear (cubic, quadratic, and concealed) trend types in satellite‐based precipitation at near‐global, continental, and climate zone scales. We found that 12.3% of pixel‐based precipitation time series across the globe have significant trend at 0.05 significance level (50% positive and 50% negative trends). In all continents except Asia, decreasing trends were found to cover larger areas than corresponding increasing trends. Regarding climate zone and precipitation trend change, our results indicate that a linear trend is dominant in the warm temperate (77.7%) and equatorial climates (80.4%) while the least linear change was detected in the polar climate (68.9%). The combined results of continental and climate zone scales indicate significant increasing trends in Asia and arid climate over the last 20 yr. Furthermore, positive trends were found to be more significant at the continental scale, particularly, in Asia relative to the climate zone scale. Linear change in precipitation (80%) was the most dominant trend observed as opposed to nonlinear (quadratic [11%] and cubic [9%]) trend types at the global scale.

Majid Kazemzadeh↗
NASA NTRS1993

Issac, Jason Cherian ses in transonic flow

Flutter analysis of a two degree of freedom airfoil in compressible flow is performed using a state-space representation of the unsteady aerodynamic behavior. Indicial response functions are used to represent the normal force and moment response of the airfoil. The structural equations of motion of the airfoil with bending and torsional degrees of freedom are coupled to the unsteady air loads and the aeroelastic system so modelled is solved as an eigenvalue problem to determine the stability. The aeroelastic equations are also directly integrated with respect to time and the time-domain results compared with the results from the eigenanalysis. A good agreement is obtained. The derivatives of the flutter speed obtained from the eigenanalysis are calculated with respect to the mass and stiffness parameters by both analytical and finite-difference methods for various transonic Mach numbers. The experience gained from the two degree of freedom model is applied to study the sensitivity of the flutter response of a wing with respect to various shape parameters. The parameters being considered are as follows: (1) aspect ratio; (2) surface area of the wing; (3) taper ratio; and (4) sweep. The wing deflections are represented by Chebyshev polynomials. The compressible aerodynamic state-space model used for the airfoil section is extended to represent the unsteady aerodynamic forces on a generally laminated tapered skewed wing. The aeroelastic equations are solved as an eigenvalue problem to determine the flutter speed of the wing. The derivatives of the flutter speed with respect to the shape parameters are calculated by both analytical and finite difference methods.

Issac, Jason Cherion↗