Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ab initio modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Ab Initio Modeling of Molecular Radiation

Radiative emission from excited states of atoms and molecules can comprise a significant fraction of the total heat flux experienced by spacecraft during atmospheric entry at hypersonic speeds. For spacecraft with ablating heat shields, some of this radiative flux can be absorbed by molecular constituents in the boundary layer that are formed by the ablation process. Ab initio quantum mechanical calculations are carried out to predict the strengths of these emission and absorption processes. This talk will describe the methods used in these calculations using, as examples, the 4th positive emission bands of CO and the 1g+ 1u+ absorption in C3. The results of these calculations are being used as input to NASA radiation modeling codes like NeqAir, HARA and HyperRad.

spectroscopy↗
NASA NTRS2009

A Model for Predicting Thermoelectric Properties of Bi2Te3

A parameterized orthogonal tight-binding mathematical model of the quantum electronic structure of the bismuth telluride molecule has been devised for use in conjunction with a semiclassical transport model in predicting the thermoelectric properties of doped bismuth telluride. This model is expected to be useful in designing and analyzing Bi2Te3 thermoelectric devices, including ones that contain such nano - structures as quantum wells and wires. In addition, the understanding gained in the use of this model can be expected to lead to the development of better models that could be useful for developing other thermoelectric materials and devices having enhanced thermoelectric properties. Bi2Te3 is one of the best bulk thermoelectric materials and is widely used in commercial thermoelectric devices. Most prior theoretical studies of the thermoelectric properties of Bi2Te3 have involved either continuum models or ab-initio models. Continuum models are computationally very efficient, but do not account for atomic-level effects. Ab-initio models are atomistic by definition, but do not scale well in that computation times increase excessively with increasing numbers of atoms. The present tight-binding model bridges the gap between the well-scalable but non-atomistic continuum models and the atomistic but poorly scalable ab-initio models: The present tight-binding model is atomistic, yet also computationally efficient because of the reduced (relative to an ab-initio model) number of basis orbitals and flexible parameterization of the Hamiltonian.

Lee, Seungwon↗
NASA NTRS1991

Effects of Mg II and Ca II ionization on ab-initio solar chromosphere models

Acoustically heated solar chromosphere models are computed considering radiation damping by (non-LTE) emission from H(-) and by Mg II and Ca II emission lines. The radiative transfer equations for the Mg II k and Ca II K emission lines are solved using the core-saturation method with complete redistribution. The Mg II k and Ca II K cooling rates are compared with the VAL model C. Several substantial improvements over the work of Ulmschneider et al. (1987) are included. It is found that the rapid temperature rises caused by the ionization of Mg II are not formed in the middle chromosphere, but occur at larger atmospheric heights. These models represent the temperature structure of the 'real' solar chromosphere much better. This result is a major precondition for the study of ab-initio models for solar flux tubes based on MHD wave propagation and also for ab-initio models for the solar transition layer.

Rammacher, W.↗
NASA NTRS2015

Integrated Physics-based Modeling and Experiments for Improved Prediction of Combustion Dynamics in Low-Emission Systems

Concurrent simulation and experiment was undertaken to assess the ability of a hybrid RANS-LES model to predict combustion dynamics in a single-element lean direct-inject (LDI) combustor showing self-excited instabilities. High frequency pressure modes produced by Fourier and modal decomposition analysis were compared quantitatively, and trends with equivalence ratio and inlet temperature were compared qualitatively. High frequency OH PLIF and PIV measurements were also taken. Submodels for chemical kinetics and primary and secondary atomization were also tested against the measured behavior. For a point-wise comparison, the amplitudes matched within a factor of two. The dependence on equivalence ratio was matched. Preliminary results from simulation using an 18-reaction kinetics model indicated instability amplitudes closer to measurement. Analysis of the simulations suggested a band of modes around 1400 Hz were due to a vortex bubble breakdown and a band of modes around 6 kHz were due to a precessing vortex core hydrodynamic instability. The primary needs are directly coupled and validated ab initio models of the atomizer free surface flow and the primary atomization processes, and more detailed study of the coupling between the 3D swirling flow and the local thermoacoustics in the diverging venturi section.

combustion↗
NASA NTRS1993

A generalized version of the Rankine-Hugoniot relations including ionization, dissociation and related phenomena

For purposes of computing shocks in stellars atmospheres and winds we have developed a generalized version of the Rankine-Hugoniot relations including ionization, dissociation, radiation and related phenomena such as excitation, rotation and vibration of molecules. The new equations are given in analytical form. They are valid as long as the internal energy E, the total pressure P, and the first adiabatic coefficient gamma(sub 1) can be evaluated. However, we have not treated shock structures. In the case of non-LTE we have to employ an approximation for gamma(sub 1) because in that case no definition exists. Our new version of the Rankine-Hugoniot relations can easily be used for many purposes including ab-initio modeling. In our derivation we introduce a parameter gamma(sub H), which is definded as the ratio of the enthalpy H (sometimes called heat function w) to the internal energy E (sometimes called U). Using this parameter we solve the equations for changing mu and (d(natural log P)/d(natural log rho))(sub ad) identically equal to gamma(sub 1) on both sides of the shock. Both gamma(sub H) and gamma(sub 1), and also mu are functions of pressure P and temperature T. We present: (1) the derivation, (2) examples of gamma(sub 1) (P,T) and gamma(sub H) (P,T) which include/exclude ionization and radiation, and (3) as an example the differences in post-shock parameters as function of the pre-shock temperature for the case with ionization and without radiation.

Nieuwenhuijzen, H.↗
NASA NTRS2018

Investigations of Spot-to-plume Mode Transition in a Hollow Cathode Discharge Using 2-D Axisymmetric Plasma Simulations

After decades of laboratory investigations that provided invaluable measurements and insight, the physics behind the transition from spot to plume modes in hollow cathodes remains one of the longest standing theoretical problems in electric propulsion. This has prohibited the development of ab initio models that allow for the prediction of the transition across different cathodes and operating conditions. Since the beginning of its development over a decade ago, simulations with the 2-D axisymmetric Orificed Cathode (OrCa2D) code have helped elucidate a wide range of processes in hollow cathode discharges. However, the code has never been used to investigate the onset of plume mode. We present results from the first OrCa2D simulations of a 25-A LaB6 cathode for a range of flow rates (5-20 sccm) in which transition from spot to plume modes is known to occur. The cathode in this study was one of the two technologies considered for the 12.5 kW Hall Effect Rocket with Magnetic Shielding (HERMeS) and operates nominally at 21 A and 14.8 sccm. The simulations capture the characteristic rise of the peak-to-peak amplitude in the keeper voltage oscillations and underscore the significance of the plume neutral gas in the transition. The plasma inside the cathode is found to be relatively quiescent throughout the transition, in agreement with previous experimental observations. The computed keeper voltage fluctuations at low flow rates (<8 sccm) are found to be driven by oscillations of the same frequency in the plasma plume with the following main characteristics: (1) they are of low frequency (<10 kHz), and associated with small longitudinal motion in the direction of the applied magnetic field, (2) they occur in a region of the plume where the neutral gas provided by the cathode has been fully depleted, and (3) they have a (small) wave velocity of about 100 m/s, which is at least ~10 smaller than the drift, thermal and acoustic speeds of the ions. At 8 sccm, when the transition to the large-amplitude oscillations begins, the ionization frequency in the neutral-depleted plume region ranges ~2-100 kHz. The frequency of the oscillations in the plasma (and keeper voltage) is found to be equal to the ionization frequency (~5 kHz) at the center of this region. The findings suggest that the transition to plume mode is driven by ionization processes in the near-plume of the cathode, in line with previous conjectures that were based solely on laboratory observations.

Polk, James E.↗
NASA NTRS1992

Ab initio atomic recombination reaction energetics on model heat shield surfaces

Ab initio quantum mechanical calculations on small hydration complexes involving the nitrate anion are reported. The self-consistent field method with accurate basis sets has been applied to compute completely optimized equilibrium geometries, vibrational frequencies, thermochemical parameters, and stable site labilities of complexes involving 1, 2, and 3 waters. The most stable geometries in the first hydration shell involve in-plane waters bridging pairs of nitrate oxygens with two equal and bent hydrogen bonds. A second extremely labile local minimum involves out-of-plane waters with a single hydrogen bond and lies about 2 kcal/mol higher. The potential in the region of the second minimum is extremely flat and qualitatively sensitive to changes in the basis set; it does not correspond to a true equilibrium structure.

Senese, Fredrick↗
NASA NTRS1994

Conformational Characteristics of Poly(tetrafluoroethylene) (PTFE) Based Upon Ab Initio Electronic Structure Calculations on Model Molecules

Conformational energy contours of perfluoroalkanes, determined from ab initio calculations, confirm the well-known spitting of trans states into two minima at plus or minus 17 degrees but also show that the gauche states split as well, with minima at plus or minus 124 degrees and plus or minus 84 in order to relieve steric crowding. The directions of such split distortions from the perfectly staggered states are strongly coupled for adjacent pairs of bonds in a manner identical to the intradyad pair for poly (isobutylene) chains. These conformational characteristics are fully represented by a six-state rotational isomeric state (RIS) model for PTFE comprised of t(+), t(-), g(sup +)+, g(sup +)-, g(sup -) + and g(sup -)-states, located at the split energy minima. The resultant 6 x 6 statistical weight matrix is described by first-order interaction parameters for the g+(+) (ca. 0.6 kcal/mol) and g+- (ca. 2.0 kcal/mol) states, and second order parameters for the g(sup +)+g(sup +)+ (ca 0.6 kcal/mol) and g(sup +)+g(sup -)+ (ca. 1.0 kcal/mol) states. This six-state RIS model, without adjustment of the geometric or energy parameters as determined from the ab initio calculations, predicts the unperturbed chain dimensions and the fraction of gauche bonds as a function of temperature for PTFE in good agreement with available experimental values.

Smith, Grant D.↗
NASA NTRS1990

Stochastic shock waves as a candidate mechanism for the formation of the He I 10830-A line in cool giant stars

The question of whether strong shocks produced in time-dependent stochastic wave models can explain the formation of the He I 10830-A line in cool giant stars is examined. The research is based on the ab initio chromosphere model for Arcturus by Cuntz (1987), showing that a stochastic distribution of wave periods leads to overtaking and merging of shocks, which occasionally produce very strong shocks with temperatures larger than 40,000 K in the postshock regions. These temperatures can easily produce a significant population in the 2s 3S state by electron collisional excitation. The 10830-A line occurs in absorption when the densities in the shocks exceed 10 million/cu cm.

Cuntz, Manfred↗
NASA NTRS2019

Warm Absorber Diagnostics of AGN Dynamics

Warm absorbers and related phenomena are some of the observable manifestations of outflows or winds from active galactic nuclei (AGNs). Warm absorbers are common in low-luminosity AGNs. They have been extensively studied observationally and are well described by simple phenomenological models. However, major open questions remain. What is the driving mechanism? What is the density and geometrical distribution? How much associated fully ionized gas is there? What is the relation to the quasi-relativistic "ultrafast outflows"? In this paper we present synthetic spectra for the observable properties of warm absorber flows and associated quantities. We use ab initio dynamical models, i.e., solutions of the equations of motion for gas in finite difference form. The models employ various plausible assumptions for the origin of the warm absorber gas and the physical mechanisms affecting its motion. The synthetic spectra are presented as an observational test of these models. In this way we explore various scenarios for warm absorber dynamics. We show that observed spectra place certain requirements on the geometrical distribution of the warm absorber gas, and that not all dynamical scenarios are equally successful at producing spectra similar to what is observed.

Kallman, Timothy R.↗
NASA NTRS1977

The molecular composition of dense interstellar clouds

Presented in this paper is an ab initio chemical model for dense interstellar clouds that incorporates 598 grain surface reactions, with small grains providing the reaction area. Gas-phase molecules are depleted through collisions with grains. The abundances of 372 chemical species are calculated as a function of time and are found to be of sufficient magnitude to explain most observations. Peak abundances are achieved on time scales of the order of 100,000 to 1 million years, depending on cloud density and kinetic temperature. The reaction rates for ion-molecule chemistry are approximately the same, indicating that surface and gas-phase chemistry may be coupled in certain regions. The composition of grain mantles is shown to be a function of grain radius. In certain grain-size ranges, large molecules containing two or more heavy atoms are more predominant than lighter 'ices' - H2O, NH3, and CH4. It is possible that absorption due to these large molecules in the mantle may contribute to the observed 3-micron band in astronomical spectra.

Allen, M.↗
NASA NTRS2014

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗
NASA NTRS2014

MHD Turbulence and Magnetic Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V↗
NASA NTRS2016

Serious Gaming for Test & Evaluation of Clean-Slate (Ab Initio) National Airspace System (NAS) Designs

Incremental approaches to air transportation system development inherit current architectural constraints, which, in turn, place hard bounds on system capacity, efficiency of performance, and complexity. To enable airspace operations of the future, a clean-slate (ab initio) airspace design(s) must be considered. This ab initio National Airspace System (NAS) must be capable of accommodating increased traffic density, a broader diversity of aircraft, and on-demand mobility. System and subsystem designs should scale to accommodate the inevitable demand for airspace services that include large numbers of autonomous Unmanned Aerial Vehicles and a paradigm shift in general aviation (e.g., personal air vehicles) in addition to more traditional aerial vehicles such as commercial jetliners and weather balloons. The complex and adaptive nature of ab initio designs for the future NAS requires new approaches to validation, adding a significant physical experimentation component to analytical and simulation tools. In addition to software modeling and simulation, the ability to exercise system solutions in a flight environment will be an essential aspect of validation. The NASA Langley Research Center (LaRC) Autonomy Incubator seeks to develop a flight simulation infrastructure for ab initio modeling and simulation that assumes no specific NAS architecture and models vehicle-to-vehicle behavior to examine interactions and emergent behaviors among hundreds of intelligent aerial agents exhibiting collaborative, cooperative, coordinative, selfish, and malicious behaviors. The air transportation system of the future will be a complex adaptive system (CAS) characterized by complex and sometimes unpredictable (or unpredicted) behaviors that result from temporal and spatial interactions among large numbers of participants. A CAS not only evolves with a changing environment and adapts to it, it is closely coupled to all systems that constitute the environment. Thus, the ecosystem that contains the system and other systems evolves with the CAS as well. The effects of the emerging adaptation and co-evolution are difficult to capture with only combined mathematical and computational experimentation. Therefore, an ab initio flight simulation environment must accommodate individual vehicles, groups of self-organizing vehicles, and large-scale infrastructure behavior. Inspired by Massively Multiplayer Online Role Playing Games (MMORPG) and Serious Gaming, the proposed ab initio simulation environment is similar to online gaming environments in which player participants interact with each other, affect their environment, and expect the simulation to persist and change regardless of any individual player's active participation.

Allen, B. Danette↗
NASA NTRS2017

Dissociation and Internal Excitation of Molecular Nitrogen Due to N + N2 Collisions Using Direct Molecular Simulation

In this work we present a molecular level study of N2+N collisions, focusing on excitation of internal energy modes and non-equilibrium dissociation. The computation technique used here is the direct molecular simulation (DMS) method and the molecular interactions have been modeled using an ab−initio potential energy surface (PES) developed at NASA's Ames Research Center. We carried out vibrational excitation calculations between 5000K and 30000K and found that the characteristic vibrational excitation time for the N + N2 process was an order of magnitude lower than that predicted by the Millikan and White correlation. It is observed that during vibrational excitation the high energy tail of the vibrational energy distribution gets over populated first and the lower energy levels get populated as the system evolves. It is found that the non-equilibrium dissociation rate coefficients for the N + N2 process are larger than those for the N2 + N2 process. This is attributed to the non-equilibrium vibrational energy distributions for the N + N2 process being less depleted than that for the N2 +N2 process. For an isothermal simulation we find that the probability of dissociation goes as 1/T(sub tr) for molecules with internal energy (epsilon(sub int)) less than approximately 9.9eV, while for molecules with epsilon (sub int) greater than 9.9eV the dissociation probability was weakly dependent on translational temperature of the system. We compared non-equilibrium dissociation rate coefficients and characteristic vibrational excitation times obtained by using the ab-initio PES developed at NASA's Ames Research Center to those obtained by using an ab-initio PES developed at the University of Minnesota. Good agreement was found between the macroscopic properties and molecular level description of the system obtained by using the two PESs.

Grover, Maninder S.↗
NASA NTRS1994

Equilibrium and Dynamics Properties of Poly(oxyethylene) Melts and Related Poly(alkylethers) from Simulations and Ab Initio Calculations

Molecular dynamics simulations of POE melts have been performed utilizing a potential force field parameterized to reproduce conformer energies and rotational energy barriers in dimethoxyethane as determined from ab initio electronic structure calculations. Chain conformations and dimensions of POE from the simulations were found to be in good agreement with predictions of a rotational isomeric state (RIS) model based upon the ab initio conformational. energies. The melt chains were found to be somewhat extended relative to chains at theta conditions. This effect will be discussed in light of neutron scattering experiments which indicate that POE chains are extended in the melt relative to theta solutions. The conformational characteristics of POE chains will also be compared with those of other poly (alkylethers), namely poly(oxymethylene), poly(oxytrimethylene) and poly(oxytetramethylene). Local conformational dynamics were found to be more rapid than in polymethylene. Calculated C-H vector correlation times were found to be in reasonable agreement with experimental values from C-13 NMR spin-lattice relaxation times. The influence of ionic salts on local conformations and dynamics will also be discussed.

Smith, Grant D.↗