Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “approximation algorithm”

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,135 records · Page 63

An efficient method for estimating neighboring steady-state numerical solutions to the Euler equations

The paper concentrates on a linear approximation method for predicting the changes occurring in steady-state numerical solutions of the Euler equations as a consequence of small changes in the independent variables which control the problem. The importance of proper boundary-condition treatment and other issues concerning the problem are covered along with the importance of proper algorithm selection for a fully supersonic inviscid flow. The method is applied to a subsonic nozzle involving variation of the pressure on the outflow boundary and to a supersonic inlet involving variation of the inflow Mach number. In the subsonic test case, the comparisons between the predicted and conventional numerical solutions are shown to be good, while in the supersonic test case, the agreement between the approximation method and conventional numerical solution starts out well but rapidly degenerates at some point in the flowfield as the perturbation of the boundary conditions is increased.

Taylor, Arthur C., III↗

Applications of an adaptive unstructured solution algorithm to the analysis of high speed flows

An upwind cell-centered scheme for the solution of steady laminar viscous high-speed flows is implemented on unstructured two-dimensional meshes. The first-order implementation employs Roe's (1981) approximate Riemann solver, and a higher-order extension is produced by using linear reconstruction with limiting. The procedure is applied to the solution of inviscid subsonic flow over an airfoil, inviscid supersonic flow past a cylinder, and viscous hypersonic flow past a double ellipse. A detailed study is then made of a hypersonic laminar viscous flow on a 24-deg compression corner. It is shown that good agreement is achieved with previous predictions using finite-difference and finite-volume schemes. However, these predictions do not agree with experimental observations. With refinement of the structured grid at the leading edge, good agreement with experimental observations for the distributions of wall pressure, heating rate and skin friction is obtained.

Thareja, R. R.↗

A Computer Program for the Calculation of Three-Dimensional Transonic Nacelle/Inlet Flowfields

A highly efficient computer analysis was developed for predicting transonic nacelle/inlet flowfields. This algorithm can compute the three dimensional transonic flowfield about axisymmetric (or asymmetric) nacelle/inlet configurations at zero or nonzero incidence. The flowfield is determined by solving the full-potential equation in conservative form on a body-fitted curvilinear computational mesh. The difference equations are solved using the AF2 approximate factorization scheme. This report presents a discussion of the computational methods used to both generate the body-fitted curvilinear mesh and to obtain the inviscid flow solution. Computed results and correlations with existing methods and experiment are presented. Also presented are discussions on the organization of the grid generation (NGRIDA) computer program and the flow solution (NACELLE) computer program, descriptions of the respective subroutines, definitions of the required input parameters for both algorithms, a brief discussion on interpretation of the output, and sample cases to illustrate application of the analysis.

Vadyak, J.↗

Linear state feedback, quadratic weights, and closed loop eigenstructures

Equations are derived for the angles of general multivariable root loci and linear quadratic optimal root loci, including angles of departure and approach. The generalized eigenvalue problem is used to compute angles of approach. Equations are also derived to find the sensitivity of closed loop eigenvalue and the directional derivatives of closed loop eigenvectors. An equivalence class of quadratic weights that produce the same asymptotic eigenstructure is defined, a canonical element is defined, and an algorithm to find it is given. The behavior of the optimal root locus in the nonasymptotic region is shown to be different for quadratic weights with the same asymptotic properties. An algorithm is presented that can be used to select a feedback gain matrix for the linear state feedback problem which produces a specified asymptotic eigenstructure. Another algorithm is given to compute the asymptotic eigenstructure properties inherent in a given set of quadratic weights. Finally, it is shown that optimal root loci for nongeneric problems can be approximated by generic ones in the nonasymptotic region.

Thompson, P. M.↗

Splitting of inviscid fluxes for real gases

Flux-vector and flux-difference splittings for the inviscid terms of the compressible flow equations are derived under the assumption of a general equation of state for a real gas in equilibrium. No necessary assumptions, approximations or auxiliary quantities are introduced. The formulas derived include several particular cases known for ideal gases and readily apply to curvilinear coordinates. Applications of the formulas in a TVD algorithm to one-dimensional shock-tube and nozzle problems show their quality and robustness.

Liou, Meng-Sing↗

Splitting of inviscid fluxes for real gases

Flux-vector and flux-difference splittings for the inviscid terms of the compressible flow equations are derived under the assumption of a general equation of state for a real gas in equilibrium. No necessary assumptions, approximations or auxiliary quantities are introduced. The formulas derived include several particular cases known for ideal gases and readily apply to curvilinear coordinates. Applications of the formulas in a TVD algorithm to one-dimensional shock-tube and nozzle problems show their quality and robustness.

Liou, Meng-Sing↗

Splitting of inviscid fluxes for real gases

Flux-vector and flux-difference splittings for the inviscid terms of the compressible flow equations are derived under the assumption of a general equation of state for a real gas in equilibrium. No necessary assumptions, approximations for auxiliary quantities are introduced. The formulas derived include several particular cases known for ideal gases and readily apply to curvilinear coordinates. Applications of the formulas in a TVD algorithm to one-dimensional shock-tube and nozzle problems show their quality and robustness.

Liou, Meng-Sing↗

The IRAS Galaxy Atlas (IGA)

In 1993 we proposed a project to NASA having the goal of producing a new infrared map of our Galaxy. In particular, we proposed to reprocess the IRAS data taken in the early 1980's using modern image processing algorithms and the large Intel parallel computers of the Center for Advanced Computing Research, (at that time called the Caltech Concurrent Supercomputing Facilities - CCSF). The rationale was simple: what took approximately 100 days on a typical workstation would take less than a day on the multi-processor parallel computers, thus making a high-resolution infrared atlas of the Galaxy feasible.

Prince, Thomas A.↗

NASA Tech Briefs, December 2009

Topics include: A Deep Space Network Portable Radio Science Receiver; Detecting Phase Boundaries in Hard-Sphere Suspensions; Low-Complexity Lossless and Near-Lossless Data Compression Technique for Multispectral Imagery; Very-Long-Distance Remote Hearing and Vibrometry; Using GPS to Detect Imminent Tsunamis; Stream Flow Prediction by Remote Sensing and Genetic Programming; Pilotless Frame Synchronization Using LDPC Code Constraints; Radiometer on a Chip; Measuring Luminescence Lifetime With Help of a DSP; Modulation Based on Probability Density Functions; Ku Telemetry Modulator for Suborbital Vehicles; Photonic Links for High-Performance Arraying of Antennas; Reconfigurable, Bi-Directional Flexfet Level Shifter for Low-Power, Rad-Hard Integration; Hardware-Efficient Monitoring of I/O Signals; Video System for Viewing From a Remote or Windowless Cockpit; Spacesuit Data Display and Management System; IEEE 1394 Hub With Fault Containment; Compact, Miniature MMIC Receiver Modules for an MMIC Array Spectrograph; Waveguide Transition for Submillimeter-Wave MMICs; Magnetic-Field-Tunable Superconducting Rectifier; Bonded Invar Clip Removal Using Foil Heaters; Fabricating Radial Groove Gratings Using Projection Photolithography; Gratings Fabricated on Flat Surfaces and Reproduced on Non-Flat Substrates; Method for Measuring the Volume-Scattering Function of Water; Method of Heating a Foam-Based Catalyst Bed; Small Deflection Energy Analyzer for Energy and Angular Distributions; Polymeric Bladder for Storing Liquid Oxygen; Pyrotechnic Simulator/Stray-Voltage Detector; Inventions Utilizing Microfluidics and Colloidal Particles; RuO2 Thermometer for Ultra-Low Temperatures; Ultra-Compact, High-Resolution LADAR System for 3D Imaging; Dual-Channel Multi-Purpose Telescope; Objective Lens Optimized for Wavefront Delivery, Pupil Imaging, and Pupil Ghosting; CMOS Camera Array With Onboard Memory; Quickly Approximating the Distance Between Two Objects; Processing Images of Craters for Spacecraft Navigation; Adaptive Morphological Feature-Based Object Classifier for a Color Imaging System; Rover Slip Validation and Prediction Algorithm; Safety and Quality Training Simulator; Supply-Chain Optimization Template; Algorithm for Computing Particle/Surface Interactions; Cryogenic Pupil Alignment Test Architecture for Aberrated Pupil Images; and Thermal Transport Model for Heat Sink Design.

Source record↗

Analytical estimates of the PP-algorithm at low number of Doppler periods per pulse length

When discussing the Doppler velocity estimators, it is of significant interest to analyze their behavior at a low number of Doppler periods n(sub D) = 2v(sub r)t(sub s)/lambda is approximately equal to 1 within the resolution cell t(sub s) (v(sub 4) is the radial velocity, lambda is the wavelength). Obviously, at n(sub D) is approximately less than 1 the velocity error is essentially increased. The problem of low n(sub D) arises in the planetary boundary layer (PBL), where higher resolutions are usually required but the signal-to-noise ratio (SNR) is relatively high. In this work analytical expression for the relative root mean square (RMS) error of the PP Doppler estimator at low number of periods for a narrowband Doppler signal and arbitrary model of the noise correlation function is obtained. The results are correct at relatively high SNR. The analysis is supported by computer simulations at various SNR's.

Angelova, M. D.↗

Evidence for Reduced Species Star Formation Rates in the Centers of Massive Galaxies at zeta = 4

We perform the first spatially-resolved stellar population study of galaxies in the early universe z equals 3.5 -6.5, utilizing the Hubble Space Telescope Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) imaging dataset over the GOODS-S field. We select a sample of 418 bright and extended galaxies at z less than or approximately equal to 3.5-6.5 from a parent sample of approximately 8000 photometric-redshift selected galaxies from Finkelstein et al. We first examine galaxies at 3.5 less than or equal to z less than or approximately equal to 4.0 using additional deep K-band survey data from the HAWK-I UDS and GOODS Survey (HUGS) which covers the 4000 Angstrom break at these redshifts. We measure the stellar mass, star formation rate, and dust extinction for galaxy inner and outer regions via spatially-resolved spectral energy distribution fitting based on a Markov Chain Monte Carlo algorithm. By comparing specific star formation rates (sSFRs) between inner and outer parts of the galaxies we find that the majority of galaxies with the high central mass densities show evidence for a preferentially lower sSFR in their centers than in their outer regions, indicative of reduced sSFRs in their central regions. We also study galaxies at z approximately equal to 5 and 6 (here limited to high spatial resolution in the rest-frame ultraviolet only), finding that they show sSFRs which are generally independent of radial distance from the center of the galaxies. This indicates that stars are formed uniformly at all radii in massive galaxies at z approximately equal to 5-6, contrary tomassive galaxies at z. less than approximately equal to 4.

early universe – galaxies: bulges – galaxies: ↗

Formulation and implementation of a practical algorithm for non-stationary adaptive state estimation

Background information on the Kalman filter is given first. A discussion of the filter parameters and their a priori determination follows. The discussion points out the need for adaptive determination of the process noise statistics. The filter innovations are presented as a means for developing the adaptive criteria. The criteria center around the estimation of the true mean and covariance of the filter innovations. A method for the numerical approximation of the mean and covariance of a locally stationary random process is presented. The definition of a local stationarity is presented. Local stationarity allows for the separation of the process statistics into a stationary component and a time-varying component. The separation method is discussed. A method for estimating the stationary and time-varying components is presented. As an example of its application to real problems, the algorithm is applied to the problem to the problem of reentry trajectory estimation for the Space Shuttle. Both the adaptive algorithm and the steady-state Kalman filter are applied to the problem. The results of the reconstructions are presented. The adaptive algorithm exhibits superior performance.

Whitemore, S. A.↗

A novel method of calculating far-field patterns of large aperture antennas

A method is described for calculation of the radiation pattern of large aperture antennas. A piece-wise linear approximation of the aperture field using overlapping pyramidal basis functions allows the radiation pattern of an aperture antenna to be calculated as though it were a two-dimensional array. The calculation of radiation pattern data versus theta and phi, suitable for 3-D or contour plot algorithms, is achieved by locating the array in the yz-plane and performing a summation over the aperture field data sampled on a square grid. A FORTRAN subroutine is provided for performing radiation pattern calculations. Numerical results are included to demonstrate the accuracy and convergence of the method. These numerical results indicate that typical accuracies of + or - 0.1 dB for Directivity, + or - dB for the 1st Sidelobe Level, and + - 2dB for the 2nd Sidelobe Level can be obtained with an aperture grid of 45x45 points and requires approximately 0.02 seconds CPU time per far-field data point on a VAX 11/750 with a floating point accelerator.

Bailey, M. C.↗

Estimating Planetary Boundary Layer Heights from NOAA Profiler Network Wind Profiler Data

An algorithm was developed to estimate planetary boundary layer (PBL) heights from hourly archived wind profiler data from the NOAA Profiler Network (NPN) sites located throughout the central United States. Unlike previous studies, the present algorithm has been applied to a long record of publicly available wind profiler signal backscatter data. Under clear conditions, summertime averaged hourly time series of PBL heights compare well with Richardson-number based estimates at the few NPN stations with hourly temperature measurements. Comparisons with clear sky reanalysis based estimates show that the wind profiler PBL heights are lower by approximately 250-500 m. The geographical distribution of daily maximum PBL heights corresponds well with the expected distribution based on patterns of surface temperature and soil moisture. Wind profiler PBL heights were also estimated under mostly cloudy conditions, and are generally higher than both the Richardson number based and reanalysis PBL heights, resulting in a smaller clear-cloudy condition difference. The algorithm presented here was shown to provide a reliable summertime climatology of daytime hourly PBL heights throughout the central United States.

Wind Profiler↗

Automated design using numerical optimization

Numerical optimization concepts are described with limited technical detail. The purpose is to provide the nonspecialist with sufficient information to judge the applicability of these methods to his particular design problem. The concepts are first described in physical terms to give a basic understanding of the iterative procedure employed by these methods. Next, the typical engineering task is presented and converted to a form ammenable to solution by numerical optimization. Basic algorithms for solving this problem are identified. Numerous applications are referenced, emphasizing the structural design discipline. The state of the art allows for the routine solution of nonlinear design problems of approximately 20 independent variables subject to 100 or more constraints. In many applications, much larger design problems may be solved. Selected references are provided which describe the methods and applications in more detail.

Vanderplaats, G. N.↗

Structural synthesis: Precursor and catalyst

More than twenty five years have elapsed since it was recognized that a rather general class of structural design optimization tasks could be properly posed as an inequality constrained minimization problem. It is suggested that, independent of primary discipline area, it will be useful to think about: (1) posing design problems in terms of an objective function and inequality constraints; (2) generating design oriented approximate analysis methods (giving special attention to behavior sensitivity analysis); (3) distinguishing between decisions that lead to an analysis model and those that lead to a design model; (4) finding ways to generate a sequence of approximate design optimization problems that capture the essential characteristics of the primary problem, while still having an explicit algebraic form that is matched to one or more of the established optimization algorithms; (5) examining the potential of optimum design sensitivity analysis to facilitate quantitative trade-off studies as well as participation in multilevel design activities. It should be kept in mind that multilevel methods are inherently well suited to a parallel mode of operation in computer terms or to a division of labor between task groups in organizational terms. Based on structural experience with multilevel methods general guidelines are suggested.

Schmit, L. A.↗

Inverse solution of Kepler's equation for hyperbolic orbits

An algorithm is presented for efficient inverse solution of Kepler's equation for hyperbolic orbits. It is shown that an expansion of Barker's equation into a bicubic polynomial provides a good approximation to obtain accurate starting values for rapid numerical solution of Kepler's equation. In the approximate equation a cubic in normalized elapsed flight time from pericenter is set equal to a cubic in a function S of eccentricity and true anomaly. The initial estimate of S to use in an iteration formula is obtained by evaluating the cubic in normalized flight time and finding in most cases the single real root of the other cubic. This initial estimate has an accuracy corresponding to values of true anomaly in error by less than 0.5 degrees generally.

Boltz, Frederick W.↗

Collision-induced gas phase dissociation rates

The Landau-Zener theory of reactive cross sections was applied to diatomic molecules dissociating from a ladder of vibrational states. The result predicts a dissociation rate that is quite well duplicated by an Arrhenius function having a preexponential temperature dependence of about T(sub -1/2), at least for inert collision partners. This relation fits experimental data reasonably well. The theory is then used to calculate the effect of vibrational nonequilibrium on dissociation rate. For Morse oscillators, the results are about the same as given by Hammerling, Kivel, and Teare in their analytic approximation for harmonic oscillators, though at very high temperature a correction for the partition function limit is included. The empirical correction for vibration nonequilibrium proposed by Park, which is a convenient algorithm for CFD calculations, is modified to prevent a drastic underestimation of dissociation rates that occurs with this method when vibrational temperature is much smaller than the kinetic temperature of the gas.

Hansen, C. Frederick↗