Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional calculation”

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 163 records · Page 9

The correlation function of galaxy ellipticities produced by gravitational lensing

The correlation of galaxy ellipticities produced by gravitational lensing is calculated as a function of the power spectrum of density fluctuations in the universe by generalizing an analytical method developed by Gunn (1967). The method is applied to a model where identical objects with spherically symmetric density profiles are randomly laid down in space, and to the cold dark matter model. The possibility of detecting this correlation is discussed. Although an ellipticity correlation can also be caused by an intrinsic alignment of the axes of galaxies belonging to a cluster or a supercluster, a method is suggested by which one type of correlation can be distinguished from another. The advantage of this ellipticity correlation is that it is one of the few astronomical observations that can directly probe large-scale mass fluctuations in the universe.

Miralda-Escude, Jordi↗

On the Reaction of FNO2 with CH3, t-butyl, and C13H21

Theoretical studies are reported for the reaction of FN02 with the radicals CH3, t-butyl, and C13H21, which are templates for the radical site of a hydrogenated diamond surface. All structures axe fully optimized using density functional theory (DFT) based on the B3LYP functional. Calibration calculations axe performed for CH3 + FNO2 using the coupled cluster approach, the the internally contracted multireference configuration. interaction method, and second order perturbation theory based upon complete active space SCF (CASSCF) reference wave function. These calibration calculations support the B3LYP approach for the calculation of bond energies, but show the B3LYP barrier is too low. Combining the calibration calculations with the larger clusters yields our best estimate of a barrier of about 10 kcal/mol. for the reaction of FNO2 with a radical site on hydrogenated diamond.

Thuemmel, H. T.↗

The effect of central starbursts on the interstellar medium of dwarf galaxies

Major starburst events can last tens of millions of years, and in the process they can deposit significant amounts of energy into the surrounding interstellar medium. This energy from supernova and stellar winds imparts enough momentum to the interstellar medium (ISM) that portions of the ISM can become unbound and leave the parent galaxy, taking the metal-enriched stellar debris along. In dwarf galaxies, starbursts can produce enough total energy to unbind most or all of the ambient ISM. Whether this actually occurs is a strong function of the ellipticity of the ISM distribution, with flat disks and spheres being the limiting cases. We calculate whether 'blow out' along the symmetry axis of 'blow away' of the entire ISM occurs during a central starburst in dwarf galaxies as a function of galactic mass, starburst energy, ISM density, and ISM ellipticity. The calculations cover a range of 10(exp 7) to 10(exp 9) solar mass for dwarf galaxies and include 'normal' galaxies of 10(exp 11) solar mass as well. No massive dark matter halos are assumed to be present. We find that for physically reasonable values of total ISM mass and starburst energy a blow out along the symmetry axis occurs in the majority of cases, though a significant fraction of small dwarf galaxies can lose most of their ISM. As no dark matter halos or clumpy ISM distributions are included, it is apparent that the ISM in most dwarf galaxies may be generally resistant to significant disruption by a central starburst event. The effects of this range of behavi or on the metallicities that would be observed in these galaxies is discussed.

De Young, David S.↗

Theoretical Studies of Liquid He-4 Near the Superfluid Transition

We performed theoretical studies of liquid helium by applying state of the art simulation and finite-size scaling techniques. We calculated universal scaling functions for the specific heat and superfluid density for various confining geometries relevant for experiments such as the confined helium experiment and other ground based studies. We also studied microscopically how the substrate imposes a boundary condition on the superfluid order parameter as the superfluid film grows layer by layer. Using path-integral Monte Carlo, a quantum Monte Carlo simulation method, we investigated the rich phase diagram of helium monolayer, bilayer and multilayer on a substrate such as graphite. We find excellent agreement with the experimental results using no free parameters. Finally, we carried out preliminary calculations of transport coefficients such as the thermal conductivity for bulk or confined helium systems and of their scaling properties. All our studies provide theoretical support for various experimental studies in microgravity.

Manousakis, Efstratios↗

Probabilistic material degradation model for aerospace materials subjected to high temperature, mechanical and thermal fatigue, and creep

A probabilistic general material strength degradation model has been developed for structural components of aerospace propulsion systems subjected to diverse random effects. The model has been implemented in two FORTRAN programs, PROMISS (Probabilistic Material Strength Simulator) and PROMISC (Probabilistic Material Strength Calibrator). PROMISS calculates the random lifetime strength of an aerospace propulsion component due to as many as eighteen diverse random effects. Results are presented in the form of probability density functions and cumulative distribution functions of lifetime strength. PROMISC calibrates the model by calculating the values of empirical material constants.

Boyce, L.↗

Auger Spectroscopy of Hydrogenated Diamond Surfaces

An energy shift and a change of the line shape of the carbon core-valence-valence Auger spectra are observed for diamond surfaces after their exposure to an electron beam, or annealing at temperatures higher then 950 C. The effect is studied for both natural diamond crystals and chemical-vapor-deposited diamond films. A theoretical model is proposed for Auger spectra of hydrogenated diamond surfaces. The observed changes of the carbon Auger line shape are shown to be related to the redistribution of the valence-band local density of states caused by the hydrogen desorption from the surface. One-electron calculation of Auger spectra of diamond surfaces with various hydrogen coverages are presented. They are based on self-consistent wave functions and matrix elements calculated in the framework of the local-density approximation and the self-consistent linear muffin-tin orbital method with static core-hole effects taken into account. The major features of experimental spectra are explained.

Krainsky, I. L.↗

Photodissociation rates of molecular oxygen in the mesosphere and lower thermosphere.

Use of recently obtained band oscillator strengths and predissociation line widths, for the Schumann-Runge band system to calculate the photodissociation rate of molecular oxygen as a function of column density and temperature by solar radiation between 1025 and 2500 A. Data are presented from which the mean transmission and photodissociation rates for given wavelength intervals between 1750 and 2025 A can be calculated for column densities between 10 to the 17th and 10 to the 23rd molecules/sq cm and for temperatures between 150 and 300 K. Also, the production of vibrationally excited oxygen molecules as a result of fluorescence in the Schumann-Runge bands is discussed.

Hudson, R. D.↗

Dislocation dynamics and the viscoplastic buckling of dendritic web type silicon ribbon

The effect of dendrites (reinforced edges) on the residual stresses, dislocation densities and buckling behavior during growth of web type silicon ribbon is studied. A viscoplastic material response function (Haasen-Sumino model) is used to calculate the stresses and the disloction density at each point in the silicon ribbon. In addition, the role of dendrites on the viscoplastic buckling behavior of the ribbon is investigated. The critical thicknesses, the corresponding deflection shapes and lateral deflection speeds are calculated. These results are then compared with similar data obtained for flat plates.

Tsai, C. T.↗

Optimized Vertex Method and Hybrid Reliability

A method of calculating the fuzzy response of a system is presented. This method, called the Optimized Vertex Method (OVM), is based upon the vertex method but requires considerably fewer function evaluations. The method is demonstrated by calculating the response membership function of strain-energy release rate for a bonded joint with a crack. The possibility of failure of the bonded joint was determined over a range of loads. After completing the possibilistic analysis, the possibilistic (fuzzy) membership functions were transformed to probability density functions and the probability of failure of the bonded joint was calculated. This approach is called a possibility-based hybrid reliability assessment. The possibility and probability of failure are presented and compared to a Monte Carlo Simulation (MCS) of the bonded joint.

Smith, Steven A.↗

A Preliminary Model of the Venus Atmosphere

“Standard’ and extreme model atmospheres of the planet Venus are calculated and presented, with temperature, density, and height given as functions of pressure, and with pressure, temperature, and density given as functions of height. Until more recent data are interpreted, the models are considered preliminary.

planetary entry↗

Plasma refractive effects in HCN lasers

Plasma refractive index effects in pulsed HCN lasers, calculating stable cavity mirror curvatures constriction as function of electron density by ray matrix approach

Mccaul, B. W.↗

Alternative communication network designs for an operational Plato 4 CAI system

The cost of alternative communications networks for the dissemination of PLATO IV computer-aided instruction (CAI) was studied. Four communication techniques are compared: leased telephone lines, satellite communication, UHF TV, and low-power microwave radio. For each network design, costs per student contact hour are computed. These costs are derived as functions of student population density, a parameter which can be calculated from census data for one potential market for CAI, the public primary and secondary schools. Calculating costs in this way allows one to determine which of the four communications alternatives can serve this market least expensively for any given area in the U.S. The analysis indicates that radio distribution techniques are cost optimum over a wide range of conditions.

Mobley, R. E., Jr.↗

Response of space shuttle insulation panels to acoustic noise pressure

The response of reusable space shuttle insulation panels to random acoustic pressure fields are studied. The basic analytical approach in formulating the governing equations of motion uses a Rayleigh-Ritz technique. The input pressure field is modeled as a stationary Gaussian random process for which the cross-spectral density function is known empirically from experimental measurements. The response calculations are performed in both frequency and time domain.

Vaicaitis, R.↗

Analysis of data from NASA B-57B gust gradient program

Statistical analysis of the turbulence measured in flight 6 of the NASA B-57B over Denver, Colorado, from July 7 to July 23, 1982 included the calculations of average turbulence parameters, integral length scales, probability density functions, single point autocorrelation coefficients, two point autocorrelation coefficients, normalized autospectra, normalized two point autospectra, and two point cross sectra for gust velocities. The single point autocorrelation coefficients were compared with the theoretical model developed by von Karman. Theoretical analyses were developed which address the effects spanwise gust distributions, using two point spatial turbulence correlations.

Frost, W.↗

Comparisons of a Three-Dimensional, Full Navier Stokes Computer Model with High Mach Number Combuster Test Data

Comparisons between scramjet combustor data and a three-dimensional full Navier-Stokes calculation have been made to verify and substantiate computational fluid dynamics (CFD) codes and application procedures. High Mach number scramjet combustor development will rely heavily on CFD applications to provide wind tunnel-equivalent data of quality sufficient to design, build and fly hypersonic aircraft. Therefore. detailed comparisons between CFD results and test data are imperative. An experimental case is presented, for which combustor wall static pressures were measured and flow-fieid interferograms were obtained. A computer model was done of the experiment, and counterpart parameters are compared with experiment. The experiment involved a subscale combustor designed and fabricated for the National Aero-Space Plane Program, and tested in the Calspan Corporation 96" hypersonic shock tunnel. The combustor inlet ramp was inclined at a 20 angle to the shock tunnel nozzle axis, and resulting combustor entrance flow conditions simulated freestream M=10. The combustor body and cowl walls were instrumented with static pressure transducers, and the combustor lateral walls contained windows through which flowfield holographic interferograms were obtained. The CFD calculation involved a three-dimensional time-averaged full Navier-Stokes code applied to the axial flow segment containing fuel injection and combustion. The full Navier-Stokes approach allowed for mixed supersonic and subsonic flow, downstream-upstream communication in subsonic flow regions, and effects of adverse pressure gradients. The code included hydrogen-air chemistry in the combustor segment which begins near fuel injection and continues through combustor exhaust. Combustor ramp and inlet segments on the combustor lateral centerline were modelled as two dimensional. Comparisons to be shown include calculated versus measured wall static pressures as functions of axial flow coordinate, and calculated path-averaged density contours versus an holographic Interferogram.

Watkins, William B.↗

Partial Model of Insulator/Insulator Contact Charging

Two papers present a two-phase equilibrium model that partly explains insulator/ insulator contact charging. In this model, a vapor of ions within a gas is in equilibrium with a submonolayer of ions of the same species that have been adsorbed on the surface of an insulator. The surface is modeled as having localized states, each with a certain energy of adsorption for an ion. In an earlier version of the model described in the first paper, the ions do not interact with each other. Using the grand canonical ensemble, the chemical potentials of both vapor and absorbed phases are derived and equated to determine the vapor pressure. If a charge is assigned to the vapor particles (in particular, if single ionization is assumed), then the surface charge density associated with adsorbed ions can be calculated as a function of pressure. In a later version of the model presented in the second paper, the submodel of the vapor phase is extended to include electrostatic interactions between vapor ions and adsorbed ones as well as the screening effect, at a given distance from the surface, of ions closer to the surface. Theoretical values of this model closely match preliminary experimental data on the discharge of insulators as a function of pressure.

Hogue, Michael↗

Calculating Free Energies Using Average Force

A new, general formula that connects the derivatives of the free energy along the selected, generalized coordinates of the system with the instantaneous force acting on these coordinates is derived. The instantaneous force is defined as the force acting on the coordinate of interest so that when it is subtracted from the equations of motion the acceleration along this coordinate is zero. The formula applies to simulations in which the selected coordinates are either unconstrained or constrained to fixed values. It is shown that in the latter case the formula reduces to the expression previously derived by den Otter and Briels. If simulations are carried out without constraining the coordinates of interest, the formula leads to a new method for calculating the free energy changes along these coordinates. This method is tested in two examples - rotation around the C-C bond of 1,2-dichloroethane immersed in water and transfer of fluoromethane across the water-hexane interface. The calculated free energies are compared with those obtained by two commonly used methods. One of them relies on determining the probability density function of finding the system at different values of the selected coordinate and the other requires calculating the average force at discrete locations along this coordinate in a series of constrained simulations. The free energies calculated by these three methods are in excellent agreement. The relative advantages of each method are discussed.

Darve, Eric↗

A model of coronal holes

It has been noted that coronal holes appear to be associated with regions of diverging magnetic field in the corona. We set out to test the hypothesis that coronal holes may be caused by an increased flow of energy into the solar wind resulting directly from this diverging magnetic field pattern. Simple models were devised to approximate the energy flow down into the transition region and up into the solar wind as a function of the temperature, density, and rate of field line divergence in the corona. By assuming the rate of mechanical energy influx into the corona to be constant, it was then possible to solve numerically for the coronal temperature and density as a function of the rate of field line divergence. The results of these calculations demonstrate that a diverging field pattern can, indeed, bring about reductions in the temperature and density at the base of the corona comparable to those observed in coronal holes.

Adams, W. M.↗