Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 415 records · Page 23

Algebraic grid generation

The numerical solution of partial differential equations about irregular geometries and with varying characteristic scales has created the need for coordinate systems and associated transformations which reflect both geometric and physical requirements. The process of finding coordinate transformations in discrete representations is called 'grid generation'. The present investigation is concerned with three algebraic grid generation methods. The methods include transfinite interpolation, the multisurface method, and the two-boundary technique. Interpolation formulas, in terms of homotopic mappings and constraints in terms of point positions and/or derivatives are the essential elements of the techniques. The methods are relatively simple to understand, they are explicit and do not require extensive computational effort, and they have a high degree of generality.

Smith, R. E.↗

A climatically-derived global soil moisture data set for use in the GLAS atmospheric circulation model seasonal cycle experiment

Algorithms for point interpolation and contouring on the surface of the sphere and in Cartesian two-space are developed from Shepard's (1968) well-known, local search method. These mapping procedures then are used to investigate the errors which appear on small-scale climate maps as a result of the all-too-common practice of of interpolating, from irregularly spaced data points to the nodes of a regular lattice, and contouring Cartesian two-space. Using mean annual air temperatures field over the western half of the northern hemisphere is estimated both on the sphere, assumed to be correct, and in Cartesian two-space. When the spherically- and Cartesian-approximted air temperature fields are mapped and compared, the magnitudes (as large as 5 C to 10 C) and distribution of the errors associated with the latter approach become apparent.

Willmott, C. J.↗

First Results of GPS Time Transfer to Australia

A global positioning system (GPS) time transfer unit was installed at Tidbinbilla Deep Space Communications Complex of the DSN in June 1983. It was used to estimate the relationship to UTC(USNO MC) of the Tidbinbilla frequency and time system TID(FTS) based on a hydrogen maser, and to estimate the performance of the Australian free-running time scale UTC(AUS). Data from the first 3 months were analyzed three ways: by two-hop common view using JPL as intermediary; by long-arc interpolation of measurements against space vehicle clocks; and by long arc interpolation of GPS-Time results. Residuals from a single quadratic fit through 3 months of UTC(USNO MC) -TID(FTS) results were white noise with standard error 15 ns, and a flying clock measurement gave 70 ns agreement. A straight line fit through results UTC(USNO MC) - UTC(AUS) gave 90 ns standard error and 120 ns agreement. It is proposed to use the GPS measurements to steer UTC(AUS) to UTC(BIH), and to rename the existing time scale TA(AUS).

Mck.luck, J.↗

Observing system simulation experiments related to space-borne Lidar wind profiling. Part 1: Forecast impacts of highly idealized observing systems

Simulation experiments comparing the relative importance of an idealized LIDAR wind profiling system with idealized temperature and pressure sounding systems on 12 h forecasts are studied for three "nature' fields representing the true evolving atmospheric states. The three fields are obtained respectively from: (1) a long integration of the GLAS 4th Order Model (4 deg lat x 5 deg lon x 9 levels), (2) a continuous sequence of NMC operational analysis and, (3) a long integration of the ECMWF high resolution (1.875 deg 1 lat x 1.875 deg lon x 15 layers) operational forecast model. These fields are interpolated to the grid of the GLAS model and used for simulating the observed global analysed fields of winds, temperature, moisture and surface pressure. The same interpolated fields are also used for verification of forecast impact. The effects of clouds, aerosol concentrations, and instrument accuracies on the simulated observation will be discussed.

Halem, M.↗

Postprocessing for the UCLA climate model

Postprocessing programs written to analyze the results of a general circulation model are described. The model history data is interpolated to pressure surfaces, converted from 64-bit Cyber words to 32 bit IBM words, and transferred to tape for further analysis. The time average program reads the interpolation output tape and produces 15-day time averages. The model outputs a very detailed record for selected grid points. This data is sent to a tape in Cyber full words. The conversion program changes the format of the tape from 64 bit Cyber words to 32 bit IBM words, processing a month of data at a time. The grid-point history sorting program inputs the converted grid point history tape and sorts it by time. The end product is a tape with a time history of each field for every point. After the PBP data has been stored, it can be plotted. The plotting package is flexible, and allows the choice of fields and grid points. The one-dimensional fields are plotted as functions of time, and a daily summary is also produced. For the two-dimensional fields, time-height cross sections and daily summaries are produced. The model plotting package produces plots for all fields on the time average and pressure tapes. These plots include world maps, polar projections, zonal means, vectors, streamlines, and latitude-height plots.

Abeles, J.↗

The use of dynamical constraints in the analysis of mesoscale rawinsonde data

A process for data interpolation which seeks to maintain dynamical consistency is applied to the analysis of mesoscale rawinsonde data obtained from the National Severe Storms Laboratory. The analysis procedure is in two steps: a successive-corrections interpolation to gridpoints, followed by a variational adjustment which uses a set of forecast equations of horizontal momentum as dynamical constraints. The objective analysis is applied in case studies of conditions on May 22, 1966 and June 8, 1966, and it is shown that the adjusted winds and geopotentials have greatly reduced residuals compared to the initial fields. Adjustments to the wind and geopotential fields have root mean square values of 1-2 m/sec and 30-50 sq m/sec per sec, respectively. In addition, the vertical velocity fields computed from the adjusted winds are found to be in closer agreement with actual weather data.

Bloom, S. C.↗

Algebraic grid generation with corner singularities

A simple noniterative algebraic procedure is presented for generating smooth computational meshes on a quadrilateral topology. Coordinate distribution and normal derivative are provided on all boundaries, one of which may include a slope discontinuity. The boundary conditions are sufficient to guarantee continuity of global meshes formed of joined patches generated by the procedure. The method extends to 3-D. The procedure involves a synthesis of prior techniques stretching functions, cubic blending functions, and transfinite interpolation - to which is added the functional form of the corner solution. The procedure introduces the concept of generalized blending, which is implemented as an automatic scaling of the boundary derivatives for effective interpolation. Some implications of the treatment at boundaries for techniques solving elliptic PDE's are discussed in an Appendix.

Vinokur, M.↗

Gravity, topography, and crustal evolution of Venus

A method for estimating radial gravity anomalies from line-of-sight acceleration data, their interpolation, and the use of iteration for improved radial anomaly estimates is outlined. Results based on spherical harmonic coefficients for Venus' gravitational field are first discussed, and the procedure is then presented for determining radial acceleration estimates from the line-of-sight vectors and for interpolating the results at preselected locations using the method of objective mapping. A preliminary gravity anomaly map of much of the surface of Venus based on the first estimation step of this technique is presented, and the results are discussed.

Bowin, C.↗

Algebraic grid generation for wing-fuselage bodies

An algebraic procedure for the generation of boundary-fitted grids about wing-fuselage configurations is presented. A wing-fuselage configuration is specified by cross sections and mathematically represented by Coons' patches. A configuration is divided into sections so that several grid blocks that either adjoin each other or partially overlap each other can be generated, and each grid has six surfaces that map into a computational cube. Grids are first determined on the six boundary surfaces and then in the interior. Grid curves that are on the surface of the configuration are derived using plane-patch intersections, and single-valued functions relating approximate arc lengths along the curves to computational coordinates define the distribution of grid points. The two-boundary technique and transfinite interpolation are used to determine the boundary surface grids that are not on the configuration, and transfinite interpolation with linear blending functions is used to determine the interior grids.

Smith, R. E.↗

On the implementation of modal filters for control of structures

The most common technique for the control of structures is modal control. In modal control, the differential equations in terms of actual coordinates are replaced by a set of ordinary differential equations in terms of the modal coordinates known as modal equations. In designing feedback controls in conjunction with the modal equations, one must know the modal states for the modes targeted for control. The sensors measure actual states, however. The modal states can be estimated by means of a Luenberger observer or modal filters. The modal filters produce estimates of the modal states from distributed measurements of the states. If distributed measurements are not available, then they can be reconstructed from measurements at discrete points via interpolation. This paper examines various questions associated with the implementation of modal filters, such as the effect of choice of interpolation functions and sensors locations, as well as of measurement errors, on the state estimation process. The method is demonstrated by means of two numerical examples.

Meirovitch, L.↗

Multigrid techniques for the numerical solution of the diffusion equation

An accurate numerical solution of diffusion problems containing large local gradients can be obtained with a significant reduction in computational time by using a multigrid computational scheme. The spatial domain is covered with sets of uniform square grids of different sizes. The finer grid patterns overlap the coarse grid patterns. The finite-difference expressions for each grid pattern are solved independently by iterative techniques. Two interpolation methods were used to establish the values of the potential function on the fine grid boundaries with information obtained from the coarse grid solution. The accuracy and computational requirements for solving a test problem by a simple multigrid and a multilevel-multigrid method were compared. The multilevel-multigrid method combined with a Taylor series interpolation scheme was found to be best.

Phillips, R. E.↗

Forecast Impact of Components of the FGGE Observing System

Experiments have been conducted to assess the summer and winter forecast impact of the FGGE system, and of its main observing components: temperature sounding data derived from the TIROS-N polar orbiting satellite, cloud track winds determined from geostationary satellite observations and drifting buoy data which were collected by satellite during FGGE. The Analysis/Forecast System used has a number of improvements upon the system utilized by Halen et al. (1982) for the FGGE Special Observing Period-1 (SOP-1). Several modifications were made in the analysis scheme, the most important being the interpolation of the analysis minus 6 h forecast deviations rather than of the analyzed fields themselves. The forecast model is still the 4 deg lat, 5 deg lon and 9 vertical levels GLAS Fourth Order GCM with several minor corrections implemented in the physics and numerics. The improved vertical interpolation in the analysis resulted in better assimilation of rawinsonde data, which has more vertical structure than satellite data. As a result, there was an improvement of the forecasts derived from conventional data only, and, consequently, a small reduction of the positive impact of satellite data from that obtained by Halem et al. (1982).

Kalnay, E.↗

First results of GPS time transfer to Australia

A Global Positioning System time transfer unit built by NBS under contract to JPL was installed at Tidbinbilla Deep Space Communications Complex of the NASA Deep Space Network in June 1983. It has been used to estimate the relationship to UTC(USNO MC) of the Tidbinbilla frequency and time system TID(FTS) based on a hydrogen maser, and thence to estimate the performance of the Australian free-running time scale UTC(AUS). Data from the first three months has been analyzed three ways: by two-hop common view using JPL as intermediary; by long-arc interpolation of measurements against space vehicle clocks; and by long-arc interpolation of GPS-Time results. Residuals from a single quadratic fit through three months of UTC(USNO MC) - TID (FTS) results were white noise with standard error 15 ns, and a flying clock measurement gave 70 ns agreement. A straight line fit through results UTC(USNO MC) - UTC (AUS) gave 90 ns standard error and 120 bns agreement. It is proposed to use the GPS measurements to steer UTS(AUS) to UTC(BIH), and to rename the existing time scale TA(AUS).

Luck, J. M.↗

Digital image processing of Seabeam bathymetric data for structural studies of seamounts near the East Pacific Rise

The problem of displaying information on the seafloor morphology is attacked by utilizing digital image processing techniques to generate images for Seabeam data covering three young seamounts on the eastern flank of the East Pacific Rise. Errors in locations between crossing tracks are corrected by interactively identifying features and translating tracks relative to a control track. Spatial interpolation techniques using moving averages are used to interpolate between gridded depth values to produce images in shaded relief and color-coded forms. The digitally processed images clarify the structural control on seamount growth and clearly show the lateral extent of volcanic materials, including the distribution and fault control of subsidiary volcanic constructional features. The image presentations also clearly show artifacts related to both residual navigational errors and to depth or location differences that depend on ship heading relative to slope orientation in regions with steep slopes.

Edwards, M. H.↗

Thermodynamic evaluation of transonic compressor rotors using the finite volume approach

The finite volume explicit time marching method was refined and improved. Previously, extension had been made to the finite volume method to improve the accuracy of the calculation of total pressure in inviscid flow, extend the method to allow the calculation of laminar and turbulent boundary layers in internal flows, and improve the shock capturing properties of the method by introducing a Mach number dependent interpolation scheme for the pressure used in the calculating the density. The current work extends these developments by using the new pressure interpolation scheme in two dimensional viscous calculations, including a more complete description of the viscous stresses, introducing a criteria for the transverse upwind differencing which is a function of the ratio of transverse and streamwise mass fluxes, and allowing the calculation of internal flow where boundary layers are present on both walls of the duct. The manner in which the viscous stresses are evaluated in the nonorthogonal, nonuniform grid is detailed. The convergence is investigated and results for calculations of laminar flow in a converging duct are presented. Results for calculations of transonic flow in a converging-diverging nozzle are presented and the results are compared with Sajben's measurements and calculations by others.

Nicholson, S.↗

An interactive technique to generate digital elevation data using a vidicon camera

A technique is proposed for the semi-automatic digitization of topographic maps using a vidicon camera interfaced with a mini-computer system. The required input is a black and white map that shows only contour lines and corresponding elevation values such as the advance prints of 7.5 minute USGS topographic sheets. A set of image processing algorithms is implemented on an image of the contour map acquired by a vidicon camera. The system developed performs extraction of contour lines, tagging elevation values and subsequent interpolation of elevations to produce, in 5 sec. x 5 sec. intervals of latitude and longitude, digital elevation matrices. The interpolation of elevations for non-contour points is achieved using a steepest descent algorithm. A communications capability allows the final data to be transferred over telephone lines to an off-site unit such as the floppy disk or a micro-computer.

Sircar, J. K.↗

Review of image processing fundamentals

Image processing through convolution, transform coding, spatial frequency alterations, sampling, and interpolation are considered. It is postulated that convolution in one domain (real or frequency) is equivalent to multiplication in the other (frequency or real), and that the relative amplitudes of the Fourier components must be retained to reproduce any waveshape. It is suggested that all digital systems may be considered equivalent, with a frequency content approximately at the Nyquist limit, and with a Gaussian frequency response. An optimized cubic version of the interpolation continuum image is derived as a set of cubic spines. Pixel replication has been employed to enlarge the visable area of digital samples, however, suitable elimination of the extraneous high frequencies involved in the visable edges, by defocusing, is necessary to allow the underlying object represented by the data values to be seen.

Billingsley, F. C.↗

Monochromatic calculations of atmospheric radiative transfer due to molecular line absorption

Sensitivity studies related to the effects of line cutoff, spectral resolution, and temperature and pressure interpolations in radiative transfer have been performed so that a data set of absorption coefficients for water vapor, CO2, and O3 may be created efficiently. Results show that computations of absorption coefficients are affected only slightly by cutting a line off at a wave number 190 times the Lorentz half width from the center, or equivalently, cutting off 0.33 percent of the line intensity from the wings. To achieve a relative cooling rate error smaller than 2 percent, it is sufficient to precompute the absorption coefficient at three temperatures (210, 250, and 290 K) and 19 pressures with Delta (log 10 p) = 0.2. The absorption coefficient at other conditions can be interpolated linearly with pressure and exponentially with a quadratic in temperature. For the spectral resolution the absorption coefficients can be adequately computed at 0.01, 0.002, 0.005, and 0.025/cm intervals in the thermal water vapor, the CO2 and O3 bands, and the solar water vapor bands, respectively, which limits the error to only a few percent in the cooling and heating rates. Using the precomputed absorption coefficients, repeated monochromatic calculations of atmospheric heating/cooling rates for radiation model developments and for comparison with less detailed calculations are no longer difficult.

Chou, M.-D.↗