Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffuse layer model”

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 199 records · Page 11

IRAS detection of very cold dust in the Lynds 134 cloud complex

IRAS maps at 60 and 100 microns of the complex of dark clouds containing L1345, L183, and L1780 are presented and discussed. Extended regions are apparent where the 60 micron emission is very low compared to the 100 micron emission. The 60-100 micron flux ratio has a constant value of 0.2 in the outer diffuse parts of the clouds but suddenly drops to less than 0.03 in the opaque cloud centers. A geometrical model is used to show that the drop occurs in a narrow transition layer that has not been resolved by IRAS. The observations provide direct evidence for dust temperatures colder than 15 K in dark clouds. A number of 60 micron-deficient regions are defined using the quantity Delta I(100) = I(60)/Theta, where Theta is the ratio I(60)/I(100) in the diffuse parts of the complex. Delta I(100) is highly proportional to extinction and (C-13)O column density. The relationships indicate that the dust temperature must be constant over a large fraction of the volume sampled by Delta I(100). An upper limit of 4.8 g/sq cm is obtained for the mass absorption coefficient of dust at 100 microns.

Laureijs, R. J.↗

Heat and mass transfer in combustion - Fundamental concepts and analytical techniques

Fundamental combustion phenomena and the associated flame structures in laminar gaseous flows are discussed on physical bases within the framework of the three nondimensional parameters of interest to heat and mass transfer in chemically-reacting flows, namely the Damkoehler number, the Lewis number, and the Arrhenius number which is the ratio of the reaction activation energy to the characteristic thermal energy. The model problems selected for illustration are droplet combustion, boundary layer combustion, and the propagation, flammability, and stability of premixed flames. Fundamental concepts discussed include the flame structures for large activation energy reactions, S-curve interpretation of the ignition and extinctin states, reaction-induced local-similarity and non-similarity in boundary layer flows, the origin and removal of the cold boundary difficulty in modeling flame propagation, and effects of flame stretch and preferential diffusion on flame extinction and stability. Analytical techniques introduced include the Shvab-Zeldovich formulation, the local Shvab-Zeldovich formulation, flame-sheet approximation and the associated jump formulation, and large activation energy matched asymptotic analysis. Potentially promising research areas are suggested.

Law, C. K.↗

Numerical simulation of compressible, turbulent, two-phase flow

A computer program for numerically simulating compressible, turbulent, two-phase flows is described and applied. Special attention is given to flows in which dust is ingested into the turbulent boundary layer behind shock waves moving over the earth's surface. it is assumed that the two phases are interpenetrating continua which are coupled by drag forces and heat transfer. The particle phase is assumed to be dilute, and turbulent effects are modeled by zero- and two-equation eddy viscosity models. An important feature of the turbulence modeling is the treatment of surface boundary conditions which control the ingestion of particles into the boundary layer by turbulent friction and diffusion. The numerical method uses second-order implicit upwind differencing of the inviscid terms of the equations and second-order central differencing of the viscous terms. A diagonal form of the implicit algorithm is used to improve efficiency, and the transformation to a curvilinear coordinate system is accomplished by the finite volume techniques. Applications to a series of representative flows include a two-phase nozzle flow, the steady flow of air over a sand bed, and the air flow behind a normal shock wave in uniform motion over a sand bed. Results of the latter two applications are compared with experimental results.

Coakley, t. J.↗

SiC Recession Due to SiO2 Scale Volatility Under Combustion Conditions: Thermodynamics and Gaseous Diffusion Model - Part 2

In combustion environments, volatilization of SiO2 to Si-O-H(g) species is a critical issue. Available thermochemical data for Si-O-H(g) species were used to calculate boundary layer controlled fluxes from SiO2. Calculated fluxes were compared to volatilization rates Of SiO2 scales grown on SiC which were measured in Part 1 of this paper. Calculated volatilization rates were also compared to those measured in synthetic combustion gas furnace tests. Probable vapor species were identified in both fuel-lean and fuel-rich combustion environments based on the observed pressure, temperature and velocity dependencies as well as the magnitude of the volatility rate. Water vapor is responsible for the degradation of SiO2 in the fuel-lean environment. Silica volatility in fuel-lean combustion environments is attributed primarily to the formation of Si(OH)4(g) with a small contribution of SiO(OH)2(g).

Opila, Elizabeth J.↗

Magnetic Flux Compression Reactor Concepts for Spacecraft Propulsion and Power (MSFC Center Director's Discretionary Fund; Project No. 99-24)

This technical publication (TP) examines performance and design issues associated with magnetic flux compression reactor concepts for nuclear/chemical pulse propulsion and power. Assuming that low-yield microfusion detonations or chemical detonations using high-energy density matter can eventually be realized in practice, various magnetic flux compression concepts are conceivable. In particular, reactors in which a magnetic field would be compressed between an expanding detonation-driven plasma cloud and a stationary structure formed from a high-temperature superconductor are envisioned. Primary interest is accomplishing two important functions: (1) Collimation and reflection of a hot diamagnetic plasma for direct thrust production, and (2) electric power generation for fusion standoff drivers and/or dense plasma formation. In this TP, performance potential is examined, major technical uncertainties related to this concept accessed, and a simple performance model for a radial-mode reactor developed. Flux trapping effectiveness is analyzed using a skin layer methodology, which accounts for magnetic diffusion losses into the plasma armature and the stationary stator. The results of laboratory-scale experiments on magnetic diffusion in bulk-processed type II superconductors are also presented.

Litchford, R. J.↗

Simple Two-Dimensional Corrections for One-Dimensional Pulse Tube Models

One-dimensional oscillating flow models are very useful for designing pulse tubes. They are simple to use, not computationally intensive, and the physical relationship between temperature, pressure and mass flow are easy to understand when used in conjunction with phasor diagrams. They do not possess, however, the ability to directly calculate thermal and momentum diffusion in the direction transverse to the oscillating flow. To account for transverse effects, lumped parameter corrections, which are obtained though experiment, must be used. Or two-dimensional solutions of the differential fluid equations must be obtained. A linear two-dimensional solution to the fluid equations has been obtained. The solution provides lumped parameter corrections for one-dimensional models. The model accounts for heat transfer and shear flow between the gas and the tube. The complex Nusselt number and complex shear wall are useful in describing these corrections, with phase relations and amplitudes scaled with the Prandtl and Valensi numbers. The calculated ratio, a, between a two-dimensional solution of the oscillating temperature and velocity and a one-dimensional solution for the same shows a scales linearly with Va for Va less than 30. In this region alpha less than 0.5, that is, the enthalpy flow calculated with a two-dimensional model is 50% of a calculation using a one-dimensional model. For Va greater than 250, alpha = 0.8, showing that diffusion is still important even when it is confined to a thing layer near the tube wall.

Lee, J. M.↗

Theory of plasma contactors in ground-based experiments and low Earth orbit

Previous theoretical work on plasma contactors as current collectors has fallen into two categories: collisionless double layer theory (describing space charge limited contactor clouds) and collisional quasineutral theory. Ground based experiments at low current are well explained by double layer theory, but this theory does not scale well to power generation by electrodynamic tethers in space, since very high anode potentials are needed to draw a substantial ambient electron current across the magnetic field in the absence of collisions (or effective collisions due to turbulence). Isotropic quasineutral models of contactor clouds, extending over a region where the effective collision frequency upsilon sub e exceeds the electron cyclotron frequency omega sub ce, have low anode potentials, but would collect very little ambient electron current, much less than the emitted ion current. A new model is presented, for an anisotropic contactor cloud oriented along the magnetic field, with upsilon sub e less than omega sub ce. The electron motion along the magnetic field is nearly collisionless, forming double layers in that direction, while across the magnetic field the electrons diffuse collisionally and the potential profile is determined by quasineutrality. Using a simplified expression for upsilon sub e due to ion acoustic turbulence, an analytic solution has been found for this model, which should be applicable to current collection in space. The anode potential is low and the collected ambient electron current can be several times the emitted ion current.

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

A simple biosphere model (SiB) for use within general circulation models

A simple realistic biosphere model for calculating the transfer of energy, mass and momentum between the atmosphere and the vegetated surface of the earth has been developed for use in atmospheric general circulation models. The vegetation in each terrestrial model grid is represented by an upper level, representing the perennial canopy of trees and shrubs, and a lower level, representing the annual cover of grasses and other heraceous species. The vegetation morphology and the physical and physiological properties of the vegetation layers determine such properties as: the reflection, transmission, absorption and emission of direct and diffuse radiation; the infiltration, drainage, and storage of the residual rainfall in the soil; and the control over the stomatal functioning. The model, with prescribed vegetation parameters and soil interactive soil moisture, can be used for prediction of the atmospheric circulation and precipitaion fields for short periods of up to a few weeks.

Sellers, P. J.↗

Thermal Pollution Mathematical Model: Verification of One-Dimensional Numerical Model at Lake Keowee - Volume 2

A one dimensional model for studying the thermal dynamics of cooling lakes was developed and verified. The model is essentially a set of partial differential equations which are solved by finite difference methods. The model includes the effects of variation of area with depth, surface heating due to solar radiation absorbed at the upper layer, and internal heating due to the transmission of solar radiation to the sub-surface layers. The exchange of mechanical energy between the lake and the atmosphere is included through the coupling of thermal diffusivity and wind speed. The effects of discharge and intake by power plants are also included. The numerical model was calibrated by applying it to Cayuga Lake. The model was then verified through a long term simulation using Lake Keowee data base. The comparison between measured and predicted vertical temperature profiles for the nine years is good. The physical limnology of Lake Keowee is presented through a set of graphical representations of the measured data base.

Lee, S. S.↗

Parametric studies with an atmospheric diffusion model that assesses toxic fuel hazards due to the ground clouds generated by rocket launches

Parametric studies were made with a multilayer atmospheric diffusion model to place quantitative limits on the uncertainty of predicting ground-level toxic rocket-fuel concentrations. Exhaust distributions in the ground cloud, cloud stabilized geometry, atmospheric coefficients, the effects of exhaust plume afterburning of carbon monoxide CO, assumed surface mixing-layer division in the model, and model sensitivity to different meteorological regimes were studied. Large-scale differences in ground-level predictions are quantitatively described. Cloud alongwind growth for several meteorological conditions is shown to be in error because of incorrect application of previous diffusion theory. In addition, rocket-plume calculations indicate that almost all of the rocket-motor carbon monoxide is afterburned to carbon dioxide CO2, thus reducing toxic hazards due to CO. The afterburning is also shown to have a significant effect on cloud stabilization height and on ground-level concentrations of exhaust products.

Stewart, R. B.↗

Similarity theory of the buoyantly interactive planetary boundary layer with entrainment

A similarity model is developed for the vertical profiles of turbulent flow variables in an entraining turbulent boundary layer of arbitrary buoyant stability. In the general formulation the vertical profiles, internal rotation of the velocity vector, discontinuities or jumps at a capping inversion and bulk aerodynamic coefficients of the boundary layer are given by solutions to a system of ordinary differential equations in the similarity variable. To close the system, a formulation for buoyantly interactive eddy diffusivity in the boundary layer is introduced which recovers Monin-Obukhov similarity near the surface and incorporates a hypothesis accounting for the observed variation of mixing length throughout the boundary layer. The model is tested in simplified versions which depend only on roughness, surface buoyancy, and Coriolis effects by comparison with planetary-boundary-layer wind- and temperature-profile observations, measurements of flat-plate boundary layers in a thermally stratified wind tunnel and observations of profiles of terms in the turbulent kinetic-energy budget of convective planetary boundary layers. On balance, the simplified model reproduced the trend of these various observations and experiments reasonably well, suggesting that the full similarity formulation be pursued further.

Hoffert, M. I.↗

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗

Numerical simulation of a supersonic reacting mixing layer

In order to arrive at physical models that can adequately describe supersonic combustion, and develop accurate and efficient numerical techniques for the solution of such models' governing equations, a computer program has been developed for the study of reacting flows which considers the multicomponent diffusion and convection of important chemical species, as well as their finite state reaction and the interaction of the fluid mechanics and the chemistry that occurred. The code employs a hybrid Chebyshev pseudospectral technique for integration of the models' resulting governing equations; the program is here used to study a spatially developing and reacting mixing layer.

Drummond, J. Philip↗

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. The Roe approximate Riemann solution scheme or the computationally less expensive Advection Upstream Splitting Method (AUSM) flux-splitting scheme are used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passage and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. Roe's approximate Riemann solution scheme or the computationally less expensive advection upstream splitting method (AUSM) flux-splitting scheme is used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passages and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗

Analysis of epitaxial drift field N on P silicon solar cells

The performance of epitaxial drift field silicon solar cell structures having a variety of impurity profiles was calculated. These structures consist of a uniformly doped P-type substrate layer, and a P-type epitaxial drift field layer with a variety of field strengths. Several N-layer structures were modeled. A four layer solar cell model was used to calculate efficiency, open circuit voltage and short circuit current. The effect on performance of layer thickness, doping level, and diffusion length was determined. The results show that peak initial efficiency of 18.1% occurs for a drift field thickness of about 30 micron with the doping rising from 10 to the 17th power atoms/cu cm at the edge of the depletion region to 10 to the 18th power atoms/cu cm in the substrate. Stronger drift fields (narrow field regions) allowed very high performance (17% efficiency) even after irradiation to 3x10 to the 14th power 1 MeV electrons/sq cm.

Baraona, C. R.↗

Study of diamond film growth and properties

The objective was to study diamond film growth and its properties in order to enhance the laser damage threshold of substrate materials. Calculations were performed to evaluate laser induced thermal stress parameter, R(sub T) of diamond. It is found that diamond has several orders of magnitude higher in value for R(sub T) compared to other materials. Thus, the laser induced damage threshold (LIDT) of diamond is much higher. Diamond films were grown using a microwave plasma enhanced chemical vapor deposition (MPECVD) system at various conditions of gas composition, pressure, temperature, and substrate materials. A 0.5 percent CH4 in H2 at 20 torr were ideal conditions for growing of high quality diamond films on substrates maintained at 900 C. The diamond films were polycrystalline which were characterized by scanning electron microscopy (SEM) and Raman scattering spectroscopy. The top surface of the growing film is always rough due to the facets of polycrystalline film while the back surface of the film replicates the substrate surface. An analytical model based on two dimensional periodic heat flow was developed to calculate the effective in-plane (face parallel) diffusivity of a two layer system. The effective diffusivity of diamond/silicon samples was measured using a laser pulse technique. The thermal conductivity of the films was measured to be 13.5 W/cm K, which is better than that of a type Ia natural diamond. Laser induced damage experiments were performed on bare Si substrates, diamond film coated Si, and diamond film windows. Significant improvements in the LIDT were obtained for diamond film coated Si compared to the bare Si.

Albin, Sacharial↗