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.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Numerical analysis of the transient response of an axisymmetric ablative char layer considering internal flow effects

The differential equations governing the transient response of the char layer of an ablating axisymmetric body, internal pyrolysis gas flow effects being considered, have been derived. These equations have been expanded into finite difference form and programed for numerical solution on a digital computer. Numerical results compare favorably with simplified exact solutions. The complete numerical analysis was used to obtain solutions for two representative body shapes subjected to a typical entry heating environment. Pronounced effects of the lateral flow of pyrolysis gases on the mass flow field within the char layer and the associated surface and pyrolysis interface recession rates are shown.

Pittman, C. M.↗

NASA NTRS1972

Numerical method and FORTRAN program for the solution of an axisymmetric electrostatic collector design problem

The numerical calculation is described of the steady-state flow of electrons in an axisymmetric, spherical, electrostatic collector for a range of boundary conditions. The trajectory equations of motion are solved alternately with Poisson's equation for the potential field until convergence is achieved. A direct (noniterative) numerical technique is used to obtain the solution to Poisson's equation. Space charge effects are included for initial current densities as large as 100 A/sq cm. Ways of dealing successfully with the difficulties associated with these high densities are discussed. A description of the mathematical model, a discussion of numerical techniques, results from two typical runs, and the FORTRAN computer program are included.

Reese, O. W.↗

NASA NTRS1972

Numerical method for the solution of large systems of differential equations of the boundary layer type

A numerical method for the solution of large systems of nonlinear differential equations of the boundary-layer type is described. The method is a modification of the technique for satisfying asymptotic boundary conditions. The present method employs inverse interpolation instead of the Newton method to adjust the initial conditions of the related initial-value problem. This eliminates the so-called perturbation equations. The elimination of the perturbation equations not only reduces the user's preliminary work in the application of the method, but also reduces the number of time-consuming initial-value problems to be numerically solved at each iteration. For further ease of application, the solution of the overdetermined system for the unknown initial conditions is obtained automatically by applying Golub's linear least-squares algorithm. The relative ease of application of the proposed numerical method increases directly as the order of the differential-equation system increases. Hence, the method is especially attractive for the solution of large-order systems. After the method is described, it is applied to a fifth-order problem from boundary-layer theory.

Green, M. J.↗

NASA NTRS1972

Numerical stability in problems of linear algebra.

Mathematical problems are introduced as mappings from the space of input data to that of the desired output information. Then a numerical process is defined as a prescribed recurrence of elementary operations creating the mapping of the underlying mathematical problem. The ratio of the error committed by executing the operations of the numerical process (the roundoff errors) to the error introduced by perturbations of the input data (initial error) gives rise to the concept of lambda-stability. As examples, several processes are analyzed from this point of view, including, especially, old and new processes for solving systems of linear algebraic equations with tridiagonal matrices. In particular, it is shown how such a priori information can be utilized as, for instance, a knowledge of the row sums of the matrix. Information of this type is frequently available where the system arises in connection with the numerical solution of differential equations.

Babuska, I.↗

NASA NTRS1973

Feasibility study of the numerical integration of shell equations using the field method

The field method is developed for arbitrary open branch domains subjected to general linear boundary conditions. Although closed branches are within the scope of the method, they are not treated here. The numerical feasibility of the method has been demonstrated by implementing it in a computer program for the linear static analysis of open branch shells of revolution under asymmetric loads. For such problems the field method eliminates the well-known numerical problem of long subintervals associated with the rapid growth of extraneous solutions. Also, the method appears to execute significantly faster than other numerical integration methods.

Cohen, G. A.↗

NASA NTRS1972

On the numerical solution of two-dimensional, laminar compressible flows with imbedded shock waves.

The complete, time-dependent Navier-Stokes equations are expressed in conservation form and solved by employing an explicit finite difference numerical technique which incorporates artificial viscosity terms of the form first suggested by Rusanov for numerical stability in the vicinity of shock waves. Surface boundary conditions are developed in a consistent and unique manner through the use of a physically oriented extrapolation procedure. From numerical experimentation an extended range for the explicit stability parameter is established. Also employed is an additional convergence parameter which relates incremental spatial steps. Convergence of the transient solution to a steady state flow was obtained after 400 to 500 time steps.

Goodrich, W. D.↗

NASA NTRS1973

Numerical integration of shell equations using the field method.

The 'field method' for the numerical solution of even-order linear boundary-value problems in ordinary differential equations is formulated. This method converts the boundary-value problem into two successive initial-value problems, which may be solved by standard forward integration techniques. The method has been implemented in a computer program to calculate the static response of ring-stiffened branched shells of revolution to asymmetric loads. For such problems the field method eliminates the well-known numerical problem of 'long subintervals' and also executes operations significantly faster than other numerical integration methods.

Cohen, G. A.↗

NASA NTRS1974

Numerical calculations of velocity and pressure distribution around oscillating airfoils

An analytical procedure based on the Navier-Stokes equations was developed for analyzing and representing properties of unsteady viscous flow around oscillating obstacles. A variational formulation of the vorticity transport equation was discretized in finite element form and integrated numerically. At each time step of the numerical integration, the velocity field around the obstacle was determined for the instantaneous vorticity distribution from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints, using the Lagrangian Multiplier Technique, at each time step of the numerical integration. The procedure was then applied for determining pressures around obstacles oscillating in unsteady flow. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗

NASA NTRS1974

Numerical calculations of two dimensional, unsteady transonic flows with circulation

The feasibility of obtaining two-dimensional, unsteady transonic aerodynamic data by numerically integrating the Euler equations is investigated. An explicit, third-order-accurate, noncentered, finite-difference scheme is used to compute unsteady flows about airfoils. Solutions for lifting and nonlifting airfoils are presented and compared with subsonic linear theory. The applicability and efficiency of the numerical indicial function method are outlined. Numerically computed subsonic and transonic oscillatory aerodynamic coefficients are presented and compared with those obtained from subsonic linear theory and transonic wind-tunnel data.

Beam, R. M.↗

NASA NTRS1973

A numerical procedure for the parametric optimization of three dimensional scramjet nozzles

A numerical procedure permitting the rapid determination of the internal performance of a class of scramjet nozzle configurations is presented. The approach developed is based on the construction of quasi two dimensional simple wave networks, where lateral expansion effects are incorporated by one dimensional approximations. A numerical procedure following this approach has has been developed and results obtained are highly comparable to those obtained employing a characteristic procedure. The numerical program developed permits the parametric variation of cowl length, turning angles on the cowl and vehicle undersurface and lateral expansion and is subject to fixed constraints such as the vehicle length and nozzle exit height. The program requires uniform initial conditions at the burner exit station and yields the location of all predominant wave zones, accounting for lateral expansion effects. In addition, the program yields the detailed pressure distribution on the cowl and vehicle undersurface and calculates the nozzle thrust, lift and pitching moment.

Dash, S.↗

NASA NTRS1973

An improved numerical process for solution of solid mechanics problems

This paper gives an overview of the development and status of an improved numerical process for the solution of solid mechanics problems. The proposed process uses a mixed formulation with the fundamental unknowns consisting of both stress and displacement parameters. The problem is formulated either by means of first-order partial differential equations or in a variational form by using a Hellinger-Reissner-type mixed variational principle. For presentation purposes, the components of a numerical process are characterized and the criteria for an ideal process are outlined. Commonly used finite-difference and finite-element procedures are examined in the light of these criteria and it is shown that they fall short in a number of ways. The proposed numerical process, on the other hand, satisfies most of the optimality criteria and appears to be particularly suited for use with the forthcoming generation computers.

Noor, A. K.↗

NASA NTRS1974

A numerical-perturbation method for the nonlinear analysis of structural vibrations

A numerical-perturbation method is proposed for the determination of the nonlinear forced response of structural elements. Purely analytical techniques are capable of determining the response of structural elements having simple geometries and simple variations in thickness and properties, but they are not applicable to elements with complicated structure and boundaries. Numerical techniques are effective in determining the linear response of complicated structures, but they are not optimal for determining the nonlinear response of even simple elements when modal interactions take place due to the complicated nature of the response. Therefore, the optimum is a combined numerical and perturbation technique. The present technique is applied to beams with varying cross sections.

Nayfeh, A. H.↗

NASA NTRS1974

Finite elasto-plastic deformation. I - Theory and numerical examples

It is demonstrated that the problem of elasto-plastic finite deformation is governed by a quasi-linear model irrespective of deformation magnitude. This feature follows from the adoption of a rate viewpoint toward finite deformation analysis in an Eulerian reference frame. Objectivity of the formulation is preserved by introduction of a frame-invariant stress rate. Equations for piece-wise linear incremental finite element analysis are developed by application of the Galerkin method to the instantaneously linear governing differential equations of the quasi-linear model. Numerical solution capability has been established for problems of plane strain and plane stress. The accuracy of the numerical analysis is demonstrated by consideration of a number of problems of homogeneous finite deformation admitting comparative analytic solution. It is shown that accurate, objective numerical solutions can be obtained for problems involving dimensional changes of an order of magnitude and rotations of a full radian.

Osias, J. R.↗

NASA NTRS1974

Use of Green's functions in the numerical solution of two-point boundary value problems

This study investigates the use of Green's functions in the numerical solution of the two-point boundary value problem. The first part deals with the role of the Green's function in solving both linear and nonlinear second order ordinary differential equations with boundary conditions and systems of such equations. The second part describes procedures for numerical construction of Green's functions and considers briefly the conditions for their existence. Finally, there is a description of some numerical experiments using nonlinear problems for which the known existence, uniqueness or convergence theorems do not apply. Examples here include some problems in finding rendezvous orbits of the restricted three body system.

Gallaher, L. J.↗

NASA NTRS1974

Automatic numerical generation of body-fitted curvilinear coordinate system for field containing any number of arbitrary two-dimensional bodies

A method for automatic numerical generation of a general curvilinear coordinate system with coordinate lines coincident with all boundaries of a general multi-connected region containing any number of arbitrarily shaped bodies is presented. With this procedure the numerical solution of a partial differential system may be done on a fixed rectangular field with a square mesh with no interpolation required regardless of the shape of the physical boundaries, regardless of the spacing of the curvilinear coordinate lines in the physical field, and regardless of the movement of the coordinate system. Numerical solutions for the lifting and nonlifting potential flow about Joukowski and Karman-Trefftz airfoils using this coordinate system generation show excellent comparison with the analytic solutions. The application to fields with multiple bodies is illustrated by a potential flow solution for multiple airfoils.

Thompson, J. F.↗

NASA NTRS1975

Numerical techniques for the solution of the compressible Navier-Stokes equations and implementation of turbulence models

The time-splitting explicit numerical method of MacCormack is applied to separated turbulent boundary layer flow problems. Modifications of this basic method are developed to counter difficulties associated with complicated geometry and severe numerical resolution requirements of turbulence model equations. The accuracy of solutions is investigated by comparison with exact solutions for several simple cases. Procedures are developed for modifying the basic method to improve the accuracy. Numerical solutions of high-Reynolds-number separated flows over an airfoil and shock-separated flows over a flat plate are obtained. A simple mixing length model of turbulence is used for the transonic flow past an airfoil. A nonorthogonal mesh of arbitrary configuration facilitates the description of the flow field. For the simpler geometry associated with the flat plate, a rectangular mesh is used, and solutions are obtained based on a two-equation differential model of turbulence.

Baldwin, B. S.↗

NASA NTRS1975

A numerical method for solving the Navier-Stokes equations with application to shock-boundary layer interactions

A numerical method for solving the compressible form of the unsteady Navier-Stokes equations is described. This method was originally presented in 1970 and has since been modified during the development of computer programs at Ames for implementing models that account for the effects of turbulence in shock-induced separated flows. Although this paper does not describe the turbulence models themselves, a complete description of the basic numerical method is given with emphasis on the choice of a computational mesh for high Reynolds number flows, finite-difference approximations for mixed partial derivatives, extension of the Courant-Friedrichs-Lewy stability condition for viscous flows, mesh boundary conditions, and numerical smoothing for strong shock-wave calculations.

Maccormack, R. W.↗

NASA NTRS1976

A method for generating numerical pilot opinion ratings using the optimal pilot model

A method for generating numerical pilot opinion ratings using the optimal pilot model is introduced. The method is contained in a rating hypothesis which states that the numerical rating which a human pilot assigns to a specific vehicle and task can be directly related to the numerical value of the index of performance resulting from the optimal pilot modeling procedure as applied to that vehicle and task. The hypothesis is tested using the data from four piloted simulations. The results indicate that the hypothesis is reasonable, but that the predictive capability of the method is a strong function of the accuracy of the pilot model itself. This accuracy is, in turn, dependent upon the parameters which define the optimal modeling problem. A procedure for specifying the parameters for the optimal pilot model in the absence of experimental data is suggested.

Hess, R. A.↗