Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

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 217 records · Page 12

Modeling and Visualizing Uncertainty in Continuous Variables Predicted using Remotely Sensed Data

The use of remotely sensed images to map continuous biophysical variables, such as those related to terrestrial vegetation amount, sea surface temperature, and many other targets of NASA s Earth Observing System (EOS), includes variable, parametric, positional, spatial support and structural sources of uncertainty. A complete description of uncertainty will lead to a probability distribution at each location, allowing the exploration of the spatial dimension of uncertainty, that is, where the field is not well quantified. To achieve this purpose, convenient visualization tools are required. We have produced such a tool, called PDFVis, that facilitates the display of probability density functions (pdfs) on a per-grid-cell basis. The density estimate from Monte-Carlo generated realizations is interactively displayed as well as parametric and non-parametric summaries of the pdf field (such as mean, median, quartiles, standard deviation, number of modes, and locations of modes). Shaded surface renderings of pdfs along a transect can also be projected onto a plane. This tool will become more useful as richer descriptions of spatial uncertainty become available.

Dungan, Jennifer L.↗

Statistical Comparison 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 paper examines the cumulative distribution function (CDF) 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 Mathematics is then used to calculate the combined value corresponding to any desired percentile along with a curve fit to this value. Another Excel macro is used to calculate the combination using a Monte Carlo simulation. 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. Also, 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 lower the design loading substantially without losing any of the identified structural reliability.

Brown, Andrew M.↗

Computation of Domain-Averaged Shortwave Irradiance by a One-Dimensional Algorithm Incorporating Correlations between Optical Thickness and Direct Incident Radiation

A one-dimensional radiative transfer algorithm that accounts for correlations between the optical thickness and the incident direct solar radiation is developed to compute the domain-averaged shortwave irradiance profile. It divides the direct irradiance into four components and treats the direct irradiance in two separate, clear and cloudy columns to account for the fact that clouds attenuate the direct irradiance more than clear-sky. The horizontal inhomogeneity of clouds in the cloudy column is treated by the gamma weighted two-stream approximation, which assumes that the optical thickness of clouds follows a gamma distribution. The algorithm inputs the cloud fraction, cumulative cloud fraction as a function of height, and a parameter expressing the shape of the probability density function of the cloud optical thickness distribution in addition to inputs required for a two-stream radiative transfer model. These cloud property inputs can be obtained using ground- and satellite-based instruments. Therefore, the algorithm can treat realistic cloud overlap features and horizontal inhomogeneity of clouds in a framework of one- dimensional radiative transfer. Heating rates computed by the algorithm using cloud fields generated by cloud resolving models agree with those computed with a Monte Carlo model. If optical properties in computational layers that divide a vertically extensive cloud are correlated, the irradiance profile computed by the algorithm further improves.

Kato, S.↗

Comparison of Ares I-X Wind-Tunnel Derived Buffet Environment with Flight Data

The Ares I-X Flight Test Vehicle (FTV), launched in October 2009, carried with it over 243 buffet verification pressure sensors and was one of the most heavily instrumented launch vehicle flight tests. This flight test represented a unique opportunity for NASA and its partners to compare the wind-tunnel derived buffet environment with that measured during the flight of Ares I-X. It is necessary to define the launch vehicle buffet loads to ensure that structural components and vehicle subsystems possess adequate strength, stress, and fatigue margins when the vehicle structural dynamic response to buffet forcing functions are considered. Ares I-X buffet forcing functions were obtained via wind-tunnel testing of a rigid buffet model (RBM) instrumented with hundreds of unsteady pressure transducers designed to measure the buffet environment across the desired frequency range. This paper discusses the comparison of RBM and FTV buffet environments, including fluctuating pressure coefficient and normalized sectional buffet forcing function root-mean-square magnitudes, frequency content of power-spectral density functions, and force magnitudes of an alternating flow phenomena. Comparison of wind-tunnel model and flight test vehicle buffet environments show very good agreement with root-mean-square magnitudes of buffet forcing functions at the majority of vehicle stations. Spectra proved a challenge to compare because of different wind-tunnel and flight test conditions and data acquisition rates. However, meaningful and promising comparisons of buffet spectra are presented. Lastly, the buffet loads resulting from the transition of subsonic separated flow to supersonic attached flow were significantly over-predicted by wind-tunnel results.

Piatak, David J.↗

Binding of an adatom to a simple metal surface

The density functional formalism of Hohenberg and Kohn is used to investigate the energies, charge densities and forces which hold an adatom on the surface of a simple metal. The valence wavefunction of the adatom is fitted to the Herman-Skillman solutions at large distance and is simplified somewhat in the core region. The field of the ion is represented by the Ashcroft pseudopotential. For the metal the jellium model is used. Detailed calculations are carried out for a sodium adatom on a sodium surface. Simply juxtaposing adatom and surface gives a binding energy of about 1/3 eV. This value is approximately twice the surface energy per atom in the close-packed plane. Charge redistributions as determined variationally increase the binding energy by about 10%. The equilibrium distance for the adatom turns out to be 1.66 A from the surface, as compared with 1.52 A, the observed value for one-half the distance between the close-packed planes.

Huntington, H. B.↗

Picturing Data With Uncertainty

NASA is in the business of creating maps for scientific purposes to represent important biophysical or geophysical quantities over space and time. For example, maps of surface temperature over the globe tell scientists where and when the Earth is heating up; regional maps of the greenness of vegetation tell scientists where and when plants are photosynthesizing. There is always uncertainty associated with each value in any such map due to various factors. When uncertainty is fully modeled, instead of a single value at each map location, there is a distribution expressing a set of possible outcomes at each location. We consider such distribution data as multi-valued data since it consists of a collection of values about a single variable. Thus, a multi-valued data represents both the map and its uncertainty. We have been working on ways to visualize spatial multi-valued data sets effectively for fields with regularly spaced units or grid cells such as those in NASA's Earth science applications. A new way to display distributions at multiple grid locations is to project the distributions from an individual row, column or other user-selectable straight transect from the 2D domain. First at each grid cell in a given slice (row, column or transect), we compute a smooth density estimate from the underlying data. Such a density estimate for the probability density function (PDF) is generally more useful than a histogram, which is a classic density estimate. Then, the collection of PDFs along a given slice are presented vertically above the slice and form a wall. To minimize occlusion of intersecting slices, the corresponding walls are positioned at the far edges of the boundary. The PDF wall depicts the shapes of the distributions very dearly since peaks represent the modes (or bumps) in the PDFs. We've defined roughness as the number of peaks in the distribution. Roughness is another useful summary information for multimodal distributions. The uncertainty of the multi-valued data can also be interpreted by the number of peaks and the widths of the peaks as shown by the PDF walls.

Kao, David↗

A dual frequency radar for ocean roughness sampling

A microwave technique for determining the roughness (rms wave height) of a randomly varying air-water interface has been developed theoretically, verified with laboratory wavetank studies and is currently being implemented for ocean surface measurements in a series of applications flight experiments. These aircraft observations will be near the Chesapeake Light Tower and will include a range of altitudes and sea conditions. The measurement concept involves cross-correlating the envelope fluctuations on two received carriers that are monochromatic when transmitted in a normal direction to the interface and are observed in a backscatter direction after reflection by the large number of randomly distributed specular points on the surface. The measured correlation coefficient (normalized covariance) as a function of carrier frequency separation will depend on the surface roughness, specifically the probability density function of the specular point height. Details of the flight system are discussed, and preliminary flight results are presented.

Weissman, D. E.↗

Weibull crack density coefficient for polydimensional stress states

A structural ceramic analysis and reliability evaluation code has recently been developed encompassing volume and surface flaw induced fracture, modeled by the two-parameter Weibull probability density function. A segment of the software involves computing the Weibull polydimensional stress state crack density coefficient from uniaxial stress experimental fracture data. The relationship of the polydimensional stress coefficient to the uniaxial stress coefficient is derived for a shear-insensitive material with a random surface flaw population.

Gross, Bernard↗

Heat transfer to and from vegetated surfaces - An analytical method for the bulk exchange coefficients

The semianalytical model outlined in a previous study (Massman, 1987) to describe momentum exchange between the atmosphere and vegetated surfaces is extended to include the exchange of heat. The methods employed are based on one-dimensional turbulent diffusivities, and use analytical solutions to the steady-state diffusion equation. The model is used to assess the influence that the canopy foliage structure and density, the wind profile structure within the canopy, and the shelter factor can have upon the inverse surface Stanton number (kB exp -1), as well as to explore the consequences of introducing a scalar displacement height which can be different from the momentum displacement height. In general, the triangular foliage area density function gives results which agree more closely with observations than that for constant foliage area density. The intended application of this work is for parameterizing the bulk aerodynamic resistances for heat and momentum exchange for use within large-scale models of plant-atmosphere exchanges.

Massman, William J.↗

A new adaptive control approach for aerospace vehicles with parameter uncertainties

A new stochastic adaptive control structure is developed for the problem of combined parameter estimation and control of aerospace vehicles with changing parameters. Parameter uncertainties are modeled as first-order Gauss-Markov processes, and are introduced to the system dynamics through a small parameter. It is assumed that an accurate inertial measurement unit gives perfect measurements of the state variables. Since the stochastic system is assumed to be Gauss-Markov, the density function of the parameters given these measurements is conditionally Gaussian. Based on this conditionally Gaussian density, the problem of minimizing a quadratic cost over an infinite time horizon can be set up within the framework of stochastic optimal control theory. The optimal feedback control law is derived from a straightforward expansion of the Hamilton-Jacobi-Bellman equation, based on the LQG solution. The resulting nonlinear controller is applied to the pitch axis control of a space platform with uncertain moments of inertia and is shown to produce marked improvement over a fixed controller.

Hahn, Yungsun↗

Studies on spectral analysis of randomly sampled signals: Application to laser velocimetry data

Spectral analysis is very useful in determining the frequency characteristics of many turbulent flows, for example, vortex flows, tail buffeting, and other pulsating flows. It is also used for obtaining turbulence spectra from which the time and length scales associated with the turbulence structure can be estimated. These estimates, in turn, can be helpful for validation of theoretical/numerical flow turbulence models. Laser velocimetry (LV) is being extensively used in the experimental investigation of different types of flows, because of its inherent advantages; nonintrusive probing, high frequency response, no calibration requirements, etc. Typically, the output of an individual realization laser velocimeter is a set of randomly sampled velocity data. Spectral analysis of such data requires special techniques to obtain reliable estimates of correlation and power spectral density functions that describe the flow characteristics. FORTRAN codes for obtaining the autocorrelation and power spectral density estimates using the correlation-based slotting technique were developed. Extensive studies have been conducted on simulated first-order spectrum and sine signals to improve the spectral estimates. A first-order spectrum was chosen because it represents the characteristics of a typical one-dimensional turbulence spectrum. Digital prefiltering techniques, to improve the spectral estimates from randomly sampled data were applied. Studies show that the spectral estimates can be increased up to about five times the mean sampling rate.

Sree, David↗

Symmetry Breaking and the B3LYP Functional

The infrared spectra of six molecules, each of which contains a five-membered ring, and their cations are determined using density functional theory (DFT); both the B3LYP and BP86 functionals are used. The computed results are compared with the experimental spectra. For the neutral molecules, both methods are in good agreement with experiment. Even the Hartree-Fock (HF) approach is qualitatively correct for the neutrals. For the cations, the HF approach fails, as found for other organic ring systems. The B3LYP and BP86 approaches are in good mutual agreement for five of the six cation spectra, and in good agreement with experiment for four of the five cations where the experimental spectra are available. It is only for the fluoranthene cation, where the BP86 and B3LYP functionals yield different results; the BP86 yields the expected C2v symmetry, while the B3LYP approach breaks symmetry. The experimental spectra supports the BP86 spectra over the B3LYP, but the quality of the experimental spectra does not allow a critical evaluation of the accuracy of the BP86 approach for this difficult system.

Bauschlicher, Charles W., Jr.↗

Remote optical measurement of wind speed spectra

The paper presents some new analytical and experimental results on the cross-correlation function and power spectrum for the flow speed by the passive optical crossed-beam correlation technique. Cross-correlation functions of crossed-beam measurements for both intersecting beams and separated beams have been reviewed for the case of varying wind speeds. The corresponding spectral density functions were derived and expressed in terms of the probability density of wind speed variations and the correlation of atmospheric turbulent eddies. The associated spectral density for a special approximate model was derived by assuming a Gaussian distribution of transit times. The theoretical predictions were then compared with actual field test data, and a reasonably good agreement in the valid frequency range was observed.

Su, M.-Y.↗

An investigation into the probabilistic combination of quasi-static and random accelerations

The development of design load factors for aerospace and aircraft components and experiment support structures, which are subject to a simultaneous vehicle dynamic vibration (quasi-static) and acoustically generated random vibration, require the selection of a combination methodology. Typically, the procedure is to define the quasi-static and the random generated response separately, and arithmetically add or root sum square to get combined accelerations. Since the combination of a probabilistic and a deterministic function yield a probabilistic function, a viable alternate approach would be to determine the characteristics of the combined acceleration probability density function and select an appropriate percentile level for the combined acceleration. The following paper develops this mechanism and provides graphical data to select combined accelerations for most popular percentile levels.

Schock, R. W.↗

Imperfection surveys on a 10-ft-diameter shell structure

The results of an extensive imperfection survey on a 10-ft-diameter integrally stiffened cylindrical shell are presented. The shape of the measured initial imperfections is clearly influenced by details of the shell construction. The modal components of the measured imperfection surface as a function of the circumferential and of the axial wave numbers are calculated. The discrete axial power spectral density functions and the corresponding root-mean-square values of the imperfections are also determined for given circumferential wave numbers. Using the Fourier coefficients of the measured initial imperfections, buckling loads are calculated by solving the nonlinear Donnell-type imperfect shell equations iteratively. The calculated lowest buckling load compares favorably with the values usually recommended for similar shell structures.

Arbocz, J.↗