Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computing methods”

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,747 records · Page 97

Measuring Positions of Objects using Two or More Cameras

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

Klinko, Steve↗

Development of a three dimensional numerical water quality model for continental shelf applications

A model to predict the distribution of water quality parameters in three dimensions was developed. The mass transport equation was solved using a non-dimensional vertical axis and an alternating-direction-implicit finite difference technique. The reaction kinetics of the constituents were incorporated into a matrix method which permits computation of the interactions of multiple constituents. Methods for the computation of dispersion coefficients and coliform bacteria decay rates were determined. Numerical investigations of dispersive and dissipative effects showed that the three-dimensional model performs as predicted by analysis of simpler cases. The model was then applied to a two dimensional vertically averaged tidal dynamics model for the Providence River. It was also extended to a steady state application by replacing the time step with an iteration sequence. This modification was verified by comparison to analytical solutions and applied to a river confluence situation.

Spaulding, M.↗

Approximate method for predicting the permanent set in a beam in vacuo and in water subject to a shock wave

An approximate method to compute the maximum deformation and permanent set of a beam subjected to shock wave laoding in vacuo and in water was investigated. The method equates the maximum kinetic energy of the beam (and water) to the elastic plastic work done by a static uniform load applied to a beam. Results for the water case indicate that the plastic deformation is controlled by the kinetic energy of the water. The simplified approach can result in significant savings in computer time or it can expediently be used as a check of results from a more rigorous approach. The accuracy of the method is demonstrated by various examples of beams with simple support and clamped support boundary conditions.

Stiehl, A. L.↗

Paneling techniques for use with the VORLAX computer program

A method is presented for determining the geometric input data required by the VORLAX computer program in order to accurately model an aircraft configuration. Techniques are described for modeling each of the major components of a configuration and for joining these individual components into a complete configuration. The effects of trailing vortex filaments and methods of avoiding their intersection with downstream panels are also discussed. The methods presented here are applicable to most conventional aircraft configurations.

Martin, G. L.↗

Steady, subsonic, lifting surface theory for wings with swept, partial span, trailing edge control surfaces

A method for computing the lifting pressure distribution on a wing with partial span, swept control surfaces is presented. This method is valid within the framework of linearized, steady, potential flow theory and consists of using conventional lifting surface theory in conjuction with a flap pressure mode. The cause of a numerical instability that can occur during the quadrature of the flap pressure mode is discussed, and an efficient technique to eliminate the instability is derived. This technique is valid for both the flap pressure mode and regular pressure modes and could be used to improve existing lifting surface methods. Examples of the use of the flap pressure mode and comparisons among this method, other theoretical methods, and experiments are given. Discrepancies with experiment are indicated and candidate causes are presented. It is concluded that the method can lead to an efficient and accurate solution of the mathematical problem when a partial span, trailing edge flap is involved.

Medan, R. T.↗

Calculation of laminar and turbulent boundary layers for two-dimensional time-dependent flows

A general method for computing laminar and turbulent boundary layers for two-dimensional time-dependent flows is presented. The method uses an eddy-viscosity formulation to model the Reynolds shear-stress term and a very efficient numerical method to solve the governing equations. The model was applied to steady two-dimensional and three-dimensional flows and was shown to give good results. A discussion of the numerical method and the results obtained by the present method for both laminar and turbulent flows are discussed. Based on these results, the method is efficient and suitable for solving time-dependent laminar and turbulent boundary layers.

Cebeci, T.↗

Global accuracy estimates of point and mean undulation differences obtained from gravity disturbances, gravity anomalies and potential coefficients

Through the method of truncation functions, the oceanic geoid undulation is divided into two constituents: an inner zone contribution expressed as an integral of surface gravity disturbances over a spherical cap; and an outer zone contribution derived from a finite set of potential harmonic coefficients. Global, average error estimates are formulated for undulation differences, thereby providing accuracies for a relative geoid. The error analysis focuses on the outer zone contribution for which the potential coefficient errors are modeled. The method of computing undulations based on gravity disturbance data for the inner zone is compared to the similar, conventional method which presupposes gravity anomaly data within this zone.

Jekeli, C.↗

A stereo model based upon mechanisms of human binocular vision

A model for stereo vision, which is based on the human-binocular vision system, is proposed. Data collected from studies of neurophysiology of the human binocular system are discussed. An algorithm for the implementation of this stereo vision model is derived. The algorithm is tested on computer-generated and real scene images. Examples of a computer-generated image and a grey-level image are presented. It is noted that the proposed method is computationally efficient for depth perception, and the results indicate accuracies that are noise tolerant.

Griswold, N. C.↗

Thin film ferroelectric electro-optic memory

An electrically programmable, optically readable data or memory cell is configured from a thin film of ferroelectric material, such as PZT, sandwiched between a transparent top electrode and a bottom electrode. The output photoresponse, which may be a photocurrent or photo-emf, is a function of the product of the remanent polarization from a previously applied polarization voltage and the incident light intensity. The cell is useful for analog and digital data storage as well as opto-electric computing. The optical read operation is non-destructive of the remanent polarization. The cell provides a method for computing the product of stored data and incident optical data by applying an electrical signal to store data by polarizing the thin film ferroelectric material, and then applying an intensity modulated optical signal incident onto the thin film material to generate a photoresponse therein related to the product of the electrical and optical signals.

Thakoor, Sarita↗

Program for calculating optimum dimensions of alpha radioisotope capsules exposed to varying stress and temperature

A method and computer program were developed for calculating the creep and optimizing the dimensions of capsules filled with alpha-emitting radioisotopes. The method solves an integral equation that was developed assuming linear accumulation of partial creep lives and relating life to time-dependent stress and temperature using the Larson-Miller parameter. The computer program, CAPSUL, is written in Fortran language for the IBM 360/75 computer. The program makes a least squares fit of the creep life function using conventional constant stress, constant temperature creep data. Dimensions of capsules having maximum thermal power per unit of weight, volume, or area are calculated for a given creep life and pressure-temperature history using a numerical Lagrange Multiplier formulation. The program also calculates the life to a prescribed strain for capsules of given dimensions and pressure-temperature history. The method has been used to analyze creep data for the alloys 304 stainless steel, Hastelloy N, Cb-1% Zr, FS-85, and T-222.

J. P. Nichols↗

Combining Ratio Estimation for Low Density Parity Check (LDPC) Coding

The Low Density Parity Check (LDPC) Code decoding algorithm make use of a scaled receive signal derived from maximizing the log-likelihood ratio of the received signal. The scaling factor (often called the combining ratio) in an AWGN channel is a ratio between signal amplitude and noise variance. Accurately estimating this ratio has shown as much as 0.6 dB decoding performance gain. This presentation briefly describes three methods for estimating the combining ratio: a Pilot-Guided estimation method, a Blind estimation method, and a Simulation-Based Look-Up table. The Pilot Guided Estimation method has shown that the maximum likelihood estimates of signal amplitude is the mean inner product of the received sequence and the known sequence, the attached synchronization marker (ASM) , and signal variance is the difference of the mean of the squared received sequence and the square of the signal amplitude. This method has the advantage of simplicity at the expense of latency since several frames worth of ASMs. The Blind estimation method s maximum likelihood estimator is the average of the product of the received signal with the hyperbolic tangent of the product combining ratio and the received signal. The root of this equation can be determined by an iterative binary search between 0 and 1 after normalizing the received sequence. This method has the benefit of requiring one frame of data to estimate the combining ratio which is good for faster changing channels compared to the previous method, however it is computationally expensive. The final method uses a look-up table based on prior simulated results to determine signal amplitude and noise variance. In this method the received mean signal strength is controlled to a constant soft decision value. The magnitude of the deviation is averaged over a predetermined number of samples. This value is referenced in a look up table to determine the combining ratio that prior simulation associated with the average magnitude of the deviation. This method is more complicated than the Pilot-Guided Method due to the gain control circuitry, but does not have the real-time computation complexity of the Blind Estimation method. Each of these methods can be used to provide an accurate estimation of the combining ratio, and the final selection of the estimation method depends on other design constraints.

Mahmoud, Saad↗

Aerodynamic interaction between vortical wakes and lifting two-dimensional bodies

Unsteady rotor wake interactions with the empennage, tail boom, and other aerodynamic surfaces of a helicopter have a significant influence on its aerodynamic performance, the ride quality, and vibration. A numerical method for computing the aerodynamic interaction between an interacting vortex wake and the viscous flow about arbitrary two-dimensional bodies was developed to address this helicopter problem. The method solves for the flow field velocities on a body-fitted computational mesh using finite-difference techniques. The interacting vortex wake is represented by an array of discrete vortices which, in turn, are represented by a finite-core model. The evolution of the interacting vortex wake is calculated by Lagrangian techniques. The viscous flow field of the two-dimensional body is calculated on an Eulerian grid. The flow around circular and elliptic cylinders in the absence of an interacting vortex wake was calculated. These results compare very well with other numerical results and with results obtained from experiment and thereby demonstrate the accuracy of the viscous solution. The interaction of a rotor wake with the flow about a 4 to 1 elliptic cylinder at 45 degree incidence was calculated for a Reynolds number of 3000. The results demonstrate the significant variations in the lift and drag on the elliptic cylinder in the presence of the interacting rotor wake.

Stremel, Paul M.↗

An Efficient Inverse Aerodynamic Design Method For Subsonic Flows

Computational Fluid Dynamics based design methods are maturing to the point that they are beginning to be used in the aircraft design process. Many design methods however have demonstrated deficiencies in the leading edge region of airfoil sections. The objective of the present research is to develop an efficient inverse design method which is valid in the leading edge region. The new design method is a streamline curvature method, and a new technique is presented for modeling the variation of the streamline curvature normal to the surface. The new design method allows the surface coordinates to move normal to the surface, and has been incorporated into the Constrained Direct Iterative Surface Curvature (CDISC) design method. The accuracy and efficiency of the design method is demonstrated using both two-dimensional and three-dimensional design cases.

Milholen, William E., II↗

Aerodynamic interaction between vortical wakes and lifting two-dimensional bodies

Unsteady rotor wake interactions with the empenage, tail boom, and other aerodynamic surfaces of a helicopter have a significant influence on its aerodynamic performance, the ride quality, and amount of vibration. A numerical method for computing the aerodynamic interaction between an interacting vortex wake and the viscous flow about arbitrary two-dimensional bodies has been developed to address this helicopter problem. The method solves for the flow field velocities on a body-fitted computational mesh using finite-difference techniques. The interaction of a rotor wake with the flow about a 4:1 elliptic cylinder at 45-deg incidence was calculated for a Reynolds number of 3000.

Stremel, Paul M.↗

Semi-Markov Unreliability-Range Evaluator

Reconfigurable, fault-tolerant systems modeled. Semi-Markov unreliability-range evaluator (SURE) computer program is software tool for analysis of reliability of reconfigurable, fault-tolerant systems. Based on new method for computing death-state probabilities of semi-Markov model. Computes accurate upper and lower bounds on probability of failure of system. Written in PASCAL.

Butler, Ricky W.↗

Body conforming grids for general unsteady airfoil motion

A method for computing a new body conforming grid at each time step of an unsteady flow calculation is presented. The method is based on a spring system analogy. It simply updates the grid from the previous time step, hence it requires that an initial grid be created using another method. The method is designed to be independent of the flow equations used, so it can be easily transferred between aerodynamic codes with the same type of initial grid (i.e., O-type). The initial development was done using an O-type grid. It was then extended to C-type grids which required only minor changes to the routine. Results are presented for both O-type and C-type grids. The quality of the updated grids were examined by observing their movement and the flow results they produced. The effect due to the shape of the prescribed grid lines for the O-type grids was examined for a NACA 0012 airfoil undergoing harmonic pitch. The time history of the pressure coefficients of this case was compared with the interpolation scheme and experimental data.

Levine, M.↗

A comparative study and validation of upwind and central-difference Navier-Stokes codes for high-speed flows

A comparative study was made using 4 different computer codes for solving the compressible Navier-Stokes equations. Three different test problems were used, each of which has features typical of high speed internal flow problems of practical importance in the design and analysis of propulsion systems for advanced hypersonic vehicles. These problems are the supersonic flow between two walls, one of which contains a 10 deg compression ramp, the flow through a hypersonic inlet, and the flow in a 3-D corner formed by the intersection of two symmetric wedges. Three of the computer codes use similar recently developed implicit upwind differencing technology, while the fourth uses a well established explicit method. The computed results were compared with experimental data where available.

Rudy, David H.↗

Viscous computations of cold air/air flow around scramjet nozzle afterbody

The flow field in and around the nozzle afterbody section of a hypersonic vehicle was computationally simulated. The compressible, Reynolds averaged, Navier Stokes equations were solved by an implicit, finite volume, characteristic based method. The computational grids were adapted to the flow as the solutions were developing in order to improve the accuracy. The exhaust gases were assumed to be cold. The computational results were obtained for the two dimensional longitudinal plane located at the half span of the internal portion of the nozzle for over expanded and under expanded conditions. Another set of results were obtained, where the three dimensional simulations were performed for a half span nozzle. The surface pressures were successfully compared with the data obtained from the wind tunnel tests. The results help in understanding this complex flow field and, in turn, should help the design of the nozzle afterbody section.

Baysal, Oktay↗