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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Intercomparison of haloe and ER-2 aircraft H2O and CH4 observations collected during the second Airborne Arctic Stratospheric Experiment (AASE-II)

HALOE observations of H2O and CH4 are compared with in situ techniques aboard the ER-2 aircraft during the northern winter of 1991/92, in particular for the dates 911208, 920108, 920217, 920222, and 920320 when the spatial coincidences are close, within +/- 1 deg latitude and +/- 12 deg longitude. The results reveal the limitations of comparing high resolution in situ aircraft data with a remote sounding limb scanner. Of the five comparison dates, three had HALOE/ER-2 coincidences which occurred near the edge of the Arctic vortex; the vertical variability in the HALOE results and the horizontal variability in the ER-2 observations on these days show that the vortex edge is not a region where exact agreement can be expected except by chance. On a 4th comparison date, there was substantial overlap away from the vortex edge, although for some species the aircraft data show considerable variability near the coincidence point. On the 5th comparison date, the ER-2 had no overlap in altitude with the lowest HALOE observations; however, a short linear interpolation over about 1 km altitude results in smooth composite profiles. Generally speaking the agreement between HALOE and the ER-2 at overlap altitudes is about 12 percent in the case of water vapor, which shows low horizontal and vertical variability. The agreement for methane, with limited data having altitude overlap, is better than 6 percent.

Tuck, A. F.↗
NASA NTRS1993

Multi-band reflector antenna with double-ring element frequency selective subreflector

Frequency selective subreflectors (FSS) are often employed in the reflector antenna system of a communication satellite or a deep space exploration vehicle for multi-frequency operations. In the past, FSS's have been designed for diplexing two frequency bands. For example, the Voyager FSS was designed to diplex S and X bands and the TDRSS FSS was designed to diplex S and Ku bands. Recently, NASA's CASSINI project requires an FSS to multiplex four frequency (S/X/Ku/Ka) bands. Theoretical analysis and experimental verifications are presented for a multi-band flat pannel FSS with double-ring elements. Both the exact formulation and the thin-ring approximation are described for analyzing and designing this multi-ring patch element FSS. It is found that the thin-ring approximation fails to predict the electrically wide ring element FSS's performance. A single screen double-ring element FSS is demonstrated for the tri-band system that reflects the X-band signal while transmitting through the S- and Ku-band signals. In addition, a double screen FSS with non-similar double-ring elements is presented for the Cassini's four-band system which reflects the X- and Ka-band signals while passing the S- and Ku-band signals. To accurately predict the FSS effects on a dual reflector antenna's radiation pattern, the FSS subreflector's transmitted/reflected field variation as functions of the polarization and incident angles with respect to the local coordinates was taken into account. An FSS transmission/reflection coefficient table is computed for TE and TM polarizations at various incident angles based on the planar FSS model. Next, the hybrid Geometric Optics (GO) and Physical Optics (PO) technique is implemented with linearly interpolating the FSS table to efficiently determine the FSS effects in a dual reflector antenna.

Wu, Te-Kao↗
NASA NTRS1994

A Lagrangian dynamic subgrid-scale model turbulence

A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

Meneveau, C.↗
NASA NTRS1994

Recalculation of the isotropic H functions

The isotropic H functions defined in radiative transfer theory by Chandrasekhar (1960) have been numerically recalculated for a larger number of single scattering albedo (w) values. The obtained H functions should be accurate to five decimal places at least. The H function values become highly nonlinear as w approaches 1, where calculations were performed with many more points than in Chandrasekhar (1960). A simple linear interpolation of Chandrasekhar's H function table could result in as much as 8% error in the H function values and 16% in the H(exp 2) values that appear in multiple scattering terms such as in Hapke (1981). On the other hand, Hapke's approximations (1981, 1993) give H values within 4.1% and 0.8% error, respectively.

Hiroi, Takahiro↗
NASA NTRS1995

Using the application visualization system to view HALOE three-dimensional satellite data

The Application Visualization System (AVS) is used to view a three-dimensional data field containing the volume mixing ratios of a chemical species in the middle atmosphere obtained by the Halogen Occultation Experiment (HALOE) aboard the Upper Atmosphere Research Satellite (UARS). Since launch in September 1991, HALOE has been collecting data on approximately 30 sunrise/sunset events in two narrow latitude bands each day. The vertical volume mixing ratio profiles are retrieved for eight species for each event. The accumulated data for approximately 30 days cover most of the globe (limited by sunlit latitudes), and this monthly data block can be described as the volume mixing ratio of a specific species in the atmosphere as a function of latitude, longitude, and height. The data were remapped using linear interpolation for pressure levels and Gaussian weighted binning from sampling locations to a three-dimensional grid. An AVS network is constructed that allows for viewing the three-dimensional field with rendered slices at constant latitudes, longitudes or pressure levels. Discussions are given on the advantages and some disadvantages learned about from experiences applying AVS to visualize HALOE three dimensional data.

Luo, Mingzhao↗
NASA NTRS1996

Variance Analysis of Unevenly Spaced Time Series Data

We have investigated the effect of uneven data spacing on the computation of delta (sub chi)(gamma). Evenly spaced simulated data sets were generated for noise processes ranging from white phase modulation (PM) to random walk frequency modulation (FM). Delta(sub chi)(gamma) was then calculated for each noise type. Data were subsequently removed from each simulated data set using typical two-way satellite time and frequency transfer (TWSTFT) data patterns to create two unevenly spaced sets with average intervals of 2.8 and 3.6 days. Delta(sub chi)(gamma) was then calculated for each sparse data set using two different approaches. First the missing data points were replaced by linear interpolation and delta (sub chi)(gamma) calculated from this now full data set. The second approach ignored the fact that the data were unevenly spaced and calculated delta(sub chi)(gamma) as if the data were equally spaced with average spacing of 2.8 or 3.6 days. Both approaches have advantages and disadvantages, and techniques are presented for correcting errors caused by uneven data spacing in typical TWSTFT data sets.

Hackman, Christine↗
NASA NTRS1998

Conversion from Engineering Units to Telemetry Counts on Dryden Flight Simulators

Dryden real-time flight simulators encompass the simulation of pulse code modulation (PCM) telemetry signals. This paper presents a new method whereby the calibration polynomial (from first to sixth order), representing the conversion from counts to engineering units (EU), is numerically inverted in real time. The result is less than one-count error for valid EU inputs. The Newton-Raphson method is used to numerically invert the polynomial. A reverse linear interpolation between the EU limits is used to obtain an initial value for the desired telemetry count. The method presented here is not new. What is new is how classical numerical techniques are optimized to take advantage of modem computer power to perform the desired calculations in real time. This technique makes the method simple to understand and implement. There are no interpolation tables to store in memory as in traditional methods. The NASA F-15 simulation converts and transmits over 1000 parameters at 80 times/sec. This paper presents algorithm development, FORTRAN code, and performance results.

Fantini, Jay A.↗
NASA NTRS1999

Frequency and Angular Variations of Land Surface Microwave Emissivities: Can we Estimate SSM/T and AMSU Emissivities from SSM/I Emissivities?

To retrieve temperature and humidity profiles from SSM/T and AMSU, it is important to quantify the contribution of the Earth surface emission. So far, no global estimates of the land surface emissivities are available at SSM/T and AMSU frequencies and scanning conditions. The land surface emissivities have been previously calculated for the globe from the SSM/I conical scanner between 19 and 85 GHz. To analyze the feasibility of deriving SSM/T and AMSU land surface emissivities from SSM/I emissivities, the spectral and angular variations of the emissivities are studied, with the help of ground-based measurements, models and satellite estimates. Up to 100 GHz, for snow and ice free areas, the SSM/T and AMSU emissivities can be derived with useful accuracy from the SSM/I emissivities- The emissivities can be linearly interpolated in frequency. Based on ground-based emissivity measurements of various surface types, a simple model is proposed to estimate SSM/T and AMSU emissivities for all zenith angles knowing only the emissivities for the vertical and horizontal polarizations at 53 deg zenith angle. The method is tested on the SSM/T-2 91.655 GHz channels. The mean difference between the SSM/T-2 and SSM/I-derived emissivities is less than or equal to 0.01 for all zenith angles with an r.m.s. difference of approx. = 0.02. Above 100 GHz, preliminary results are presented at 150 GHz, based on SSM/T-2 observations and are compared with the very few estimations available in the literature.

Prigent, Catherine↗
NASA NTRS2000

On Fast Post-Processing of Global Positioning System Simulator Truth Data and Receiver Measurements and Solutions Data

Post-Processing of data related to a Global Positioning System (GPS) simulation is an important activity in qualification of a GPS receiver for space flight. Because a GPS simulator is a critical resource it is desirable to move off the pertinent simulation data from the simulator as soon as a test is completed. The simulator data files are usually moved to a Personal Computer (PC), where the post-processing of the receiver logged measurements and solutions data and simulated data is performed. Typically post-processing is accomplished using PC-based commercial software languages and tools. Because of commercial software systems generality their general-purpose functions are notoriously slow and more than often are the bottleneck problem even for short duration experiments. For example, it may take 8 hours to post-process data from a 6-hour simulation. There is a need to do post-processing faster, especially in order to use the previous test results as feedback for a next simulation setup. This paper demonstrates that a fast software linear interpolation algorithm is applicable to a large class of engineering problems, like GPS simulation data post-processing, where computational time is a critical resource and is one of the most important considerations. An approach is developed that allows to speed-up post-processing by an order of magnitude. It is based on improving the post-processing bottleneck interpolation algorithm using apriori information that is specific to the GPS simulation application. The presented post-processing scheme was used in support of a few successful space flight missions carrying GPS receivers. A future approach to solving the post-processing performance problem using Field Programmable Gate Array (FPGA) technology is described.

Kizhner, Semion↗
NASA NTRS2000

TAILSIM Users Guide

The TAILSIM program uses a 4th order Runge-Kutta method to integrate the standard aircraft equations-of-motion (EOM). The EOM determine three translational and three rotational accelerations about the aircraft's body axis reference system. The forces and moments that drive the EOM are determined from aerodynamic coefficients, dynamic derivatives, and control inputs. Values for these terms are determined from linear interpolation of tables that are a function of parameters such as angle-of-attack and surface deflections. Buildup equations combine these terms and dimensionalize them to generate the driving total forces and moments. Features that make TAILSIM applicable to studies of tailplane stall include modeling of the reversible control System, modeling of the pilot performing a load factor and/or airspeed command task, and modeling of vertical gusts. The reversible control system dynamics can be described as two hinged masses connected by a spring. resulting in a fifth order system. The pilot model is a standard form of lead-lag with a time delay applied to an integrated pitch rate and/or airspeed error feedback. The time delay is implemented by a Pade approximation, while the commanded pitch rate is determined by a commanded load factor. Vertical gust inputs include a single 1-cosine gust and a continuous NASA Dryden gust model. These dynamic models. coupled with the use of a nonlinear database, allow the tailplane stall characteristics, elevator response, and resulting aircraft response, to be modeled. A useful output capability of the TAILSIM program is the ability to display multiple post-run plot pages to allow a quick assessment of the time history response. There are 16 plot pages currently available to the user. Each plot page displays 9 parameters. Each parameter can also be displayed individually. on a one plot-per-page format. For a more refined display of the results the program can also create files of tabulated data. which can then be used by other plotting programs. The TAILSIM program was written straightforwardly assuming the user would want to change the database tables, the buildup equations, the output parameters. and the pilot model parameters. A separate database file and input file are automatically read in by the program. The use of an include file to set up all common blocks facilitates easy changing of parameter names and array sizes.

Hiltner, Dale W.↗
NASA NTRS2002

Origin of the Bandgap Bowing in ZNSe(1-x)Te(x) Alloys

Photomodulated reflection and optical absorption spectroscopies were used to measure the composition and hydrostatic pressure dependencies of interband optical transitions in ZnSe(1-x)Te(x) alloys. On the Te-rich side, the bandgap reduction and the reduction in the pressure dependence of the bandgap are well explained by the band anticrossing interaction between the ZnTe conduction band states and the Se localized state located at 2.85 eV above the valence band top of ZnTe. On the Se-rich side, an interaction between the degenerate Gamma valence bands of ZnSe and the localized Te state located at 0.1 eV above the valence band top of ZnSe is proposed to generate the bandgap reduction and the rapid increase of the spin-orbit splitting with increasing Te concentration. A linear interpolation between these two models accounts for the composition and pressure dependencies of the bandgap in the entire composition range.

Walukiewicz, W.↗
NASA NTRS2002

Origin of the Large Bandgap Bowing in Highly Mismatched Semiconductor Alloys

We use a combination of optical and soft x-ray fluorescence techniques to investigate the composition dependence of the band edge energies in ZnSe(1-x)Te(x) and ZnS(1-x)Te(x) alloys. The results reveal entirely different origins of the large bandgap bowing for alloys with small and large Te content. On the Te-rich side, the reduction of the bandgap is well explained by the band anticrossing interaction between the Se or S localized states and the ZnTe conduction band states. On the Se or S-rich side, an interaction between the localized Te states and the degenerate GAMMA valence bands of ZnSe or ZnS is responsible for the bandgap reduction and the rapid increase of the spin-orbit splitting with increasing Te concentration. Results of soft x-ray fluorescence experiments provide direct evidence for the valence band anticrossing interaction. The bandgap bowing over the entire composition range is accounted for by a linear interpolation between the conduction band anticrossing and valence band anticrossing models.

Wu, J.↗
NASA NTRS2001

Global Surface Solar Energy Anomalies Including El Nino and La Nina Years

This paper synthesizes past events in an attempt to define the general magnitude, duration, and location of large surface solar anomalies over the globe. Surface solar energy values are mostly a function of solar zenith angle, cloud conditions, column atmospheric water vapor, aerosols, and surface albedo. For this study, solar and meteorological parameters for the 10-yr period July 1983 through June 1993 are used. These data were generated as part of the Release 3 Surface meteorology and Solar Energy (SSE) activity under the NASA Earth Science Enterprise (ESE) effort. Release 3 SSE uses upgraded input data and methods relative to previous releases. Cloud conditions are based on recent NASA Version-D International Satellite Cloud Climatology Project (ISCCP) global satellite radiation and cloud data. Meteorological inputs are from Version-I Goddard Earth Observing System (GEOS) reanalysis data that uses both weather station and satellite information. Aerosol transmission for different regions and seasons are for an 'average' year based on historic solar energy data from over 1000 ground sites courtesy of Natural Resources Canada (NRCan). These data are input to a new Langley Parameterized Shortwave Algorithm (LPSA) that calculates surface albedo and surface solar energy. That algorithm is an upgraded version of the 'Staylor' algorithm. Calculations are performed for a 280X280 km equal-area grid system over the globe based on 3-hourly input data. A bi-linear interpolation process is used to estimate data output values on a 1 X 1 degree grid system over the globe. Maximum anomalies are examined relative to El Nino and La Nina events in the tropical Pacific Ocean. Maximum year-to-year anomalies over the globe are provided for a 10-year period. The data may assist in the design of systems with increased reliability. It may also allow for better planning for emergency assistance during some atypical events.

Whitlock, C. H.↗
NASA NTRS2002

Interface Conditions for Wave Propagation Through Mesh Refinement Boundaries

We study the propagation of waves across fixed mesh refinement boundaries in linear and nonlinear model equations in 1-D and 2-D, and in the 3-D Einstein equations of general relativity. We demonstrate that using linear interpolation to set the data in guard cells leads to the production of reflected waves at the refinement boundaries. Implementing quadratic interpolation to fill the guard cells eliminates these spurious signals.

Choi, Dae-II↗
NASA NTRS2003

Composition Dependence of the Hydrostatic Pressure Coefficients of the Bandgap of ZnSe(1-x)Te(x) Alloys

Optical absorption experiments have been performed using diamond anvil cells to measure the hydrostatic pressure dependence of the fundamental bandgap of ZnSe(sub 1-xTe(sub x) alloys over the entire composition range. The first and second-order pressure coefficients are obtained as a function of composition. Starting from the ZnSe side, the magnitude of both coefficients increases slowly until x approx. 0.7, where the ambient-pressure bandgap reaches a minimum. For larger values of x the coefficients rapidly approach the values of ZnTe. The large deviations of the pressure coefficients from the linear interpolation between ZnSe and ZnTe are explained in terms of the band anticrossing model.

Wu, J.↗
NASA NTRS1988

Variability in surface ECG morphology: signal or noise?

Using data collected from canine models of acute myocardial ischemia, we investigated two issues of major relevance to electrocardiographic signal averaging: ECG epoch alignment, and the spectral characteristics of the beat-to-beat variability in ECG morphology. With initial digitization rates of 1 kHz, an iterative a posteriori matched filtering alignment scheme, and linear interpolation, we demonstrated that there is sufficient information in the body surface ECG to merit alignment to a precision of 0.1 msecs. Applying this technique to align QRS complexes and atrial pacing artifacts independently, we demonstrated that the conduction delay from atrial stimulus to ventricular activation may be so variable as to preclude using atrial pacing as an alignment mechanism, and that this variability in conduction time be modulated at the frequency of respiration and at a much lower frequency (0.02-0.03Hz). Using a multidimensional spectral technique, we investigated the beat-to-beat variability in ECG morphology, demonstrating that the frequency spectrum of ECG morphological variation reveals a readily discernable modulation at the frequency of respiration. In addition, this technique detects a subtle beat-to-beat alternation in surface ECG morphology which accompanies transient coronary artery occlusion. We conclude that physiologically important information may be stored in the variability in the surface electrocardiogram, and that this information is lost by conventional averaging techniques.

NASA Discipline Cardiopulmonary↗
NASA NTRS2006

3-D Mixed Mode Delamination Fracture Criteria - An Experimentalist's Perspective

Many delamination failure criteria based on fracture toughness have been suggested over the past few decades, but most only covered the region containing mode I and mode II components of loading because that is where toughness data existed. With new analysis tools, more 3D analyses are being conducted that capture a mode III component of loading. This has increased the need for a fracture criterion that incorporates mode III loading. The introduction of a pure mode III fracture toughness test has also produced data on which to base a full 3D fracture criterion. In this paper, a new framework for visualizing 3D fracture criteria is introduced. The common 2D power law fracture criterion was evaluated to produce unexpected predictions with the introduction of mode III and did not perform well in the critical high mode I region. Another 2D criterion that has been shown to model a wide range of materials well was used as the basis for a new 3D criterion. The new criterion is based on assumptions that the relationship between mode I and mode III toughness is similar to the relation between mode I and mode II and that a linear interpolation can be used between mode II and mode III. Until mixed-mode data exists with a mode III component of loading, 3D fracture criteria cannot be properly evaluated, but these assumptions seem reasonable.

Reeder, James R.↗
NASA NTRS2010

Characteristics of Whipple Shield Performance in the Shatter Regime

Ballistic limit equations define the failure of metallic Whipple shields in three parts: low velocity, shatter, and hypervelocity. Failure limits in the shatter regime are based on a linear interpolation between the onset of projectile fragmentation, and impulsive rupture of the shield rear wall. A series of hypervelocity impact tests have been performed on aluminum alloy Whipple shields to investigate failure mechanisms and performance limits in the shatter regime. Test results demonstrated a more rapid increase in performance than predicted by the latest iteration of the JSC Whipple shield BLE following the onset of projectile fragmentation. This increase in performance was found to level out between 4.0-5.0 km/s, with a subsequent decrease in performance for velocities up to 6.0 km/s. For a detached spall failure criterion, the failure limit was found to continually decrease up to a velocity of 7.0 km/s, substantially varying from the BLE, while for perforation-based failure an increase in performance was observed. An existing phenomenological ballistic limit curve was found to provide a more accurate reproduction of shield behavior that the BLE, however a number of underlying assumptions such as the occurrence of complete projectile fragmentation and the effect on performance of incipient projectile melt were found to be inaccurate. A cratering relationship based on the largest residual fragment size has been derived for application at velocities between 3.0-4.0 km/s, and was shown to accurately reproduce the trends of the experimental data. Further investigation is required to allow a full analytical description of shatter regime performance for metallic Whipple shields.

Ryan, S.↗