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 271 records · Page 15

Theoretical Prediction of Thermal Expansion Anisotropy for Y 2 Si 2 O 7 Environmental Barrier Coatings Using a Deep Neural Network Potential and Comparison to Experiment

Environmental barrier coatings (EBCs) are an enabling technology for silicon carbide (SiC)-based ceramic matrix composites (CMCs) in extreme environments such as gas turbine engines. However, development of new coating systems is hindered by the large design space and difficulty in predicting properties for these materials. Density Functional Theory (DFT) has successfully been used to model and predict some thermodynamic and thermo-mechanical properties of high-temperature ceramics for EBCs, although these calculations are challenging due to their high computational costs. In this work, we use machine learning to train a deep neural network potential (DNP) for Y 2 Si 2 O 7 , which is then applied to calculate thermodynamic and thermo-mechanical properties at near-DFT accuracy much faster and using less computational resources than DFT. We use this DNP to predict phonon-based thermodynamic properties of Y 2 Si 2 O 7 with good agreement to DFT and experiments. We also utilize the DNP to calculate the anisotropic, lattice direction-dependent coefficients of thermal expansion (CTEs) for Y 2 Si 2 O 7 . Molecular dynamics trajectories using the DNP correctly demonstrate accurate prediction of the anisotropy of the CTE in good agreement with diffraction experiments. In the future, this DNP could be applied to accelerate additional property calculations for Y 2 Si 2 O 7 compared to DFT or experiments.

rare earth silicates↗

Molecular Dynamics of a Water-Lipid Bilayer Interface

We present results of molecular dynamics simulations of a glycerol 1-monooleate bilayer in water. The total length of analyzed trajectories is 5ns. The calculated width of the bilayer agrees well with the experimentally measured value. The interior of the membrane is in a highly disordered fluid state. Atomic density profile, orientational and conformational distribution functions, and order parameters indicate that disorder increases toward the center of the bilayer. Analysis of out-of-plane thermal fluctuations of the bilayer surfaces occurring at the time scale of the present calculations reveals that the distribution of modes agrees with predictions of the capillary wave model. Fluctuations of both bilayer surfaces are uncorrelated, yielding Gaussian distribution of instantaneous widths of the membrane. Fluctuations of the width produce transient thinning defects in the bilayer which occasionally span almost half of the membrane. The leading mechanism of these fluctuations is the orientational and conformational motion of head groups rather than vertical motion of the whole molecules. Water considerably penetrates the head group region of the bilayer but not its hydrocarbon core. The total net excess dipole moment of the interfacial water points toward the aqueous phase, but the water polarization profile is non-monotonic. Both water and head groups significantly contribute to the surface potential across the interface. The calculated sign of the surface potential is in agreement with that from experimental measurements, but the value is markedly overestimated. The structural and electrical properties of the water-bilayer system are discussed in relation to membrane functions, in particular transport of ions and nonelectrolytes across membranes.

Wilson, Michael A.↗

Short-term frequency stability of the Rb-87 maser.

Measurements of the short-term stability of the Rb-87 maser as a function of the maser power output and the receiver cutoff frequency are reported. The experimental data are compared to theoretical results obtained from an approximate theory. In this theory the transfer function of the maser for thermal noise is derived, and the spectral density of the phase fluctuations is calculated. An analytical expression for the 'Allan variance' is also given. A comparison of the stability of the Rb-87 maser with existing frequency standards shows its superiority for averaging times less than 1 sec.

Tetu, M.↗

Probability density functions in turbulent channel flow

The probability density functions (pdf's) of the fluctuating velocity components, as well as their first and second derivatives, are calculated using data from the direct numerical simulations (DNS) of fully developed turbulent channel flow. It is observed that, beyond the buffer region, the pdf of each of these quantities is independent of the distance from the channel wall. It is further observed that, beyond the buffer region, the pdf's for all the first derivatives collapse onto a single universal curve and those of the second derivatives also collapse onto another universal curve, irrespective of the distance from the wall. The kinetic-energy dissipation rate exhibits log normal behavior.

Dinavahi, Surya P. G.↗

Calculations on the competition between association and reaction for C3H(+) + H2

The ab initio results presently obtained for the potential energy surface of competing associative and reactive ion-molecule processes of the reactants H2 and C3H(+) show that the linear ion is able to directly access the deep potential well of the propargyl ion H2CCCH(+). Using the minimum energy potential pathway and properties of the stationary point structures determined via the ab initio methods, calculation results are obtained for both the association-rate coefficient for C3H3(+) production as a function of density, and the branching ratio between the propargyl and cyclic structures of the ion.

Maluendes, Sergio A.↗

Energetic electrons in the magnetophere of Saturn

The energy spectra and angular distributions of electrons observed by Pioneer 11 as a function of radial distance in the inner magnetosphere of Saturn are reanalyzed and phase space densities are then calculated. The radial dependence of phase space density requires a distributed loss process. The loss is greatest in the region of the E ring (5.5 less than L less than 8.5) and is attributed to collisions with the ring particles in agreement with earlier work by Van Allen et al. (1980b). Quantitative analysis yields the following properties of the E ring: the particle radii are in the range of 4 x 10(exp -5) to 3.2 x 10(exp -4) cm and the thickness of the ring is approximately 3 R(sub s). Between the inner edge of the E ring (5.5 R(sub s)) and the outer edge of the A ring (2.3 R(sub s)) there are more energetic electrons than can be supplied by radial diffusion from an external source. Detailed calculations show that a cosmic ray albedo neutron decay (CRAND) source in the A and B rings is a plausible source for this excess. The radial diffusion coefficient required to explain the E ring absorption and CRAND source for electrons is 1 x 10(exp -12) greater than D(sub 0) greater than 3 x 10(exp -12) R(sub s)(exp 2)/s, assuming that D(sub LL) = D(sub 0)L(exp 3). As part of the reanalysis program, a method for the deconvolution of pitch angle distributions observed by simple detectors on a rotating spacecraft is developed. This process removes the instrumental response and rotational smear due to finite sampling periods and yields true angular distributions.

Randall, B. A.↗

Energetic electrons in the magnetosphere of Saturn

The energy spectra and angular distributions of electrons observed by Pioneer 11 as a function of radial distance in the inner magnetosphere of Saturn are reanalyzed and phase space densities are then calculated. The radial dependence of phase space density requires a distributed loss process. The loss is greatest in the region of the E ring (5.5 less than L less than 8.5) and is attributed to collisions with the ring particles in agreement with earlier work by Van Allen et al. (1980). Quantitative analysis yields the following properties of the E ring: the particle radii are in the range of 4 x 10(exp -5) to 3.2 x 10(exp -4) cm and the thickness of the ring is approximately 3 R(sub s). Between the inner edge of the E ring (5.5 R(sub s)) and the outer edge of the A ring (2.3 R(sub s)) there are more energetic electrons than can be supplied by radial diffusion from an external source. Detailed calculations show that a cosmic ray albedo neutron decay (CRAND) source in the A and B rings is a plausible source for this excess. The radial diffusion coefficient required to explain the E ring absorption and CRAND source for electrons is 1 x 10(exp -12) greater than D(sub 0) greater than 3 x 10(exp -12) R(exp 2, sub s)/s, assuming that D(sub LL) = D(sub 0)L(exp 3). As part of the reanalysis program, a method for the deconvolution of pitch angle distributions observed by simple detectors on a rotating spacecraft is developed. This process removes the instrumental response and rotational smear due to finite sampling periods and yields true angular distributions.

Randall, B. A.↗

Statistical Short-Range Guidance for Peak Wind Speed Forecasts at Edwards Air Force Base, CA

The peak winds near the surface are an important forecast element for Space Shuttle landings. As defined in the Shuttle Flight Rules (FRs), there are peak wind thresholds that cannot be exceeded in order to ensure the safety of the shuttle during landing operations. The National Weather Service Spaceflight Meteorology Group (SMG) is responsible for weather forecasts for all shuttle landings. They indicate peak winds are a challenging parameter to forecast. To alleviate the difficulty in making such wind forecasts, the Applied Meteorology Unit (AMTJ) developed a personal computer based graphical user interface (GUI) for displaying peak wind climatology and probabilities of exceeding peak-wind thresholds for the Shuttle Landing Facility (SLF) at Kennedy Space Center. However, the shuttle must land at Edwards Air Force Base (EAFB) in southern California when weather conditions at Kennedy Space Center in Florida are not acceptable, so SMG forecasters requested that a similar tool be developed for EAFB. Marshall Space Flight Center (MSFC) personnel archived and performed quality control of 2-minute average and 10-minute peak wind speeds at each tower adjacent to the main runway at EAFB from 1997- 2004. They calculated wind climatologies and probabilities of average peak wind occurrence based on the average speed. The climatologies were calculated for each tower and month, and were stratified by hour, direction, and direction/hour. For the probabilities of peak wind occurrence, MSFC calculated empirical and modeled probabilities of meeting or exceeding specific 10-minute peak wind speeds using probability density functions. The AMU obtained and reformatted the data into Microsoft Excel PivotTables, which allows users to display different values with point-click-drag techniques. The GUT was then created from the PivotTables using Visual Basic for Applications code. The GUI is run through a macro within Microsoft Excel and allows forecasters to quickly display and interpret peak wind climatology and likelihoods in a fast-paced operational environment. A summary of how the peak wind climatologies and probabilities were created and an overview of the GUT will be presented.

Dreher, Joseph↗

Analysis of Shockwave Radiation Data in Nitrogen

Data from a pure nitrogen test series in the Electric Arc Shock Tube Facility were previously reported for velocities spanning 6-12 km/s at a free-stream pressure of 0.2 Torr. This test series provides validation data for a range of physical phenomena to investigate, including vibrational relaxation, molecular radiation, nitrogen dissociation and ionization, and atomic radiation and ionization. This paper details analysis of data obtained at a nominal velocity of 10.3 km/s. The spectra are analyzed to extract temperatures and the densities of excited states as a function of position behind the shock. The effect of different methods for calculating state populations and ionization processes is assessed, as is a rigorous assessment of the atomic line lists, with both missing and extra lines identified.

Nitrogen↗

MIRO Computational Model

A computational model calculates the excitation of water rotational levels and emission-line spectra in a cometary coma with applications for the Micro-wave Instrument for Rosetta Orbiter (MIRO). MIRO is a millimeter-submillimeter spectrometer that will be used to study the nature of cometary nuclei, the physical processes of outgassing, and the formation of the head region of a comet (coma). The computational model is a means to interpret the data measured by MIRO. The model is based on the accelerated Monte Carlo method, which performs a random angular, spatial, and frequency sampling of the radiation field to calculate the local average intensity of the field. With the model, the water rotational level populations in the cometary coma and the line profiles for the emission from the water molecules as a function of cometary parameters (such as outgassing rate, gas temperature, and gas and electron density) and observation parameters (such as distance to the comet and beam width) are calculated.

Broderick, Daniel↗

Vibration-rotation line shifts for 1 sigma g + H2/V,J/-1S/0/ He computed via close coupling - Temperature dependence

The density shifting of vibration-rotation transitions of H2 perturbed by He was computed (as a function of temperature) with no adjustable parameters. The calculation was carried out using the framework of the impact theory of Baranger with S-matrix elements obtained via close coupling calculations which incorporated the ab initio H2-H2 system potential of Tsapline et al.(1977). Vibrational and rotational inelasticity were neglected in the calculations; nevertheless good agreement with experimental data was obtained, up to moderate temperatures, for the density shift. A much poorer comparison was obtained for the density broadening.

Hahne, G. E.↗

Prediction of crosslink density of solid propellant binders

A quantitative theory is outlined which allows calculation of crosslink density of solid propellant binders from a small number of predetermined parameters such as the binder composition, the functionality distributions of the ingredients, and the extent of the curing reaction. The parameter which is partly dependent on process conditions is the extent of reaction. The proposed theoretical model is verified by independent measurement of effective chain concentration and sol and gel fractions in simple compositions prepared from model compounds. The model is shown to correlate tensile data with composition in the case of urethane-cured polyether and certain solid propellants. A formula for the branching coefficient is provided according to which if one knows the functionality distributions of the ingredients and the corresponding equivalent weights and can measure or predict the extent of reaction, he can calculate the branching coefficient of such a system for any desired composition.

Marsh, H. E., Jr.↗

Frequency modulation television analysis: Threshold impulse analysis

A computer program is developed to calculate the FM threshold impulse rates as a function of the carrier-to-noise ratio for a specified FM system. The system parameters and a vector of 1024 integers, representing the probability density of the modulating voltage, are required as input parameters. The computer program is utilized to calculate threshold impulse rates for twenty-four sets of measured probability data supplied by NASA and for sinusoidal and Gaussian modulating waveforms. As a result of the analysis several conclusions are drawn: (1) The use of preemphasis in an FM television system improves the threshold by reducing the impulse rate. (2) Sinusoidal modulation produces a total impulse rate which is a practical upper bound for the impulse rates of TV signals providing the same peak deviations. (3) As the moment of the FM spectrum about the center frequency of the predetection filter increases, the impulse rate tends to increase. (4) A spectrum having an expected frequency above (below) the center frequency of the predetection filter produces a higher negative (positive) than positive (negative) impulse rate.

Hodge, W. H.↗

Mass flow in loop type coronal transients

Coronal transients having characteristics of a well-defined loop structure are examined, particularly with respect to temporal changes in the density and mass per unit length along the loop over periods of several days after the initial eruption. Measurements of mass distributions as a function of time are presented for eight transients, and one particular transient with a fairly simple configuration is investigated in more detail. Theoretical calculations are combined with the masses and densities derived from the observations to obtain estimates of the material flow in the transients; this flow is modeled on the assumption that magnetic forces drive and confine the loop. The flow field is found to be diverging everywhere, indicating that the density decreases in time. It is inferred that the transient legs are approximately in hydrostatic equilibrium and that most of the mass of the transient is lost from the sun during the initial phase.

Anzer, U.↗

Trajectory Software With Upper Atmosphere Model

The Trajectory Software Applications 6.0 for the Dec Alpha platform has an implementation of the Jacchia-Lineberry Upper Atmosphere Density Model used in the Mission Control Center for International Space Station support. Previous trajectory software required an upper atmosphere to support atmosphere drag calculations in the Mission Control Center. The Functional operation will differ depending on the end-use of the module. In general, the calling routine will use function-calling arguments to specify input to the processor. The atmosphere model will then compute and return atmospheric density at the time of interest.

Barrett, Charles↗

Measurement of the lunar photoelectron layer in the geomagnetic tail

Stable photoelectron fluxes (with energies between 40 and 200 eV) were observed by the Apollo 15 Charged Particle Lunar Environment Experiment (CPLEE). These observations, made in the magnetotail under near vacuum conditions, are compared with numerically calculated photoemission spectra to determine the approximate potential difference between ground and CPLEE'S apertures (a distance of 26 cm). Numerically calculated density and potential distributions are compared with the measured values to provide an estimate of the photoelectron yield function of the dust layer covering the moon.

Reasoner, D. L.↗

Brown dwarfs as dark galactic halos

The possibility that the dark matter in galactic halos can consist of brown dwarf stars is considered. The radiative signature for such halos consisting solely of brown dwarfs is calculated, and the allowed range of brown dwarf masses, the initial mass function (IMF), the stellar properties, and the density distribution of the galactic halo are discussed. The prediction emission from the halo is compared with existing observations. It is found that, for any IMF of brown dwarfs below the deuterium burning limit, brown dwarf halos are consistent with observations. Brown dwarf halos cannot, however, explain the recently observed near-IR background. It is shown that future satellite missions will either detect brown dwarf halos or place tight constraints on the allowed range of the IMF.

Adams, Fred C.↗