Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear interpolation”

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 37 records · Page 2

Frequency-shaped large-angle maneuvers

The paper considers the problem of maneuvering a flexible spacecraft through large angles in finite time. The basic control problem is divided into two parts. The first part consists of generating a frequency-shaped open-loop solution for the nonlinear rigid body as the nominal solution. The resulting two-point boundary-value problem is solved by introducing a continuation method for altering the mass distribution and boundary conditions for the spacecraft. For the second part, a feedback control is designed by linearizing the flexible body response about several points along the rigid body nominal solution. The perturbation gains are designed by using a frequency-shaped cost functional approach. The gains are linearly interpolated to produce smooth control time-histories as the linear piecewise constant plant models change during the maneuver.

Chun, Hon M.↗

MODIS Solar Diffuser Degradation Determination and Its Spectral Dependency

This study presents a modeling approach to improve solar diffuser (SD) degradation determination from SD stability monitor (SDSM) measurements. The MODIS instrument uses a SDto calibrate its reflective solar bands (RSBs) on-orbit. Due to the imperfectly designed SDSM sun view screen, the SD reflectance tracked by SDSM has large noise. The SDSM measurements noise is spectrally coherent and can be minimized by normalizing measurements to the least degradeddetector 9 (936 nm). In this study, a SD degradation model is used to determine the SDdegradation's wavelength dependency and the detector 9 degradation is estimated by the model solution.The results show the SD degradations measured at 6 SDSM detectors (554 _ 936 nm) have stable relationships, where the degradation is inversely proportion to 1/wavelength^4. The model estimated SD degradation at SDSM detector 9 wavelength (936 nm) is ~0.9% from 2002 to 2018.Based on the SD degradation model solution, the SD degradation at short/mid wave bands are estimated to improve short/mid wave bands calibration. The model can also be used to improve interpolating SD degradation at SDSM detectors to RSB wavelengths. Compared to linear interpolation, bands 9 and 10 show the largest differences of up to 0.3 and 0.4% respectively. These differences directly impact the calibration coefficients of these bands.

radiometric calibration↗

Secondary Sonic Boom Ray Path Prediction using sBOOM

Sonic boom ray path equations are solved numerically to show the extent of secondary booms compared to primary sonic boom carpet. Secondary booms are shown to reach altitudes of ~125 km before turning toward the ground and impinging at lateral distances of ~50 km and beyond. Exact closed-form solutions of the ray paths and their partial derivatives with respect to certain relevant parameters are derived for the simplified case of windless atmospheres and compared against numerical solutions. Ray tube areas are computed using the Jacobian of a single ray as opposed to the traditional technique of using 4 rays that are temporally and azimuthally perturbed. The ray tube areas corresponding to secondary booms can get 1000 times larger than those from primary booms. The numerically computed ray tube areas are verified using exact expressions for the special case of windless atmospheres. Finally, pressure scaling due to ray tube area changes and atmospheric stratification is presented and the shortcomings associated with the current atmospheric interpolation techniques are discussed. It is determined that piece-wise linear interpolation of the atmospheric temperature may be inadequate and could lead to numerical noise in the prediction of secondary booms.

Sriram K Rallabhandi↗

Input-output characterization of an ultrasonic testing system by digital signal analysis

Ultrasonic test system input-output characteristics were investigated by directly coupling the transmitting and receiving transducers face to face without a test specimen. Some of the fundamentals of digital signal processing were summarized. Input and output signals were digitized by using a digital oscilloscope, and the digitized data were processed in a microcomputer by using digital signal-processing techniques. The continuous-time test system was modeled as a discrete-time, linear, shift-invariant system. In estimating the unit-sample response and frequency response of the discrete-time system, it was necessary to use digital filtering to remove low-amplitude noise, which interfered with deconvolution calculations. A digital bandpass filter constructed with the assistance of a Blackman window and a rectangular time window were used. Approximations of the impulse response and the frequency response of the continuous-time test system were obtained by linearly interpolating the defining points of the unit-sample response and the frequency response of the discrete-time system. The test system behaved as a linear-phase bandpass filter in the frequency range 0.6 to 2.3 MHz. These frequencies were selected in accordance with the criterion that they were 6 dB below the maximum peak of the amplitude of the frequency response. The output of the system to various inputs was predicted and the results were compared with the corresponding measurements on the system.

Williams, J. H., Jr.↗

Triangle based TVD schemes for hyperbolic conservation laws

A triangle based total variation diminishing (TVD) scheme for the numerical approximation of hyperbolic conservation laws in two space dimensions is constructed. The novelty of the scheme lies in the nature of the preprocessing of the cell averaged data, which is accomplished via a nearest neighbor linear interpolation followed by a slope limiting procedures. Two such limiting procedures are suggested. The resulting method is considerably more simple than other triangle based non-oscillatory approximations which, like this scheme, approximate the flux up to second order accuracy. Numerical results for linear advection and Burgers' equation are presented.

Durlofsky, Louis J.↗

Triangle based adaptive stencils for the solution of hyperbolic conservation laws

A triangle based total variation diminishing (TVD) scheme for the numerical approximation of hyperbolic conservation laws in two space dimensions is constructed. The novelty of the scheme lies in the nature of the preprocessing of the cell averaged data, which is accomplished via a nearest neighbor linear interpolation followed by a slope limiting procedures. Two such limiting procedures are suggested. The resulting method is considerably more simple than other triangle based non-oscillatory approximations which, like this scheme, approximate the flux up to second order accuracy. Numerical results for linear advection and Burgers' equation are presented.

Durlofsky, Louis J.↗

A New Machine Learning Based Analysis for Improving Satellite Retrieved Atmospheric Composition Data: OMI SO2 as an Example

Despite recent progress, satellite retrievals of anthropogenic SO2 still suffer from relatively low signal-tonoise ratios. In this study, we demonstrate a new machine learning data analysis method to improve the quality of satellite SO2 products. In the absence of large ground-truth datasets for SO2, we start from SO2 slant column densities (SCDs) retrieved from the Ozone Monitoring Instrument (OMI) using a data-driven, physically based algorithm and calculate the ratio between the SCD and the root mean square (rms) of the fitting residuals for each pixel. To build the training data, we select presumably clean pixels with small SCD / rms ratios (SRRs) and set their target SCDs to zero. For polluted pixels with relatively large SRRs, we set the target to the original retrieved SCDs. We then train neural networks (NNs) to reproduce the target SCDs using predictors including SRRs for individual pixels, solar zenith, viewing zenith and phase angles, scene reflectivity, and O3 column amounts, as well as the monthly mean SRRs. For data analysis, we employ two NNs: (1) one trained daily to produce analyzed SO2 SCDs for polluted pixels each day and (2) the other trained once every month to produce analyzed SCDs for less polluted pixels for the entire month. Test results for 2005 show that our method can significantly reduce noise and artifacts over background regions. Over polluted areas, the monthly mean NN-analyzed and original SCDs generally agree to within ±15 %, indicating that our method can retain SO2 signals in the original retrievals except for large volcanic eruptions. This is further confirmed by running both the NN-analyzed and original SCDs through a topdown emission algorithm to estimate the annual SO2 emissions for ∼ 500 anthropogenic sources, with the two datasets yielding similar results. We also explore two alternative approaches to the NN-based analysis method. In one, we employ a simple linear interpolation model to analyze the original SCD retrievals. In the other, we develop a PCA–NN algorithm that uses OMI measured radiances, transformed and dimension-reduced with a principal component analysis (PCA) technique, as inputs to NNs for SO2 SCD retrievals. While the linear model and the PCA–NN algorithm can reduce retrieval noise, they both underestimate SO2 over polluted areas. Overall, the results presented here demonstrate that our new data analysis method can significantly improve the quality of existing OMI SO2 retrievals. The method can potentially be adapted for other sensors and/or species and enhance the value of satellite data in air quality research and applications.

Can Li↗

Polynomial Expressions for Estimating Elastic Constants From the Resonance of Circular Plates

Two approaches were taken to make convenient spread sheet calculations of elastic constants from resonance data and the tables in ASTM C1259 and E1876: polynomials were fit to the tables; and an automated spread sheet interpolation routine was generated. To compare the approaches, the resonant frequencies of circular plates made of glass, hardened maraging steel, alpha silicon carbide, silicon nitride, tungsten carbide, tape cast NiO-YSZ, and zinc selenide were measured. The elastic constants, as calculated via the polynomials and linear interpolation of the tabular data in ASTM C1259 and E1876, were found comparable for engineering purposes, with the differences typically being less than 0.5 percent. Calculation of additional v values at t/R between 0 and 0.2 would allow better curve fits. This is not necessary for common engineering purposes, however, it might benefit the testing of emerging thin structures such as fuel cell electrolytes, gas conversion membranes, and coatings when Poisson s ratio is less than 0.15 and high precision is needed.

Salem, Jonathan A.↗

Spectral Analysis Of Linear, Shift-Invariant Interpolants

Method of analysis provides quantitative measure of reconstruction and interpolation performances of linear, shift-invariant interpolants. Criterion of performance based upon mean-square error of difference between sampled and reconstructed functions. Applicable to reconstruction algorithms used in processing of signals and images and to types of interpolants used in numerical analysis, computer-aided design, and computer graphics.

Lansing, Donald L.↗

Fully energy-dependent HZETRN (a galactic cosmic-ray transport code)

For extended manned space missions, the radiation shielding design requires efficient and accurate cosmic-ray transport codes that can handle the physics processes in detail. The Langley Research Center galactic cosmic-ray transport code (HZETRN) is currently under development for such design use. The cross sections for the production of secondary nucleons in the existing HZETRN code are energy dependent only for nucleon collisions. The approximation of energy-independent, heavy-ion fragmentation cross section is now removed by implementing a mathematically simplified energy-dependent stepping formalism for heavy ions. The cross section at each computational grid is obtained by linear interpolation from a few tabulated data to minimize computing time. Test runs were made for galactic cosmic-ray transport through a liquid hydrogen shield and a water shield at solar minimum. The results show no appreciable change in total fluxes or computing time compared with energy-independent calculations. Differences in high LET (linear energy transfer) spectra are noted, however, because of the large variation in cross sections at the low-energy region. The high LET components are significantly higher in the new code and have important implications on biological risk estimates for heavy-ion exposure.

Shinn, Judy L.↗

GMI-IPS: Python Processing Software for Aircraft Campaigns

NASA's Atmospheric Tomography Mission (ATom) seeks to understand the impact of anthropogenic air pollution on gases in the Earth's atmosphere. Four flight campaigns are being deployed on a seasonal basis to establish a continuous global-scale data set intended to improve the representation of chemically reactive gases in global atmospheric chemistry models. The Global Modeling Initiative (GMI), is creating chemical transport simulations on a global scale for each of the ATom flight campaigns. To meet the computational demands required to translate the GMI simulation data to grids associated with the flights from the ATom campaigns, the GMI ICARTT Processing Software (GMI-IPS) has been developed and is providing key functionality for data processing and analysis in this ongoing effort. The GMI-IPS is written in Python and provides computational kernels for data interpolation and visualization tasks on GMI simulation data. A key feature of the GMI-IPS, is its ability to read ICARTT files, a text-based file format for airborne instrument data, and extract the required flight information that defines regional and temporal grid parameters associated with an ATom flight. Perhaps most importantly, the GMI-IPS creates ICARTT files containing GMI simulated data, which are used in collaboration with ATom instrument teams and other modeling groups. The initial main task of the GMI-IPS is to interpolate GMI model data to the finer temporal resolution (1-10 seconds) of a given flight. The model data includes basic fields such as temperature and pressure, but the main focus of this effort is to provide species concentrations of chemical gases for ATom flights. The software, which uses parallel computation techniques for data intensive tasks, linearly interpolates each of the model fields to the time resolution of the flight. The temporally interpolated data is then saved to disk, and is used to create additional derived quantities. In order to translate the GMI model data to the spatial grid of the flight path as defined by the pressure, latitude, and longitude points at each flight time record, a weighted average is then calculated from the nearest neighbors in two dimensions (latitude, longitude). Using SciPya's Regular Grid Interpolator, interpolation functions are generated for the GMI model grid and the calculated weighted averages. The flight path points are then extracted from the ATom ICARTT instrument file, and are sent to the multi-dimensional interpolating functions to generate GMI field quantities along the spatial path of the flight. The interpolated field quantities are then written to a ICARTT data file, which is stored for further manipulation. The GMI-IPS is aware of a generic ATom ICARTT header format, containing basic information for all flight campaigns. The GMI-IPS includes logic to edit metadata for the derived field quantities, as well as modify the generic header data such as processing dates and associated instrument files. The ICARTT interpolated data is then appended to the modified header data, and the ICARTT processing is complete for the given flight and ready for collaboration. The output ICARTT data adheres to the ICARTT file format standards V1.1. The visualization component of the GMI-IPS uses Matplotlib extensively and has several functions ranging in complexity. First, it creates a model background curtain for the flight (time versus model eta levels) with the interpolated flight data superimposed on the curtain. Secondly, it creates a time-series plot of the interpolated flight data. Lastly, the visualization component creates averaged 2D model slices (longitude versus latitude) with overlaid flight track circles at key pressure levels. The GMI-IPS consists of a handful of classes and supporting functionality that have been generalized to be compatible with any ICARTT file that adheres to the base class definition. The base class represents a generic ICARTT entry, only defining a single time entry and 3D spatial positioning parameters. Other classes inherit from this base class; several classes for input ICARTT instrument files, which contain the necessary flight positioning information as a basis for data processing, as well as other classes for output ICARTT files, which contain the interpolated model data. Utility classes provide functionality for routine procedures such as: comparing field names among ICARTT files, reading ICARTT entries from a data file and storing them in data structures, and returning a reduced spatial grid based on a collection of ICARTT entries. Although the GMI-IPS is compatible with GMI model data, it can be adapted with reasonable effort for any simulation that creates Hierarchical Data Format (HDF) files. The same can be said of its adaptability to ICARTT files outside of the context of the ATom mission. The GMI-IPS contains just under 30,000 lines of code, eight classes, and a dozen drivers and utility programs. It is maintained with GIT source code management and has been used to deliver processed GMI model data for the ATom campaigns that have taken place to date.

Damon, M. R.↗

Optical Flow for Intermediate Frame Interpolation of Multispectral Geostationary Satellite Data

Applications in areas such as weather tracking and modeling, ecosystem monitoring, wildfire detection, and land-cover change are heavily dependent on spatial and temporal resolutions of satellite observations. However, there are typically trade-offs between spatial and temporal resolutions in dataset selection. For instance, geostationary weather tracking satellites are designed to take snapshots many times throughout the day but sensor hardware limits data collection. In this work we tackle this limitation, developing a method for temporal upsampling of multi-spectral satellite imagery using optical flow video interpolation deep convolutional neural networks. The presented model, extends Super SloMo (SSM) from single optical flow estimates to multichannel where flows are computed per band. We apply this technique on 8 multi-spectral bands of NOAA/NASA's GOES-16 mesoscale dataset to temporally enhance full disk hemispheric snapshots from 15 minutes to 1 minute. Through extensive experimentation, we show SSM vastly outperforms the linear interpolation baseline and that multichannel optical flows improves performance on GOES-16. A visual analysis of optical flow vectors clearly identifies hurricanes and large-scale atmospheric dynamics. Furthermore, we discuss challenges and open questions related to optical flow and temporal interpolation of multispectral geostationary satellite imagery.

Optical Flow↗

Computer programs for thermodynamic and transport properties of hydrogen (tabcode-II)

The thermodynamic and transport properties of para and equilibrium hydrogen have been programmed into a series of computer routines. Input variables are the pair's pressure-temperature and pressure-enthalpy. The programs cover the range from 1 to 5000 psia with temperatures from the triple point to 6000 R or enthalpies from minus 130 BTU/lb to 25,000 BTU/lb. Output variables are enthalpy or temperature, density, entropy, thermal conductivity, viscosity, at constant volume, the heat capacity ratio, and a heat transfer parameter. Property values on the liquid and vapor boundaries are conveniently obtained through two small routines. The programs achieve high speed by using linear interpolation in a grid of precomputed points which define the surface of the property returned.

Roder, H. M.↗

An interative solution of an integral equation for radiative transfer by using variational technique

An effective iterative technique is introduced to solve a nonlinear integral equation frequently associated with radiative transfer problems. The problem is formulated in such a way that each step of an iterative sequence requires the solution of a linear integral equation. The advantage of a previously introduced variational technique which utilizes a stepwise constant trial function is exploited to cope with the nonlinear problem. The method is simple and straightforward. Rapid convergence is obtained by employing a linear interpolation of the iterative solutions. Using absorption coefficients of the Milne-Eddington type, which are applicable to some planetary atmospheric radiation problems. Solutions are found in terms of temperature and radiative flux. These solutions are presented numerically and show excellent agreement with other numerical solutions.

Yoshikawa, K. K.↗

Recent developments at JPL in the application of digital image processing techniques to astronomical images

Several techniques of a digital image-processing nature are illustrated which have proved useful in visual analysis of astronomical pictorial data. Processed digital scans of photographic plates of Stephans Quintet and NGC 4151 are used as examples to show how faint nebulosity is enhanced by high-pass filtering, how foreground stars are suppressed by linear interpolation, and how relative color differences between two images recorded on plates with different spectral sensitivities can be revealed by generating ratio images. Analyses are outlined which are intended to compensate partially for the blurring effects of the atmosphere on images of Stephans Quintet and to obtain more detailed information about Saturn's ring structure from low- and high-resolution scans of the planet and its ring system. The employment of a correlation picture to determine the tilt angle of an average spectral line in a low-quality spectrum is demonstrated for a section of the spectrum of Uranus.

Lorre, J. J.↗

Comment on 'Drag coefficient of spheres in continuum and rarefied flows'

A paper by Henderson (1976) provides a method of predicting experimental sphere drag data. This approach uses two equations for the drag coefficient, one for relative Mach number less than one, one for relative Mach number greater than 1.75. For relative Mach numbers between these limits a linear interpolation procedure is followed. In a comment on this paper, it is claimed, on the basis of comparing predictions with experimental results, that a method proposed by Walsh (1975) gives better predictions of the drag coefficient for relative Mach numbers less than 1.75, provided that a modification of the procedure is made for relative Mach numbers less than 0.1. For values over 1.75, both methods are considered equally accurate. In a reply to this comment, it is agreed that the Walsh method is more accurate when Reynolds numbers are within a range between 20 and 200, and Mach numbers are between 0.5 and 1.25. Presumed errors and possible limitations in the Walsh procedure for predicting drag coefficients are discussed.

Walsh, M. J.↗