Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy 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 109 records · Page 6

Lower solar chromosphere-corona transition region. III - Implications of the observed quiet-sun emission measure including wave pressure effects

The observed form of the emission measure (EM) is used as a function of temperature to infer the wave energy flux density and pressure throughout the lower transition region (TR). This procedure eliminates the need for specifying how the wave energy flux density is damped and addresses the question of whether there is any form of the mechanical heating associated with the degradation of an upward traveling wave energy flux density which is consistent with the observed EM and other observational constraints for the quiet sun. It is found that the observed form of the EM curve is incompatible with waves traveling vertically at the sound speed, regardless of any filling factor arguments. The same conclusion also applies to waves traveling at the Alfven speed, unless it is assumed that the emission in lower TR lines originates solely from small, spatially unresolved regions of large magnetic field strength (100 G), which cover a small fraction (filling factors of 1 percent) of the solar surface.

Woods, D. Tod↗

Estimating Resolution Lengths of Hybrid Turbulence Models

A two-stage procedure has been devised for estimating the spatial resolution achievable in the simulation of a given flow on a given computational grid by a computational fluid dynamics (CFD) code that incorporates a hybrid model of turbulence. The hybrid models to which this procedure is especially relevant are those of the Reynolds-averaged Navier-Stokes (RANS) and the partial-averaged Navier-Stokes (PANS) approaches. This procedure represents the first step toward adding variable-resolution turbulence-modeling capabilities to CFD codes as part of a continuing effort to increase the accuracy and robustness of CFD simulations of unsteady flows. Some background information is prerequisite to a meaningful summary of the procedure. Among experts in CFD, it is well known that combination of the Reynolds-averaged Navier-Stokes (RANS) approach and eddy-viscosity turbulence models offers limited capability for simulating unsteady and complex flows. The RANS approach includes an assumption that most of the energy in a given flow is modeled through turbulence-transport equations and is resolved in a computational grid used to simulate the flow. RANS also overpredicts eddy viscosity, thereby yielding excessive damping of unsteady motion. The eddy viscosity attains an unphysically large value because of unresolved scales, and suppresses most temporal and spatial fluctuations in the resolved flow field. One approach used to overcome this deficiency is to provide a mechanism for the RANS equations to resolve motion only on the largest scales and to use a hybrid model to represent effects at smaller scales. The RANS approach involves the use of a standard two-equation turbulence model in which the effect of turbulence is summarized by a viscosity that is a function of (1) the time-averaged kinetic- energy density (k) associated with the local fluctuating (turbulent) component of flow and (2) the time-averaged rate of dissipation of the turbulent-kinetic- energy density ( ). In PANS, which was developed to overcome the grid dependency associated with other hybrid turbulence models (including that of RANS), the standard two-equation turbulence model is replaced by another two-equation model in which one solves for the previously unresolved k and , which are now allowed to vary in space and time.

Abdol-Hamid, Khaled S.↗

The Role of Higher-Order Modes on the Electromagnetic Whistler-Cyclotron Wave Fluctuations of Thermal and Non-Thermal Plasmas

Two fundamental challenging problems of laboratory and astrophysical plasmas are the understanding of the relaxation of a collisionless plasmas with nearly isotropic velocity distribution functions and the resultant state of nearly equipartition energy density with electromagnetic plasma turbulence. Here, we present the results of a study which shows the role that higher-order-modes play in limiting the electromagnetic whistler-like fluctuations in a thermal and non-thermal plasma. Our main results show that for a thermal plasma the magnetic fluctuations are confined by regions that are bounded by the least-damped higher order modes. We further show that the zone where the whistler-cyclotron normal modes merges the electromagnetic fluctuations shifts to longer wavelengths as the beta(sub e) increases. This merging zone has been interpreted as the beginning of the region where the whistler-cyclotron waves losses their identity and become heavily damped while merging with the fluctuations. Our results further indicate that in the case of nonthermal plasmas, the higher-order modes do not confine the fluctuations due to the effective higher-temperature effects and the excess of suprathermal plasma particles. The analysis presented here considers the second-order theory of fluctuations and the dispersion relation of weakly transverse fluctuations, with wave vectors parallel to the uniform background magnetic field, in a finite temperature isotropic bi-Maxwellian and Tsallis-kappa-like magnetized electron-proton plasma. Our results indicate that the spontaneously emitted electromagnetic fluctuations are in fact enhanced over these quasi modes suggesting that such modes play an important role in the emission and absorption of electromagnetic fluctuations in thermal or quasi-thermal plasmas.

Vinas, Adolfo F.↗

C2 Fragmentation Energy of C60 Revisited: Theory Disagrees with Most Experiments

Following our earlier work on the subject, we have carried out density functional theory (DFT) and second-order Moller-Plesset perturbation theory (MP2) calculations of the dissociation energy of the reaction C60 yields C58 + C2 using polarized basis sets and geometries optimized with DFT methods. The present theoretical results support an electronic fragmentation energy D(sub e) around 10-11 eV in disagreement with most experimental results that place the dissociation energy D(sub o) (including zero point energy) around 7-8 eV. The plausible errors remaining in the theoretical calculations are unlikely to account for this big difference (2-4 eV).

Boese, A. Daniel↗

Theory of dense hydrogen

Equations of state for molecular and metallic hydrogen are calculated in an accurate and comparable manner by using the density-functional method applied to static lattices. By isolating the electronic energies in this way, attention can be drawn to the relative importance of the protonic degrees of freedom. From the equation-of-state results it can be concluded that a remnant molecular pairing is preferred in a band-overlap metallic state, complete dissociation occurring only at very high densities. An accurate determination of the pressure needed to achieve complete dissociation is shown to require a self-consistent treatment of electron and proton degrees of freedom.

Chakravarty, S.↗

The initiation and growth of delaminations induced by matrix microcracks in laminated composites

A recent variational mechanics analysis of microcracking damage in cross-ply laminates of the form /(S)/90n/s, where (S) is any orthotropic sublaminate much stiffer than /90n/, has been extended to account for the presence of delaminations emanating from the tips of microcracks in the /90 2n/T sublaminate. The new two-dimensional stress analysis is used to calculate the total strain energy, effective modulus, and longitudinal thermal expansion coefficient for a laminate having microcracks and delaminations. These results are used to calculate the energy release rate for the initiation and growth of a delamination induced by a matrix microcrack. At low crack densities, /(S)/90n/s laminates are expected to fail by microcracking and to show little or no delamination. At some critical crack density, which is a function of laminate structure and material properties, the energy release rate for delamination exceeds that for microcracking and delamination is predicted to dominate over microcracking. A quasi-three-dimensional model is used to predict the propagation of arbitrarily shaped delamination fronts. All predictions agree with experimental observations.

Nairn, J. A.↗

Pressure-strain-rate events in homogeneous turbulent shear flow

A detailed study of the intercomponent energy transfer processes by the pressure-strain-rate in homogeneous turbulent shear flow is presented. Probability density functions (pdf's) and contour plots of the rapid and slow pressure-strain-rate show that the energy transfer processes are extremely peaky, with high-magnitude events dominating low-magnitude fluctuations, as reflected by very high flatness factors of the pressure-strain-rate. A concept of the energy transfer class was applied to investigate details of the direction as well as magnitude of the energy transfer processes. In incompressible flow, six disjoint energy transfer classes exist. Examination of contours in instantaneous fields, pdf's and weighted pdf's of the pressure-strain-rate indicates that in the low magnitude regions all six classes play an important role, but in the high magnitude regions four classes of transfer processes, dominate. The contribution to the average slow pressure-strain-rate from the high magnitude fluctuations is only 50 percent or less. The relative significance of high and low magnitude transfer events is discussed.

Brasseur, James G.↗

Electron escape depth variation in thin SiO2 films measured with variable photon energy

A double crystal monochromator at the Stanford Synchrotron Radiation Laboratory is used to study the Si/SiO2 interface, using photon energies of hv = 1950-3700 eV. This photon energy range allows interfaces to be observed through oxide layers 50 A thick or more. Variations in electron escape depth and/or oxide density as a function of distance from the interface are observed over the entire kinetic energy range (100-3600 eV). These differences are attributed to a strained oxide layer near the interface.

Hecht, M. H.↗

Electron beam experiments at high altitudes

Experiments with the electron gun on the SCATHA satellite produced evidence of beam-plasma interactions, and heating of the low energy electrons around the satellite. These experiments were conducted near geosynchronous orbit, in the dusk bulge, and plasma sheet, with one short operation in the lobe regions, providing a range of ambient plasma densities. The electron gun was operated at 50 eV, with beam currents of 1, 10, and 100 micro-A. Data from electrostatic analyzers and the dc electric field experiment show that the satellite charged to near the beam energy in sunlight, if the beam current was sufficient. Higher ambient densities required higher beam currents. The electrostatic analyzers showed distribution functions which had peaks, or plateaus, at energies greater than the satellite potential. These measurements indicate heating of the ambient plasma at several Debye lengths from the satellite, with the heated plasma then accelerated into the satellite. It is likely that the ambient plasma is in fact the photoelectron sheath generated by the satellite.

Olsen, R. C.↗

Electron beam experiments at high altitudes

Experiments with the electron gun on the SCATHA satellite produced evidence of beam-plasma interactions, and heating of the low energy electrons around the satellite. These experiments were conducted near geosynchronous orbit, in the dusk bulge, and plasma sheet, with one short operation in the lobe regions, providing a range of ambient plasma densities. The electron gun was operated at 50 eV, with beam currents of 1, 10, and 100 micro-A. Data from electrostatic analyzers and the DC electric field experiment show that the satellite charged to near the beam energy in sunlight, if the beam current was sufficient. Higher ambient densities required higher beam currents. The electrostatic analyzers showed distribution functions which had peaks, or plateaus, at energies greater than the satellite potential. These measurements indicate heating of the ambient plasma at several Debye lengths from the satellite, with the heated plasma then accelerated into the satellite. It is likely that the ambient plasma is in fact the photoelectron sheath generated by the satellite.

Olsen, R. C.↗

Clustering correlations and limits on cosmological gravitational waves

A cosmological background of gravitational waves induces angular deviations in the propagation of light traversing it. All observed astrophysical sources might therefore have varying apparent positions, with the time dependence set by the wave period. Wavelengths greater than a kiloparsec are examined, so the positions are frozen but in general correlated. Comparison with observed galaxy-galaxy n-point correlation functions provide a spectrum-independent limit on the energy density of gravitational waves with wavelengths between a few tens of kiloparsecs and a few hundred megaparsecs of Omega(GW) less than 0.001.

Linder, Eric V.↗

An Ab Initio Study of CuCO

Modified coupled-pair functional (MCPF) calculations and coupled cluster singles and doubles calculations, which include a perturbational estimate of the connected triples [CCSD(T)], yield a bent structure for CuCO, thus, supporting the prediction of a nonlinear structure based on density functional (DF) calculations. Our best estimate for the binding energy is 4.9 +/- 1.4 kcal/mol; this is in better agreement with experiment (6.0 +/- 1.2 kcal/mol) than the DF approach which yields a value (19.6 kcal/mol) significantly larger than experiment.

Bauschlicher, Charles W., Jr.↗

Application of the embedded atom method to Pb and Be

We have derived the embedding energy functional and two-body potential of the embedded atom method using decreasing exponentials for both the electron density and the two body potential. The embedding function was obtained from the equation of state given by Rose et al. (1984). Because of the form of the embedding function, the equilibrium lattice constant, cohesive energy, and bulk modulus are automatically satisfied. The two parameters phi(e) and gamma of the two-body potential were determined by fitting to shear modulus and the single vacancy formation energy. Contributions of up to the third nearest neighbors were included in the evaluation of the charge density rho and the two-body potential phi. The stability and anisotropy of each structure were estimated and compared with the available experimental data.

Karimi, M.↗

Equatorial electron energy and number densities in the Jovian magnetosphere

A synchrotron model with a Maxwellian energy distribution of the form e to the (-E/E sub 0) power is used in a comparison with spatially resolved radio interferometric measurements of the Jovian emission. The observations of the decimeter radiation as a function of equatorial distance at 10.4 and 21 cm wavelength were reduced to source emission/cc of source electrons in each of 16 concentric rings. The peak energies for isotropically distributed electrons exceeded the maximum energy for flat orbiting electrons, and the peaks were generally located from 2.25 to 3 Jupiter radii. Beyond 3 radii, the order of magnitude on number density became a sensitive function of pitch angle distribution. The total equatorial intensities at 75 cm wavelength were computed for (E sub 0)(r) and n(r) at different values of B sub 0. The radiative half life for electrons of initial energy E sub 0 in a dipole field was calculated and found to be nearly constant at one year or less for altitudes at and below the position in peak energy.

Luthey, J. L.↗

A redshift survey of IRAS galaxies. III - Reconstruction of the velocity and density fields in N-body model universes

N-body simulations of a 'cold dark matter' universe are presently used to calibrate the accuracy, and assess the limitations, of the procedure previously employed to predict the velocity field within 8000 km/sec of the Local Group through the application of linear gravitational theory to a full-sky, flux-limited sample of IRAS galaxies. The rms difference between the one-dimensional acceleration and velocity of field particles is an increasing function of local density; linear theory can in this way account for all but one-sixth of kinetic energy. A series of artificial IRAS catalogs closely matching the real sample in space density and clustering amplitude is constructed. Velocity correlation functions are used to demonstrate that the predicted velocity fields are in good agreement with the true velocity fields on large scales.

Davis, Marc↗

Some strategies for quantitative scanning Auger electron microscopy

The general applicability of power law forms of the background in electron spectra is pointed out and exploited for background removal from under Auger peaks. This form of B(E) is found to be extremely sensitive to instrumental alignment and to fault-free construction - an observation which can be used to set up analyser configurations in an accurate way. Also, differences between N(E) and B(E) can be used to derive a spectrometer transmission function T(E). The questions of information density in an energy-analysing spatially-resolving instrument are addressed after reliable instrumental characterization has been established. Strategies involving ratio histograms, showing the population distribution of the ratio of a pair of Auger peak heights, composition scatter diagrams and windowed imaging are discussed and illustrated.

Browning, R.↗

Turbulent small-scale neutral and ion density fluctuations as measured during MAC/Epsilon

During the MAC/Epsilon campaign (Fall 1987, from Andoya, Northern Norway, 69 N, 16 E) a total of four altitude profiles of neutral gas number densities and six profiles of ion number densities were measured with high spatial resolution in the height range from 60 to 120 km. First results of these rocket-borne experiments are presented with emphasis on small scale turbulent density variations and related turbulent parameter as structure function constants and energy dissipation rates.

Luebken, F.-J.↗

Collisionless Plasma Modeling in an Arbitrary Potential Energy Distribution

A new technique for calculating a collisionless plasma along a field line is presented. The primary feature of the new model is that it can handle an arbitrary (including nonmonotonic) potential energy distribution. This was one of the limiting constraints on the existing models in this class, and these constraints are generalized for an arbitrary potential energy composition. The formulation for relating current density to the field-aligned potential as well as formulas for density, temperature and energy flux calculations are presented for several distribution functions, ranging from a bi-Lorentzian with a loss cone to an isotropic Maxwellian. A comparison of these results with previous models shows that the formulation reduces.to the earlier models under similar assumptions.

Liemohn, M. W.↗