Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional approximations”

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 235 records · Page 13

Low-energy ion distribution functions on a magnetically quiet day at geostationary altitude /L = 7/

Ion energy and pitch angle distribution functions are examined for a magnetically quiet day using averaged data from ATS 6. For both field-aligned and perpendicular fluxes, the populations have a mixture of characteristic energies, and the distribution functions can be fairly well approximated by Maxwellian distributions over three different energy bands in the range 3-600 eV. Pitch angle distributions varying with local time, and energy distributions are used to compute total ion density. Pitch angle scattering mechanisms responsible for the observed transformation of pitch angle distribution are examined, and it is found that a magnetic noise of a certain power spectral density belonging to the electromagnetic ion cyclotron mode near the ion cyclotron frequency can be effective in trapping the field aligned fluxes by pitch angle scattering.

Singh, N.↗

Radio and UV Observations of High Latitude HI

Interstellar Lyman alpha observations toward high latitude stars, which give the column density of H I between Earth and the star, and 21 cm radio observations, which give the total H I column density to and beyond the star, can be combined to map out the amount of neutral gas as a function of distance from the Sun. Observations of the 21 cm and Lyman alpha lines made toward approximately 50 high latitude OB stars are presented. These data are used to discuss the distribution of H I in the solar neighborhood.

Lockman, F. J.↗

Solutions to bi-Maxwellian transport equations for radial solar wind beyond 28 R(S)

This paper presents solar wind solutions for radial flow between 28 solar radii and 1 AU using the bi-Maxwellian-based 16-moment set of transport equations. In addition to the number density, drift velocity, and parallel and perpendicular temperatures, the 16-moment equations account for the transport of both longitudinal and transverse thermal energies as well as stress. Also, using the 16-moment approximation for the distribution function and assuming plasma parameter values characteristic of the solar wind, contour plots are generated for the proton velocity distribution function. It is shown how the shape of these plots depends on various macroscopic plasma parameters.

Demars, H. G.↗

Bidirectional Reflectance of a Macroscopically Flat, High-Albedo Particulate Surface: An Efficient Radiative Transfer Solution and Applications to Regoliths

Many remote sensing applications rely on accurate knowledge of the bidirectional reflection function (BRF) of surfaces composed of discrete, randomly positioned scattering particles. Theoretical computations of BRFs for plane-parallel particulate layers are usually reduced to solving the radiative transfer equation (RTE) using one of existing exact or approximate techniques. Since semi-empirical approximate approaches are notorious for their low accuracy, violation of the energy conservation law, and ability to produce unphysical results, the use of numerically exact solutions of RTE has gained justified popularity. For example, the computation of BRFs for macroscopically flat particulate surfaces in many geophysical publications is based on the adding-doubling (AD) and discrete ordinate (DO) methods. A further saving of computer resources can be achieved by using a more efficient technique to solve the plane-parallel RTE than the AD and DO methods. Since many natural particulate surfaces can be well represented by the model of an optically semi-infinite, homogeneous scattering layer, one can find the BRF directly by solving the Ambartsumian's nonlinear integral equation using a simple iterative technique. In this way, the computation of the internal radiation field is avoided and the computer code becomes highly efficient and very accurate and compact. Furthermore, the BRF thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. In this paper, we discuss numerical aspects and the computer implementation of this technique, examine the applicability of the Henyey-Greenstein phase function and the sigma-Eddington approximation in BRF and flux calculations, and describe sample applications demonstrating the potential effect of particle shape on the bidirectional reflectance of flat regolith surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web, and can be applied to a wide range of remote sensing problems. BRF computations for undulated (macroscopically rough) surfaces are more complicated and often rely on time consuming Monte Carlo procedures. This approach is especially inefficient for optically thick, weakly absorbing media (e.g., snow and desert surfaces at visible wavelengths since a photon may undergo many internal scattering events before it exists the medium or is absorbed. However, undulated surfaces can often be represented as collections of locally flat tilted facets characterized by the BRF found from the traditional plane parallel RTE. In this way the MOnte Carlo procedure could be used only to evaluate the effects of surface shadowing and multiple surface reflections, thereby bypassing the time-consuming ray tracing inside the medium and providing a great savings of CPU time.

Mishchenko, Michael I.↗

Computation of Domain-Averaged Shortwave Irradiance by a One-Dimensional Algorithm Incorporating Correlations between Optical Thickness and Direct Incident Radiation

A one-dimensional radiative transfer algorithm that accounts for correlations between the optical thickness and the incident direct solar radiation is developed to compute the domain-averaged shortwave irradiance profile. It divides the direct irradiance into four components and treats the direct irradiance in two separate, clear and cloudy columns to account for the fact that clouds attenuate the direct irradiance more than clear-sky. The horizontal inhomogeneity of clouds in the cloudy column is treated by the gamma weighted two-stream approximation, which assumes that the optical thickness of clouds follows a gamma distribution. The algorithm inputs the cloud fraction, cumulative cloud fraction as a function of height, and a parameter expressing the shape of the probability density function of the cloud optical thickness distribution in addition to inputs required for a two-stream radiative transfer model. These cloud property inputs can be obtained using ground- and satellite-based instruments. Therefore, the algorithm can treat realistic cloud overlap features and horizontal inhomogeneity of clouds in a framework of one- dimensional radiative transfer. Heating rates computed by the algorithm using cloud fields generated by cloud resolving models agree with those computed with a Monte Carlo model. If optical properties in computational layers that divide a vertically extensive cloud are correlated, the irradiance profile computed by the algorithm further improves.

Kato, S.↗

Light Echos in Kerr Geometry: A Source of High Frequency QPOs from Random X-ray Bursts

We propose that high frequency quasi-periodic oscillations (HFQPOs) can be produced from randomly-formed X-ray bursts (flashes) by plasma interior to the ergosphere of a rapidly-rotating black hole. We show by direct computation of their orbits that the photons comprising the observed X-ray light curves, if due to a multitude of such flashes, are affected significantly by the black hole's dragging of inertial frames; the photons of each such burst arrive to an observer at infinity in multiple (double or triple), distinct 'bunches' separated by a roughly constant time lag of t/M approximately equal to 14, regardless of the bursts' azimuthal position. We argue that every other such 'bunch' represents photons that follow trajectories with an additional orbit around the black hole at the photon circular orbit radius (a photon 'echo'). The presence of this constant lag in the response function of the system leads to a QPO feature in its power density spectra, even though the corresponding light curve consists of a totally stochastic signal. This effect is by and large due to the black hole spin and is shown to gradually diminish as the spin parameter a decreases or the radial position of the burst moves outside the static limit surface (ergosphere). Our calculations indicate that for a black hole with Kerr parameter of a/M=0.99 and mass of M=10*Msun the QPO is expected at a frequency of approximately 1.3-1.4 kHz. We discuss the plausibility and observational implications of our model/results as well as its limitations.

Fukumura, K.↗

Far-infrared rotational emission by carbon monoxide

Accurate theoretical collisional excitation rates are used to determine the emissivities of CO rotational lines 10 to the 4th power/cu cm n(H2), 100 K T 2000 K, and J 50. An approximate analytic expression for the emissitivities which is valid over most of this region is obtained. Population inversions in the lower rotational levels occur for densities n(H2) approximately 10 (to the 3rd to 5th power)/cu cm and temperatures T approximately 50 K. Interstellar shocks observed edge on are a potential source of millimeter wave CO maser emission. The CO rotational cooling function suggested by Hollenbach and McKee (1979) is verified, and accurate numerical values given. Application of these results to other linear molecules should be straightforward.

Mckee, C. F.↗

Theory of inhomogeneous fluids

A theory of nonuniform liquids is presented which is based on the Yvon-Born-Green equation for the one-particle density and the Ornstein-Zernike equation of an inhomogeneous system. The necessary closure is effected by exploiting the solution of a modified hypernetted-chain equation and making a local approximation on the (highly universal) bridge function. The theory is successfully applied to model problems that have also been studied by direct-simulation methods. A brief generalization to quantum liquids is also given.

Nieminen, R. M.↗

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

EUPDF: Eulerian Monte Carlo Probability Density Function Solver for Applications With Parallel Computing, Unstructured Grids, and Sprays

The success of any solution methodology used in the study of gas-turbine combustor flows depends a great deal on how well it can model the various complex and rate controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as convective and radiative heat transfer and other phenomena. The phenomena to be modeled, which are controlled by these processes, often strongly interact with each other at different times and locations. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. The influence of turbulence in a diffusion flame manifests itself in several forms, ranging from the so-called wrinkled, or stretched, flamelets regime to the distributed combustion regime, depending upon how turbulence interacts with various flame scales. Conventional turbulence models have difficulty treating highly nonlinear reaction rates. A solution procedure based on the composition joint probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices (such as extinction, blowoff limits, and emissions predictions) because it can account for nonlinear chemical reaction rates without making approximations. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on the PDF method to unstructured grids, parallel computing, and sprays. EUPDF, which was developed by M.S. Raju of Nyma, Inc., was designed to be massively parallel and could easily be coupled with any existing gas-phase and/or spray solvers. EUPDF can use an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements. The application of the PDF method showed favorable results when applied to several supersonic-diffusion flames and spray flames. The EUPDF source code will be available with the National Combustion Code (NCC) as a complete package.

Raju, M. S.↗

Blade Motion and Bouncing Tests of KD-1 Autogiro

There are forwarded herewith the results of blade motion and bouncing tests on the Kellett KD-1 three-bladed autogiro. Motion picture records and two-component accelerometer records were taken in flight during glides at air speeds from 30 miles per hour to 100 miles per hour indicator readings. Calibration curves of correct indicated air speed and rotor speed as functions of air speed meter reading were established with a trailing pitot-static head and a rotoscope, at 2,000 ft. altitude and an air density of 0.00231 slug/ cu. ft., all tests being made at approximately that density.

Wheatley, John B.↗

QSO evolution and the intergalactic medium

A model is developed of a cold, photoionized intergalactic medium in which only discrete, evolving sources (quasars) produce the observed extragalactic hard X-ray background. Constraints on the model are imposed by observed optical depth limits to various redshifts, far UV and visible flux measurements, X-ray background and source counts in the region 2-6 keV, and it makes use of results of the Palomar Bright Quasar Survey and the Braccesi Quasar survey to obtain expressions for the quasar evolution function, absolute luminosity, present number density and switch-on redshift. Results of computer modeling reveal an intergalactic medium with a density with respect to closure less than 0.35, temperature of 16,000 K and fractional ion abundances of approximately 0.000001 for H I and 0.001 for He II, under the assumption of a Hubble constant of 50 km/sec per Mpc. In addition, spectral indices of -0.4 and -1.3 are obtained for the regions of the emission spectrum above and below 4 x 10 to the 18th Hz.

Sherman, R. D.↗

The origin and evolution of atmospheric oxygen

This paper discusses the chemical processes involved in the evolution of the earth's atmospheric oxygen and ozone, as well as the sources, sinks, and transfer rates of oxygen in the present atmosphere. Special attention is given the evolution of atmospheric O3 as a function of the buildup of O2, with the results of calculations presented as the vertical profiles of O3, in terms of the present atmospheric level (PAL) oxygen values. Calculations show that the total O3 column density that is approximately half of the present level was reached when atmospheric oxygen level reached 0.1 PAL. At this level of ozone, the biological shielding of the earth's surface from the UV radiation is believed to have been achieved.

Levine, Joel S.↗

The mass function of Seyfert 1 nuclei

The first mass function of Seyfert 1 nuclei is derived from optical spectra of the complete CfA sample of Seyfert galaxies by estimating the mass for each object from a dynamical relation. An independent estimate is also derived using a complete infrared-selected sample. The two mass functions are indistinguishable. The mean mass of Seyfert 1 nuclei is about 2 x 10 to the 7th solar masses, and the integrated mass density is about 6 x 10 to the 11th solar masses/cu Gpc. This is approximately two orders of magnitude less than the value inferred from the energetics associated with quasar counts. A careful analysis of the various parameters and assumptions involved suggests that this large difference is not due to systematic errors in the determinations. Therefore, the bulk of mass related to the accretion processes connected with past quasar activity does not reside in Seyfert 1 nuclei. Instead, the remnants of past activity must be present in a much larger number of galaxies, and a one-to-one relation between distant and local active galactic nuclei seems then to be excluded.

Padovani, P.↗

Hot gas in the cold dark matter scenario: X-ray clusters from a high-resolution numerical simulation

A new, three-dimensional, shock-capturing hydrodynamic code is utilized to determine the distribution of hot gas in a standard cold dark matter (CDM) model of the universe. Periodic boundary conditions are assumed: a box with size 85 h(exp -1) Mpc having cell size 0.31 h(exp -1) Mpc is followed in a simulation with 270(exp 3) = 10(exp 7.3) cells. Adopting standard parameters determined from COBE and light-element nucleosynthesis, sigma(sub 8) = 1.05, omega(sub b) = 0.06, and assuming h = 0.5, we find the X-ray-emitting clusters and compute the luminosity function at several wavelengths, the temperature distribution, and estimated sizes, as well as the evolution of these quantities with redshift. We find that most of the total X-ray emissivity in our box originates in a relatively small number of identifiable clusters which occupy approximately 10(exp -3) of the box volume. This standard CDM model, normalized to COBE, produces approximately 5 times too much emission from clusters having L(sub x) is greater than 10(exp 43) ergs/s, a not-unexpected result. If all other parameters were unchanged, we would expect adequate agreement for sigma(sub 8) = 0.6. This provides a new and independent argument for lower small-scale power than standard CDM at the 8 h(exp -1) Mpc scale. The background radiation field at 1 keV due to clusters in this model is approximately one-third of the observed background, which, after correction for numerical effects, again indicates approximately 5 times too much emission and the appropriateness of sigma(sub 8) = 0.6. If we have used the observed ratio of gas to total mass in clusters, rather than basing the mean density on light-element nucleosynthesis, then the computed luminosity of each cluster would have increased still further, by a factor of approximately 10. The number density of clusters increases to z approximately 1, but the luminosity per typical cluster decreases, with the result that evolution in the number density of bright clusters is moderate in this redshift range, showing a broad peak near z = 0.7, and then a rapid decline above redshift z = 3. Detailed computations of the luminosity functions in the range L(sub x) = 10(exp 40) - 10(exp 44) ergs/s in various energy bands are presented for both cluster central regions and total luminosities to be used in comparison with ROSAT and other observational data sets. The quantitative results found disagree significantly with those found by other investigators using semianalytic techniques. We find little dependence of core radius on cluster luminosity and a dependence of temperature on luminosity given by log kT(sub x) = A + B log L(sub x), which is slightly steeper (B = 0.38) than is indicated by observations. Computed temperatures are somewhat higher than observed, as expected, in that COBE-normalized CDM has too much power on the relevant scales. A modest average temperature gradient is found, with temperatures dropping to 90% of central values at 0.4 h(exp -1) Mpc and 70% of central values at 0.9 h(exp -1) Mpc. Examining the ratio of gas to total mass in the clusters normalized to Omega(sub B) h(exp 2) = 0.015, and comparing with observations, we conclude, in agreement with White (1991), that the cluster observations argue for an open universe.

Kang, Hyesung↗

Ground state energy of electrons in a static point-ion lattice

The ground state energy of a neutral collection of protons and electrons was investigated under the assumption that in the ground state configuration, static protons occupy the sites of a rigid Bravais lattice. The Wigner-Seitz method was used in conjunction with three postulated potentials: bare Coulomb, Thomas-Fermi screening, and screening by a uniform bare background charge. Within these approximations, the exact band-minimum energy and wave functions are derived. For each of the three potentials, the approximate minimum ground state energy per proton (relative to isolated electrons and protons) is, respectively, -1.078 Ry, -1.038 Ry, and -1.052 Ry. These three minima all fall at a density of about 0.60 gm/cu cm, which is thus an approximate lower bound on the density of metallic hydrogen at its transition pressure.

Styer, D. F.↗

Statistics of Sxy estimates

The statistics of Sxy estimates derived from orthogonal-component measurements are examined. Based on results of Goodman (1957), the probability density function (pdf) for Sxy(f) estimates is derived, and a closed-form solution for arbitrary moments of the distribution is obtained. Characteristic functions are used to derive the exact pdf of Sxy(tot). In practice, a simple Gaussian approximation is found to be highly accurate even for relatively few degrees of freedom. Implications for experiment design are discussed, and a maximum-likelihood estimator for a posterior estimation is outlined.

Freilich, M. H.↗

Influence of aerodynamic roughness length on aeolian processes: Earth, Mars, Venus

The aerodynamic roughness length (z sub 0) is the height at which a wind profile assumes a zero velocity. The lower part of the atmospheric boundary layer will be impeded by friction with the surface. An increase in surface roughness will also increase the shear stress required to initiate particle entrainment by the wind. Bagnold (1941) estimated z sub 0 as being 1/30 of the mean particle size. In Nature, surface roughness is composed of nonerodible elements as well as sand-size erodible particles. To assess z sub 0 values as a function of terrain, field experiments were conducted to obtain wind profiles monitored over natural surfaces at 15 sites in the Mojave Desert, Death Valley, and Nye County, Nevada. These sites span a variety of arid-land terrains, including smooth playas, alluvial fans, and lava flows; z sub 0 values ranged from 0.0001 cm to 1 cm. These values were incorporated in a threshold model and a flux model to assess transport efficiency over such terrains in three planetary environments (Venus, Earth, and Mars), and for particle sizes ranging from 60-500 micron. Threshold and flux are a function of planetary environment, particle density and size (Dp), and z sub 0, and the shear velocity of 1.2 x U sub *t (for Dp = 250 micron and z sub 0 = 0.84). Results show that flux on Mars is approximately 14 g/(cm x s), on Earth it is approximately 3 g/(cm x s), and on Venus 0.5 g/(cm x s). Under all planetary environments, the results also show a dramatic decrease in the flux for particles greater than 200 microns when z sub 0 increases above 0.0085 cm (corresponding to sites consisting of alluvium). When z sub 0 approaches 0.03 cm (corresponding to a mantled pahoehoe lava), the flux diminishes.

Blumberg, Dan G.↗