Search NASA⌕ Search

SEARCH · Search NASA

Results for “Accuracy and stability”

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 289 records · Page 16

Study of second order upwind differencing in a recirculating flow

The accuracy and stability of the second order upwind differencing scheme was investigated. The solution algorithm employed is based on a coupled solution of the nonlinear finite difference equations by the multigrid technique. Calculations have been made of the driven cavity flow for several Reynolds numbers and finite difference grids. In comparison with the hybrid differencing, the second order upwind differencing is somewhat more accurate but it is not monotonically accurate with mesh refinement. Also, the convergence of the solution algorithm deteriorates with the use of the second order upwind differencing.

Vanka, S. P.↗

Analysis of some fully-discrete algorithms for the one-dimensional heat equation

The present investigation is concerned with a fully discrete accuracy and stability analysis of the one-dimensional heat equation, taking into account the evaluation of two-pass explicit schemes which simultaneously employ lumped and coupled capacity matrices. Schemes of the considered characteristics are not amenable to uncoupled semidiscrete and ordinary differential equation analyses. The obtained results illustrate that superior behavior may be achieved by schemes of the employed type when compared with the performance of the standard one-pass explicit schemes. The key idea in the considered approach is related to the utilization of a reduced-quadrature capacity matrix in the evaluation of the right-hand-side residual.

Hughes, T. J. R.↗

Variable-Conductance Heat Pipes

In response to need to accurately and efficiently predict performance of variable-conductance heat pipes (VCHP's) incorporated in spacecraft thermalcontrol systems, computer code VCHPDA developed to interact with thermal analyzer programs such as SINDA (Systems Improved Numerical Differencing Analyzer). Calculates length of gas-blocked region and vapor temperature in active portion. Advantages of VCHPDA over prior programs improved accuracy, unconditional stability, and increased efficiency of solution resulting from novel approach and use of state-of-the-art numerical techniques for solving VCHP mathematical model. Code valuable tool in design and evaluation of advanced thermal-control systems using variable-conductance heat pipes. Written in FORTRAN IV for use on CDC 600 computers.

Antoniuk, D.↗

On the solution of problems involving impact type loading

The response of beam and plate members to pulse and impact loading is investigated, and numerical methods are analyzed and compared. The equations of motion, the finite difference and finite element methods, the time integration/Runge-Kutta techniques, and material modeling are discussed in detail. It is found that both the finite difference and the finite element methods could be accurately employed to discretize the spatial variation in the displacements. Central differences or a fourth order Runge-Kutta algorithm could be used for the time integration. The total energy of the system would give the stability and accuracy of the solution. Results showed that the finite element method provided better efficiency in obtaining accurate solution than the finite difference method when a scalar processor is used. However, the finite difference method was more efficient on the vector processor. Therefore, the most efficient method of solution depends on the type of computer present for the analysis.

Moyer, E. T., Jr.↗

El Nino and outgoing longwave radiation - Observations from Nimbus-7 ERB

Nimbus-7 longwave radiation data are analyzed with respect to temporal and spatial changes in the anomalies induced by the El Nino event of 1982-1983. The components and functions of the earth radiation budget (ERB) instruments, which record satellite-altitude irradiances, are described. Local-noon analyses of outgoing longwave radiation (OLR) are presented; interannual variability is observed. El Nino-induced perturbations in the planetary-scale fields which affect OLR are studied. The onset, intensification, expansion, and decay of the OLR associated with El Nino are examined. A comparison of the ERS data with a three month average OLR anomaly map reveals that the ERS data have sufficient accuracy and stability with which to observe OLR anomalies induced by El Nino.

Ardanuy, P. E.↗

Numerical advection by conservation of second-order moments

A new, accurate, and nondiffusive method for three-dimensional advection of trace species is presented. The method preserves tracer structures by conserving the second-order moments of the spatial distribution of tracer during advection. Upstream transport and second-order tracer distribution are described, and the moments of the tracer distribution about the center of a grid box are formally defined and related to the polynomial distribution. Formulas are presented which describe how the moments of a grid box are decomposed into a unique set of moments centered about each subbox and how they are reassembled into new grid boxes. A one-dimensional example if tracer transport is given, and limits necessary to maintain positive tracer concentrations are derived. The accuracy and stability of this method are analytically examined, and numerical experiments testing the effective resolution are presented. The method is compared with other methods for numerical advection of tracers.

Prather, M. J.↗

Future payload isolation and pointing system technology

An overview is presented of the technology associated with the precision pointing of future spaceborne science instruments. High-level architectural options for dealing with the precision pointing problem are considered, and a representative NASA mission set for the 1990s is given. Pointing accuracy and stability requirements are extracted from the mission requirements, with the sub-0.1 arcsec stability requirement emerging as the primary driver. The state of the art of current technology is assessed, including an evaluation of gimbal systems, suspension systems, and actuator and sensor component technology. Areas where the technology needs to be pushed to satisfy future requirements are identified, and some promising design options are proposed.

Laskin, R. A.↗

Mathematical model for the Bridgman-Stockbarger crystal growing system

In a major technical breakthrough, a computer model for Bridgman-Stockbarger crystal growth was developed. The model includes melt convection, solute effects, thermal conduction in the ampule, melt, and crystal, and the determination of the curved moving crystal-melt interface. The key to the numerical method is the use of a nonuniform computational mesh which moves with the interface, so that the interface is a mesh surface. In addition, implicit methods are used for advection and diffusion of heat, concentration, and vorticity, for interface movement, and for internal gracity waves. This allows large time-steps without loss of stability or accuracy. Numerical results are presented for the interface shape, temperature distribution, and concentration distribution, in steady-state crystl growth. Solutions are presented for two test cases using water, with two different salts in solution. The two diffusivities differ by a factor of ten, and the concentrations differ by a factor of twenty.

Roberts, G. O.↗

Accuracy and long-term stability of amorphous-silicon measurements

The measurement system requirements to obtain accurate electrical performance measurements of amorphous silicon cells and modules were described. The progress achieved in modifying the Jet Propulsion Laboratory (JPL) system toward that objective were reviewed.

Mueller, R.↗

Decoupled direct method for sensitivity analysis in combustion kinetics

An efficient, decoupled direct method for calculating the first order sensitivity coefficients of homogeneous, batch combustion kinetic rate equations is presented. In this method the ordinary differential equations for the sensitivity coefficients are solved separately from , but sequentially with, those describing the combustion chemistry. The ordinary differential equations for the thermochemical variables are solved using an efficient, implicit method (LSODE) that automatically selects the steplength and order for each solution step. The solution procedure for the sensitivity coefficients maintains accuracy and stability by using exactly the same steplengths and numerical approximations. The method computes sensitivity coefficients with respect to any combination of the initial values of the thermochemical variables and the three rate constant parameters for the chemical reactions. The method is illustrated by application to several simple problems and, where possible, comparisons are made with exact solutions and those obtained by other techniques.

Radhakrishnan, Krishnan↗

Finite element solution techniques for large-scale problems in computational fluid dynamics

Element-by-element approximate factorization, implicit-explicit and adaptive implicit-explicit approximation procedures are presented for the finite-element formulations of large-scale fluid dynamics problems. The element-by-element approximation scheme totally eliminates the need for formation, storage and inversion of large global matrices. Implicit-explicit schemes, which are approximations to implicit schemes, substantially reduce the computational burden associated with large global matrices. In the adaptive implicit-explicit scheme, the implicit elements are selected dynamically based on element level stability and accuracy considerations. This scheme provides implicit refinement where it is needed. The methods are applied to various problems governed by the convection-diffusion and incompressible Navier-Stokes equations. In all cases studied, the results obtained are indistinguishable from those obtained by the implicit formulations.

Liou, J.↗

Higher-order flux difference splitting schemes for the Euler equations using upstream interpolations

A class of explicit two time-level, 2p + 1 space-point, (2p 1)-th order, upwind-biased flux difference splitting schemes are proposed for the numerical advection based on Lagrange's interpolation, and the method is an accord with the physical domain of dependence. A normalized Jacobian coefficient matrix is introduced to convert the schemes to hyperbolic systems of conservation laws, and approaches to make the higher-order schemes total variation stable are discussed. Accuracy and stability of the present schemes are examined, and implicit total variation diminishing schemes are developed for steady-state calculations.Application to gasdynamic problems for both steady and unsteady flows covering a wide range of Mach numbers is considered, and results for a blast wave passing a cylinder, and head-on collision of two blast waves over a circular arc, are presented. The flow patterns were found to be symmetric, and good resolution of flow structures was obtained.

Yang, J. Y.↗

Enhancements to Program LAURA for computation of three-dimensional hypersonic flow

Changes to Program Laura (Langley Aerothermodynamic Upwind Relaxation Algorithm) are presented which enhance both stability and accuracy of the algorithm. A discussion of iteration/sweeping strategies and their relation to computer architectures is included to best exploit the capabilities of serial, vector, and parallel processor machines. Test cases for Mach 10 perfect gas flow and Mach 32 real gas flow in chemical nonequilibrium over a blunt, raked elliptic cone using the thin-layer Navier-Stokes equations are presented in order to demonstrate the current improved capabilities. Algorithm changes include the use of volume averaging, application of a symmetric total variation diminishing (TVD) scheme, and stronger interaction between the grid/shock alignment routine and the relaxation algorithm. Good comparisons with heat transfer and pitching moment data at three different angles of attack for the Mach 10 tests serve to further validate the present algorithm. Parameters are defined which control the coupling of the specie continuity equations with the solution of the mixture conservation equations. A discussion of the consequences involved in the choice of strong versus weak coupling is presented, and a sample nonequilibrium calculation on a fine grid over a full scale model of the Aeroassist Flight Experiment (AFE) demonstrates current capabilities.

Gnoffo, Peter A.↗

Algorithm developments for the Euler equations with calculations of transonic flows

A new algorithm has been developed for the Euler equations that uses flux vector splitting in combination with the concept of rotating the coordinate system to the local streamwise direction. Flux vector biasing is applied along the local streamwise direction and central differencing is used transverse to the flow direction. The flux vector biasing is switched from upwind for supersonic flow to downwind-biased for subsonic flow. This switching is based on the Mach number; hence the proper domain of dependence is used in the supersonic regions and the switching occurs across shock waves. The theoretical basis and the development of the formulas for flux vector splitting are presented. Then several one-dimensional calculations are presented of steady and unsteady transonic flows, which demonstrate the stability and accuracy of the algorithm. Finally results are shown for unsteady transonic flow over an airfoil. The pressure coefficient plots show sharp transonic shock profiles, and the Mach contour plots show smoothly varying contours.

Goorjian, Peter M.↗

Numerical simulation of three dimensional transonic flows

The three-dimensional flow over a projectile has been computed using an implicit, approximately factored, partially flux-split algorithm. A simple composite grid scheme has been developed in which a single grid is partitioned into a series of smaller grids for applications which require an external large memory device such as the SSD of the CRAY X-MP/48, or multitasking. The accuracy and stability of the composite grid scheme has been tested by numerically simulating the flow over an ellipsoid at angle of attack and comparing the solution with a single grid solution. The flowfield over a projectile at M = 0.96 and 4 deg angle-of-attack has been computed using a fine grid, and compared with experiment.

Sahu, Jubaraj↗

Second-order upwind differencing in a recirculating flow

The accuracy and stability of the second order upwind differencing scheme was investigated. The solution algorithm employed is based on a coupled solution of the nonlinear finite difference equations by the multigrid technique. Calculations have been made of the driven cavity flow for several Reynolds numbers and finite difference grids. In comparison with the hybrid differencing, the second order upwind differencing is somewhat more accurate but it is not monotonically accurate with mesh refinement. Also, the convergence of the solution algorithm deteriorates with the use of the second order upwind differencing.

Vanka, S. P.↗

A forecast of new test capabilities using Magnetic Suspension and Balance Systems

This paper outlines the potential of Magnetic Suspension and Balance System (MSBS) technology to solve existing problems related to support interference in wind tunnels. Improvement of existing test techniques and exciting new techniques are envisioned as a result of applying MSBS. These include improved data accuracy, dynamic stability testing, two-body/stores release testing, and pilot/designer-in-the-loop tests. It also discusses the use of MSBS for testing exotic configurations such as hybrid hypersonic vehicles. A new facility concept that combines features of ballistic tubes, magnetic suspension, and cryogenic tunnels is described.

Lawing, Pierce L.↗

Aeroelastic modeling for the FIT team F/A-18 simulation

Some details of the aeroelastic modeling of the F/A-18 aircraft done for the Functional Integration Technology (FIT) team's research in integrated dynamics modeling and how these are combined with the FIT team's integrated dynamics model are described. Also described are mean axis corrections to elastic modes, the addition of nonlinear inertial coupling terms into the equations of motion, and the calculation of internal loads time histories using the integrated dynamics model in a batch simulation program. A video tape made of a loads time history animation was included as a part of the oral presentation. Also discussed is work done in one of the areas of unsteady aerodynamic modeling identified as needing improvement, specifically, in correction factor methodologies for improving the accuracy of stability derivatives calculated with a doublet lattice code.

Zeiler, Thomas A.↗