Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “horizontal 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

G-Function Library for Modeling Vertical Bore Ground Heat Exchanger

This library contains g-functions (thermal response functions) for standard, regularly spaced vertical borehole ground heat exchangers. In total, it contains 34, 321 configurations. To permit interpolation, each configuration has g-functions for heights of 24, 48, 96, 192, and 384 m. All the g-functions were calculated with burial depths of 2m, and borehole diameters of 15 to 17.5 cm, depending on height. In configurations with uniform spacing, the spacing between the boreholes is set to 5m, though it can be scaled to other horizontal spacings.

15 GEOTHERMAL ENERGY↗
NASA NTRS1981

Low-level circulations induced by an upper-level jet in a fine mesh numerical model

Attention is given to varying the initial conditions in the Penn State three-dimensional numerical model in the development of the low-level jet over the Ohio Valley during the 10 May 1973 severe weather event. The results of two forecasts are summarized: (1) the control forecast which interpolated the National Meteorological Center's operational analysis to the mesoscale model grid; and (2) modifying the control forecast by inserting artificial soundings and surface observations. Enhancements are made on the basis of subsynoptic scale analyses by Uccellini and surface data. The results indicate that the low-level jet was in part induced by pressure changes associated with deceleration and horizontal divergence of air in the exit region of the upper-level jet. The enhancement of the initial moisture field using surface reports of clouds and precipitation has a substantial effect on precipitation.

Burkhart, R. P.↗
NASA NTRS2018

File Specification for M2AMIP Products: Version Number - 1.0

The Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) is an atmospheric reanalysis computed with the Goddard Earth Observing (EOS) System, Version 5.12.4 (GEOS) data assimilation system (Gelaro et al., 2017). To supplement the reanalysis, the GEOS General Circulation Model (GCM) used in MERRA-2 has been used to generate a 10-member ensemble of simulations, configured following the convention of the Atmospheric Model Intercomparison Project (AMIP; Gates et al., 1992). Each ensemble member was initialized using meteorological fields from a different date in November 1979. The AMIP simulations used the sea-surface temperature (SST) and sea-ice boundary conditions that were used in MERRA-2 (Bosilovich et al., 2016). This 10-member ensemble of AMIP simulations, denoted M2AMIP, is available for download in a group of self-describing files, which are documented in this office note. All data collections are provided on the same horizontal grid as MERRA-2. This grid has 576 points in the longitudinal direction and 361 points in the latitudinal direction, corresponding to a resolution of 0.625 degrees by 0.5 degrees. Although data collections are available at this grid, all fields are computed on a cubed-sphere grid with an approximate resolution of 50 km by 50 km and are then spatially interpolated to the latitude-longitude grid. There are no changes in the vertical grids used: variables are provided on either the native vertical grid of 72 model layers, or interpolated to 42 standard pressure levels. Unlike MERRA, no data collections are available at the vertical layer edges. More details on the grid are provided in Section 4. MERRA-2 introduced observation-based precipitation forcing for the land surface parameterization and the corresponding variable PRECTOTCORR in the MERRA-2 FLX (surface turbulent fluxes and related quantities) and LFO (land-surface forcing) collections (see Section 6; Reichle et al., 2017). While this variable is still available for M2AMIP, there was no observation-based forcing, making the value identical to the model derived precipitation, PRECTOT. Similarly, without data assimilation, the values for the analysis increments, D*DTANA, in the tendency and vertically integrated file collections are zero. The M2AMIP data are available for download online through the NASA Center for Climate Simulation (NCCS) DataPortal (https://portal.nccs.nasa.gov/datashare/gmao_m2amip/). Data are arranged in subdirectories based on ensemble member, followed by year and month. Control files that are compatible with the Grid Analysis and Display System (GrADS) are available in the ctl_daily and ctl_monthly directories for the hourly, three hourly, and monthly mean data. Control files for the monthly mean diurnal cycle can be found in the ctl_diurnal subdirectory within the directory for each individual ensemble member.

M2AMIP↗
NASA NTRS2020

Improved LOLA Elevation Maps for South Pole Landing Sites: Error Estimates and Their Impact on Illumination Conditions

We present new high-resolution topographic models of 4 high-priority lunar south pole landing sites based exclusively on the laser altimetry data acquired by the Lunar Orbiter Laser Altimeter (LOLA) onboard the Lunar Reconnaissance Orbiter. By iteratively adjusting the LOLA tracks to the LOLA-based digital elevation model (LDEM) in a self-consistent fashion, we reduce the orbital geolocation errors by over a factor of 10 such that the new ground track geolocation uncertainty is ~10–20 ​cm horizontally and ~2–4 ​cm vertically over each 16 ​× ​16 km region. These new and improved 5 ​m/pix LDEMs will be useful to constrain higher-resolution topographic models derived from imagery, which are not as well controlled geodetically and which can be hindered by shadows. We developed a method to estimate surface height uncertainty in the new LDEMs, which accounts for the reduced orbital errors and interpolation errors by assuming a fractal behavior for the short-scale topography. The LDEM surface height and slope uncertainties have typical RMS values of ~0.30–0.50 ​m and ~1.5–2.5°, respectively. Finally, we examine how height uncertainties propagate to variations in horizon elevation and thus the predicted illumination conditions at these polar latitudes, and we show how this error characterization can inform landing site studies.

Michael K Barker↗
NASA NTRS2005

Applications and Improvement of a Coupled, Global and Cloud-Resolving Modeling System

Recently Grabowski (2001) and Khairoutdinov and Randall (2001) have proposed the use of 2D CFWs as a "super parameterization" [or multi-scale modeling framework (MMF)] to represent cloud processes within atmospheric general circulation models (GCMs). In the MMF, a fine-resolution 2D CRM takes the place of the single-column parameterization used in conventional GCMs. A prototype Goddard MMF based on the 2D Goddard Cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM) is now being developed. The prototype includes the fvGCM run at 2.50 x 20 horizontal resolution with 32 vertical layers from the surface to 1 mb and the 2D (x-z) GCE using 64 horizontal and 32 vertical grid points with 4 km horizontal resolution and a cyclic lateral boundary. The time step for the 2D GCE would be 15 seconds, and the fvGCM-GCE coupling frequency would be 30 minutes (i.e. the fvGCM physical time step). We have successfully developed an fvGCM-GCE coupler for this prototype. Because the vertical coordinate of the fvGCM (a terrain-following floating Lagrangian coordinate) is different from that of the GCE (a z coordinate), vertical interpolations between the two coordinates are needed in the coupler. In interpolating fields from the GCE to fvGCM, we use an existing fvGCM finite- volume piecewise parabolic mapping (PPM) algorithm, which conserves the mass, momentum, and total energy. A new finite-volume PPM algorithm, which conserves the mass, momentum and moist static energy in the z coordinate, is being developed for interpolating fields from the fvGCM to the GCE. In the meeting, we will discuss the major differences between the two MMFs (i.e., the CSU MMF and the Goddard MMF). We will also present performance and critical issues related to the MMFs. In addition, we will present multi-dimensional cloud datasets (i.e., a cloud data library) generated by the Goddard MMF that will be provided to the global modeling community to help improve the representation and performance of moist processes in climate models and to improve our understanding of cloud processes globally (the software tools needed to produce cloud statistics and to identify various types of clouds and cloud systems from both high-resolution satellite and model data will be also presented).

Tao, W.-K.↗
NASA NTRS2009

The Impact of the Assimilation of AIRS Radiance Measurements on Short-term Weather Forecasts

Advanced spaceborne instruments have the ability to improve the horizontal and vertical characterization of temperature and water vapor in the atmosphere through the explicit use of hyperspectral thermal infrared radiance measurements. The incorporation of these measurements into a data assimilation system provides a means to continuously characterize a three-dimensional, instantaneous atmospheric state necessary for the time integration of numerical weather forecasts. Measurements from the National Aeronautics and Space Administration (NASA) Atmospheric Infrared Sounder (AIRS) are incorporated into the gridpoint statistical interpolation (GSI) three-dimensional variational (3D-Var) assimilation system to provide improved initial conditions for use in a mesoscale modeling framework mimicking that of the operational North American Mesoscale (NAM) model. The methodologies for the incorporation of the measurements into the system are presented. Though the measurements have been shown to have a positive impact in global modeling systems, the measurements are further constrained in this system as the model top is physically lower than the global systems and there is no ozone characterization in the background state. For a study period, the measurements are shown to have positive impact on both the analysis state as well as subsequently spawned short-term (0-48 hr) forecasts, particularly in forecasted geopotential height and precipitation fields. At 48 hr, height anomaly correlations showed an improvement in forecast skill of 2.3 hours relative to a system without the AIRS measurements. Similarly, the equitable threat and bias scores of precipitation forecasts of 25 mm (6 hr)-1 were shown to be improved by 8% and 7%, respectively.

McCarty, Will↗
NASA NTRS1991

Concept for the subresolution measurement of earthquake strain fields using SPOT panchromatic imagery

The use of satellite-derived imagery for measuring subresolution horizontal terrain displacments associated with present-day earthquakes is discussed with reference to data from the French SPOT satellite whose sensor array provides 10-m panchromatic imagery. The measured terrain displacements can be up to several meters, but usually no more than 6-8 m even for major earthquakes. The general approach is to spatially match the after image to the before image at each point on a half-kilometer grid by iteratively interpolating one and testing its correlation with the other. The discussion covers the basic algorithm, results of initial tests, error types and limitations, and future work.

Crippen, Robert E.↗
NASA NTRS1975

Large-scale kinetic energy spectra from Eulerian analysis of EOLE wind data

A data set of 56,000 winds determined from the horizontal displacements of EOLE balloons at the 200 mb level in the Southern Hemisphere during the period October 1971-February 1972 is utilized for the computation of planetary- and synoptic-scale kinetic energy space spectra. However, the random distribution of measurements in space and time presents some problems for the spectral analysis. Two different approaches are used, i.e., a harmonic analysis of daily wind values at equi-distant points obtained by space-time interpolation of the data, and a correlation method using the direct measurements. Both methods give similar results for small wavenumbers, but the second is more accurate for higher wavenumbers (k above or equal to 10). The spectra show a maximum at wavenumbers 5 and 6 due to baroclinic instability and then decrease for high wavenumbers up to wavenumber 35 (which is the limit of the analysis), according to the inverse power law k to the negative p, with p close to 3.

Desbois, M.↗
NASA NTRS1971

Bayesian recursive image estimation.

A procedure for recursively estimating images that are characterized statistically by the mean and correlation functions associated with the random process representing the brightness level is proposed for the case where the images are corrupted by additive noise. First, a dynamic model is developed with a response characteristic which matches that of the scanner output (the input of the estimator is the output of a horizontal line scanner) in a statistical sense. Such models have the form of an ordinary differential or difference equation with white noise input. An insignificant approximation is introduced by using a constant-coefficient model. The appropriate model is a vector valued difference equation with the solution representing a vector Markov process. The next step is to obtain the minimum mean square estimate of the image by using a Kalman filter. Since the image estimation is an interpolation problem, two successive runs over the observation are performed in opposite directions and the resultant estimates are averaged. Examples are included for illustration.

Nahi, N. E.↗
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 NTRS1985

Noise correction on LANDSAT images using a spline-like algorithm

Many applications using LANDSAT images face a dilemma: the user needs a certain scene (for example, a flooded region), but that particular image may present interference or noise in form of horizontal stripes. During automatic analysis, this interference or noise may cause false readings of the region of interest. In order to minimize this interference or noise, many solutions are used, for instane, that of using the average (simple or weighted) values of the neighboring vertical points. In the case of high interference (more than one adjacent line lost) the method of averages may not suit the desired purpose. The solution proposed is to use a spline-like algorithm (weighted splines). This type of interpolation is simple to be computer implemented, fast, uses only four points in each interval, and eliminates the necessity of solving a linear equation system. In the normal mode of operation, the first and second derivatives of the solution function are continuous and determined by data points, as in cubic splines. It is possible, however, to impose the values of the first derivatives, in order to account for shapr boundaries, without increasing the computational effort. Some examples using the proposed method are also shown.

Vijaykumar, N. L.↗
NASA NTRS2007

Large and Small Droplet Impingement Data on Airfoils and Two Simulated Ice Shapes

Water droplet impingement data were obtained at the NASA Glenn Icing Research Tunnel (IRT) for four wings and one wing with two simulated ice shapes. The wings tested include three 36-in. chord wings (MS(1)-317, GLC-305, and a NACA 652-415) and a 57-in. chord Twin Otter horizontal tail section. The simulated ice shapes were 22.5- and 45-min glaze ice shapes for the Twin Otter horizontal tail section generated using the LEWICE 2.2 ice accretion program. The impingement experiments were performed with spray clouds having median volumetric diameters of 11, 21, 79, 137, and 168 mm. Comparisons to the experimental data were generated which showed good agreement for the clean wings and ice shapes at lower drop sizes. For larger drop sizes LEWICE 2.2 over predicted the collection efficiencies due to droplet splashing effects which were not modeled in the program. Also for the more complex glaze ice shapes interpolation errors resulted in the over prediction of collection efficiencies in cove and shadow regions of ice shapes.

Papadakis, Michael↗
DOE OSTI2023

A phased workflow to define permit‐ready locations for large volume CO 2 injection and storage

Abstract To‐date, only two UIC Class VI permits have been issued by the US Environmental Protection Agency. We illustrate a four‐phase workflow to first identify regional storage resources and then down‐select sites to yield permit‐ready locations that can accept and store large volumes of CO 2 . Specific permit requirements should guide objectives and define deliverables of respective workflow phases. In the first phase we used available regional data and screened structure and injection zones to locate resources that match CO 2 volumes planned to be captured. Available data were also used to assess presence and depth of usable groundwater, the key resource being protected via permitting. We then used advanced, closed‐form, analytical solutions (EASiTool) to estimate CO 2 injectivity into each hydrologically connected injection compartment. In the second phase we acquired and conditioned additional wireline logs and leased available seismic datasets. We interpreted the depositional systems from wireline well‐log character and mapped sandbody geometry to interpolate injection and confining‐zone distribution. Using available data, we mapped faults and locations of freshwater and overpressure (or other capacity‐limiting geologic parameters) in more detail. In the third phase, we used the augmented geologic data to develop a static model for the selected area, extracted the areas of highest interest, and generated and ran dynamic (flow) models. In a fourth phase, we reduced major uncertainties identified in earlier phases. Our case study indicates that to complete preparation of a permit application requires (1) improved lithologic characterization information (thicknesses and horizontal and vertical connectivity) and (2) better definition of poorly defined local faults. © 2023 The Authors. Greenhouse Gases: Science and Technology published by Society of Chemical Industry and John Wiley & Sons Ltd.

58 GEOSCIENCES↗
NASA NTRS1986

Detailed three-dimensional isentropic analysis using an objective cross-sectional approach

An approach is presented which is designed to give detailed, computationally-efficient meteorological objective analyses over data-rich regions, with results independent of numerical model influences. Observed station data and interstation gradients are used to interpolate pressure, wind, and moisture information from vertical cross-sections onto isentropci surfaces by means of cubic polynomials. Three-dimensional grid-point analyses of mass and moisture are obtained by first constructing a series of nearly parallel vertical cross sections over a specified area, and then filling in the space between adjacent vertical analyses by employing a second set of nearly orthogonal cross sections. Gradients obtained from the mass analyses are then combined with observed wind data to obtain thermally-consistent observed wind analyses. It is concluded that the resulting analyses are useful for preserving significnat horizontal and vertical gradients, such as those associated with jet-streak systems. Examples are given which verify the applicability of this analysis technique to radiosonde data.

Petersen, R. A.↗
NASA NTRS2001

A Two-Dimensional Variational Analysis Method for NSCAT Ambiguity Removal: Methodology, Sensitivity, and Tuning

In this study, we apply a two-dimensional variational analysis method (2d-VAR) to select a wind solution from NASA Scatterometer (NSCAT) ambiguous winds. 2d-VAR determines a "best" gridded surface wind analysis by minimizing a cost function. The cost function measures the misfit to the observations, the background, and the filtering and dynamical constraints. The ambiguity closest in direction to the minimizing analysis is selected. 2d-VAR method, sensitivity and numerical behavior are described. 2d-VAR is compared to statistical interpolation (OI) by examining the response of both systems to a single ship observation and to a swath of unique scatterometer winds. 2d-VAR is used with both NSCAT ambiguities and NSCAT backscatter values. Results are roughly comparable. When the background field is poor, 2d-VAR ambiguity removal often selects low probability ambiguities. To avoid this behavior, an initial 2d-VAR analysis, using only the two most likely ambiguities, provides the first guess for an analysis using all the ambiguities or the backscatter data. 2d-VAR and median filter selected ambiguities usually agree. Both methods require horizontal consistency, so disagreements occur in clumps, or as linear features. In these cases, 2d-VAR ambiguities are often more meteorologically reasonable and more consistent with satellite imagery.

Hoffman, R. N.↗
NASA NTRS1997

Recovery of motion parameters from distortions in scanned images

Scanned images, such as those produced by the scanning-laser ophthalmoscope (SLO), show distortions when there is target motion. This is because pixels corresponding to different image regions are acquired sequentially, and so, in essence, are slices of different snapshots. While these distortions create problems for image registration algorithms, they are potentially useful for recovering target motion parameters at temporal frequencies above the frame rate. Stetter, Sendtner and Timberlake measured large distortions in SLO images to recover the time course of rapid horizontal saccadic eye movements. Here, this work is extended with the goal of automatically recovering small eye movements in two dimensions. Eye position during the frame interval is modeled using a low dimensional parametric description, which in turn is used to generate predicted distortions of a reference template. The input image is then registered to the distorted template using normalized cross correlation. The motion parameters are then varied, and the correlation recomputed, to find the motion which maximizes the peak value of the correlation. The location and value of the correlation maximum are determined with sub-pixel precision using biquadratic interpolation, yielding eye position resolution better than 1 arc minute. This method of motion parameter estimation is tested using actual SLO images as well as simulated images. Motion parameter estimation might also be applied to individual video lines in order to reduce pipeline delays for a near real-time system.

Mulligan, Jeffrey B.↗
NASA NTRS1989

Grid generation and inviscid flow computation about aircraft geometries

Grid generation and Euler flow about fighter aircraft are described. A fighter aircraft geometry is specified by an area ruled fuselage with an internal duct, cranked delta wing or strake/wing combinations, canard and/or horizontal tail surfaces, and vertical tail surfaces. The initial step before grid generation and flow computation is the determination of a suitable grid topology. The external grid topology that has been applied is called a dual-block topology which is a patched C (exp 1) continuous multiple-block system where inner blocks cover the highly-swept part of a cranked wing or strake, rearward inner-part of the wing, and tail components. Outer-blocks cover the remainder of the fuselage, outer-part of the wing, canards and extend to the far field boundaries. The grid generation is based on transfinite interpolation with Lagrangian blending functions. This procedure has been applied to the Langley experimental fighter configuration and a modified F-18 configuration. Supersonic flow between Mach 1.3 and 2.5 and angles of attack between 0 degrees and 10 degrees have been computed with associated Euler solvers based on the finite-volume approach. When coupling geometric details such as boundary layer diverter regions, duct regions with inlets and outlets, or slots with the general external grid, imposing C (exp 1) continuity can be extremely tedious. The approach taken here is to patch blocks together at common interfaces where there is no grid continuity, but enforce conservation in the finite-volume solution. The key to this technique is how to obtain the information required for a conservative interface. The Ramshaw technique which automates the computation of proportional areas of two overlapping grids on a planar surface and is suitable for coding was used. Researchers generated internal duct grids for the Langley experimental fighter configuration independent of the external grid topology, with a conservative interface at the inlet and outlet.

Smith, Robert E.↗
NASA NTRS1996

Four-Dimensional Oceanic and Atmosperic Data Assimilation with Tropical Rainfall Measuring Mission Data

An oceanic data assimilation system which allows to utilize the forthcoming Tropical Rainfall Measuring Mission (TRMM) data has been developed and applied to the Pacific Ocean to produce the velocity field. The assimilated data will be indispensable to examine the effects of rainfall and its variability on the structure and circulation of the tropical oceans and to assess the impact of global warming due to the increase of carbon dioxide on the ocean circulation system and the marine pollution caused by oil spill and ocean damping of radionuclide. The data will also provide the verification for the oceanic and ocean-atmosphere coupled General Circulation Models (GCM's). The system consists of oceanic GCM, analysis scheme and data. In the system the flow field has been determined to be physically consistent with the observed density field and the sea surface winds derived from the Special Sensor Microwave Imagery (SSM/I) data which drive the ocean current. The time integration has been performed for five years until the flow field near the surface attained the steady state starting from the rest ocean with observed temperature and salinity fields, and the SSM/I surface wind velocity. The resultant flow field showed high producibility of the system. Especially the flow near the ocean surface agreed well with available observed data. The system, for the first time, succeeded to produce the eastward subtropical current which has been discovered in the joint investigation on Kuroshio current (CSK) in the 1960s. To verify the quality of the flow field a trajectory analysis has been carried out and compared with the Algos buoy data. BRIEF DESCRIPTION OF THE DATA ASSIMILATION SYSTEM ## Oceanic GCM and analysis scheme--The basic equations are much the same as used for the GCM's, except for the Newtonian damping terms introduced into the prediction equations for the potential temperature and salinity to maintain these fields as observed. The C grid of 2'lat. by 2'long. in horizontal and the 11 vertical levels are applied to the entire Pacific Ocean. At the east and west ocean boundaries the periodic boundary conditions are applied creating fictitious ocean there. The SMAC Method is used to increase the accuracy of mass conservation. * Data--The JODC temperature and salinity data obtained from 1906 to 1988 are used in the system between Long.100'E. and 60'W. The surface wind data are derived from the SSM/I data by Dr-R. Atlas of NASA/GSFC. The data set contains every 6 hours data from July 1987 to June 1989 on the grid of 2'lat. by 2.5'long. The averaged for the whole period and then interpolated into the 2'lat. by 2'long. grid data are used to force the system. The sea bottom topography data was based on the General Bathymetric Chart of the Ocean (GEBCO) supplied by the Canadian Hydrographic Service under contract with the International Hydrographic Organization and International Oceanographic Commission of UNESCO.

Takano, Kenji↗