Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerics”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

On the numerical solution of time-dependent viscous incompressible fluid flows involving solid boundaries

An inherent numerical problem associated with the fully explicit pseudospectral numerical simulation of the incompressible Navier-Stokes equation for viscous flows with no-slip walls is described. A semi-implicit scheme which circumvents this numerical difficulty is presented. In this algorithm the equation of continuity rather than the Poisson equation for pressure is solved directly. Pseudospectral formulation of the channel flow problem using Fourier series and Chebyshev polynomials expansions is given for this scheme. An example demonstrating the applicability of the method is given.

Moin, P.↗
NASA NTRS1980

Viscous compressible flow about blunt bodies using a numerically generated orthogonal coordinate system

A numerical solution to the Navier-Stokes equations was obtained for blunt axisymmetric entry bodies of arbitrary shape in supersonic flow. These equations are solved on a finite difference mesh obtained from a simple numerical technique which generates orthogonal coordinates between arbitrary boundaries. The governing equations are solved in time dependent form using Stetter's improved stability three step predictor corrector method. For the present application, the metric coefficients were obtained numerically using fourth order accurate, finite difference relations and proved to be totally reliable for the highly stretched mesh used to resolve the thin viscous boundary layer. Solutions are obtained for a range of blunt body nose shapes including concavities.

Graves, R. A., Jr.↗
NASA NTRS1981

Numerical simulation of natural convection in a spherical container due to cooling at the center (idealization of the Lal/Kroes experiment)

Natural convection in a spherical container with cooling at the center was numerically simulated using a numerical fluid dynamics computer program. The numerical analysis was simplified by assuming axisymmetric flow in the spherical container, with the symmetry axis being a sphere diagonal parallel to the gravity vector. This axisymmetric spherical geometry was intended as an idealization of the proposed Lal/Kroes crystal growing experiment to be performed on Spacelab. Results were obtained for a range of Rayleigh numbers from 25 to 10,000. The computed velocities were found to be approximately proportional to the Rayleigh number over the range of Rayleigh numbers investigated.

Robertson, S. J.↗
NASA NTRS1981

Numerical Boundary Condition Procedures

Topics include numerical procedures for treating inflow and outflow boundaries, steady and unsteady discontinuous surfaces, far field boundaries, and multiblock grids. In addition, the effects of numerical boundary approximations on stability, accuracy, and convergence rate of the numerical solution are discussed.

Source record↗
NASA NTRS1981

Stability analysis of numerical boundary conditions and implicit difference approximations for hyperbolic equations

Implicit, noniterative, finite difference schemes were recently developed by several authors for multidimensional systems of nonlinear hyperbolic partial differential equations. When applied to linear model equations with periodic boundary conditions those schemes are unconditionally stable (A-stable). As applied in practice the algorithms often face a severe time step restriction. A major source of the difficulty is the treatment of the numerical boundary conditions. One conjecture was that unconditional stability requires implicit numerical boundary conditions. An apparent counter example was the space time extrapolation considered by Gustafsson, Kreiss, and Sunstrom. Spatial (implicit) and space time (explicit) extrapolation using normal mode analysis for a finite and infinite number of spatial mesh intervals are examined. The results indicate that for unconditional stability with a finite number of spatial mesh intervals, the numerical boundary conditions must be implicit.

Beam, R. M.↗
NASA NTRS1981

Numerical model of the moon's rotation

The differential equations for the Euler angles that define the moon's orientation with respect to an inertial coordinate frame have been integrated numerically. This numerical model of the moon's rotation, when used in conjunction with current lunar-orbit and earth-rotation models, fits lunar laser ranging observations over a five-year period within 28 cm (rms). The present model is also compared with the numerical model of Williams (1975) and Eckhardt's (1981) semianalytical model.

Cappalo, R. J.↗
NASA NTRS1981

A numerical study of two-dimensional shock vortex interaction

The numerical analysis of shock vortex interaction in a two-dimensional channel is discussed in this paper. The Euler equations and the two-step MacCormack scheme are employed in this numerical analysis. The comparison of the numerical results with the theoretical analyses and experimental results in the literature has provided unique insights into the dynamics of the transient interaction process and the needs for further theoretical developments. The results are useful in applications to problems in aeroacoustics and the transient aerodynamics of high performance aircraft.

Pao, S. P.↗
NASA NTRS1981

Manufacturing in space: Fluid dynamics numerical analysis

Natural convection in a spherical container with cooling at the center was numerically simulated using the Lockheed-developed General Interpolants Method (GIM) numerical fluid dynamic computer program. The numerical analysis was simplified by assuming axisymmetric flow in the spherical container, with the symmetry axis being a sphere diagonal parallel to the gravity vector. This axisymmetric spherical geometry was intended as an idealization of the proposed Lal/Kroes growing experiments to be performed on board Spacelab. Results were obtained for a range of Rayleigh numbers from 25 to 10,000. For a temperature difference of 10 C from the cooling sting at the center to the container surface, and a gravitional loading of 0.000001 g a computed maximum fluid velocity of about 2.4 x 0.00001 cm/sec was reached after about 250 sec. The computed velocities were found to be approximately proportional to the Rayleigh number over the range of Rayleigh numbers investigated.

Robertson, S. J.↗
NASA NTRS1982

The numerical analysis of the rotational theory for the formation of lunar globules

The morphology of lunar globules is studied through the application of a numerical analysis of their rotation in space during cooling. It is assumed that molten rock is shot from the surface of the moon, solidifies in space above the moon and then falls back to the surface. The rotational theory studied makes the following assumptions: the volume of the molten rock does not change during cooling; the angular momentum is conserved; there are no internal motions because of the high viscosity of the molten rock, i.e., in equilibrium the globule is rotating as a rigid body; finally, the kinetic reaction of the globule to the forces is fast relative to the rate of cooling, i.e., the globule reaches equilibrium at constant energy. These assumptions are subjected to numerical analysis yielding good agreement between the actual globule shapes and the numerical results, but leaving some doubt as to the validity of the rotational theory due to the failure to establish the existence of true local minima and an incomplete understanding of the thermokentics.

Ross, J.↗
NASA NTRS1982

Engine dynamic analysis with general nonlinear finite element codes. Part 2: Bearing element implementation overall numerical characteristics and benchmaking

Finite element codes are used in modelling rotor-bearing-stator structure common to the turbine industry. Engine dynamic simulation is used by developing strategies which enable the use of available finite element codes. benchmarking the elements developed are benchmarked by incorporation into a general purpose code (ADINA); the numerical characteristics of finite element type rotor-bearing-stator simulations are evaluated through the use of various types of explicit/implicit numerical integration operators. Improving the overall numerical efficiency of the procedure is improved.

Padovan, J.↗
NASA NTRS1981

Simulation of a numerical filter for enhancing earth radiation budget measurements

The Earth Radiation Budget Experiment has the objective to collect the radiation budget data which are needed to determine the radiation budget at the top of the atmosphere (TOA) on a regional scale. A second objective is to determine the accuracy of the results. Three satellites will carry wide and medium field of view radiometers which measure the longwave and shortwave components of radiation. Scanning radiometers will be included to detect small spatial features. A proposal has been made to employ for the nonscanning radiometers a one-dimensional numerical filter which reduces satellite measurements to TOA radiant excitances. The numerical filter was initially formulated by House (1980). It enhances the resolution of the radiation budget along the satellite groundtrack. The accuracy of the numerical filter estimate is studied by simulating the data gathering and measurement inversion process. The results of the study are discussed, taking into account two error sources.

Green, R. N.↗
NASA NTRS1982

Numerical simulation of viscous-inviscid interactions on indented nose tips

An implicit numerical algorithm to solve the unsteady thin-layer Navier-Stokes equations in a strong conservative form has been used to compute the viscous flow over indented nose tips placed in a supersonic free stream. Numerical solutions are presented for axisymmetric and three-dimensional indented configurations for laminar flow conditions. Results demonstrate the capability of the present numerical procedure to predict flow fields that contain strong viscous-inviscid interactions, including boundary-layer separation, reattachment, and embedded discontinuities. Comparisons with available experimental data for the shock shape and surface pressure are also given.

Rizk, Y. M.↗
NASA NTRS1981

International Conference on Numerical Methods in Fluid Dynamics, 7th, Stanford University, Stanford and Moffett Field, CA, June 23-27, 1980, Proceedings

Topics discussed include polygon transformations in fluid mechanics, computation of three-dimensional horseshoe vortex flow using the Navier-Stokes equations, an improved surface velocity method for transonic finite-volume solutions, transonic flow calculations with higher order finite elements, the numerical calculation of transonic axial turbomachinery flows, and the simultaneous solutions of inviscid flow and boundary layer at transonic speeds. Also considered are analytical solutions for the reflection of unsteady shock waves and relevant numerical tests, reformulation of the method of characteristics for multidimensional flows, direct numerical simulations of turbulent shear flows, the stability and separation of freely interacting boundary layers, computational models of convective motions at fluid interfaces, viscous transonic flow over airfoils, and mixed spectral/finite difference approximations for slightly viscous flows.

Reynolds, W. C.↗
NASA NTRS1982

Error propagation in the numerical solutions of the differential equations of orbital mechanics

The relationship between the eigenvalues of the linearized differential equations of orbital mechanics and the stability characteristics of numerical methods is presented. It is shown that the Cowell, Encke, and Encke formulation with an independent variable related to the eccentric anomaly all have a real positive eigenvalue when linearized about the initial conditions. The real positive eigenvalue causes an amplification of the error of the solution when used in conjunction with a numerical integration method. In contrast an element formulation has zero eigenvalues and is numerically stable.

Bond, V. R.↗
NASA NTRS1982

Numerical simulation of magnetohydrodynamic shock propagation in the corona

Recent developments in the field of numerical simulation models for the study of shock wave propagation in the corona are presented. These models are based on gasdynamic (GD) and ideal (that is, dissipationless, except at shocks) magnetohydrodynamic (MHD) theories. The characteristics and physical interpretations of the results derived from these models are discussed in some detail. The most significant physical results obtained to date are provided by the two-dimensional non-planar, time-dependent, MHD numerical simulation model. In this model, the non-linear interaction among the three essential MHD waves, i.e., fast-, slow-, and Alfven waves are demonstrated. Finally, the physical relevance of these numerical simulation models in relation to observed solar activity is presented.

Wu, S. T.↗
NASA NTRS1982

On the energy dependence of the radial diffusion coefficient and spectra of inner radiation belt particles - Analytic solutions and comparison with numerical results

A theoretical method by which the energy dependence of the radial diffusion coefficient may be deduced from spectral observations of the particle population at the inner edge of the earth's radiation belts is presented. This region has previously been analyzed with numerical techniques; in this report an analytical treatment that illustrates characteristic limiting cases in the L shell range where the time scale of Coulomb losses is substantially shorter than that of radial diffusion (L approximately 1-2) is given. It is demonstrated both analytically and numerically that the particle spectra there are shaped by the energy dependence of the radial diffusion coefficient regardless of the spectral shapes of the particle populations diffusing inward from the outer radiation zone, so that from observed spectra the energy dependence of the diffusion coefficient can be determined. To insure realistic simulations, inner zone data obtained from experiments on the DIAL, AZUR, and ESRO 2 spacecraft have been used as boundary conditions. Excellent agreement between analytic and numerical results is reported.

Westphalen, H.↗
NASA NTRS1983

Numerical computation of transonic flow governed by the full-potential equation

Numerical solution techniques for solving transonic flow fields governed by the full potential equation are discussed. In a general sense relaxation schemes suitable for the numerical solution of elliptic partial differential equations are presented and discussed with emphasis on transonic flow applications. The presentation can be divided into two general categories: An introductory treatment of the basic concepts associated with the numerical solution of elliptic partial differential equations and a more advanced treatment of current procedures used to solve the full potential equation for transonic flow fields. The introductory material is presented for completeness and includes a brief introduction (Chapter 1), governing equations (Chapter 2), classical relaxation schemes (Chapter 3), and early concepts regarding transonic full potential equation algorithms (Chapter 4).

Holst, T. L.↗
NASA NTRS1983

A numerical method based on the Fourier-Fourier transform approach for modeling 1-D electron plasma evolution

A numerical method is presented for studying one-dimensional electron plasma evolution under typical interplanetary conditions. The method applies the Fourier-Fourier transform approach to a plasma model that is a generalization of the electrostatic Vlasov-Poisson system of equations. Conservation laws that are modified to include the plasma model generalization and also the boundary effects of nonperiodic solutions are given. A new conservation law for entropy in the transformed space is then introduced. These conservation laws are used to verify the numerical solutions. A discretization error analysis is presented. Two numerical instabilities and the methods used for their suppression are treated. It is shown that in interplanetary plasma conditions, the bump-on-tail instability produces significant excitation of plasma oscillations at the Bohm-Gross frequency and its second harmonic. An explanation of the second harmonic excitation is given in terms of wave-wave coupling during the growth phase of the instability.

Klimas, A. J.↗