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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

NASA NTRS1979

Numerical spatial marching techniques in duct acoustics

Direct calculation of the internal structure of a ducted noise source from farfield pressure measurements is regarded as an initial value problem, where the pressure and pressure gradient (farfield impedance) are assumed to be known along a line in the farfield. If pressure and impedance are known at the boundary of the farfield, the pressure can be uniquely determined in the vicinity of the inlet and inside the inlet ducting. A marching procedure is developed which, with this information obtained from measurements, enables a description of a ducted noise source. The technique uses a finite difference representation of the homogeneous Helmholtz equation.

Baumeister, K. J.↗

NASA NTRS1980

Numerical grid generation techniques

The state of the art in topology and flow geometry is presented. Solution techniques for partial differential equations are reviewed and included developments in coordinate transformations, conformal mapping, and invariant imbeddings. Applications of these techniques in fluid mechanics, flow geometry, boundary value problems, and fluidics are presented.

Source record↗

NASA NTRS1986

On substructuring algorithms and solution techniques for the numerical approximation of partial differential equations

Substructuring methods are in common use in mechanics problems where typically the associated linear systems of algebraic equations are positive definite. Here these methods are extended to problems which lead to nonpositive definite, nonsymmetric matrices. The extension is based on an algorithm which carries out the block Gauss elimination procedure without the need for interchanges even when a pivot matrix is singular. Examples are provided wherein the method is used in connection with finite element solutions of the stationary Stokes equations and the Helmholtz equation, and dual methods for second-order elliptic equations.

Gunzburger, M. D.↗

NASA NTRS1989

Aspects of efficient and reliable multibody system simulation

Multibody system equations can be generated in various forms. All of these may be interpreted as results of two basic approaches, the augmentation- and the elimination-method. The former method yields the descriptor form of the system motion, a set of differential-algebraic equations (DAE), and the latter the state space representation, a minimal set of ordinary differential equations (ODE). Both of these methods are surveyed. Particular emphasis is on the discussion of recursive computational schemes, generating the equations of motion with a number of operations, which is proportional to the number N of system bodies (O(N)-formulations). For simulation purposes one would like to create that set of system equations, which can be generated most efficiently and for which the most efficient and reliable solution techniques are available. Numerical solution techniques for ODE have been studied in great detail and they are well-developed. By contrast, DAE have not been investigated for such a long time. In view of new developments in the latter field the generation of all the equations required for an efficient and reliable solution of DAE describing multibody system motion is discussed. These methods, i.e., an O(N)-formulation and new techniques for solving DAE, are implemented in the SIMPACK code. Its capabilities are illustrated by simulation of multibody robot models.

Schwertassek, R.↗

NASA NTRS1996

Dynamics of Nuclear Regions of Galaxies

Current research carried out with the help of the ASEE-NASA Summer Faculty Program, at NASA-Ames, is concentrated on the dynamics of nuclear regions of galaxies. From a dynamical point of view a galaxy is a collection of around 10(sup 11) stars like our Sun, each of which moves in the summed gravitational field of all the remaining stars. Thus galaxy dynamics becomes a self-consistent n-body problem with forces given by Newtonian gravitation. Strong nonlinearity in the gravitational force and the inherent nonlinearity of self-consistent problems both argue for a numerical approach. The technique of numerical experiments consis of constructing an environment in the computer that is as close as possible to the physical conditions in a real galaxy and then carrying out experiments much like laboratory experiments in physics or engineering, in this environment. Computationally, an experiment is an initial value problem, and a good deal of thought and effort goes into the design of the starting conditions that serve as initial values. Experiments are run at Ames because all the 'equipment' is in place-the programs, the necessary computational power, and good facilities for post-run analysis. Our goal for this research program is to study the nuclear regions in detail and this means replacing most of the galaxy by a suitable boundary condition to allow the full capability of numerical experiments to be brought to bear on a small region perhaps 1/1000 of the linear dimensions of an entire galaxy. This is an extremely delicate numerical problem, one in which some small feature overlook, can easily lead to a collapse or blow-up of the entire system. All particles attract each other in gravitational problems, and the 1/r(sup 2) force is: (1) nonlinear; (2) strong at short range; (3) long-range, and (4) unscreened at any distance.

Miller, Richard H.↗

NASA NTRS2003

Analytic Method for Computing Instrument Pointing Jitter

A new method of calculating the root-mean-square (rms) pointing jitter of a scientific instrument (e.g., a camera, radar antenna, or telescope) is introduced based on a state-space concept. In comparison with the prior method of calculating the rms pointing jitter, the present method involves significantly less computation. The rms pointing jitter of an instrument (the square root of the jitter variance shown in the figure) is an important physical quantity which impacts the design of the instrument, its actuators, controls, sensory components, and sensor- output-sampling circuitry. Using the Sirlin, San Martin, and Lucke definition of pointing jitter, the prior method of computing the rms pointing jitter involves a frequency-domain integral of a rational polynomial multiplied by a transcendental weighting function, necessitating the use of numerical-integration techniques. In practice, numerical integration complicates the problem of calculating the rms pointing error. In contrast, the state-space method provides exact analytic expressions that can be evaluated without numerical integration.

Bayard, David↗

NASA NTRS2003

A New Objective Technique for Verifying Mesoscale Numerical Weather Prediction Models

This report presents a new objective technique to verify predictions of the sea-breeze phenomenon over east-central Florida by the Regional Atmospheric Modeling System (RAMS) mesoscale numerical weather prediction (NWP) model. The Contour Error Map (CEM) technique identifies sea-breeze transition times in objectively-analyzed grids of observed and forecast wind, verifies the forecast sea-breeze transition times against the observed times, and computes the mean post-sea breeze wind direction and speed to compare the observed and forecast winds behind the sea-breeze front. The CEM technique is superior to traditional objective verification techniques and previously-used subjective verification methodologies because: It is automated, requiring little manual intervention, It accounts for both spatial and temporal scales and variations, It accurately identifies and verifies the sea-breeze transition times, and It provides verification contour maps and simple statistical parameters for easy interpretation. The CEM uses a parallel lowpass boxcar filter and a high-order bandpass filter to identify the sea-breeze transition times in the observed and model grid points. Once the transition times are identified, CEM fits a Gaussian histogram function to the actual histogram of transition time differences between the model and observations. The fitted parameters of the Gaussian function subsequently explain the timing bias and variance of the timing differences across the valid comparison domain. Once the transition times are all identified at each grid point, the CEM computes the mean wind direction and speed during the remainder of the day for all times and grid points after the sea-breeze transition time. The CEM technique performed quite well when compared to independent meteorological assessments of the sea-breeze transition times and results from a previously published subjective evaluation. The algorithm correctly identified a forecast or observed sea-breeze occurrence or absence 93% of the time during the two- month evaluation period from July and August 2000. Nearly all failures in CEM were the result of complex precipitation features (observed or forecast) that contaminated the wind field, resulting in a false identification of a sea-breeze transition. A qualitative comparison between the CEM timing errors and the subjectively determined observed and forecast transition times indicate that the algorithm performed very well overall. Most discrepancies between the CEM results and the subjective analysis were again caused by observed or forecast areas of precipitation that led to complex wind patterns. The CEM also failed on a day when the observed sea- breeze transition affected only a very small portion of the verification domain. Based on the results of CEM, the RAMS tended to predict the onset and movement of the sea-breeze transition too early and/or quickly. The domain-wide timing biases provided by CEM indicated an early bias on 30 out of 37 days when both an observed and forecast sea breeze occurred over the same portions of the analysis domain. These results are consistent with previous subjective verifications of the RAMS sea breeze predictions. A comparison of the mean post-sea breeze winds indicate that RAMS has a positive wind-speed bias for .all days, which is also consistent with the early bias in the sea-breeze transition time since the higher wind speeds resulted in a faster inland penetration of the sea breeze compared to reality.

Case, Jonathan L.↗

NASA NTRS1988

The solution of the Elrod algorithm for a dynamically loaded journal bearing using multigrid techniques

A numerical solution to a theoretical model of vapor cavitation in a dynamically loaded journal bearing is developed utilizing a multigrid iteration technique. The method is compared with a noniterative approach in terms of computational time and accuracy. The computational model is based on the Elrod algorithm, a control volume approach to the Reynolds equation which mimics the Jakobsson-Floberg and Olsson cavitation theory. Besides accounting for a moving cavitation boundary and conservation of mass at the boundary, it also conserves mass within the cavitated region via a smeared mass or striated flow extending to both surfaces in the film gap. The mixed nature of the equations (parabolic in the full film zone and hyperbolic in the cavitated zone) coupled with the dynamic aspects of the problem create interesting difficulties for the present solution approach. Emphasis is placed on the methods found to eliminate solution instabilities. Excellent results are obtained for both accuracy and reduction of computational time.

Woods, Claudia M.↗

NASA NTRS1989

The solution of the Elrod algorithm for a dynamically loaded journal bearing using multigrid techniques

Woods, C. M.↗

NASA NTRS1973

A comparison of digital computer programs for the numerical solution of ordinary differential equations

Recently the determination of the best technique for numerically solving systems of ordinary differential equations on a digital computer has received much attention. The use of these formulas in conjunction with a stepsize control developed is explained, and one of the formulas is chosen for comparison with other integration techniques. This comparison of one of the best of Fehlberg's formulas with the different numerical techniques described in previous studies on a variety of test problems clearly shows the superiority of Fehlberg's formula. That is, on each of the test problems, the chosen Fehlberg formula is able to achieve a given accuracy in less computer time than any of the other techniques tested. Also, the computer program for the chosen Fehlberg formula is less complex and easier to use than the computer programs for most of the other techniques. To illustrate the use of the chosen Fehlberg formula, a computer listing of its application to several example problems is included along with a sample of the computer output from these applications.

Ingram, H. L.↗

NASA NTRS1988

Euler equation computations for the flow over a hovering helicopter rotor

A numerical solution technique is developed for computing the flow field around an isolated helicopter rotor in hover. The flow is governed by the compressible Euler equations which are integrated using a finite volume approach. The Euler equations are coupled to a free wake model of the rotary wing vortical wake. This wake model is incorporated into the finite volume solver using a prescribed flow, or perturbation, technique which eliminates the numerical diffusion of vorticity due to the artificial viscosity of the scheme. The work is divided into three major parts: (1) comparisons of Euler solutions to experimental data for the flow around isolated wings show good agreement with the surface pressures, but poor agreement with the vortical wake structure; (2) the perturbation method is developed and used to compute the interaction of a streamwise vortex with a semispan wing. The rapid diffusion of the vortex when only the basic Euler solver is used is illustrated, and excellent agreement with experimental section lift coefficients is demonstrated when using the perturbation approach; and (3) the free wake solution technique is described and the coupling of the wake to the Euler solver for an isolated rotor is presented. Comparisons with experimental blade load data for several cases show good agreement, with discrepancies largely attributable to the neglect of viscous effects. The computed wake geometries agree less well with experiment, the primary difference being that too rapid a wake contraction is predicted for all the cases.

Roberts, Thomas Wesley↗

NASA NTRS1985

Simulation analysis of a novel high efficiency silicon solar cell

It is recognized that crystalline silicon photovoltaic module efficiency of 15 percent or more is required for cost-effective photovoltaic energy utilization. This level of module efficiency requires large-area encapsulated production cell efficiencies in the range of 18 to 20 percent. Though the theoretical maximum of silicon solar cell efficiency for an idealized case is estimated to be around 30 percent, practical performance of cells to-date are considerably below this limit. This is understood to be largely a consequence of minority carrier losses in the bulk as well as at all surfaces including those under the metal contacts. In this paper a novel device design with special features to reduce bulk and surface recombination losses is evaluated using numerical analysis technique. Details of the numerical model, cell design, and analysis results are presented.

Mokashi, Anant R.↗

NASA NTRS1976

Automated optimization techniques for aircraft synthesis

Application of numerical optimization techniques to automated conceptual aircraft design is examined. These methods are shown to be a general and efficient way to obtain quantitative information for evaluating alternative new vehicle projects. Fully automated design is compared with traditional point design methods and time and resource requirements for automated design are given. The NASA Ames Research Center aircraft synthesis program (ACSYNT) is described with special attention to calculation of the weight of a vehicle to fly a specified mission. The ACSYNT procedures for automatically obtaining sensitivity of the design (aircraft weight, performance and cost) to various vehicle, mission, and material technology parameters are presented. Examples are used to demonstrate the efficient application of these techniques.

Vanderplaats, G. N.↗

NASA NTRS1987

Wave envelope technique for multimode wave guide problems

A numerical solution technique for two-dimensional EM waveguide problems is developed analytically and demonstrated. The frequency-domain waveguide problem is transformed to a time-dependent problem and treated following the limiting-amplitude principle; the wave-envelope technique is applied; and the resulting problem is solved numerically using a finite-difference scheme based on a variation of the finite-area method (Baumeister, 1977 and 1984). Results for several sample problems, obtained in 10 min or less of execution time on PCs, are presented in tables.

Hariharan, S. I.↗

NASA NTRS1984

Grid flexibility and patching techniques

The numerical determination of combustor flowfields is of great value to the combustor designer. An a priori knowledge of the flow behavior can speed the combustor design process and reduce the number of experimental test rigs required to arrive at an optimal design. Even 2-D steady incompressible isothermal flow predictions are of use; many codes of this kind are available, each employing different techniques to surmount the difficulties arising from the nonlinearity of the governing equations and from typically irregular combustor geometries. Mapping techniques (algebraic and elliptic PDE), and adaptive grid methods (both multi-grid and grid embedding) as applied to axisymmetric combustors are discussed.

Keith, T. G.↗

NASA NTRS1972

Computation of supersonic flow fields about bodies in coning motion using a shock-capturing finite-difference technique.

A numerical method for computing the nonlinear inviscid flow field surrounding a body performing coning motion is described. The method permits accurate computation of the aerodynamic moment due to one of the four motions characterizing an arbitrary nonplanar motion. Results of computations for a slender circular cone in coning motion are presented, and show good agreement with experiment for angles of attack up to twice the cone half angle. The computational results display significant departure of the side moment from the linear theory value with increasing angle of attack, but agree well with experimental measurements. This indicates that the initial nonlinear behavior of the aerodynamic moment is determined primarily by the inviscid flow.

Schiff, L. B.↗

NASA NTRS1972

Three dimensional nozzle-exhaust flow field analysis by a reference plane technique.

A numerical method based on reference plane characteristics has been developed for the calculation of highly complex supersonic nozzle-exhaust flow fields. The difference equations have been developed for three coordinate systems. Local reference plane orientations are employed using the three coordinate systems concurrently thus catering to a wide class of flow geometries. Discontinuities such as the underexpansion shock and contact surfaces are computed explicitly for nonuniform vehicle external flows. The nozzles considered may have irregular cross-sections with swept throats and may be stacked in modules using the vehicle undersurface for additional expansion. Results are presented for several nozzle configurations.

Dash, S. M.↗

NASA NTRS1976

Development of a computer code for calculating the steady super/hypersonic inviscid flow around real configurations. Volume 1: Computational technique

A numerical procedure has been developed to compute the inviscid super/hypersonic flow field about complex vehicle geometries accurately and efficiently. A second order accurate finite difference scheme is used to integrate the three dimensional Euler equations in regions of continuous flow, while all shock waves are computed as discontinuities via the Rankine Hugoniot jump conditions. Conformal mappings are used to develop a computational grid. The effects of blunt nose entropy layers are computed in detail. Real gas effects for equilibrium air are included using curve fits of Mollier charts. Typical calculated results for shuttle orbiter, hypersonic transport, and supersonic aircraft configurations are included to demonstrate the usefulness of this tool.

Marconi, F.↗