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,441 records · Page 80

Structure Identification Within a Transitioning Swept-Wing Boundary Layer

Extensive measurements are made in a transitioning swept-wing boundary layer using hot-film, hot-wire and cross-wire anemometry. The crossflow-dominated flow contains stationary vortices that breakdown near mid-chord. The most amplified vortex wavelength is forced by the use of artificial roughness elements near the leading edge. Two-component velocity and spanwise surface shear-stress correlation measurements are made at two constant chord locations, before and after transition. Streamwise surface shear stresses are also measured through the entire transition region. Correlation techniques are used to identify stationary structures in the laminar regime and coherent structures in the turbulent regime. Basic techniques include observation of the spatial correlations and the spatially distributed auto-spectra. The primary and secondary instability mechanisms are identified in the spectra in all measured fields. The primary mechanism is seen to grow, cause transition and produce large-scale turbulence. The secondary mechanism grows through the entire transition region and produces the small-scale turbulence. Advanced techniques use Linear Stochastic Estimation (LSE) and Proper Orthogonal Decomposition (POD) to identify the spatio-temporal evolutions of structures in the boundary layer. LSE is used to estimate the instantaneous velocity fields using temporal data from just two spatial locations and the spatial correlations. Reference locations are selected using maximum RMS values to provide the best available estimates. POD is used to objectively determine modes characteristic of the measured flow based on energy. The stationary vortices are identified in the first laminar modes of each velocity component and shear component. Experimental evidence suggests that neighboring vortices interact and produce large coherent structures with spanwise periodicity at double the stationary vortex wavelength. An objective transition region detection method is developed using streamwise spatial POD solutions which isolate the growth of the primary and secondary instability mechanisms in the first and second modes, respectively. Temporal evolutions of dominant POD modes in all measured fields are calculated. These scalar POD coefficients contain the integrated characteristics of the entire field, greatly reducing the amount of data to characterize the instantaneous field. These modes may then be used to train future flow control algorithms based on neural networks.

Chapman, Keith↗

Analytical Verifications in Cryogenic Testing of NGST Advanced Mirror System Demonstrators

Ground based testing is a critical and costly part of component, assembly, and system verifications of large space telescopes. At such tests, however, with integral teamwork by planners, analysts, and test personnel, segments can be included to validate specific analytical parameters and algorithms at relatively low additional cost. This paper opens with strategy of analytical verification segments added to vacuum cryogenic testing of Advanced Mirror System Demonstrator (AMSD) assemblies. These AMSD assemblies incorporate material and architecture concepts being considered in the Next Generation Space Telescope (NGST) design. The test segments for workmanship testing, cold survivability, and cold operation optical throughput are supplemented by segments for analytical verifications of specific structural, thermal, and optical parameters. Utilizing integrated modeling and separate materials testing, the paper continues with support plan for analyses, data, and observation requirements during the AMSD testing, currently slated for late calendar year 2002 to mid calendar year 2003. The paper includes anomaly resolution as gleaned by authors from similar analytical verification support of a previous large space telescope, then closes with draft of plans for parameter extrapolations, to form a well-verified portion of the integrated modeling being done for NGST performance predictions.

Cummings, Ramona↗

Comparison of Modal Analysis Methods Applied to a Vibro-Acoustic Test Article

Modal testing of a vibro-acoustic test article referred to as the Aluminum Testbed Cylinder (ATC) has provided frequency response data for the development of validated numerical models of complex structures for interior noise prediction and control. The ATC is an all aluminum, ring and stringer stiffened cylinder, 12 feet in length and 4 feet in diameter. The cylinder was designed to represent typical aircraft construction. Modal tests were conducted for several different configurations of the cylinder assembly under ambient and pressurized conditions. The purpose of this paper is to present results from dynamic testing of different ATC configurations using two modal analysis software methods: Eigensystem Realization Algorithm (ERA) and MTS IDEAS Polyreference method. The paper compares results from the two analysis methods as well as the results from various test configurations. The effects of pressurization on the modal characteristics are discussed.

Pritchard, Jocelyn↗

Unsteady Aerodynamic Force Sensing from Strain Data

A simple approach for computing unsteady aerodynamic forces from simulated measured strain data is proposed in this study. First, the deflection and slope of the structure are computed from the unsteady strain using the two-step approach. Velocities and accelerations of the structure are computed using the autoregressive moving average model, on-line parameter estimator, low-pass filter, and a least-squares curve fitting method together with analytical derivatives with respect to time. Finally, aerodynamic forces over the wing are computed using modal aerodynamic influence coefficient matrices, a rational function approximation, and a time-marching algorithm.

aerodynamic load sensing↗

Baryon–Antibaryon Photoproduction off the Proton

The GlueX Experiment at Jefferson Lab uses linearly polarized photon beam to study the photoproduction of light mesons and baryons off a proton target with unprecedented precision. We have observed p¯p photoproduction and, for the first time, ?¯? photoproduction from thresholds up to E = 11.6 GeV. We reproduce the observed kinematic distributions with a phenomenological model that includes t-channel Regge exchange and double Regge exchange. The model uses few parameters to match all the kinematic spectra via an optimization algorithm developed specifically for this purpose. We report the total cross sections of the two reactions, as well as differential cross sections with respect to various momentum transfers and invariant mass spectra of p¯p, ?¯?, p¯? systems which show enhancements right above production thresholds. Furthermore, the beam asymmetries of the p¯p system are extracted from the GlueX data to further probe the nature of the t-channel Regge exchanges. We also report preliminary measurements of the spin polarizations of the individual hyperons. No prominent resonance structure was found in these reaction channels, but the angular distributions and the correlation of longitudinal momenta indicate the existence of multiple production mechanisms. Besides a dominant forward peak in the polar angle of the baryon-antibaryon pair in t-channel production, we see wide-angle antibaryon distributions different from that of the baryon. In the hyperon channel, we see a clear separation between the photoproduction of the ?¯?, p¯? systems.

Li, Hao↗

Experimental Characterization of Hysteresis in a Revolute Joint for Precision Deployable Structures

Recent studies of the micro-dynamic behavior of a deployable telescope metering truss have identified instabilities in the equilibrium shape of the truss in response to low-energy dynamic loading. Analyses indicate that these micro-dynamic instabilities arise from stick-slip friction within the truss joints (e.g., hinges and latches). The present study characterizes the low-magnitude quasi-static load cycle response of the precision revolute joints incorporated in the deployable telescope metering truss, and specifically, the hysteretic response of these joints caused by stick-slip friction within the joint. Detailed descriptions are presented of the test setup and data reduction algorithms, including discussions of data-error sources and data-filtering techniques. Test results are presented from thirteen specimens, and the effects of joint preload and manufacturing tolerances are investigated. Using a simplified model of stick-slip friction, a relationship is made between joint load-cycle behavior and micro-dynamic dimensional instabilities in the deployable telescope metering truss.

Lake, Mark S.↗

Scalable implementation of polynomial filtering for density functional theory calculation in PARSEC

In this work, we present an efficient implementation of polynomial filtering methods in PARSEC, a real-space pseudopotential based Kohn–Sham density functional theory solver. The implementation described here improves upon a Chebyshev-filtered subspace iteration algorithm used in the previous version of PARSEC. We present a hybrid polynomial filtering scheme that combines Chebyshev-filtered subspace iteration and a spectrum slicing method that partitions the spectrum into several spectral slices and uses bandpass-filtered subspace iteration to compute approximate eigenpairs within each interior slice simultaneously. We describe a procedure to partition a spectrum and construct polynomial filters. We also discuss a number of practical issues such as the use of appropriate data layouts for carrying out the computation on a two-dimensional process grid and how to achieve good load balance by allocating an appropriate number of process groups to each spectral slice. Numerical examples are presented to demonstrate the effectiveness of the hybrid polynomial filtering method as well as the superior parallel scalability of spectrum slicing in comparison to that of Chebyshev-filtered subspace iteration.

97 MATHEMATICS AND COMPUTING↗

Extended dynamic mode decomposition for model reduction in fluid dynamics simulations

High computational cost and storage/memory requirements of fluid dynamics simulations constrain their usefulness as a predictive tool. Reduced-order models (ROMs) provide a viable solution to this challenge by extracting the key underlying dynamics of a complex system directly from data. We investigate the efficacy and robustness of an extended dynamic mode decomposition (xDMD) algorithm in constructing ROMs of three-dimensional cardiovascular computations. Focusing on the ROMs' accuracy in representation and interpolation, we relate these metrics to the truncation rank of singular value decomposition, which underpins xDMD and other approaches to ROM construction. Our key innovation is to relate the truncation rank to the singular values of the original flow problem. This result establishes a priori guidelines for the xDMD deployment and its likely success as a means of data compression and reconstruction of the system's dynamics from dominant spatiotemporal structures present in the data.

Mechanics↗

Stereo image visualization of numerically simulated turbulence

Stereo imaging techniques are used in the current study to analyze the results of a direct Navier-Stokes simulation of a turbulent boundary layer. Coherent structures embedded within the unsteady, highly three-dimensional turbulent flow fields are significantly enhanced through the use of stereo computer graphic renderings. In the paper, the geometrical foundations of two different methods for generating stereo images are reviewed. Sufficient detail is provided to enable readers to code stereo display algorithms from scratch. The effects of varying the geometric parameters of the stereo algorithm are demonstrated. Sample image pairs from animations of the turbulence numerical simulation are presented in color. Stereo imagery, especially when used with time-evolving data, provides new insight into the physical processes responsible for the generation an maintenance of turbulence within boundary layers.

Robinson, Stephen K.↗

Practical aspects of on-orbit modal identification using free-decay data

Future large space structures such as Space Station Freedom can be tested as a total structure on Earth. Size, zero g design, and the existing atmosphere conflict with traditional structural dynamic testing. On-orbit modal identification will become necessary. Practical aspects of performing such tests using time domain analysis of free decay data are discussed. The effects of environmental constraints, structural characteristics, and excitation and sensing are reviewed. A recent laboratory application and sensing are reviewed. In this test, an on-orbit experiment is simulated using free decay data and a limited number of excitation and measurement points. The test article is dynamically similar to the truss/solar arrays of Space Station. The identified modal parameters from this simulation correlate well, though not uniquely, with those obtained in a complete modal survey. Practical difficulties in performing the correlation are illustrated. Model parameters are identified with the Eigensystem Realization Algorithm.

Schenk, Axel↗

Ice tracking techniques, implementation, performance, and applications

Present techniques of ice tracking make use both of cross-correlation and of edge tracking, the former being more successful in heavy pack ice, the latter being critical for the broken ice of the pack margins. Algorithms must assume some constraints on the spatial variations of displacements to eliminate fliers, but must avoid introducing any errors into the spatial statistics of the measured displacement field. We draw our illustrations from the implementation of an automated tracking system for kinematic analyses of ERS-1 and JERS-1 SAR imagery at the University of Alaska - the Alaska SAR Facility's Geophysical Processor System. Analyses of the ice kinematic data that might have some general interest to analysts of cloud-derived wind fields are the spatial structure of the fields, and the evaluation and variability of average deformation and its invariants: divergence, vorticity and shear. Many problems in sea ice dynamics and mechanics can be addressed with the kinematic data from SAR.

Rothrock, D. A.↗

Optimized System Identification

In system identification, one usually cares most about finding a model whose outputs are as close as possible to the true system outputs when the same input is applied to both. However, most system identification algorithms do not minimize this output error. Often they minimize model equation error instead, as in typical least-squares fits using a finite-difference model, and it is seen here that this distinction is significant. Here, we develop a set of system identification algorithms that minimize output error for multi-input/multi-output and multi-input/single-output systems. This is done with sequential quadratic programming iterations on the nonlinear least-squares problems, with an eigendecomposition to handle indefinite second partials. This optimization minimizes a nonlinear function of many variables, and hence can converge to local minima. To handle this problem, we start the iterations from the OKID (Observer/Kalman Identification) algorithm result. Not only has OKID proved very effective in practice, it minimizes an output error of an observer which has the property that as the data set gets large, it converges to minimizing the criterion of interest here. Hence, it is a particularly good starting point for the nonlinear iterations here. Examples show that the methods developed here eliminate the bias that is often observed using any system identification methods of either over-estimating or under-estimating the damping of vibration modes in lightly damped structures.

Juang, Jer-Nan↗

Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to demonstrate the use of design and analysis of computer experiments (DACE) methods in Sandia’s DAKOTA software package for surrogate modeling and optimization. These methods were applied to a flow- path fueled with an interdigitated flushwall injector suitable for scramjet applications at hyper- velocity conditions and ascending along a constant dynamic pressure flight trajectory. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. Because the RAS of this case are computationally expensive, surrogate models are used for optimization. To build a surrogate model a RAS database is created. The sequence of the design variables comprising the database were generated using a Latin hypercube sampling (LHS) method. A methodology was also developed to automatically build geometries and generate structured grids for each design point. The ensuing RAS analysis generated the simulation database from which the two objective functions were computed using a one-dimensionalization (1D) of the three-dimensional simulation data. The data were fitted using four surrogate models: an artificial neural network (ANN), a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model predicted an optimal solution set that exhibited high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts may be required to lower the surrogate model errors and perform more accurate surrogate-model-based optimization.

Shenoy, Rajiv R.↗

Dense 3D-Reconstruction from Monocular Image Sequences for Computationally Constrained UAS

The ability to find safe landing sites over complex 3D terrain is an essential safety feature for fully autonomous small unmanned aerial systems (UAS), which requires on-board perception for 3D reconstruction and terrain analysis if the overflown terrain is unknown. This is a challenge for UAS that are limited in size, weight and computational power, such as small rotorcrafts executing autonomous missions on Earth, or in planetary applications such as the Mars Helicopter. For such a computationally constraint system, we propose a structure from motion approach that uses inputs from a single downward facing camera to produce dense point clouds of the overflown terrain in real time. In contrast to existing approaches, our method uses metric pose information from a visual-inertial odometry algorithm as camera pose priors, which allows deploying a fast pose refinement step to align camera frames such that a conventional stereo algorithm can be used for dense 3D reconstruction. We validate the performance of our approach with extensive evaluations in simulation, and demonstrate the feasibility with data from UAS flights.

Brockers, Roland↗

Integrated Modeling for Payload Test of the Roman Space Telescope

The Nancy Grace Roman Telescope (RST) is a NASA observatory designed to unravel the secrets of dark energy and dark matter, search for and image exoplanets, and explore many topics in infrared optics. Scheduled to launch in no earlier than October 2026, this 2.4 meter aperture telescope has a field of view 100 times greater than the Hubble Space Telescope. The mission is currently in its construction phase, where integrated modeling between thermal, structural, and optical models of the observatory is necessary to demonstrate science quality images over the range of operational parameters. This presentation discusses the most recent integrated modeling analysis cycle for Roman, including model correlation with our instrument level testing. We include a discussion on improved processes of the handling of the various flows of data between the modeling disciplines and discipline specific monte-carlo analysis predictions. We will finish with the predicted uncertainties and expected performance for our upcoming observatory alignment verification test using machine learning algorithms.

telescope↗

A model of asynchronous iterative algorithms for solving large, sparse, linear systems

Solving large, sparse, linear systems of equations is one of the fundamental problems in large scale scientific and engineering computation. A model of a general class of asynchronous, iterative solution methods for linear systems is developed. In the model, the system is solved by creating several cooperating tasks that each compute a portion of the solution vector. This model is then analyzed to determine the expected intertask data transfer and task computational complexity as functions of the number of tasks. Based on the analysis, recommendations for task partitioning are made. These recommendations are a function of the sparseness of the linear system, its structure (i.e., randomly sparse or banded), and dimension.

Reed, D. A.↗

Model-based Reconstruction for Single Particle Cryo-Electron Microscopy

Single particle cryo-electron microscopy is a vital tool for 3D characterization of protein structures. A typical workflow involves acquiring projection images of a collection of randomly oriented particles, picking and classifying individual particle projections by orientation, and finally using the individual particle projections to reconstruct a 3D map of the electron density profile. The reconstruction is challenging because of the low signal-to-noise ratio of the data, the unknown orientation of the particles, and the sparsity of data especially when dealing with flexible proteins where there may not be sufficient data corresponding to each class to obtain an accurate reconstruction using standard algorithms. In this paper we present a model-based image reconstruction technique that uses a regularized cost function to reconstruct the 3D density map by assuming known orientations for the particles. Our method casts the reconstruction as minimizing a cost function involving a novel forward model term that accounts for the contrast transfer function of the microscope, the orientation of the particles and the center of rotation offsets. We combine the forward model term with a regularizer that enforces desirable properties in the volume to be reconstructed. Using simulated data, we demonstrate how our method can significantly improve upon the typically used approach.

Venkatakrishnan, Singanallur↗

Gravitation and Mesh Adaption

The Gravitation and Mesh Adaptation (GaMA) toolbox consists of a set of Matlab classes for modeling the environment around asteroids and comets. A variety of gravitational models and supporting algorithms from are consolidated alongside a custom meshing utility tailored for the application. This allows the user to work within a single streamlined environment to import and manipulate surface definitions, probe the dynamical environment, integrate trajectories, and post-process results. For other applications, modified surface meshes can be exported. Gravity Models: 1) Werner's analytic polyhedron model 2) Mascon model 3) Gottlieb's spherical harmonic model 4) Approximate polyhedron models 5) Curvilinear surface models 6) Custom composite models Meshing Features: 1) Array-based half-edge data structure 2) Supports curvilinear surface definitions up to degree 4 3) Ray-tracing 4) Coarsening 5) Feature-based refinement 6) Projection 7) Smoothing 8) Mesh quality validity tests 9) Mesh repair Additional Features: 1) Solar radiation pressure model 2) Distant 3rd body model 3) Collision detection 4) Trajectory integration and post-processing 5) Visualization of surface fields

Pearl, Jason↗