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 problem is considered which arises in the numerical analysis of the degenerate partial differential equations of stochastic control theory. The relationship between diffusions and finite state Markov chains is considered. The role of probability theory in numerical analysis is illustrated.
The nonlinear evolution of unstable sound waves in a uranium plasma has been calculated using a multiple time-scale asymptotic expansion scheme. The fluid equations used include the fission power density, radiation diffusion, and the effects of the changing degree of ionization of the uranium atoms. The nonlinear growth of unstable waves is shown to be limited by mode coupling to shorter wavelength waves which are damped by radiation diffusion. This mechanism limits the wave pressure fluctuations to values of order delta P/P approximates 0.00001 in the plasma of a typical gas-core nuclear rocket engine. The instability is thus not expected to present a control problem for this engine.
The nonlinear evolution of unstable sound waves in a uranium plasma has been calculated using a multiple time-scale asymptotic expansion scheme. The fluid equations used include the fission power density, radiation diffusion, and the effects of the changing degree of ionization of the uranium atoms. The nonlinear growth of unstable waves is shown to be limited by mode coupling to shorter wavelength waves which are damped by radiation diffusion. This mechanism limits the wave pressure fluctuations to values of order delta P/P equal to about .00001 in the plasma of a typical gas-core nuclear rocket engine. The instability is thus not expected to present a control problem for this engine.
The equations giving the response of a slightly ionized plasma with monatomic components to sinusoidal perturbations have been formulated. Included in the model equations were the electron Hall effect, electron thermal diffusion, radiation, and electron-atom rate processes. Plasma conditions were limited to those where viscous effects, the induced magnetic field, ion slip, and atom-atom inelastic processes can be neglected. Presented are results of numerical calculations for MHD generators with a working fluid of potassium seeded argon.
The method for calculating the pitch angle diffusion coefficient is investigated for cosmic rays in a static random magnetic field, using the resonance integral method. The pitch angle diffusion coefficient may be derived from the Vlasov equation via ensemble averaging.
A FORTRAN program for calculating the off-design performance of centrifugal compressors with channel diffusers is presented. Use of the program requires complete knowledge of the overall impeller and diffuser geometries. Individual losses are computed using analytical equations and empirical correlations which relate loss levels to velocity diagram characteristics and overall geometry. On a given speed line compressor performance is calculated for a range of inlet velocity levels. At flow rates between surge and choke, individual efficiency decrements, compressor overall efficiency, and compressor total pressure ratio are tabulated. An example case of performance comparison with a compressor built by a commercial engine manufacturer is presented to demonstrate the correlation with limited experimental data.
When the distribution function for energetic particles moving in interplanetary space is expressed in terms of eigenfunctions of an operator which describes pitch angle scattering by random magnetic fields, the familiar phenomenon of diffusion is the dominant solution of the transport equations provided that the spectrum of magnetic fluctuations is not too steep. The present paper describes a novel mode of particle transport when the latter condition is not satisfied (i.e., when the exponent for the spectrum of magnetic fluctuations is greater than 2). In this mode, the density inhomogeneities propagate coherently, in wavelike fashion, and with very little dispersion along the guiding field at a characteristic velocity equal to half the particle speed. The most striking manifestation of this phenomenon is the observed scatter-free propagation of solar flare electrons below 1 MeV.-
When the distribution function for energetic particles moving in interplanetary space is expressed in terms of eigenfunctions of an operator which describes pitch-angle scattering by random magnetic fields, the familiar phenomenon of diffusion is the dominant solution of the transport equations provided that the spectrum of magnetic fluctuations is not too steep. However, if the power-versus-wavenumber spectrum of the random fields has a spectral index greater than 2, the dominant solution of the same equations is a qualitatively different mode of transport in which density inhomogeneities propagate coherently along the guiding field with a characteristic velocity equal to half the particle speed. The present work, which extends classical transport theory to give a detailed treatment of the new mode, presents explicit formulae which describe the spatial and temporal structure of the coherent disturbance. It is shown that this disturbance evolves asymptotically into a moving Gaussian pulse whose width increases with time.
The present work describes some multilayer diffusion models and a computer program for these models developed to predict the impact of ground clouds formed during Space Shuttle launches on ambient air quality. The diffusion models are based on the Gaussian plume equation for an instantaneous volume source. Cloud growth is estimated on the basis of measurable meteorological parameters: standard deviation of the wind azimuth angle, standard deviation of wind elevation angle, vertical wind-speed shear, vertical wind-direction shear, and depth of the surface mixing layer. Calculations using these models indicate that Space Shuttle launches under a variety of meteorological regimes at Kennedy Space Center and Vandenberg AFB are unlikely to endanger the exposure standards for HCl; similar results have been obtained for CO and Al2O3. However, the possibility that precipitation scavenging of the ground cloud might result in an acidic rain that could damage vegetation has not been investigated.
The hydrodynamic properties and the reaction kinetics of axisymmetric coaxial turbulent jets having steady mean quantities are investigated. From the analysis, limited to free turbulent boundary layer mixing of such jets, it is found that the two-equation model of turbulence is adequate for most nonreactive flows. For the reactive flows, where an allowance must be made for second order correlations of concentration fluctuations in the finite rate chemistry for initially inhomogeneous mixture, an equation similar to the concentration fluctuation equation of a related model is suggested. For diffusion limited reactions, the eddy breakup model based on concentration fluctuations is found satisfactory and simple to use. The theoretical results obtained from these various models are compared with some of the available experimental data.
The paper examines the mathematical and physical foundations for the kinetic theory of reactive turbulent flows, discussing the differences and relation between the kinetic and averaged equations, and comparing some solutions of the kinetic equations obtained by the Green's function method with those obtained by the approximate bimodal method. The kinetic method described consists essentially in constructing the probability density functions of the chemical species on the basis of solutions of the Langevin stochastic equation for the influence of eddies on the behavior of fluid elements. When the kinetic equations are solved for the structure of the diffusion flame established in a shear layer by the bimodal method, discontinuities in gradients of the mean concentrations at the two flame edges appear. This is a consequence of the bimodal approximation of all distribution functions by two dissimilar half-Maxwellian functions, which is a very crude approximation. These discontinuities do not appear when the solutions are constructed by the Green's function method described here.
The Reynolds stress equation, pressure strain correlation, and dissipative terms and diffusion are discussed in relation to turbulence modelling using the Reynolds stress model. Algebraic modeling of Reynolds stresses and calculation of the boundary layer over an axial cylinder are examined with regards to the kinetic energy model for turbulence modelling. The numerical analysis of blade and hub wall boundary layers, and an experimental study of rotor blade boundary layer in an axial flow compressor rotor are discussed. The Patankar-Spalding numerical method for two dimensional boundary layers is included.
To study the nonlinear physics of inhomogeneous turbulent shear flow, the unaveraged Navier-Stokes equations are solved numerically. For initial conditions a three-dimensional cosine velocity fluctuation and a mean-velocity profile with a step are used. Although the initial conditions are nonrandom. The flow soon becomes turbulent. Concentrated turbulent energy develops near the plane where the mean velocity gradient is initially infinite. The terms in the one-point correlation equation for turbulent energy, including those for the diffusion and production of turbulence, are calculated, the diffusion terms tend to make the turbulence more homogeneous.
This paper describes numerical simulations of self-excited oscillations in a two-dimensional transonic inlet-diffuser flow by solving the Navier-Stokes equations with a two-equation turbulence model. The calculated amplitudes of oscillations for the terminal shock and the velocity fields compare well with experimental measurements; however, the predicted frequency of oscillations is about 50 percent higher. The formation of a pair of downstream-traveling, counter-rotating vortices at each cycle of velocity fluctuations, as reported experimentally, is vividly revealed by the numerical results.
Regarding new bipolar magnetic regions as sources of flux, the evolution of the radial component of the solar photospheric magnetic field is simulated during 1976-1984, and the corresponding evolution of the line-of-sight polar fields as seen from earth is derived. The observed timing and strength of the polar-field reversal during cycle 21 can be accounted for by supergranular diffusion alone, for a diffusion coefficient of 800 sq km/sec. For an assumed 300 sq km/sec rate of diffusion, on the other hand, a poleward meridional flow with a moderately broad profile and a peak speed of 10 m/s reached at about 5 deg latitude is required to obtain agreement between the simulated and observed fields. Such a flow accelerates the transport of following-polarity flux to the polar caps, but also inhibits the diffusion of leading-polarity flux across the equator. For flows faster than about 10 m/s, the latter effect dominates, and the simulated polar fields reverse increasingly later and more weakly than the observed fields.
This note discusses the consequences of the interaction of radiatively driven mass loss with magnetic fields in upper-main-sequence helium-peculiar stars. The field constrains mass loss to occur in narrow cones above the magnetic poles. These jets can explain the observed behavior of UV resonance lines in the helium-peculiar stars, especially the sn and helium-rich stars. Suppression of mass outflow at the magnetic equator produces a region in which radiatively driven diffusive mass transport can levitate helium in the hottest stars, producing helium enrichment of the atmosphere. In the cooler stars, the surface becomes progressively more helium poor.
An explicit solution of the spectral radiance leaving an arbitrary point on the wall of a spherical cavity with diffuse reflectivity is obtained. The solution is applicable to spheres with an arbitrary number of openings of any size and shape, an arbitrary number of light sources with possible non-diffuse characteristics, a non-uniform sphere wall temperature distribution, non-uniform and non-diffuse sphere wall emissivity and non-uniform but diffuse sphere wall spectral reflectivity. A general measurement equation describing the output of a sensor with a given field of view, angular and spectral response measuring the sphere output is obtained. The results are applied to the Earth Radiation Budget Experiment (ERBE) integrating sphere. The sphere wall radiance uniformity, loading effects and non-uniform wall temperature effects are investigated. It is shown that using appropriate interpretation and processing, a high-accuracy short-wave calibration of the ERBE sensors can be achieved.
The Eulerian mapping closure approach is developed for uncertainty propagation in computational fluid mechanics. The approach is used to study the Probability Density Function (PDF) for the concentration of species advected by a random shear flow. An analytical argument shows that fluctuation of the concentration field at one point in space is non-Gaussian and exhibits stretched exponential form. An Eulerian mapping approach provides an appropriate approximation to both convection and diffusion terms and leads to a closed mapping equation. The results obtained describe the evolution of the initial Gaussian field, which is in agreement with direct numerical simulations.