Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density Functional calculations”

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 325 records · Page 18

Synchrotron radiation as a probe of the inner magnetosphere of Jupiter

A short review is given on the characteristics of Jupiter's inner magnetosphere derived from radio observations in the decimetric wavelength range. A comparison of the data with sophisticated model calculations yields information on the magnetic field configuration and the electron distribution, its density, energy spectrum, and pitch angle dependence as a function of spatial coordinates. The latter information can be used to derive, e.g., the radial diffusion parameters plus the effects of the satellites, Jupiter's ring, and wave-particle interactions upon the electron distribution.

De Pater, I.↗

Gravitational microlensing of high-redshift supernovae by compact objects

An analysis of the effect of microlensing by a cosmologically dominant density of compact objects is performed, using high-redshift Type Ia supernovae (SN Ia's) as probes. The compact objects are modeled as a three-dimensional distribution of point masses, and Monte Carlo simulations are done to calculate the resulting amplification probability distributions for several column densities and cosmologies. By combining these distributions with the intrinsic SN Ia luminosity function and comparing with the results for a perfectly smooth universe, estimates are made of the number of supernovae that would need to be observed to confirm or rule out this lensing scenario. It is found that about 1000 SN Ia's with redshifts of z = 1 would be needed to perform this test, which is beyond what current searches can hope to accomplish. Observations of many fewer high-redshift supernovae, used merely as standard candles, appears a promising way of distinguishing between different cosmological models.

Rauch, Kevin P.↗

Statistical Evaluation and Improvement of Methods for Combining Random and Harmonic Loads

Structures in many environments experience both random and harmonic excitation. A variety of closed-form techniques has been used in the aerospace industry to combine the loads resulting from the two sources. The resulting combined loads are then used to design for both yield/ultimate strength and high- cycle fatigue capability. This Technical Publication examines the cumulative distribution percentiles obtained using each method by integrating the joint probability density function of the sine and random components. A new Microsoft Excel spreadsheet macro that links with the software program Mathematica to calculate the combined value corresponding to any desired percentile is then presented along with a curve tit to this value. Another Excel macro that calculates the combination using Monte Carlo simulation is shown. Unlike the traditional techniques. these methods quantify the calculated load value with a consistent percentile. Using either of the presented methods can be extremely valuable in probabilistic design, which requires a statistical characterization of the loading. Additionally, since the CDF at high probability levels is very flat, the design value is extremely sensitive to the predetermined percentile; therefore, applying the new techniques can substantially lower the design loading without losing any of the identified structural reliability.

Brown, A. M.↗

Computational Investigation of the Near-Field Plasma Plume in Ion-Ion Propulsion

A two-fluid numerical model of plasma flows was developed to investigate the plume of an ion-ion propulsion system. The densities of positive and negative ions, and the associated values of net charge, electric field, and electric potential were calculated as a function of time throughout the domain. The computational domain was chosen to be large enough (25 thruster diameters downstream of the exit plane) to allow for examining the neutralization of the plume. The resulting plasma electric potential and charge neutrality at the downstream end of the domain are shown and they indicate that it is possible to alternatively accelerate oppositely charged ions without the need for an electron-emitting neutralizer and without facing any electric potential hills that could cause stagnation. However, compared to existing literature on ion-ion plasma thrusters, the results from this simulation predict a longer length-scale for voltage decay.

Sankaran, K.↗

Inverse Compton Cooling in the Coronae of Simulated Black Hole Accretion Flows

We present a formulation for a local cooling function to be employed in the diffuse, hot corona region of 3D GRMHD simulations of accreting black holes. This new cooling function calculates the cooling rate due to inverse Compton scattering by considering the relevant microphysics in each cell in the corona and approximating the radiation energy density and Compton temperature there by integrating over the thermal seed photon flux from the disk surface. The method either assumes ion and electron tem- peratures are equal (1T), or calculates them separately (2T) using an instantaneous equilibrium approach predicated on the actual relevant rate equations (Coulomb and Compton). The method is shown to be consistent with a more detailed ray-tracing calculation where the bulk of the cooling occurs, but is substantially less costly to per- form. As an example, we apply these methods to a harm3d simulation of a 10M⊙, non-spinning black hole, accreting at nominally 1% the Eddington value. Both 1T and 2T approaches lead to increased radiative efficiency and a larger fraction of total cooling in the corona as compared to the original target-temperature cooling function used by harm3d, especially in the 1T case. Time-averaged post-processing reveals that the continuum spectral observations predicted from these simulations are qualitatively similar to actual X-ray binary data, especially so for the 1T approach which yields a harder power-law component (Γ = 2.25) compared to the 2T version (Γ = 2.53).

Brooks E. Kinch↗

Scalar dissipation and mixing in turbulent reacting flows

The mixing term that appears in probability density function-based formulations of the statistical behavior of turbulent reacting flows was studied with the use of direct numerical simulations of an irreversible, second-order, isothermal chemical reaction between initially segregated reactants in decaying, homogeneous turbulence. Three-dimensional, time-dependent pseudospectral calculations with 64-cubed Fourier modes were made as well as were Monte Carlo simulations of the joint concentration pdf equations for coalescence-dispersion (C/D) models. It was found that the development of the scalar microscales depends largely on the initial concentration distributions (both isotropic and nonisotropic distributions were used) and is relatively independent of the reaction rate constant. Joint concentration probability density functions were compared to predictions of D/D models that make use of parameters evaluated from the nonreacting or inert system.

Leonard, Andy D.↗

Ab initio Screening of Refractory Nitrides and Carbides for High Temperature Hydrogen Permeation Barriers

Density functional theory was used to screen eleven refractory materials – two pure metals, six nitrides, and three carbides–as high-temperature hydrogen permeation barriers to prevent hydrogen embrittlement. Activation energies were calculated for atomic hydrogen (H) diffusion into the first subsurface layer from the lowest energy surface of the high-temperature phase of each candidate material. The candidate barrier materials with the highest activation energies are h-BN, c-BN, HfN, and ZrN with predicted barriers of 3.25 eV, 3.23 eV, 3.14 eV, and 2.76 eV, respectively. Strain energies, Bader charges, and density of states were calculated for the diffusing H at the relaxed initial state and the transition state to provide insight into contributing factors to high energy barriers. The diffusing H atom in materials with the highest predicted barriers are protic. In addition, interstitial H atoms induce mid-gap states in the density of states of both BN polymorphs. The nitrogen retention of each nitride material at high temperatures was predicted using nitrogen vacancy formation energies with respect to gaseous nitrogen. Experimental evaluation of nitrogen retention in h-BN, ZrN, and TiN confirmed their resistance to nitrogen loss at 1773 K. However, of these nitrides, TiN is predicted to be the least stable. This work identifies multiple promising materials that are predicted to be effective hydrogen barriers at high temperatures and that are stable at temperatures above 2700 K, with BN predicted to perform best.

Density Functional Theory (DFT)↗

CNO abundances and temperature fluctuations in the Orion Nebula

The largest UV-optical data base to date of the Orion Nebula has been assembled, and the CNO abundances their distribution at a statistically significant number of positions to the north and west of the Trapezium and out to a distance of 5 arcmin has been examined. The forbidden O III temperature increases with distance from Theta1 Ori C while the forbidden N II temperature decreases with increasing distance. The density decreases as a function of increasing distance from Theta1 Ori C. The C/H, N/H, O/H, and C/O distributions are constant across the inner 5 arcmin of the northwestern quadrant of the nebula. The empirically determined mean C/H abundance calculated from UV lines is almost a factor of two lower in the number density than earlier estimates. Photoionization model calculations are consistent with a solar abundance for N/H and O/H.

Walter, Donald K.↗

Geoid anomalies in a dynamic earth

Potential and surface deformation Love numbers for internal loading have been calculated in order to obtain a dynamically consistent relationship between the geoid and the earth's response to internal buoyancy forces. These quantities depend on the depth and harmonic degree of loading, and can be integrated as Green functions to obtain the dynamic response due to an arbitrary distribution of internal density contrasts. Constructing a series of spherically symmetric, self-gravitating flow models for a variety of radial Newtonian viscosity variations and flow configurations, and calculating relaxation times for spherically symmetric viscous earth models, it is demonstrated that boundary deformation due to internal loading reaches its equilibrium value on the same time scale as postglacial rebound; this is much less time than the time scale for significant change in the convective flow pattern.

Richards, M. A.↗

Atomic Data and Spectral Line Intensities for Ne III

Electron impact collision strengths, energy levels, oscillator strengths, and spontaneous radiative decay rates are calculated for Ne III. The configurations used are 2s(sup 2)2p(sup 4), 2s2p(sup 5), 2p(sup 6), 2s(sup 2)2p(sup 3)3(sub s), 2s(sup 2)2p(sup 3)3p, and 2s(sup 2)2p(sup 3)3d giving rise to 86 fine-structure levels in intermediate coupling. Collision strengths are calculated at five incident energies, 5.2, 10, 15, 20 and 25 Ry. Excitation rate coefficients are calculated as a function of electron temperature by assuming a Maxwellian electron velocity distribution. Using the excitation rate coefficients and the radiative transition rates, statistical equilibrium equations for level populations are solved at electron densities covering the range of 10(exp 8) - 10(exp 14)/cc at an electron temperature of logTe/K = 5.0, corresponding to maximum abundance of Ne III. Relative spectral line intensities are calculated.

Landi, E.↗

Analysis of electromagnetic scattering from irregularly shaped, thin, metallic flat plates

This report describes an application of the method of moments to calculate the electromagnetic scattering from irregularly shaped, thin, metallic flat plates in free space. In the present technique, an irregularly shaped plate is enclosed by a rectangle on which the surface-current density is then expressed in terms of subdomain functions by dividing the rectangle into subsections. A shape function is introduced to ensure zero current outside the patch. The surface-current density is determined using the electric field integral equation (EFIE) approach in conjunction with the method of moments, and from a knowledge of the surface-current density, the electromagnetic scattering from a plate is calculated. Using this technique, the electromagnetic scattering from a hexagonal plate; an equilateral triangular plate; an equilateral triangular plate with a concentric, equilateral triangular hole and an inverted, equilateral triangular hole; and a diamond-shaped plate is computed and compared with the numerical results obtained by using the Electromagnetic Surface Patch (ESP) code developed by Ohio State University. The numerical results compare favorably with the measurements performed on these shapes in the Langley Experimental Test Range facility.

Deshpande, Manohar D.↗

Statistical time-dependent model for the interstellar gas

We present models for temperature and ionization structure of low, uniform-density (approximately 0.3 per cu cm) interstellar gas in a galactic disk which is exposed to soft X rays from supernova outbursts occurring randomly in space and time. The structure was calculated by computing the time record of temperature and ionization at a given point by Monte Carlo simulation. The calculation yields probability distribution functions for ionized fraction, temperature, and their various observable moments. These time-dependent models predict a bimodal temperature distribution of the gas that agrees with various observations. Cold regions in the low-density gas may have the appearance of clouds in 21-cm absorption. The time-dependent model, in contrast to the steady-state model, predicts large fluctuations in ionization rate and the existence of cold (approximately 30 K), ionized (ionized fraction equal to about 0.1) regions.

Gerola, H.↗

Mars - VLA observations of the northern hemisphere and the north polar region at wavelengths of 2 and 6 cm

Calculations based on 2- and 6-cm observations of Mars with the A configuration of the VLA have yielded a whole-disk effective dielectric constant of 2.34 + or - 0.05, implying a subsurface density of 1.24 + or - 0.11 g/cu cm at 2 cm, as well as 1.45 + or - 0.10 g/cu cm effective density and 2.70 + or - 0.10 dielectric constant at 6 cm. These parameters have also been estimated as a function of latitude over the 15 deg S - 60 deg N range; subsurface radio absorption length was estimated to be about 15 wavelengths at most of these latitudes. Most of the subsurface density calculations yielded results in the 1-2-g/cu cm range, implying that the subsurface is not very different from the surface observed by Viking and Mariner spacecraft; the decrease in correlation with depth is in keeping with slow variation of the subsurface in the near-subsurface region.

Rudy, Donald J.↗

The Solar Wind

Parker's model of a spherically expanding corona, the "solar wind," is compared with D. E. Blackwell's observations of the 1954 minimum equatorial corona. A significant discrepancy is found between the predicted and the observed electron densities at distances from the sun greater than 20 solar radii. Blackwell's data are found to be consistent with a model in which the corona expands mostly within a disk less than 25 solar radii thick, lying within the sun's equatorial plane. The thickness of the disk as a function of distance from the sun is qualitatively explained in terms of magnetic pressure. The solar wind is found to have a considerable effect on the lunar atmosphere. First, the calculated density of the lunar atmosphere is greatly reduced by collisions with protons in the solar wind. If the flux of particles in this wind has the conventional values ranging between 10(exp 10) to 10(exp 11) per sq cm-sec, the calculations yield a lunar pressure of 10(exp -13) atmosphere of argon, in agreement with the value predicted by Elsmore and Whitfield on the basis of observations on the occultation of radio stars. Second, following a suggestion by Gold, it was found that the collisions of solar-wind protons with the lunar surface produce an atmosphere of cold neutral hydrogen with a density of 10(exp 5) per cu cm at the lunar surface. The density falls off at greater distances in accordance with the inverse-square law. Estimates indicate that the interaction of solar particles with the neutral hydrogen will produce an extended lunar ionosphere with a density of the order of 400 protons/cu cm in the vicinity of the moon.

Herring, J. R.↗

High zonal harmonics of rapidly rotating planets

A new perturbation expansion is derived for the structure of rotating bodies in hydrostatic equilibrium. The method uses an expansion of the density on Legendre polynomial functions of angle, and can be developed analytically in a manner analogous to the standard level-surface perturbation theory. The new theory proceeds from a prescribed pressure-density relation rather than from a prescribed density distribution, and is both simpler and more physically transparent than the level-surface approach. High zonal harmonics are shown to arise via a transfer function involving derivatives of the interior sound velocity, and via mixing of multipole density components in the outer shell of the planet. Sample calculations for polytropic sequences are presented, as well as standard gravity models for Jupiter and Saturn. Mathematical subleties of the theory are discussed in an appendix.

Hubbard, W. B.↗

An analysis of auroral E region neutral winds based on incoherent scatter radar observations at Chatanika

The physical origin of auroral E-region neutral winds during geomagnetic disturbances is investigated by comparing incoherent-backscatter radar observations made at Chatanika, Alaska, during such disturbances with detailed theoretical calculations of neutral velocities for these conditions. In the analysis, an ion-drag collisional-coupling model is adopted as the primary driving force for the neutral winds. Electric fields and electron number densities observed by the radar facility are used in the calculations, and the coupled ion and neutral momentum equations are solved numerically as functions of altitude and time in the ion-drag approximation. The altitude structure of the velocity field is examined, and the detailed response of the neutral gas to the ion-drag mechanism is delineated. Weighted altitude averages of the calculated velocities are obtained, the accuracy of the procedure for deriving neutral winds from incoherent-scatter radar observations is evaluated, and the comparison is carried out. It is concluded that both ion drag and pressure gradients appear to be required to explain the neutral winds observed in the auroral E-region during geomagnetic disturbances.

Comfort, R. H.↗

Tearing at the dayside magnetopause

Initial observations of the dayside magnetopause thickness using magnetometers on board ISEE 1 (International Sun-Earth Explorer) and ISEE 2 indicate that this boundary may only be a few ion gyro-radii thick during periods of observed southward interplanetary magnetic field (IMF). Such a thin current layer can destabilize the collisionless tearing mode, a probable first step in the reconnection sequence. The linear growth rate is calculated assuming typical magnetopause parameters and it is found that it is a sensitive function of IMF orientation, magnetopause thickness, and electron number density. It is also shown that the finite size of the dayside magnetopause necessitates the existence of a two-dimensional tearing mode wave vector spectrum, a consequence of the inability of a thermal electron to maintain Landau resonance with the wave for a growth period. Implications regarding reconnection are discussed.

Quest, K. B.↗

Hexagonal and nematic phases of chains. I - Correlation functions. II - Phase transitions

The statistical mechanics of a system of semiflexible chains, which can represent polymer liquid crystals, long-chain biomolecules, stiff wormlike micelles, or columns of discotic liquid crystals, are examined. A continuum theory is used to calculate static correlation functions in the hexagonal and nematic phases. Two correlation functions are considered: (1) the structure factor which describes fluctuations in the density; and (2) the director fluctuation spectrum, which describes fluctuations in the local optical axis. In addition, a model is developed for the phase transitions of a system of infinitely long, semiflexible chains which interact through a steric, excluded-volume repulsion. The model yields generic phase diagrams in terms of pressure or density vs. persistence length or temperature.

Selinger, Jonathan V.↗