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 181 records · Page 10

Energy Landscape, and Structural and Spectroscopic Characterization of Diazirine and Its Cyclic Isomers

Identifying new nitrogenated hydrocarbon molecules in the interstellar medium (ISM) is challenging due to the lack of comprehensive spectroscopic data from experiments. In this computational work, we focus on investigating the structures, relative energies, spectroscopic constants, and energy landscape of the cyclic isomers of diazirine (c-CH2N2) using ab initio quantum chemical methods. Density functional theory (DFT) methods and coupled cluster theory with singles and doubles including perturbative triples [CCSD(T)] and CCSD(T) with the explicitly correlated F12 b correction [CCSD(T)-F12b] were employed for this purpose along with large correlation consistent cc-pVTZ, cc-pVQZ, and cc-pV5Z basis sets. Harmonic vibrational frequencies, infrared vibrational intensities, rotational constants, and dipole moments are reported. An harmonic vibrational fundamentals along with centrifugal distortion constants, and vibration-rotation interaction constants are also reported for all the cyclic isomers. The energies computed with the CCSD(T) and CCSD(T)-F12b methods were extrapolated to the one-particle complete basis set (CBS) limit following a three-point formula. At the CCSD(T)-F12b/CBS level of theory the 3,3H-diazirine (c-CH2N2) is the lowest energy cyclic isomer followed by 1,3H-diazirine, (E)-1,2H-Diazirine, and (Z)-1,2H-diazirine which are 20.1, 47.8, and 51.3kcal mol-1above the 3,3H-Diazirine, respectively. Accurate structures and spectroscopic constants that are reported here could be useful for future identification of these cyclic nitrogenated organic molecules in the interstellar medium or circumstellar disks.

Energy Landscape↗

Geometric model from microscopic theory for nuclear absorption

A parameter-free geometric model for nuclear absorption is derived herein from microscopic theory. The expression for the absorption cross section in the eikonal approximation, taken in integral form, is separated into a geometric contribution that is described by an energy-dependent effective radius and two surface terms that cancel in an asymptotic series expansion. For collisions of light nuclei, an expression for the effective radius is derived from harmonic oscillator nuclear density functions. A direct extension to heavy nuclei with Woods-Saxon densities is made by identifying the equivalent half-density radius for the harmonic oscillator functions. Coulomb corrections are incorporated, and a simplified geometric form of the Bradt-Peters type is obtained. Results spanning the energy range from 1 MeV/nucleon to 1 GeV/nucleon are presented. Good agreement with experimental results is obtained.

John, Sarah↗

Ponderomotive effects on distributions of O(+) ions in the auroral zone

Test particle calculations are used to compute the effects of gravity and ponderomotive acceleration by shear Alfven wave oscillations on the distribution function of O(+) ions along auroral field lines, assuming an ionospheric Maxwellian source of the ions at 2000 km altitude with approximately 0.5 eV of thermal energy in the parallel component of velocity. The electric field model corresponds to a standing wave oscillation with a frequency approximately 1 Hz in the azimuthal direction superimposed on the background dipole field, in which the wave amplitude is either increasing or decreasing in time. The electric field is taken to be primarily in the perpendicular direction. The time varying wave produces broad distributions with widths of 2 to 10 times the initial 0.5-eV thermal energy of the Maxwellian source, and the density and flux of upward going O(+) ions at one Earth radius are both enhanced in this model. The oxygen ion distribution functions at 1 R(sub E) altitude resulting from interaction with waves whose amplitudes are increasing in time have a more gradual lower energy cutoff than do the distribution functions resulting from decaying waves. The high-energy part of the distribution functions in growing waves reflects the temperature of the Maxwellian source, while the high-energy part of the distributions resulting from decaying waves steepens with time, independent of the source temperature.

Witt, E.↗

Galactic propagation models consistent with the cosmic ray lifetime derived from Be-10 measurements

Using a propagation calculation with energy dependent parameters, including the depletion of short pathlengths, and incorporating experimental nuclear excitation functions, the variation of the Be-10/Be9 ratio with the matter densities in two nested confinement regions is investigated. It is shown that there is no unique correspondence between a Be-10/Be9 measurement at low energy and the density of matter in the galaxy. Be-10/Be9 measurements at both low and high energy are needed to fully specify the matter densities.

Guzik, T. G.↗

Evolution of interstellar pickup ions in the solar wind

A model is constructed for the evolution of an interstellar pickup ion distribution in the solar wind. The model assumes that the ions are immediately isotropized at ionization and follows the subsequent development of the distribution function as the particles are convected with the solar wind. The effects of energy diffusion in an ambient wave field with a power law spectrum, adiabatic deceleration in the expanding solar wind, and continual addition of newly ionized particles are all included in the model. An analytical expression describing the evolution of the distribution function in phase space velocity and heliocentric radius is obtained. The distribution quickly approaches an asymptotic shape in phase space which depends on the relative efficiency of the energy diffusion process compared to that of adiabatic deceleration. At large distances from the sun the density of pickup ions falls as 1/r in this model. An expression for the distribution function at large distances and for large particle speed is presented. The asymptotic shape should describe the distribution of pickup ions in the outer heliosphere and could be used as an input distribution for a model of the anomalous component of heliosphere and could be used as an input distribution for a model of the anomalous component of cosmic rays. Comparison of this work with the recent observation of He(+) at 1 AU implies that the energy diffusion process is very weak inside 1 AU.

Isenberg, Philip A.↗

AladynPi – Adaptive Neural Network Molecular Dynamics Simulation Code with Physically Informed Potential: Computational Materials Mini-Application

This report provides an overview and description of commands used in the Computational Materials mini-application, AladynPi. AladynPi is an extension of a previously released mini-application, Aladyn (https://github.com/nasa/aladyn; Yamakov, V.I., and Glaessgen, E.H., NASA/TM-2018-220104). Aladyn and AladynPi are basic molecular dynamics codes written in FORTRAN 2003, which are designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory method. An input for the ANN is a set of structure coefficients, characterizing the local atomic environment of each atom, for which the atomic energy is obtained in the ANN inference process. In Aladyn, the ANN gives directly the energy of interatomic interactions. In AladynPi, the ANN gives optimized parameters for a predefined empirical function, known as bond-order-potential (BOP). The parameterized BOP function is then used to calculate the energy. AladynPi code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing manycore architectures based on graphic processing units (GPUs). The effort is supported by the High Performance Computing incubator (HPCi) project at NASA Langley Research Center.

Yamakov, Vesselin I.↗

A comparison of two types of structural optimization procedures for satisfying flutter requirements

A comparison is made of results obtained from mathematical programming and optimality criteria procedures for the minimum-mass design of typical aircraft wing structures to satisfy prescribed flutter requirements. The mathematical programming method is based on an interior penalty function approach. A rigorous optimality criterion and an intuitive optimality criterion based on uniform strain energy density are considered. An intuitive resizing procedure is used for both optimality criterion solutions. All results are calculated using the same computer program, changing only the optimization procedure. Both high- and low-aspect-ratio wings are examined. Finite elements are used for structural modeling, and the generalized coordinates for the flutter solution are based on the natural vibration modes of the structure.

Haftka, R. T.↗

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

Kapton charging characteristics: Effects of material thickness and electron-energy distribution

Charging characteristics of polyimide (Kapton) of varying thicknesses under irradiation by a very-low-curent-density electron beam, with the back surface of the sample grounded are reported. These charging characteristics are in good agreement with a simple analytical model which predicts that in thin samples at low current density, sample surface potential is limited by conduction leakage through the bulk material. The charging of Kapton in a low-current-density electron beam in which the beam energy was modulated to simulate Maxwellian and biMaxwellian distribution functions is measured.

Williamson, W. S.↗

Direct numerical simulation of sheared turbulent flow

The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.

Harris, Vascar G.↗

Geometric model for nuclear absorption from microscopic theory

A parameter-free geometric model for nuclear absorption is derived from microscopic theory. The expression for the absorption cross section in the eikonal approximation taken in integral form is separated into a geometric contribution, described by an energy-dependent effective radius, and two surface terms which are shown to cancel in an asymptotic series expansion. For collisions of light nuclei, an expression for the effective radius is derived using harmonic-oscillator nuclear density functions. A direct extension to heavy nuclei with Woods-Saxon densities is made by identifying the equivalent half density radius for the harmonic-oscillator functions. Coulomb corrections are incorporated and a simplified geometric form of the Bradt-Peters type obtained. Results spanning the energy range of 1 MeV/nucleon to 1 GeV/nucleon are presented. Good agreement with experimental results are obtained.

John, S.↗

Constraints on the Energy Density Content of the Universe Using Only Clusters of Galaxies

We demonstrate that it is possible to constrain the energy content of the Universe with high accuracy using observations of clusters of galaxies only. The degeneracies in the cosmological parameters are lifted by combining constraints from different observables of galaxy clusters. We show that constraints on cosmological parameters from galaxy cluster number counts as a function of redshift and accurate angular diameter distance measurements to clusters are complementary to each other and their combination can constrain the energy density content of the Universe well. The number counts can be obtained from X-ray and/or SZ (Sunyaev-Zeldovich effect) surveys, the angular diameter distances can be determined from deep observations of the intra-cluster gas using their thermal bremsstrahlung X-ray emission and the SZ effect (X-SZ method). In this letter we combine constraints from simulated cluster number counts expected from a 12 deg2 SZ cluster survey and constraints from simulated angular diameter distance measurements based on using the X-SZ method assuming an expected accuracy of 7% in the angular diameter distance determination of 70 clusters with redshifts less than 1.5. We find that R, can be determined within about 25%, A within 20%, and w within 16%. Any cluster survey can be used to select clusters for high accuracy distance measurements, but we assumed accurate angular diameter distance measurements for only 70 clusters since long observations are necessary to achieve high accuracy in distance measurements. Thus the question naturally arises: How to select clusters of galaxies for accurate diameter distance determinations? In this letter, as an example, we demonstrate that it is possible to optimize this selection changing the number of clusters observed, and the upper cut off of their redshift range. We show that constraints on cosmological parameters from combining cluster number counts and angular diameter distance measurements, as opposed to general expectations, will not improve substantially selecting clusters with redshifts higher than one. This important conclusion allow us to restrict our cluster sample to clusters closer than one, in a range where the observational time for accurate distance measurements are more manageable. Subject headings: cosmological parameters - cosmology: theory - galaxies: clusters: general - X-rays: galaxies: clusters

Molnar, Sandor M.↗

The electronic Hamiltonian for cuprates

A realistic many-body Hamiltonian for the cuprate superconductors should include both copper d and oxygen p states, hopping matrix elements between them, and Coulomb energies, both on-site and inter-site. We have developed a novel computational scheme for deriving the relevant parameters ab initio from a constrained occupation local density functional. The scheme includes numerical calculation of appropriate Wannier functions for the copper and oxygen states. Explicit parameter values are given for La2CuO4. These parameters are generally consistent with other estimates and with the observed superexchange energy. Secondly, we address whether this complicated multi-band Hamiltonian can be reduced to a simpler one with fewer basis states per unit cell. We propose a mapping onto a new two-band effective Hamiltonian with one copper d and one oxygen p derived state per unit cell. This mapping takes into account the large oxygen-oxygen hopping given by the ab initio calculations.

Annett, James F.↗

Evolution of the Suprathermal Proton Population at Interplanetary Shocks

We investigate the evolution of the suprathermal (ST) proton population as interplanetary shocks cross 1 au. The variability of the ST proton intensities and energy spectra upstream of the shocks is analyzed in terms of the shock parameters, upstream magnetic field configurations, and preexisting upstream populations. Propitious conditions for the observation of ST particles at distances far upstream from the shock occur in parallel shock configurations when particles can easily escape from the shock vicinity. In this situation, ST intensity enhancements show onsets characterized by velocity dispersion effects and energy spectra that develop into a “hump” profile peaking around ∼10 keV just before the arrival of the shock. The observation of field-aligned proton beams at low energies (5–10 keV) is possible under conditions that facilitate the scatter-free propagation of the particles streaming out of the shock. Upstream of perpendicular shocks, ST intensity enhancements are only observed in close proximity to the shock. Power-law proton spectra develop downstream of the shocks. The functional form for the downstream phase-space density proportional to v(exp -5) is observed only over a limited range of ST energies. The absence of ST populations observed far upstream of interplanetary shocks raises questions about whether ST protons contribute as a seed particle population in the processes of particle acceleration at shocks.

Shock waves↗

Acceleration of low-energy protons and alpha particles at interplanetary shock waves

The low-energy protons and alpha particles in the energy range 30 keV/charge to 150 keV/charge associated with three different interplanetary shock waves in the immediate preshock and postshock region are studied using data obtained by the ISEE 3. The spatial distributions in the preshock and postshock medium are presented, and the dependence of the phase space density at different energies on the distance from the shock and on the form of the distribution function of both species immediately at the shock is examined. It is found that in the preshock region the particles are flowing in the solar wind frame of reference away from the shock and in the postshock medium the distribution is more or less isotropic in this frame of reference. The distribution function in the postshock region can be represented by a power law in energy which has the same spectral exponent for both protons and alpha particles. It is concluded that the first-order Fermi acceleration process can consistently explain the data, although the spectra of diffuse bow shock associated particles are different from the spectra of the interplanetary shock-associated particles in the immediate vicinity of the shock. In addition, the mean free path of the low energy ions in the preshock medium is found to be considerably smaller than the mean free path determined by the turbulence of the background interplanetary medium.

Scholer, M.↗

On the pulsational overstability in narrowly confined viscous rings

Viscous accretion disks and rings are subject to pulsational overstability if the effective viscous stress is a rapidly increasing function of the disks' (or rings') surface density. Such a phenomenon is induced by the transfer of energy from shear to normal mode oscillations through viscous stress. The paper examines the global response of viscous accretion disks or rings to pulsational overstabilities. The analytic and preliminary numerical results indicate that all modes with wavelength longer than the scale height of the disk or ring are overstable. From the results, some examples of the complex time-dependent structure which may be induced by this overstability are illustrated. Finally, results are applied to narrow planetary rings to address the issues of maintenance of uniform precession and origin of eccentricity.

Papaloizou, J. C. B.↗

Langmuir Probe Measurements of Inductively Coupled Plasma in CF4/AR/O2 Mixtures

Inductively coupled plasmas of CF4:Ar:O2, which have been of importance to material processing, were studied in the GEC cell at 80:10:10, 60:20:20, and 40:30:30 mixture ratios. Radial distributions of plasma potential (V(sub p)), electron and ion number densities (n(sub e) and n(sub i)), electron temperature (T(sub e)), and electron energy distribution functions (EEDFs) were measured in the mid-plane of plasma across the electrodes in the pressure range of 10-50 mTorr, and RF power of 200 and 300 W. V(sub p), n(sub e) and n(sub i), which peak in the center of the plasma, increase with decrease of pressure. T(sub e) also increases with pressure but peaks toward the electrode edge. Both V(sub p) and T(sub e) remain nearly independent of RF power, whereas n(sub e) and n(sub i) increase with power. In all conditions the EEDFs exhibit non-Maxwellian shape and are more like Druyvesteyn form at higher energies. They exhibit a broad lip in the energy range 0-10 eV suggesting an electron loss mechanism, which could be due to via resonance electron attachment processes producing negative ions in this rich electronegative gas mixture. This behavior is more prominent towards the electrode edge.

Rao, M. V. V. S.↗

Langmuir Probe Measurements of Inductively Coupled Plasmas in CF4/Ar/O2 Mixtures

Inductively coupled plasmas of CF4:Ar:O2, which have been of importance to material processing, were studied in the GEC cell at 80:10:10, 60:20:20, and 40:30:30 mixture ratios. Radial distributions of plasma potential (V(sub p)), electron and ion number densities (n(sub e) and n(sub i), electron temperature (T(sub e)), and electron energy distribution functions (EEDFs) were measured in the mid-plane of plasma across the electrodes in the pressure range of 10-50 mTorr, and RF (radio frequency) power of 200 and 300 W. V(sub p), n(sub e) and n(sub i), which peak in the center of the plasma, increase with decrease of pressure. T(sub e) also increases with pressure but peaks toward the electrode edge. Both V(sub p) and T(sub e) remain nearly independent of RF power, whereas n(sub e) and n(sub i) increase with power. In all conditions the EEDFs exhibit non-Maxwellian shape and are more like Druyvesteyn form at higher energies. They exhibit a broad dip in the energy range 0-10 eV suggesting an electron loss mechanism, which could be due to via resonance electron attachment processes producing negative ions in this rich electronegative gas mixture. This behavior is more prominent towards the electrode edge.

Rao, M. V. V. S.↗