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 253 records · Page 14

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

Nonlinear structure formation in flat cosmological models

This paper describes the formation of nonlinear structure in flat (zero curvature) Friedmann cosmological models. We consider models with two components: the usual nonrelativistic component that evolves under gravity and eventually forms the large-scale structure of the universe, and a uniform dark matter component that does not clump under gravity, and whose energy density varies with the scale factor a(t) like a(t)(sup -n), where n is a free parameter. Each model is characterized by two parameters: the exponent n and the present density parameter Omega(sub 0) of the nonrelativistic component. The linear perturbation equations are derived and solved for these models, for the three different cases n = 3, n is greater than 3, and n is less than 3. The case n = 3 is relevant to model with massive neutrinos. The presence of the uniform component strongly reduces the growth of the perturbation compared with the Einstein-de Sitter model. We show that the Meszaros effect (suppression of growth at high redshift) holds not only for n = 4, radiation-dominated models, but for all models with n is greater than 3. This essentially rules out any such model. For the case n is less than 3, we numerically integrate the perturbation equations from the big bang to the present, for 620 different models with various values of Omega(sub 0) and n. Using these solutions, we show that the function f(Omega(sub 0), n) = (a/delta(sub +))d(delta)(sub +)/da, which enters in the relationship between the present density contrast delta(sub 0) and peculiar velocity field u(sub 0) is essentially independent of n. We derive approximate solutions for the second-order perturbation equations. These second-order solutions are tested against the exact solutions and the Zel'dovich approximation for spherically symmetric perturbations in the marginally nonlinear regime (the absolute value of delta is less than or approximately 1). The second-order and Zel'dovich solutions have comparable accuracy, significantly higher than the accuracy of the linear solutions. We then investigate the dependence of the delta(sub 0) - u(sub 0) relationship upon the value of n in the nonlinear regime using the second-order solutions for marginally nonlinear, general perturbations, and the exact solutions for strongly nonlinear, spherically symmetric perturbations. In both cases, we find that the delta(sub 0) - u(sub 0) relationship remains independent of n. We speculate that this result extends to strongly nonlinear, general perturbations as well. This eliminates any hope to determine the presence of the uniform component or the value of n using dynamical methods. Finally, we compute the nonlinear evolution of the skewness of the distribution of values of delta, assuming Gaussian initial conditions. We find that the skewness is not only independent of n, but also of Omega(sub 0). Thus the skewness cannot be used to discriminate among various models with Gaussian initial conditions. However, it can be used for testing the Gaussianity of the initial conditions themselves. We conclude that the uniform component leaves no observable signature in the present large-scale structure of the universe. To determine its presence and nature, we must investigate the relationship between the past and present universe, using redshift-dependent tests.

Martel, Hugo↗

Pitch Angle Distributions of Solar Energetic Particles and the Local Scattering Properties of the Interplanetary Magnetic Field

An approximate solution of the Fokker-Planck equation containing pitch angle scattering and adiabatic focusing is discussed. For modest focusing effects the omnidirectional density obeys an ordinary diffusion equation with a modified diffusion coefficient. The anisotropic part of the distribution function is properly normalized and split into an even and an odd part. The even part is determined by the ratio between the scattering mean free path and the focusing length and by the degree of polarization of the magnetic field fluctuations. The odd part is determined by the deviation of the pitch angle scattering from isotropic scattering. The method supplies a powerful tool to obtain the local characteristics of pitch angle scattering. It is insensitive to long lasting solar injections and to moderate radial variations of the mean free path. The method is applied to solar particle events observed on Helio-1 and -2.

Beeck, J.↗

Convection and Easterly Wave Structure Observed in the Eastern Pacific Warm-Pool during EPIC-2001

During September-October 2001, the East Pacific Investigation of Climate Processes in the Coupled Ocean-Atmosphere System (EPIC-2001) ITCZ field campaign focused on studies of deep convection in the warm-pool region of the East Pacific. In addition to the TAO mooring array, observational platforms deployed during the field phase included the NOAA ship RN Ronald H. Brown, the NSF ship RN Horizon, and the NOAA P-3 and NCAR C-130 aircraft. This study combines C-band Doppler radar, rawinsonde, and surface heat flux data collected aboard the RN Brown to describe ITCZ convective structure and rainfall statistics in the eastern Pacific as a function of 3-5 day easterly wave phase. Three distinct easterly wave passages occurred during EPIC-2001. Wind and thermodynamic data reveal that the wave trough axes exhibited positively correlated U and V winds and a slight westward phase tilt with height. A relatively strong (weak) northeasterly deep tropospheric shear followed the trough (ridge) axis. Temperature and humidity perturbations exhibited mid-to upper level cooling (warming) and drying (moistening) in the northerly (trough and southerly) phase. At low levels warming (cooling) occurred in the northerly (southerly) phase with little change in the relative humidity, though mixed layer mixing ratios were larger during the northerly phase. When composited, radar, sounding, lightning and surface heat flux observations suggest the following systematic behavior as a function of wave phase: approximately zero to one quarter wavelength ahead of (behind) the wave trough in northerly (southerly) flow, larger (smaller) CAPE, lower (higher) CIN, weaker (stronger) tropospheric shear, higher (lower) conditional mean rain rates, higher (lower) lightning flash densities, and more (less) robust convective vertical structure occurred. Latent and sensible heat fluxes reached a minimum in the northerly phase and then increased through the trough, reaching a peak during the ridge phase (leading the peak in CAPE). From a radar echo coverage perspective, larger areas of light rain and slightly larger (10%) area averaged rain rates occurred in the vicinity of, and just behind, the trough axes in southerly flow. Importantly, the transition in convective structure observed across the trough axis when considered with the relatively small change in area mean rain rates suggests the presence of a transition in the vertical structure of diabatic heating across the easterly waves examined. The inferred transition in heating structure is supported by radar diagnosed divergence profiles that exhibit convective (stratiform) characteristics ahead of (behind) the trough.

Peterson, Walter A.↗

Complex Permittivity for NUW-LHT-5M at 2.466 GHz Over a Range of Densities

Microwave radiation near 2.45 GHz are being considered as a mechanism to heat lunar regolith. Because microwave absorption is in part a function of particle packing, the effect of bulk density of the lunar highland simulant, NUW-LHT-5M, was investigated. It was computed that over the practically achievable range of compaction, 1.4 g/cm 3 to 2.094 g/cm 3 , tanδ changed approximately 9%/0.2 g/cm 3 and half power depth changed is approximately 14%/0.2 g/cm 3 . Energy consumption is tied to these values.

NUW-LHT-5M↗

Properties of a B0 I stellar wind and interstellar grains derived from Ginga observations of the binary X-ray pulsar 4U 1538-52

From measurements of the X-ray eclipse phenomena of the binary X-ray pulsar 4U 1538-52, we derive properties of the stellar wind of its B0 I companion, QV Nor, and a constraint on models of interstellar grains. Estimates of the wind density as a function of the distance from the center of QV Nor are obtained from an analysis of the variation of X-ray attenuation during an eclipse egress. A Monte Carlo computation of the absorption and scattering of X-rays in the X-ray-ionized wind accounts for approximately two-thirds of the spectrum of X-rays with energies above 4.5 keV observed during the eclipse. Our upper limit on R(sub XV) is 0.06/mag, which implies that the X-ray scattering efficiency of interstellar dust is less than expected for solid grains with a size distribution of the form n(sub g)(a) approximately a(exp -3.5) in the range from 0.005 to 0.25 microns and composed of silicate (R(sub XV) = 0.22/mag) or a silicate-graphite mixture (R(sub XV) = 0.11/mag) as derived from the calculations of Martin & Rouleau (1991). This lends support to the idea (Mathis & Whiffen 1989) that interstellar grains are 'fluffy' aggregates with an average bulk density less than that of their constitutent particles. Such aggregates would have a smaller ratio of X-ray scattering efficiency to optical extinction efficiency compared with solid grains of the same material.

Clark, George W.↗

Pair production rates in mildly relativistic, magnetized plasmas

Electron-positron pairs may be produced by either one or two photons in the presence of a strong magnetic field. In magnetized plasmas with temperatures kT approximately sq mc, both of these processes may be important and could be competitive. The rates of one-photon and two-photon pair production by photons with Maxwellian, thermal bremsstrahlung, thermal synchrotron and power law spectra are calculated as a function of temperature or power law index and field strength. This allows a comparison of the two rates and a determination of the conditions under which each process may be a significant source of pairs in astrophysical plasmas. It is found that for photon densities n(gamma) or = 10 to the 25th power/cu cm and magnetic field strengths B or = 10 to the 12th power G, one-photon pair production dominates at kT approximately sq mc for a Maxwellian, at kT approximately 2 sq mc for a thermal bremsstrahlung spectrum, at all temperatures for a thermal synchrotron spectrum, and for power law spectra with indices s approximately 4.

Burns, M. L.↗

Pair production rates in mildly relativistic magnetized plasmas

Electron-positron pairs may be produced by either one or two photons in the presence of a strong magnetic field. In magnetized plasmas with temperatures kT approximately sq mc, both of these processes may be important and could be competitive. The rates of one-photon and two-photon pair production by photons with Maxwellian, thermal bremsstrahlung, thermal synchrotron and power law spectra are calculated as a function of temperature or power law index and field strength. This allows a comparison of the two rates and a determination of the conditions under which each process may be a significant source of pairs in astrophysical plasmas. It is found that for photon densities n(gamma) or = 10 to the 25th power/cu cm and magnetic field strengths B or = 10 to the 12th power G, one-photon pair production dominates at kT approximately sq mc for a Maxwellian, at kT approximately 2 sq mc for a thermal bremsstrahlung spectrum, at all temperatures for a thermal synchrotron spectrum, and for power law spectra with indices s approximately 4.

Burns, M. L.↗

Some influences of approximate values for velocity, density and total temperature sensitivities on hot wire anemometer results

There is a renewed interest in hot wire anemometry at transonic speeds. Recent results were published which indicate that at transonic speeds a heated wire is sensitive only to mass flow and total temperature, results similar to those obtained for supersonic flows. Other results were obtained to show that the sensitivity is a function of velocity, density, and total temperature, results in agreement with many of those obtained in the 1950s. An analysis of anemometry results was made to evaluate possible errors when various assumptions were made concerning the sensitivity of a heated wire to fluid flow variables.

Stainback, P. C.↗

Application of a Reynolds stress turbulence model to a supersonic radiating hydrogen-air diffusion flame

A second-order differential Reynolds Stress turbulence model has been applied to the Favre-averaged Navier-Stokes equations for the study of supersonic flows with finite-rate chemistry and radiation. An assumed Beta Probability Density Function is applied to account for the chemical source terms and the radiative flux terms in the conservation equations. A seven-species, seven-reaction finite rate chemistry mechanism is used to simulate the combustion process. The tangent slab approximation is used in radiative flux formulation. A pseudo-gray gas model is used to represent the absorption-emission characteristics of the participating species. The turbulence/radiation interaction is achieved via a new formulation. The resulting formulation is validated by comparison with experimental data on reacting supersonic axisymmetric jets. Results obtained for specific conditions indicate that the effect of chemical reaction on the turbulence is significant. Also, the radiative heat transfer is enhanced by the turbulence.

Chandrasekhar, R.↗

Modeling of SAR returns from a red pine stand

Bright P-band radar returns from red pine forests have been observed on synthetic aperture radar (SAR) images in Bangor, Maine. A plot of red pine trees was selected for the characterization and modeling to understand the cause of the high P-band returns. The red pine stand under study consisted of mature trees. Diameter at breast height (DBH) measurements were made to determine stand density as a function of tree diameter. Soil moisture and bulk density measurements were taken along with ground rough surface profiles. Detailed biomass measurements of the needles, shoots, branches, and trunks were also taken. These site statistics have been used in a distorted Born approximation model of the forest. Computations indicate that the direct-reflected or the double-bounce contributions from the ground are responsible for the high observed P-band returns for HH polarization.

Lang, R. H.↗

The Gaseous Content of the Universe at Zeta less than 1.6

Together with graduate student Hsiao-Wen Chen, I have measured and analyzed structural and morphological parameters of 38 galaxies in eight fields for which sensitive measurements of corresponding Ly(alpha) absorption toward background QSOs are available. These measurements are based on Wide Field Planetary Camera 2 (WFPC2) observations obtained with the Hubble Space Telescope (HST) and provide a first look at how the incidence and extent of tenuous gas around galaxies depends on galaxy luminosity, size, and morphological type and on geometry of the impact. The primary result of the analysis is that the amount of gas encountered along the line of sight depends on the galaxy impact parameter and B-band luminosity but does not depend strongly on the galaxy average surface brightness, disk-to-bulge ratio, or redshift. This result confirms and improves upon an anti-correlation between Ly(alpha) absorption equivalent width and galaxy impact parameter found previously. More importantly, this result provides the first quantitative means of relating statistics of faint galaxies to statistics of Ly(alpha) absorption systems. which we plan to exploit to constrain the luminosity function of galaxies beyond the realm of current surveys. Results have been submitted for publication and will greatly improve our statistical conclusions . Together with graduate student Noriaki Yahata. I have measured and classified spectral properties of over 1000 faint galaxies and stars obtained in our low-resolution spectroscopic survey. The goal of this project is two-fold: (1) to exhaustively characterize the spectral properties of all faint galaxies that comprise our current survey, and (2) to gain experience with our measurement and classification code. which ultimately will be used on a data base of 20,000 galaxies to be obtained with the Two-Degree Field (2df) spectrograph at the Anglo-Australian Telescope (AAT). The results will ultimately be used for many goals, but so far we have concentrated on using the results to make a binary classification of the galaxies (i.e. early type versus late type) and to then exploit the density-morphology relationship to obtain a crude density indicator. The primary result of the analysis is that the incidence and extent of tenuous gas around galaxies shows no strong preference for local galaxy environment, at least over the range of densities spanned by the current observations. Along similar lines, two instances of Ly(alpha) absorption lines that arise in groups or clusters were examined. Analysis demonstrates that some can produce corresponding absorption lines and that LY(alpha) absorption lines do not avoid a high-density environment. A new measure of the galaxy-absorber cross-correlation function defines the statistical criterion by which galaxies and absorber pairs are to be matched. I have identified a damped Ly(alpha) absorption system at redshift z equals approximately 0.16, the lowest redshift confirmed to date. The most important results of the analysis are learning that the metal abundances of the absorption system are less than 10 percent of the solar metal abundance and that the absorbing gas is not rotating with the galaxy disk.

Source record↗

Two-Dimensional Ffowcs Williams/Hawkings Equation Solver

FWH2D is a Fortran 90 computer program that solves a two-dimensional (2D) version of the equation, derived by J. E. Ffowcs Williams and D. L. Hawkings, for sound generated by turbulent flow. FWH2D was developed especially for estimating noise generated by airflows around such approximately 2D airframe components as slats. The user provides input data on fluctuations of pressure, density, and velocity on some surface. These data are combined with information about the geometry of the surface to calculate histories of thickness and loading terms. These histories are fast-Fourier-transformed into the frequency domain. For each frequency of interest and each observer position specified by the user, kernel functions are integrated over the surface by use of the trapezoidal rule to calculate a pressure signal. The resulting frequency-domain signals are inverse-fast-Fourier-transformed back into the time domain. The output of the code consists of the time- and frequency-domain representations of the pressure signals at the observer positions. Because of its approximate nature, FWH2D overpredicts the noise from a finite-length (3D) component. The advantage of FWH2D is that it requires a fraction of the computation time of a 3D Ffowcs Williams/Hawkings solver.

Lockard, David P.↗

Nonlocal screening in metal surfaces

Due to the effect of the nonuniform environment on the static screening of the Coulomb potential, the local-density approximation for the particle-hole interaction is found to be inadequate to determine the surface energy of simple metals. Use of the same set of single-particle states, and thus the same one-body density and the same work function, has eliminated the single-electron states in favor of the structure of the short-ranged correlations as the basis of this effect. A posteriori simplifications of the Fermi hypernetted-chain theory may be found to allow the same calculational accuracy with simpler computational tools.

Krotscheck, E.↗

Fanaroff-Riley I galaxies as the parent population of BL Lacertae objects. II - Optical constraints

The hypothesis is investigated that BL Lacertae objects are a subset of Fanaroff-Riley (FR) I galaxies dominated by beamed emission from a relativistic jet aligned with the line of sight. It is found that FR I galaxies, although of low radio luminosity, are bright optical galaxies. If the beaming hypothesis is correct, the optical luminosity function of BL Lac objects, including the effect of the host galaxy, can be approximated by a double power law with a given differential index. New estimates of the number counts of BL Lac objects are made which constrain their local number density for L(B) greater than about 6 x 10 to the 43rd ergs/s to be in the range 110-720/cu Gpc. The corresponding ratio between the total number densities of BL Lac objects and FR I galaxies is 0.003 to 0.02. This limits the Lorentz factor of the optical emission to the range 8-20, larger than that of the X-ray emitting plasma and roughly of the same order as the radio Lorentz factor.

Padovani, Paolo↗

Optimal design of an unsupervised adaptive classifier with unknown priors

An adaptive detection scheme for M hypotheses was analyzed. It was assumed that the probability density function under each hypothesis was known, and that the prior probabilities of the M hypotheses were unknown and sequentially estimated. Each observation vector was classified using the current estimate of the prior probabilities. Using a set of nonlinear transformations, and applying stochastic approximation theory, an optimally converging adaptive detection and estimation scheme was designed. The optimality of the scheme lies in the fact that convergence to the true prior probabilities is ensured, and that the asymptotic error variance is minimum, for the class of nonlinear transformations considered. An expression for the asymptotic mean square error variance of the scheme was also obtained.

Kazakos, D.↗