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 73 records · Page 4

Interpolating Fields of Carbon Monoxide Data Using a Hybrid Statistical-Physical Model

Atmospheric Carbon Monoxide (CO) is a pollutant gas of which the US congress has mandated regular monitoring, and satellite sensors can be used to retrieve regional concentrations of CO over several vertical layers. However, CO at cloudy locations cannot be observed and have to be estimated from the observed data set, resulting in an interpolation problem. The current state-of-the-art solution is to combine prior information, computed by a deterministic physical model, with observations. However, the deterministic model may introduce uncertainties that do not derive from the data. While sharing certain features with the physical model, this paper presents a Bayesian hierarchical model for interpolating CO on a 3-dimensional spatial grid, across time. To our knowledge such a model has not been considered before. The model is applied to a hypothetical air-quality monitoring scenario, and is compared to existing interpolation methods. The results provide motivation for the use of the statistical model for regional to local applications.

Arellano, A. A.↗

Status of the NASA GMAO Observing System Simulation Experiment

An Observing System Simulation Experiment (OSSE) is a pure modeling study used when actual observations are too expensive or difficult to obtain. OSSEs are valuable tools for determining the potential impact of new observing systems on numerical weather forecasts and for evaluation of data assimilation systems (DAS). An OSSE has been developed at the NASA Global Modeling and Assimilation Office (GMAO, Errico et al 2013). The GMAO OSSE uses a 13-month integration of the European Centre for Medium- Range Weather Forecasts 2005 operational model at T511/L91 resolution for the Nature Run (NR). Synthetic observations have been updated so that they are based on real observations during the summer of 2013. The emulated observation types include AMSU-A, MHS, IASI, AIRS, and HIRS4 radiance data, GPS-RO, and conventional types including aircraft, rawinsonde, profiler, surface, and satellite winds. The synthetic satellite wind observations are colocated with the NR cloud fields, and the rawinsondes are advected during ascent using the NR wind fields. Data counts for the synthetic observations are matched as closely as possible to real data counts, as shown in Figure 2. Errors are added to the synthetic observations to emulate representativeness and instrument errors. The synthetic errors are calibrated so that the statistics of observation innovation and analysis increments in the OSSE are similar to the same statistics for assimilation of real observations, in an iterative method described by Errico et al (2013). The standard deviations of observation minus forecast (xo-H(xb)) are compared for the OSSE and real data in Figure 3. The synthetic errors include both random, uncorrelated errors, and an additional correlated error component for some observational types. Vertically correlated errors are included for conventional sounding data and GPS-RO, and channel correlated errors are introduced to AIRS and IASI (Figure 4). HIRS, AMSU-A, and MHS have a component of horizontally correlated error. The forecast model used by the GMAO OSSE is the Goddard Earth Observing System Model, Version 5 (GEOS-5) with Gridpoint Statistical Interpolation (GSI) DAS. The model version has been updated to v. 5.13.3, corresponding to the current operational model. Forecasts are run on a cube-sphere grid with 180 points along each edge of the cube (approximately 0.5 degree horizontal resolution) with 72 vertical levels. The DAS is cycled at 6-hour intervals, with 240 hour forecasts launched daily at 0000 UTC. Evaluation of the forecasting skill for July and August is currently underway. Prior versions of the GMAO OSSE have been found to have greater forecasting skill than real world forecasts. It is anticipated that similar forecast skill will be found in the updated OSSE.

OSSE↗

Modeling Martian Dust Using Mars-GRAM

Engineering-level atmospheric model widely used for diverse mission applications. Mars-GRAM s perturbation modeling capability is commonly used, in a Monte-Carlo mode, to perform high fidelity engineering end-to-end simulations for entry, descent, and landing (EDL). From the surface to 80 km altitude, Mars-GRAM is based on NASA Ames Mars General Circulation Model (MGCM). Mars-GRAM and MGCM use surface topography from Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA), with altitudes referenced to the MOLA areoid, or constant potential surface. Traditional Mars-GRAM options for representing the mean atmosphere along entry corridors include: TES Mapping Years 1 and 2, with Mars-GRAM data coming from MGCM model results driven by observed TES dust optical depth TES Mapping Year 0, with user-controlled dust optical depth and Mars-GRAM data interpolated from MGCM model results driven by selected values of globally-uniform dust optical depth. Mars-GRAM 2005 has been validated against Radio Science data, and both nadir and limb data from the Thermal Emission Spectrometer (TES).

Justh, Hilary L.↗

Utilizing Mars Global Reference Atmospheric Model (Mars-GRAM 2005) to Evaluate Entry Probe Mission Sites

Engineering-level atmospheric model widely used for diverse mission applications. Mars-GRAM s perturbation modeling capability is commonly used, in a Monte-Carlo mode, to perform high fidelity engineering end-to-end simulations for entry, descent, and landing (EDL)1. Traditional Mars-GRAM options for representing the mean atmosphere along entry corridors include: a) TES Mapping Years 1 and 2, with Mars-GRAM data coming from MGCM model results driven by observed TES dust optical depth; and b) TES Mapping Year 0, with user-controlled dust optical depth and Mars-GRAM data interpolated from MGCM model results driven by selected values of globally-uniform dust optical depth. From the surface to 80 km altitude, Mars-GRAM is based on NASA Ames Mars General Circulation Model (MGCM). Mars-GRAM and MGCM use surface topography from Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA), with altitudes referenced to the MOLA areoid, or constant potential surface. Mars-GRAM 2005 has been validated2 against Radio Science data, and both nadir and limb data from the Thermal Emission Spectrometer (TES)

Justh, Hilary L.↗

Evaluating Mars Science Laboratory Landing Sites with the Mars Global Reference Atmospheric Model (Mars-GRAM 2005)

The Mars Global Reference Atmospheric Model (Mars-GRAM) is an engineering-level atmospheric model widely used for diverse mission applications. Mars-GRAM s perturbation modeling capability is commonly used, in a Monte-Carlo mode, to perform high fidelity engineering end-to-end simulations for entry, descent, and landing (EDL) [1]. From the surface to 80 km altitude, Mars-GRAM is based on the NASA Ames Mars General Circulation Model (MGCM). Mars-GRAM and MGCM use surface topography from Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA), with altitudes referenced to the MOLA areoid, or constant potential surface. Traditional Mars-GRAM options for representing the mean atmosphere along entry corridors include: (1) Thermal Emission Spectrometer (TES) mapping years 1 and 2, with Mars-GRAM data coming from NASA Ames Mars General Circulation Model (MGCM) results driven by observed TES dust optical depth or (2) TES mapping year 0, with user-controlled dust optical depth and Mars-GRAM data interpolated from MGCM model results driven by selected values of globally-uniform dust optical depth. Mars-GRAM 2005 has been validated [2] against Radio Science data, and both nadir and limb data from TES [3]. There are several new features included in Mars-GRAM 2005. The first is the option to use input data sets from MGCM model runs that were designed to closely simulate conditions observed during the first two years of TES observations at Mars. The TES Year 1 option includes values from April 1999 through January 2001. The TES Year 2 option includes values from February 2001 through December 2002. The second new feature is the option to read and use any auxiliary profile of temperature and density versus altitude. In exercising the auxiliary profile Mars-GRAM option, values from the auxiliary profile replace data from the original MGCM databases. Some examples of auxiliary profiles include data from TES nadir or limb observations and Mars mesoscale model output at a particular location and time. The final new feature is the addition of two Mars-GRAM parameters that allow standard deviations of Mars-GRAM perturbations to be adjusted. The parameter rpscale can be used to scale density perturbations up or down while rwscale can be used to scale wind perturbations.

Justh, H. L.↗

Kamodo’s Model-Agnostic Satellite Flythrough: Lowering the Utilization Barrier for Heliophysics Model Outputs

Heliophysics model outputs are increasingly accessible, but typically are not usable by the majority of the community unless directly collaborating with the relevant model developers. Prohibitive factors include complex file output formats, cryptic metadata, unspecified and often customized coordinate systems, and non-linear coordinate grids. Some pockets of progress exist, giving interfaces to various simulation outputs, but only for a small set of outputs and typically not with open-source, freely available packages. Additionally, the increasing array of tools built upon these sporadic interfaces are typically model-specific. We present Kamodo’s model-agnostic satellite flythrough capabilities as the solution to the utilization barrier for heliophysics model outputs. Developed at the Community Coordinated Modeling Center, these flythrough capabilities are built in Python upon a network of model-agnostic interfaces developed in collaboration with model developers, providing interpolation results the community can trust. Kamodo’s flythrough capabilities present the user with a growing variety of flythrough tools based upon a rapidly expanding library of heliophysics model outputs in several domains, currently including a variety of Ionosphere-Thermosphere-Mesosphere and global magnetosphere model outputs. Each capability is designed to be easily accessible via simplistic model-agnostic syntax, with the entire package freely available in the cloud on Github. Here, we describe the tools developed, include several sample applications for common science questions, demonstrate interoperability with selected packages, and summarize ongoing developments.

Software↗

Reconstructing irregularly sampled images by neural networks

Neural-network-like models of receptor position learning and interpolation function learning are being developed as models of how the human nervous system might handle the problems of keeping track of the receptor positions and interpolating the image between receptors. These models may also be of interest to designers of image processing systems desiring the advantages of a retina-like image sampling array.

Ahumada, Albert J., Jr.↗

Theory of the lattice Boltzmann Method: Dispersion, Dissipation, Isotropy, Galilean Invariance, and Stability

The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyper-viscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters which optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some 2D LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long wave-length limit (wave vector k = 0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.

Lallemand, Pierre↗

Finite element algorithms for compressible flow computation on a supercomputer

Methods of applying computational fluid dynamics techniques to engineering problems are discussed. The linking factors between theoretical and applied research are the shape and weight functions which produce the Galerkin finite element schemes. The flow is unsteady, viscous, heat conducting and compressible, and steady-state flows are the asymptotic limit of unsteady flows. All flows are described by the time-averaged Navier-Stokes equations (NSE) with mass and energy conservation. Local curvilinear intrinsic coordinates are applied to discretize the NSE in arbitrary geometric domains. A modified weighted residuals approach defines the discrete analogs of the physical systems modeled and a general interpolants method is used to derive families of numerical models, both implicit and explicit, finite difference and finite elemnt, from a single point of departure. Time-dependent element approximations are achieved with a progressive assembly of generalized elements method. Hyperbolic steady-state Euler algorithms and quasi-parabolic are spatial marching algorithms for solving the discretized equations. Application of the techniques is illustrated through calculation of the compressible flow around a hypersonic flight vehicle at 100,000 ft at zero angle of attack at speeds of Mach 4-10.

Spradley, L. W.↗

Estimation of the Ocean Skin Temperature using the NASA GEOS Atmospheric Data Assimilation System

This report documents the status of the development of a sea surface temperature (SST) analysis for the Goddard Earth Observing System (GEOS) Version-5 atmospheric data assimilation system (ADAS). Its implementation is part of the steps being taken toward the development of an integrated earth system analysis. Currently, GEOS-ADAS SST is a bulk ocean temperature (from ocean boundary conditions), and is almost identical to the skin sea surface temperature. Here we describe changes to the atmosphere-ocean interface layer of the GEOS-atmospheric general circulation model (AGCM) to include near surface diurnal warming and cool-skin effects. We also added SST relevant Advanced Very High Resolution Radiometer (AVHRR) observations to the GEOS-ADAS observing system. We provide a detailed description of our analysis of these observations, along with the modifications to the interface between the GEOS atmospheric general circulation model, gridpoint statistical interpolation-based atmospheric analysis and the community radiative transfer model. Our experiments (with and without these changes) show improved assimilation of satellite radiance observations. We obtained a closer fit to withheld, in-situ buoys measuring near-surface SST. Evaluation of forecast skill scores corroborate improvements seen in the observation fits. Along with a discussion of our results, we also include directions for future work.

AGCM↗

Postprocessing for the UCLA climate model

Postprocessing programs written to analyze the results of a general circulation model are described. The model history data is interpolated to pressure surfaces, converted from 64-bit Cyber words to 32 bit IBM words, and transferred to tape for further analysis. The time average program reads the interpolation output tape and produces 15-day time averages. The model outputs a very detailed record for selected grid points. This data is sent to a tape in Cyber full words. The conversion program changes the format of the tape from 64 bit Cyber words to 32 bit IBM words, processing a month of data at a time. The grid-point history sorting program inputs the converted grid point history tape and sorts it by time. The end product is a tape with a time history of each field for every point. After the PBP data has been stored, it can be plotted. The plotting package is flexible, and allows the choice of fields and grid points. The one-dimensional fields are plotted as functions of time, and a daily summary is also produced. For the two-dimensional fields, time-height cross sections and daily summaries are produced. The model plotting package produces plots for all fields on the time average and pressure tapes. These plots include world maps, polar projections, zonal means, vectors, streamlines, and latitude-height plots.

Abeles, J.↗

GMI-IPS: Processing & Visualization Software Used in ATom DC-8 Aircraft Studies

NASA's Atmospheric Tomography Mission (ATom) deployed in each of the four seasons during 2016-2018, the DC-8 aircraft in order to establish global-scale datasets intended to improve the representation of chemically reactive gases in global atmospheric chemistry models (ACMs). The Global Modeling Initiative (GMI) executed simulations for each ATom flight using the GMI Chemistry Transport Model (GMI-CTM) to provide species concentrations of chemical gases along the DC-8 flight transects. To solve the problem of translating the GMI-CTM simulation data to the unique spatial resolutions of each ATom flight, the GMI ICARTT Processing Software (GMI-IPS) was developed.The GMI-IPS is written in Python and provides data processing, flight extraction, and visualization support for aircraft research projects using ICARTT format, which is a standard format for airborne instrument data. Additionally, the GMI-IPS interpolates global gridded model data from Hierarchical Data Format (HDF) to ICARTT compatible flight transects. Software classes for instruments and collections provided by the ATom DC-8 aircraft such as MER10, MMS, etc. are derived from a common base class. Other functionality provided by the GMI-IPS are: deriving missing flight entries along a transect, reading ICARTT entries from file, and providing Python data structures for storing flight and model information, and more.The GMI-IPS is GIT source controlled, has approximately 30,000 lines of code, and supports parallelization across data collections. It delivered GMI-CTM data for more than forty distinct DC-8 aircraft flights that took place under ATom. The output ICARTT files adhere to format standard V1.1, and pass the scan utility provided by NASA LaRC Airborne Science Data for Atmospheric Composition. This presentation will include a software and methods overview, and results from ATom, including assessments using the GMI-CTM showing how well observations from ATom flight transects represent a broader region.

Damon, M. R.↗

Photogrammetric and Global Positioning System Measurements of Active Pahoehoe Lava Lobe Emplacement on Kilauea, Hawaii

Basalt is the most common rock type on the surface of terrestrial bodies throughout the solar system and -- by total volume and areal coverage -- pahoehoe flows are the most abundant form of basaltic lava in subaerial and submarine environments on Earth. A detailed understanding of pahoehoe emplacement processes is necessary for developing accurate models of flow field development, assessing hazards associated with active lava flows, and interpreting the significance of lava flow morphology on Earth and other planetary bodies. Here, we examine the active emplacement of pahoehoe lobes along the margins of the Hook Flow from Pu'u 'O'o on Kilauea, Hawaii. Topographic data were acquired between 21 and 23 February 2006 using stereo-imaging and differential global positing system (DGPS) measurements. During this time, the average discharge rate for the Hook Flow was 0.01-0.05 cubic m/s. Using stereogrammetric point clouds and interpolated digital terrain models (DTMs), active flow fronts were digitized at 1 minute intervals. These areal spreading maps show that the lava lobe grew by a series of breakouts tha t broadly fit into two categories: narrow (0.2-0.6 m-wide) toes that grew preferentially down-slope, and broad (1.4-3.5 m-wide) breakouts that formed along the sides of the lobe, nearly perpendicular to the down-flow axis. These lobes inflated to half of their final thickness within approx 5 minutes, with a rate of inflation that generally deceased with time. Through a combination of down-slope and cross-slope breakouts, lobes developed a parabolic cross-sectional shape within tens of minutes. We also observed that while the average local discharge rate for the lobe was generally constant at 0.0064 +/- 0.0019 cubic m/s, there was a 2 to 6 fold increase in the areal coverage rate every 4.1 +/- 0.6 minutes. We attribute this periodicity to the time required for the dynamic pressurization of the liquid core of the lava lobe to exceed the cooling-induced strength of the lobe margins. Using DGPS-derived DTMs of the topography before and after pahoehoe lobe emplacement, we observed that the lava typically concentrated within existing topographic lows, with the lobe reaching a maximum thickness of approx 1.2 m above the lowest points of the initial topography and above reverse-facing slopes. Lobe margins were typically controlled by high-standing topography, with the zone directly adjacent to the final flow margin having average relief that is approx 4 cm higher than the lava-inundated region. This suggests that irregularities approx 25% of the height of the smallest breakout elements (i.e., toes) can exert a strong control on the paths of low-discharge pahoehoe lobes, with stagnated toes forming confining margins that allow interior portions of flow to topographically invert the landscape by inflation.

Hamilton, Christopher W.↗

Detailed Uncertainty Analysis of the Ares I A106 Liftoff/Transition Database

The Ares I A106 Liftoff/Transition Force and Moment Aerodynamics Database describes the aerodynamics of the Ares I Crew Launch Vehicle (CLV) from the moment of liftoff through the transition from high to low total angles of attack at low subsonic Mach numbers. The database includes uncertainty estimates that were developed using a detailed uncertainty quantification procedure. The Ares I Aerodynamics Panel developed both the database and the uncertainties from wind tunnel test data acquired in the NASA Langley Research Center s 14- by 22-Foot Subsonic Wind Tunnel Test 591 using a 1.75 percent scale model of the Ares I and the tower assembly. The uncertainty modeling contains three primary uncertainty sources: experimental uncertainty, database modeling uncertainty, and database query interpolation uncertainty. The final database and uncertainty model represent a significant improvement in the quality of the aerodynamic predictions for this regime of flight over the estimates previously used by the Ares Project. The maximum possible aerodynamic force pushing the vehicle towards the launch tower assembly in a dispersed case using this database saw a 40 percent reduction from the worst-case scenario in previously released data for Ares I.

Hanke, Jeremy L.↗

A Hybrid CFD/Engineering Model for Predicting Plume Induced Erosion and Cratering

The Descent Interpolated Gas-Granular Erosion Model (DIGGEM) is a hybrid computational fluid dynamics (CFD) / engineering model calibrated with flight data from Apollo. It is used for the prediction of plume induced viscous erosion and cratering during the descent of landing vehicles. The model assumes a functional relationship between local shear stress and eroded mass flux. It is used along with a suite of other tools to predict complex plume induced environments for Human Landing System (HLS) and Commercial Lunar Payload Services (CLPSS) vehicles.

PSI↗

A Final Approach Trajectory Model for Current Operations

Predicting accurate trajectories with limited intent information is a challenge faced by air traffic management decision support tools in operation today. One such tool is the FAA's Terminal Proximity Alert system which is intended to assist controllers in maintaining safe separation of arrival aircraft during final approach. In an effort to improve the performance of such tools, two final approach trajectory models are proposed; one based on polynomial interpolation, the other on the Fourier transform. These models were tested against actual traffic data and used to study effects of the key final approach trajectory modeling parameters of wind, aircraft type, and weight class, on trajectory prediction accuracy. Using only the limited intent data available to today's ATM system, both the polynomial interpolation and Fourier transform models showed improved trajectory prediction accuracy over a baseline dead reckoning model. Analysis of actual arrival traffic showed that this improved trajectory prediction accuracy leads to improved inter-arrival separation prediction accuracy for longer look ahead times. The difference in mean inter-arrival separation prediction error between the Fourier transform and dead reckoning models was 0.2 nmi for a look ahead time of 120 sec, a 33 percent improvement, with a corresponding 32 percent improvement in standard deviation.

Gong, Chester↗

Interpolation schemes for orbit determination with the global positioning system

This paper demonstrates that the Global Positioning System (GPS) and simple interpolation schemes can be utilized to satisfy typical orbit determination demands of applications satellites. The complete GPS consists of 24 satellites and permits a position fix of a user satellite at any arbitrary instant. The interpolation formulae used in this report fit generalized Keplerian orbits through a number of position fixes. For a given accuracy level the telemetry requirements as measured by the time intervals between position fixes and the computational load of an interpolation as measured by the number of points through which a generalized Keplerian orbit is fitted vary inversely. A set of possible compromises between these two factors is presented. The Phase I GPS consists of 6 satellites and permits a position fix of a user satellite just 36% of the time. It is shown that with this system more sophisticated interpolation schemes which model atmospheric drag and higher degree terms of the earth's gravity field must be employed to obtain accurate orbits.

Argentiero, P.↗

Mixed models and reduction method for dynamic analysis of anisotropic shells

A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.

Noor, A. K.↗