Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Methods”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Space shuttle auxiliary propulsion system design study. Phase C report: Oxygen-hydrogen RCS/OMS integration study

A comparison of the concepts of auxiliary propulsion systems proposed for the space shuttle vehicle is discussed. An evaluation of the potential of integration between the reaction control system and the orbit maneuvering system was conducted. Numerous methods of implementing the various levels of integration were evaluated. Preferred methods were selected and design points were developed for two fully integrated systems, one partially integrated system, and one separate system.

Bruns, A. E.↗
NASA NTRS1973

Behavior of nitric oxide formed by the space shuttle in the mesosphere

A numerical model of eddy diffusive transport of emitted gases from the space shuttle wake, including chemical reactions between the emitted constituents and the ambient atmosphere, has been constructed for 75 km altitude. The numerical methods involve explicit solution of the diffusion equation and Runge-Kutta method for the chemical reactions. The time required to reach background levels of nitric oxide concentration of 7 x 10 to the 7th power molecules/cc has been calculated. This relaxation time depends strongly on atmospheric conditions.

Eberstein, I. J.↗
NASA NTRS1972

A nonlinear analysis of pulsatile flow in arteries.

An approximate numerical method for calculating flow profiles in arteries is developed. The theory takes into account the nonlinear terms of the Navier-Stokes equations as well as the nonlinear behaviour and large deformations of the arterial wall. Through the locally measured values of the pressure, pressure gradient, and pressure-radius function, the velocity distribution and wall shear at a given location along the artery can be determined. The computed results agree well with the corresponding experimental data.

Ling, S. C.↗
NASA NTRS1973

Direct numerical solution of three-dimensional equations containing elliptic operators.

A direct three-dimensional elliptic solver is presented for application in a wide class of numerical methods for solving partial differential equations in physics and engineering. The derived algorithm and FORTRAN code implement Buzbee, Golub and Nielson's proposed extension of Buneman's Cyclic-Reduction Poisson solver to three dimensions. Both a 'most direct' cyclic reduction and a revised method (to eliminate roundoff error difficulties) are derived. Tests on an IBM 360/67 computer, using various optional combinations of subroutines, showed significant differences in accuracy and computing time, with the optimum subroutine combination depending on mesh size.

Martin, E. D.↗
NASA NTRS1973

Energy partitioning of gaseous ions in an electric field.

The partitioning of ion energy among thermal energy, drift energy, and random-field energy is studied by solution of the Boltzmann equation. An expansion in powers of the square of the electric field strength is obtained by Kihara's method. Numerical calculations for several ion-neutral force laws show that Wannier's constant mean-free-time model gives a reasonable first approximation. The formal extension to multicomponent mixtures is also given. The matrix elements obtained are tabulated, and can be used to study the field dependence of other moments of the ion-distribution function.

Hahn, H.-S.↗
NASA NTRS1972

Analysis and design of three dimensional supersonic nozzles. Volume 1: Nozzle-exhaust flow field analysis by a reference plane characteristics technique

A second order numerical method employing reference plane characteristics has been developed for the calculation of geometrically complex three dimensional nozzle-exhaust flow fields, heretofore uncalculable by existing methods. The nozzles may have irregular cross sections with swept throats and may be stacked in modules using the vehicle undersurface for additional expansion. The nozzles may have highly nonuniform entrance conditions, the medium considered being an equilibrium hydrogen-air mixture. The program calculates and carries along the underexpansion shock and contact as discrete discontinuity surfaces, for a nonuniform vehicle external flow.

Dash, S.↗
NASA NTRS1973

Mach Reflection, Mach Disc, and the Associated Nozzle Free Jet Flows

The numerical method involving both the method of integral relations and the method of characteristics have been applied to investigate the steady flow phenomena associated with the accurrence of Mach reflection and Mach disc from nozzle flows. The solutions of triple-shock intersection are presented. The regime where Mach configuration appears is defines for the inviscid analysis. The method of integral relations developed for the blunt body problem is modified and extended to the attached shock wave and to internal nozzle flow problems.

Chang, I.↗
NASA NTRS1974

Stress intensity factors for deep cracks emanating from the corner formed by a hole intersecting a plate surface

A technique consisting of a marriage between stress freezing photoelasticity and a numerical method was used to obtain stress intensity factors for natural cracks emanating from the corner at which a hole intersects a plate surface. Geometrics studied were: crack depth to thickness ratios of approximately 0.2, 0.5, and 0.75; crack depth to crack length ratios of approximately 1.0 to 2.0. All final crack geometries were grown under monotonic loading and growth was not self similar with most of the growth occurring through the thickness under remote extension. Stress intensity plate surface K sub s factors were determined at the intersection of the flaw border with the plate surface K sub s and with the edge of the hole K sub h. Results showed that for the relatively shallow flaws K sub h approximately equal to 1.5 K sub s, for the moderately deep flaws K sub h approximately equal to K sub s, and for the deep flaws K sub h approximately equal to 0.5 K sub s, revealing a severe sensitivity of K to flaw geometry.

Mcgowan, J. J.↗
NASA NTRS1973

Transverse cracks in a strip with reinforced surfaces

The symmetrical problem of two transverse cracks in an elastic strip with reinforced surfaces is formulated in terms of a singular integral equation. The special cases of one central crack or two edge cracks are discussed. Numerical methods for solving the problems with internal cracks are outlined and stress intensity factors are presented for various geometrics and degrees of surface reinforcement.

Krenk, S.↗
NASA NTRS1974

Radiative property investigation

A detailed formulation of radiative heat transfer problems involving the use of bidirectional reflectance is discussed. Except for very simple systems this formulation is very intricate. For this reason and the fact that bidirectional data are scarce, such a formulation is not in common use. However, with the development of the digital computer there have been numerical methods developed for detailed radiative investigations using the bidirectional reflectance. Since computations using bidirectional reflectance are coming into use for spacecraft radiative studies, it is necessary to have reflectance data on the materials involved. Materials that are frequently studied are coatings which are used in controlling the thermal environment of spacecrafts. This study is primarily concerned with the bidirectional reflectance of zinc oxide which is one constituent commonly used for coatings.

Scott, R. L., Jr.↗
NASA NTRS1974

Stress intensity factors for deep cracks emanating from the corner formed by a hole intersecting a plate surface

A technique consisting of freezing photo-elasticity and a numerical method was used to obtain stress intensity factors for natural cracks emanating from the corner at which a hole intersects a plate surface. Geometries studied were: (1) crack depth to thickness ratios of approximately 0.2, (2) 0.5 and 0.75; (3) crack depth to crack length ratios of approximately 1.0 to 2.0; and (4) crack length to hole radius ratios of about 0.5 to 2.0. All final crack geometries were grown under monotonic loading and growth was not self similar, with most of the growth occuring through the thickness under remote extension. Stress intensity factors were determined at the intersection of the flaw border.

Mcgowan, J. J.↗
NASA NTRS1974

Experimental determination of turbulence in a GH2-GOX rocket combustion chamber

The intensity of turbulence and the Lagrangian correlation coefficient for a gaseous rocket combustion chamber have been determined from the experimental measurements of the tracer gas diffusion. A combination of Taylor's turbulent diffusion theory and Spalding's numerical method for solving the conservation equations of fluid mechanics was used to calculate these quantities. Taylor's theory was extended to consider the inhomogeneity of the turbulence field in the axial direction of the combustion chamber. An exponential function was used to represent the Lagrangian correlation coefficient. The results indicate that the maximum value of the intensity of turbulence is about 15% and the Lagrangian correlation coefficient drops to about 0.12 in one inch of the chamber length.

Tou, P.↗
NASA NTRS1974

Lee surface flow phenomena over space shuttle at large angles of attack at M sub infinity equal 6

Surface pressure and heat transfer, flow separation, flow field, and oil flow patterns on the leeward side of a space shuttle orbiter model are investigated at a free stream Mach number of 6. The free stream Reynolds numbers are between 1.64 times 10 to the 7th power and 1.31 times 10 to the 8th power per meter, and the angle of attack is varied between 0 deg and 40 deg for the present experiments. The stagnation temperatures for the tests are approximately 500 K and the wall temperature is maintained at 290 K. Existing numerical methods of three-dimensional inviscid supersonic flow theory and compressible boundary layer theory are used to predict the present experimental measurements. Results of the present tests indicate two distinct types of flow separation and surface peak heating depending on the angle of attack.

Zakkay, V.↗
NASA NTRS1974

Analysis of dynamic stall using unsteady boundary-layer theory

The unsteady turbulent boundary layer and potential flow about a pitching airfoil are analyzed using numerical methods to determine the effect of pitch rate on the delay in forward movement of the rear flow reversal point. An explicit finite difference scheme is used to integrate the unsteady boundary layer equations, which are coupled at each instant of time to a fully unsteady and nonlinear potential flow analysis. A substantial delay in forward movement of the reversal point is demonstrated with increasing pitch rate, and it is shown that the delay results partly from the alleviation of the gradients in the potential flow, and partly from the effects of unsteadiness in the boundary layer itself. The predicted delay in flow-reversal onset, and its variation with pitch rate, are shown to be in reasonable agreement with experimental data relating to the delay in dynamic stall. From the comparisons it can be concluded (a) that the effects of time-dependence are sufficient to explain the failure of the boundary layer to separate during the dynamic overshoot, and (b) that there may be some link between forward movement of the reversal point and dynamic stall.

Scruggs, R. M.↗
NASA NTRS1973

Numerical analysis of stiffened shear webs in the postbuckling range

The postbuckling behavior of shear webs divided into rectangular panels by stiffeners (uprights) was studied numerically, using the STAGS program in which two-dimensional finite differences are used to solve buckling and nonlinear problems. Universal nondimensional parameters, suggested by linear buckling data, are found to hold for the postbuckling range. Results indicate that the postbuckling stiffness of shear webs with isotropic panels is roughly two thirds of the prebuckling stiffness. The postbuckling behavior of shear webs with isotropic and orthotropic material properties is compared. Some practical aspects of solving nonlinear problems of this type by the numerical method employed are examined.

Stein, M.↗
NASA NTRS1973

Computational fracture mechanics

Numerical procedures for accurate determination of elastic stress intensity factors for the general two-dimensional crack problem are reviewed. The elastic perfectly plastic state of crack tip deformation is studied by a finite element procedure. Elastic-plastic fields in the immediate vicinity of a crack tip are determined numerically by finite element procedures based on asymptotic studies of crack tip singularities in plastic materials. The small-scale yielding problem is modeled, and expressions for crack tip opening displacement, shear singularity amplitude, and plastic zone extent are derived. A finite element solution to the large-scale yielding of a circumferentially cracked round tension bar is obtained. The three-dimensional aspects of flawed structures and numerical methods of treating them are studied. Ductile fracture mechanisms, in particular crack tip fracture on the microscale, are discussed.

Rice, J. R.↗
NASA NTRS1974

Analysis of flow-reversal delay for a pitching airfoil

The coupled unsteady turbulent boundary layer and potential flow about a pitching airfoil are analyzed using numerical methods to determine the effect of pitch rate on the delay in forward movement of the rear flow reversal point. An explicit finite difference scheme is used to integrate the unsteady boundary layer equations, which are coupled at each instant of time to a fully unsteady and nonlinear potential flow analysis. A substantial delay in forward movement of the reversal point is demonstrated with increasing pitch rate, and it is shown that the delay results partly from the alleviation of the gradients in the potential flow, and partly from the effects of unsteadiness in the boundary layer itself.

Scruggs, R. M.↗