SEARCH · Engineering Papers
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.
Programs manual for one-dimensional and two- dimensional steady state models of catalyzed hydrazine decomposition reaction chambers
An approach for the solution of transient-forced convection film boiling on an isothermal flat plate using the boundary layer model is developed. The similarity variables are used to convert the governing partial differential equations to ordinary ones. The results of numerical solutions of these ordinary equations indicate that the transient process can be classified as one-dimensional conduction, intermediate, and the steady state regions. The time required for the one-dimensional conduction and the time necessary to attain a steady state condition are obtained. The influence of interfacial shear is seen to be negligible while the Prandtl Number and the ratio (C sub p delta T divided by h sub fg times Pr) have major influence. The use of local similarity approximations for the intermediate regime facilitates prediction of complete boundary layer growth. Using the ratio of time at any instant to the steady state time as abscissa, the curves representing the boundary layer growth can be merged into a single mean curve within 5 percent. Further, the analysis shows that the average rate of heat transfer during transient is 50 to 100 percent higher than those at steady state. The average rate of vapor convected away is 10 to 15 percent lower than at steady state while the average rate of accumulation to form the vapor layer is 1 to 14 times larger. Further, the total heat transferred during transient increases and the evaporation decreases for increasing values of C sub p delta T divided by h sub fg times Pr
Development of a new approach for the solution of transient-forced convection film boiling on an isothermal flat plate using the boundary layer model. The similarity variables are used to convert the governing partial differential equations to ordinary ones. The results of numerical solutions of these ordinary equations indicate that the transient process can be classified as one-dimensional conduction, intermediate, and the steady-state regions. The time required for the one-dimensional conduction and the time necessary to attain a steady-state condition are obtained. The use of local similarity approximations for the intermediate regime facilitates prediction of complete boundary layer growth. Using the ratio of time at any instant to the steady-state time as abscissa, the curves representing the boundary layer growth can be merged into a single mean curve within 5%. Further, the analysis shows that the average rate of heat transfer during transient is 50 to 100% higher than those at steady state. The average rate of vapor convected away is 10 to 15% lower than at steady state while the average rate of accumulation to form the vapor layer is 1 to 14 times larger.
One-dimensional systems of N = 500 and 250 particles in equilibrium are numerically simulated utilizing the method of molecular dynamics. Periodic boundary conditions are imposed. The classical two-body interaction potential is short range, repulsive and has a corresponding finite force. The equations of state are determined for densities both less and greater than one. Corresponding theoretical isochores are determined from models based on nearest-neighbor interactions and on a truncated virial expansion, and a comparison is made with the experimental isochores. Time independent radial distributions are constructed numerically and discussed. A change of state from a solidlike state to a fluid-gas state based on the penetrability of the particles is predicted. The transition temperatures are estimated from the radial distribution functions and the nearest-neighbor model. Self-diffusion is observed and the corresponding constants are determined from the velocity autocorrelation functions.
A finite element algorithm for solving unsteady, three-dimensional high velocity impact problems is presented. A computer program was developed based on the Eulerian hydroelasto-viscoplastic formulation and the utilization of the theorem of weak solutions. The equations solved consist of conservation of mass, momentum, and energy, equation of state, and appropriate constitutive equations. The solution technique is a time-dependent finite element analysis utilizing three-dimensional isoparametric elements, in conjunction with a generalized two-step time integration scheme. The developed code was demonstrated by solving one-dimensional as well as three-dimensional impact problems for both the inviscid hydrodynamic model and the hydroelasto-viscoplastic model.
The computer simulation of the electrical behavior of a Gate-Controlled Diode (GCD) fabricated in Silicon-On-Sapphire (SOS) was described. A procedure for determining lifetime profiles from capacitance and reverse current measurements on the GCD was established. Chapter 1 discusses the SOS structure and points out the need of lifetime profiles to assist in device design for GCD's and bipolar transistors. Chapter 2 presents the one-dimensional analytical formula for electrostatic analysis of the SOS-GCD which are useful for data interpretation and setting boundary conditions on a simplified two-dimensional analysis. Chapter 3 gives the results of a two-dimensional analysis which treats the field as one-dimensional until the silicon film is depleted and the field penetrates the sapphire substrate. Chapter 4 describes a more complete two-dimensional model and gives results of programs implementing the model.
A nonlinear analysis is developed for sound propagation in a variable area duct in which the mean flow approaches choking conditions. A quasi-one-dimensional model is used; results of the standard linear theory are compared with the nonlinear results to assess the significance of the nonlinear terms. The nonlinear analysis represents the acoustic disturbance as a sum of interacting harmonics. Numerical results show that the basic signal is unaffected by the presence of higher harmonics if the throat Mach number is not too large, but as the Mach number approaches unity more harmonics are needed to describe the acoustic propagation. The strong interactions among harmonics in the numerical results occur in a region which is generally consistent with the nonlinear inner-expansion region of Callegari and Myers.
Results are presented from an experimental investigation of perpendicular, hydrogen fuel injection and combustion from opposing walls in a scramjet combustor model using a longitudinally staged laterally inline step-injection configuration. The model represents a portion of the flow in the Langley integrated modular scramjet engine combustor operating at a flight Mach number of 7. When operating at a ratio of jet pressure to free-stream dynamic pressure of 3, the injectors produce a bulk equivalence ratio of unity. This investigation represents part of a continuing study of the modular engine fuel injectors and is specifically designed to eliminate the adverse lateral pressure gradient observed at the injector location in a previous test. Flow survey contours at three axial locations, ranging from one-third of the engine combustor length to the total engine combustor length, are used to determine mixing efficiency and fuel distribution. Wall static pressures are analyzed by using one-dimensional theory to determine the combustion efficiency. Results show a significant improvement over previous injector designs tested in this duct geometry.
The problem of flutter of a tensioned elastic plate is reconsidered on a model of two-dimensional aerodynamic flow over a one-dimensional, infinitely wide plate. The purpose is to clarify the effect of variations in plate bending stiffness on the flutter of a plate under a given tension, and vice versa. On the basis that the flutter mode for the one-dimensional infinitely wide plate is essentially the first plate natural mode for subsonic and low supersonic Mach numbers, and that the dynamic pressure required for flutter is proportional to the square of the first natural mode frequency, it is shown that the membrane paradox disappears for Mach numbers 0 and 1.3, so that reduced bending stiffness reduces the dynamic pressure required for flutter.
An optimal solution of the ascent trajectory of the Space Shuttle for the first orbital flight test is presented; the optimization is a minimum propellant, four-control problem in yaw angle, roll angle, pitch angle and vacuum thrust of each Space Shuttle main engine. Piecewise linear segments with juncture points treated as parameters are employed to model the controls. Equations of motion for a three-dimensional flight with pitch plane moment balance about an oblate are integrated numerically with a fourth-order Runge-Kutta method; two- and one-dimensional cubic spline function curve fits of aerodynamic coefficients are used during the first and second stages, respectively. The constraint minimization problem is solved with the Davidon-Fletcher-Powell function method.
This paper deals with general inferences about the low-density phase of interstellar gas having temperatures well above 200,000 K (the 'coronal gas') which can be drawn from O VI absorption data for 72 stars. Attention is given to the behavior of radial velocities, possible evidence for circumstellar O VI, the space distribution of the O VI gas, crowding of normal interstellar gas, and temperature distributions for the coronal gas. A model is adopted in which the coronal gas is contained within randomly distributed and nonoverlapping parcels, each with a size, pressure, and internal temperature distribution that do not vary markedly from one unit to the next. It is shown that the one-dimensional velocity dispersion for O VI regions (26 km/s) is substantially higher than the value for ordinary interstellar clouds (6.4 km/s).
Critical to the advancement of computational aerodynamics capability is the ability to simulate flows about three-dimensional configurations that contain both compressible and viscous effects, including turbulence and flow separation at high Reynolds numbers. Analyses were conducted of two solution techniques for solving the Reynolds averaged Navier-Stokes equations describing the mean motion of a turbulent flow with certain terms involving the transport of turbulent momentum and energy modeled by auxiliary equations. The first solution technique is an implicit approximate factorization finite-difference scheme applied to three-dimensional flows that avoids the restrictive stability conditions when small grid spacing is used. The approximate factorization reduces the solution process to a sequence of three one-dimensional problems with easily inverted matrices. The second technique is a hybrid explicit/implicit finite-difference scheme which is also factored and applied to three-dimensional flows. Both methods are applicable to problems with highly distorted grids and a variety of boundary conditions and turbulence models.
Several FORTRAN program modules for performing one-dimensional and two-dimensional discrete Fourier transforms, Mellin, and Fourier-Bessel transforms are described along with programs that realize the algebra of high speed Fourier transforms on a computer. The programs can perform numerical harmonic analysis of functions, synthesize complex optical filters on a computer, and model holographic image processing methods.
Possible applications of the principle of similarity to fluid mechanics is described and illustrated. In correlating thermophysical properties of fluids, the similarity principle transcends the traditional corresponding states principle. In fluid mechanics the similarity principle is useful in correlating flow processes that can be modeled adequately with one independent variable (i.e., one-dimensional flows). In this paper we explore the concept of transforming the conservation equations by combining similarity principles for thermophysical properties with those for fluid flow. We illustrate the usefulness of the procedure by applying such a transformation to calculate two phase critical mass flow through a nozzle.
A survey of current concepts regarding electrostatic ion-cyclotron turbulence (theory and experiment), and regarding inclined double potential layers in the magnetospheric plasma is presented. Anomalous resistance governed by electrostatic ion-cyclotron turbulence, and one-dimensional and two-dimensional models of double electrostatic layers in the magnetospheric plasma are examined.
An approximate method to determine the index of refraction of a dielectric obstacle is presented. For simplicity one dimensional models of electromagnetic scattering are treated. The governing equations yield a second order boundary value problem, in which the index of refraction appears as a functional parameter. The availability of reflection coefficients yield two additional boundary conditions. The index of refraction by a k-th order spline which can be written as a linear combination of B-splines is approximated. For N distinct reflection coefficients, the resulting N boundary value problems yield a system of N nonlinear equations in N unknowns which are the coefficients of the B-splines.
Solute redistribution during directional solidification of HgCdTe is addressed. Both one-dimensional and two-dimensional models for solute redistribution are treated and model results compared to experiment. The central problem studied is the cause of radial inhomogeneities found in directionally solidified HgCdTe. A large scale gravity-driven interface instability, termed shape instability, is postulated to be the cause of radial inhomogeneities. Recommendations for future work, along with appropriate computer programs, are included.
An approximate method to determine the index of refraction of a dielectric obstacle is presented. For simplicity one dimensional models of electromagnetic scattering are treated. The governing equations yield a second order boundary value problem, in which the index of refraction appears as a functional parameter. The availability of reflection coefficients yield two additional boundary conditions. The index of refraction by a k-th order spline which can be written as a linear combination of B-splines is approximated. For N distinct reflection coefficients, the resulting N boundary value problems yield a system of N nonlinear equations in N unknowns which are the coefficients of the B-splines.