Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

The application of Encke's method to long arc orbit determination solutions

The Laser Geodynamics Satellite (LAGEOS) was launched on May 4, 1976 to provide geophysical measurements by means of laser ranging techniques. To date, over twelve years of laser range measurements have been collected from various tracking stations located around the world. Laser range measurements to LAGEOS have contributed to studies of earth rotation, plate tectonics, global baseline, and the gravity field as well as many other areas. Some of these studies are based upon the determination of a single, continuous orbit for LAGEOS for time spans on the order of several years. Current studies at the University of Texas Center for Space Research include the precision orbit determination of LAGEOS for arc lengths of up to 12.8 years which represents over 31,000 orbital revolutions. These long arc studies have led to the implementation of Encke's method to improve the convergence of the batch filter while reducing numerical integration errors. While the technique has been successfully applied to arc lengths of up to 12.8 years, the results presented focus on the solution of a six-year arc.

Lundberg, J. B.↗
NASA NTRS2013

Extended Analytic Device Optimization Employing Asymptotic Expansion

Analytic optimization of a thermoelectric junction often introduces several simplifying assumptionsincluding constant material properties, fixed known hot and cold shoe temperatures, and thermallyinsulated leg sides. In fact all of these simplifications will have an effect on device performance,ranging from negligible to significant depending on conditions. Numerical methods, such as FiniteElement Analysis or iterative techniques, are often used to perform more detailed analysis andaccount for these simplifications. While numerical methods may stand as a suitable solution scheme,they are weak in gaining physical understanding and only serve to optimize through iterativesearching techniques. Analytic and asymptotic expansion techniques can be used to solve thegoverning system of thermoelectric differential equations with fewer or less severe assumptionsthan the classic case. Analytic methods can provide meaningful closed form solutions and generatebetter physical understanding of the conditions for when simplifying assumptions may be valid.In obtaining the analytic solutions a set of dimensionless parameters, which characterize allthermoelectric couples, is formulated and provide the limiting cases for validating assumptions.Presentation includes optimization of both classic rectangular couples as well as practically andtheoretically interesting cylindrical couples using optimization parameters physically meaningful toa cylindrical couple. Solutions incorporate the physical behavior for i) thermal resistance of hot andcold shoes, ii) variable material properties with temperature, and iii) lateral heat transfer through legsides.

Thermoelectric Generators↗
NASA NTRS1974

A review of averaging techniques and their application to orbit determination systems

The theory of numerical averaging and analytical averaging will be reviewed and the application of these techniques to orbit and parameter estimation problems will be presented. Comparisons will be made between utilizing mean elements versus tracking data as the observation types. Results will be presented comparing the accuracy and efficiency of the combined orbit estimation and orbit prediction problem using averaged equations of motion, the Cowell equations of motion and the Brouwer general perturbation theory. The problem of converting the averaged element space back to osculating element space for orbit operations will also be discussed.

Velez, C. E.↗
NASA NTRS1976

Triangular covariance factorizations for

An improved computational form of the discrete Kalman filter is derived using an upper triangular factorization of the error covariance matrix. The covariance P is factored such that P = UDUT where U is unit upper triangular and D is diagonal. Recursions are developed for propagating the U-D covariance factors together with the corresponding state estimate. The resulting algorithm, referred to as the U-D filter, combines the superior numerical precision of square root filtering techniques with an efficiency comparable to that of Kalman's original formula. Moreover, this method is easily implemented and involves no more computer storage than the Kalman algorithm. These characteristics make the U-D method an attractive realtime filtering technique. A new covariance error analysis technique is obtained from an extension of the U-D filter equations. This evaluation method is flexible and efficient and may provide significantly improved numerical results. Cost comparisons show that for a large class of problems the U-D evaluation algorithm is noticeably less expensive than conventional error analysis methods.

Thornton, C. L.↗
NASA NTRS1980

Analytical and design techniques for drag reduction studies on wavy surfaces

Numerical models for two dimensional turbulent boundary layers over wavy surfaces were investigated. Computations for wavy wall boundary layers indicate possibilities of overall drag reduction in a parameter range of the geometry of the wall. The correction technique using integral methods for analyzing arbitrary surfaces was found to be unsuitable for some cases of interest in drag reduction; a Navier-Stokes solver for wavy walls was built to test these problems. Test results of the Navier-Stokes solver indicate that the solution techniques are accurate enough to handle complex geometries and steep variations in fluid properties.

Balasubramanian, R.↗
NASA NTRS2002

Analysis of High Order Difference Methods for Multiscale Complex Compressible Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes with incremental studies was initiated. Here we further refine the analysis on, and improve the understanding of the adaptive numerical dissipation control strategy. Basically, the development of these schemes focuses on high order nondissipative schemes and takes advantage of the progress that has been made for the last 30 years in numerical methods for conservation laws, such as techniques for imposing boundary conditions, techniques for stability at shock waves, and techniques for stable and accurate long-time integration. We concentrate on high order centered spatial discretizations and a fourth-order Runge-Kutta temporal discretizations as the base scheme. Near the bound-aries, the base scheme has stable boundary difference operators. To further enhance stability, the split form of the inviscid flux derivatives is frequently used for smooth flow problems. To enhance nonlinear stability, linear high order numerical dissipations are employed away from discontinuities, and nonlinear filters are employed after each time step in order to suppress spurious oscillations near discontinuities to minimize the smearing of turbulent fluctuations. Although these schemes are built from many components, each of which is well-known, it is not entirely obvious how the different components be best connected. For example, the nonlinear filter could instead have been built into the spatial discretization, so that it would have been activated at each stage in the Runge-Kutta time stepping. We could think of a mechanism that activates the split form of the equations only at some parts of the domain. Another issue is how to define good sensors for determining in which parts of the computational domain a certain feature should be filtered by the appropriate numerical dissipation. For the present study we employ a wavelet technique introduced in as sensors. Here, the method is briefly described with selected numerical experiments.

Sjoegreen, Bjoern↗
NASA NTRS1987

Geoid undulation computations at laser tracking stations

Geoid undulation computations were performed at 29 laser stations distributed around the world using a combination of terrestrial gravity data within a cap of radius 2 deg and a potential coefficient set up to 180 deg. The traditional methods of Stokes' and Meissl's modification together with the Molodenskii method and the modified Sjoberg method were applied. Performing numerical tests based on global error assumptions regarding the terrestrial data and the geopotential set it was concluded that the modified Sjoberg method is the most accurate and promising technique for geoid undulation computations. The numerical computations for the geoid undulations using all the four methods resulted in agreement with the ellipsoidal minus orthometric value of the undulations on the order of 60 cm or better for most of the laser stations in the eastern United States, Australia, Japan, Bermuda, and Europe. A systematic discrepancy of about 2 meters for most of the western United States stations was detected and verified by using two relatively independent data sets. For oceanic laser stations in the western Atlantic and Pacific oceans that have no terrestrial data available, the adjusted GEOS-3 and SEASAT altimeter data were used for the computation of the geoid undulation in a collocation method.

Despotakis, Vasilios K.↗
NASA NTRS1968

NASA Contributions to Development of Special-Purpose Thermocouples. A Survey

The thermocouple has been used for measuring temperatures for more than a century, but new materials, probe designs, and techniques are continually being developed. Numerous contributions have been made by the National Aeronautics and Space Administration and its contractors in the aerospace program. These contributions have been collected by Midwest Research Institute and reported in this publication to enable American industrial engineers to study them and adapt them to their own problem areas. Potential applications are suggested to stimulate ideas on how these contributions can be used.

Moeller, C. Eugene↗
NASA NTRS1971

Efficient digital comparison technique for logic circuits

Tolerance compare technique indicates discompare only when numerical difference value exceeds prescribed limit. Algorithm involving binary number properties is defined, in lieu of arithmetic operation which requires relatively complex circuitry. Extension of algorithm may be made to encompass tolerances other than one unit.

Mccarthy, C. E.↗
NASA NTRS1973

Numerical studies of unsteady two dimensional subsonic flows using the ICE method

A numerical program was developed to compute transient compressible and incompressible laminar flows in two dimensions with multicomponent mixing and chemical reaction. The algorithm used the Los Alamos Scientific Laboratory ICE (Implicit Continuous-Fluid Eulerian) method as its base. The program can compute both high and low speed compressible flows. The numerical program incorporating the stabilization techniques was quite successful in treating both old and new problems. Detailed calculations of coaxial flow very close to the entry plane were possible. The program treated complex flows such as the formation and downstream growth of a recirculation cell. An implicit solution of the species equation predicted mixing and reaction rates which compared favorably with the literature.

Wieber, P. R.↗
NASA NTRS1974

Solutions of the Navier-Stokes equations for vortex breakdown

Steady solutions of the Navier-Stokes equations, in terms of velocity and pressure, for breakdown in an unconfined viscous vortex are obtained numerically using the artificial compressibility technique of Chorin combined with an ADI finite-difference scheme. Axisymmetry is assumed and boundary conditions are carefully applied at the boundaries of a large finite region in an axial plane while resolution near the axis is maintained by a coordinate transformation. The solutions, which are obtained for Reynolds numbers up to 200 based on the free-stream axial velocity and a characteristic core radius, show that breakdown results from the diffusion and convection of vorticity away from the vortex core which, because of the strong coupling between the circumferential and axial velocity fields in strongly swirling flows, can lead to stagnation and reversal of the axial flow near the axis.

Grabowski, W. J.↗
NASA NTRS1974

A study to determine the usefulness of interval analysis in solving problems in celestial mechanics

This investigation was undertaken to determine the usefulness of interval analysis to numerical integration and matrix inversion techniques and to combine these results to determine the value of interval analysis in bounding computational errors in the two-body problem. Conclusions were that interval analysis may be worthwhile in certain small scale isolated problems, but its usefulness in any large scale problem is doubtful.

Walling, D.↗
NASA NTRS1977

Contouring randomly spaced data

Computer program using triangulation contouring technique contours data points too numerous to fit into rectangular grid. Using random access procedures, program can handle up to 56,000 data points and provides up to 20 contour intervals for multiple number of parameters.

Kibler, J. F.↗
NASA NTRS1975

Amplitude effects on the dynamic performance of a hydrostatic gas thrust bearing

The Reynolds' equation is applied to a strip gas thrust bearing to analyze amplitude disturbance effects on its dynamic performance. The Reynolds' equation is numerically approximated using finite difference techniques. The time dependent load carrying capacity is represented by a Fourier series up to and including the third harmonics. Design curves for the load capacity and the linear stiffness and damping are presented as a function of inlet location, restrictor coefficient, supply pressure, amplitude of oscillation, and squeeze number. For the range of amplitudes investigated the dimensionless load capacity, stiffness and damping does not exhibit an appreciable change in magnitude; thus, only one design curve is needed to represent each relationship. A design methodology is presented.

Stiffler, A. K.↗
NASA NTRS1976

Design analogy between optimal time-fuel and rate-ledge relay controllers

An empirically obtained reaction control jet relay control law with deadband is used as the basis for determining an equivalent weighted time-fuel optimal switching curve according to a least-squares criterion. The derived transformation from the empirical to the optimal law is found to be reversible and to yield a unique transformed control law. The proposed method provides a basis for determining the behavior of an easily implemented relay control law using well-known optimal control results, as well as determining the equivalent relay law corresponding to an analytically determined optimal control law. A numerical example illustrates the transformation technique and simulation results are presented to compare the two control laws.

White, L. R.↗
NASA NTRS1977

Development and application of a three dimensional numerical model for predicting pollutant and sediment transport using an Eulerian-Lagrangian marker particle technique

A computer coded Lagrangian marker particle in Eulerian finite difference cell solution to the three dimensional incompressible mass transport equation, Water Advective Particle in Cell Technique, WAPIC, was developed, verified against analytic solutions, and subsequently applied in the prediction of long term transport of a suspended sediment cloud resulting from an instantaneous dredge spoil release. Numerical results from WAPIC were verified against analytic solutions to the three dimensional incompressible mass transport equation for turbulent diffusion and advection of Gaussian dye releases in unbounded uniform and uniformly sheared uni-directional flow, and for steady-uniform plug channel flow. WAPIC was utilized to simulate an analytic solution for non-equilibrium sediment dropout from an initially vertically uniform particle distribution in one dimensional turbulent channel flow.

Pavish, D. L.↗
NASA NTRS1978

Inviscid, nonadiabatic flow fields over blunt, sonic corner bodies for outer planet entry conditions by a method of integral relations

An investigation has been made into the ability of a method of integral relations to calculate inviscid zero degree angle of attack, radiative heating distributions over blunt, sonic corner bodies for some representative outer planet entry conditions is investigated. Comparisons have been made with a more detailed numerical method, a time asymptotic technique, using the same equilibrium chemistry and radiation transport subroutines. An effort to produce a second order approximation (two-strip) method of integral relations code to aid in this investigation is also described and a modified two-strip routine is presented. Results indicate that the one-strip method of integral relations cannot be used to obtain accurate estimates of the radiative heating distribution because of its inability to resolve thermal gradients near the wall. The two-strip method can sometimes be used to improve these estimates; however, the two-strip method has only a small range of conditions over which it will yield significant improvement over the one-strip method.

Gnoffo, P. A.↗