Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ATOMIC MODELS”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Pairing preferences of the model mono-valence mono-atomic ions investigated by molecular simulation

We carried out a series of potential of mean force calculations to study the pairing preferences of a series of model mono-atomic 1:1 ions with evenly varied sizes. The probabilities of forming the contact ion pair (CIP) and the single water separate ion pair (SIP) were presented in the two-dimensional plots with respect to the ion sizes. The pairing preferences reflected in these plots largely agree with the empirical rule of matching ion sizes in the small and big size regions. In the region that the ion sizes are close to the size of the water molecule; however, a significant deviation from this conventional rule is observed. Our further analysis indicated that this deviation originates from the competition between CIP and the water bridging SIP state. The competition is mainly an enthalpy modulated phenomenon in which the existing of the water bridging plays a significant role.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Structure and aggregation in model tetramethylurea solutions

The structure of model aqueous tetramethylurea (TMU) solutions is investigated employing large-scale (32 000, 64 000 particles) molecular dynamics simulations. Results are reported for TMU mole fractions, X{sub t}, ranging from infinite dilution up to 0.07, and for two temperatures, 300 and 330 K. Two existing force fields for TMU-water solutions are considered. These are the GROMOS 53A6 united-atom TMU model combined with SPC/E water [TMU(GROMOS-UA)/W(SPC/E)], and the more frequently employed AMBER03 all-atom force field for TMU combined with the TIP3P water model [TMU(AMBER-AA)/W(TIP3P)]. It is shown that TMU has a tendency towards aggregation for both models considered, but the tendency is significantly stronger for the [TMU(AMBER-AA)/W(TIP3P)] force field. For this model signs of aggregation are detected at X{sub t} = 0.005, aggregation is a well established feature of the solution at X{sub t} = 0.02, and the aggregates increase further in size with increasing concentration. This is in agreement with at least some experimental studies, which report signals of aggregation in the low concentration regime. The TMU aggregates exhibit little structure and are simply loosely ordered, TMU-rich regions of solution. The [TMU(GROMOS-UA)/W(SPC/E)] model shows strong signs of aggregation only at higher concentrations (X{sub t} ≳ 0.04), and the aggregates appear more loosely ordered, and less well-defined than those occurring in the [TMU(AMBER-AA)/W(TIP3P)] system. For both models, TMU aggregation increases when the temperature is increased from 300 to 330 K, consistent with an underlying entropy driven, hydrophobic interaction mechanism. At X{sub t} = 0.07, the extra-molecular correlation length expected for microheterogeneous solutions has become comparable with the size of the simulation cell for both models considered, indicating that even the systems simulated here are sufficiently large only at low concentrations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Insights from high-fidelity modeling of industrial rotary bell atomization

The global automotive industry sprayed over 2.6 billion liters of paint in 2018, much of which through electrostatic rotary bell atomization, a highly complex process involving the fluid mechanics of rapidly rotating thin films tearing apart into micrometer-thin filaments and droplets. Coating operations account for 65% of the energy usage in a typical automotive assembly plant, representing 10,000s of gigawatt-hours each year in the United States alone. Optimization of these processes would allow for improved robustness, reduced material waste, increased throughput, and significantly reduced energy usage. Here, we introduce a high-fidelity mathematical and algorithmic framework to analyze rotary bell atomization dynamics at industrially relevant conditions. Our approach couples laboratory experiment with the development of robust non-Newtonian fluid models; devises high-order accurate numerical methods to compute the coupled bell, paint, and gas dynamics; and efficiently exploits high-performance supercomputing architectures. These advances have yielded insight into key dynamics, including i) parametric trends in film, sheeting, and filament characteristics as a function of fluid rheology, delivery rates, and bell speed; ii) the impact of nonuniform film thicknesses on atomization performance; and iii) an understanding of spray composition via primary and secondary atomization. These findings result in coating design principles that are poised to improve energy- and cost-efficiency in a wide array of industrial and manufacturing settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Multisite Proton–Coupled Electron Transfer at a Keggin-Type Polyoxotungstate

Proton−coupled electron transfer (PCET) governs many redox transformations, but is thermodynamically constrained when proton and electron transfer occur at a single site. Here, we introduce a new multisite PCET (MSPCET) platform, based on the Keggin-type polyoxotungstate, [VW 12 O 40 ] 3− (VW 12 ). Pairing VW 12 with either Brønsted bases or acids yields reagent pairs with tunable effective bond dissociation free energies (BDFE eff ) over 15 kcal mol −1 , enabling both oxidative and reductive H atom transfer reactions. Kinetic studies on the oxidative pathway by using 2,4,6- t Bu 3 PhOH as a model hydrogen atom (H atom) donor reveal a product-like, entropy-dominated concerted proton−electron transfer (CPET) pathway from a preorganized hydrogen-bonded complex. By contrast, reductive H atom transfer reactions exhibit larger ΔH ‡ values, measurable kinetic isotope effects, and balanced Brønsted slope, consistent with synchronous CPET-type mechanism. Extension to N−H, O−H, and C−H substrates demonstrates the versatility of the VW 12 MS-PCET platform for tunable (de)hydrogenation.

Charge transfer↗
DOE OSTI2014

Thermal transpiration: A molecular dynamics study

Thermal transpiration is a phenomenon where fluid molecules move from the cold end towards the hot end of a channel under the influence of longitudinal temperature gradient alone. Although the phenomenon of thermal transpiration is observed at rarefied gas conditions in macro systems, the phenomenon can occur at atmospheric pressure if the characteristic dimensions of the channel is less than 100 nm. The flow through these nanosized channels is characterized by the free molecular flow regimes and continuum theory is inadequate to describe the flow. Thus a non-continuum method like molecular dynamics (MD) is necessary to study such phenomenon. In the present work, MD simulations were carried out to investigate the occurance of thermal transpiration in copper and platinum nanochannels at atmospheric pressure conditions. The mean pressure of argon gas confined inside the nano channels was maintained around 1 bar. The channel height is maintained at 2nm. The argon atoms interact with each other and with the wall atoms through the Lennard-Jones potential. The wall atoms are modelled using an EAM potential. Further, separate simulations were carried out where a Harmonic potential is used for the atom-atom interaction in the platinum channel. A thermally insulating wall was introduced between the low and high temperature regions and those wall atoms interact with fluid atoms through a repulsive potential. A reduced cut off radius were used to achieve this. Thermal creep is induced by applying a temperature gradient along the channel wall. It was found that flow developed in the direction of the increasing temperature gradient of the wall. An increase in the volumetric flux was observed as the length of the cold and the hot regions of the wall were increased. The effect of temperature gradient and the wall-fluid interaction strength on the flow parameters have been studied to understand the phenomenon better.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Features of Calculation of the Equation of State, Composition, and Conductivity for a Plasma of Dense, Supercritical Metal Vapors—a Plasma Fluid

The caloric and thermal equations of state, composition, and conductivity have been calculated for a supercritical aluminum plasma fluid. A previously proposed chemical plasma model called the “3+”-component one was used for the calculations. The model includes atoms, electrons, ions, and an electron jellium. The thermodynamic functions have been calculated for the first time within the “3+” model for a plasma fluid. The magnification and compensation of intercharge and interatomic interactions when calculating the equation of state and composition are analyzed. The introduction of the jellium leads to an increase in conductivity under compression, while the compensation of interactions when calculating the composition leads to a virtually ideal-gas behavior of the equation of state. Comparison with the data from physical and numerical experiments has confirmed our conclusions and demonstrated that the hypothesis about jellium, a new gas–plasma component, is constructive.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Nuclear Structure Studies at the Limits of Stability

On the quest for a predictive model of atomic nuclei, the properties of rare isotopes have emerged as critical observables that benchmark extrapolations, identify missing physics ingredients in models, or that enable discovery science. Formerly NSCL and now FRIB facilitate world-leading experiments due to the availability of fast, stopped, and reaccelerated rare-isotope beams in conjunction with unique instruments. DOE-SC NP Award DE-SC0020451 supported from 2019 to 2023 MSU’s nuclear structure and reactions program that tackles compelling questions in the field with complementary techniques, probing single-particle and collective degrees of freedom in relevant regions of the nuclear chart and enabled developments that positioned the PIs and their research groups for first FRIB science.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1973

X-ray spectrum of Cassiopeia A - Evidence for iron line emission.

A sensitive measurement by rocket-borne detectors of the X-ray flux from Cas A has revealed a steep continuum and a broad spectral feature in the region where line radiation from iron nuclei would be expected. The flux in this feature is 0.0122 plus or minus 0.0017 photons per sq cm per sec. The presence of broad iron lines is consistent with a model in which about 13 MeV-per-nucleon iron nuclei charge exchange with surrounding interstellar oxygen and other heavy atoms. The model suggests that a substantial fraction of the energy from the outburst has gone into low-energy cosmic rays which produce the observed H II region surrounding the remnant.

Serlemitsos, P. J.↗
NASA NTRS2019

Physically-Informed Artificial Neural Networks for Atomistic Modeling of Materials

A new approach is presented for the development of classical interatomic potentials using physically-informed neural networks (PINN) combined with an analytical bond-order atomic interaction model. Due to the strong physical underpinnings, the PINN potentials demonstrate much better transferability than the existing machine-learning potentials while drastically improving the accuracy in comparison with traditional potentials. PINN potentials can be constructed for both metallic and covalent materials in a unified manner. A number of applications of PINN potentials to large-scale molecular dynamics and Monte Carlo simulations and calculation of thermal and mechanical properties of diverse materials are demonstrated. Some of the specific materials systems include silicon and aluminum, as well as alloys and compounds. Computational aspects of PINN potentials are discussed and future developments in this field are outlined.

J Hickman↗
DOE OSTI2021

Machine Learning Predictions of Transition Probabilities in Atomic Spectra

Forward modeling of optical spectra with absolute radiometric intensities requires knowledge of the individual transition probabilities for every transition in the spectrum. In many cases, these transition probabilities, or Einstein A-coefficients, quickly become practically impossible to obtain through either theoretical or experimental methods. Complicated electronic orbitals with higher order effects will reduce the accuracy of theoretical models. Experimental measurements can be prohibitively expensive and are rarely comprehensive due to physical constraints and sheer volume of required measurements. Due to these limitations, spectral predictions for many element transitions are not attainable. In this work, we investigate the efficacy of using machine learning models, specifically fully connected neural networks (FCNN), to predict Einstein A-coefficients using data from the NIST Atomic Spectra Database. For simple elements where closed form quantum calculations are possible, the data-driven modeling workflow performs well but can still have lower precision than theoretical calculations. For more complicated nuclei, deep learning emerged more comparable to theoretical predictions, such as Hartree–Fock. Unlike experiment or theory, the deep learning approach scales favorably with the number of transitions in a spectrum, especially if the transition probabilities are distributed across a wide range of values. It is also capable of being trained on both theoretical and experimental values simultaneously. In addition, the model performance improves when training on multiple elements prior to testing. The scalability of the machine learning approach makes it a potentially promising technique for estimating transition probabilities in previously inaccessible regions of the spectral and thermal domains on a significantly reduced timeline.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Analysis of the crystal lattice instability for cage–cluster systems using the superatom model

We have investigated the lattice dynamics for a number of rare-earth hexaborides based on the superatom model within which the boron octahedron is substituted by one superatom with a mass equal to the mass of six boron atoms. Phenomenological models have been constructed for the acoustic and lowenergy optical phonon modes in RB{sub 6} (R = La, Gd, Tb, Dy) compounds. Using DyB{sub 6} as an example, we have studied the anomalous softening of longitudinal acoustic phonons in several crystallographic directions, an effect that is also typical of GdB{sub 6} and TbB{sub 6}. The softening of the acoustic branches is shown to be achieved through the introduction of negative interatomic force constants between rare-earth ions. We discuss the structural instability of hexaborides based on 4f elements, the role of valence instability in the lattice dynamics, and the influence of the number of f electrons on the degree of softening of phonon modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

X-ray heating and electron temperature of laboratory photoionized plasmas

We discuss experimental and modeling results for the x-ray heating and temperature of laboratory photoionized plasmas. A method is used to extract electron temperature based on the analysis of transmission spectroscopy data that is independent of atomic kinetics modeling. Boltzmann kinetic and radiation-hydrodynamics simulations emphasized the critical interplay between atomic physics and heating in photoionized plasmas. They also demonstrated the dramatic impact of photoexcitation on excited state populations, line emissivity, and radiation cooling. Temperature calculations performed with astrophysical codes significantly overestimated the measurement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Characterization of organic aerosol across the global remote troposphere: a comparison of ATom measurements and global chemistry models

The spatial distribution and properties of submicron organic aerosol (OA) are among the key sources of uncertainty in our understanding of aerosol effects on climate. Uncertainties are particularly large over remote regions of the free troposphere and Southern Ocean, where very few data have been available and where OA predictions from AeroCom Phase II global models span 2 to 3 orders of magnitude, greatly exceeding the model spread over source regions. The (nearly) pole-to-pole vertical distribution of non-refractory aerosols was measured with an aerosol mass spectrometer onboard the NASA DC-8 aircraft as part of the Atmospheric Tomography (ATom) mission during the Northern Hemisphere summer (August 2016) and winter (February 2017). This study presents the first extensive characterization of OA mass concentrations and their level of oxidation in the remote atmosphere. OA and sulfate are the major contributors by mass to submicron aerosols in the remote troposphere, together with sea salt in the marine boundary layer. Sulfate was dominant in the lower stratosphere. OA concentrations have a strong seasonal and zonal variability, with the highest levels measured in the lower troposphere in the summer and over the regions influenced by biomass burning from Africa (up to 10 µg sm -3 ). Lower concentrations (~0.1–0.3 µg sm -3 ) are observed in the northern middle and high latitudes and very low concentrations (<0.1 µg sm -3 ) in the southern middle and high latitudes. The ATom dataset is used to evaluate predictions of eight current global chemistry models that implement a variety of commonly used representations of OA sources and chemistry, as well as of the AeroCom-II ensemble. The current model ensemble captures the average vertical and spatial distribution of measured OA concentrations, and the spread of the individual models remains within a factor of 5. These results are significantly improved over the AeroCom-II model ensemble, which shows large overestimations over these regions. However, some of the improved agreement with observations occurs for the wrong reasons, as models have the tendency to greatly overestimate the primary OA fraction and underestimate the secondary fraction. Measured OA in the remote free troposphere is highly oxygenated, with organic aerosol to organic carbon (OA/OC) ratios of ~2.2–2.8, and is 30 %–60 % more oxygenated than in current models, which can lead to significant errors in OA concentrations. The model–measurement comparisons presented here support the concept of a more dynamic OA system as proposed by Hodzic et al. (2016), with enhanced removal of primary OA and a stronger production of secondary OA in global models needed to provide better agreement with observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Quantum effect on the nucleation of plastic deformation carriers and destruction in crystals

New concepts on the irreversible crystal deformation as a structure transformation caused by a change in interatomic interactions at fluctuations of the electron density under loading are described. The change in interatomic interactions lead to the excitation of dynamical displacements of atoms. A model and a theory of a deformable pristine crystal taking into account the excitation of thermally activated and dynamical displacements of atoms are suggested. New mechanisms of the nucleation of plastic deformation carriers and destruction in pristine crystals at the real value of the deforming stress are studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2008

Uncertainty Analysis of Air Radiation for Lunar Return Shock Layers

By leveraging a new uncertainty markup technique, two risk analysis methods are used to compute the uncertainty of lunar-return shock layer radiation predicted by the High temperature Aerothermodynamic Radiation Algorithm (HARA). The effects of epistemic uncertainty, or uncertainty due to a lack of knowledge, is considered for the following modeling parameters: atomic line oscillator strengths, atomic line Stark broadening widths, atomic photoionization cross sections, negative ion photodetachment cross sections, molecular bands oscillator strengths, and electron impact excitation rates. First, a simplified shock layer problem consisting of two constant-property equilibrium layers is considered. The results of this simplified problem show that the atomic nitrogen oscillator strengths and Stark broadening widths in both the vacuum ultraviolet and infrared spectral regions, along with the negative ion continuum, are the dominant uncertainty contributors. Next, three variable property stagnation-line shock layer cases are analyzed: a typical lunar return case and two Fire II cases. For the near-equilibrium lunar return and Fire 1643-second cases, the resulting uncertainties are very similar to the simplified case. Conversely, the relatively nonequilibrium 1636-second case shows significantly larger influence from electron impact excitation rates of both atoms and molecules. For all cases, the total uncertainty in radiative heat flux to the wall due to epistemic uncertainty in modeling parameters is 30% as opposed to the erroneously-small uncertainty levels (plus or minus 6%) found when treating model parameter uncertainties as aleatory (due to chance) instead of epistemic (due to lack of knowledge).

Kleb, Bil↗
DOE OSTI2015

United polarizable multipole water model for molecular mechanics simulation

We report the development of a united AMOEBA (uAMOEBA) polarizable water model, which is computationally 3–5 times more efficient than the three-site AMOEBA03 model in molecular dynamics simulations while providing comparable accuracy for gas-phase and liquid properties. In this coarse-grained polarizable water model, both electrostatic (permanent and induced) and van der Waals representations have been reduced to a single site located at the oxygen atom. The permanent charge distribution is described via the molecular dipole and quadrupole moments and the many-body polarization via an isotropic molecular polarizability, all located at the oxygen center. Similarly, a single van der Waals interaction site is used for each water molecule. Hydrogen atoms are retained only for the purpose of defining local frames for the molecular multipole moments and intramolecular vibrational modes. The parameters have been derived based on a combination of ab initio quantum mechanical and experimental data set containing gas-phase cluster structures and energies, and liquid thermodynamic properties. For validation, additional properties including dimer interaction energy, liquid structures, self-diffusion coefficient, and shear viscosity have been evaluated. The results demonstrate good transferability from the gas to the liquid phase over a wide range of temperatures, and from nonpolar to polar environments, due to the presence of molecular polarizability. The water coordination, hydrogen-bonding structure, and dynamic properties given by uAMOEBA are similar to those derived from the all-atom AMOEBA03 model and experiments. Thus, the current model is an accurate and efficient alternative for modeling water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1973

Rates and mechanisms of the atomic oxygen reaction with nickel at elevated temperatures

The oxidation of nickel by atomic oxygen at pressure from 1 to 45 N/sq m between 1050 and 1250 K was investigated. In these ranges, the oxidation was found to follow the parobolic rate law, viz., K sub p = 0.0000114 exp(-13410/T) g squared/cm4/sec for films of greater than 1 micron thickness and was pressure independent. The activation enthalpy for the oxidation reaction was 112 + or - 11 kj/mole (27 + or - 3 kcal/mole). Of a number of possible mechanisms and defect structures considered, it was shown that the most likely was a saturated surface defect model for atomic oxidation, based on reaction activation enthalpies, impurity effects, pressure independence, and magnitudes of rates. A model judged somewhat less likely was one having doubly ionized cationic defects rate controlling in both atomic and molecular oxygen. From comparisons of the appropriate processes, the following enthalpy values were derived: enthalpy of activation (Ni diffusion in Ni0) = 110 + or - 30 kj/mole and standard enthalpy change for reaction formation (doubly ionized cation vacancies in Ni0 from atomic oxygen)= -9 + or - 25 kj/mole.

Christian, J. D.↗
DOE OSTI2017

The Influence of Atomic Diffusion on Stellar Ages and Chemical Tagging

In the era of large stellar spectroscopic surveys, there is an emphasis on deriving not only stellar abundances but also the ages for millions of stars. In the context of Galactic archeology, stellar ages provide a direct probe of the formation history of the Galaxy. We use the stellar evolution code MESA to compute models with atomic diffusion—with and without radiative acceleration—and extra mixing in the surface layers. The extra mixing consists of both density-dependent turbulent mixing and envelope overshoot mixing. Based on these models we argue that it is important to distinguish between initial, bulk abundances (parameters) and current, surface abundances (variables) in the analysis of individual stellar ages. In stars that maintain radiative regions on evolutionary timescales, atomic diffusion modifies the surface abundances. We show that when initial, bulk metallicity is equated with current, surface metallicity in isochrone age analysis, the resulting stellar ages can be systematically overestimated by up to 20%. The change of surface abundances with evolutionary phase also complicates chemical tagging, which is the concept that dispersed star clusters can be identified through unique, high-dimensional chemical signatures. Stars from the same cluster, but in different evolutionary phases, will show different surface abundances. We speculate that calibration of stellar models may allow us to estimate not only stellar ages but also initial abundances for individual stars. In the meantime, analyzing the chemical properties of stars in similar evolutionary phases is essential to minimize the effects of atomic diffusion in the context of chemical tagging.

79 ASTRONOMY AND ASTROPHYSICS↗