Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordinate 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 919 records · Page 51

Algorithmic considerations of integrated design for CSI on a hypercube architecture

An approach is presented to the integrated design problem for actively controlled large, flexible mechanical systems for which Control Structure Interaction (CSI) problems are of concern. The two coupled design problems were identified as the optimal Structural Design problem the optimal Controller Design problem. These two problems can be addressed within a decision making loop that would consider each separately, and then sequentially analyze the effects of one on the other. Embedded in such a loop would be the simulation and coordination tasks as part of the decision tools required in a total (software) package. All of the above are compute-intensive tasks. In any such task, possible decompositions and gains due to the inherent parallelism have to be exploited. The problems under consideration, as applied to large flexible mechanical structures, are particularly suited to be mapped onto multicomputer systems in a hypercube topology.

Oezguener, UE.↗

Evaluation of lens distortion errors in video-based motion analysis

In an effort to study lens distortion errors, a grid of points of known dimensions was constructed and videotaped using a standard and a wide-angle lens. Recorded images were played back on a VCR and stored on a personal computer. Using these stored images, two experiments were conducted. Errors were calculated as the difference in distance from the known coordinates of the points to the calculated coordinates. The purposes of this project were as follows: (1) to develop the methodology to evaluate errors introduced by lens distortion; (2) to quantify and compare errors introduced by use of both a 'standard' and a wide-angle lens; (3) to investigate techniques to minimize lens-induced errors; and (4) to determine the most effective use of calibration points when using a wide-angle lens with a significant amount of distortion. It was seen that when using a wide-angle lens, errors from lens distortion could be as high as 10 percent of the size of the entire field of view. Even with a standard lens, there was a small amount of lens distortion. It was also found that the choice of calibration points influenced the lens distortion error. By properly selecting the calibration points and avoidance of the outermost regions of a wide-angle lens, the error from lens distortion can be kept below approximately 0.5 percent with a standard lens and 1.5 percent with a wide-angle lens.

Poliner, Jeffrey↗

Computational Materials Research

Computational Materials aims to model and predict thermodynamic, mechanical, and transport properties of polymer matrix composites. This workshop, the second coordinated by NASA Langley, reports progress in measurements and modeling at a number of length scales: atomic, molecular, nano, and continuum. Assembled here are presentations on quantum calculations for force field development, molecular mechanics of interfaces, molecular weight effects on mechanical properties, molecular dynamics applied to poling of polymers for electrets, Monte Carlo simulation of aromatic thermoplastics, thermal pressure coefficients of liquids, ultrasonic elastic constants, group additivity predictions, bulk constitutive models, and viscoplasticity characterization.

Hinkley, Jeffrey A.↗

Improvements in Block-Krylov Ritz Vectors and the Boundary Flexibility Method of Component Synthesis

A method of dynamic substructuring is presented which utilizes a set of static Ritz vectors as a replacement for normal eigenvectors in component mode synthesis. This set of Ritz vectors is generated in a recurrence relationship, proposed by Wilson, which has the form of a block-Krylov subspace. The initial seed to the recurrence algorithm is based upon the boundary flexibility vectors of the component. Improvements have been made in the formulation of the initial seed to the Krylov sequence, through the use of block-filtering. A method to shift the Krylov sequence to create Ritz vectors that will represent the dynamic behavior of the component at target frequencies, the target frequency being determined by the applied forcing functions, has been developed. A method to terminate the Krylov sequence has also been developed. Various orthonormalization schemes have been developed and evaluated, including the Cholesky/QR method. Several auxiliary theorems and proofs which illustrate issues in component mode synthesis and loss of orthogonality in the Krylov sequence have also been presented. The resulting methodology is applicable to both fixed and free- interface boundary components, and results in a general component model appropriate for any type of dynamic analysis. The accuracy is found to be comparable to that of component synthesis based upon normal modes, using fewer generalized coordinates. In addition, the block-Krylov recurrence algorithm is a series of static solutions and so requires significantly less computation than solving the normal eigenspace problem. The requirement for less vectors to form the component, coupled with the lower computational expense of calculating these Ritz vectors, combine to create a method more efficient than traditional component mode synthesis.

Carney, Kelly Scott↗

Format for Interchange and Display of 3D Terrain Data

Visible Scalable Terrain (ViSTa) is a software format for production, interchange, and display of three-dimensional (3D) terrain data acquired by stereoscopic cameras of robotic vision systems. ViSTa is designed to support scalability of data, accuracy of displayed terrain images, and optimal utilization of computational resources. In a ViSTa file, an area of terrain is represented, at one or more levels of detail, by coordinates of isolated points and/or vertices of triangles derived from a texture map that, in turn, is derived from original terrain images. Unlike prior terrain-image software formats, ViSTa includes provisions to ensure accuracy of texture coordinates. Whereas many such formats are based on 2.5-dimensional terrain models and impose additional regularity constraints on data, ViSTa is based on a 3D model without regularity constraints. Whereas many prior formats require external data for specifying image-data coordinate systems, ViSTa provides for the inclusion of coordinate-system data within data files. ViSTa admits highspeed loading and display within a Java program. ViSTa is designed to minimize file sizes and maximize compressibility and to support straightforward reduction of resolution to reduce file size for Internet-based distribution.

Backes, Paul↗

Tunguska Workshop: Applying Modern Tools to Understand the 1908 Tunguska Impact

A one-day workshop was held at NASA Ames Research Center, January 16, 2018, to re-examine the 1908 Tunguska impact using modern computational tools, many of them developed in response to the 2013 Chelyabinsk airburst. Twelve international experts gave presentations, with another 40 attending in-person or remotely. The most likely models for Tunguska converged on an energy of 10-20 Megatons, released in an airburst at a height of about 10 km. If the Tunguska impactor was a stony asteroid similar to Chelyabinsk, the diameter was roughly 50-80m. A comparison with current understanding of the population of asteroids in this size range indicates that the interval between such events is millennia, not centuries as had been concluded previously. The primary constraints on our understanding of Tunguska are the dearth of quantitative data, not weakness of the computational models. The workshop was sponsored by the NASA Ames Asteroid Threat Assessment Project and supported the NASA Planetary Defense Coordination Office.

Morrison, David↗

Coordinated observations of X-ray bright BL Lacertae objects

No new International Ultraviolet Explorer (IUE) observations were scheduled. The work on the spectral extraction routine was continued. The computer code was adapted from the RDAF PDP 11/44 computer to a VAX 11/750, in order that it run more quickly. A chapter on ultraviolet observations of blazars for a book celebrating the scientific accomplishment of IUE was written in colloboration. A copy of the final draft is attached.

Urry, C. M.↗

Algorithm for Wavefront Sensing Using an Extended Scene

A recently conceived algorithm for processing image data acquired by a Shack-Hartmann (SH) wavefront sensor is not subject to the restriction, previously applicable in SH wavefront sensing, that the image be formed from a distant star or other equivalent of a point light source. That is to say, the image could be of an extended scene. (One still has the option of using a point source.) The algorithm can be implemented in commercially available software on ordinary computers. The steps of the algorithm are the following: 1. Suppose that the image comprises M sub-images. Determine the x,y Cartesian coordinates of the centers of these sub-images and store them in a 2xM matrix. 2. Within each sub-image, choose an NxN-pixel cell centered at the coordinates determined in step 1. For the ith sub-image, let this cell be denoted as si(x,y). Let the cell of another subimage (preferably near the center of the whole extended-scene image) be designated a reference cell, denoted r(x,y). 3. Calculate the fast Fourier transforms of the sub-sub-images in the central NxN portions (where N < N and both are preferably powers of 2) of r(x,y) and si(x,y). 4. Multiply the two transforms to obtain a cross-correlation function Ci(u,v), in the Fourier domain. Then let the phase of Ci(u, v) constitute a phase function, phi(u,v). 5. Fit u and v slopes to phi (u,v) over a small u,v subdomain. 6. Compute the fast Fourier transform, Si(u,v) of the full NxN cell si(x,y). Multiply this transform by the u and phase slopes obtained in step 4. Then compute the inverse fast Fourier transform of the product. 7. Repeat steps 4 through 6 in an iteration loop, cumulating the u and slopes, until a maximum iteration number is reached or the change in image shift becomes smaller than a predetermined tolerance. 8. Repeat steps 4 through 7 for the cells of all other sub-images.

Sidick, Erkin↗

Latitude dependence of eddy variances

The eddy variance of a meteorological field must tend to zero at high latitudes due solely to the nature of spherical polar coordinates. The zonal averaging operator defines a length scale: the circumference of the latitude circle. When the circumference of the latitude circle is greater than the correlation length of the field, the eddy variance from transient eddies is the result of differences between statistically independent regions. When the circumference is less than the correlation length, the eddy variance is computed from points that are well correlated with each other, and so is reduced. The expansion of a field into zonal Fourier components is also influenced by the use of spherical coordinates. As is well known, a phenomenon of fixed wavelength will have different zonal wavenumbers at different latitudes. Simple analytical examples of these effects are presented along with an observational example from satellite ozone data. It is found that geometrical effects can be important even in middle latitudes.

Bowman, Kenneth P.↗

New Rover Conops with High-Performance Onboard Computing: Give Up Raw Data to Reduce Ops Cost and Do More Science

A major portion of time during the tactical operation of Mars rovers is spent for selecting, prioritizing, and coordinating sciences and engineering activities such that they fit within resource constraints, including the downlink data volume, energy, and time. In particular, the downlink data volume constraint is getting particularly tighter in recent missions because modern instruments produce increasingly high data volume while the communication bandwidth is essentially bounded by the law of physics. Tactical operation would be substantially simplified, hence the operation cost could be reduced, if the data volume constraint is relaxed or even removed. In this abstract, we propose a new operation paradigm for achieving this goal. The key observation is that, both in science and engineering applications, the bit size of raw data is typically much greater than the volume of processed information that is needed for scientific or engineering analysis. For example, a full-resolution image from Mastcam-Z, the main science camera on Perseverance, is about 700 kB in volume and we downlinked 29,685 images up to Sol 243, totaling ~20 GB of data. But of course, scientists do not use every pixel of these images; what they really look for in the images are geological features, typically represented by specific geometric configurations or textures. An end product after processing hundreds of Mascam-Z images could be a single geological map summarizing the spatial distribution of the features. For another example, a 100-meter drive of Perseverance produces 7-12 MB of drive telemetry, which records every detail of the rover's motion at 8 Hz, including position, attitude, steering angles, encoder readings, motor currents and many other information. But what the ground engineers eventually pay attention to is the signs of anomaly, such as excessive motor currents or high slip; if a drive is nominal, the vast majority of this data is unused. What if, then, we process the raw data onboard and only downlink the processed data that is relevant to scientific or engineering analyses, such as a list of detected science features (with cropped images) or a list of potential signs of anomaly while driving? A major roadblock for such onboard, high-level information processing has been the onboard computational resource. RAD750, the main onboard computer of Perseverance, is obviously not sufficient for performing complex image or signal processing such as object detection, semantic segmentation, or anomaly detection. Interestingly, RAD750 is not the best processor that Perseverance has; Qualcomm's Snapdragon 801, a modern mobile processor, is on her Heli Base Station, a device for communicating with Mars Helicopter Ingenuity; also, Intel's Atom E3845 processors are on engineering cameras. In the reminder of this paper, we will introduce two particular uses cases of these high-performance co-processors (meaning auxiliary CPU, GPU, or other types of processors that are separate from the main processor that runs the main flight software) for lowering operation cost and accommodating more science activities for a given communication constraint.

Didier, A.↗

Three-dimensional unsteady viscous flow analysis over airfoil sections

A three-dimensional solution procedure for the approximate form of the Navier-Stokes equation was exercised in the two- and three-dimensional modes to compute the unsteady turbulent boundary layer on a flat plate corresponding to the data of Karlsson. The procedure is based on the use of a consistently split Linearized Block Implicit technique in conjunction with a QR operator scheme. New time-dependent upstream boundary conditions were developed that yielded realistic solutions for the interior in the vicinity of the upstream boundary. Comparisons of the computation employing these boundary conditions with the data indicate that both qualitative and quantitative agreement was obtained for the mean velocity and the in phase and out of phase components of the first harmonic of the velocity. In addition, the calculation gave results for the skin friction phase angle that had expected physical behavior for large distances downstream of the inflow boundary. For the three-dimensional case, the two-dimensional data of Karlsson was considered, but in a coordinate system skewed at 45 deg to the free stream direction. The results of the calculations were in excellent agreement with the data and the two-dimensional computations.

Weinberg, B. C.↗

Computational study of duct and pipe flows using the method of pseudocompressibility

A viscous, three-dimensional, incompressible, Navier-Stokes Computational Fluid Dynamics code employing pseudocompressibility is used for the prediction of laminar primary and secondary flows in two 90-degree bends of constant cross section. Under study are a square cross section duct bend with 2.3 radius ratio and a round cross section pipe bend with 2.8 radius ratio. Sensitivity of predicted primary and secondary flow to inlet boundary conditions, grid resolution, and code convergence is investigated. Contour and velocity versus spanwise coordinate plots comparing prediction to experimental data flow components are shown at several streamwise stations before, within, and after the duct and pipe bends. Discussion includes secondary flow physics, computational method, computational requirements, grid dependence, and convergence rates.

Williams, Robert W.↗

Development of a three-dimensional supersonic inlet flow analysis

A method for computing three dimensional flow in supersonic inlets is described. An approximate set of governing equations is given for viscous flows which have a primary flow direction. The governing equations are written in general orthogonal coordinates. These equations are modified in the subsonic region of the flow to prevent the phenomenon of branching. Results are presented for the two sample cases: a Mach number equals 2.5 flow in a square duct, and a Mach number equals 3.0 flow in a research jet engine inlet. In the latter case the computed results are compared with the experimental data. A users' manual is included.

Buggeln, R. C.↗

Analysis and computation of three-dimensional flow in strongly curved ducts

Numerical solutions are presented for three-dimensional laminar and turbulent flow in curved ducts of rectangular cross section and significant curvature. The analysis is based on a primary-secondary velocity decomposition in a given coordinate system, and leads to approximate governing equations which correct an a priori inviscid solution for viscous effects, secondary flows, total pressure distortion, heat transfer, and internal flow blockage and losses. Solution of the correction equations is accomplished as an initial-value problem in space using an implicit forward-marching technique. The overall solution procedure requires significantly less computational effort than Navier-Stokes algorithms. The present solution procedure is effective even with the extreme local mesh resolution which is necessary to resolve near-wall sublayer regions in turbulent flow calculations. Computed solutions for both laminar and turbulent flow compare very favorably with available analytical and experimental results.

Kreskovsky, J. P.↗

Asymmetric hypersonic flow

A general method for the analysis of the inviscid asymmetric hypersonic flow fields enveloping smooth bodies of general shape is given. The method is based on the assumption of a thin shock layer which yields an expression for pressure in generalized Mises coordinates. Numerical results for elliptic cones at angle of attack are shown to compare well with experiments and other theories. The computing logic for a blunt body is described, and a limiting solution at the stagnation point is presented.

Maslen, S. H.↗

Control point measurements on Mariner 9 pictures

Measurements are given of the control points on the Mariner 9 pictures which were used in the computation of the June 1973 control net of Mars. The method of making the measurements is discussed along with the picture coordinate system and the removal of distortions. Table are presented of the 9804 measurements of 1645 points on 660 pictures.

Davies, M. E.↗

A system to geometrically rectify and map airborne scanner imagery and to estimate ground area

A system of computer programs were developed which performs geometric rectification and line-by-line mapping of airborne multispectral scanner data to ground coordinates and estimates ground area. The system requires aircraft attitude and positional information furnished by ancillary aircraft equipment, as well as ground control points. The geometric correction and mapping procedure locates the scan lines, or the pixels on each line, in terms of map grid coordinates. The area estimation procedure gives ground area for each pixel or for a predesignated parcel specified in map grid coordinates. The results of exercising the system with simulated data showed the uncorrected video and corrected imagery and produced area estimates accurate to better than 99.7%.

Spencer, M. M.↗