Search NASA⌕ Search

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.

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

At least 163 records · Page 9

Autonomous antenna tracking system for mobile symphonie ground stations

The implementation of a satellite tracking and antenna control system is described. Due to the loss of inclination control for the symphonie satellites, it became necessary to equip the parabolic antennas of the mobile Symphonie ground station with tracking facilities. For the relatively low required tracking accuracy of 0.5 dB, a low cost, step track system was selected. The step track system developed for this purpose and tested over a long period of time in 7 ground stations is based on a search step method with subsequent parabola interpolation. As compared with the real search step method, the system has the advantage of a higher pointing angle resolution, and thus a higher tracking accuracy. When the pilot signal has been switched off for a long period of time, as for instance after the eclipse, the antenna is repointed towards the satellite by an automatically initiated spiral search scan. The function and design of the tracking system are detailed, while easy handling and tracking results.

Ernsberger, K.↗

Assimilation of satellite surface wind speed data using the GLA analysis/forecast system

Research at the Goddard Laboratory for Atmospheres to evaluate methods to assign directions to both real and simulated surface wind speed data is discussed. Surface wind speed measurements are obtained from satellites, including Seasat, Nimbus-7, and Geosat. The methods include the interpolation of modeled forecast winds to the wind speed datum location, the use of surface pressures with a balance relation, and a variational analysis method. It is found that the best estimates may be obtained using a multipass approach to perform a surface analysis incorporating all conventional surface data.

Bloom, S. C.↗

Probabilistic Machine Learning Estimation of Ocean Mixed Layer Depth from Dense Satellite and Sparse In-Situ Observations

The ocean mixed layer plays an important role in the coupling between the upper ocean and atmosphere across a wide range of time scales. Estimation of the variability of the ocean mixed layer is therefore important for atmosphere-ocean prediction and analysis. The increasing coverage of in situ Argo profile data allows for an increasingly accurate analysis of the mixed layer depth (MLD) variability associated with deviations from the seasonal climatology. However, sampling rates are not sufficient to fully resolve subseasonal (<90 day) MLD variability. Yet, many multivariate observations-based analyses include implicit modeled subseasonal MLD variability. One analysis method is optimal interpolation of in situ data, but the interior analysis can be improved by leveraging surface data with regression or variational approaches. Here, we demonstrate how machine learning methods and satellite sea surface temperature, salinity, and height facilitate MLD estimation in a pilot study of two regions: the mid-latitude southern Indian and the eastern equatorial Pacific Oceans. We construct multiple machine learning architectures to produce weekly 1/2° gridded MLD anomaly fields (relative to a monthly climatology) with uncertainty estimates. We test multiple traditional and probabilistic machine learning techniques to compare both accuracy and probabilistic calibration. We validate our methodology by applying it to ocean model simulations. We find that incorporating sea surface data through a machine learning model improves the performance of spatiotemporal MLD variability estimation compared to optimal interpolation of Argo observations alone. These preliminary results are a promising first step for the application of machine learning to MLD prediction.

Machine Learning↗

A second-order accurate kinetic-theory-based method for inviscid compressible flows

An upwind method for the numerical solution of the Euler equations is presented. This method, called the kinetic numerical method (KNM), is based on the fact that the Euler equations are moments of the Boltzmann equation of the kinetic theory of gases when the distribution function is Maxwellian. The KNM consists of two phases, the convection phase and the collision phase. The method is unconditionally stable and explicit. It is highly vectorizable and can be easily made total variation diminishing for the distribution function by a suitable choice of the interpolation strategy. The method is applied to a one-dimensional shock-propagation problem and to a two-dimensional shock-reflection problem.

Deshpande, Suresh M.↗

Discrete Fourier transforms of nonuniformly spaced data

Time series or spatial series of measurements taken with nonuniform spacings have failed to yield fully to analysis using the Discrete Fourier Transform (DFT). This is due to the fact that the formal DFT is the convolution of the transform of the signal with the transform of the nonuniform spacings. Two original methods are presented for deconvolving such transforms for signals containing significant noise. The first method solves a set of linear equations relating the observed data to values defined at uniform grid points, and then obtains the desired transform as the DFT of the uniform interpolates. The second method solves a set of linear equations relating the real and imaginary components of the formal DFT directly to those of the desired transform. The results of numerical experiments with noisy data are presented in order to demonstrate the capabilities and limitations of the methods.

Swan, P. R.↗

Elastic Model Transitions Using Quadratic Inequality Constrained Least Squares

A technique is presented for initializing multiple discrete finite element model (FEM) mode sets for certain types of flight dynamics formulations that rely on superposition of orthogonal modes for modeling the elastic response. Such approaches are commonly used for modeling launch vehicle dynamics, and challenges arise due to the rapidly time-varying nature of the rigid-body and elastic characteristics. By way of an energy argument, a quadratic inequality constrained least squares (LSQI) algorithm is employed to e ect a smooth transition from one set of FEM eigenvectors to another with no requirement that the models be of similar dimension or that the eigenvectors be correlated in any particular way. The physically unrealistic and controversial method of eigenvector interpolation is completely avoided, and the discrete solution approximates that of the continuously varying system. The real-time computational burden is shown to be negligible due to convenient features of the solution method. Simulation results are presented, and applications to staging and other discontinuous mass changes are discussed

Orr, Jeb S.↗

Interpolating for the location of remote sensor data

An interpolation algorithm is presented as a practical alternative to common interpolation and approximation methods when applied to the problem of determining the location of remote sensor data. This algorithm is based upon knowledge of the geometry of the problem and is shown to be inherently more accurate than common interpolation schemes which may be applied to all types of data. A practical location problem is used to demonstrate its accuracy and computational cost.

Puccinelli, E. F.↗

Analysis of Doppler lidar wind measurements

New analysis methods were developed and refined for the synthesis of wind fields from multiple Doppler radar data. Considerable software was written and debugged for wind field synthesis. The analysis methods were applied to radar data collected during the CCOPE experiment in which lidar wind measurements were also made. The analysis method differs from existing methods in that interpolation (using filters with assigned weights) and advection of data, to account for time differential between data acquisition at different points, are not used. Instead, functional fits were used in space and time coordinates to transfer radial velocity data to a space-time grid. Subsequently the wind synthesis is performed. This method gives better space resolution in the synthesized winds, compared to other methods in common use, by better preserving the gradients of velocities. The wind fields show the cellular structure characteristic of boundary-layer convection. Also evident in the analysis are the sharp transitions and gradients in the wind field which were preserved by the synthesis.

Srivastava, R. C.↗

Comparison of Response Surface Construction Methods for Derivative Estimation Using Moving Least Squares, Kriging and Radial Basis Functions

Response construction methods using Moving Least Squares (MLS), Kriging and Radial Basis Functions (RBF) are compared with the Global Least Squares (GLS) method in three numerical examples for derivative generation capability. Also, a new Interpolating Moving Least Squares (IMLS) method adopted from the meshless method is presented. It is found that the response surface construction methods using the Kriging and RBF interpolation yields more accurate results compared with MLS and GLS methods. Several computational aspects of the response surface construction methods also discussed.

Krishnamurthy, Thiagarajan↗

The piecewise-parabolic method in curvilinear coordinates

We derive interpolation formulae for a third-order finite difference method in curvilinear, orthogonal coordinate systems. These formulae serve as a supplement to Colella and Woodward's PPM scheme for problems where the coordinate origin is included in the computational domain. Numerical examples of the improved accuracy of the advection scheme near coordinate singularities are shown.

Blondin, John M.↗

A case study of the sensitivity of forecast skill to data and data analysis techniques

A series of experiments have been conducted to examine the sensitivity of forecast skill to various data and data analysis techniques for the 0000 GMT case of January 21, 1979. These include the individual components of the FGGE observing system, the temperatures obtained with different satellite retrieval methods, and the method of vertical interpolation between the mandatory pressure analysis levels and the model sigma levels. It is found that NESS TIROS-N infrared retrievals seriously degrade a rawinsonde-only analysis over land, resulting in a poorer forecast over North America. Less degradation in the 72-hr forecast skill at sea level and some improvement at 500 mb is noted, relative to the control with TIROS-N retrievals produced with a physical inversion method which utilizes a 6-hr forecast first guess. NESS VTPR oceanic retrievals lead to an improved forecast over North America when added to the control.

Baker, W. E.↗

Global/local stress analysis of composite panels

A method for performing a global/local stress analysis is described, and its capabilities are demonstrated. The method employs spline interpolation functions which satisfy the linear plate bending equation to determine displacements and rotations from a global model which are used as boundary conditions for the local model. Then, the local model is analyzed independent of the global model of the structure. This approach can be used to determine local, detailed stress states for specific structural regions using independent, refined local models which exploit information from less-refined global models. The method presented is not restricted to having a priori knowledge of the location of the regions requiring local detailed stress analysis. This approach also reduces the computational effort necessary to obtain the detailed stress state. Criteria for applying the method are developed. The effectiveness of the method is demonstrated using a classical stress concentration problem and a graphite-epoxy blade-stiffened panel with a discontinuous stiffener.

Ransom, Jonathan B.↗

Global/local stress analysis of composite structures

A method for performing a global/local stress analysis is described and its capabilities are demonstrated. The method employs spline interpolation functions which satisfy the linear plate bending equation to determine displacements and rotations from a global model which are used as boundary conditions for the local model. Then, the local model is analyzed independent of the global model of the structure. This approach can be used to determine local, detailed stress states for specific structural regions using independent, refined local models which exploit information from less-refined global models. The method presented is not restricted to having a priori knowledge of the location of the regions requiring local detailed stress analysis. This approach also reduces the computational effort necessary to obtain the detailed stress state. Criteria for applying the method are developed. The effectiveness of the method is demonstrated using a classical stress concentration problem and a graphite-epoxy blade-stiffened panel with a discontinuous stiffener.

Ransom, Jonathan B.↗

Global/local stress analysis of composite panels

A method for performing a global/local stress analysis is described and its capabilities are demonstrated. The method employs spline interpolation functions which satisfy the linear plate bending equation to determine displacements and rotations from a global model which are used as 'boundary conditions' for the local model. Then, the local model is analyzed independently of the global model of the structure. This approach can be used to determine local, detailed stress states for specific structural regions using independent, refined local models which exploit information from less-refined global models. The method presented is not restricted to having a priori knowledge of the location of the regions requiring local detailed stress analysis. This approach also reduces the computational effort necessary to obtain the detailed stress state. Criteria for applying the method are developed. The effectiveness of the method is demonstrated using a classical stress concentration problem and a graphite-epoxy blade-stiffened panel with a discontinuous stiffener.

Ransom, J. B.↗

Estimation of total ozone from satellite measurements of backscattered ultraviolet earth radiance.

Total ozone is estimated from Nimbus IV satellite measurements of the attenuation of backscattered radiances at wavelengths between 3100 and 3400 A. A measurement of the backscattered radiance at 3800 A, outside the ozone absorption band, is used to determine an equivalent Lambert albedo for the cloud-ground-haze surface viewed by the instrument. The measured relative attenuation at two wavelengths is compared with such values precomputed for a series of standard ozone profiles and corrected for the equivalent Lambert albedo. Total ozone is obtained by interpolation. Two alternative methods are used to assign an equivalent Lambert albedo at the absorbing wavelengths. Total ozone values estimated by these methods are compared with the Dobson (ground-truth) data by linear regression. The available evidence suggests that the true standard error of the satellite data may be 0.015 atm-cm or less for solar zenith angles smaller than 60 deg.

Mateer, C. L.↗

Higher-order flux difference splitting schemes for the Euler equations using upstream interpolations

A class of explicit two time-level, 2p + 1 space-point, (2p 1)-th order, upwind-biased flux difference splitting schemes are proposed for the numerical advection based on Lagrange's interpolation, and the method is an accord with the physical domain of dependence. A normalized Jacobian coefficient matrix is introduced to convert the schemes to hyperbolic systems of conservation laws, and approaches to make the higher-order schemes total variation stable are discussed. Accuracy and stability of the present schemes are examined, and implicit total variation diminishing schemes are developed for steady-state calculations.Application to gasdynamic problems for both steady and unsteady flows covering a wide range of Mach numbers is considered, and results for a blast wave passing a cylinder, and head-on collision of two blast waves over a circular arc, are presented. The flow patterns were found to be symmetric, and good resolution of flow structures was obtained.

Yang, J. Y.↗

Retrieval of aerosol size distribution moments from multiwavelength particulate extinction data

An iterative method for deriving the particle size distribution from multiwavelength satellite remote sensing data is presented and its use illustrated with an application to SAGE II data. The compression size distribution method (CSD) features interpolation and look-up features for relating the coefficients of a least-squares fit as a function of wavelength to parameters of the compression function that have been obtained by Mie calculations. Results are reported of comparisons between particulate radii predicted for the El Chichon volcano aerosol based on SAGE II data and 14 radii measured with NASA flights through the aerosol with aircraft carrying wire impactors.

Livingston, J. M.↗

Automatic computation of Euler-marching and subsonic grids for wing-fuselage configurations

Algebraic procedures are described for the automatic generation of structured, single-block flow computation grids for relatively simple configurations (wing, fuselage, and fin). For supersonic flows, a quasi two-dimensional grid for Euler-marching codes is developed, and some sample results in graphical form are included. A type of grid for subsonic flow calculation is also described. The techniques are algebraic and are based on a generalization of the method of transfinite interpolation.

Barger, Raymond L.↗