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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
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

A Kalman filter for a two-dimensional shallow-water model

A two-dimensional Kalman filter is described for data assimilation for making weather forecasts. The filter is regarded as superior to the optimal interpolation method because the filter determines the forecast error covariance matrix exactly instead of using an approximation. A generalized time step is defined which includes expressions for one time step of the forecast model, the error covariance matrix, the gain matrix, and the evolution of the covariance matrix. Subsequent time steps are achieved by quantifying the forecast variables or employing a linear extrapolation from a current variable set, assuming the forecast dynamics are linear. Calculations for the evolution of the error covariance matrix are banded, i.e., are performed only with the elements significantly different from zero. Experimental results are provided from an application of the filter to a shallow-water simulation covering a 6000 x 6000 km grid.

Parrish, D. F.↗
NASA NTRS2015

Spectral Analysis of Forecast Error Investigated with an Observing System Simulation Experiment

The spectra of analysis and forecast error are examined using the observing system simulation experiment (OSSE) framework developed at the National Aeronautics and Space Administration Global Modeling and Assimilation Office (NASAGMAO). A global numerical weather prediction model, the Global Earth Observing System version 5 (GEOS-5) with Gridpoint Statistical Interpolation (GSI) data assimilation, is cycled for two months with once-daily forecasts to 336 hours to generate a control case. Verification of forecast errors using the Nature Run as truth is compared with verification of forecast errors using self-analysis; significant underestimation of forecast errors is seen using self-analysis verification for up to 48 hours. Likewise, self analysis verification significantly overestimates the error growth rates of the early forecast, as well as mischaracterizing the spatial scales at which the strongest growth occurs. The Nature Run-verified error variances exhibit a complicated progression of growth, particularly for low wave number errors. In a second experiment, cycling of the model and data assimilation over the same period is repeated, but using synthetic observations with different explicitly added observation errors having the same error variances as the control experiment, thus creating a different realization of the control. The forecast errors of the two experiments become more correlated during the early forecast period, with correlations increasing for up to 72 hours before beginning to decrease.

OSSE↗
NASA NTRS2024

Impact Response of Thermoplastic Composites – Experiments and Modeling

This study investigates the impact response of thermoplastic composite materials. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK reinforced with T700G (T700/LMPAEK ). Several High Energy Dynamic Impact (HEDI) tests were performed wherein an aluminum projectile impacted T700/LMPAEK panels at different velocities. In the impact experiment studied herein, the projectile impacted the T700/LMAPEK panel at 73.2 m/s and rebounded. The ensuing damage to the T700/LMPAEK panel was characterized via ultrasonic C-scans, visualizing damage modes including matrix cracking and delaminations. Since the projectile rebounded, the fiber breakage damage mode was not evident. The impact event was modeled using a combined continuum damage mechanics and cohesive zone analysis approach with the commercially available hydrocode LS-DYNA®. The dynamic deformation and damage of the continuum, i.e., the plies in the T700/LMPAEK laminate, was modeled using a Deformation Gradient Decomposition (DGD) material model, i.e., MAT299 available within LS-DYNA. The interlaminar regions were treated as zero-thickness cohesive zones, wherein the dynamic delamination mechanics was modeled using a traction-separation formulation. Mixed-mode delamination growth was captured using the Benzegaggh-Kenane (B-K) traction-separation interpolation law. Strain-rate sensitivity of interlaminar fracture toughness was included using approximations based on published literature. The outcomes of the modeling effort are discussed herein which demonstrate good correlations with test data regarding panel deflection, panel-projectile interactions, damage modes and damage extent.

Impact modeling↗
NASA NTRS1990

Analytical aerodynamic model of a high alpha research vehicle wind-tunnel model

A 6 DOF analytical aerodynamic model of a high alpha research vehicle is derived. The derivation is based on wind-tunnel model data valid in the altitude-Mach flight envelope centered at 15,000 ft altitude and 0.6 Mach number with Mach range between 0.3 and 0.9. The analytical models of the aerodynamics coefficients are nonlinear functions of alpha with all control variable and other states fixed. Interpolation is required between the parameterized nonlinear functions. The lift and pitching moment coefficients have unsteady flow parts due to the time range of change of angle-of-attack (alpha dot). The analytical models are plotted and compared with their corresponding wind-tunnel data. Piloted simulated maneuvers of the wind-tunnel model are used to evaluate the analytical model. The maneuvers considered are pitch-ups, 360 degree loaded and unloaded rolls, turn reversals, split S's, and level turns. The evaluation finds that (1) the analytical model is a good representation at Mach 0.6, (2) the longitudinal part is good for the Mach range 0.3 to 0.9, and (3) the lateral part is good for Mach numbers between 0.6 and 0.9. The computer simulations show that the storage requirement of the analytical model is about one tenth that of the wind-tunnel model and it runs twice as fast.

Cao, Jichang↗
NASA NTRS1994

The scale and persistence of soil moisture anomalies as simulated in a global model

Short term variability of climate is intimately connected with soil moisture variability. Soil moisture provides the storage and subsequent return to the atmosphere, through evaporation and transpiration, of precipitation anomalies over land. Global Circulation Model (GCM) simulations enable consistent identification of correlations and dynamical connections between the hydrologic variables, many of which are incompletely observed. One way to facilitate understanding with these increasingly intricate models is to perform sensitivity studies in which a boundary condition or process is prescribed. In this study we will report on a sensitivity study in which a GCM with a sophisticated land surface representation is used to investigate soil moisture variability in the model climate. The simulations to be used in this study were made at R15 resolution (approximately 4.5 deg latitude x 7.5 deg longitude) with prescribed sea surface temperatures (SST) in the GENESIS model (Thompson and Pollard, 1994), which is coupled to a Land Surface Transfer model (LSX) at 2 deg x 2 deg resolution (Pollard and Thompson, 1994). All the results represented here were taken from the monthly averages of the model results. The LSX model accounts for the physical effects of vegetation with two layers specified at each grid point. Vegetation attributes such as leaf area indices, fractional cover, leaf albedos, etc., were taken from the global dataset in Dorman and Sellers (1989). A six-layer soil model extends from the surface to 4.25 m depth. SST's were prescribed in two ten year experiments using monthly SST values with the daily value being interpolated from the nearest two months. In the first experiment monthly climatological values were used, and in the second, the Atmospheric Model Intercomparison Project (AMIP) observed SST's for the years 1979 through 1988 were used (Gates, 1992). Thus, the former experiment gives a measure of the intrinsic model variability, to be compared with that of the latter experiment, which includes month-to-month variability due to ocean forcing.

Fitzjarrald, D.↗
NASA NTRS1989

Six-degree-of-freedom aircraft simulation with mixed-data structure using the applied dynamics simulation language, ADSIM

A realistic simulation of an aircraft in the flight using the AD 100 digital computer is presented. The implementation of three model features is specifically discussed: (1) a large aerodynamic data base (130,00 function values) which is evaluated using function interpolation to obtain the aerodynamic coefficients; (2) an option to trim the aircraft in longitudinal flight; and (3) a flight control system which includes a digital controller. Since the model includes a digital controller the simulation implements not only continuous time equations but also discrete time equations, thus the model has a mixed-data structure.

Savaglio, Clare↗
NASA NTRS1984

Interpolation problems in meteor radar analysis, part 7.6A

Meteor echoes come from random points in the observation volume, and are irregularly spaced in time. This precludes the use of fast fourier transformations (FFT) techniques on the raw data to give the spectrum of waves that are present. One way around this obstacle is to restrict our interest to a particular class of waves, and fit a corresponding model to the raw data. It is assumed that there is no horizontal variation across the observation volume for tides, but in the vertical this is certainly not the case. If, in addition, we are interested in other types of waves which may be present and whose periods are unknown, then examining the raw line-of-sight velocities does not tell us how to modify the model, since the line-of-sight direction is not fixed. This is the motivation for interpolation. Interpolation takes a temporal series of line-of-sight velocities, and transforms it to a temporal series of wind velocities for each orthogonal direction. The velocities along a given direction can then be examined readily for any waves in addition to tides.

Tetenbaum, D.↗
NASA NTRS2012

Projections of Declining Surface-Water Availability for the Southwestern United States

16 of the CMIP5 models had all the data needed for this work for at least one simulation that was continuous from 1950 to 2040. Details of the models analyzed here are provided in Table S1. The model data analyzed here are available at http://strega.ldeo.columbia.edu:81/expert/home/.naomi/.AR5/.v2/.historical:rcp85/.mmm16/ a. Assessing the climatology of the models Despite increases in horizontal resolution of many models compared to their CMIP3 counterparts none of these models can adequately resolve the topography of the south west United States, such as the Sierra Nevada and Rocky Mountains and the associated orographic precipitation. This requires that caution be used when interpreting the results presented here. To assess the ability of the models to simulate the current hydroclimate, in Figure S1 we show the observed (from the Global Precipitation Climatology Centre gridded rain gauge data, (1)) monthly climatology of precipitation and the same for all the models and the multimodel mean for the California-Nevada, Colorado headwaters and Texas regions. The GPCC data uses rain gauges only and interpolates to regular grids of which we used the 1◦ by 1◦ one. Details of the data set can be found in (2). While the models apparently overestimate precipitation in California and Nevada the seasonal cycle with wet winters and dry summers is very well represented. It is also possible that the rain gauge observations are biased low by inadequately sampling the higher mountain regions. How ever the models might also be expected to underestimate orographic precipitation due to inadequate horizontal resolution. The 25 models are also too wet in the Colorado headwaters region but correctly represent the quite even distribution though the year. The bimodal distribution of precipitation in Texas, with peaks in May and September, and the absolute amounts, are well modeled but with the September peak too weak. The positive precipitation bias translates into a positive runoff bias for the Colorado headwaters as also shown in Figure S1. Here the observed runoff values are taken from simulations of the Variable Infiltration Capacity (VIC) land surface-hydrology model (3) forced by observed meteorology (5) that were conducted as part of the North American Land Data Assimilation System project phase 2 ( (NLDAS-2), http://www.emc.ncep.noaa.gov/mmb/nldas/. Runoff for California-Nevada is better simulated but there is a positive bias over Texas despite no strong precipitation bias. To check whether regional climate models better simulate P and runoff in these regions we analyzed the historical simulation with the Regional Climate Model version 3 driven by the National Centers for Environmental Prediction-Department of Energy Reanalysis 2 available from the North American Regional Climate Change Assessment Program (http://www.narccap.ucar.edu). This model configuration retained these biases in P and runoff although they were reduced in amplitude. Given these varying biases we plot P and P − E changes in actual values but apply the simplest bias correction possible to the runoff and soil moisture values and show the modeled changes in terms of percentages of the 20th Century model climatologies. A thorough assessment of the simulation of North American climate in CMIP5 models is conducted in Sheffield at al. (North American Climate in CMIP5 Experiments. Part I: Evaluation of 20th Century Continental and Regional Climatology, manuscript submit ted to J. Climate, available at http://www.climate.noaa.gov/index.jsp?pg=./cpo pa/ mapp/cmip5 publications.html). Sheffield et al. analyze the climatology of precipitation, surface air temperature, low level winds, moisture fluxes, runoff etc. and conclude that the main features of the hydrological cycle, including characteristics of the atmospheric moisture balance and its seasonality, are captured in the CMP5 models subject to biases in total precipitation amounts. We chose to use all available models instead of selecting some and rejecting others based on an assessment of model realism. This is in accord with the suggestions of Mote et al. for CMIP3 (4) but future work needs to revisit this matter for the case of the CMIP5 ensemble.

Seager, Richard↗
NASA NTRS1982

Observation of additional low-degree 5-min modes of solar oscillation

High-order solar oscillations with degrees l=3, 4, and 5 could be detected. The observations were made by measuring the difference between the shifts in the Fe 5,124 spectrum line from light integrated from a central circular portion of the solar disk and from an annular portion exterior to it. The frequencies of the octupole modes agree well with the values obtained from whole-disk measurements at the South Pole. A least-squares fit of the observed frequencies to values interpolated between and extrapolated from the predictions of a sequence of solar models with different chemical compositions selects two models. One, a helium-rich solution, agrees with that of similar analyses of whole-disk data. The extrapolated solution has a relatively deep convection zone, and is thus consistent with analyses of 5-min oscillations of high degree.

Scherrer, P. H.↗
NASA NTRS1983

Propeller performance and weight predictions appended to the Navy/NASA engine program

The Navy/NASA Engine Performance (NNEP) is a general purpose computer program currently employed by government, industry and university personnel to simulate the thermodynamic cycles of turbine engines. NNEP is a modular program which has the ability to evaluate the performance of an arbitrary engine configuration defined by the user. In 1979, a program to calculate engine weight (WATE-2) was developed by Boeing's Military Division under NASA contract. This program uses a preliminary design approach to determine engine weights and dimensions. Because the thermodynamic and configuration information required by the weight code was available in NNEP, the weight code was appended to NNEP. Due to increased emphasis on fuel economy, a renewed interest has developed in propellers. This report describes the modifications developed by NASA to both NNEP and WATE-2 to determine the performance, weight and dimensions of propellers and the corresponding gearbox. The propeller performance model has three options, two of which are based on propeller map interpolation. Propeller and gearbox weights are obtained from empirical equations which may easily be modified by the user.

Plencner, R. M.↗
NASA NTRS1988

Nonlinear modeling, identification and estimation of slew induced structural deformation

This proposal discusses a new nonlinear, nonparametric method for off-line modeling and on-line estimation of the deformation of a flexible structure undergoing rapid retargeting maneuvers. In these circumstances, the structural stiffness and damping coefficients depend on the angular acceleration omega(dot), the angular rate omega, and the square of the angular rate omega. In the single axis case, the excitation of the structure is represented by the vector u(exp T) = (omega(dot), omega(exp 2), 2(omega)), to which the structural dynamics responds as a 'bilinear' (i.e., parametrically excited) system. A similar technique for multiaxial rotations yields a bilinear model with respect to matrix valued excitations. Three methods of estimation and modeling are described in this proposal to achieve deformation state determination: (1) a method based on a feedback linearized procedure which gives an estimate of the state by means of observers installed in the deformable body; (2) off-line modeling of the deformation state of the structure by means of topological interpolators; and (3) an on-line structural state estimation method based on a combination of the two previous techniques.

Dwyer, Thomas A., III↗
NASA NTRS1993

A finite-element-based perturbation model for the rotordynamic analysis of shrouded pump impellers: Part 2: User's guide

This report describes the computational steps involved in executing a finite-element-based perturbation model for computing the rotor dynamic coefficients of a shrouded pump impeller or a simple seal. These arise from the fluid/rotor interaction in the clearance gap. In addition to the sample cases, the computational procedure also applies to a separate category of problems referred to as the 'seal-like' category. The problem, in this case, concerns a shrouded impeller, with the exception that the secondary, or leakage, passage is totally isolated from the primary-flow passage. The difference between this and the pump problem is that the former is analytically of the simple 'seal-like' configuration, with two (inlet and exit) flow-permeable stations, while the latter constitutes a double-entry / double-discharge flow problem. In all cases, the problem is that of a rotor clearance gap. The problem here is that of a rotor excitation in the form of a cylindrical whirl around the housing centerline for a smooth annular seal. In its centered operation mode, the rotor is assumed to give rise to an axisymmetric flow field in the clearance gap. As a result, problems involving longitudinal or helical grooves, in the rotor or housing surfaces, go beyond the code capabilities. Discarding, for the moment, the pre- and post-processing phases, the bulk of the computational procedure consists of two main steps. The first is aimed at producing the axisymmetric 'zeroth-order' flow solution in the given flow domain. Detailed description of this problem, including the flow-governing equations, turbulence closure, boundary conditions, and the finite-element formulation, was covered by Baskharone and Hensel. The second main step is where the perturbation model is implemented, with the input being the centered-rotor 'zeroth-order' flow solution and a prescribed whirl frequency ratio (whirl frequency divided by the impeller speed). The computational domain, in the latter case, is treated as three dimensional, with the number of computational planes in the circumferential direction being specified a priori. The reader is reminded that the deformations in the finite elements are all infinitesimally small because the rotor eccentricity itself is a virtual displacement. This explains why we have generically termed the perturbation model the 'virtually' deformable finite-element category. The primary outcome of implementing the perturbation model is the tangential and radial components, F(sub theta)(sup *) and F(sub r)(sup *) of the fluid-exerted force on the rotor surface due to the whirling motion. Repetitive execution of the perturbation model subprogram over a sufficient range of whirl frequency ratios, and subsequent interpolation of these fluid forces, using the least-square method, finally enable the user to compute the impeller rotor dynamic coefficients of the fluid/rotor interaction. These are the direct and cross-coupled stiffness, damping, and inertia effects of the fluid/rotor interaction.

Baskharone, Erian A.↗
NASA NTRS1995

An Experimental and Analytical Investigation of Stirling Space Power Converter Heater Head

NASA has identified the Stirling power converter as a prime candidate for the next generation power system for space applications requiring 60000 hr of operation. To meet this long-term goal, several critical components of the power converter have been analyzed using advanced structural assessment methods. Perhaps the most critical component, because of its geometric complexity and operating environment, is the power converter's heater head. This report describes the life assessment of the heater head which includes the characterization of a viscoplastic material model, the thermal and structural analyses of the heater head, and the interpolation of fatigue and creep test results of a nickel-base superalloy, Udimet 720 LI (Low Inclusions), at several elevated temperatures for life prediction purposes.

Abdul-Aziz, Ali↗
NASA NTRS1996

A Multigrid Algorithm for Immersed Interface Problems

Many physical problems involve interior interfaces across which the coefficients in the problem, the solution, its derivatives, the flux, or the source term may have jumps. These interior interfaces may or may not align with a underlying Cartesian grid. Zhilin Li, in his dissertation, showed how to discretize such elliptic problems using only a Cartesian grid and the known jump conditions to second order accuracy. In this paper, we describe how to apply the full multigrid algorithm in this context. In particular, the restriction, interpolation, and coarse grid problem will be described. Numerical results for several model problems are given to demonstrate that good rates can be obtained even when jumps in the coefficients are large and do not align with the grid.

Adams, Loyce↗
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↗
NASA NTRS2007

Water Droplet Impingement on Simulated Glaze, Mixed, and Rime Ice Accretions

Water droplet impingement data were obtained at the NASA Glenn Icing Research Tunnel (IRT) for a 36-in. chord NACA 23012 airfoil with and without simulated ice using a dye-tracer method. The simulated ice shapes were defined with the NASA Glenn LEWICE 2.2 ice accretion program and including one rime, four mixed and five glaze ice shapes. The impingement experiments were performed with spray clouds having median volumetric diameters of 20, 52, 111, 154, and 236 micron. Comparisons to the experimental data were generated which showed good agreement for the rime and mixed shapes at lower drop sizes. For larger drops 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 or shadow regions of ice shapes.

Papadakis, Michael↗
NASA NTRS2021

Temporal Interpolation of Geostationary Satellite Imagery With Optical Flow

Applications of satellite data in areas such as weather tracking and modeling, ecosystem monitoring, wildfire detection, and land-cover change are heavily dependent on the tradeoffs to spatial, spectral, and temporal resolutions of observations. In weather tracking, high-frequency temporal observations are critical and used to improve forecasts, study severe events, and extract atmospheric motion, among others. However, while the current generation of geostationary (GEO) satellites has hemispheric coverage at 10-15-min intervals, higher temporal frequency observations are ideal for studying mesoscale severe weather events. In this work, we present a novel application of deep learning-based optical flow to temporal upsampling of GEO satellite imagery. We apply this technique to 16 bands of the GOES-R/Advanced Baseline Imager mesoscale dataset to temporally enhance full-disk hemispheric snapshots of different spatial resolutions from 10 to 1 min. Experiments show the effectiveness of task-specific optical flow and multiscale blocks for interpolating high-frequency severe weather events relative to bilinear and global optical flow baselines. Finally, we demonstrate strong performance in capturing variability during convective precipitation events.

Image processing↗