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 343 records · Page 19

Measuring Positions of Objects using Two or More Cameras

An improved method of computing positions of objects from digitized images acquired by two or more cameras (see figure) has been developed for use in tracking debris shed by a spacecraft during and shortly after launch. The method is also readily adaptable to such applications as (1) tracking moving and possibly interacting objects in other settings in order to determine causes of accidents and (2) measuring positions of stationary objects, as in surveying. Images acquired by cameras fixed to the ground and/or cameras mounted on tracking telescopes can be used in this method. In this method, processing of image data starts with creation of detailed computer- aided design (CAD) models of the objects to be tracked. By rotating, translating, resizing, and overlaying the models with digitized camera images, parameters that characterize the position and orientation of the camera can be determined. The final position error depends on how well the centroids of the objects in the images are measured; how accurately the centroids are interpolated for synchronization of cameras; and how effectively matches are made to determine rotation, scaling, and translation parameters. The method involves use of the perspective camera model (also denoted the point camera model), which is one of several mathematical models developed over the years to represent the relationships between external coordinates of objects and the coordinates of the objects as they appear on the image plane in a camera. The method also involves extensive use of the affine camera model, in which the distance from the camera to an object (or to a small feature on an object) is assumed to be much greater than the size of the object (or feature), resulting in a truly two-dimensional image. The affine camera model does not require advance knowledge of the positions and orientations of the cameras. This is because ultimately, positions and orientations of the cameras and of all objects are computed in a coordinate system attached to one object as defined in its CAD model.

Klinko, Steve↗

A simplified scheme for computing radiation transfer in the troposphere

A scheme is presented, for the heating of clear and cloudy air by solar and infrared radiation transfer, designed for use in tropospheric general circulation models with coarse vertical resolution. A bulk transmission function is defined for the infrared transfer. The interpolation factors, required for computing the bulk transmission function, are parameterized as functions of such physical parameters as the thickness of the layer, the pressure, and the mixing ratio at a reference level. The computation procedure for solar radiation is significantly simplified by the introduction of two basic concepts. The first is that the solar radiation spectrum can be divided into a scattered part, for which Rayleigh scattering is significant but absorption by water vapor is negligible, and an absorbed part for which absorption by water vapor is significant but Rayleigh scattering is negligible. The second concept is that of an equivalent cloud water vapor amount which absorbs the same amount of radiation as the cloud.

Katayama, A.↗

Comparison of techniques for approximating ocean bottom topography in a wave-refraction computer model

A study of the effects of using different methods for approximating bottom topography in a wave-refraction computer model was conducted. Approximation techniques involving quadratic least squares, cubic least squares, and constrained bicubic polynomial interpolation were compared for computed wave patterns and parameters in the region of Saco Bay, Maine. Although substantial local differences can be attributed to use of the different approximation techniques, results indicated that overall computed wave patterns and parameter distributions were quite similar.

Poole, L. R.↗

Design and testing of the navigation model for three axis stabilized earth oriented satellites applied to the ATS-6 satellite image data base

An earth edge methodology has been developed to account for the relative attitude changes between successive ATS-6 images which allows reasonable high quality wind sets to be produced. The method consists of measuring the displacements of the right and left infrared earth edges between successive ATS-6 images as a function of scan line; from these measurements the attitude changes can be deduced and used to correct the apparent cloud displacement measurements. The wind data sets generated from ATS-6 using the earth-edge methodology were compared with those derived from the SMS-1 images (and model) covering the same time period. Quantitative comparisons for low level trade cumuli were made at interpolated uniformly spaced grid points and for selected individual comparison clouds. Selected individual comparison clouds, the root-mean-square differences for the U and V components were 1.0 and 1.2 meters per second with a maximum wind direction difference of 15 deg.

Kuhlow, W. W.↗

Analysis methods and preliminary design study

The shuttle orbiter body flap was studied to define loads and temperatures and to assess various structural concepts. Laminated-structure analysis capability was required. An appraisal of the available structural analysis programs, NASTRAN and SPAR, indicates that neither had all the capabilities required. The thermal analysis program, MITAS solved the problem once a model was generated but model generation and verification was laborious and transfer of temperatures to the structural program usually required interpolation. Therefore it was decided to incorporate both capabilities in the SPAR finite element program.

Blackburn, C. L.↗

Forecast experiments with the NASA/GLA stratospheric/tropospheric data assimilation system

For the first time, a four-dimensional stratospheric/tropospheric data assimilation system with a top analysis level at 0.4 mb has been developed and used to produce physically consistent gridded analyses for the stratosphere as well as the troposphere for a period during the First GARP Global Experiment (FGGE) and Limb Infrared Monitor of the Stratosphere (LIMS) (November 1978-May 1979). The system consists of a two-dimensional optimum interpolation analysis with 18 mandatory pressure levels and a 19-level fourth order stratospheric/tropospheric general circulation model with a horizontal resolution of 4 (latitude) by 5 deg (longitude) and a top at 0.3 mb. The system allows the utilization of stratospheric data including LIMS, Tiros-N retrievals, rocketsondes and vertical temperature profile radiometer soundings in addition to the other FGGE level 2b data. These data are analyzed every six hours. In order to examine the quality of the analyzed data, forecast experiments starting from different analyses are performed for the period of the stratospheric sudden warming of late February 1979. The results indicate that by employing the present four-dimensional assimilation approach, the medium-range forecast skill for this event is improved.

Takano, Kenji↗

Multi-band reflector antenna with double-ring element frequency selective subreflector

Frequency selective subreflectors (FSS) are often employed in the reflector antenna system of a communication satellite or a deep space exploration vehicle for multi-frequency operations. In the past, FSS's have been designed for diplexing two frequency bands. For example, the Voyager FSS was designed to diplex S and X bands and the TDRSS FSS was designed to diplex S and Ku bands. Recently, NASA's CASSINI project requires an FSS to multiplex four frequency (S/X/Ku/Ka) bands. Theoretical analysis and experimental verifications are presented for a multi-band flat pannel FSS with double-ring elements. Both the exact formulation and the thin-ring approximation are described for analyzing and designing this multi-ring patch element FSS. It is found that the thin-ring approximation fails to predict the electrically wide ring element FSS's performance. A single screen double-ring element FSS is demonstrated for the tri-band system that reflects the X-band signal while transmitting through the S- and Ku-band signals. In addition, a double screen FSS with non-similar double-ring elements is presented for the Cassini's four-band system which reflects the X- and Ka-band signals while passing the S- and Ku-band signals. To accurately predict the FSS effects on a dual reflector antenna's radiation pattern, the FSS subreflector's transmitted/reflected field variation as functions of the polarization and incident angles with respect to the local coordinates was taken into account. An FSS transmission/reflection coefficient table is computed for TE and TM polarizations at various incident angles based on the planar FSS model. Next, the hybrid Geometric Optics (GO) and Physical Optics (PO) technique is implemented with linearly interpolating the FSS table to efficiently determine the FSS effects in a dual reflector antenna.

Wu, Te-Kao↗

A Kernel-Free Particle-Finite Element Method for Hypervelocity Impact Simulation

An improved hybrid particle-finite element method has been developed for the simulation of hypervelocity impact problems. Unlike alternative methods, the revised formulation computes the density without reference to any kernel or interpolation functions, for either the density or the rate of dilatation. This simplifies the state space model and leads to a significant reduction in computational cost. The improved method introduces internal energy variables as generalized coordinates in a new formulation of the thermomechanical Lagrange equations. Example problems show good agreement with exact solutions in one dimension and good agreement with experimental data in a three dimensional simulation.

Park, Young-Keun↗

Hybrid Data Assimilation without Ensemble Filtering

The Global Modeling and Assimilation Office is preparing to upgrade its three-dimensional variational system to a hybrid approach in which the ensemble is generated using a square-root ensemble Kalman filter (EnKF) and the variational problem is solved using the Grid-point Statistical Interpolation system. As in most EnKF applications, we found it necessary to employ a combination of multiplicative and additive inflations, to compensate for sampling and modeling errors, respectively and, to maintain the small-member ensemble solution close to the variational solution; we also found it necessary to re-center the members of the ensemble about the variational analysis. During tuning of the filter we have found re-centering and additive inflation to play a considerably larger role than expected, particularly in a dual-resolution context when the variational analysis is ran at larger resolution than the ensemble. This led us to consider a hybrid strategy in which the members of the ensemble are generated by simply converting the variational analysis to the resolution of the ensemble and applying additive inflation, thus bypassing the EnKF. Comparisons of this, so-called, filter-free hybrid procedure with an EnKF-based hybrid procedure and a control non-hybrid, traditional, scheme show both hybrid strategies to provide equally significant improvement over the control; more interestingly, the filter-free procedure was found to give qualitatively similar results to the EnKF-based procedure.

Kalman Filter↗

Impacts of aerosol size-dependent below-cloud scavenging on tropospheric aerosol in the NASA GEOS model

Cloud scavenging is the major sink for a suite of tropospheric aerosols, and thus largely affects the global distribution and lifetime of aerosols as well as the formation of clouds. Simplified model parameterizations of cloud scavenging of aerosols substantially contribute to large uncertainties in the simulated aerosol burdens, spatial distributions, lifetimes, as well as aerosol radiative forcing. Here we evaluate the wet scavenging scheme in the current NASA GEOS model, develop more physically-based parameterizations of below-cloud scavenging, and assess the impacts of the new scavenging parameterizations on the GEOS simulations of aerosols. The current below-cloud scavenging (BCS; washout) scheme uses a first-order removal with a constant scavenging coefficient. We replaced the current scheme with a size-dependent parameterization of Croft et al.(Atmos. Chem. Phys., 2009; hereafter referred to as C09). The C09 BCS parameterizations calculate aerosol below-cloud scavenging coefficients based on a size-dependent collision efficiency between aerosol particles and raindrops (or snow crystals). A look-up table of compiled collision efficiencies provides values for sixty aerosol sizes and nine rain flux rates. Snow-aerosol collision efficiencies are a function of the sixty aerosol sizes only. We assume a log-normal size distribution for aerosols (bulk) in GOCART. The collision efficiencies are then calculated online using a bilinear interpolation from the look-up table for the given aerosol sizes and precipitation flux rates in the model.

aerosols↗

Toward Improving Short-Term Predictions of Fine Particulate Matter Over the United States Via Assimilation of Satellite Aerosol Optical Depth Retrievals

This study develops a new approach to improve simulations of the particulate matter of aerodynamic diameter smaller than 2.5μm (PM2.5) in the Community Multiscale Air Quality (CMAQ) model via assimilation of Moderate Resolution Imaging Spectroradiometer (MODIS) aerosol optical depth(AOD) retrievals using the Grid point Statistical Interpolation (GSI) system. In contrast to previous studies that only consider errors due to transport, our computation of the background error covariance matrix incorporates uncertainties in anthropogenic emissions. To understand the impact of this approach, three experiments (one background and two assimilations) are performed over the contiguous United States(CONUS) from 15 July to 14 August 2014. The background CMAQ experiment significantly underestimates both the MODIS AOD and surface PM2.5levels. MODIS AOD assimilation pushes both the CMAQ AOD and surface PM2.5distributions toward the observed distributions, but CMAQ still underestimates the observations. Averaged over CONUS, the two assimilation experiments with and without including the anthropogenic emission uncertainties improve the correlation coefficient between the model and independent observations of PM2.5 by ~67% and ~48%, respectively, and reduces the mean bias by ~38% and~10%, respectively. The assimilation improves the model performance everywhere over CONUS, except the New York and Wisconsin, where CMAQ overestimates the observed PM2.5during nighttime after assimilation likely because of overcorrection of aerosol mass concentrations by the AOD assimilation. Future work should incorporate uncertainties in other processes (biomass burning and biogenic emissions, deposition, chemistry, transport, and boundary conditions) to further enhance the value of assimilating spaceborne AOD retrievals.

Assimilation↗

An operator calculus for surface and volume modeling

The mathematical techniques which form the foundation for most of the surface and volume modeling techniques used in practice are briefly described. An outline of what may be termed an operator calculus for the approximation and interpolation of functions of more than one independent variable is presented. By considering the linear operators associated with bivariate and multivariate interpolation/approximation schemes, it is shown how they can be compounded by operator multiplication and Boolean addition to obtain a distributive lattice of approximation operators. It is then demonstrated via specific examples how this operator calculus leads to practical techniques for sculptured surface and volume modeling.

Gordon, W. J.↗

Recommendations for Using Noise Monitors to Estimate Noise Exposure During X-59 Community Tests

A low fidelity simulation approach is used to explore how to place and use noise monitors during X-59 QueSST community tests, where people’s annoyance to the noise produced by the X-59 aircraft will be gathered. Several recommendations are provided including: 1) the desired number of sparsely spaced noise monitor sites within the survey area, 2) whether to group and average measurements across multiple noise monitors located at a site, 3) what spacing should be used if grouped noise monitors are used, 4) an approach to mitigate ambient noise contamination at the measurement sites, 5) a method to combine empirical and predicted dose estimates to provide a single dose estimate for respondents, and 6) assessing how changes in turbulence intensity and array configuration affect dose uncertainty. To make these recommendations, the error that is expected when fitting contrived, smoothly varying sonic boom “reference exposure surfaces” is studied when a spatially sparse and scattered set of samples is used as responses for the fit. The reference exposure surfaces mimic the sonic boom exposure at ground level that might be expected in the X-59 survey area in the absence of atmospheric turbulence, ambient noise, and other localized effects. The spatial extent of these surfaces varies and is representative of the different survey area sizes that might be expected during future X-59 community overflight tests. These contrived reference surfaces are sampled, and those reference samples are then perturbed to mimic atmospheric turbulence, ambient noise and other localized effects that might affect noise monitor measurements within overflown communities. Two different surface fitting methods are investigated when fitting these perturbed samples to approximate the reference surface. The first method uses interpolation between the perturbed data at the scattered sites to compute the fit. The second method fits a polynomial surface model to the perturbed data using ordinary least squares regression analysis. For both fitting methods, the root mean square fit error is computed from the pointwise difference between the fit surface and the reference surface as the count and configuration of the sites is varied while also averaging the error across many different realizations of both the smooth variation of the reference exposure surface and the random, localized perturbations at the sample sites. Different site configurations are compared using this error statistic to make the recommendations noted above. Additionally, the two fitting approaches (interpolation vs linear regression) are compared based on the fit error observed in these simulations. These analyses, comparisons, and recommendations should inform future decisions on the noise monitor placement and the methods used to analyze the noise monitor data that is collected during X-59 community overflights.

sonic boom↗

A Flexible Method for Producing F.E.M. Analysis of Bone Using Open-Source Software

This project, performed in support of the NASA GRC Space Academy summer program, sought to develop an open-source workflow methodology that segmented medical image data, created a 3D model from the segmented data, and prepared the model for finite-element analysis. In an initial step, a technological survey evaluated the performance of various existing open-source software that claim to perform these tasks. However, the survey concluded that no single software exhibited the wide array of functionality required for the potential NASA application in the area of bone, muscle and bio fluidic studies. As a result, development of a series of Python scripts provided the bridging mechanism to address the shortcomings of the available open source tools. The implementation of the VTK library provided the most quick and effective means of segmenting regions of interest from the medical images; it allowed for the export of a 3D model by using the marching cubes algorithm to build a surface mesh. To facilitate the development of the model domain from this extracted information required a surface mesh to be processed in the open-source software packages Blender and Gmsh. The Preview program of the FEBio suite proved to be sufficient for volume filling the model with an unstructured mesh and preparing boundaries specifications for finite element analysis. To fully allow FEM modeling, an in house developed Python script allowed assignment of material properties on an element by element basis by performing a weighted interpolation of voxel intensity of the parent medical image correlated to published information of image intensity to material properties, such as ash density. A graphical user interface combined the Python scripts and other software into a user friendly interface. The work using Python scripts provides a potential alternative to expensive commercial software and inadequate, limited open-source freeware programs for the creation of 3D computational models. More work will be needed to validate this approach in creating finite-element models.

gravitational physiology↗

Impact of Lidar Wind Sounding on Mesoscale Forecast

An Observing System Simulation Experiment (OSSE) was conducted to study the impact of airborne lidar wind sounding on mesoscale weather forecast. A wind retrieval scheme, which interpolates wind data from a grid data system, simulates the retrieval of wind profile from a satellite lidar system. A mesoscale forecast system based on the PSU/NCAR MM5 model is developed and incorporated the assimilation of the retrieved line-of-sight wind. To avoid the "identical twin" problem, the NCEP reanalysis data is used as our reference "nature" atmosphere. The simulated space-based lidar wind observations were retrieved by interpolating the NCEP values to the observation locations. A modified dataset obtained by smoothing the NCEP dataset was used as the initial state whose forecast was sought to be improved by assimilating the retrieved lidar observations. Forecasts using wind profiles with various lidar instrument parameters has been conducted. The results show that to significantly improve the mesoscale forecast the satellite should fly near the storm center with large scanning radius. Increasing lidar firing rate also improves the forecast. Cloud cover and lack of aerosol degrade the quality of the lidar wind data and, subsequently, the forecast.

Miller, Timothy L.↗

A mesoscale sixth-order numerical modelling system

A numerical simulation system is currently under development for NASA which is intended to improve the modeling of subsynoptic and mesoscale adjustments associated with cyclogenesis, severe storm development and atmospheric transport processes. The model utilizes a standard hydrostatic sigma-p coordinate primitive equation set, with x,y-space differencing accurate to eighth order. A three-step dynamic initialization procedure is employed between the analysis of real-time data and grid interpolation. Results of an 18-hour simulation during which synoptic scale cyclogenesis, subsynoptic scale jet streak adjustments, mesoscale convergence zones and tornadic storms were observed have shown the present model to have the potential for simulating the fine-scale structure of features associated with cyclogenesis and intense squall-line development. The mesoscale model was also found to produce less truncation than the NWS LFM model, although a frictionless version of the mesoscale model somewhat overdeepens and overaccelerates features.

Kaplan, M. L.↗

Sensitivity of Age-of-Air Calculations to the Choice of Advection Scheme

The age of air has recently emerged as a diagnostic of atmospheric transport unaffected by chemical parameterizations, and the features in the age distributions computed in models have been interpreted in terms of the models' large-scale circulation field. This study shows, however, that in addition to the simulated large-scale circulation, three-dimensional age calculations can also be affected by the choice of advection scheme employed in solving the tracer continuity equation, Specifically, using the 3.0deg latitude X 3.6deg longitude and 40 vertical level version of the Geophysical Fluid Dynamics Laboratory SKYHI GCM and six online transport schemes ranging from Eulerian through semi-Lagrangian to fully Lagrangian, it will be demonstrated that the oldest ages are obtained using the nondiffusive centered-difference schemes while the youngest ages are computed with a semi-Lagrangian transport (SLT) scheme. The centered- difference schemes are capable of producing ages older than 10 years in the mesosphere, thus eliminating the "young bias" found in previous age-of-air calculations. At this stage, only limited intuitive explanations can be advanced for this sensitivity of age-of-air calculations to the choice of advection scheme, In particular, age distributions computed online with the National Center for Atmospheric Research Community Climate Model (MACCM3) using different varieties of the SLT scheme are substantially older than the SKYHI SLT distribution. The different varieties, including a noninterpolating-in-the-vertical version (which is essentially centered-difference in the vertical), also produce a narrower range of age distributions than the suite of advection schemes employed in the SKYHI model. While additional MACCM3 experiments with a wider range of schemes would be necessary to provide more definitive insights, the older and less variable MACCM3 age distributions can plausibly be interpreted as being due to the semi-implicit semi-Lagrangian dynamics employed in the MACCM3. This type of dynamical core (employed with a 60-min time step) is likely to reduce SLT's interpolation errors that are compounded by the short-term variability characteristic of the explicit centered-difference dynamics employed in the SKYHI model (time step of 3 min). In the extreme case of a very slowly varying circulation, the choice of advection scheme has no effect on two-dimensional (latitude-height) age-of-air calculations, owing to the smooth nature of the transport circulation in 2D models. These results suggest that nondiffusive schemes may be the preferred choice for multiyear simulations of tracers not overly sensitive to the requirement of monotonicity (this category includes many greenhouse gases). At the same time, age-of-air calculations offer a simple quantitative diagnostic of a scheme's long-term diffusive properties and may help in the evaluation of dynamical cores in multiyear integrations. On the other hand, the sensitivity of the computed ages to the model numerics calls for caution in using age of air as a diagnostic of a GCM's large-scale circulation field.

Eluszkiewicz, Janusz↗

Numerical prediction of the Mid-Atlantic states cyclone of 18-19 February 1979

A series of forecast experiments was conducted to assess the accuracy of the GLAS model, and to determine the importance of large scale dynamical processes and diabatic heating to the cyclogenesis. The GLAS model correctly predicted intense coastal cyclogenesis and heavy precipitation. Repeated without surface heat and moisture fluxes, the model failed to predict any cyclone development. An extended range forecast, a forecast from the NMC analysis interpolated to the GLAS grid, and a forecast from the GLAS analysis with the surface moisture flux excluded predicted weak coastal low development. Diabatic heating resulting from oceanic fluxes significantly contributed to the generation of low level cyclonic vorticity and the intensification and slow rate of movement of an upper level ridge over the western Atlantic. As an upper level short wave trough approached this ridge, diabatic heating associated with the release of latent heat intensified, and the gradient of vorticity, vorticity advection and upper level divergence in advance of the trough were greatly increased, providing strong large scale forcing for the surface cyclogenesis.

Atlas, R.↗