Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY FUNCTIONAL METHOD”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Cosmic temperature fluctuations from two years of COBE differential microwave radiometers observations

The first two years of Cosmic Background Explorer (COBE) Differential Microwave Radiometers (DMR) observations of the cosmic microwave background (CMB) anisotropy are analyzed and compared with our previously published first year results. The results are consistent, but the addition of the second year of data increases the precision and accuracy detected CMB temperature fluctuations. The 2 yr 53 GHz data are characterized by rms temperature fluctuations of (delta-T)(sub rms) (7 deg) = 44 +/- 7 micro-K and (delta-T)(sub rms) (10 deg) = 30.5 +/- 2.7 micro-K at 7 deg and 10 deg angular resolution, respectively. The 53 x 90 GHz cross-correlation amplitude at zero lag is C(0)(sup 1/2) = 36 +/- 5 micro-K (68% CL) for the unsmoothed (7 deg resolution) DMR data. We perform a likelihood analysis of the cross-correlation function, with Monte Carlo simulations to infer biases of the method, for a power-law model of initial density fluctuations, P(k) proportional to R(exp n). The Monte Carlo simulations indicate that derived estimates of n are biased by +0.11 +/- 0.01, while the subset of simulations with a low quadrupole (as observed) indicate a bias of +0.31+/- 0.04. Derived values for 68% confidence intervals are given corrected (and not corrected) for our estimated biases. Including the quadrupole anisotropy, the most likely quadrupole-normalized amplitude is Q(sub rms-PS) = 14.3(sup + 5.2 sub -3.3) micro-K (12.8(sup + 5.2 sub -3.3) micro-K0 with a spectral index n = 1.42(sup + 0.49 sub -0.55)(n = 1.53(sup + 0.49 sub -0.55). With n fixed to 1.0 the most likely amplitude is 18.2 +/- 11.5 micro-K (17.4 +/- 1.5 micro-K). The marginal likelihood of n is 1.42 +/- 0.37 (1.53 +/- 0.37). Excluding the quadrupole anisotropy, the most likely quadrupole-normalized amplitude is Q(sub rms-PS) = 17.4(sup + 7.5 sub -5.2) micro-K (15.8(sup + 7.5 sub -5.2) micro-K) with a spectral index n = 1.11(sup + 0.60 sub -0.55) (n = 1.22(sup + 0.60 sub -0.55). With n fixed to 1.0 the most likely amplitude is 18.6 +/- 1.6 micro-K (18.2 +/- 1.6 micro-K). The marginal likelihood of n is 1.11 +/- 0.40 (1.22 +/- 0.40). Our best estimate of the dipole from the 2 yr DMR data is 3.363 +/- 0.024 mK toward Galactic coordinates (l, b) = (264.4 deg +/- 0.2 deg, 48.1 deg +/- 0.4 deg), and our best estimate of the rms quadrupole amplitude in our sky is 6 +/- 3 micro-K (68% CL).

Bennett, C. L.↗
NASA NTRS2004

Probabilistic Aeroelastic Analysis of Turbomachinery Components

A probabilistic approach is described for aeroelastic analysis of turbomachinery blade rows. Blade rows with subsonic flow and blade rows with supersonic flow with subsonic leading edge are considered. To demonstrate the probabilistic approach, the flutter frequency, damping and forced response of a blade row representing a compressor geometry is considered. The analysis accounts for uncertainties in structural and aerodynamic design variables. The results are presented in the form of probabilistic density function (PDF) and sensitivity factors. For subsonic flow cascade, comparisons are also made with different probabilistic distributions, probabilistic methods, and Monte-Carlo simulation. The approach shows that the probabilistic approach provides a more realistic and systematic way to assess the effect of uncertainties in design variables on the aeroelastic instabilities and response.

Reddy, T. S. R.↗
NASA NTRS2013

Towards Probablistic Assessment of Hypobaric Decompression Sickness Treatment

INTRODUCTION: Pressure, oxygen (O2), and time are the pillars to effective treatment of decompression sickness (DCS). The NASA DCS Treatment Model links a decrease in computed bubble volume to the resolution of a symptom. The decrease in volume is realized in two stages: a) during the Boyle's Law compression and b) during subsequent dissolution of the gas phase by the O2 window. METHODS: The cumulative distribution of 154 symptoms that resolved during repressurization was described with a log-logistic density function of pressure difference (deltaP as psid) associated with symptom resolution and two other explanatory variables. The 154 symptoms originated from 119 cases of DCS during 969 exposures in 47 different altitude tests. RESULTS: The probability of symptom resolution [P(symptom resolution)] = 1 / (1+exp(- (ln(deltaP) - 1.682 + 1.089×AMB - 0.00395×SYMPTOM TIME) / 0.633)), where AMB is 1 when the subject ambulated as part of the altitude exposure or else 0 and SYMPTOM TIME is the elapsed time in min from start of the altitude exposure to recognition of a DCS symptom. The P(symptom resolution) was estimated from computed deltaP from the Tissue Bubble Dynamics Model based on the "effective" Boyle's Law change: P2 - P1 (deltaP, psid) = P1×V1/V2 - P1, where V1 is the computed volume of a spherical bubble in a unit volume of tissue at low pressure P1 and V2 is computed volume after a change to a higher pressure P2. V2 continues to decrease through time at P2, at a faster rate if 100% ground level O2 was breathed. The computed deltaP is the effective treatment pressure at any point in time as if the entire deltaP was just from Boyle's Law compression. DISCUSSION: Given the low probability of DCS during extravehicular activity and the prompt treatment of a symptom with options through the model it is likely that the symptom and gas phase will resolve with minimum resources and minimal impact on astronaut health, safety, and productivity.

Conkin, J.↗
NASA NTRS1987

A comparative study of some mathematical models of the mean wind structure and aerodynamic drag of plant canopies

A semianalytical method for describing the mean wind profile and shear stress within plant canopies and for estimating the roughness length and the displacement height is presented. This method incorporates density and vertical structure of the canopy and includes simple parameterizations of the roughness sublayer and shelter factor. Some of the wind profiles examined are consistent with first-order closure techniques while others are consistent with second-order closure techniques. Some profiles show a shearless region near the base of the canopy; however, none displays a secondary maximum there. Comparing several different analytical expressions for the canopy wind profile against observations suggests that one particular type of profile (an Airy function which is associated with the triangular foliage surface area density distribution) is superior to the others. Because of the numerical simplicity of the methods outlined, it is suggested that they may be profitably used in large-scale models of plant-atmosphere exchanges.

Massman, William↗
NASA NTRS2013

Understanding How Kurtosis Is Transferred from Input Acceleration to Stress Response and Its Influence on Fatigue Llife

High cycle fatigue of metals typically occurs through long term exposure to time varying loads which, although modest in amplitude, give rise to microscopic cracks that can ultimately propagate to failure. The fatigue life of a component is primarily dependent on the stress amplitude response at critical failure locations. For most vibration tests, it is common to assume a Gaussian distribution of both the input acceleration and stress response. In real life, however, it is common to experience non-Gaussian acceleration input, and this can cause the response to be non-Gaussian. Examples of non-Gaussian loads include road irregularities such as potholes in the automotive world or turbulent boundary layer pressure fluctuations for the aerospace sector or more generally wind, wave or high amplitude acoustic loads. The paper first reviews some of the methods used to generate non-Gaussian excitation signals with a given power spectral density and kurtosis. The kurtosis of the response is examined once the signal is passed through a linear time invariant system. Finally an algorithm is presented that determines the output kurtosis based upon the input kurtosis, the input power spectral density and the frequency response function of the system. The algorithm is validated using numerical simulations. Direct applications of these results include improved fatigue life estimations and a method to accelerate shaker tests by generating high kurtosis, non-Gaussian drive signals.

Kihm, Frederic↗
NASA NTRS1977

The analysis and interpretation of solar X-ray photographs

The paper examines methods of analyzing and interpreting solar X-ray photographs obtained by soft X-ray telescopes such as those used during the Skylab mission, to obtain information on the physical conditions that give rise to coronal X-ray emitting features revealed by inspection or microdensitometry of the X-ray photographs. First, the differential emission measure function is defined. By choosing models for this function and calculating the integral for the irradiance at the film for each filter, one calculates the effective temperature and emission measure. This method gives correct results only when applied to isothermal regions. For heterothermal regions, a modeling approach which predicts the run of temperature and density, and thus the form of the differential measure function, is used.

Underwood, J. H.↗
NASA NTRS1984

Radiative transfer theory for active remote sensing of a layer of nonspherical particles

The radiative transfer theory is applied to calculate the scattering by a layer of randomly positioned and oriented nonspherical particles. The scattering amplitude functions of each individual particle are calculated with Waterman's T matrix method, which utilizes vector spherical wave functions for expansion of incident, scattered, and surface fields. The orientation of the particles is described by a probability density function of the Eulerian angles of rotation. A rotation matrix is used to relate the T matrix of the principal frame to that of the natural frame of the particle. The extinction matrix and phase matrix of the radiative transfer equations are expressed in terms of the T matrix elements. The extinction matrix for nonspherical particles is generally nondiagonal. There are only two attenuation rates in a specified direction of propagation. The radiative transfer equations are solved by an iterative method to first order in albedo. Numerical results are illustrated as functions of incidence angle and frequency with applications to active remote sensing.

Tsang, L.↗
NASA NTRS2006

Learning User Preferences for Sets of Objects

Most work on preference learning has focused on pairwise preferences or rankings over individual items. In this paper, we present a method for learning preferences over sets of items. Our learning method takes as input a collection of positive examples--that is, one or more sets that have been identified by a user as desirable. Kernel density estimation is used to estimate the value function for individual items, and the desired set diversity is estimated from the average set diversity observed in the collection. Since this is a new learning problem, we introduce a new evaluation methodology and evaluate the learning method on two data collections: synthetic blocks-world data and a new real-world music data collection that we have gathered.

preferences↗
NASA NTRS1984

Use of shape-preserving interpolation methods in surface modeling

In many large-scale scientific computations, it is necessary to use surface models based on information provided at only a finite number of points (rather than determined everywhere via an analytic formula). As an example, an equation of state (EOS) table may provide values of pressure as a function of temperature and density for a particular material. These values, while known quite accurately, are typically known only on a rectangular (but generally quite nonuniform) mesh in (T,d)-space. Thus interpolation methods are necessary to completely determine the EOS surface. The most primitive EOS interpolation scheme is bilinear interpolation. This has the advantages of depending only on local information, so that changes in data remote from a mesh element have no effect on the surface over the element, and of preserving shape information, such as monotonicity. Most scientific calculations, however, require greater smoothness. Standard higher-order interpolation schemes, such as Coons patches or bicubic splines, while providing the requisite smoothness, tend to produce surfaces that are not physically reasonable. This means that the interpolant may have bumps or wiggles that are not supported by the data. The mathematical quantification of ideas such as physically reasonable and visually pleasing is examined.

Ftitsch, F. N.↗
NASA NTRS1977

Probing the solar wind with radio measurements of the second moment field

The dependence of spectral broadening on anisotropic electron-density irregularities and velocity fluctuations in the solar wind is investigated on the basis of results for a second moment mutual coherence function obtained by a parabolic-equation method. It is found that although anisotropic irregularities and wind-velocity fluctuations decrease the bandwidth of the spectrum relative to that for isotropic irregularities and no velocity fluctuations, the shape of the spectrum remains unchanged. Examination of the frequency dependence and shape of a computed angular spectrum indicates that the electron-density spectrum is consistent with a power law having a spectral index close to the value appropriate for the Kolmogorov spectrum. It is proposed that multiple-station observations of some total field be employed to measure the solar wind.

Woo, R.↗
NASA NTRS1974

Estimations of the smoothing operator response characteristics

The mean response characteristic of the graphical smoothing method is discussed. The method is illustrated by analysis of latitude observations at Washington from 1915.9 to 1941.0. Spectral density, frequency distribution, and distribution functions are also discussed.

Yatskiv, Y. S.↗
NASA NTRS1975

Tunnel injection into gate oxide traps

An experimental method is described for measuring the density of oxide traps in the gate oxide of an MOS transistor as a function of energy and position near the silicon interface. Measurements are obtained from different oxide growth processes and after Co-60 irradiation. The results are related to long-term drift of threshold voltage.

Maserjian, J.↗
NASA NTRS1992

Dielectronic recombination as a function of electric field strength

Dielectronic recombination (DR) is the dominant recombination mechanism at coronal temperatures and densities. We present a procedure for calculating DR rate coefficients as a function of electric field strength and apply this method to carbon ions. We focus on the competing effects of enhancement by plasma microfields and rate decrease through collisional excitation and ionization. We find that, in the case of C(3+), a significant rate enhancement results, leading to a reinterpretation of C IV emission-line intensities in the sun and late-type stars. We further consider how macroscopic electric fields, in particular motional electric fields, can affect DR rate coefficients, demonstrating dramatic rate increases for a number of the carbon ions.

Reisenfeld, Daniel B.↗
NASA NTRS1999

Phase Diagram of the Two-Chain Hubbard Model

We have calculated the charge gap and spin gap for the two-chain Hubbard model as a function of the on-site Coulomb interaction and the interchain hopping amplitude. We used the density matrix renormalization group method and developed a method to calculate separately the gaps numerically for the symmetric and antisymmetric modes with respect to the exchange of the chain indices. We have found very different behaviors for the weak and strong interaction cases. Our calculated phase diagram is compared to the one obtained by Balents and Fisher using the weak coupling renormalization group technique.

Park, Youngho↗
NASA NTRS1994

The Molecular Structure of cis-FONO

The molecular structure of cis-FONO has been determined with the CCSD(T) correlation method using an spdf quality basis set. In agreement with previous coupled-cluster calculations but in disagreement with density functional theory, cis-FONO is found to exhibit normal bond distances. The quadratic and cubic force fields of cis-FONO have also been determined in order to evaluate the effect of vibrational averaging on the molecular geometry. Vibrational averaging is found to increase bond distances, as expected, but it does not affect the qualitative nature of the bonding. The CCSD(T)/spdf harmonic frequencies of cis-FONO support our previous assertion that a band observed at 1200 /cm is a combination band (upsilon(sub 3) + upsilon(sub 4)), and not a fundamental.

Lee, Timothy J.↗
NASA NTRS1997

Formation of Carbon Nanotube Based Gears: Quantum Chemistry and Molecular Mechanics Study of the Electrophilic Addition of o-Benzyne to Fullerenes, Graphene, and Nanotubes

Considerable progress has been made in recent years in chemical functionalization of fullerene molecules. In some cases, the predominant reaction products are different from those obtained (using the same reactants) from polycyclic aromatic hydrocarbons (PAHs). One such example is the cycloaddition of o-benzyne to C60. It is well established that benzyne adds across one of the rings in naphthalene, anthracene and other PAHs forming the [2+4] cycloaddition product (benzobicyclo[2.2.2.]-octatriene with naphthalene and triptycene with anthracene). However, Hoke et al demonstrated that the only reaction path for o-benzyne with C60 leads to the [2+2] cycloaddition product in which benzyne adds across one of the interpentagonal bonds (forming a cyclobutene ring in the process). Either reaction product results in a loss of aromaticity and distortion of the PAH or fullerene substrate, and in a loss of strain in the benzyne. It is not clear, however, why different products are preferred in these cases. In the current paper, we consider the stability of benzyne-nanotube adducts and the ability of Brenner's potential energy model to describe the structure and stability of these adducts. The Brenner potential has been widely used for describing diamondoid and graphitic carbon. Recently it has also been used for molecular mechanics and molecular dynamics simulations of fullerenes and nanotubes. However, it has not been tested for the case of functionalized fullerenes (especially with highly strained geometries). We use the Brenner potential for our companion nanogear simulations and believe that it should be calibrated to insure that those simulations are physically reasonable. In the present work, Density Functional theory (DFT) calculations are used to determine the preferred geometric structures and energetics for this calibration. The DFT method is a kind of ab initio quantum chemistry method for determining the electronic structure of molecules. For a given basis set expansion, it is comparable in accuracy to the MP2 method (better than Hartree Fock, but less accurate than more extensive electron correlation methods such as MP4 or CCSD). However, for systems with large numbers of basis functions it more efficient than any other methods that include electron correlation effects. In this presentation we show the results of DFT calculations for the reaction of benzyne with naphthalene, C60, and nanotube models. We compare energies for [2+2] and [2+4] cycloaddition products. The preferred products for the naphthalene and C60 reactions have been determined by experiment and, thus, these cases serve as a validation of our quantum chemical approach. We also compare the DFT and Brenner potential results. Finally we can predict the likelihood of reaction between benzyne and nanotubes.

Jaffe, Richard↗
NASA NTRS2011

Ab initio Computations of the Electronic, Mechanical, and Thermal Properties of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2

Refractory materials such as metallic borides, often considered as ultra high temperature ceramics (UHTC), are characterized by high melting point, high hardness, and good chemical inertness. These materials have many applications which require high temperature materials that can operate with no or limited oxidation. Ab initio, first principles methods are the most accurate modeling approaches available and represent a parameter free description of the material based on the quantum mechanical equations. Using these methods, many of the intrinsic properties of these material can be obtained. We performed ab initio calculations based on density functional theory for the UHTC materials ZrB2 and HfB2. Computational results are presented for structural information (lattice constants, bond lengths, etc), electronic structure (bonding motifs, densities of states, band structure, etc), thermal quantities (phonon spectra, phonon densities of states, specific heat), as well as information about point defects such as vacancy and antisite formation energies.

Lawson, John W.↗
NASA NTRS1971

Techniques of orbital decay and long-term ephemeris prediction for satellites in earth orbit

In the special perturbation method, Cowell and variation-of-parameters formulations of the motion equations are implemented and numerically integrated. Variations in the orbital elements due to drag are computed using the 1970 Jacchia atmospheric density model, which includes the effects of semiannual variations, diurnal bulge, solar activity, and geomagnetic activity. In the general perturbation method, two-variable asymptotic series and automated manipulation capabilities are used to obtain analytical solutions to the variation-of-parameters equations. Solutions are obtained considering the effect of oblateness only and the combined effects of oblateness and drag. These solutions are then numerically evaluated by means of a FORTRAN program in which an updating scheme is used to maintain accurate epoch values of the elements. The atmospheric density function is approximated by a Fourier series in true anomaly, and the 1970 Jacchia model is used to periodically update the Fourier coefficients. The accuracy of both methods is demonstrated by comparing computed orbital elements to actual elements over time spans of up to 8 days for the special perturbation method and up to 356 days for the general perturbation method.

Barry, B. F.↗