SEARCH · Search NASA

Results for “interpolating methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Algebraic grid generation for complex geometries

An efficient computer program called GRID2D/3D has been developed to generate single and composite grid systems within geometrically complex two- and three-dimensional (2D and 3D) spatial domains that can deform with time. GRID2D/3D generates single grid systems by using algebraic grid generation methods based on transfinite interpolation. The distribution of grid points within the spatial domain is controlled by stretching functions and grid lines can intersect boundaries of the spatial domain orthogonally. GRID2D/3D generates composite grid systems by patching together two or more single grid systems. The patching can be discontinuous or continuous. For 2D spatial domains the boundary curves are constructed by using either cubic or tension spline interpolation. For 3D spatial domains the boundary surfaces are constructed by using a new technique, developed in this study, referred to as 3D bidirectional Hermite interpolation.

Shih, T. I.-P.↗

NASA NTRS1993

Use of Fourier transforms for asynoptic mapping: Applications to the Upper Atmosphere Research Satellite microwave limb sounder

Fourier analysis has been applied to data obtained from limb viewing instruments on the Upper Atmosphere Research Satellite. A coordinate system rotation facilitates the efficient computation of Fourier transforms in the temporal and longitudinal domains. Fields such as ozone (O3), chlorine monoxide (ClO), temperature, and water vapor have been transformed by this process. The transforms have been inverted to provide maps of these quantities at selected times, providing a method of accurate time interpolation. Maps obtained by this process show evidence of both horizontal and vertical transport of important trace species such as O3 and ClO. An examination of the polar regions indicates that large-scale planetary variations are likely to play a significant role in transporting midstratospheric O3 into the polar regions. There is also evidence that downward transport occurs, providing a means of moving O3 into the polar vortex at lower altitudes. The transforms themselves show the structure and propagation characteristics of wave variations.

Elson, Lee S.↗

NASA NTRS1996

Propagation and Breaking at High Altitudes of Gravity Waves Excited by Tropospheric Forcing

An anelastic approximation is used with a time-variable coordinate transformation to formulate a two-dimensional numerical model that describes the evolution of gravity waves. The model is solved using a semi-Lagrangian method with monotone (nonoscillatory) interpolation of all advected fields. The time-variable transformation is used to generate disturbances at the lower boundary that approximate the effect of a traveling line of thunderstorms (a squall line) or of flow over a broad topographic obstacle. The vertical propagation and breaking of the gravity wave field (under conditions typical of summer solstice) is illustrated for each of these cases. It is shown that the wave field at high altitudes is dominated by a single horizontal wavelength; which is not always related simply to the horizontal dimension of the source. The morphology of wave breaking depends on the horizontal wavelength; for sufficiently short waves, breaking involves roughly one half of the wavelength. In common with other studies, it is found that the breaking waves undergo "self-acceleration," such that the zonal-mean intrinsic frequency remains approximately constant in spite of large changes in the background wind. It is also shown that many of the features obtained in the calculations can be understood in terms of linear wave theory. In particular, linear theory provides insights into the wavelength of the waves that break at high altitudes, the onset and evolution of breaking. the horizontal extent of the breaking region and its position relative to the forcing, and the minimum and maximum altitudes where breaking occurs. Wave breaking ceases at the altitude where the background dissipation rate (which in our model is a proxy for molecular diffusion) becomes greater than the rate of dissipation due to wave breaking, This altitude, in effect, the model turbopause, is shown to depend on a relatively small number of parameters that characterize the waves and the background state.

Prusa, Joseph M.↗

NASA NTRS2010

ISS Solar Array Management

The International Space Station (ISS) Solar Array Management (SAM) software toolset provides the capabilities necessary to operate a spacecraft with complex solar array constraints. It monitors spacecraft telemetry and provides interpretations of solar array constraint data in an intuitive manner. The toolset provides extensive situational awareness to ensure mission success by analyzing power generation needs, array motion constraints, and structural loading situations. The software suite consists of several components including samCS (constraint set selector), samShadyTimers (array shadowing timers), samWin (visualization GUI), samLock (array motion constraint computation), and samJet (attitude control system configuration selector). It provides high availability and uptime for extended and continuous mission support. It is able to support two-degrees-of-freedom (DOF) array positioning and supports up to ten simultaneous constraints with intuitive 1D and 2D decision support visualizations of constraint data. Display synchronization is enabled across a networked control center and multiple methods for constraint data interpolation are supported. Use of this software toolset increases flight safety, reduces mission support effort, optimizes solar array operation for achieving mission goals, and has run for weeks at a time without issues. The SAM toolset is currently used in ISS real-time mission operations.

Williams, James P.↗

NASA NTRS1967

Poisson equation on a sphere. II.

Interpolation of Poisson equation by method adjusting equations to insure existence of discrete solution

DISCRETE FUNCTION↗

NASA NTRS1990

Global convergence of inexact Newton methods for transonic flow

In computational fluid dynamics, nonlinear differential equations are essential to represent important effects such as shock waves in transonic flow. Discretized versions of these nonlinear equations are solved using iterative methods. In this paper an inexact Newton method using the GMRES algorithm of Saad and Schultz is examined in the context of the full potential equation of aerodynamics. In this setting, reliable and efficient convergence of Newton methods is difficult to achieve. A poor initial solution guess often leads to divergence or very slow convergence. This paper examines several possible solutions to these problems, including a standard local damping strategy for Newton's method and two continuation methods, one of which utilizes interpolation from a coarse grid solution to obtain the initial guess on a finer grid. It is shown that the continuation methods can be used to augment the local damping strategy to achieve convergence for difficult transonic flow problems. These include simple wings with shock waves as well as problems involving engine power effects. These latter cases are modeled using the assumption that each exhaust plume is isentropic but has a different total pressure and/or temperature than the freestream.

Young, David P.↗

NASA NTRS1981

Optimal interpolation and the Kalman filter

The estimation theory of stochastic-dynamic systems is described and used in a numerical study of optimal interpolation. The general form of data assimilation methods is reviewed. The Kalman-Bucy, KB filter, and optimal interpolation (OI) filters are examined for effectiveness in performance as gain matrices using a one-dimensional form of the shallow-water equations. Control runs in the numerical analyses were performed for a ten-day forecast in concert with the OI method. The effects of optimality, initialization, and assimilation were studied. It was found that correct initialization is necessary in order to localize errors, especially near boundary points. Also, the use of small forecast error growth rates over data-sparse areas was determined to offset inaccurate modeling of correlation functions near boundaries.

Cohn, S.↗

NASA NTRS1983

Spectral multigrid methods for elliptic equations 2

A detailed description of spectral multigrid methods is provided. This includes the interpolation and coarse-grid operators for both periodic and Dirichlet problems. The spectral methods for periodic problems use Fourier series and those for Dirichlet problems are based upon Chebyshev polynomials. An improved preconditioning for Dirichlet problems is given. Numerical examples and practical advice are included.

Zang, T. A.↗

NASA NTRS1984

Spectral multigrid methods for elliptic equations II

Zang, T. A.↗

NASA NTRS2007

Smooth Phase Interpolated Keying

Smooth phase interpolated keying (SPIK) is an improved method of computing smooth phase-modulation waveforms for radio communication systems that convey digital information. SPIK is applicable to a variety of phase-shift-keying (PSK) modulation schemes, including quaternary PSK (QPSK), octonary PSK (8PSK), and 16PSK. In comparison with a related prior method, SPIK offers advantages of better performance and less complexity of implementation. In a PSK scheme, the underlying information waveform that one seeks to convey consists of discrete rectangular steps, but the spectral width of such a waveform is excessive for practical radio communication. Therefore, the problem is to smooth the step phase waveform in such a manner as to maintain power and bandwidth efficiency without incurring an unacceptably large error rate and without introducing undesired variations in the amplitude of the affected radio signal. Although the ideal constellation of PSK phasor points does not cause amplitude variations, filtering of the modulation waveform (in which, typically, a rectangular pulse is converted to a square-root raised cosine pulse) causes amplitude fluctuations. If a power-efficient nonlinear amplifier is used in the radio communication system, the fluctuating-amplitude signal can undergo significant spectral regrowth, thus compromising the bandwidth efficiency of the system. In the related prior method, one seeks to solve the problem in a procedure that comprises two major steps: phase-value generation and phase interpolation. SPIK follows the two-step approach of the related prior method, but the details of the steps are different. In the phase-value-generation step, the phase values of symbols in the PSK constellation are determined by a phase function that is said to be maximally smooth and that is chosen to minimize the spectral spread of the modulated signal. In this step, the constellation is divided into two groups by assigning, to information symbols, phase values that result in equal numbers of clockwise and counter-clockwise phase rotations for equally likely symbols. The purpose served by assigning phase values in this way is to prevent unnecessary generation of spectral lines and prevent net shifts of the carrier signal. In the phase-interpolation step, the smooth phase values are interpolated over a number, n, of consecutive symbols (including the present symbol) by means of an unconventional spline curve fit.

Borah, Deva K.↗

NASA NTRS2016

Simple Test Functions in Meshless Local Petrov-Galerkin Methods

Two meshless local Petrov-Galerkin (MLPG) methods based on two different trial functions but that use a simple linear test function were developed for beam and column problems. These methods used generalized moving least squares (GMLS) and radial basis (RB) interpolation functions as trial functions. These two methods were tested on various patch test problems. Both methods passed the patch tests successfully. Then the methods were applied to various beam vibration problems and problems involving Euler and Beck's columns. Both methods yielded accurate solutions for all problems studied. The simple linear test function offers considerable savings in computing efforts as the domain integrals involved in the weak form are avoided. The two methods based on this simple linear test function method produced accurate results for frequencies and buckling loads. Of the two methods studied, the method with radial basis trial functions is very attractive as the method is simple, accurate, and robust.

Raju, Ivatury S.↗

NASA NTRS2014

On Formulations of Discontinuous Galerkin and Related Methods for Conservation Laws

A formulation for the discontinuous Galerkin (DG) method that leads to solutions using the differential form of the equation (as opposed to the standard integral form) is presented. The formulation includes (a) a derivative calculation that involves only data within each cell with no data interaction among cells, and (b) for each cell, corrections to this derivative that deal with the jumps in fluxes at the cell boundaries and allow data across cells to interact. The derivative with no interaction is obtained by a projection, but for nodal-type methods, evaluating this derivative by interpolation at the nodal points is more economical. The corrections are derived using the approximate (Dirac) delta functions. The formulation results in a family of schemes: different approximate delta functions give rise to different methods. It is shown that the current formulation is essentially equivalent to the flux reconstruction (FR) formulation. Due to the use of approximate delta functions, an energy stability proof simpler than that of Vincent, Castonguay, and Jameson (2011) for a family of schemes is derived. Accuracy and stability of resulting schemes are discussed via Fourier analyses. Similar to FR, the current formulation provides a unifying framework for high-order methods by recovering the DG, spectral difference (SD), and spectral volume (SV) schemes. It also yields stable, accurate, and economical methods.

Huynh, H. T.↗

NASA NTRS2012

Method for Pre-Conditioning a Measured Surface Height Map for Model Validation

This software allows one to up-sample or down-sample a measured surface map for model validation, not only without introducing any re-sampling errors, but also eliminating the existing measurement noise and measurement errors. Because the re-sampling of a surface map is accomplished based on the analytical expressions of Zernike-polynomials and a power spectral density model, such re-sampling does not introduce any aliasing and interpolation errors as is done by the conventional interpolation and FFT-based (fast-Fourier-transform-based) spatial-filtering method. Also, this new method automatically eliminates the measurement noise and other measurement errors such as artificial discontinuity. The developmental cycle of an optical system, such as a space telescope, includes, but is not limited to, the following two steps: (1) deriving requirements or specs on the optical quality of individual optics before they are fabricated through optical modeling and simulations, and (2) validating the optical model using the measured surface height maps after all optics are fabricated. There are a number of computational issues related to model validation, one of which is the "pre-conditioning" or pre-processing of the measured surface maps before using them in a model validation software tool. This software addresses the following issues: (1) up- or down-sampling a measured surface map to match it with the gridded data format of a model validation tool, and (2) eliminating the surface measurement noise or measurement errors such that the resulted surface height map is continuous or smoothly-varying. So far, the preferred method used for re-sampling a surface map is two-dimensional interpolation. The main problem of this method is that the same pixel can take different values when the method of interpolation is changed among the different methods such as the "nearest," "linear," "cubic," and "spline" fitting in Matlab. The conventional, FFT-based spatial filtering method used to eliminate the surface measurement noise or measurement errors can also suffer from aliasing effects. During re-sampling of a surface map, this software preserves the low spatial-frequency characteristic of a given surface map through the use of Zernike-polynomial fit coefficients, and maintains mid- and high-spatial-frequency characteristics of the given surface map by the use of a PSD model derived from the two-dimensional PSD data of the mid- and high-spatial-frequency components of the original surface map. Because this new method creates the new surface map in the desired sampling format from analytical expressions only, it does not encounter any aliasing effects and does not cause any discontinuity in the resultant surface map.

Sidick, Erkin↗

NASA NTRS1991

An evaluation of HEMT potential for millimeter-wave signal sources using interpolation and harmonic balance techniques

A large-signal analysis method based on an harmonic balance technique and a 2-D cubic spline interpolation function has been developed and applied to the prediction of InP-based HEMT oscillator performance for frequencies extending up to the submillimeter-wave range. The large-signal analysis method uses a limited number of DC and small-signal S-parameter data and allows the accurate characterization of HEMT large-signal behavior. The method has been validated experimentally using load-pull measurement. Oscillation frequency, power performance, and load requirements are discussed, with an operation capability of 300 GHz predicted using state-of-the-art devices (fmax is approximately equal to 450 GHz).

Kwon, Youngwoo↗

NASA NTRS1978

Interpolation and Extrapolation of Creep Rupture Data by the Minimum Commitment Method. Part 3: Analysis of Multiheats

The Minimum Commitment Method was applied to two sets of data for which multiple heat information was available. For one alloy, a 304 stainless steel studied in Japan, data on nine well characterized heats were used, while for a proprietary low alloy carbon steel studied in the United Kingdom data were available on seven heats - in many cases to very long rupture times. For this preliminary study no instability factors were used. It was discovered that heat-to-heat variations would be accounted for by introducing heat identifiers in the form A + B log sigma where sigma is the stress and the constants A and B depend only on the heat. With these identifiers all the data could be collapsed onto a single master curve, even though there was considerable scatter among heats. Using these identifiers together with the average behavior of all heats made possible the determination of an accurate constitutive equation for each individual heat. Two basic approaches are discussed for applying the results of the analysis.

Manson, S. S.↗

NASA NTRS1978

Interpolation and extrapolation of creep rupture data by the minimum commitment method. II - Oblique translation

An outline is presented of a new procedure, termed the oblique translation method, which emerged in the development of the focal point convergence method. Approaches for implementing the oblique translation method are discussed. It is shown that the new method is essentially a minimum commitment method when manually-graphically implemented, in the sense that the form of the functions involved are not forced into particular analytical forms. The individual constants and functions are independently determined. The minimum commitment concept is extended to the analysis of creep rupture data wherein each isothermal is to be generated by an oblique translation of the 'master curve' when plotted. Attention is given to a manual-graphical analysis, the preassessment of data, and an analysis by computer code.

Manson, S. S.↗

NASA NTRS1978

Interpolation and extrapolation of creep rupture data by the minimum commitment method. I - Focal-point convergence

The minimum commitment method has been applied to the analysis of creep-rupture data. The method is based on a parameter representing the focal point of convergence of all isothermals when extended to the long or short times necessary for such convergence to occur. The technique may be applied by manual-graphic analysis on computer code. It is illustrated for the nickel-base alloy Astroloy.

Manson, S. S.↗

NASA NTRS1978

Interpolation and extrapolation of creep rupture data by the minimum commitment method. III - Analysis of multiheats

An outline is presented of approaches for treating multiheats on the basis of the focal point convergence method. The method has been employed in the case of two highly characterized multiheats, including a 304 stainless steel studied in Japan and a low alloy carbon steel studied in England. The method makes use of the same functional form for all materials. Only the constants are varied for each multiheat. Completely computerized procedures are employed for the determination of the constants. Once the basic analysis has been performed, the representation of various members in the same system is achieved by adding linear expressions of log stress, changing only two constants in the equations to represent a selected heat.

Manson, S. S.↗