Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “polynomial expansions”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Kinetic theory of turbulent compressible flows and comparison with classical theory

A complete set of turbulent correlations is given in terms of expansion coefficients in a double series of Hermite polynomials of a two-particle correlation in phase space. Only two of these coefficients, corresponding to Reynolds stress and turbulent heat flux, are shown to appear in gasdynamic equations of turbulence. This is accomplished by incorporating the phase-space correlation in the collision integral of Boltzmann's equation and by deriving generalized Navier-Stokes and Fourier transport relations. In contrast to the classical formalism, the simple expression remains invariant in form whether the flow is compressible or not. Nevertheless, a bilinear transformation of fluctuating quantities shows that the two formalisms are identical with regard to terms of double correlations. The kinetic theory justifies nonexistence of the higher-order correlations characteristic of compressible-turbulence equations in the classical regime.

Tsuge, S.↗
NASA NTRS1987

Inverse solution of Kepler's equation for hyperbolic orbits

An algorithm is presented for efficient inverse solution of Kepler's equation for hyperbolic orbits. It is shown that an expansion of Barker's equation into a bicubic polynomial provides a good approximation to obtain accurate starting values for rapid numerical solution of Kepler's equation. In the approximate equation a cubic in normalized elapsed flight time from pericenter is set equal to a cubic in a function S of eccentricity and true anomaly. The initial estimate of S to use in an iteration formula is obtained by evaluating the cubic in normalized flight time and finding in most cases the single real root of the other cubic. This initial estimate has an accuracy corresponding to values of true anomaly in error by less than 0.5 degrees generally.

Boltz, Frederick W.↗
NASA NTRS2013

RUSHMAPS: Real-Time Uploadable Spherical Harmonic Moment Analysis for Particle Spectrometers

RUSHMAPS is a new onboard data reduction scheme that gives real-time access to key science parameters (e.g. moments) of a class of heliophysics science and/or solar system exploration investigation that includes plasma particle spectrometers (PPS), but requires moments reporting (density, bulk-velocity, temperature, pressure, etc.) of higher-level quality, and tolerates a lowpass (variable quality) spectral representation of the corresponding particle velocity distributions, such that telemetry use is minimized. The proposed methodology trades access to the full-resolution velocity distribution data, saving on telemetry, for real-time access to both the moments and an adjustable-quality (increasing quality increases volume) spectral representation of distribution functions. Traditional onboard data storage and downlink bandwidth constraints severely limit PPS system functionality and drive cost, which, as a consequence, drives a limited data collection and lower angular energy and time resolution. This prototypical system exploit, using high-performance processing technology at GSFC (Goddard Space Flight Center), uses a SpaceCube and/or Maestro-type platform for processing. These processing platforms are currently being used on the International Space Station as a technology demonstration, and work is currently ongoing in a new onboard computation system for the Earth Science missions, but they have never been implemented in heliospheric science or solar system exploration missions. Preliminary analysis confirms that the targeted processor platforms possess the processing resources required for realtime application of these algorithms to the spectrometer data. SpaceCube platforms demonstrate that the target architecture possesses the sort of compact, low-mass/power, radiation-tolerant characteristics needed for flight. These high-performing hybrid systems embed unprecedented amounts of onboard processing power in the CPU (central processing unit), FPGAs (field programmable gate arrays), and DSP (digital signal processing) elements. The fundamental computational algorithm de constructs 3D velocity distributions in terms of spherical harmonic spectral coefficients (which are analogous to a Fourier sine-cosine decomposition), but uses instead spherical harmonics Legendre polynomial orthogonal functions as a basis for the expansion, portraying each 2D angular distribution at every energy or, geometrically, spherical speed-shell swept by the particle spectrometer. Optionally, these spherical harmonic spectral coefficients may be telemetered to the ground. These will provide a smoothed description of the velocity distribution function whose quality will depend on the number of coefficients determined. Successfully implemented on the GSFC-developed processor, the capability to integrate the proposed methodology with both heritage and anticipated future plasma particle spectrometer designs is demonstrated (with sufficiently detailed design analysis to advance TRL) to show specific science relevancy with future HSD (Heliophysics Science Division) solar-interplanetary, planetary missions, sounding rockets and/or CubeSat missions.

Figueroa-Vinas, Adolfo↗
DOE OSTI2022

Extended kinetic lattice grand canonical Monte Carlo simulation method for transport of multicomponent ion mixtures through a model nanopore system

Abstract An extended version of the original kinetic lattice grand canonical Monte Carlo simulation method combined with mean field theory (KLGCMC/MF) ( J. Chem. Phys. 2007 , 127, 024706) is presented for the study of transport of multicomponent ion mixtures through a model nanopore. Comparison of the extended KLGCMC/MF (eKLGCMC/MF) simulation results with Poisson–Nernst–Planck (PNP) calculations is also made to confirm the validity of the extended simulation approach. Unlike the original version of KLGCMC/MF simulation method that treats only a binary ionic solution with one cation and one anion species, this extended version can deal with a system that includes ternary ion mixtures. A diffusion probability algorithm is also added to the extended version of the simulation method to describe the inhomogeneous diffusivity of ions that is often observed in the ion permeation through nanopores. Both Legendre and Chebyshev polynomials of the second kind were tested as a basis set for the basis set expansion (BSE) method with which to calculate the reaction field energy in the eKLGCMC/MF simulation. It turned out that the Legendre polynomials perform better than the Chebyshev polynomials, and as a result, the Legendre polynomials were implemented in the current version of eKLGCMC/MF simulation algorithm. The presented eKLGCMC/MF simulation method with new features finds its potential applications in nanopore systems where the correlation between ion species with the same sign of charges plays a key role such as oscillating ion currents or anomalous mole fraction effects.

Choi, Inhyeok↗
DOE OSTI2023

Photon production rate from Transverse-Longitudinal ( T – L) mesonic correlator on the lattice

Thermal photons from the QGP provide important information about the interaction among plasma constituents. The photon production rate from a thermally equilibrated system is proportional to the transverse spectral function ρ T (ω = |k → |,k → ). One can also calculate the photon production rate from the difference between ρ T (ω, k → ) (transverse) and ρ L (ω, k → ) (longitudinal) projections, as ρL vanishes on the photon point. Because the UV part of ρ T – ρ L is suppressed, the corresponding Euclidean correlator receives most of its contribution from the IR part. We calculate the T – L correlator on N f = 2 + 1 flavour HISQ configurations with m l = m s /5 at temperature of about 1.15T ρc (220 MeV). We have used two ansätze for the spectral function: 1) A polynomial connected to the UV region consistent with OPE expansion and 2) a hydro-inspired spectral function. We have also applied the Backus-Gilbert method to estimate the spectral function. All these different approaches are combined to estimate the photon production rate.

Bala, Dibyendu↗
NASA NTRS1984

Aberrations of varied line-space grazing incidence gratings in converging light beams

Analyses of the imaging properties of several designs for varied-line space gratings in converging beams of light in grazing-incidence spectrometers are presented. An explicit model is defined for the case of a plane-reflection grating intercepting light that converges and is reflected to a stigmatic point associated with the zero-order image of the grating. Smooth spatial variation of the grating constant then permits aberration correction. The aberrations are expressed as polynomials in the grating lens coordinates using power series expansions. Application of the model is illustrated in terms of aberrations experienced with the short wavelength spectrometer on the EUVE satellite. Attention is given to straight and parallel in-plane grooves, curved groove in-plane designs and off-plane grooves. Aberrations due to dispersions and misalignment are also considered.

Hettrick, M. C.↗
NASA NTRS1994

An Efficient Spectral Method for Ordinary Differential Equations with Rational Function Coefficients

We present some relations that allow the efficient approximate inversion of linear differential operators with rational function coefficients. We employ expansions in terms of a large class of orthogonal polynomial families, including all the classical orthogonal polynomials. These families obey a simple three-term recurrence relation for differentiation, which implies that on an appropriately restricted domain the differentiation operator has a unique banded inverse. The inverse is an integration operator for the family, and it is simply the tridiagonal coefficient matrix for the recurrence. Since in these families convolution operators (i.e. matrix representations of multiplication by a function) are banded for polynomials, we are able to obtain a banded representation for linear differential operators with rational coefficients. This leads to a method of solution of initial or boundary value problems that, besides having an operation count that scales linearly with the order of truncation N, is computationally well conditioned. Among the applications considered is the use of rational maps for the resolution of sharp interior layers.

Coutsias, Evangelos A.↗
DOE OSTI2022

Analysis of sparse recovery for Legendre expansions using envelope bound

We provide novel sufficient conditions for the uniform recovery of sparse Legendre expansions using ℓ 1 minimization, where the sampling points are drawn according to orthogonalization (uniform) measure. So far, conditions of the form m ≳ Θ 2 s x log factors have been relied on to determine the minimum number of samples m that guarantees successful reconstruction of s-sparse vectors when the measurement matrix is associated to an orthonormal system. However, in case of sparse Legendre expansions, the uniform bound Θ of Legendre systems is so high that these conditions are unable to provide meaningful guarantees. Here, in this paper, we present an analysis which employs the envelop bound of all Legendre polynomials instead, and prove a new recovery guarantee for s-sparse Legendre expansions, m ≳ Θs 2 x log factors, which is independent of Θ. Arguably, this is the first recovery condition established for orthonormal systems without assuming the uniform boundedness of the sampling matrix. The key ingredient of our analysis is an extension of chaining arguments, recently developed in Bourgain and Chkifa et al., to handle the envelope bound. Furthermore, our recovery condition is proved via restricted eigenvalue property, a less demanding replacement of restricted isometry property which is perfectly suited to the considered scenario. Along the way, we derive simple criteria to detect good sample sets. Our numerical tests show that sets of uniformly sampled points that meet these criteria will perform better recovery on average.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Systematic, Polynomial-Cost Approach to Exact Correlation Energies (Final Technical Report)

The full configuration interaction (FCI) wave function provides the exact solution to the Schrödinger equation in a given basis set. While FCI is intractable due to exponential computational costs, the many-body expansion (or the method of increments) can reduce scaling to a low-order polynomial with system size. This project entailed advances in the incremental FCI (iFCI) approach, designed to allow iFCI to reach larger system sizes and maintain its intrinsic high accuracy. This document describes advances in solvers for iFCI, strategies to treat multiple charge and spin states, and virtual state management methods to reduce memory requirements. Overall, this project allows iFCI to correlate (for the first time) 142 valence electrons in 444 orbitals in a realistic model of a transition metal complex.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2004

Validation of an Impedance Education Method in Flow

This paper reports results of a research effort to validate a method for educing the normal incidence impedance of a locally reacting liner, located in a grazing incidence, nonprogressive acoustic wave environment with flow. The results presented in this paper test the ability of the method to reproduce the measured normal incidence impedance of a solid steel plate and two soft test liners in a uniform flow. The test liners are known to be locally react- ing and exhibit no measurable amplitude-dependent impedance nonlinearities or flow effects. Baseline impedance spectra for these liners were therefore established from measurements in a conventional normal incidence impedance tube. A key feature of the method is the expansion of the unknown impedance function as a piecewise continuous polynomial with undetermined coefficients. Stewart's adaptation of the Davidon-Fletcher-Powell optimization algorithm is used to educe the normal incidence impedance at each Mach number by optimizing an objective function. The method is shown to reproduce the measured normal incidence impedance spectrum for each of the test liners, thus validating its usefulness for determining the normal incidence impedance of test liners for a broad range of source frequencies and flow Mach numbers. Nomenclature

Watson, Willie R.↗
DOE OSTI2024

The linear thermal expansion of 11 polymers from approximately -100 to +100°C

The linear coefficient of thermal expansion has been measured for 11 polymers from approximately -100 to +100°C. The polymers tested were poly(tetrafluoroethylene) (PTFE), poly(chlorotrifluoroethylene) (PCTFE), three polyethylenes; high density (HDPE), ultra-high molecular weight (UHMWPE) and cross-linked (XPE). The commercial fluoropolymers Kel-F 800 (also known as FK-800®) and THV 500®were tested together with samples of poly(ether ether ketone) (PEEK), poly(methyl methyl methacrylate) (PMMA), polycarbonate (PC) and poly(vinylidene fluoride) (PVDF). Polynomial fits were made to the data over appropriate ranges allowing the expansivity (α) to be estimated.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Efficient Probabilistic Prediction and Uncertainty Quantification of Tropical Cyclone–Driven Storm Tides and Inundation

Abstract This study proposes and assesses a methodology to obtain high-quality probabilistic predictions and uncertainty information of near-landfall tropical cyclone–driven (TC-driven) storm tide and inundation with limited time and resources. Forecasts of TC track, intensity, and size are perturbed according to quasi-random Korobov sequences of historical forecast errors with assumed Gaussian and uniform statistical distributions. These perturbations are run in an ensemble of hydrodynamic storm tide model simulations. The resulting set of maximum water surface elevations are dimensionality reduced using Karhunen–Loève expansions and then used as a training set to develop a polynomial chaos (PC) surrogate model from which global sensitivities and probabilistic predictions can be extracted. The maximum water surface elevation is extrapolated over dry points incorporating energy head loss with distance to properly train the surrogate for predicting inundation. We find that the surrogate constructed with third-order PCs using elastic net penalized regression with leave-one-out cross validation provides the most robust fit across training and test sets. Probabilistic predictions of maximum water surface elevation and inundation area by the surrogate model at 48-h lead time for three past U.S. landfalling hurricanes (Irma in 2017, Florence in 2018, and Laura in 2020) are found to be reliable when compared to best track hindcast simulation results, even when trained with as few as 19 samples. The maximum water surface elevation is most sensitive to perpendicular track-offset errors for all three storms. Laura is also highly sensitive to storm size and has the least reliable prediction. Significance Statement The purpose of this study is to develop and evaluate a methodology that can be used to provide high-quality probabilistic predictions of hurricane-induced storm tide and inundation with limited time and resources. This is important for emergency management purposes during or after the landfall of hurricanes. Our results show that sampling forecast errors using quasi-random sequences combined with machine learning techniques that fit polynomial functions to the data are well suited to this task. The polynomial functions also have the benefit of producing exact sensitivity indices of storm tide and inundation to the forecasted hurricane properties such as path, intensity, and size, which can be used for uncertainty estimation. The code implementing the presented methodology is publicly available on GitHub.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1989

A numerically exact full wave packet approach to molecule-surface scattering

A numerically exact spectral method for solving the time-dependent Schroedinger equation in spherical coordinates is described. The angular dependence of the wave function is represented on a two-dimensional grid of evenly spaced points. The fast Fourier transform algorithm is used to transform between the angle space representation of the wave function and its conjugate representation in momentum space. The time propagation of the wave function is evaluated using an expansion of the time evolution operator as a series of Chebyshev polynomials. Calculations performed for a model system representing H2 scattering from a rectangular corrugated surface yield transition probabilities that are in excellent agreement with those obtained using the close-coupling wave packet (CCWP) method. However, the new method is found to require substantially more computation time than the CCWP method because of the large number of grid points needed to represent the angular dependence of the wave function and the variation in the number of terms required in the Chebyshev representation of the time evolution operator.

Mowrey, R. C.↗
NASA NTRS1989

The effect of an unsteady three-dimensional wake on elastic blade-flapping eigenvalues in hover

This paper describes the formulation of a finite-state inflow model based on an acceleration potential and a helical wake geometry. The states of the model are coefficients of an inflow expansion in terms of a Fourier series (azimuthally) and of special polynomials (radially). The integrals over the wake are done in closed-form to obtain a set of ordinary differential equations for the inflow coefficients. The forcing functions for these equations are generalized forces which are integrals of the blade loading exactly as in structural dynamics. This model implicitly includes (for the hover case) Prandtl-Goldstein tip losses, dynamic inflow, and Theodorsen/Loewy lift deficiency. Thus, it is a fully three-dimensional unsteady wake model. Here, this model is coupled with elastic-blade equations in hover and eigenvalues are found. The results show that the three-dimensional wake has a large effect on the flap damping of all modes.

Peters, David A.↗
NASA NTRS1998

Global Nonlinear Parametric Modeling with Application to F-16 Aerodynamics

A global nonlinear parametric modeling technique is described and demonstrated. The technique uses multivariate orthogonal modeling functions generated from the data to determine nonlinear model structure, then expands each retained modeling function into an ordinary multivariate polynomial. The final model form is a finite multivariate power series expansion for the dependent variable in terms of the independent variables. Partial derivatives of the identified models can be used to assemble globally valid linear parameter varying models. The technique is demonstrated by identifying global nonlinear parametric models for nondimensional aerodynamic force and moment coefficients from a subsonic wind tunnel database for the F-16 fighter aircraft. Results show less than 10% difference between wind tunnel aerodynamic data and the nonlinear parameterized model for a simulated doublet maneuver at moderate angle of attack. Analysis indicated that the global nonlinear parametric models adequately captured the multivariate nonlinear aerodynamic functional dependence.

Morelli, Eugene A.↗
NASA NTRS1997

Global Nonlinear Parametric Modeling with Application to F-16 Aerodynamics

A global nonlinear parametric modeling technique is described and demonstrated. The technique uses multivariate orthogonal modeling functions generated from the data to determine nonlinear model structure, then expands each retained modeling function into an ordinary multivariate polynomial. The final model form is a finite multivariate power series expansion for the dependent variable in terms of the independent variables. Partial derivatives of the identified models can be used to assemble globally valid linear parameter varying models. The technique is demonstrated by identifying global nonlinear parametric models for nondimensional aerodynamic force and moment coefficients from a subsonic wind tunnel database for the F-16 fighter aircraft. Results show less than 10% difference between wind tunnel aerodynamic data and the nonlinear parameterized model for a simulated doublet maneuver at moderate angle of attack. Analysis indicated that the global nonlinear parametric models adequately captured the multivariate nonlinear aerodynamic functional dependence.

Morelli, Eugene A.↗
NASA NTRS2017

Autonomous Guidance Navigation and Control for Agile Quadrotors Using Polynomial Trajectory Planning and L1 Adaptive Control

We address the challenge to allow efficient autonomous flight in real world environments, both indoor and outdoor. We use a straight-line SE-SCP (Spherical Expansion and Sequential Convex Programming) [algorithm] to find an initial route through the environment and minimum snap trajectory generation using piecewise polynomials. Then, we implement an adaptive robust control able to address some robustness issues for quadrotors in outdoor flight, such as mass variation and wind disturbances. Coupling these techniques we allow high-speed and aggressive autonomous flight through obstacle-dense indoor environments, as well as address outdoor disturbances.

Landolfi, Mattia↗