Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data intensive computing”

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 253 records · Page 14

Introduction to Analysis Methods for Big Earth Data

Big Earth Data are too big to be tractable to simple data inspection. Thus, they typically require models to make sense of all the data. Useful models for Big Earth Data may be physical, statistical, or machine learning based. While physical models are ideal for understanding the data, they are not always feasible, particularly when our ability to observe at finer scales exceeds our ability to incorporate the physics. Statistical models are more generalized, but computationally intensive for many Earth Observation datasets. Machine Learning models generally scale well but are sometimes limited in the physical understanding they can offer. Hybrid models combine attributes—and advantages—of two or more of these types.

Christopher Lynnes↗

Automatic Incorporation of Flaws into Cad and FEA Models Based On Measured NDE Data

Automatic incorporation of non-destructive evaluation (NDE) data into CAD and finite element (FE) models for structural analysis is essential for accurate residual-life and strength analyses of composite structures. Currently, two significant but independent challenges exist for incorporating NDE data in such analyses: interpretation of the NDE signals to locate and characterize the structural defects and the incorporation of that information into structural analysis tools. The first challenge will be addressed through linear inversion of calculated synthetic responses generated by a simulation engine to evaluate flaw size, location, and other characteristics. The second challenge will be addressed by building on prior and parallel defect modeling efforts to integrate suitable flaw models into the CAD/finite-element representation of a specimen. Ultrasonic NDE information is usually represented using C-Scans of time-of-flight or intensity information. Other modalities such as thermography and computed tomography (CT) could also be used to characterize the structure. This information will first be projected into material coordinates and then aligned with the specimen’s finite element mesh. The NDE data can then be incorporated by inserting flaw geometry into the CAD and FE analysis model and automatic modification of local material properties and boundary conditions. This process will enable the use of structural analysis tools even for practical modalities such as ultrasound and thermography which may not directly give accurate-enough geometric information to fully define the flaw. The full flaw characteristics will be inferred using prior expert and modeling knowledge to categorize flaw signatures. Outputs for this project will include demonstrated methods for creating structural analysis geometric models that incorporate flaw characterization information; these models can then be used directly in popular commercial FEA software. In addition, the primary deliverable for the project will include a software that will be able to automatically incorporate different kinds of flaws into an existing FE models. The models will be generated in standard file formats to enable direct import into the analysis software.

Adarsh Krishnamurthy↗

The QuakeSim Project: Numerical Simulations for Active Tectonic Processes

In order to develop a solid earth science framework for understanding and studying of active tectonic and earthquake processes, this task develops simulation and analysis tools to study the physics of earthquakes using state-of-the art modeling, data manipulation, and pattern recognition technologies. We develop clearly defined accessible data formats and code protocols as inputs to the simulations. these are adapted to high-performance computers because the solid earth system is extremely complex and nonlinear resulting in computationally intensive problems with millions of unknowns. With these tools it will be possible to construct the more complex models and simulations necessary to develop hazard assessment systems critical for reducing future losses from major earthquakes.

tectonic processes↗

Pioneer Venus polarimetry and haze optical thickness

The Pioneer Venus mission provided us with high-resolution measurements at four wavelengths of the linear polarization of sunlight reflected by the Venus atmosphere. These measurements span the complete phase angle range and cover a period of more than a decade. A first analysis of these data by Kawabata et al. confirmed earlier suggestions of a haze layer above and partially mixed with the cloud layer. They found that the haze exhibits large spatial and temporal variations. The haze optical thickness at a wavelength of 365 nm was about 0.06 at low latitudes, but approximately 0.8 at latitudes from 55 deg poleward. Differences between morning and evening terminator have also been reported by the same authors. Using an existing cloud/haze model of Venus, we study the relationship between the haze optical thickness and the degree of linear polarization. Variations over the visible disk and phase angle dependence are investigated. For that purpose, exact multiple scattering computations are compared with Pioneer Venus measurements. To get an impression of the variations over the visible disk, we have first studied scans of the polarization parallel to the intensity equator. After investigating a small subset of the available data we have the following results. Adopting the haze particle characteristics given by Kawabata et al., we find a thickening of the haze at increasing latitudes. Further, we see a difference in haze optical thickness between the northern and southern hemispheres that is of the same order of magnitude as the longitudinal variation of haze thickness along a scan line. These effects are most pronounced at a wavelength of 935 nm. We must emphasize the tentative nature of the results, because there is still an enormous amount of data to be analyzed. We intend to combine further polarimetric research of Venus with constraints on the haze parameters imposed by physical and chemical processes in the atmosphere.

Knibbe, W. J. J.↗

Remote sensing data applied to the evaluation of soil erosion caused by land-use. Ribeirao Anhumas Basin Area: A case study

Formerly covered by a tropical forest, the study area was deforested in the early 40's for coffee plantation and cattle raising, which caused intense gully erosion problems. To develop a method to analyze the relationship between land use and soil erosion, visual interpretations of aerial photographs (scale 1:25.000), MSS-LANDSAT imagery (scale 1:250,000), as well as automatic interpretation of computer compatible tapes by IMAGE-100 system were carried out. From visual interpretation the following data were obtained: land use and cover tapes, slope classes, ravine frequency, and a texture sketch map. During field work, soil samples were collected for texture and X-ray analysis. The texture sketch map indicate that the areas with higher slope angles have a higher susceptibilty to the development of gullies. Also, the over carriage of pastureland, together with very friable lithologies (mainly sandstone) occuring in that area, seem to be the main factors influencing the catastrophic extension of ravines in the study site.

Dejesusparada, N.↗

Phase Time and Envelope Time in Time-Distance Analysis and Acoustic Imaging

Time-distance analysis and acoustic imaging are two related techniques to probe the local properties of solar interior. In this study, we discuss the relation of phase time and envelope time between the two techniques. The location of the envelope peak of the cross correlation function in time-distance analysis is identified as the travel time of the wave packet formed by modes with the same w/l. The phase time of the cross correlation function provides information of the phase change accumulated along the wave path, including the phase change at the boundaries of the mode cavity. The acoustic signals constructed with the technique of acoustic imaging contain both phase and intensity information. The phase of constructed signals can be studied by computing the cross correlation function between time series constructed with ingoing and outgoing waves. In this study, we use the data taken with the Taiwan Oscillation Network (TON) instrument and the Michelson Doppler Imager (MDI) instrument. The analysis is carried out for the quiet Sun. We use the relation of envelope time versus distance measured in time-distance analyses to construct the acoustic signals in acoustic imaging analyses. The phase time of the cross correlation function of constructed ingoing and outgoing time series is twice the difference between the phase time and envelope time in time-distance analyses as predicted. The envelope peak of the cross correlation function between constructed ingoing and outgoing time series is located at zero time as predicted for results of one-bounce at 3 mHz for all four data sets and two-bounce at 3 mHz for two TON data sets. But it is different from zero for other cases. The cause of the deviation of the envelope peak from zero is not known.

Chou, Dean-Yi↗

Photoionization and Recombination

Theoretically self-consistent calculations for photoionization and (e + ion) recombination are described. The same eigenfunction expansion for the ion is employed in coupled channel calculations for both processes, thus ensuring consistency between cross sections and rates. The theoretical treatment of (e + ion) recombination subsumes both the non-resonant recombination ("radiative recombination"), and the resonant recombination ("di-electronic recombination") processes in a unified scheme. In addition to the total, unified recombination rates, level-specific recombination rates and photoionization cross sections are obtained for a large number of atomic levels. Both relativistic Breit-Pauli, and non-relativistic LS coupling, calculations are carried out in the close coupling approximation using the R-matrix method. Although the calculations are computationally intensive, they yield nearly all photoionization and recombination parameters needed for astrophysical photoionization models with higher precision than hitherto possible, estimated at about 10-20% from comparison with experimentally available data (including experimentally derived DR rates). Results are electronically available for over 40 atoms and ions. Photoionization and recombination of He-, and Li-like C and Fe are described for X-ray modeling. The unified method yields total and complete (e+ion) recombination rate coefficients, that can not otherwise be obtained theoretically or experimentally.

Nahar, Sultana N.↗

Scientific Programming Using Java: A Remote Sensing Example

This paper presents results of a project to port remote sensing code from the C programming language to Java. The advantages and disadvantages of using Java versus C as a scientific programming language in remote sensing applications are discussed. Remote sensing applications deal with voluminous data that require effective memory management, such as buffering operations, when processed. Some of these applications also implement complex computational algorithms, such as Fast Fourier Transformation analysis, that are very performance intensive. Factors considered include performance, precision, complexity, rapidity of development, ease of code reuse, ease of maintenance, memory management, and platform independence. Performance of radiometric calibration code written in Java for the graphical user interface and of using C for the domain model are also presented.

Prados, Don↗

First-Order Model Management With Variable-Fidelity Physics Applied to Multi-Element Airfoil Optimization

First-order approximation and model management is a methodology for a systematic use of variable-fidelity models or approximations in optimization. The intent of model management is to attain convergence to high-fidelity solutions with minimal expense in high-fidelity computations. The savings in terms of computationally intensive evaluations depends on the ability of the available lower-fidelity model or a suite of models to predict the improvement trends for the high-fidelity problem, Variable-fidelity models can be represented by data-fitting approximations, variable-resolution models. variable-convergence models. or variable physical fidelity models. The present work considers the use of variable-fidelity physics models. We demonstrate the performance of model management on an aerodynamic optimization of a multi-element airfoil designed to operate in the transonic regime. Reynolds-averaged Navier-Stokes equations represent the high-fidelity model, while the Euler equations represent the low-fidelity model. An unstructured mesh-based analysis code FUN2D evaluates functions and sensitivity derivatives for both models. Model management for the present demonstration problem yields fivefold savings in terms of high-fidelity evaluations compared to optimization done with high-fidelity computations alone.

Alexandrov, N. M.↗

RECOVER: An Automated Cloud-Based Decision Support System for Post-fire Rehabilitation Planning

RECOVER is a site-specific decision support system that automatically brings together in a single analysis environment the information necessary for post-fire rehabilitation decision-making. After a major wildfire, law requires that the federal land management agencies certify a comprehensive plan for public safety, burned area stabilization, resource protection, and site recovery. These burned area emergency response (BAER) plans are a crucial part of our national response to wildfire disasters and depend heavily on data acquired from a variety of sources. Final plans are due within 21 days of control of a major wildfire and become the guiding document for managing the activities and budgets for all subsequent remediation efforts. There are few instances in the federal government where plans of such wide-ranging scope and importance are assembled on such short notice and translated into action more quickly. RECOVER has been designed in close collaboration with our agency partners and directly addresses their high-priority decision-making requirements. In response to a fire detection event, RECOVER uses the rapid resource allocation capabilities of cloud computing to automatically collect Earth observational data, derived decision products, and historic biophysical data so that when the fire is contained, BAER teams will have a complete and ready-to-use RECOVER dataset and GIS analysis environment customized for the target wildfire. Initial studies suggest that RECOVER can transform this information-intensive process by reducing from days to a matter of minutes the time required to assemble and deliver crucial wildfire-related data.

cloud computing↗

Second moment closure analysis of the backstep flow database

A Second Moment Closure computation (SMC) is compared in detail with the Direct Numerical Simulation (DNS) data of Le and Moin for the backstep flow at Re = 5,000 in an attempt to understand why the intensity of the backflow and, consequently, the friction coefficient in the recirculation bubble are severely underestimated. The data show that this recirculation bubble is far from being laminar except in the very near wall layer. A novel 'differential a priori' procedure was used, in which the full transport equation for one isolated component of the Reynolds stress tensor was solved using DNS data as input. Conclusions are then different from what would have been deduced by comparing a full simulation to a DNS. One cause of discrepancy was traced back to insufficient transfer of energy to the normal stress by pressure strain, but was not cured. A significant finding, confirmed by the DNS data in the core region of a channel flow, is that the coefficient that controls destruction of dissipation, C epsilon(sub 2), should be decreased by a factor of 2 when production is vanishing. This is also the case in the recirculation bubble, and a new formulation has cured 25% of the backflow discrepancy.

Parneix, S.↗

Efficient algorithms for robust feature matching

One of the basic building blocks in any point-based registration scheme involves matching feature points that are extracted from the sensed image to their counterparts in the reference image. This leads to the fundamental problem of point matching: given two sets of points, find the affine transformation that transforms one point set so that its distance from the other point set is minimized. Because of measurement errors and the presence of outlying data points, it is important that the distance measure between two point sets be robust to these effects. We measure distances using the generalized Hausdorff distance. Point matching can be a computationally intensive task, and there have been a number of algorithms and approaches proposed for solving this problem both theoretical and applied. We present two approaches to the point matching problem, in an attempt to reduce the computational complexity of the problem, while still providing guarantees on the quality of the final match. Our first method is an approximation algorithm, which is loosely based on a branch-and-bound approach due to Huttenlocher and Rucklidge. We show that by varying the approximation error bounds, it is possible to achieve a tradeoff between the quality of the match and the running time of the algorithm. Our second method involves a Monte Carlo method for accelerating the search process used in the first algorithm. With high probability this method succeeds in finding an approximately optimal match. We establish the efficiency of our approaches empirically.

Mount, David M.↗

Satellite Orbital Drag During Magnetic Storms

We investigate satellite orbital drag effects at lowEarth orbit associated with thermosphere heating during magnetic storms caused by coronal mass ejections. CHAllenge Minisatellite Payload (CHAMP) and Gravity Recovery And Climate Experiment (GRACE) neutral density data are used to compute orbital drag. Stormtoquiet density comparisons are performed with background densities obtained by the JacchiaBowman 2008 (JB2008) empirical model. Our storms are grouped in different categories regarding their intensities as indicated by minimum values of the SYMH index. We then perform superposed epoch analyses with storm main phase onset as zero epoch time. In general, we find that orbital drag effects are larger for CHAMP (lower altitudes) in comparison to GRACE (higher altitudes). Results show that storm time drag effects manifest first at high latitudes, but for extreme storms, particularly observed by GRACE, stronger orbital drag effects occur during early main phase at low/equatorial latitudes, probably due to heating propagation from high latitudes. We find that storm time orbital decay along the satellites' path generally increases with storm intensity, being stronger and faster for the most extreme events. For these events, orbital drag effects decrease faster probably due to elevated cooling effects caused by nitric oxide, which introduce modeled density uncertainties during storm recovery phase. Errors associated with total orbit decay introduced by JB2008 are generally the largest for the strongest storms and increase during storm times, particular during recovery phases. We discuss the implication of these uncertainties for the prediction of collision between space objects at lowEarth orbit during magnetic storms.

Oliveira, D. M.↗

Comparison of ERBS orbit determination accuracy using batch least-squares and sequential methods

The Flight Dynamics Div. (FDD) at NASA-Goddard commissioned a study to develop the Real Time Orbit Determination/Enhanced (RTOD/E) system as a prototype system for sequential orbit determination of spacecraft on a DOS based personal computer (PC). An overview is presented of RTOD/E capabilities and the results are presented of a study to compare the orbit determination accuracy for a Tracking and Data Relay Satellite System (TDRSS) user spacecraft obtained using RTOS/E on a PC with the accuracy of an established batch least squares system, the Goddard Trajectory Determination System (GTDS), operating on a mainframe computer. RTOD/E was used to perform sequential orbit determination for the Earth Radiation Budget Satellite (ERBS), and the Goddard Trajectory Determination System (GTDS) was used to perform the batch least squares orbit determination. The estimated ERBS ephemerides were obtained for the Aug. 16 to 22, 1989, timeframe, during which intensive TDRSS tracking data for ERBS were available. Independent assessments were made to examine the consistencies of results obtained by the batch and sequential methods. Comparisons were made between the forward filtered RTOD/E orbit solutions and definitive GTDS orbit solutions for ERBS; the solution differences were less than 40 meters after the filter had reached steady state.

Oza, D. H.↗

SP mountain data analysis

An analysis of synthetic aperture radar data of SP Mountain was undertaken to demonstrate the use of digital image processing techniques to aid in geologic interpretation of SAR data. These data were collected with the ERIM X- and L-band airborne SAR using like- and cross-polarizations. The resulting signal films were used to produce computer compatible tapes, from which four-channel imagery was generated. Slant range-to-ground range and range-azimuth-scale corrections were made in order to facilitate image registration; intensity corrections were also made. Manual interpretation of the imagery showed that L-band represented the geology of the area better than X-band. Several differences between the various images were also noted. Further digital analysis of the corrected data was done for enhancement purposes. This analysis included application of an MSS differencing routine and development of a routine for removal of relief displacement. It was found that accurate registration of the SAR channels is critical to the effectiveness of the differencing routine. Use of the relief displacement algorithm on the SP Mountain data demonstrated the feasibility of the technique.

Rawson, R. F.↗

Predicting Strength and Fatigue for Suited and Unsuited Conditions from Empirical Data

The need for longer and more labor-intensive extra-vehicular activities (EVA) is required for construction and maintenance of the International Space Station (ISS). Issues pertaining to human performance while wearing a space suit (EMU) for prolonged periods have become more important. This project was conducted to investigate how a pressurized Extra-vehicular Mobility Unit (EMU) affects human upper body joint strength and fatigue and how to predict it from computer models based on the data collected.

Maida, James C.↗

Computation of the tip vortex flowfield for advanced aircraft propellers

The tip vortex flowfield plays a significant role in the performance of advanced aircraft propellers. The flowfield in the tip region is complex, three-dimensional and viscous with large secondary velocities. An analysis is presented using an approximate set of equations which contains the physics required by the tip vortex flowfield, but which does not require the resources of the full Navier-Stokes equations. A computer code was developed to predict the tip vortex flowfield of advanced aircraft propellers. A grid generation package was developed to allow specification of a variety of advanced aircraft propeller shapes. Calculations of the tip vortex generation on an SR3 type blade at high Reynolds numbers were made using this code and a parametric study was performed to show the effect of tip thickness on tip vortex intensity. In addition, calculations of the tip vortex generation on a NACA 0012 type blade were made, including the flowfield downstream of the blade trailing edge. Comparison of flowfield calculations with experimental data from an F4 blade was made. A user's manual was also prepared for the computer code (NASA CR-182178).

Tsai, Tommy M.↗

Precision processing of earth image data

Precise corrections of Landsat data are useful for generating land-use maps, detecting various crops and determining their acreage, and detecting changes. The paper discusses computer processing and visualization techniques for Landsat data so that users can get more information from the imagery. The elementary unit of data in each band of each scene is the integrated value of intensity of reflected light detected in the field of view by each sensor. To develop the basic mathematical approach for precision correction of the data, differences between positions of ground control points on the reference map and the observed control points in the scene are used to evaluate the coefficients of cubic time functions of roll, pitch, and yaw, and a linear time function of altitude deviation from normal height above local earth's surface. The resultant equation, termed a mapping function, corrects the warped data image into one that approximates the reference map. Applications are discussed relative to shade prints, extraction of road features, and atlas of cities.

Bernstein, R.↗