Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “approximate 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 451 records · Page 25

Space Station Active Thermal Control System modeling

The Space Station Active Thermal Control System (ATCS) has been modeled using modified SINDA/SINFLO programs to solve two-phase Thermo-fluid problems. The modifications include changes in several subroutines to incorporate implicit solution which allows larger time step as compared to that for explicit solutions. Larger time step saves computer time but involves larger computational error. Several runs were made using various time steps for the ATCS model. It has been found that for a reasonable approach, three times larger time step as compared to that used in explicit method is a good value which will reduce the computer time by approximately 50 percent and still maintain the accuracy of the output data to within 90 percent of the explicit values.

Hye, Abdul↗

Resonant decay of global MHD modes at 'thick' interfaces

A global surface mode can exist on a 'surface' which is not a true discontinuity, but without sufficient dissipation it is not a normal mode and it decays in time via resonance absorption. The decay rate can be calculated analytically when the 'surface' is thin. The goal of this paper is to numerically estimate the decay rate when the surface is not thin. A cold plasma in a uniform magnetic field is considered in which density varies linearly across the 'surface'. In this linearized calculation, the global surface mode is driven in steady state by an antenna located in one of the uniform regions external to the surface. The frequency is a free parameter, and resonance curves are computed numerically without approximations. The widths of the resonance curves are used to estimate the free decay times of undriven surface modes, via the uncertainty principle. When applied to the active solar corona, the decay rates are large enough to account for the coronal heating, but it should be kept in mind that the role of nonlinearity in resonance absorption is still undetermined.

Hollweg, Joseph V.↗

Probabilistic structural analysis: Introductory remarks

The development of probabilistic structural analysis methodology consists of the following program elements: (1) composite load spectra models, (2) computational probabilistic structural analysis methods, and (3) probabilistic constitutive relationships. The development of the probabilistic structural analysis methodology is a joint program of NASA Lewis in-house and sponsored research. The objective of this session is to illustrate recent progress on the application of this methodology to determine the reliability of structural components for rocket propulsion systems. The session contains descriptions of and progress reports on the following specific activities: (1) The NESSUS computer code, (2) approximate methods, (3) advanced methods, (4) composite load spectra applications, (5) probabilistic fracture mechanisms, and (6) probability of failure and risk analysis. Collectively, the progress to date demonstrates that the structural durability of hot engine structural components can be effectively evaluated in a formal probabilistic/reliability framework.

Chamis, Christos C.↗

Navier-Stokes calculation of transonic flow past the NTF 65-deg delta wing

Viscous flow past a wind tunnel model of a 65-degree swept angle Delta wing at transonic speeds is being studied. The model was tested in the 8-foot cryogenic transonic wind tunnel at the National Transonic Facility. Aerodynamic forces and wing surface pressure data were obtained at various angles of attack, Mach numbers, and Reynold's numbers for four different leading edges of the wing. The objectives of the present investigation are: (1) to perform numerical modeling of the flow around the wing; (2) to validate the experimental data with a Navier-Stokes computational fluid dynamics code and vice versa; (3) to investigate the effects of the sting mount of the wing; (4) to evaluate the effects of leading edge radius on the flow; and (5) to explain the Reynold's number effect as indicated by the test data. Several computer programs were developed to define the surfaces of the wing, the four leading edges, and the sting mount. Based on these geometric databases, the surface grids of a single-block computational domain was generated interactively on the IRIS workstation using the GRIDGEN2D module of GRIDGEN. To refine the grids and to avoid excessive loss of grid points due to collapsed edges, a 9-block computational domain containing approximately 750,000 grid points was developed with the GRIDBLOCK module to replace the single-block grid.

Wu, Chivey↗

Confidence bounds on structural reliability

Different approaches for quantifying physical, statistical, and model uncertainties associated with the distribution parameters which are aimed at determining structural reliability are described. Confidence intervals on the distribution parameters of the input random variables are estimated using four algorithms to evaluate uncertainty of the response. Design intervals are evaluated using either Monte Carlo simulation or an iterative approach. A first order approach can be used to compute a first approximation of the design interval, but its accuracy is not satisfactory. The regression approach which combines the iterative approach with Monte Carlo simulation is capable of providing good results if the performance function can be accurately represented using regression analysis. It is concluded that the design interval-based approach seems to be quite general and takes into account distribution and model uncertainties.

Mehta, S. R.↗

L-band radar scattering from grass

A radar system based on a network analyzer has been developed to study the backscatter from vegetation. The radar is operated at L-band. Radar measurements of a grass field were made in 1991. The radar returns from the grass were measured at three incidence angles. Ground truth and canopy parameters such as blade and stem dimensions, moisture content of the grass and the soil, and blade and stem density, were measured. These parameters are used in a distorted Born approximation model to compute the backscatter coefficients from the grass layer. The model results are compared with the radar data.

Chauhan, N.↗

Computation of the two-dimensional flow in a laminar boundary layer

A comparison is made of the boundary-layer flow computed by the approximate method developed by K. Pohlhausen with the exact solutions which have been published for several special cases. A modification of Pohlhausen's method has been developed which extends the range of application at the expense of some decrease in the accuracy of the approximations.

Dryden, Hugh L↗

Feasibility study of a single, elliptical heliocentric Earth-Mars trajectory

The initial intent of this design project was to evaluate the existence and feasibility of a single elliptical heliocentric Earth/Mars trajectory. This trajectory was constrained to encounter Mars twice in its orbit, within a time interval of 15 to 180 Earth days between encounters. The single ellipse restriction was soon found to be prohibitive for reasons shown later. Therefore, the approach taken in the design of the round-trip mission to Mars was to construct single-leg trajectories which connected two planets on two prescribed dates. Three methods of trajectory design were developed. Method 1 is an eclectic approach and employs Gaussian Orbit Determination (Method 1A) and Lambert-Euler Preliminary Orbit Determination (Method 1B) in conjunction with each other. Method 2 is an additional version of Lambert's Solution to orbit determination, and both a coplanar and a noncoplanar solution were developed within Method 2. In each of these methods, the fundamental variables are two position vectors and the time between the position vectors. In all methods, the motion was considered Keplerian motion and the reference frame origin was located at the sun. Perturbative effects were not considered in Method 1. The feasibility study of round-trip Earth/Mars trajectories maintains generality by considering only heliocentric trajectory parameters and planetary approach conditions. The coordinates and velocity components of the planets, for the standard epoch J2000, were computed from an approximate set of osculating elements by the procedure outlined in an ephemeris of coordinates.

Blake, M.↗

A sensitivity equation approach to shape optimization in fluid flows

A sensitivity equation method to shape optimization problems is applied. An algorithm is developed and tested on a problem of designing optimal forebody simulators for a 2D, inviscid supersonic flow. The algorithm uses a BFGS/Trust Region optimization scheme with sensitivities computed by numerically approximating the linear partial differential equations that determine the flow sensitivities. Numerical examples are presented to illustrate the method.

Borggaard, Jeff↗

Algebraic Multigrid by Smoothed Aggregation for Second and Fourth Order Elliptic Problems

Multigrid methods are very efficient iterative solvers for system of algebraic equations arising from finite element and finite difference discretization of elliptic boundary value problems. The main principle of multigrid methods is to complement the local exchange of information in point-wise iterative methods by a global one utilizing several related systems, called coarse levels, with a smaller number of variables. The coarse levels are often obtained as a hierarchy of discretizations with different characteristic meshsizes, but this requires that the discretization is controlled by the iterative method. To solve linear systems produced by existing finite element software, one needs to create an artificial hierarchy of coarse problems. The principal issue is then to obtain computational complexity and approximation properties similar to those for nested meshes, using only information in the matrix of the system and as little extra information as possible. Such algebraic multigrid method that uses the system matrix only was developed by Ruge. The prolongations were based on the matrix of the system by partial solution from given values at selected coarse points. The coarse grid points were selected so that each point would be interpolated to via so-called strong connections. Our approach is based on smoothed aggregation introduced recently by Vanek. First the set of nodes is decomposed into small mutually disjoint subsets. A tentative piecewise constant interpolation (in the discrete sense) is then defined on those subsets as piecewise constant for second order problems, and piecewise linear for fourth order problems. The prolongation operator is then obtained by smoothing the output of the tentative prolongation and coarse level operators are defined variationally.

Vanek, Petr↗

A Users Guide for the NASA ANOPP Propeller Analysis System

The purpose of this report is to document improvements to the Propeller Analysis System of the Aircraft Noise Prediction Program (PAS-ANOPP) and to serve as a users guide. An overview of the functional modules and modifications made to the Propeller ANOPP system are described. Propeller noise predictions are made by executing a sequence of functional modules through the use of ANOPP control statements. The most commonly used ANOPP control statements are discussed with detailed examples demonstrating the use of each control statement. Originally, the Propeller Analysis System included the angle-of-attack only in the performance module. Recently, modifications have been made to also include angle-of-attack in the noise prediction module. A brief description of PAS prediction capabilities is presented which illustrate the input requirements necessary to run the code by way of ten templates. The purpose of the templates are to provide PAS users with complete examples which can be modified to serve their particular purposes. The examples include the use of different approximations in the computation of the noise and the effects of synchrophasing. Since modifications have been made to the original PAS-ANOPP, comparisons of the modified ANOPP and wind tunnel data are also included. Two appendices are attached at the end of this report which provide useful reference material. One appendix summarizes the PAS functional modules while the second provides a detailed discussion of the TABLE control statement.

Nguygen, L. Cathy↗

Heat Transfer Computations of Internal Duct Flows With Combined Hydraulic and Thermal Developing Length

This study investigated the Navier-Stokes computations of the surface heat transfer coefficients of a transition duct flow. A transition duct from an axisymmetric cross section to a non-axisymmetric cross section, is usually used to connect the turbine exit to the nozzle. As the gas turbine inlet temperature increases, the transition duct is subjected to the high temperature at the gas turbine exit. The transition duct flow has combined development of hydraulic and thermal entry length. The design of the transition duct required accurate surface heat transfer coefficients. The Navier-Stokes computational method could be used to predict the surface heat transfer coefficients of a transition duct flow. The Proteus three-dimensional Navier-Stokes numerical computational code was used in this study. The code was first studied for the computations of the turbulent developing flow properties within a circular duct and a square duct. The code was then used to compute the turbulent flow properties of a transition duct flow. The computational results of the surface pressure, the skin friction factor, and the surface heat transfer coefficient were described and compared with their values obtained from theoretical analyses or experiments. The comparison showed that the Navier-Stokes computation could predict approximately the surface heat transfer coefficients of a transition duct flow.

Wang, C. R.↗

Nano-Plasticity of Single-Wall Carbon Nanotubes Under Uniaxial Compression

Nano-plasticity of thin single-wall carbon nanotubes under uniaxial compression is investigated through generalized tight-binding molecular dynamics (GTBMD) and ab-initio electronic structure methods. A novel mechanism of nano-plasticity of carbon nanotubes under uniaxial compression is observed in which bonding geometry collapses from a graphitic (sp(sup 2)) to a localized diamond like (sp(sup 3)) reconstruction. The computed critical stress (approximately equals 153 G Pa) and the shape of the resulting plastic deformation is in good agreement with recent experimental observation of collapse and fracture of compressed carbon nanotubes in polymer composites.

Srivastava, Deepak↗

RF Device for Acquiring Images of the Human Body

A safe, non-invasive method for forming images through clothing of large groups of people, in order to search for concealed weapons either made of metal or not, has been developed. A millimeter wavelength scanner designed in a unique, ring-shaped configuration can obtain a full 360 image of the body with a resolution of less than a millimeter in only a few seconds. Millimeter waves readily penetrate normal clothing, but are highly reflected by the human body and concealed objects. Millimeter wave signals are nonionizing and are harmless to human tissues when used at low power levels. The imager (see figure) consists of a thin base that supports a small-diameter vertical post about 7 ft (=2.13 m) tall. Attached to the post is a square-shaped ring 2 in. (=5 cm) wide and 3 ft (=91 cm) on a side. The ring is oriented horizontally, and is supported halfway along one side by a connection to a linear bearing on the vertical post. A planar RF circuit board is mounted to the inside of each side of the ring. Each circuit board contains an array of 30 receivers, one transmitter, and digitization electronics. Each array element has a printed-circuit patch antenna coupled to a pair of mixers by a 90 coupler. The mixers receive a reference local oscillator signal to a subharmonic of the transmitter frequency. A single local oscillator line feeds all 30 receivers on the board. The resulting MHz IF signals are amplified and carried to the edge of the board where they are demodulated and digitized. The transmitted signal is derived from the local oscillator at a frequency offset determined by a crystal oscillator. One antenna centrally located on each side of the square ring provides the source illumination power. The total transmitted power is less than 100 mW, resulting in an exposure level that is completely safe to humans. The output signals from all four circuit boards are fed via serial connection to a data processing computer. The computer processes the approximately 1-MB data set into a three-dimensional image in a matter of seconds. The innovation is to configure the receiver array in a ring topology surrounding the scanned object. The ring is then scanned vertically to cover the necessary two-dimensional surface. This fabrication of the ring is made possible by using planar antenna and circuit technology. A planar circuit board serves as a medium for both antennas and signal processing components. Using this technique, parts counts are kept low, and the cost per element is a small fraction of a waveguide-based system.

Gaier, Todd C.↗

Optimizing Integrated Terminal Airspace Operations Under Uncertainty

In the terminal airspace, integrated departures and arrivals have the potential to increase operations efficiency. Recent research has developed geneticalgorithm- based schedulers for integrated arrival and departure operations under uncertainty. This paper presents an alternate method using a machine jobshop scheduling formulation to model the integrated airspace operations. A multistage stochastic programming approach is chosen to formulate the problem and candidate solutions are obtained by solving sample average approximation problems with finite sample size. Because approximate solutions are computed, the proposed algorithm incorporates the computation of statistical bounds to estimate the optimality of the candidate solutions. A proof-ofconcept study is conducted on a baseline implementation of a simple problem considering a fleet mix of 14 aircraft evolving in a model of the Los Angeles terminal airspace. A more thorough statistical analysis is also performed to evaluate the impact of the number of scenarios considered in the sampled problem. To handle extensive sampling computations, a multithreading technique is introduced.

air traffic optimization↗

Electromagnetic Scattering by Discrete Random Media. IV: Coherent Backscattering

The problem of backscattering of light by a discrete random medium illuminated by an obliquely incident plane electromagnetic wave is considered.The analysis is performed in a linear-polarization basis and includes a complete derivation of the cross reflection matrix for a layer with densely and sparsely distributed particles, the design of an approximate method for computing the ladder and cross reflection matrices in the case of a semi-infinite medium with a sparse distribution of particles, the derivation of the relations between the elements of the ladder and cross reflection matrices in the exact backscattering direction for dense and sparse media, and the development of practical algorithms for solving the underlying integral equations by the method of Picard iterations and the discrete ordinate method. Simulation results for particles with large size parameters are also presented.

Adrian Doicu↗

Optimal Control Prediction Method for Control Allocation

This paper proposes a novel prediction method for online optimal control allocation that extends the volume of moments achievable with the Moore-Penrose generalized inverse to the entire Attainable Moment Set. This method formulates the control allocation problem using selected basis vectors and associated gains which reduces the optimization problem dimensions and provides physical insight into the resulting optimal solutions. The proposed algorithm finds the entire family of unique optimal control solutions along the desired moment vector from the origin to the boundary of the Attainable Moment Set. Numerical results for the Moore-Penrose prediction method show that the unique minimal controls obtained yield the desired moment with near machine precision accuracy while maintaining control effectors within specified position limits. This method has been fully validated against the unique solution obtained on the boundary of the Attainable Moment Set using the Durham Direct Allocation method. Minimal control solutions obtained for moments in the interior of the Attainable Moment Set, similarly yield the desired moment to near machine precision while providing control solutions that are smaller (i.e. 2-norm) than solutions found with traditional control allocation algorithms (e.g. interior point methods) applied to the minimal control problem. Numerical simulations using a Matlab® autocoded executable (MEX) for the representative real world problem of 3-moments with 20 individual control effectors and prescribed control position limits show a mean computation speed of approximately 125 Hz which is sufficient to enable real-time flight allocation.

Acheson, Michael J.↗

Towards Formalization of Advanced Linear Algebra with Applications to Dynamical Systems using PVS

Linear Algebra is essential for numerous aerospace problems of interest. Formal reasoning about hybrid systems that contain variables modeled by differential equations rely on concepts from Linear Algebra such as eigenvalues, matrix decompositions, and matrix valued functions. For example, the long-term dynamics of a system of differential equations depend on the stability/instability of its equilibrium points, which often reduces to an eigenvalue problem. This talk will embark on a quest to formalize theorems and results about eigenvalues and eigenvectors using PVS. We shall start our journey with 2 x 2 complex matrices, where we will apply our PVS code to a simple example of a dynamical system. Since it can be difficult or impossible to give simple expressions of eigenvalues for larger matrices (i.e. 5 x 5 or higher), we then move towards specifying the power method for verified computation of eigenvalue approximations in PVS. This effort requires development of multivariate complex arithmetic. At the end of the day, having such additions to the PVS NASA libraries will help move towards the use of formal methods to verify concepts of control theory and system level verification.

Linear Algebra↗