Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “algorithms and data structure”

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 1,351 records · Page 75

A heuristic finite-state model of the human driver in a car-following situation

An approach to modeling human driver behavior in single-lane car following which is based on a finite-state decision structure is considered. The specific strategy at each point in the decision tree was obtained from observations of typical driver behavior. The synthesis of the decision logic is based on position and velocity thresholds and four states defined by regions in the phase plane. The performance of the resulting assumed intuitively logical model was compared with actual freeway data. The match of the model to the data was optimized by adapting the model parameters using a modified PARTAN algorithm. The results indicate that the heuristic model behavior matches actual car-following performance better during deceleration and constant velocity phases than during acceleration periods.

Burnham, G. O.↗

Prediction of pilot reserve attention capacity during air-to-air target tracking

Reserve attention capacity of a pilot was calculated using a pilot model that allocates exclusive model attention according to the ranking of task urgency functions whose variables are tracking error and error rate. The modeled task consisted of tracking a maneuvering target aircraft both vertically and horizontally, and when possible, performing a diverting side task which was simulated by the precise positioning of an electrical stylus and modeled as a task of constant urgency in the attention allocation algorithm. The urgency of the single loop vertical task is simply the magnitude of the vertical tracking error, while the multiloop horizontal task requires a nonlinear urgency measure of error and error rate terms. Comparison of model results with flight simulation data verified the computed model statistics of tracking error of both axes, lateral and longitudinal stick amplitude and rate, and side task episodes. Full data for the simulation tracking statistics as well as the explicit equations and structure of the urgency function multiaxis pilot model are presented.

Onstott, E. D.↗

Reflection Coefficients on Surfaces of Different Periodic Structure

Diffraction properties of lossy periodic gratings with the metal base were investigated by solving Maxwell's equations numerically using the differential method. Two periodic surfaces were employed in the simulation: triangle structure and tilting triangle structure. Based on the numerical solution and in conjunction with the algorithm of Adams-Moulton, we computed reflection coefficients of plane waves with different wavelengths and different incident angles. The dielectric properties were also explored using various dielectric constants. The results show that the reflection coefficients of both TE and TM waves are quite sensitive to the incident angles of the plane waves when the metal sheet exists, which is in good agreement with the experimental data.

Niu, Pengfei↗

The Final Approach Spacing Tool

A system for assisting terminal area air traffic controllers in the management and control of arrival traffic, referred to as the Final Approach Spacing Tool (FAST), is being developed at NASA Ames Research Center. In a cooperative program, NASA and FAA have efforts underway to install and evaluate the system at the Dallas/Fort Worth Terminal Radar Approach Control facility. This paper will review the software architecture, the algorithms components, and the human-machine interface. The system is based on continuous updates of a detailed trajectory analyses of all arrival aircraft. FAST interprets the results of these trajectory analyses to build an efficient and procedurally acceptable plan for the arrival traffic that consists of a sequence, schedule, and runway assignment. The system utilizes a heuristically-based conflict resolution algorithm to build a solution trajectory that satisfies the plan, It extracts a series of speed and heading advisories from the solution trajectory to assist the controller in efficiently managing and controlling the arrival traffic down to the runway. The advisories are displayed in a graphical format to the controller. In addition to the radar tracking data, the system also relies on a series of data bases. These data bases contain aircraft performance models, airline preferred operational procedures, airspace structure, air traffic procedural models, and a three dimensional wind model. Field evaluation of FAST is expected to begin in 1994.

Davis, Thomas J.↗

A Parameterized Inversion Model for Soil Moisture and Biomass from Polarimetric Backscattering Coefficients

A semi-empirical algorithm for the retrieval of soil moisture, root mean square (RMS) height and biomass from polarimetric SAR data is explained and analyzed in this paper. The algorithm is a simplification of the distorted Born model. It takes into account the physical scattering phenomenon and has three major components: volume, double-bounce and surface. This simplified model uses the three backscattering coefficients ( sigma HH, sigma HV and sigma vv) at low-frequency (P-band). The inversion process uses the Levenberg-Marquardt non-linear least-squares method to estimate the structural parameters. The estimation process is entirely explained in this paper, from initialization of the unknowns to retrievals. A sensitivity analysis is also done where the initial values in the inversion process are varying randomly. The results show that the inversion process is not really sensitive to initial values and a major part of the retrievals has a root-mean-square error lower than 5% for soil moisture, 24 Mg/ha for biomass and 0.49 cm for roughness, considering a soil moisture of 40%, roughness equal to 3cm and biomass varying from 0 to 500 Mg/ha with a mean of 161 Mg/ha

soil moisture↗

First Image Products from EcoSAR - Osa Peninsula, Costa Rica

Designed especially for forest ecosystem studies, EcoSAR employs state-of-the-art digital beamforming technology to generate wide-swath, high-resolution imagery. EcoSARs dual antenna single-pass imaging capability eliminates temporal decorrelation from polarimetric and interferometric analysis, increasing the signal strength and simplifying models used to invert forest structure parameters. Antennae are physically separated by 25 meters providing single pass interferometry. In this mode the radar is most sensitive to topography. With 32 active transmit and receive channels, EcoSARs digital beamforming is an order of magnitude more versatile than the digital beamforming employed on the upcoming NISAR mission. EcoSARs long wavelength (P-band, 435 MHz, 69 cm) measurements can be used to simulate data products for ESAs future BIOMASS mission, allowing scientists to develop algorithms before the launch of the satellite. EcoSAR can also be deployed to collect much needed data where BIOMASS satellite wont be allowed to collect data (North America, Europe and Arctic), filling in the gaps to keep a watchful eye on the global carbon cycle. EcoSAR can play a vital role in monitoring, reporting and verification schemes of internationals programs such as UN-REDD (United Nations Reducing Emissions from Deforestation and Degradation) benefiting global society. EcoSAR was developed and flown with support from NASA Earth Sciences Technology Offices Instrument Incubator Program.

InSAR↗

Finite Element Based Model of Fiber-optic Strain Sensor Mounted on a Semi-monocoque Structure

Fiber-optic strain sensing (FOSS) provides an opportunity for additional awareness and can enable other mitigations such as flutter suppression or gust load alleviation but requires an accurate computational model as part of a design process that integrates structures and controls. A modeling approach for fiber-optic strain sensors has been developed using a shell element formulation in a detailed finite element model. The new FOSS modeling approach is designed to better scale to more complex realistic structures and fiber layouts than previous approaches. The modeling method was applied to the X-56A vehicle and results were compared to experimental static ground-test data. Additionally, the experimental sensor noise of these fiber optic sensors was assessed for the correlations in time and correlations along the fiber. The noise does show significant correlations both in time and along the fibers. These correlations of the noise should also be considered in modeling and assessing algorithms that will use FOSS.

Jeffrey Ouellette↗

Finite Element Based Model of Fiber-optic Strain Sensor Mounted on a Semi-monocoque Structure

Fiber-optic strain sensing (FOSS) provides an opportunity for additional awareness and can enable other mitigations such as flutter suppression or gust load alleviation but requires an accurate computational model as part of a design process that integrates structures and controls. A modeling approach for fiber-optic strain sensors has been developed using a shell element formulation in a detailed finite element model. The new FOSS modeling approach is designed to better scale to more complex realistic structures and fiber layouts than previous approaches. The modeling method was applied to the X-56A vehicle and results were compared to experimental static ground-test data. Additionally, the experimental sensor noise of these fiber optic sensors was assessed for the correlations in time and correlations along the fiber. The noise does show significant correlations both in time and along the fibers. These correlations of the noise should also be considered in modeling and assessing algorithms that will use FOSS.

Jeffrey Ouellette↗

Validation of the EDGES Low-band Antenna Beam Model

The response of the antenna is a source of uncertainty in measurements with the Experiment to Detect the Global Epoch of Reionization Signature (EDGES). We aim to validate the electromagnetic beam model of the low-band (50–100 MHz) dipole antenna with comparisons between models and against data. We find that simulations of a simplified model of the antenna over an infinite perfectly conducting ground plane are, with one exception, robust to changes in the numerical electromagnetic solver code or algorithm. For simulations of the antenna with the actual finite ground plane and realistic soil properties, we find that two out of three numerical solvers agree well. Applying our analysis pipeline to a simulated drift-scan observation from an early EDGES low-band instrument that had a 10 m × 10 m ground plane, we find residual levels after fitting and removing a five-term foreground model from the simulated data binned in local sidereal time (LST) average about 250 mK with ±40 mK variation between numerical solvers. A similar analysis of the primary 30 m × 30 m sawtooth ground plane reduced the LST-averaged residuals to about 90 mK with ±10 mK between the two viable solvers. More broadly we show that larger ground planes generally perform better than smaller ground planes. Simulated data have a power that is within 4% of real observations, a limitation of net accuracy of the sky and beam models. We observe that residual spectral structures after foreground model fits match qualitatively between simulated data and observations, suggesting that the frequency dependence of the beam is reasonably represented by the models. We find that a soil conductivity of 0.02 S m{sup -1} and relative permittivity of 3.5 yield good agreement between simulated spectra and observations. This is consistent with the soil properties reported by Sutinjo et al. for the Murchison Radio-astronomy Observatory, where EDGES is located.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Unsteady Euler algorithm with unstructured dynamic mesh for complex-aircraft aeroelastic analysis

A finite-volume unstructured-grid FEM scheme with multistage Runge-Kutta time stepping is applied to the three-dimensional time-dependent Euler equations for inviscid flows on complex aircraft configurations undergoing structural deformation. The derivation of the model, the solution procedure, and the computer implementation are described, and results are presented graphically for a NASA Langley supersonic fighter aircraft model in steady and unsteady (harmonic oscillation in complete-vehicle bending mode) flow regimes. Good agreement between FEM predictions and experimental data is demonstrated.

Batina, John T.↗

Radiation image reconstruction and uncertainty quantification using a Gaussian process prior

We propose a complete framework for Bayesian image reconstruction and uncertainty quantification based on a Gaussian process prior (GPP) to overcome limitations of maximum likelihood expectation maximization (ML-EM) image reconstruction algorithm. The prior distribution is constructed with a zero-mean Gaussian process (GP) with a choice of a covariance function, and a link function is used to map the Gaussian process to an image. Unlike many other maximum a posteriori approaches, our method offers highly interpretable hyperparamters that are selected automatically with the empirical Bayes method. Furthermore, the GP covariance function can be modified to incorporate a priori structural priors, enabling multi-modality imaging or contextual data fusion. Lastly, we illustrate that our approach lends itself to Bayesian uncertainty quantification techniques, such as the preconditioned Crank–Nicolson method and the Laplace approximation. The proposed framework is general and can be employed in most radiation image reconstruction problems, and we demonstrate it with simulated free-moving single detector radiation source imaging scenarios. We compare the reconstruction results from GPP and ML-EM, and show that the proposed method can significantly improve the image quality over ML-EM, all the while providing greater understanding of the source distribution via the uncertainty quantification capability. Furthermore, significant improvement of the image quality by incorporating a structural prior is illustrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A vectorized Lanczos eigensolver for high-performance computers

The computational strategies used to implement a Lanczos-based-method eigensolver on the latest generation of supercomputers are described. Several examples of structural vibration and buckling problems are presented that show the effects of using optimization techniques to increase the vectorization of the computational steps. The data storage and access schemes and the tools and strategies that best exploit the computer resources are presented. The method is implemented on the Convex C220, the Cray 2, and the Cray Y-MP computers. Results show that very good computation rates are achieved for the most computationally intensive steps of the Lanczos algorithm and that the Lanczos algorithm is many times faster than other methods extensively used in the past.

Bostic, Susan W.↗

Cost - The challenge for advanced materials and structures

Information is presented on the cost of various aircraft structures, together with methods for predicting and reducing cost. The need for the development of cost models, and of a comparative cost algorithm which could function as an engineering design tool to evaluate different design concepts, is emphasized. Efforts are underway to develop cost models that establish building-block unit cell elements that represent different material forms, geometric shapes, fabrication processes, and methods of assembly, with the purpose of expressing cost per pound or labor per pound data, with physical design and manufacture variables that a designer can visualize.

Davis, John G., Jr.↗

TRMM Field Campaigns: Objectives and Status Report

The Tropical Rainfall Measuring Mission (TRMM) satellite has been sending valuable data since launch in November 1997. Some of the key goals of the joint NASA (US) and NASDA (Japan) mission are: (1) to estimate the four-dimensional diabatic heating in the tropical and subtropical atmosphere, (2) understand the role of latent heating in driving tropical and extratropical circulations, (3) obtain monthly area-averaged estimates of rainfall over the data-sparse oceans, and (4) estimate the relative contribution of convective and stratiform precipitation over different regions during different seasons. The overarching scientific objective is to understand and improve estimates of rainfall and latent heating profiles throughout the global tropics. This requires observations for fundamental understanding of cloud dynamics and microphysics, as well as for validation, testing assumptions and error estimates of cloud-resolving models, forward radiative transfer models, algorithms used to estimate rainfall statistics and vertical structure of precipitation from surface-based radar, and from satellites. Field experiments designed to contribute to this understanding have been conducted in Texas and the South China Sea in spring of 1998, Florida in summer of 1998, and interior Brazil in (boreal) winter 1999. In summer 1999, a major oceanic campaign will be based at Kwajalein Atoll. Some early results will be highlighted, noting some significant contrasts between oceanic and continental convective systems.

Zipser, Edward I.↗

Clustering of the Diffuse Infrared Light from the COBE DIRBE Maps: Power Spectrum Analysis and Excess Isotropic Component of Fluctuations - 3

The cosmic infrared background (CIB) radiation is the cosmic repository for energy release throughout the history of the universe. The spatial fluctuations of the CIB resulting from galaxy clustering are expected to be at least a few percent on scales of a degree, depending on the luminosity and clustering history of the early universe. Using the all-sky data from the COBE DIRBE instrument at wavelengths 1.25 - 100 microns we attempt to measure the CIB fluctuations. In the near-IR, foreground emission is dominated by small scale structure due to stars in the Galaxy. There we find a strong correlation between the amplitude of the fluctuations and Galactic latitude after removing bright foreground stars. Using data outside the Galactic plane (absolute value of b > 20 deg) and away from the center (90 deg < l < 270 deg) we extrapolate the amplitude of the fluctuations to cosec absolute value of b = 0. We find a positive intercept of delta.F(sub rms) = 15.5(sup +3.7, sub -7.0), 5.9(sup +1.6, sub -3.7), 2.4(sup +0.5, sub -0.9), 2.0(sup +0.25, sub -0.5) nW/sq m.sr at 1.25, 2.2, 3.5 and 4.9 microns respectively, where the errors are the range of 92% confidence limits. For color subtracted maps between band 1 and 2 we find the isotropic part of the fluctuations at 7.6(sup +1.2, sub -2.4) nW/sq m.sr. Based on detailed numerical and analytic models, this residual is not likely to originate from the Galaxy, our clipping algorithm, or instrumental noise. We demonstrate that the residuals from the fit used in the extrapolation are distributed isotropically and suggest that this extra variance may result from structure in the CIB. We also obtain a positive intercept from a linear combination of maps at 1.25 and 2.2 microns. For 2 deg < theta < 15 deg, a power-spectrum analysis yields limits of (theta/5 deg) x delta.F(sub rms)(theta) < 6, 2.5, 0.8, 0.5 nW/sq m.sr at 1.25, 2.2, 3.5 and 4.9 microns respectively. From 10 - 100 microns, the dominant foregrounds are emission by dust in the Solar system and the Galaxy. There the upper limits on the CIB fluctuations are below 1 nW/sq m.sr and are lowest (< equal 0.5 nW/sq m.sr) at 25 microns.

Kashlinsky, A.↗

Multi-GPU immersed boundary method hemodynamics simulations

Large-scale simulations of blood flow that resolve the 3D deformation of each comprising cell are increasingly popular owing to algorithmic developments in conjunction with advances in compute capability. Among different approaches for modeling cell-resolved hemodynamics, fluid structure interaction (FSI) algorithms based on the immersed boundary method are frequently employed for coupling separate solvers for the background fluid and the cells within one framework. GPUs can accelerate these simulations; however, both current pre-exascale and future exascale CPU-GPU heterogeneous systems face communication challenges critical to performance and scalability. In this paper, we describe, to our knowledge, the largest distributed GPU-accelerated FSI simulations of high hematocrit cell-resolved flows with over 17 million red blood cells. We compare scaling on a fat node system with six GPUs per node and on a system with a single GPU per node. Through comparison between the CPU- and GPU-based implementations, we identify the costs of data movement in multiscale multi-grid FSI simulations on heterogeneous systems and show it to be the greatest performance bottleneck on the GPU.

97 MATHEMATICS AND COMPUTING↗

CP decomposition for tensors via alternating least squares with QR decomposition

The CP tensor decomposition is used in applications such as machine learning and signal processing to discover latent low-rank structure in multidimensional data. Computing a CP decomposition via an alternating least squares (ALS) method reduces the problem to several linear least squares problems. The standard way to solve these linear least squares subproblems is to use the normal equations, which inherit special tensor structure that can be exploited for computational efficiency. However, the normal equations are sensitive to numerical ill-conditioning, which can compromise the results of the decomposition. In this paper, we develop versions of the CP-ALS algorithm using the QR decomposition and the singular value decomposition, which are more numerically stable than the normal equations, to solve the linear least squares problems. Our algorithms utilize the tensor structure of the CP-ALS subproblems efficiently, have the same complexity as the standard CP-ALS algorithm when the input is dense and the rank is small, and are shown via examples to produce more stable results when ill-conditioning is present. Our MATLAB implementation achieves the same running time as the standard algorithm for small ranks, and we show that the new methods can obtain lower approximation error.

97 MATHEMATICS AND COMPUTING↗

Achieving precise multiparameter measurements with distributed optical fiber sensor using wavelength diversity and deep neural networks

The development of advanced distributed optical fiber sensing systems that are capable of performing accurate and spatially resolved multiparameter measurements is of great interest to a wide range of scientific and industrial applications. Here, in this paper, we propose and experimentally demonstrate a wavelength diversity based advanced distributed optical fiber sensor system to accomplish multiparameter sensing while greatly enhancing measurement accuracy. A suite of deep neural network (DNN) algorithms are developed and verified for data denoising, rapid Brillouin frequency shift estimation, and vibration data event classification. As a proof-of-concept, we demonstrate the effectiveness of the proposed advanced wavelength diversity distributed fiber sensor system assisted by DNN for simultaneous, independent measurements of static strain, temperature, and acoustic vibrations over a 25 km long sensing fiber at 3 m spatial resolution. These results suggest the potential for an intelligent multiparameter monitoring system with enhanced performance in advanced structural health monitoring applications.

47 OTHER INSTRUMENTATION↗