Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Polynomial models”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Cutting Corners: Curvilinear-Surface-Based Gravity Models for Asteroids and Comets

We report that asteroids and comets are often highly irregular in shape and can contain large density inhomogenities. For such bodies, the simple point-mass gravity model does not capture the dynamics of the environment and this has driven a need for higher fidelity alternatives. A number of approaches have been developed with one of the most popular being the analytic polyhedral model. The analytic polyhedral model approximates the irregular shaped body as a polyhedron, typically consisting of triangular facets, with some assumed internal density distribution. Several analytic formulations have been derived for the gravitational fields of a constant-density polyhedron. A thorough history of the model is provided in Ref. [6] as well as a comparison of the implementations of Refs. [2–4]. The three versions provide near identical accuracy with the implementation of Werner requiring the fewest transcendental function evaluations suggesting its superior efficiency. Implementations with linear density contrasts and arbitrary polynomial density contrasts have also been developed.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1981

Atmospheric corrections to satellite radiometric data over rugged terrain

Radiometric measurements from satellites in the solar portion of the electromagnetic spectrum can be converted to measurements of surface exitance. Over rugged terrain, the satellite image must be precisely registered to a terrain data set. For small areas a first-order polynomial interpolation scheme is generally satisfactory for the geometric rectification. If there are saturated pixels, a nearest-neighbor procedure is used for the interpolated satellite radiance numbers; otherwise a cubic-convolution algorithm is used. Path radiance and path transmittance are calculated with a simple spectral model, which requires an estimate of the water vapor and aerosol content of the atmosphere.

Dozier, J.↗
NASA NTRS1986

Simplified curve fits for the thermodynamic properties of equilibrium air

New improved curve fits for the thermodynamic properties of equilibrium air were developed. The curve fits are for p = p(e,rho), a = a(e,rho), T = T(e,rho), s = s(e,rho), T = T(p,rho), h = h(p,rho), rho = rho(p,s), e = e(p,s) and a = a(p,s). These curve fits can be readily incorporated into new or existing Computational Fluid Dynamics (CFD) codes if real-gas effects are desired. The curve fits were constructed using Grabau-type transition functions to model the thermodynamic surfaces in a piecewise manner. The accuracies and continuity of these curve fits are substantially improved over those of previous curve fits appearing in NASA CR-2470. These improvements were due to the incorporation of a small number of additional terms in the approximating polynomials and careful choices of the transition functions. The ranges of validity of the new curve fits are temperatures up to 25,000 K and densities from 10 to the minus 7th to 100 amagats (rho/rho sub 0).

Srinivasan, S.↗
NASA NTRS2006

Analytical Derivation of the Vegetation Optical Depth from the Microwave Polarization Difference Index

A numerical solution for the canopy optical depth in an existing microwave-based land surface parameter retrieval model is presented. The optical depth is derived from the microwave polarization difference index and the dielectric constant of the soil. The original procedure used an approximation in the form of a logarithmic decay function to define this relationship, and was derived through a series of lengthy polynomials. These polynomials had to be recalculated when the scattering albedo or antenna incidence angle changes. The new procedure is computationally more efficient and accurate.

Meesters, Antoon G. C. A.↗
NASA NTRS2018

Influence of Core and Hexapod Geometry, and Local Reinforcement on the Performance of Ultra-Lightweight ULE Mirror

The Habitable Exoplanet Observatory (HabEx) mission has unique optical performance requirements which drive the mirror design process beyond the traditional criteria. While mass and stiffness are still important, the response to inertia loading (expressed in terms of Zernike coefficients) to omni-directional excitation dominates the effort. While a Zerodur mirror is the current baseline, as mass budgets change, a ULE design is being studied as a potential alternative. This trade study looked at over 264 design variations using the Arnold Mirror Modeler and ANSYS(c) to investigate the influence of various design elements, including: substrate thickness, core cell size, hexapod geometry and local reinforcement. Design 'goodness' was evaluated based on the mirror's inertial deformation response to omni-directional input. This response was calculated via RSSing Zernike polynomial responses to (XYZ) accelerations.

Space Telescope Technology↗
NASA NTRS2000

Motion Cueing Algorithm Development: Initial Investigation and Redesign of the Algorithms

In this project four motion cueing algorithms were initially investigated. The classical algorithm generated results with large distortion and delay and low magnitude. The NASA adaptive algorithm proved to be well tuned with satisfactory performance, while the UTIAS adaptive algorithm produced less desirable results. Modifications were made to the adaptive algorithms to reduce the magnitude of undesirable spikes. The optimal algorithm was found to have the potential for improved performance with further redesign. The center of simulator rotation was redefined. More terms were added to the cost function to enable more tuning flexibility. A new design approach using a Fortran/Matlab/Simulink setup was employed. A new semicircular canals model was incorporated in the algorithm. With these changes results show the optimal algorithm has some advantages over the NASA adaptive algorithm. Two general problems observed in the initial investigation required solutions. A nonlinear gain algorithm was developed that scales the aircraft inputs by a third-order polynomial, maximizing the motion cues while remaining within the operational limits of the motion system. A braking algorithm was developed to bring the simulator to a full stop at its motion limit and later release the brake to follow the cueing algorithm output.

Telban, Robert J.↗
NASA NTRS1992

Polynomial coefficients for calculating O2 Schumann-Runge cross sections at 0.5/cm resolution

O2 cross sections from 49,000 to 57,000/cm have been fitted with temperature dependent polynomial expressions, providing an accurate and efficient means of determining Schumann-Runge band cross sections for temperatures between 130 and 500 K. The least squares fits were carried out on a 0.5/cm spectral grid, using cross sections obtained from a Schumann-Runge line-by-line model that incorporates the most recent spectroscopic data. The O2 cross sections do not include the underlying Herzberg continuum, but they do contain contributions from the temperature dependent Schumann-Runge continuum. The cross sections are suitable for use in UV transmission calculations at high spectral resolution. They should also prove useful for updating existing parameterizations of ultraviolet transmission and O2 photolysis.

Minschwaner, K.↗
DOE OSTI2021

Time-Resolved X-ray Operando Observations of Lithiation Gradients across the Cathode Matrix and Individual Oxide Particles during Fast Cycling of a Li-Ion Cell

Lithiated transition metal oxides serve as active materials in the positive electrode (cathode) of lithium-ion cells. During electrochemical cycling, lithium ions intercalate and deintercalate into these oxide particles. This behavior causes two types of lithiation gradients to emerge: (i) a bulk gradient across the depth of the cathode matrix (averaged over individual oxide particles) and (ii) a microscopic gradient across the particles themselves, which also depends on their location in the electrode. Here we show how both gradients can be studied using operando X-ray diffraction during 4C charge and 4C discharge. The oxide (de)lithiation is estimated from the unit cell parameters by indexing the X-ray diffraction spectra. By fitting the lithiation profiles with orthogonal polynomials, the bulk gradients across the electrode thickness are quantified. These gradients develop as the current flows through the cell and dissipate during open-circuit and potentiostatic-hold periods. Further details of lithiation dynamics can be obtained through shape analysis of the Bragg peaks. In particular, from electrochemical model simulations, we show that the width and skewness of the (003) peak track (de)lithiation fronts moving across the individual oxide particles.

(003) peak↗
NASA NTRS1990

Analysis of aircraft tires via semianalytic finite elements

A computational procedure is presented for the geometrically nonlinear analysis of aircraft tires. The tire was modeled by using a two-dimensional laminated anisotropic shell theory with the effects of variation in material and geometric parameters included. The four key elements of the procedure are: (1) semianalytic finite elements in which the shell variables are represented by Fourier series in the circumferential direction and piecewise polynomials in the meridional direction; (2) a mixed formulation with the fundamental unknowns consisting of strain parameters, stress-resultant parameters, and generalized displacements; (3) multilevel operator splitting to effect successive simplifications, and to uncouple the equations associated with different Fourier harmonics; and (4) multilevel iterative procedures and reduction techniques to generate the response of the shell.

Noor, Ahmed K.↗
NASA NTRSDate not supplied

Jig Twist Optimization of Mach 0.8 Transonic Truss-Braced Wing Aircraft

This paper presents a jig twist optimization study of Mach 0.8 Transonic Truss-Braced Wing (TTBW) aircraft using FUN3D. The solver has been validated against wind tunnel test data. An aeroelastic model deformation correction method is developed for the wind tunnel model to improve the validation. The flight condition corresponding to Mach 0.8 and design lift coefficient 0.695 is selected for the TTBW aircraft jig twist optimization to reduce the drag coefficient. The static aeroelastic simulation module in FUN3D is used for the aeroelastic simulation. Three different orders of the Chebyshev polynomials are used as shape function to represent the variations of the jig twist distribution along the wing span. The preliminary optimization results show that the Mach 0.8 TTBW aircraft with the optimized jig twist achieves a drag reduction of about 5.7 counts at the design condition when 5 𝑡ℎ order Chebyshev polynomials are used as shape function.

TTBW↗
DOE OSTI2023

Applying Quantum Computing to Simulate Power System Dynamics

Power system dynamics are generally modeled by high dimensional nonlinear differential-algebraic equations due to a large number of generators, loads, and transmission lines. Thus, its computational complexity grows exponentially with the system size. This paper demonstrates the potential use of quantum computing algorithms to model the power system dynamics. Leveraging a symbolic programming framework, we equivalently convert the power system dynamics’ differential algebraic equations (DAEs) into ordinary differential equations (ODEs), where the data of the state vector can be encoded into quantum computers via amplitude encoding. The system's nonlinearity is captured by Taylor polynomial expansion, the quantum state tensor, and Hamiltonian simulation, whereas state variables can be updated by a quantum linear equation solver. Our results show that quantum computing can simulate the dynamics of the power system with high accuracy, whereas its complexity is polynomial in the logarithm of the system dimension. Our work also illustrates the use of scientific machine learning tools for implementing scientific computing concepts, e.g., Taylor expansion, DAEs/ODEs transform, and quantum computing solver, in the field of power engineering.

Tran, Huynh↗
DOE OSTI2021

Complexity growth in integrable and chaotic models

We use the SYK family of models with N Majorana fermions to study the complexity of time evolution, formulated as the shortest geodesic length on the unitary group manifold between the identity and the time evolution operator, in free, integrable, and chaotic systems. Initially, the shortest geodesic follows the time evolution trajectory, and hence complexity grows linearly in time. We study how this linear growth is eventually truncated by the appearance and accumulation of conjugate points, which signal the presence of shorter geodesics intersecting the time evolution trajectory. By explicitly locating such “shortcuts” through analytical and numerical methods, we demonstrate that: (a) in the free theory, time evolution encounters conjugate points at a polynomial time; consequently complexity growth truncates at O($\sqrt{N}$), and we find an explicit operator which “fast-forwards” the free N-fermion time evolution with this complexity, (b) in a class of interacting integrable theories, the complexity is upper bounded by O(poly(N)), and (c) in chaotic theories, we argue that conjugate points do not occur until exponential times O(e N ), after which it becomes possible to find infinitesimally nearby geodesics which approximate the time evolution operator. Finally, we explore the notion of eigenstate complexity in free, integrable, and chaotic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Results on elastic cross sections in proton–proton collisions at $\sqrt{s}$ = 510 GeV with the STAR detector at RHIC

We report results on an elastic cross section measurement in proton–proton collisions at a center-of-mass energy $\sqrt{s}$ = 510 GeV, obtained with the Roman Pot setup of the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The elastic differential cross section is measured in the four-momentum transfer squared range 0.23 ≤ −t ≤ 0.67 GeV 2 . This is the only measurement of the proton-proton elastic cross section in this t range for collision energies above the Intersecting Storage Rings (ISR) and below the Large Hadron Collider (LHC) colliders. We find that a constant slope B does not fit the data in the aforementioned t range, and we obtain a much better fit using a second-order polynomial for B(t). This is the first measurement below the LHC energies for which the non-constant behavior B(t) is observed. The t dependence of B is also determined using six subintervals of t in the STAR measured t range, and is in good agreement with the phenomenological models. The measured elastic differential cross section dσ/dt agrees well with the results obtained at $\sqrt{s}$ = 540 GeV for proton–antiproton collisions by the UA4 experiment. We also determine that the integrated elastic cross section within the STAR t-range is σ$^{fid}_{el}$ = 462.1 ± 0.9(stat.) ± 1.1(syst.) ± 11.6(scale) μb.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Coherency-Constrained Spectral Clustering for Power Network Reduction

This paper presents a methodology for reducing the complexity of large-scale power network models using spectral clustering, aggregation of electrical components, and cost function approximation. Two approaches are explored using unconstrained and constrained spectral clustering to determine areas for effective system reduction. Once the system areas are determined, both loads and generators by type are aggregated, and their new cost function is approximated through polynomial curve-fitting or statistical methods. The performance of reduced networks is evaluated in terms of their ability to follow the true daily cost of the original system over a 24-hour period considering a set of several days. Two test systems are taken as test beds. Application of the methodology to a modified version of the IEEE 39-bus system reduces it from 17 generators to a 4-bus system and 9 generators with about 93% of accuracy. Similarly, the IEEE 118-bus system is reduced from 19 generators to a 3-bus system with three aggregated units achieving over 99% of accuracy. These findings address scalability challenges and enhance accuracy for high and mid-loading level conditions, and by aggregating thermal units with similar cost functions.

42 ENGINEERING↗
DOE OSTI2021

Comprehensive Material Characterization and Simultaneous Model Calibration for Improved Computational Simulation Credibility

Computational simulation is increasingly relied upon for high-consequence engineering decisions, and a foundational element to solid mechanics simulations is a credible material model. Our ultimate vision is to interlace material characterization and model calibration in a real-time feedback loop, where the current model calibration results will drive the experiment to load regimes that add the most useful information to reduce parameter uncertainty. The current work investigated one key step to this Interlaced Characterization and Calibration (ICC) paradigm, using a finite load-path tree to incorporate history/path dependency of nonlinear material models into a network of surrogate models that replace computationally-expensive finite-element analyses. Our reference simulation was an elastoplastic material point subject to biaxial deformation with a Hill anisotropic yield criterion. Training data was generated using either a space-filling or adaptive sampling method, and surrogates were built using either Gaussian process or polynomial chaos expansion methods. Surrogate error was evaluated to be on the order of 10 ⁻5 and 10 ⁻3 percent for the space-filling and adaptive sampling training data, respectively. Direct Bayesian inference was performed with the surrogate network and with the reference material point simulator, and results agreed to within 3 significant figures for the mean parameter values, with a reduction in computational cost over 5 orders of magnitude. These results bought down risk regarding the surrogate network and facilitated a successful FY22-24 full LDRD proposal to research and develop the complete ICC paradigm.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Line shape of the non-thermal 6300 A O/1D/ emission

The two-population model of Schmitt, Abreu and Hays (1981) is used to calculate the line shape of the atomic oxygen metastable state, nonthermal O(1D) 6300 A emission, in order to simulate observations made from a space platform at different zenith angles and altitudes. The Addition theorem, for spherical harmonics of a Legendre polynomial expansion of the nonthermal population distribution function, is used to obtain nonthermal line shapes observed at zenith angles other than the local vertical one.

Schmitt, G. A.↗
NASA NTRS1995

Estimating errors in least-squares fitting

While least-squares fitting procedures are commonly used in data analysis and are extensively discussed in the literature devoted to this subject, the proper assessment of errors resulting from such fits has received relatively little attention. The present work considers statistical errors in the fitted parameters, as well as in the values of the fitted function itself, resulting from random errors in the data. Expressions are derived for the standard error of the fit, as a function of the independent variable, for the general nonlinear and linear fitting problems. Additionally, closed-form expressions are derived for some examples commonly encountered in the scientific and engineering fields, namely ordinary polynomial and Gaussian fitting functions. These results have direct application to the assessment of the antenna gain and system temperature characteristics, in addition to a broad range of problems in data analysis. The effects of the nature of the data and the choice of fitting function on the ability to accurately model the system under study are discussed, and some general rules are deduced to assist workers intent on maximizing the amount of information obtained form a given set of measurements.

Richter, P. H.↗