Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional approximations”

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 91 records · Page 5

Numerical models for Halley dust environments

The algebra, parameter values, and results of a reference model for the dust environment of Comet P/Halley are presented. The model employs simplified treatments of the icy conglomerate nucleus and of the hydrodynamic outflow of gas and dust. The model includes a gas production rate derived from observations of Halley's visual magnitude in 1910, gas and dust flux nonuniformly distributed over the nuclear surface, a dust distribution function consistent with large particle inferences from comet antitails, dust density variable with mass, and an approximation to the distribution of particles near the dust envelope. Many parameter values are derived by analogy with the results for other comets. The features are incorporated in a computer program which calculates the concentration, flux, and fluence of dust particles as functions of time, position and particle mass.

Divine, N.↗

A mathematical analysis of the theory of interplanetary scintillation in the weak scattering approximation

A simplified analytical technique is presented for modeling the interplanetary scintillation of radio sources of finite angular size with a power-law electron-density-fluctuation power spectrum. The simplification results from the representation of the scintillation spectrum in confluent hypergeometric functions. The approximations presented allow fast numerical evaluation of a spectrum for a weakly scattering but extended medium with less than 10% error over the entire spectrum. Parameters describing anisotropic electron irregularities as well as anisotropic source structure are included, and the dependence of the spectrum normalization on the scales of the medium is derived explicitly. The parametric description of the domains of convergence of the approximate expansions also provides a simple conceptualization of the relative contributions of the scattered radiation along the line of sight to the observed spectrum. This is particularly useful for sources of finite angular size. This technique is applied to previously published observations.

Mitchell, D. G.↗

Structural Reliability Analysis and Optimization: Use of Approximations

This report is intended for the demonstration of function approximation concepts and their applicability in reliability analysis and design. Particularly, approximations in the calculation of the safety index, failure probability and structural optimization (modification of design variables) are developed. With this scope in mind, extensive details on probability theory are avoided. Definitions relevant to the stated objectives have been taken from standard text books. The idea of function approximations is to minimize the repetitive use of computationally intensive calculations by replacing them with simpler closed-form equations, which could be nonlinear. Typically, the approximations provide good accuracy around the points where they are constructed, and they need to be periodically updated to extend their utility. There are approximations in calculating the failure probability of a limit state function. The first one, which is most commonly discussed, is how the limit state is approximated at the design point. Most of the time this could be a first-order Taylor series expansion, also known as the First Order Reliability Method (FORM), or a second-order Taylor series expansion (paraboloid), also known as the Second Order Reliability Method (SORM). From the computational procedure point of view, this step comes after the design point identification; however, the order of approximation for the probability of failure calculation is discussed first, and it is denoted by either FORM or SORM. The other approximation of interest is how the design point, or the most probable failure point (MPP), is identified. For iteratively finding this point, again the limit state is approximated. The accuracy and efficiency of the approximations make the search process quite practical for analysis intensive approaches such as the finite element methods; therefore, the crux of this research is to develop excellent approximations for MPP identification and also different approximations including the higher-order reliability methods (HORM) for representing the failure surface. This report is divided into several parts to emphasize different segments of the structural reliability analysis and design. Broadly, it consists of mathematical foundations, methods and applications. Chapter I discusses the fundamental definitions of the probability theory, which are mostly available in standard text books. Probability density function descriptions relevant to this work are addressed. In Chapter 2, the concept and utility of function approximation are discussed for a general application in engineering analysis. Various forms of function representations and the latest developments in nonlinear adaptive approximations are presented with comparison studies. Research work accomplished in reliability analysis is presented in Chapter 3. First, the definition of safety index and most probable point of failure are introduced. Efficient ways of computing the safety index with a fewer number of iterations is emphasized. In chapter 4, the probability of failure prediction is presented using first-order, second-order and higher-order methods. System reliability methods are discussed in chapter 5. Chapter 6 presents optimization techniques for the modification and redistribution of structural sizes for improving the structural reliability. The report also contains several appendices on probability parameters.

Grandhi, Ramana V.↗

Solar albedo and spectral reflectance for Apollo 15 and 16 lunar fines

The spectral directional reflectance of Apollo 15 and 16 fines was obtained for bulk densities of approximately 1000 and 1600 kg/cu m. The solar albedo as a function of illumination was calculated from these results. Comparison of solar albedos shows that Apollo 11, 12, and 15 soils fall into one group and the Apollo 14 and 16 soil results into a second group corresponding to mare or lunar highland materials.

Birkebak, R. C.↗

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.↗

AutoClass: A Bayesian Approach to Classification

We describe a Bayesian approach to the untutored discovery of classes in a set of cases, sometimes called finite mixture separation or clustering. The main difference between clustering and our approach is that we search for the "best" set of class descriptions rather than grouping the cases themselves. We describe our classes in terms of a probability distribution or density function, and the locally maximal posterior probability valued function parameters. We rate our classifications with an approximate joint probability of the data and functional form, marginalizing over the parameters. Approximation is necessitated by the computational complexity of the joint probability. Thus, we marginalize w.r.t. local maxima in the parameter space. We discuss the rationale behind our approach to classification. We give the mathematical development for the basic mixture model and describe the approximations needed for computational tractability. We instantiate the basic model with the discrete Dirichlet distribution and multivariant Gaussian density likelihoods. Then we show some results for both constructed and actual data.

Stutz, John↗

Quantiles, parametric-select density estimation, and bi-information parameter estimators

A quantile-based approach to statistical analysis and probability modeling of data is presented which formulates statistical inference problems as functional inference problems in which the parameters to be estimated are density functions. Density estimators can be non-parametric (computed independently of model identified) or parametric-select (approximated by finite parametric models that can provide standard models whose fit can be tested). Exponential models and autoregressive models are approximating densities which can be justified as maximum entropy for respectively the entropy of a probability density and the entropy of a quantile density. Applications of these ideas are outlined to the problems of modeling: (1) univariate data; (2) bivariate data and tests for independence; and (3) two samples and likelihood ratios. It is proposed that bi-information estimation of a density function can be developed by analogy to the problem of identification of regression models.

Parzen, E.↗

Determination of Particle Sources for a Geotail Distribution Function Observed on May 23, 1995

On May 23, 1995, the Comprehensive Plasma Instrumentation (CPI) onboard the Geotail spacecraft observed a complex and structured ion distribution function near the magnetotail midplane at x approximately -10 R(sub E). On the same day, the Wind spacecraft observed a very high density (approximately 40/cubic cm) solar wind and an interplanetary magnetic field (IMF) that was predominantly northward but had several southward turnings. We have inferred the sources of the ions in this distribution function by following approximately 90,000 ion trajectories backward in time using time-dependent electric and magnetic fields obtained from a global MHD (magnetohydrodynamic) simulation. Wind data were used as input for the MHD model. We found that three sources contributed to this distribution: the ionosphere, the plasma mantle which had near-Earth and distant tail components, and the low latitude boundary layer (LLBL). Moreover, distinct structures in the low energy part of the distribution function were found to be associated with individual sources. Structures near 0 deg pitch angle were made up of either ionospheric or plasma mantle ions, while structures near 90 deg pitch angle were dominated by ions from the LLBL source. Particles that underwent nonadiabatic acceleration were numerous in the higher energy part of the ion distribution function, whereas ionospheric and LLBL ions were mostly adiabatic. A large proportion of the near-Earth mantle ions underwent adiabatic acceleration, while most of the distant mantle ions experienced nonadiabatic acceleration.

Ashour-Abdalla, M.↗

Thermal characteristics of Apollo 16 lunar fines

The vacuum thermal conductivity of the Apollo 16 fines is presented as a function of temperature for the approximate range of diurnal temperatures on the moon. The density used is 1500 kg/cu m, which is approximately the density of the core tube samples returned from the landing site. The thermal diffusivity was calculated from the conductivity, density and specific heat and is also presented. These properties along with the reflectance and emittance are used in a calculation of temperature in the lunar surface layer.

Cremers, C. J.↗

Thermal conductivity of Apollo 16 lunar fines

The vacuum thermal conductivity of the Apollo 16 fines is presented as a function of temperature for the approximate range of diurnal temperatures on the moon. The density used is 1500 kg/cu m, which is approximately the density of the core-tube samples returned from the landing site. The thermal diffusivity was calculated from the conductivity, density, and specific heat and is also presented.

Cremers, C. J.↗

Critical currents in sputtered PbMo6S8

Critical currents in sputtered Chevrel-phase lead molybdenum sulfide have been measured at several temperatures as a function of applied magnetic field up to 19 T. A critical current density of approximately 50 million A/sq m was obtained at 15 T and 1.8 K, and the effective upper critical field was estimated to be 30 T. The pinning force at low fields was found to be dependent on the amount of free lead.

Alterovitz, S. A.↗

New explicit equations for the accurate calculation of the growth and evaporation of hydrometeors by the diffusion of water vapor

The traditional explicit growth equation has been widely used to calculate the growth and evaporation of hydrometeors by the diffusion of water vapor. This paper reexamines the assumptions underlying the traditional equation and shows that large errors (10-30 percent in some cases) result if it is used carelessly. More accurate explicit equations are derived by approximating the saturation vapor-density difference as a quadratic rather than a linear function of the temperature difference between the particle and ambient air. These new equations, which reduce the error to less than a few percent, merit inclusion in a broad range of atmospheric models.

Srivastava, R. C.↗

Molecular Static Third-Order Polarizabilities of Carbon-Cage Fullerene and Their Correlation with Three Geometric Properties: Symmetry, Aromaticity, and Size

The static third-order polarizabilities (gamma) of C60, C70, five isomers of C78 and two isomers of C84 were analyzed in terms of three properties, from a geometric point of view: symmetry, aromaticity and size. The polarizability values were based on the finite field approximation using a semiempirical Hamiltonian (AM1) and applied to molecular structures obtained from density functional theory calculations. Symmetry was characterized by the molecular group order. The selection of 6-member rings as aromatic was determined from an analysis of bond lengths. Maximum interatomic distance and surface area were the parameters considered with respect to size. Based on triple linear regression analysis, it was found that the static linear polarizability (alpha) and gamma in these molecules respond differently to geometrical properties: alpha depends almost exclusively on surface area while gamma is affected by a combination of number of aromatic rings, length and group order, in decreasing importance. In the case of alpha, valence electron contributions provide the same information as all-electron estimates. For gamma, the best correlation coefficients are obtained when all-electron estimates are used and when the dependent parameter is ln(gamma) instead of gamma.

Moore, C. E.↗

Constituent Effects on the Stress-Strain Behavior of Woven Melt-Infiltrated SiC Composites

The stress-strain behavior of 2D woven SiC fiber reinforced, melt-infiltrated SiC matrix composites with BN interphases were studied for composites fabricated with different fiber tow ends per unit length, different composite thickness, and different numbers of plies. In general, the stress-strain behavior, i.e., the 'knee' in the curve and the final slope of the stress-strain curve, was controlled by the volume fraction of fibers. Some of the composites exhibited debonding and sliding in between the interphase and the matrix rather than the more common debonding and sliding interface between the fiber and the interphase. Composites that exhibited this 'outside debonding' interface, in general, had lower elastic moduli and higher ultimate strains as well as longer pull-out lengths compared to the 'inside debonding' interface composites. Stress-strain curves were modeled where matrix crack formation as a function of stress was approximated from the acoustic emission activity and the measured crack density from the failed specimens. Interfacial shear strength measurements from individual fiber push-in tests were in good agreement with the interfacial shear strength values used to model the stress-strain curves.

Morscher, Gregory N.↗

Heats of Segregation of BCC Metals Using Ab Initio and Quantum Approximate Methods

Many multicomponent alloys exhibit surface segregation, in which the composition at or near a surface may be substantially different from that of the bulk. A number of phenomenological explanations for this tendency have been suggested, involving, among other things, differences among the components' surface energies, molar volumes, and heats of solution. From a theoretical standpoint, the complexity of the problem has precluded a simple, unified explanation, thus preventing the development of computational tools that would enable the identification of the driving mechanisms for segregation. In that context, we investigate the problem of surface segregation in a variety of bcc metal alloys by computing dilute-limit heats of segregation using both the quantum-approximate energy method of Bozzolo, Ferrante and Smith (BFS), and all-electron density functional theory. In addition, the composition dependence of the heats of segregation is investigated using a BFS-based Monte Carlo procedure, and, for selected cases of interest, density functional calculations. Results are discussed in the context of a simple picture that describes segregation behavior as the result of a competition between size mismatch and alloying effects

Good, Brian↗

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗

Analytical methods for performance evaluation of nonlinear filters.

In the investigation, the filtering problem is considered in the continuous time domain. The postulated simple suboptimal nonlinear filter structure closely parallels the structure of the Kalman-Bucy optimal linear filter algorithm. Two filter performance evaluation methods are developed based on the Kolmogorov equations for the transition density of Markov processes. The expansions in the approximations for the nonlinear system and observation functions are in effect carried out up to second-order terms in both methods. The description of the filter's performance is sought in terms of second-order statistics in both methods.

Bejczy, A. K.↗