Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical method”

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 1,405 records · Page 78

Three dimensional nozzle-exhaust flow field analysis by a reference plane technique.

A numerical method based on reference plane characteristics has been developed for the calculation of highly complex supersonic nozzle-exhaust flow fields. The difference equations have been developed for three coordinate systems. Local reference plane orientations are employed using the three coordinate systems concurrently thus catering to a wide class of flow geometries. Discontinuities such as the underexpansion shock and contact surfaces are computed explicitly for nonuniform vehicle external flows. The nozzles considered may have irregular cross-sections with swept throats and may be stacked in modules using the vehicle undersurface for additional expansion. Results are presented for several nozzle configurations.

Dash, S. M.↗

Buckling and vibration analysis for stiffened orthotropic shells of revolution.

Development of a numerical method, using the multistep integration approach, in which the buckling and vibration analyses are formulated as a succession of linear eigenvalue problems. The method does not require an estimate of the eigenvector, and once an eigenvalue has been converged good estimates for other eigenvalues are automatically available. This is accomplished through the use of an in-core Householder scheme for solution of the eigenvalue problem. Furthermore, since the method uses an eigenvalue solution, the possibility of missing modes is eliminated.

Svalbonas, V.↗

The radiative transfer equation and environmental effects in the upper atmosphere.

A highly effective approach to the solution of the radiative transfer equation, including the effects of emission, absorption, and multiply scattered radiation, is described. The transfer equation is cast in its unsteady hyperbolic form which is amenable to several different numerical methods of solution. Comparisons of solutions with emission, absorption, and scattering of radiation, are presented to demonstrate the validity of the present approach.

Callis, L. B.↗

Aeroelastic optimization of a panel in high Mach number supersonic flow.

Solution for a least-weight skin thickness distribution for a panel with a flutter parameter constraint. This panel weighs less than any similar constant thickness panel, but has the same critical supersonic panel flutter parameter. The panel rests on simple supports and is of sandwich construction. The span to chord ratio is large enough that the inertial, elastic, and aerodynamic behavior is one-dimensional. The Mach number is great enough that the aerodynamic forces acting on the upper panel surface may be accurately described by quasi-steady, linearized, supersonic aerodynamic theory. The final optimum design is obtained from theoretical and numerical methods adapted from optimal control theory. The results of this investigation show that the optimal panel thickness distribution is symmetric about the panel chord midpoint. Compared to a reference panel with constant thickness, optimum panels are found to be nearly 12% lighter.

Weisshaar, T. A.↗

Saturn systems holddown acoustic efficiency and normalized acoustic power spectrum.

Saturn systems field acoustic data are used to derive mid- and far-field prediction parameters for rocket engine noise. The data were obtained during Saturn vehicle launches at the Kennedy Space Center. The data base is a sorted set of acoustic data measured during the period 1961 through 1971 for Saturn system launches SA-1 through AS-509. The model assumes hemispherical radiation from a simple source located at the intersection of the longitudinal axis of each booster and the engine exit plane. The model parameters are evaluated only during vehicle holddown. The acoustic normalized power spectrum and efficiency for each system are isolated as a composite from the data using linear numerical methods. The specific definitions of each allows separation. The resulting power spectra are nondimensionalized as a function of rocket engine parameters. The nondimensional Saturn system acoustic spectrum and efficiencies are compared as a function of Strouhal number with power spectra from other systems.

Gilbert, D. W.↗

Theoretical and experimental investigations of antennas and waves in plasma

The near field characteristics of an antenna when it is covered by a plasma sheath were investigated. The areas of investigation were: (1) current distribution, (2) input admittance, and (3) radiation pattern as functions of the electrical length of the antenna. Diagnostic measurements were also made to determine the properties of the plasma. A numerical method was evolved to solve the problem of a finite cylindrical antenna in an infinite plasma column.

Source record↗

Numerical calculations of flow fields

The solutions to the equations of motion for inviscid fluid flow around a pointed elliptic cone at incidence are presented. The numerical method used, MacCormack's second order preferential predictor-corrector finite difference approximation, is applied to the fluid flow equations derived in conservation-law form. The entropy boundary condition, hitherto unused for elliptic cone problems, is investigated and compared to reflection boundary condition solutions. The stagnation streamline movement of the inclined elliptic cone is noted and surface pressure coefficients are plotted. Also presented are solutions for an elliptic cone and a circular cone at zero incidence and a circular cone at a small angle of attack. Comparisons are made between these present solutions and previously published theory.

Anderson, D. M.↗

Tidal perturbations on the orbits of GEOS-1 and GEOS-2

Analysis of the luni-solar tidal perturbations of the inclination of GEOS-1 (1965-89A) and GEOS-2 (1968-002A) yielded the values k2 = 0.22 (sigma = 0.02) and 0.31 (sigma = 0.01) respectively for the second degree Love number. For GEOS-1 a new, purely numerical method involving osculating elements was employed. For GEOS-2 it was necessary to analyze the variations of the mean elements because of the very long period (450d) of the dominant solar tidal perturbation. An additional analysis of the variation of the mean elements of GEOS-1 confirmed the value of k2 obtained from the osculating elements.

Douglas, B. C.↗

Performance Analysis and Design Synthesis (PADS) computer program. Volume 2: Program description, part 2

The QL module of the Performance Analysis and Design Synthesis (PADS) computer program is described. Execution of this module is initiated when and if subroutine PADSI calls subroutine GROPE. Subroutine GROPE controls the high level logical flow of the QL module. The purpose of the module is to determine a trajectory that satisfies the necessary variational conditions for optimal performance. The module achieves this by solving a nonlinear multi-point boundary value problem. The numerical method employed is described. It is an iterative technique that converges quadratically when it does converge. The three basic steps of the module are: (1) initialization, (2) iteration, and (3) culmination. For Volume 1 see N73-13199.

Source record↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗