Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interpolation models”

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 379 records · Page 21

Dynamical Interpolation of Mesoscale Flows in the TOPEX/Poseidon Diamond Surrounding the U.S. Joint Global Ocean Flux Study Bermuda Atlantic Time-Series Study Site

An open boundary ocean model is configured in a domain bounded by the four TOPEX/Poseidon (T/P) ground tracks surrounding the US Joint Global Ocean Flux Study Bermuda Atlantic Time-Series Study (BATS) site. This implementation facilitates prescription of model boundary conditions directly from altimetric measurements (both TIP and ERS-2). The expected error characteristics for a domain of this size with periodically updated boundary conditions are established with idealized numerical experiments using simulated data. A hindcast simulation is then constructed using actual altimetric observations during the period October 1992 through September 1998. Quantitative evaluation of the simulation suggests significant skill. The correlation coefficient between predicted sea level anomaly and ERS observations in the model interior is 0.89; that for predicted versus observed dynamic height anomaly based on hydrography at the BATS site is 0.73. Comparison with the idealized experiments suggests that the main source of error in the hindcast is temporal undersampling of the boundary conditions. The hindcast simulation described herein provides a basis for retrospective analysis of BATS observations in the context of the mesoscale eddy field.

McGillicuddy, Dennis J., Jr.↗

Finite element solver for 3-D compressible viscous flows

The space shuttle main engine (SSME) has extremely complex internal flow structure. The geometry of the flow domain is three-dimensional with complicated topology. The flow is compressible, viscous, and turbulent with large gradients in flow quantities and regions of recirculations. The analysis of the flow field in SSME involves several tedious steps. One is the geometrical modeling of the particular zone of the SSME being studied. Accessing the geometry definition, digitalizing it, and developing surface interpolations suitable for an interior grid generator require considerable amount of manual labor. There are several types of grid generators available with some general-purpose finite element programs. An efficient and robust computational scheme for solving 3D Navier-Stokes equations has to be implemented. Post processing software has to be adapted to visualize and analyze the computed 3D flow field. The progress made in a project to develop software for the analysis of the flow is discussed. The technical approach to the development of the finite element scheme and the relaxation procedure are discussed. The three dimensional finite element code for the compressible Navier-Stokes equations is listed.

Reddy, K. C.↗

Computation Of Flow About A Helicopter Rotor

Improved method developed to simulate numerically flow about helicopter rotor with multiple blades. Computes vortical wake beneath rotor without using simplified ad hoc mathematical models of effects in wake. Patched grids used in numerical simulation of flow about rotor blades. Cylindrical configuration simplifies interpolation at interface between rotating and stationary grids.

Chen, C. L.↗

Parameter estimation in space systems using recurrent neural networks

The identification of time-varying parameters encountered in space systems is addressed, using artificial neural systems. A hybrid feedforward/feedback neural network, namely a recurrent multilayer perception, is used as the model structure in the nonlinear system identification. The feedforward portion of the network architecture provides its well-known interpolation property, while through recurrency and cross-talk, the local information feedback enables representation of temporal variations in the system nonlinearities. The standard back-propagation-learning algorithm is modified and it is used for both the off-line and on-line supervised training of the proposed hybrid network. The performance of recurrent multilayer perceptron networks in identifying parameters of nonlinear dynamic systems is investigated by estimating the mass properties of a representative large spacecraft. The changes in the spacecraft inertia are predicted using a trained neural network, during two configurations corresponding to the early and late stages of the spacecraft on-orbit assembly sequence. The proposed on-line mass properties estimation capability offers encouraging results, though, further research is warranted for training and testing the predictive capabilities of these networks beyond nominal spacecraft operations.

Parlos, Alexander G.↗

Modification of a Successive Corrections Objective Analysis for Improved Derivative Calculations

The use of objectively analyzed fields of meteorological data for complex diagnostic studies and for the initialization of numerical prediction models places the requirements upon the objective method that derivatives of the gridded fields be accurate and free from interpolation error. A modification of an objective analysis developed by Barnes provides improvements in analyses of both the field and its derivatives. Theoretical comparisons, comparisons between analyses of analytical monochromatic waves, and comparisons between analyses of actual weather data are used to show the potential of the new method. The new method restores more of the amplitudes of desired wavelengths while simultaneously filtering more of the amplitudes of undesired wavelengths. These results also hold for the first and second derivatives calculated from the gridded fields. Greatest improvements were for the Laplacian of the height field; the new method reduced the variance of undesirable very short wavelengths by 72 pct. Other improvements were found in the divergence of the gridded wind field and near the boundaries of the field data.

Achtemeier, Gary L.↗

An Integrated Approach to Swept Wing Icing Simulation

This presentation describes the various elements of a simulation approach used to develop a database of ice shape geometries and the resulting aerodynamic performance data for a representative commercial transport wing model exposed to a variety of icing conditions. Methods for capturing full three-dimensional ice shape geometries, geometry interpolation along the span of the wing, and creation of artificial ice shapes based upon that geometric data were developed for this effort. The icing conditions used for this effort were representative of actual ice shape encounter scenarios and run the gamut from ice roughness to full three-dimensional scalloped ice shapes.

aircraft icing↗

First-Order-hold interpolation digital-to-analog converter with application to aircraft simulation

Those who design piloted aircraft simulations must contend with the finite size and speed of the available digital computer and the requirement for simulation reality. With a fixed computational plant, the more complex the model, the more computing cycle time is required. While increasing the cycle time may not degrade the fidelity of the simulated aircraft dynamics, the larger steps in the pilot cue feedback variables (such as the visual scene cues), may be disconcerting to the pilot. The first-order-hold interpolation (FOHI) digital-to-analog converter (DAC) is presented as a device which offers smooth output, regardless of cycle time. The Laplace transforms of these three conversion types are developed and their frequency response characteristics and output smoothness are compared. The FOHI DAC exhibits a pure one-cycle delay. Whenever the FOHI DAC input comes from a second-order (or higher) system, a simple computer software technique can be used to compensate for the DAC phase lag. When so compensated, the FOHI DAC has (1) an output signal that is very smooth, (2) a flat frequency response in frequency ranges of interest, and (3) no phase error. When the input comes from a first-order system, software compensation may cause the FOHI DAC to perform as an FOHE DAC, which, although its output is not as smooth as that of the FOHI DAC, has a smoother output than that of the ZOH DAC.

Cleveland, W. B.↗

Comparison of three atmospheric correction models for a vegetated airborne visible/infrared imaging spectrometer (AVIRIS) scene

Current atmospheric correction models applied to imaging spectroscopy data include such methods as residual or scene average, flat field correction, regression method or empirical line algorithm, the continuum interpolated band ratio (CIBR) derivation and the LOWTRAN 7 method. Due to the limitations of using residual and flat field corrections on vegetated scenes, three methods will be compared: regression, CIBR derivation and LOWTRAN 7. Field-measured bright and dark targets taken at the time of the 13 April, 1989 AVIRIS overflight of Jasper Ridge, California were used to formulate the regression method atmospheric correction. Using this corrected scene as 'ground truth', the CIBR derivation and the LOWTRAN 7 method with both input models are compared on the vegetated Jasper Ridge scene. Although representing a qualitative approach, this is a first approximation and shows the need for more quantitative analysis.

Van Den Bosch, J. M.↗

Dynamical Interpolation of Mesoscale Flows in the TOPEX/ Poseidon Diamond Surrounding the U.S. Joint Global Ocean Flux Study Bermuda Atlantic Time-Series Study Site

An open boundary ocean model is configured in a domain bounded by the four TOPEX/Poseidon (TIP) ground tracks surrounding the U.S. Joint Global Ocean Flux Study Bermuda Atlantic Time-series Study (BATS) site. This implementation facilitates prescription of model boundary conditions directly from altimetric measurements (both TIP and ERS-2). The expected error characteristics for a domain of this size with periodically updated boundary conditions are established with idealized numerical experiments using simulated data. A hindcast simulation is then constructed using actual altimetric observations during the period October 1992 through September 1998. Quantitative evaluation of the simulation suggests significant skill. The correlation coefficient between predicted sea level anomaly and ERS observations in the model interior is 0.89; that for predicted versus observed dynamic height anomaly based on hydrography at the BATS site is 0.73. Comparison with the idealized experiments suggests that the main source of error in the hindcast is temporal undersampling of the boundary conditions. The hindcast simulation described herein provides a basis for retrospective analysis of BATS observations in the context of the mesoscale eddy field.

McGillicuddy, D. J.↗

In Situ Trace Element Measurements on Roda and the Origin of Diogenites

The origin of diogenites remains poorly understood. A recent model interprets many diogenites to have been formed from melts that were derived by remelting initial magma ocean cumulates, and these penultimate parent melts were then contaminated by melts derived from remelting of the basaltic (eucritic) crust to form the ultimate diogenite parent melts [1] (hereafter the remelting model). This is a very complicated petrogenesis that has profound implications for the geological evolution of 4 Vesta if correct. This model was developed based on trace element analyses of bulk rock samples that had been leached in acids to remove phosphates; the compositions of the residues were interpreted to be close to those of cumulus orthopyroxenes plagioclase, chromite and olivine [1]. In situ measurements of phases in diogenites can be used to test this model. We have begun a campaign of laser ablation ICP-MS of orthopyroxene grains in diogenites for this purpose. Here we report our first results on one diogenite, Roda. We have determined a suite of trace lithophile elements on nine, mm-sized pyroxene grains separated from Roda that have previously been studied [2, 3]. A key observation supporting the remelting model is the very low Eu/Eu* of leached residues; values too low to represent orthopyroxene that crystallized from melts with chondritic Sm/Eu and Gd/Eu [1]. (Eu* = Eu interpolated from REE diagrams.) Crustal remelts have low Sm/Eu and Gd/Eu, and orthopyroxenes that crystallized from parent melts contaminated by them would have very low Eu/Eu* [1]. Roda grains have Eu/Eu* of 0.243 to 0.026; the latter a value lower than any measured on bulk diogenite leached residues (0.041) [1]. There is a general negative correlation between Eu/Eu* and some incompatible elements (Zr, Nb, Hf), but not others (LREE). This appears inconsistent with the remelting model as it would suggest an evolving parent melt with La de-creasing as Zr increased and Eu/Eu* decreased. Grain R-15 includes trace-element-rich trapped melt phases [2, 3]. This grain has the highest Eu/Eu* and LREE contents, indicating that the trapped melt had a high Eu/Eu*. Thus, our first data on one diogenite do not provide support for the remelting model [1]. Roda is unusual in that its orthopyroxene grains show wide ranges in trace element contents [4]. Previous in situ REE analyses of grain R-15 did not reveal evidence for subsolidus equilibration with trace-element-rich trapped melt phases, and led to the suggestion that Roda may be polymict, with different grains representing different lithologies of diverse compositions [3]. Thus, based on our results on Roda, it is perhaps premature to abandon the remelting model. In situ measurements on a suite of diogenites is planned to further address this issue.

Mittlefehldt, David W.↗

A calibration model for a stellar photometer using a SEC vidicon.

Description of a SEC vidicon response model that was used to calculate intensities of stars in the ultraviolet. The tube was used only as a stellar filter photometer. The model proposed is applicable for point sources. Experiments with both curve- or surface-fitting techniques and with extrapolation-interpolation procedures were conducted. The curve-fitting techniques were not successful, and the target gain was calculated with the aid of interpolation techniques.

Deutschman, W. A.↗

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.↗

High Accuracy In-Flight Wavelength Calibration of Imaging Spectrometry Data

Accurate wavelength calibration of imaging spectrometer data is essential if proper atmospheric transmission corrections are to be made to obtain apparent surface reflectance. Accuracies of 0.1 nm are necessary for a 10 nm-sampling instrument in order to match the slopes of the deep atmospheric water vapor features that dominate the 0.7-2.3 micrometer regions. The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) is calibrated in the laboratory to determine the wavelength position and full-width-half-maximum (FWHM) response for each of the 224 channels. The accuracies are limited by the quality of the monochromator used as a source. The accuracies vary from plus or minus to plus or minus 1.5 nm depending on the wavelength region, in general decreasing with increasing wavelength. Green et al. make corrections to the wavelength calibrations by using the known positions of 14 atmospheric absorption features throughout the 0.4-2.5 micrometer wavelength region. These features, having varying width and intensity, were matched to the MODTRAN model with a non-linear least squares fitting algorithm. A complete calibration was developed for all 224 channels by interpolation. Instrument calibration cannot be assumed to be stable to 0.1 nm over a flight season given the rigors of thermal cycling and launch and landing loads. The upcoming sensor HYDICE will require a means for in-flight spectral calibration of each scene because the calibration is both temperature and pressure sensitive. In addition, any sensor using a two-dimensional array has the potential for systematic wavelength shifts as a function of cross-track position, commonly called 'smile'. Therefore, a rapid means for calibrating complete images will be required. The following describes a method for determining instrument wavelength calibration using atmospheric absorption features that is efficient enough to be used for entire images on workstations. This study shows the effect of the surface reflectance on the calibration accuracy and the calibration history for the AVIRIS B spectrometer over the 1992 flight season.

Goetz, Alexander F. H.↗

Measurement of Geopotential Heights by GPS Radio Occultation

Geopotential heights of constant pressure surfaces are retrieved from global positioning system (GPS) radio occultation data. In order to assess accuracy a subset of data obtained by GPS/MET during spring 1995 and summer 1995 are compared to the output of the European Centre for Medium-Range Weather Forecasts (ECMWF) global model. The root-mean-square measurement error is 20 m throughout the upper troposphere and lower stratosphere. Furthermore, the ECMWF global model contains enhanced errors in the southeast Pacific. In probing the data for potential utility in climate studies, a Bayesian interpolation technique is used to map the geopotential height fields in the upper troposphere during the summer. Despite limitations of the GPS/MET data set the global average 300-mbar geopotential height over a 2-week period in summer 1995 is determined with an accuracy of 7 m. By obtaining greater coverage and partially resolving synoptic variability, a future constellation of 16 orbiting receivers could obtain global average geopotential height estimates in the upper troposphere with an accuracy of 1 m each day. Accuracy would be somewhat worse for regional studies, except in the tropics where synoptic variability is depressed.

Leroy, Stephen S.↗

Discovery of a Similar to 5 Day Characteristic Timescale in the Kepler Power Spectrum of Zw 229-15

We present time series analyses of the full Kepler dataset of Zw 229− 15. This Kepler light curve- with a baseline greater than three years, composed of virtually continuous, evenly sampled 30-minute measurements - is unprecedented in its quality and precision. We utilize two methods of power spectral analysis to investigate the optical variability and search for evidence of a bend frequency associated with a characteristic optical variability timescale. Each method yields similar results. The first interpolates across data gaps to use the standard Fourier periodogram. The second, using the CARMA-based time-domain modeling technique of Kelly et al., does not need evenly-sampled data. Both methods find excess power at high frequencies that may be due to Kepler instrumental effects. More importantly both also show strong bends (delta alpha is approx. 2) at timescales of approx. 5 days, a feature similar to those seen in the X-ray PSDs of AGN but never before in the optical. This observed approx. 5 day timescale may be associated with one of several physical processes potentially responsible for the variability. A plausible association could be made with light -crossing, dynamical or thermal timescales, depending on the assumed value of the accretion disk size and on unobserved disk parameters such as alpha and H¬R. This timescale is not consistent with the viscous timescale, which would be years in a approx. 10(exp7) solar mass AGN such as Zw 229− 15. However there must be a second bend on long (& 1 year) timescales, and that feature could be associated with the viscous timescale.

timescale↗

Flight Trajectory Prediction Based on Hybrid-Recurrent Networks

The development of future technologies for the National Airspace System (NAS) will be reliant on a new communications infrastructure capable of managing the limited available spectrum for communications among aircraft and ground systems. Emerging approaches to autonomous allocation of aviation spectrum mostlyrely on machine learning techniques, where 4D (longitude, latitude, altitude, time) trajectory prediction is an important data input to enable real-time resource allocation. This study explores and evaluates effective data sources and deep recurrent neural network techniques when determining flight trajectories. Specifically, data are collected and evaluated in a 100-day and 14-day period. Sources of data include NASA Sherlock Data Warehouse, MIT Lincoln Labs Corridor Integrated Weather Service (CIWS), and assorted NOAA weather datasets. Deep learning models for 4D predictions all utilize a hybrid-recurrent technique. A baseline model is considered via the convolutional-LSTM design from the existing literature. The modified design considers Gated Recurrent Units (GRU), Independently Recurrent Neural Networks (IndRNN), and stand-alone self-attention layers. Results indicatethe effectiveness of LSTM and GRUcells for state-of-the-art data processing (interpolation). Additionally, GRUs may be quickly trained with limited data, allowing for exacting improvements with optimizer selection. Attention mechanisms provide notable performance improvements to convolutional layers and may extend dimensional capabilities of a learning model. Finally, NOAA measurements provide only a supplemental value, requiring support from tailored measurements for Air Traffic Management.

Nathan Schimpf↗

Simulating Virtual Terminal Area Weather Data Bases for Use in the Wake Vortex Avoidance System (Wake VAS) Prediction Algorithm

During the research project, sounding datasets were generated for the region surrounding 9 major airports, including Dallas, TX, Boston, MA, New York, NY, Chicago, IL, St. Louis, MO, Atlanta, GA, Miami, FL, San Francico, CA, and Los Angeles, CA. The numerical simulation of winter and summer environments during which no instrument flight rule impact was occurring at these 9 terminals was performed using the most contemporary version of the Terminal Area PBL Prediction System (TAPPS) model nested from 36 km to 6 km to 1 km horizontal resolution and very detailed vertical resolution in the planetary boundary layer. The soundings from the 1 km model were archived at 30 minute time intervals for a 24 hour period and the vertical dependent variables as well as derived quantities, i.e., 3-dimensional wind components, temperatures, pressures, mixing ratios, turbulence kinetic energy and eddy dissipation rates were then interpolated to 5 m vertical resolution up to 1000 m elevation above ground level. After partial validation against field experiment datasets for Dallas as well as larger scale and much coarser resolution observations at the other 8 airports, these sounding datasets were sent to NASA for use in the Virtual Air Space and Modeling program. The application of these datasets being to determine representative airport weather environments to diagnose the response of simulated wake vortices to realistic atmospheric environments. These virtual datasets are based on large scale observed atmospheric initial conditions that are dynamically interpolated in space and time. The 1 km nested-grid simulated datasets providing a very coarse and highly smoothed representation of airport environment meteorological conditions. Details concerning the airport surface forcing are virtually absent from these simulated datasets although the observed background atmospheric processes have been compared to the simulated fields and the fields were found to accurately replicate the flows surrounding the airport where coarse verification data were available as well as where airport scale datasets were available.

Kaplan, Michael L.↗

Applications of remote sensing to stream discharge predictions

A feasibility study has been initiated on the use of remote earth observations for augmenting stream discharge prediction for the design and/or operation of major reservoir systems, pumping systems and irrigation systems. The near-term objectives are the interpolation of sparsely instrumented precipitation surveillance networks and the direct measurement of water loss by evaporation. The first steps of the study covered a survey of existing reservoir systems, stream discharge prediction methods, gage networks and the development of a self-adaptive variation of the Kentucky Watershed model, SNOPSET, that includes snowmelt. As a result of these studies, a special three channel scanner is being built for a small aircraft, which should provide snow, temperature and water vapor maps for the spatial and temporal interpolation of stream gages.

Krause, F. R.↗