Search NASA⌕ Search

SEARCH · Search NASA

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 469 records · Page 26

Comparison of approximate and numerical analyses of nonlinear combustion instability

At the present time, there are three general analytical techniques available to study problems of unsteady motions in rocket motors: linear stability analysis; approximate nonlinear analysis, founded on examining the behavior of coupled normal modes; and numerical calculations based on the conservation equations for one-dimensional flows. The last two yield the linear results as a limit. It is the main purpose of this paper to check the accuracy of the approximate analysis against the numerical analysis for some special cases. The results provide some justification for using the approximate analysis to study three-dimensional problems.

Culick, F. E. C.↗

Recent advances in numerical analysis of structural eigenvalue problems

A wide range of eigenvalue problems encountered in practical structural engineering analyses is defined, in which the structures are assumed to be discretized by any suitable technique such as the finite-element method. A review of the usual numerical procedures for the solution of such eigenvalue problems is presented and is followed by an extensive account of recently developed eigenproblem solution procedures. Particular emphasis is placed on the new numerical algorithms and associated computer programs based on the Sturm sequence method. Eigenvalue algorithms developed for efficient solution of natural frequency and buckling problems of structures are presented, as well as some eigenvalue procedures formulated in connection with the solution of quadratic matrix equations associated with free vibration analysis of structures. A new algorithm is described for natural frequency analysis of damped structural systems.

Gupta, K. K.↗

Numerical evaluation of electromagnetic fields due to dipole antennas in the presence of stratified media

Two numerical methods are used to evaluate the integrals that express the em fields due to dipole antennas radiating in the presence of a stratified medium. The first method is a direct integration by means of Simpson's rule. The second method is indirect and approximates the kernel of the integral by means of the fast Fourier transform. In contrast to previous analytical methods that applied only to two-layer cases the numerical methods can be used for any arbitrary number of layers with general properties.

Tsang, L.↗

The numerical calculation of inviscid plume flow fields

A numerical method is presented for the computation of inviscid, axisymmetric, underexpanded plumes. The numerical techniques developed by Moretti (1969, 1971, 1972) are used in conjunction with Abbett's (1970) theory for locating the Mach disk. Abbett's theory and three other prediction methods are compared to experimental results. Results are presented for jets exhausting into static ambients and supersonic free stream. Detailed results from a flow field with multiple Mach disks are also presented. Finally, some problems associated with the computation of very small Mach disks are discussed.

Salas, M. D.↗

Numerical methods for the design and analysis of wings at supersonic speeds

Numerical methods for the design and analysis of arbitrary-planform wings at supersonic speeds are reviewed. Certain deficiencies are revealed, particularly in application to wings with slightly subsonic leading edges. Recently devised numerical techniques which overcome the major part of these deficiencies are presented. The original development as well as the more recent revisions are subjected to a thorough review.

Carlson, H. W.↗

Numerical solution for unsteady sonic flow over thin wings

A numerical solution precedure of a simplified unsteady transonic equation which is fast, reasonably accurate, and takes into account many of the effects of the steady flow field is described. The numeric solution of this equation is accurate and is accomplished on an IBM 360/65 computer. Arbitrary planform shape is accommodated and variable local Mach number effects from the steady flow are easily handled.

Kimble, K. R.↗

Numerical computation of viscous blunt body flows with a planar impinging shock

Two- and three-dimensional, viscous blunt body flows with planar impinging shocks are computed using an explicit, time-dependent, finite-difference method to solve the complete set of Navier-Stokes equations. The bow shock is treated as a discontinuity, while all interior shock layer detail such as shear layers, shock waves, jets and the wall boundary layer are automatically captured in the solution. Numerical results are presented for cases in which planar shock waves of different strengths and orientations are allowed to impinge on the flow field surroundings an infinite cylinder resulting in two- and three-dimensional shock interference patterns. The numerical results are compared with experiment.

Holst, T. L.↗

Numerical solution of the Navier-Stokes equations for arbitrary two-dimensional multi-element airfoils

Abstracts are presented on a method of numerical solution of the Navier-Stokes equation for the flow about arbitrary airfoils, using a numerically generated curvilinear coordinate system having a coordinate line coincident with the body contour. Results of continuing research are reported and include: application of the Navier-Stokes solution in the vorticity-stream function formulation to a number of single airfoils at Reynolds numbers up to 2000; programming of the Navier-Stokes solution for multiple airfoils in the primitive variable formulation; testing of the potential flow solution of multiple bodies; and development of a generalized coordinate system program.

Thompson, J. F.↗

Numerical computation of three-dimensional blunt body flow fields with an impinging shock

A time-marching finite-difference method was used to solve the compressible Navier-Stokes equations for the three-dimensional wing-leading-edge shock impingement problem. The bow shock was treated as a discontinuity across which the exact shock jump conditions were applied. All interior shock layer detail such as shear layers, shock waves, jets, and the wall boundary layer were automatically captured in the solution. The impinging shock was introduced by discontinuously changing the freestream conditions across the intersection line at the bow shock. A special storage-saving procedure for sweeping through the finite-difference mesh was developed which reduces the required amount of computer storage by at least a factor of two without sacrificing the execution time. Numerical results are presented for infinite cylinder blunt body cases as well as the three-dimensional shock impingement case. The numerical results are compared with existing experimental and theoretical results.

Holst, T. L.↗

An improved numerical procedure for the parametric optimization of three dimensional scramjet nozzles

A parametric numerical procedure permitting the rapid determination of the performance of a class of scramjet nozzle configurations is presented. The geometric complexity of these configurations ruled out attempts to employ conventional nozzle design procedures. The numerical program developed permitted the parametric variation of cowl length, turning angles on the cowl and vehicle undersurface and lateral expansion, and was subject to fixed constraints such as the vehicle length and nozzle exit height. The program required uniform initial conditions at the burner exit station and yielded the location of all predominant wave zones, accounting for lateral expansion effects. In addition, the program yielded the detailed pressure distribution on the cowl, vehicle undersurface and fences, if any, and calculated the nozzle thrust, lift and pitching moments.

Dash, S.↗

High-order numerical solutions using cubic splines

The cubic spline collocation procedure for the numerical solution of partial differential equations was reformulated so that the accuracy of the second-derivative approximation is improved and parallels that previously obtained for lower derivative terms. The final result is a numerical procedure having overall third-order accuracy for a nonuniform mesh and overall fourth-order accuracy for a uniform mesh. Application of the technique was made to the Burger's equation, to the flow around a linear corner, to the potential flow over a circular cylinder, and to boundary layer problems. The results confirmed the higher-order accuracy of the spline method and suggest that accurate solutions for more practical flow problems can be obtained with relatively coarse nonuniform meshes.

Rubin, S. G.↗

Higher order accurate partial implicitization: An unconditionally stable fourth-order-accurate explicit numerical technique

The previously obtained second-order-accurate partial implicitization numerical technique used in the solution of fluid dynamic problems was modified with little complication to achieve fourth-order accuracy. The Von Neumann stability analysis demonstrated the unconditional linear stability of the technique. The order of the truncation error was deduced from the Taylor series expansions of the linearized difference equations and was verified by numerical solutions to Burger's equation. For comparison, results were also obtained for Burger's equation using a second-order-accurate partial-implicitization scheme, as well as the fourth-order scheme of Kreiss.

Graves, R. A., Jr.↗

Short-periodic variations and second-order numerical averaging

The principal disadvantage of the method of numerical averaging is that it provides only the average time history of the orbital elements and yields no information about the high-frequency or short-periodic variations that occur inside the averaging interval. This paper contains a description of a technique for recovering the short-periodic variations by minor modifications to the averaging process so as to permit the construction of a Fourier series for the osculating elements. The availability of this series permits the extension of the averaging technique to higher order and allows us to account for short-periodic coupling of the perturbations. Comparisons of the results with numerically integrated solutions are presented for three distinct orbit prediction problems.

Lutzky, D.↗

Recent progress in numerical techniques for flow simulation

Recent developments in the use of numerical methods for fluid flow simulation show an increasing tendency to use numerical operators that can, in various ways, be factored. Use of methods having this property often increases the accuracy and efficiency of computer codes. Tracing the factorization property provides a unity to the basic concepts involved in the development of cyclic reduction, predictor-corrector, splitting, fast Fourier transform, and pseudospectral methods.-

Lomax, H.↗

Higher-order numerical solutions using cubic splines

A cubic spline collocation procedure has recently been developed for the numerical solution of partial differential equations. In the present paper, this spline procedure is reformulated so that the accuracy of the second-derivative approximation is improved and parallels that previously obtained for lower derivative terms. The final result is a numerical procedure having overall third-order accuracy for a non-uniform mesh and overall fourth-order accuracy for a uniform mesh. Solutions using both spline procedures, as well as three-point finite difference methods, will be presented for several model problems.-

Rubin, S. G.↗

Use of numerically generated body-fitted coordinate systems for solution of the Navier-Stokes equations

A procedure for numerical solution of the time-dependent, two-dimensional incompressible Navier-Stokes equations that can treat the unsteady laminar flow about bodies of arbitrary shape, such as two-dimensional airfoils, multiple airfoils, and submerged hydrofoils, as naturally as it can deal with the flow about simple bodies. The solution is based on a method of automatic numerical generation of a general curvilinear coordinate system with coordinate lines coincident with all boundaries of a general multiconnected region containing any number of arbitrarily shaped bodies. The curvilinear coordinates are generated as the solution of two elliptical partial differential equations with Dirichlet boundary conditions, one coordinate being specified to be constant on each of the boundaries, and a distribution of the other being specified along the boundaries. The solution compares excellently with the Blasius boundary layer solution for the flow past a semiinfinite flat plate.

Thompson, J. F.↗

Evolution of rotating interstellar clouds. I - Numerical techniques

A method is described for the numerical calculation of the hydrodynamic evolution of a self-gravitating configuration in two space dimensions with assumed axial symmetry. The calculation is formulated in cylindrical coordinates with respect to a moving Eulerian grid and is solved using explicit hydrodynamics combined with implicit radiative transfer. The physics included is appropriate for calculation of the collapse of a rotating protostellar cloud. The gravitational field is obtained by means of an alternating-direction iterative technique. Numerical tests to demonstrate the correctness of the method are presented for special cases.

Black, D. C.↗

Numerical procedure for analyzing Langmuir probe data

A numerical procedure is proposed for analyzing Langmuir probe data in the presence of a two-group plasma containing both primary and Maxwellian electrons. The procedure is known as a least-squares differential-correction technique for determining the unknown coefficients of the governing equation. It is shown that for a given set of input data the results of the analysis are unique and independent of the initial estimate of electron temperature, that convergence is fastest when electron temperature is overestimated, that the results are sensitive to the region of the curve used as input and also to the voltage increment, and that plasma properties determined by the proposed numerical procedure are either consistent with those determined graphically or closer to expected values. With a suitable data acquisition system, the Fortran IV program worked out for this procedure could be used to provide real-time plasma diagnostic information for an operating ion thruster.

Beattie, J. R.↗