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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Numerical methods for stiff systems of two-point boundary value problems

Numerical procedures are developed for constructing asymptotic solutions of certain nonlinear singularly perturbed vector two-point boundary value problems having boundary layers at one or both endpoints. The asymptotic approximations are generated numerically and can either be used as is or to furnish a general purpose two-point boundary value code with an initial approximation and the nonuniform computational mesh needed for such problems. The procedures are applied to a model problem that has multiple solutions and to problems describing the deformation of thin nonlinear elastic beam that is resting on an elastic foundation.

Flaherty, J. E.↗
NASA NTRS1983

Numerical viscosity and the entropy condition for conservative difference schemes

Consider a scalar, nonlinear conservative difference scheme satisfying the entropy condition. It is shown that difference schemes containing more numerical viscosity will necessarily converge to the unique, physically relevant weak solution of the approximated conservation equation. In particular, entropy satisfying convergence follows for E schemes - those containing more numerical viscosity than Godunov's scheme.

Tadmor, E.↗
NASA NTRS1983

Numerical solution of a coupled pair of elliptic equations from solid state electronics

Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem.

Phillips, T. N.↗
NASA NTRS1983

Numerical computations of turbulence amplification in shock wave interactions

Numerical computations are presented which illustrate and test various effects pertinent to the amplification and generation of turbulence in shock wave turbulent boundary layer interactions. Several fundamental physical mechanisms are identified. Idealizations of these processes are examined by nonlinear numerical calculations. The results enable some limits to be placed on the range of validity of existing linear theories.

Zang, T. A.↗
NASA NTRS1983

Numerical solutions of Navier-Stokes equations for compressible turbulent two/three dimensional flows in terminal shock region of an inlet/diffuser

The multidimensional ensemble averaged compressible time dependent Navier Stokes equations in conjunction with mixing length turbulence model and shock capturing technique were used to study the terminal shock type of flows in various flight regimes occurring in a diffuser/inlet model. The numerical scheme for solving the governing equations is based on a linearized block implicit approach and the following high Reynolds number calculations were carried out: (1) 2 D, steady, subsonic; (2) 2 D, steady, transonic with normal shock; (3) 2 D, steady, supersonic with terminal shock; (4) 2 D, transient process of shock development and (5) 3 D, steady, transonic with normal shock. The numerical results obtained for the 2 D and 3 D transonic shocked flows were compared with corresponding experimental data; the calculated wall static pressure distributions agree well with the measured data.

Liu, N. S.↗
NASA NTRS1983

A comparative study of time-marching and space-marching numerical methods

Menees (1981) has conducted an evaluation of three different flowfield codes for the Jupiter entry conditions. However, a comparison of the codes has been made difficult by the fact that the three codes use different solution procedures, different computational mesh sizes, and a different convergence criterion. There are also other differences. For an objective evaluation of the different numerical solution methods employed by the codes, it would be desirable to select a simple no-blowing perfect-gas flowfield case for which the turbulent models are well established. The present investigation is concerned with the results of such a study. It is found that the choice of the numerical method is rather problem dependent. The time-marching and the space-marching method provide both comparable results if care is taken in selecting the appropriate mesh size near the body surface.

Gupta, R. N.↗
NASA NTRS1982

Implicit numerical integration for periodic solutions of autonomous nonlinear systems

A change of variables that stabilizes numerical computations for periodic solutions of autonomous systems is derived. Computation of the period is decoupled from the rest of the problem for conservative systems of any order and for any second-order system. Numerical results are included for a second-order conservative system under a suddenly applied constant load. Near the critical load for the system, a small increment in load amplitude results in a large increase in amplitude of the response.

Thurston, G. A.↗
NASA NTRS1982

Numerical evaluation of the Rayleigh integral for planar radiators using the FFT

Rayleigh's integral formula is evaluated numerically for planar radiators of any shape, with any specified velocity in the source plane using the fast Fourier transfrom algorithm. The major advantage of this technique is its speed of computation - over 400 times faster than a straightforward two-dimensional numerical integration. The technique is developed for computation of the radiated pressure in the nearfield of the source and can be easily extended to provide, with little computation time, the vector intensity in the nearfield. Computations with the FFT of the nearfield pressure of baffled rectangular plates with clamped and free boundaries are compared with the 'exact' solution to illuminate any errors. The bias errors, introduced by the FFT, are investigated and a technique is developed to significantly reduce them.

Williams, E. G.↗
NASA NTRS1982

A numerical study of the effects of ambipolar diffusion on the collapse of magnetic gas clouds

The gravitational collapse of isothermal, nonrotating magnetic gas clouds have been calculated numerically, including the effects of ambipolar diffusion. The fractional ionization in the clouds is approximated by a power-law function of the gas density, f = K/n to the q-power, where K and q are adjustable parameters. Eleven numerical experiments were run, and the results indicate that the asymptotic character of collapse is determined mainly by the value of q and is largely independent of the other parameters characterizing a cloud (e.g., K, cloud mass). In particular, there is nearly a one-to-one correspondence between q and the slope, x, of the central magnetic field strength-gas density relationship. If q is no more than 0.8, a cloud collapses asymptotically, as though the magnetic field were 'frozen' to the neutral matter. The magnetic field strength at the center of a collapsing cloud is strongly amplified during collapse even for values of q of about 1, despite extremely low values of fractional ionization. A discussion of the theoretical basis for this unexpected behavior is given. Possible implications of our results for the problems of magnetic braking of rotating protostars and star formation in general are also presented.

Black, D. C.↗
NASA NTRS1983

Numerical simulation of a disk-shaped electron accelerating electrostatic probe

This paper presents a method utilizing the numerical plasma simulation particle-in-cell technique for studying electron accelerating probe characteristics. A probe configuration other than the commonly used plane, cylindrical, or spherical case will be studied; the probe configuration used is that of a disk. Cylindrical geometry is used in the simulation model. One of the advantages of a numerical simulation is the ease with which the details of the system can be examined. Therefore, in addition to the current-voltage characteristics of the static probe, other results are presented such as the time evolution of the plasma properties and studies of the sheath region surrounding the probe.

Nonnast, J. H.↗
NASA NTRS1983

A numerical simulation of transition in plane channel flow

This paper involves a numerical simulation of the final stages of transition to turbulence in plane channel flow at a Reynolds number of 7500. Three-dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time-evolution of two- and three-dimensional finite-amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32 x 33 x 32 grid. Solutions indicate the existence of structures similar to those observed in the laboratory and which are characteristic of various stages of transition that lead to final breakdown. Details of the resulting flow field after breakdown indicate the evolution of streak-like formations found in turbulent flows. Although the flow field does approach a steady-state (turbulent channel flow), implementation of subgrid-scale terms are necesary to obtain proper turbulent statistics.

Biringen, S.↗
NASA NTRS1982

A numerical study of the two-dimensional Navier-Stokes equations in vorticity-velocity variables

The application of solution methods for compact finite-difference schemes to a vorticity-velocity form of the two-dimensional unsteady Navier-Stokes equations is described. An account is also given of numerical experiments for stagnation point and driven cavity flows. One experiment treats the flow impinging on a flat plate, and the numerical results can be compared with the analytical steady state solution (Batchelor, 1967). The other deals with the driven cavity problem, the primary purpose being to describe the time evolution of this classical problem (Donovan, 1970).

Gatski, T. B.↗
NASA NTRS1983

Importance of inlet boundary conditions for numerical simulation of combustor flows

Fluid dynamic computer codes for the mathematical simulation of problems in gas turbine engine combustion systems are required as design and diagnostic tools. To eventually achieve a performance standard with these codes of more than qualitative accuracy it is desirable to use benchmark experiments for validation studies. Typical of the fluid dynamic computer codes being developed for combustor simulations is the TEACH (Teaching Elliptic Axisymmetric Characteristics Heuristically) solution procedure. It is difficult to find suitable experiments which satisfy the present definition of benchmark quality. For the majority of the available experiments there is a lack of information concerning the boundary conditions. A standard TEACH-type numerical technique is applied to a number of test-case experiments. It is found that numerical simulations of gas turbine combustor-relevant flows can be sensitive to the plane at which the calculations start and the spatial distributions of inlet quantities for swirling flows.

Sturgess, G. J.↗
NASA NTRS1983

An experimental and numerical investigation of the impingement of an oblique shock wave on a body of revolution

Attention is given to the experimental study and numerical simulation of the impingement of oblique shock wave on a cylinder, in order both to document the complex three-dimensional shock wave and boundary layer interaction occurring in practical problems (such as stores carriage interference in a supersonic tactical aircraft) and to conduct a critical comparison of experimental measurements and numerical computations for such complex flows. A thin layer approximation of the Navier-Stokes equations was solved by means of a mixed explicit-implicit scheme. Experimental measurements reveal a highly complex flow field with two distinct adjacent separation zones, regions of high cross flow, and multiply reflected shocks and expansion fans.

Brosh, A.↗
NASA NTRS1983

Numerical solution of the Navier-Stokes equations for compressible turbulent two/three dimensional flows in the terminal shock region of an inlet/diffuser

The multidimensional, ensemble-averaged, compressible, time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique have been used to study the terminal shock type of flows in various flight regimes occurring in a diffuser/inlet model. The numerical scheme for solving the governing equations is based on a linearized block implicit approach and the following high Reynolds number calculations have been carried out: (1) 2-D, steady, subsonic; (2) 2-D, steady, transonic with normal shock, (3) 2-D, steady, supersonic with terminal shock, (4) 2-D, transient process of shock development and (5) 3-D, steady, transonic with normal shock. The numerical results obtained for the 2-D and 3-D transonic shocked flows have been compared with corresponding experimental data; the calculated wall static pressure distributions agree well with the measured data. The predicted transient responses of the flow to externally applied impulsive disturbance are consistent with results obtained via asymptotic analysis.

Liu, N.-S.↗
NASA NTRS1982

Numerical methods and computers used in elastohydrodynamic lubrication

Some of the methods of obtaining approximate numerical solutions to boundary value problems that arise in elastohydrodynamic lubrication are reviewed. The highlights of four general approaches (direct, inverse, quasi-inverse, and Newton-Raphson) are sketched. Advantages and disadvantages of these approaches are presented along with a flow chart showing some of the details of each. The basic question of numerical stability of the elastohydrodynamic lubrication solutions, especially in the pressure spike region, is considered. Computers used to solve this important class of lubrication problems are briefly described, with emphasis on supercomputers.

Hamrock, B. J.↗
NASA NTRS1983

Numerical modeling of turbulent flow in a channel

Two-dimensional incompressible turbulent flow in a channel with a backward-facing step was studied numerically by Chorin's Random Vortex Method (RVM), an algorithm capable of tracing the action of elementary turbulent eddies and their cumulative effects without imposing any restrictions upon their motions. The step occurs in one side of a channel with otherwise flat, parallel walls; its height equals 1/3, 1/4 or 1/5 the width of the channel downstream. The main objective was to investigate the behavior of the large-scale turbulent eddies in a flow and the flow characteristics in the separated shear layer, the reattached zone, and the rebuilding boundary layer after reattachment. The unsteady vorticity field and the distribution of time-averaged turbulent statistics were obtained. The effects of expansion step height and initial boundary layer state were also studied. Comparisons were made with the available experimental results. The agreement is satisfactory in the velocity profiles and in the reattachment length, and fairly good in the turbulence profiles. Also a mechanism of the development of the reattaching turbulent flow was suggested by the numerical results.

Dai, Y. W.↗