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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Galactic propagation models consistent with the cosmic ray lifetime derived from Be-10 measurements

Using a propagation calculation with energy dependent parameters, including the depletion of short pathlengths, and incorporating experimental nuclear excitation functions, the variation of the Be-10/Be9 ratio with the matter densities in two nested confinement regions is investigated. It is shown that there is no unique correspondence between a Be-10/Be9 measurement at low energy and the density of matter in the galaxy. Be-10/Be9 measurements at both low and high energy are needed to fully specify the matter densities.

Guzik, T. G.↗
NASA NTRS1990

Recognition and characterization of hierarchical interstellar structure. I - Correlation function

The problem of the quantitative description of multiscale structure in interstellar cloud complexes and gravitational collapse calculations is considered, emphasizing the recognition and characterization of hierarchical fragmentation structure. The response of the two-point correlation function to a variety of analytical models for density structure is discussed for simple clustering of pointlike clouds to more complex models involving clouds with a distribution of sizes and densities and hierarchical substructure. By expressing the density distribution as the superposition of individual clouds, it is shown that the correlation function generates two types of terms: those involving each cloud's density convolved with itself and those involving pairs of different clouds. Major distortion of the correlation function are introduced by the presence of any image features with size scales a significant fraction of the image size.

Houlahan, Padraig↗
NASA NTRS1996

Use of Density Functional Method to Study Molecular Vibrations

Use of a scale factor reduces the errors in computed vibrational frequencies (errors due to anharmonicities, basis set deficiencies, and approximate treatment of electron correlation). Additional errors come from modelling solids with small molecules. Nevertheless, we hope that the present calculations add to the understanding of surface IR measurements of the heat tiles.

Chong, Delano P.↗
NASA NTRS1974

The effect of ultrasonics on the strength properties of carbon steel processed by cold plastic deformation

A study was made of the influence of ultrasounds on the mechanical properties of OLT 35 carbon steel tubes cold-drawn on a plug ultrasonically activated by longitudinal waves. Experimental results indicate that: 1. The reduction in the values of the flow limit and tensile strength is proportional to the increase in acoustic energy introduced into the material subjected to deformation. 2. The diminution in influence of ultrasounds on tensile strength and flow rate that is due to an increased degree of deformation is explained by a reduction in specific density of the acoustic energy at the focus of deformation. 3. The relations calculated on the basis of the variation in the flow limit and tensile strength as a function of acoustic energy intensity was verified experimentally.

Atanasiu, N.↗
NASA NTRS1980

Simulation of spacecraft charging environments by monoenergetic beams

To simulate the distributed spectra of plasmas in space, the differential current density spectrum of the plasma was divided into a number of energy bands and the beam energy and current were calculated for each band to provide a piecewise reproduction of the distributed spectrum. Beam energies and current densities were chosen to match the velocity moments of the plasma distribution function. The velocity moments are averages related to physical quantities such as particle density, flux, pressure, and energy flux, and have been used extensively to characterize the measured properties of plasmas in space. Combinations of one, two, and three beams were found to match two to six velocity moments of Maxwellian distributions. A computational model was used to compare the charging of a spacecraft by plasmas with distributed spectra and by monoenergetic beams.

Halverson, W.↗
NASA NTRS1981

Experimental and Theoretical Studies of Interstellar Grains

Steady state vibrational populations of SiO and CO in dilute black body radiation fields were calculated as a function of total pressure, kinetic temperature and chemical composition of the gas. Approximate calculations for polyatomic molecules are presented. Vibrational disequilibrium becomes increasingly significant as total pressure and radiation density decrease. Many regions of postulated grain formation are found to be far from thermal equilibrium before the onset of nucleation. Calculations based upon classical nucleation theory or equilibrium thermodynamics are expected to be of dubious value in such regions. Laboratory measurements of the extinction of small iron and magnetite grains were made from 195 nm to 830 nm and found to be consistent with predictions based upon published optical constants. This implies that small iron particles are not responsible for the 220 nm interstellar extinction features. Additional measurements are discussed.

Nuth, J. A., III↗
NASA NTRS2006

Computer Modeling of Thermal Convection in Melts to Explain Glass Formation in Low Gravity and on Earth

Experiments conducted up to this time on glass forming melts in the low gravity environment of space show that glasses prepared in low-g are more chemically homogeneous and more resistant to crystallization than the comparable glasses prepared at 1-g on Earth. This result is somewhat surprising and opposite to the accepted concept on glass formation for a melt. A hypothesis based on "shear thinning" of a melt, a decrease in viscosity with increasing shear stress, is proposed as an explanation for the observed low-gravity results. This paper describes detailed simulation procedures to test the role of thermal convection in introducing shear stress in glass forming melts, using a lithium disilcate melt as a model. The simulation system in its idealized version consists of a cylinder that is heated at one end and cooled at the other with gravity acting in a transverse direction to the thermal gradient. The side wall of the cylinder is assumed to be insulating. The governing equations of motion and energy are solved using variable properties for viscosity (Arrehenius and non-Arrehenius behaviors) and density (constant and temperature dependent). Other parametric variables in the calculations include gravity level and gravity vector orientation. The shear stress in the system are then computed as a function of gravity from the calculated values of maximum melt velocity, and its effect on melt viscosity (shear thinning) is predicted. Also included and discussed are the modeling efforts related to other potential convective processes in glass forming melts and their possible effects on melt viscosity.

Ray, Chandra S.↗
NASA NTRS2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λc, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λc-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

Amanda Stricklan↗
NASA NTRS2016

Modified Spectral Fatigue Methods for S-N Curves With MIL-HDBK-5J Coefficients

The rainflow method is used for counting fatigue cycles from a stress response time history, where the fatigue cycles are stress-reversals. The rainflow method allows the application of Palmgren-Miner's rule in order to assess the fatigue life of a structure subject to complex loading. The fatigue damage may also be calculated from a stress response power spectral density (PSD) using the semi-empirical Dirlik, Single Moment, Zhao-Baker and other spectral methods. These methods effectively assume that the PSD has a corresponding time history which is stationary with a normal distribution. This paper shows how the probability density function for rainflow stress cycles can be extracted from each of the spectral methods. This extraction allows for the application of the MIL-HDBK-5J fatigue coefficients in the cumulative damage summation. A numerical example is given in this paper for the stress response of a beam undergoing random base excitation, where the excitation is applied separately by a time history and by its corresponding PSD. The fatigue calculation is performed in the time domain, as well as in the frequency domain via the modified spectral methods. The result comparison shows that the modified spectral methods give comparable results to the time domain rainflow counting method.

Irvine, Tom↗
NASA NTRS1995

Numerical investigations of low-density nozzle flow by solving the Boltzmann equation

A two-dimensional finite-difference code to solve the BGK-Boltzmann equation has been developed. The solution procedure consists of three steps: (1) transforming the BGK-Boltzmann equation into two simultaneous partial differential equations by taking moments of the distribution function with respect to the molecular velocity u(sub z), with weighting factors 1 and u(sub z)(sup 2); (2) solving the transformed equations in the physical space based on the time-marching technique and the four-stage Runge-Kutta time integration, for a given discrete-ordinate. The Roe's second-order upwind difference scheme is used to discretize the convective terms and the collision terms are treated as source terms; and (3) using the newly calculated distribution functions at each point in the physical space to calculate the macroscopic flow parameters by the modified Gaussian quadrature formula. Repeating steps 2 and 3, the time-marching procedure stops when the convergent criteria is reached. A low-density nozzle flow field has been calculated by this newly developed code. The BGK Boltzmann solution and experimental data show excellent agreement. It demonstrated that numerical solutions of the BGK-Boltzmann equation are ready to be experimentally validated.

Deng, Zheng-Tao↗
NASA NTRS1997

Analysis of Stub Loaded Microstrip Patch Antennas

A microstrip patch antenna fed by a coaxial probe and reactively loaded by a open circuited microstrip line has been used previously to produce circular polarization and also as a building block for a series fed microstrip patch array. Rectangular and circular patch antennas loaded with a microstrip stub were previously analyzed using the generalized Thevenin theorem. In the Thevenin theorem approach, the mutual coupling between the patch current and the surface current on the stub was not taken into account. Also, the Thevenin theorem approach neglects continuity of current at the patch-stub junction. The approach in this present paper includes the coupling between the patch and stub currents as well as continuity at the patch-stub junction. The input impedance for a stub loaded microstrip patch is calculated by the general planar dielectric dyadic Green's function approach in the spectral domain, as was initiated much earlier and has been extensively expanded upon and utilized successfully throughout the literature for microstrip antenna configurations. Using the spectral domain dyadic Green s function derived earlier with the electric field integral equation (EFIE), the problem is formulated by using entire domain basis functions to represent the surface current densities on the patch, the loading stub and the attachment mode at the junction. Galerkin's procedure is used to reduce the EFIE to a matrix equation, which is then solved to obtain the amplitudes of the surface currents. These surface currents are then used for calculating the input impedance of stub loaded rectangular and circular microstrip patches. Numerical results are compared with measured results and with previous results calculated by the Thevenin's theorem approach.

Deshpande, M. D.↗
NASA NTRS1999

Cloud Liquid Water, Mean Droplet Radius and Number Density Measurements Using a Raman Lidar

A new technique for measuring cloud liquid water, mean droplet radius and droplet number density is outlined. The technique is based on simultaneously measuring Raman and Mie scattering from cloud liquid droplets using a Raman lidar. Laboratory experiments on liquid micro-spheres have shown that the intensity of Raman scattering is proportional to the amount of liquid present in the spheres. This fact is used as a constraint on calculated Mie intensity assuming a gamma function particle size distribution. The resulting retrieval technique is shown to give stable solutions with no false minima. It is tested using Raman lidar data where the liquid water signal was seen as an enhancement to the water vapor signal. The general relationship of retrieved average radius and number density is consistent with traditional cloud physics models. Sensitivity to the assumed maximum cloud liquid water amount and the water vapor mixing ratio calibration are tested. Improvements to the technique are suggested.

Whiteman, David N.↗
NASA NTRS1982

Thermophysical properties of germanium for thermal analysis of growth from the melt

The thermal diffusivity of Ge has been measured over a temperature range from 300 C to 1010 C which includes values for the melt. Specific heat has been measured from room temperature to 727 C. Thermal conductivity has been calculated over the same temperature range as the diffusivity measurements. These data are reported along with the best values from the literature for the other parameters which are required to calculate the temperature and convective fields for the growth of germanium by the Bridgman method. These parameters include the specific heat, the viscosity, the emissivity, and the density as a function of temperature.

Crouch, R. K.↗
NASA NTRS1984

Line emission from H II blister models

Numerical techniques to calculate the thermal and geometric properties of line emission from H II 'blister' regions are presented. It is assumed that the density distributions of the H II regions are a function of two dimensions, with rotational symmetry specifying the shape in three-dimensions. The thermal and ionization equilibrium equations of the problem are solved by spherical modeling, and a spherical sector approximation is used to simplify the three-dimensional treatment of diffuse ionizing radiation. The global properties of H II 'blister' regions near the edges of a molecular cloud are simulated by means of the geometry/density distribution, and the results are compared with observational data. It is shown that there is a monotonic increase of peak surface brightness from the i = 0 deg (pole-on) observational position to the i = 90 deg (edge-on) position. The enhancement of the line peak intensity from the edge-on to the pole-on positions is found to depend on the density, stratification, ionization, and electron temperature weighting. It is found that as i increases, the position of peak line brightness of the lower excitation species is displaced to the high-density side of the high excitation species.

Rubin, R. H.↗
NASA NTRS2010

Packing Optimization of an Intentionally Stratified Sorbent Bed Containing Dissimilar Media Types

The Fire Cartridge is a packed bed air filter with two different and separate layers of media designed to provide respiratory protection from combustion products after a fire event on the International Space Station (ISS). The first layer of media is a carbon monoxide catalyst made from gold nanoparticles dispersed on iron oxide. The second layer of media is universal carbon, commonly used in commercial respirator filters. Each layer must be optimally packed to effectively remove contaminants from the air. Optimal packing is achieved by vibratory agitations. However, if post-packing movement of the media within the cartridge occurs, mixing of the bed layers, air voids, and channeling could cause preferential air flow and allow contaminants to pass. Several iterations of prototype fire cartridges were developed to reduce post-packing movement of the media within each layer (settling), and to prevent mixing of the two media types. Both types of movement of the media contribute to decreased fire cartridge performance. Each iteration of the fire cartridge design was tested to demonstrate mechanical loads required to cause detrimental movement within the bed, and resulting level of functionality of the media beds after movement was detected. In order to optimally pack each layer, vertical, horizontal, and orbital agitations were tested and a final packed bulk density was calculated for each method. Packed bulk density must be calculated for each lot of catalyst to accommodate variations in particle size, shape, and density. In addition, a physical divider sheet between each type of media was added within the fire cartridge design to further inhibit intermixing of the bed layers.

Kidd, Jessica↗
NASA NTRS2013

Current Flow and Pair Creation at Low Altitude in Rotation-Powered Pulsars' Force-Free Magnetospheres: Space Charge Limited Flow

We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of rotation-powered pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We apply a new 1D hybrid plasma simulation code to the dynamical problem, using Particle-in-Cell methods for the dynamics of the charged particles, including a determination of the collective electrostatic fluctuations in the plasma, combined with a Monte Carlo treatment of the high-energy gamma-rays that mediate the formation of the electron-positron pairs.We assume the electric current flowing through the pair creation zone is fixed by the much higher inductance magnetosphere, and adopt the results of force-free magnetosphere models to provide the currents which must be carried by the accelerator. The models are spatially one dimensional, and designed to explore the physics, although of practical relevance to young, high-voltage pulsars. We observe novel behaviour (a) When the current density j is less than the Goldreich-Julian value (0 < j/j(sub GJ) < 1), space charge limited acceleration of the current carrying beam is mild, with the full Goldreich-Julian charge density comprising the charge densities of the beam and a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are of the order of mc(sup 2)/e, and pair creation is absent. (b) When the current density exceeds the Goldreich-Julian value (j/j(sub GJ) > 1), the system develops high voltage drops (TV or greater), causing emission of curvature gamma-rays and intense bursts of pair creation. The bursts exhibit limit cycle behaviour, with characteristic time-scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). (c) In return current regions, where j/j(sub GJ) < 0, the system develops similar bursts of pair creation. These discharges are similar to those encountered in previous calculations by Timokhin of pair creation when the surface has a high work function and cannot freely emit charge. In cases (b) and (c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady (stationary with small fluctuations in the rotating frame), finding that such steady flows can occupy only a small fraction of the current density parameter space exhibited by the force-free magnetospheric model. The generic polar flow dynamics and pair creation are strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage (with pair creation, when the voltage drops are sufficiently large; without, when they are small) as a function of the applied current.

Timokhin, A. N.↗
NASA NTRS1967

Subsonic and Supersonic Flutter Analysis of a Highly Tapered Swept-Wing Planform, Including Effects of Density Variation and Finite Wing Thickness, and Comparison with Experiments

The flutter characteristics of several wings with an aspect-ratio of 4.0, a taper ratio of 0.2, and a quarter-chord sweepback of 45 deg. have been investigated analytically for Mach numbers up to 2.0. The calculations were based on the modified-strip-analysis method, the subsonic-kernel-function method, piston theory, and quasi-steady second-order theory. Results of t h e analysis and comparisons with experiment indicated that: (1) Flutter speeds were accurately predicted by the modified strip analysis, although accuracy at t h e highest Mach numbers required the use of nonlinear aerodynamic theory (which accounts for effects of wing thickness) for the calculation of the aerodynamic parameters. (2) An abrupt increase of flutter-speed coefficient with increasing Mach number, observed experimentally in the transonic range, was also indicated by the modified strip analysis. (3) In the low supersonic range for some densities, a discontinuous variation of flutter frequency with Mach number was indicated by the modified strip analysis. An abrupt change of frequency appeared experimentally in the transonic range. (4) Differences in flutter-speed-coefficient levels obtained from tests at low supersonic Mach numbers in two wind tunnels were also predicted by the modified strip analysis and were shown to be caused primarily by differences in mass ratio. (5) Flutter speeds calculated by the subsonic-kernel-function method were in good agreement with experiment and with the results of the modified strip analysis. (6) Flutter speed obtained from piston theory and from quasi-steady second-order theory were higher than experimental values by at least 38 percent.

Yates, Carson, Jr.↗
NASA NTRS1993

Relaxation of molecular velocities and spectral line profiles of astrophysical masers

The rates at which the velocities of water molecules relax to a Maxwellian distribution are calculated under conditions that are representative of those in astrophysical water masers. Populations are calculated as a function of velocity for the lowest 40 rotational energy levels of ortho-water in order to make a detailed determination of the velocity relaxation rate. Relaxation rates for the 22 GHz maser at H2 densities of 10 exp 9/cu cm are 2.0/s and 5.1/s for temperatures of 400 K and 1000 K, respectively. Under these and other conditions, as well as for other maser transitions, the relaxation rate far exceeds the loss rate - usually by more than a factor of 10. Thus maser lines can remain narrow up to intensities that are at least a factor of 10 greater than the saturation intensity, and to path lengths that are at least a factor of 10 greater than the path length at which saturation occurs. Relaxation is dominated by the effects of trapped infrared radiation, though elastic and inelastic collisions are also included in the calculations. These relaxation rates, together with the observed line breadths of the 22 GHz masers, provide limits on maser luminosities which indicate that the brightest 22 GHz masers are highly beamed, with solid angles smaller than 10 exp -5 to 10 exp -4 Sr.

Anderson, Nels↗