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.
A three-dimensional, time-dependent mathematical model to represent Mobile Bay was developed. Computer programs were developed which numerically solve the appropriate conservation equations for predicting bay and estuary flow fields. The model is useful for analyzing the dispersion of sea water into fresh water and the transport of sediment, and for relating field and physical model data.
Conservation equations are developed in the three-fluid approximation for general application problems of ionospheric dynamics in the altitude region 90 km to 800 km for all geographic locations. These equations are applied to a detailed study of auroral E region neutral winds and their relationship to ionospheric plasma motions.
An analytical study of nonsimilar jet mixing is made for compressible, nonisoenergetic flows. The conservation equations are solved for each of the streams above and below the dividing streamline by using Meksyn's asymptotic method of integration for solving boundary-layer problems. The problem of laminar mixing between two parallel streams is investigated for the case of a constant pressure gradient. It is found that the velocity and temperature profiles from the exact solution to the nonsimilar governing equations can be well approximated by the locally similar solution.
Wave propagation in a quasi-chemical equilibrium plasma is studied. The plasma is infinite and without external fields. The chemical reactions are assumed to result from the ionization and recombination processes. When the gas is near equilibrium, the dominant role describing the evolution of a reacting plasma is played by the global conservation equations. These equations are first derived and then used to study the small amplitude wave motion for a near-equilibrium situation. Nontrivial damping effects have been obtained by including the conduction current terms.
A theoretical investigation was undertaken to develop a relatively simple model of a two-dimensional, turbulent wall jet in a coflowing stream. The incompressible jet flow was modeled by using an integral method which includes turbulent shear stress, entrainment, and heat transfer. The method solves the conservation equations for the average jet flow properties and uses the velocity profile suggested by Escudier and Nicoll to obtain detailed characteristics of the jet on a flat plate. The analytical results compare favorably with experimental data for a range of injection velocities, which demonstrates the usefulness of the theory for estimating jet growth, velocity decay, and wall skin friction. The theory, which was applied to a Coanda jet on a circular cylinder, provided estimates of suction pressures aft of the jet exit that were in close agreement with experimental values.
The intensity of turbulence and the Lagrangian correlation coefficient for a LOX-GH2 rocket combustion chamber was determined from experimental measurements of tracer gas diffusion. A combination of Taylor's turbulent diffusion theory and a 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, and an exponential function was used to represent the Lagrangian correlation coefficient. The results indicate that the value of the intensity of turbulence reaches a maximum of 14% at a location about 7" downstream from the injector. The Lagrangian correlation coefficient associated with this value is given by the above exponential expression where alpha = 10,000/sec.
The results of a theoretical investigation of thermodynamic effects on developed cavitation are presented. An approximate solution to the conservation equations for a two-phase laminar boundary layer is obtained. This analysis produces an expression for the temperature difference between the liquid and vapor phases which can be applied to developed cavity flows. Experimental data of cavity temperature depressions are correlated using this result. In addition, a theoretical estimate of the Nusselt number for the cavity is made using a turbulent boundary layer cavity model proposed by Brennen. The result agrees in part with empirically determined expressions for the cavity Nusselt number.
The nonlinear growth and limiting amplitude of acoustic waves in a combustion chamber are considered. A formal framework is provided within which practical problems can be treated with a minimum of effort and expense. The general conservation equations were expanded in two small parameters, one characterizing the mean flow field and one measuring the amplitude of oscillations, and then combined to yield a nonlinear inhomogeneous wave equation. The unsteady pressure and velocity fields were expressed as syntheses of the normal modes of the chamber, but with unknown time-varying amplitudes. This procedure yielded a representation of a general unsteady field as a system of coupled nonlinear oscillators. The system of nonlinear equations was treated by the method of averaging to produce a set of coupled nonlinear first order differential equations for the amplitudes and phases of the modes. The analysis is applicable to any combustion chamber. The most interesting applications are probably to solid rockets, liquid rockets, or thrust augmentors on jet engines.
Approximate analytical expressions are derived for the velocity and temperature distributions in steady state coal slag deposits flowing over MHD generator walls. Effects of slag condensation and Joule heating are included in the analysis. The transport conditions and the slag temperature at the slag-gas interface are taken to be known parameters in the formulation. They are assumed to have been predetermined either experimentally or from the slag properties and the gas dynamic calculations of the free stream flow. The analysis assumes a power law velocity profile for the slag and accounts for the coupling between the energy and momentum conservation equations. Comparisons are made with the more exact numerical solutions to verify the accuracy of the results.
The underlying bases and developments in two techniques of detailed turbulence modeling are described where the flow is treated in the Eulerian sense, and one technique where the Lagrangian motions of vortices are followed. First, a technique is described for solving the single-point statistically averaged conservation equations. The Reynolds stresses that appear in these equations are evaluated by solving supplemental differential equations which contain terms that are modeled. A sequence of increasingly complex, but also increasingly general, modeling equations is described and computations based on these equations are compared with experimental data. The hierarchy of models described terminates with equations for the individual components of the Reynolds stress tensor. The second Eulerian technique approach to turbulence modeling is the direct numerical simulation of turbulent fields. In this approach, all three dimensional eddies between a predetermined range of sizes are computed in time within a specified volume of flow. Present day computers require a tradeoff between the size of the volume that can be considered and the degree of resolution of the turbulent eddies. Techniques of modeling the smallest eddies are described that permit enlarging the volume, or Reynolds number, that can be considered.
A semi-empirical fluid mechanical model of the acoustic behavior of Helmholtz resonators is presented which predicts impedance as a function of the amplitude and frequency of the incident sound pressure field and resonator geometry. The model assumes that the particle velocity approaches the orifice in a spherical manner. The incident and cavity sound fields are connected by solving the governing oscillating mass and momentum conservation equations. The model is in agreement with the Rayleigh slug-mass model at low values of incident sound pressure level. At high values, resistance is predicted to be independent of frequency, proportional to the square root of the amplitude of the incident sound pressure field, and virtually independent of resonator geometry. Reactance is predicted to depend in a very complicated way upon resonator geometry, incident sound pressure level, and frequency. Nondimensional parameters are defined that divide resonator impedance into three categories corresponding to low, moderately low, and intense incident sound pressure amplitudes. The two-microphone method was used to measure the impedance of a variety of resonators. The data were used to refine and verify the model.
The general problem of the nonlinear growth and limiting amplitude of acoustic waves in a combustion chamber is treated in three parts: (1) the general conservation equations are expanded in two small parameters, and then combined to yield a nonlinear inhomogeneous wave equation, (2) the unsteady pressure and velocity fields are expressed as a synthesis of the normal modes of the chamber, but with unknown time-varying amplitudes, and (3) the system of nonlinear equations is treated by the method of averaging to produce a set of coupled nonlinear first order differential equations for the amplitudes and phases of the modes. This approximate analysis is applied to the investigation of the unstable motions in a solid propellant rocket engine and in a T burner.
A numerical program was developed to compute transient laminar flows in two dimensions including multicomponent mixing and chemical reaction. The program can compute both incompressible flows and compressible flows at all speeds, and it is applied to describe transient and steady state solutions for low subsonic, coaxial entry, tue flows. Single component, nonreacting flows comprise most of the solutions, but one steady state solution is presented for trace concentration constituents engaging in a second order reaction. Numerical stability was obtained by adding at each calculation point a correction for numerical diffusion errors caused by truncation of the Taylor series used to finite difference the conservation equations. Transient computations were made for fluids initially at rest, then subjected to step velocity inputs that were uniform across each region of the entry plane and were held constant throughout the computation period. For center tube to annulus velocity ratios of 0.5 and 2.0, the bulk fluid in the tube initially moved in plug flow, but strong radial flows developed near the injection plane which moved the fluid into the high shear region between the jets and away from the tube wall.
Earlier studies have shown that interplanetary disturbances observed by spacecraft at large heliocentric distances from the sun can be directly associated with various solar activities such as surges, sprays, eruptive prominences, and flares. A digital computer code to simulate these phenomena is developed. This numerical code solves a set of conservation equations for an adiabatic time-dependent spherically-symmetric fluid flow. The solar events are simulated by perturbations in the dependent variables at the lower boundary of the quiet sun corona, and the time-dependent numerical solution simulates the coronal response. Three physical cases are simulated by using the proposed numerical code: (1) part of a surge, spray and prominence; (2) a single flare; and (3) two flares whose interplanetary consequences interact in space. Results for the first two cases are compared with both ground-based and space probe observations. Fairly good agreement with observations is found. A comparison of the third case with observations is presently under way.
The response of a lunar-sized object to the impact of meteoroids no more than about 100 km in radius is studied by means of a numerical model. The partitioning of impact energy into the kinetic and internal energy of the ejecta is obtained by using the conservation of mass, momentum, and energy conservation equations in finite-difference form within an Eulerian framework with approximate equations of state. The calculations are performed for a 15 km/sec impact of an iron object 5 cm in radius on a gabbroic anorthosite surface. Ejecta ballistic analysis is then performed. Most of the material lost escaping the moon is lunar crust material. Only 0.2% of the meteoroid escapes, all in the vapor phase.
Photoionization of the upper atmosphere of Titan by sunlight is expected to produce a substantial ionospheric layer. One-dimensional forms of the mass, momentum, and energy conservation equations for ions and electrons have been solved along with electron number densities of about 1000/cu cm, using various model atmospheres. The significant ions in a CH4-H2 atmosphere are H(+), H3(+), CH5(+), CH3(+), and C2H5(+). Electron temperatures may be as high as 1000 K, depending on the abundance of hydrogen in the high atmosphere. Interaction of the solar wind with the ionosphere is also discussed.
In their most general form, the Reynolds averaged conservation equations result from ensemble or time averages of the instantaneous Navier-Stokes equations or their compressible counterparts. For these averaging processes to be consistent, the averaging time period must exceed the periods identified with the largest time scales of the turbulence, and yet be shorter than the characteristic times of the flow field. With these equations long period variations in the flow fields are deterministic, provided initial conditions are known. The average dependent variables are sufficiently smooth to be resolvable by finite difference techniques consistent with the size and speed of modern computers.
Flame speed and mean degree of reaction profiles are determined for an ensemble mean one-dimensional premixed turbulent flame. Elementary closure models are employed to treat turbulence closure problems arising from the nonlinear convective and chemical reaction rate terms in the time-dependent fuel conservation equation; a constant eddy diffusivity is useed to model correlation terms derived from the convective term. An expression developed for fuel consumption is solved to obtain the turbulent flame speed and the variation of mean fuel mass fraction across the flame. Comparison of the numerical analysis with experimental data indicates the importance of the quenching of the combustion reaction in high-wavenumber eddies.