Search NASA⌕ Search

Engineering topics

Oran, E. S.

Publications and source records attributed to Oran, E. S..

Unsteady numerical simulations of the stability and dynamics of flames

In this report we describe the research performed at the Naval Research Laboratory in support of the NASA Microgravity Science and Applications Program over the past three years (from Feb. 1992) with emphasis on the work performed since the last microgravity combustion workshop. The primary objective of our research is to develop an understanding of the differences in the structure, stability, dynamics and extinction of flames in earth gravity and in microgravity environments. Numerical simulations, in which the various physical and chemical processes can be independently controlled, can significantly advance our understanding of these differences. Therefore, our approach is to use detailed time-dependent, multi-dimensional, multispecies numerical models to perform carefully designed computational experiments. The basic issues we have addressed, a general description of the numerical approach, and a summary of the results are described in this report. More detailed discussions are available in the papers published which are referenced herein. Some of the basic issues we have addressed recently are (1) the relative importance of wall losses and gravity on the extinguishment of downward-propagating flames; (2) the role of hydrodynamic instabilities in the formation of cellular flames; (3) effects of gravity on burner-stabilized flames, and (4) effects of radiative losses and chemical-kinetics on flames near flammability limits. We have also expanded our efforts to include hydrocarbon flames in addition to hydrogen flames and to perform simulations in support of other on-going efforts in the microgravity combustion sciences program. Modeling hydrocarbon flames typically involves a larger number of species and a much larger number of reactions when compared to hydrogen. In addition, more complex radiation models may also be needed. In order to efficiently compute such complex flames recent developments in parallel computing have been utilized to develop a state-of-the-art parallel flame code. This is discussed below in some detail after a brief discussion of the numerical models.

Kailasanath, K.↗

A barely implicit correction for flux-corrected transport

In virtue of adding one elliptic equation that must be implicitly solved, the 'barely implicit correction' method presented removes the stringent sidestep limit imposed by sound speed in the explicit methods. The barely implicit method is presently combined with a flux-corrected transport algorithm, which has been rendered fourth-order by the appropriate subtraction of corrected fluxes, in order to accurately represent the sharp gradients in subsonic flows. The multigrid MGRID technique (DeVore, 1984) is used to efficiently solve the elliptic pressure equation.

Patnaik, G.↗

Implicit flux-corrected transport scheme for low-speed flow

A new high-order nonlinear Eulerian convection algorithm, the barely implicit correction (BIC)-FCT method, is presented. The BIC-FCT makes it possible to perform accurate, multidimensional, slow-flow calculations in which the fluid expansion is important. Owing to its low cost even in two dimensions, fairly detailed chemistry models and other physical processes can be included. Areas of BIC-FCT applicability include premixed flames, diffusion flames, and turbulent jet flames.

Patnaik, G.↗

Systematic development of reduced reaction mechanisms for dynamic modeling

A method for systematically developing a reduced chemical reaction mechanism for dynamic modeling of chemically reactive flows is presented. The method is based on the postulate that if a reduced reaction mechanism faithfully describes the time evolution of both thermal and chain reaction processes characteristic of a more complete mechanism, then the reduced mechanism will describe the chemical processes in a chemically reacting flow with approximately the same degree of accuracy. Here this postulate is tested by producing a series of mechanisms of reduced accuracy, which are derived from a full detailed mechanism for methane-oxygen combustion. These mechanisms were then tested in a series of reactive flow calculations in which a large-amplitude sinusoidal perturbation is applied to a system that is initially quiescent and whose temperature is high enough to start ignition processes. Comparison of the results for systems with and without convective flow show that this approach produces reduced mechanisms that are useful for calculations of explosions and detonations. Extensions and applicability to flames are discussed.

Frenklach, M.↗

Detailed studies of the dynamics and energetics of coronal bullets

Coronal bullets are small ejecta of cool, dense plasma observed to accelerate through the solar atmosphere from 20 to 450 km/s. The NRL Dynamic Flux Tube Model has been used to simulate the evolving physical properties of these dynamic events. The present calculations utilize an adaptive-gridding technique to resolve the fine structure within and around the bullets. In this work, an identification was made of a component of shocked plasma which piles up ahead of the bullet and eventually dominates both the dynamics and heating of the original bullet mass. The observational consequences of this shocked component are discussed in terms of the available HRTS EUV data, and suggestions are made for optimizing future observations of this phenomenon. An investigation has also been conducted of the structure of the bullet material visible in EUV spectral lines and the observable characteristics of the EUV-emitting plasma. Finally, the most likely mechanisms for accelerating the bullets, as well as favorable sites of origin are evaluated.

Karpen, J. T.↗

Numerical simulations of fuel droplet flows using a Lagrangian triangular mesh

The incompressible, Lagrangian, triangular grid code, SPLISH, was converted for the study of flows in and around fuel droplets. This involved developing, testing and incorporating algorithms for surface tension and viscosity. The major features of the Lagrangian method and the algorithms are described. Benchmarks of the algorithms are given. Several calculations are presented for kerosene droplets in air. Finally, extensions which make the code compressible and three dimensional are discussed.

Fritts, M. J.↗

Numerical Simulation for Droplet Combustion Using Lagrangian Hydrodynamics

A predictive model of spray combustion must incorporate models for the wide variety of physical environments in a practical combustor. In regions where droplets are closely spaced, combustion resembles a diffusion flame; where they are well separated, an envelope or wake flame results. The relative velocity field between the fuel droplets and oxidizer in influences boundary layer development about the droplet, recirculating flow patterns, and droplet shape and stability. A model must encompass these interacting temporal and spatial effects as well as complicated combustor boundaries. The objective of the current work is to develop the triangular gridding method for describing the individual and collective properties of vaporizing and burning fuel droplets.

Fritts, M. J.↗

Numerical simulations of loops heated to solar flare temperatures. I - Gasdynamics. II - X-ray and UV spectroscopy

The NRL's Dynamic Flux Tube Model is used to numerically simulate the dynamic response of a coronal magnetic loop to an energy input of the order encountered in solar flares. The coronal plasma is heated by the deposition of flare energy at the top of the loop to more than 10 million K, yielding a conduction front that moves toward the chromosphere, where the plasma is heated by the large downward conductive flux and ablates upward to the coronal part of the loop at velocities of a few hundred km/sec. The conduction front simultaneously produces chromospheric ablation and compresses the material ahead of it. With the aid of compressional instabilities, the compressed plasma grows throughout the flare heating phase, presenting a possible source of the flare optical continuum emission which is correlated with soft X-ray radiation. The observational consequences of rapidly heated loop gas dynamic processes are discussed. In the second part of this presentation, the dynamical calculation results previously obtained are used to predict the spectral line intensities, profiles and wavelengths of several X-ray lines and the UV line of Fe XXI at 1354.1 A. Three different viewing orientations of the loop are considered.

Cheng, C.-C.↗

The dynamics of accelerating coronal bullets

Results are presented of computer simulations of the jets that accelerate through the corona at velocities of 50 to 400 km/s. Particular emphasis is placed on the sensitivity of the induced acceleration to the form in which energy is put into the system. A comparison is made between the observed and predicted physical characteristics of the high-velocity bullets; the potential contribution of the bullets to the mass and energy balance of the solar corona is considered. It is found that the velocity and temperature evolution of the bullets can be modeled successfully by assuming energy input in the form of an external force, pushing continuously on the ejected material. From the physical characteristics of the model bullets and the energy input required to reproduce the observations, it is concluded that the bullets may constitute a significant fraction of the coronal mass flux but only a negligible component of the coronal energy budget.

Karpen, J. T.↗

Solar transition region response to variations in the heating rate

The response of a numerical model for the upper chromosphere, transition region, and corona to variations in the energy input has been examined. The numerical model solves the set of one-dimensional two-fluid hydrodynamic equations in a simple vertical magnetic flux tube. The atmosphere responds to both the increase and decrease in energy deposition by smoothly readjusting the temperature gradient and the amount of material in the region of peak radiating efficiency to radiate away energy being deposited. At no time during this readjustment is a departure from a thin laminar transition region structure seen. In addition, a time-dependent description of the nonequilibrium ionization of all of the ionization stages of oxygen has been included. This calculation is coupled with the self-consistent calculations of the dynamical variables. It is found that the nonequilibrium ionization balance calculations for both heating and cooling small loops in the quiet sun predict relative ionic abundances which differ substantially from those which would be predicted by an equilibrium calculation

Mariska, J. T.↗

The condensational instability in the solar transition region and corona

The stability of plasmas at temperatures and densities typical of the solar transition region and corona is investigated using both a linear analysis and nonlinear time-dependent numerical simulations. Growth rates, decay rates, and oscillation frequencies of the perturbations determined from the linear analysis are in excellent agreement with the simulations. The nonlinear regime is characterized by a bifurcation of the plasma into a cool dense condensation surrounded by a hot tenuous corona. The condensation may then be accelerated by forces in the plasma such as those arising from gravity or differential heating. Finally, the results of the detailed simulation show that the transition region is a dynamically stable structure which is the result of the nonlinear evolution of the condensational instability.

Oran, E. S.↗

F region above Kauai - Measurement, model, modification

Results of the Lagopedo II experiment conducted from Kauai, Hawaii to investigate the ionospheric modification that occurs when rocket combustion products are introduced into the O(+)-rich F region are presented. The experiment involved the detonation of a chemical explosion in the F2 peak accompanied by rocket-borne measurements of ion composition and electron content in the vicinity of the explosion. The experimental data is found to be in good agreement with the predictions of a model of the nighttime ion densities in the midlatitude laminar ionosphere, with the exception of N2(+) densities before the explosion. H2O(+) and H3O(+) currents produced by considerable H2O outgassing from the rocket are used to determine a H3O(+)/H2O(+) dissociative recombination rate averaging 1.6 to 1.08, depending on model assumptions. At the time of the explosion, an ionic void 1 km in radius is observed, the boundary of which is characterized by a steep gradient in ionic densities. Evidence of variations in the concentrations of ambient ion species, new reactant species and ionic depletion by sweeping is also obtained.

Johnson, C. Y.↗

High-altitude atomic nitrogen densities

Theoretical calculations of the seasonal and diurnal variations of atomic nitrogen are compared with measurements made by the open source neutral mass spectrometer on the AE-C satellite. With the simultaneous measurements of molecular nitrogen and atomic oxygen densities as input, model calculations of odd nitrogen densities predict the same trends in atomic nitrogen as those observed. From these comparisons it is inferred that horizontal transport significantly reduces the diurnal variation of atomic nitrogen. Estimates are given of the sensitivity of atomic nitrogen densities to variations in the photoelectron flux, the neutral temperatures, and the neutral winds.

Oran, E. S.↗

Measurement of the nitric oxide altitude distribution in the mid-latitude mesosphere

A simple rocket-borne system has been developed for deriving nitric oxide concentration in the mesosphere from measurements of resonantly scattered solar ultraviolet radiation in the NO gamma band near 2150 A. The NO signal is extracted from a background of Rayleigh-scattered light by a direct measurement of this scattered light by use of an NO absorption cell. This system was utilized to measure the NO concentration above White Sands, New Mexico, in the altitude range from 69 to 101 km. The derived NO profile shows a concentration of about 30 million per cu cm at 67 and 100 km with a deep valley (more than an order of magnitude less) at about 84 km. Comparison with model calculations indicates that the time constant for NO chemistry is smaller than or comparable with that for transport.

Baker, K. D.↗

The aeronomy of odd nitrogen in the thermosphere. II - Twilight emissions

A model developed for the aeronomy of odd nitrogen in the thermosphere is used to analyze rocket measurements of N(4S) and NO densities. Data from Atmosphere Explorer were used to develop a consistent reaction kinetics model for odd nitrogen chemistry. It is concluded that most NO(+) dissociative recombination events must produce N(2D), that N(2D) is quenched by O at a rate of 1 trillionth cu cm per sec, and that the atmospheric O2 quenching rate of N(2D) is consistent with the laboratory rate. The major quenching agent of N(2D) between 140 and 220 km is atomic oxygen, and this reaction is the major source of N(4S). Peak N(4S) densities of about (20-60) million per cu cm at 140-150 km are predicted, with the variability being indicative of the model sensitivity to a factor of 2 change in the O/O2 ratio in the thermosphere.

Strobel, D. F.↗