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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Mapping the body’s battle with Ebola and Zika

The viruses that cause Ebola and Zika, daunting diseases that inspire concern at every outbreak, share a strong similarity in how they first infiltrate a host’s cells. Through the computer modeling capabilities of Los Alamos National Laboratory, the molecular calisthenics involved in invading a cell are visually documented, an essential step toward vaccine and therapeutic medicine development. “The two very different viruses both use a surface protein to initiate a process called ‘membrane fusion,’ which allows the viral entry into the host cell,” said Chang-Shung Tung, the leading and corresponding author of the recent work in the journal Biomolecules. “Before the membrane fusion started, these two proteins went through large conformational changes that was not fully understood, as the proteins would take on the physical shape that would allow them to fuse with the cell and begin the infection. In this work, we modeled the structures of the two proteins in the fusion-initiation state. Then we used the structures to explain how four different antibodies in our bodies can fight the infection and block the viral entry, based on how they bind to these two proteins,” Tung said. The aim of the team’s work was to understand the specific structure-function relationship of the Ebola glycoprotein (EBOV GP) and Zika envelope (ZIKA E) proteins, the ones that are instrumental in the fusion with the host cell. “This computer modeling approach allowed us to generate atomic-level models of EBOV GP and ZIKV E proteins in the fusion state for the first time,” said Anna Lappala, a postdoctoral fellow at the Center for Nonlinear Studies and the Theoretical Biology & Biophysics Group at Los Alamos. “These models, together with the knowledge of the envelope proteins in the pre-fusion state, allowed us to study the function of neutralizing antibodies as they bind to the viral envelope proteins.” Read the published scientific paper here: http://www.mdpi.com/2218-273X/8/2/25/pdf

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2004

Studies of Pressure-Broadening of Alkali Atom Resonance Lines for Modeling Atmospheres of Extrasolar Giant Planets and Brown Dwarfs

In L-dwarfs and T-dwarfs the resonance lines of sodium and potassium are so profoundly pressure-broadened that their wings extend several hundred nanometers from line center. With accurate knowledge of the line profiles as a function of temperature and pressure: such lines can prove to be valuable diagnostics of the atmospheres of such objects. We have initiated a joint program of theoretical and experimental research to study the line-broadening of alkali atom resonance lines due to collisions with species such as helium and molecular hydrogen. Although potassium and sodium are the alkali species of most interest in the atmospheres of cool brown dwarfs and extrasolar giant planets, some of our theoretical focus this year has involved the calculation of pressure-broadening of lithium resonance lines by He, as a test of a newly developed suite of computer codes. In addition, theoretical calculations have been carried out to determine the leading long range van der Waals coefficients for the interactions of ground and excited alkali metal atoms with helium atoms, to within a probable error of 2%. Such data is important in determining the behavior of the resonance line profiles in the far wings. Important progress has been made on the experimental aspects of the program since the arrival of a postdoctoral fellow in September. A new absorption cell has been designed, which incorporates a number of technical improvements over the previous cell, including a larger cell diameter to enhance the signal, and fittings which allow for easier cleaning, thereby significantly reducing the instrument down-time.

Kirby, Kate↗
DOE OSTI2022

Model development for flash boiling spray and validations with isooctane, hexane, ethanol and their binary mixtures

The objective of this study is to develop a comprehensive model for spray atomization under the consideration of flash boiling, which includes nucleation, bubble growth, and droplet breakup. Although flash boiling has been previously studied, there is still no well-established analytical and numerical model of their mechanisms, or general guidelines on the implementation of this atomization strategy. In this study, detailed formulations are elaborated for the unified spray atomization model under both non-flash and flash boiling conditions. Next, the model was implemented into CFD software and 3-D simulations were performed for model validations. Two fuel blends are tested: ethanol-isooctane blends and hexane-isooctane blends. Simulation results were compared with the experimental measurements conducted using two methodologies, DBI (Diffuse Backlight Imaging) system for spray images captured by a high-speed camera, and a PDA (Phase Doppler Anemometry) system for droplet size measurement during the spray injection process. It was found that the current flash boiling model can accurately capture the spray characteristics under different environment conditions, such as the spray plume interactions and liquid penetration. The averaged Sauter mean diameter (SMD) calculations match the measurement sufficiently well. The spray structures for both single- and multi-component fuels during the injection process are well represented by current model.

42 ENGINEERING↗
DOE OSTI2013

Study on glass-forming ability and hydrogen storage properties of amorphous Mg{sub 60}Ni{sub 30}La{sub 10−x}Co{sub x} (x = 0, 4) alloys

Mg{sub 60}Ni{sub 30}La{sub 10−x}Co{sub x} (x = 0, 4) amorphous alloys were prepared by rapid solidification, using a melt-spinning technique. X-ray diffraction and differential scanning calorimetry analysis were employed to measure their microstructure, thermal stability and glass-forming ability, and hydrogen storage properties were studied by means of PCTPro2000. Based on differential scanning calorimetry results, their glass-forming ability and thermal stability were investigated by Kissinger method, Lasocka curves and atomic cluster model, respectively. The results indicate that glass-forming ability, thermal properties and hydrogen storage properties in the Mg-rich corner of Mg–Ni–La–Co system alloys were enhanced by Co substitution for La. It can be found that the smaller activation energy (ΔΕ) and frequency factor (υ{sub 0}), the bigger value of B (glass transition point in Lasocka curves), and higher glass-forming ability of Mg–Ni–La–Co alloys would be followed. In addition, atomic structure parameter (λ), deduced from atomic cluster model is valuable in the design of Mg–Ni–La–Co system alloys with good glass-forming ability. With an increase of Co content from 0 to 4, the hydrogen desorption capacity within 4000 s rises from 2.25 to 2.85 wt.% at 573 K. - Highlights: • Amorphous Mg{sub 60}Ni{sub 30}La{sub 10−x}Co{sub x} (x = 0 and 4) alloys were produced by melt spinning. • The GFA and hydrogen storage properties were enhanced by Co substitution for La. • With an increase of Co content, the hydrogen desorption capacity rises at 573 K.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Single Atomic Iron Catalysts for Oxygen Reduction in Acidic Media: Particle Size Control and Thermal Activation

It remains a grand challenge to replace platinum group metal (PGM) catalysts with earth-abundant materials for the oxygen reduction reaction (ORR) in acidic media, which is crucial for large-scale deployment of proton exchange membrane fuel cells (PEMFCs). Here, we report a high-performance atomic Fe catalyst derived from chemically Fe-doped zeolitic imidazolate frameworks (ZIFs) by directly bonding Fe ions to imidazolate ligands within 3D frameworks. Although the ZIF was identified as a promising precursor, the new synthetic chemistry enables the creation of well-dispersed atomic Fe sites embedded into porous carbon without the formation of aggregates. The size of catalyst particles is tunable through synthesizing Fe-doped ZIF nanocrystal precursors in a wide range from 20 to 1000 nm followed by one-step thermal activation. Similar to Pt nanoparticles, the unique size control without altering chemical properties afforded by this approach is able to increase the number of PGM-free active sites. The best ORR activity is measured with the catalyst at a size of 50 nm. Further size reduction to 20 nm leads to significant particle agglomeration, thus decreasing the activity. Using the homogeneous atomic Fe model catalysts, we elucidated the active site formation process through correlating measured ORR activity with the change of chemical bonds in precursors during thermal activation up to 1100 °C. The critical temperature to form active sites is 800 °C, which is associated with a new Fe species with a reduced oxidation number (from Fe3+ to Fe2+) likely bonded with pyridinic N (FeN4) embedded into the carbon planes. Further increasing the temperature leads to continuously enhanced activity, linked to the rise of graphitic N and Fe–N species. The new atomic Fe catalyst has achieved respectable ORR activity in challenging acidic media (0.5 M H2SO4), showing a half-wave potential of 0.85 V vs RHE and leaving only a 30 mV gap with Pt/C (60 μgPt/cm2). Enhanced stability is attained with the same catalyst, which loses only 20 mV after 10 000 potential cycles (0.6–1.0 V) in O2 saturated acid. The high-performance atomic Fe PGM-free catalyst holds great promise as a replacement for Pt in future PEMFCs.

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

Monte Carlo Computational Modeling of the Energy Dependence of Atomic Oxygen Undercutting of Protected Polymers

A Monte Carlo computational model has been developed which simulates atomic oxygen attack of protected polymers at defect sites in the protective coatings. The parameters defining how atomic oxygen interacts with polymers and protective coatings as well as the scattering processes which occur have been optimized to replicate experimental results observed from protected polyimide Kapton on the Long Duration Exposure Facility (LDEF) mission. Computational prediction of atomic oxygen undercutting at defect sites in protective coatings for various arrival energies was investigated. The atomic oxygen undercutting energy dependence predictions enable one to predict mass loss that would occur in low Earth orbit, based on lower energy ground laboratory atomic oxygen beam systems. Results of computational model prediction of undercut cavity size as a function of energy and defect size will be presented to provide insight into expected in-space mass loss of protected polymers with protective coating defects based on lower energy ground laboratory testing.

Banks, Bruce A.↗
DOE OSTI2015

Development of high-fidelity models for liquid fuel spray atomization and mixing processes in transportation and energy systems

Significant inadequacies of current models for multiphase flows present a major barrier to rapid development of advanced high-efficiency low-emissions combustion devices. Liquid spray atomization processes largely determine fuel-air mixture formation, which subsequently govern combustion and controls performance, emissions, and durability of a device. The current study presents a fundamentally-consistent framework to model the effects of breakup processes, liquid drop deformations, and internal flow dynamics on mass, momentum, and energy exchange functions. This framework builds on the Taylor Analogy Breakup (TAB) model which naturally quantifies local drop deformation dynamics. Real-fluid multicomponent thermodynamic property modeling and Gradient Theory simulations facilitate accurate calculations of molecular two-phase interface exchange functions, surface tensions forces, drop oscillations, and breakup processes. The analysis establishes that these drop dynamics, along with finite liquid viscosity effects, significantly alter gas-liquid exchange functions. The study quantifies these effects for the resulting drag coefficients of liquid drops and demonstrates significant deviations from the classic dynamic drag model, which is widely applied in modern simulations performed in academia and industry. This work also quantifies effects, which originate from gas-liquid coupling dynamics, on evaporation and heating rates. The analysis establishes that the consideration of these coupling dynamics modify mass and energy transfer rates even more significantly than the corresponding drag forces from momentum exchange. This physical complexity, however, is largely neglected in modern studies from academia and industry. A new set of equations is presented to improve the modeling of drop breakup processes to address the current shortcomings in the prediction of resulting drops proper- ties over the full range of relevant breakup conditions. The framework is based on a refined energy balance equation which explicitly enforces drop momentum conservation during the breakup process. The introduced modeling framework is entirely derived from conservation equations for mass, momentum, and energy and does not, as a consequence, introduce new modeling constants. The significance of the developed modeling advances to fuel injection processes is demonstrated using Large Eddy Simulation (LES) with a Lagrangian-Eulerian modeling approach.

42 ENGINEERING↗
DOE OSTI2026

Towards Spatial Control of Reaction Selectivity on Photocatalysts using Area Selective Atomic Layer Deposition on Model Dual Site Electrocatalyst Platform

Photocatalytic water splitting is a promising route to low-cost H2; however, this approach is currently limited by solar-to-hydrogen (STH) conversion efficiencies in the sub-10% range. Z-scheme water splitting, in which H2 and O2 evolving particles are operated in separate compartments and electronically coupled by a soluble redox mediator, offers improvements in STH efficiency but introduces high rates of undesired side and back reactions. Nanoscopic oxide coatings (e.g., CrOx, and SiOx) have previously been used to selectively block undesired reactants from reaching active sites; however, a coating encapsulating the entire photocatalyst particle limits activity as the particle can no longer facilitate both half reactions. In contrast to photodeposition, which may produce non-uniform overlayers due to the difficulty of controlling local electrochemical reactivity on the surface, area-selective atomic layer deposition (AS-ALD) can be used to deposit conformal ultrathin oxide coatings while preserving access to non-growth areas of a substrate surface. To develop this technique for Z-scheme photocatalysts, we performed AS-ALD of TiO2 on a dual site planar electrocatalyst based on interdigitated arrays of Pt and Au. Self-assembled monolayers of 1-octadecanethiol (ODT) were used to block growth on the Au array, resulting in a patterned surface in which only the Pt sites were encapsulated with TiO2. These electrocatalysts demonstrated localized reaction selectivity for hydrogen evolution reaction (HER) over the TiO2/Pt sites, while the uncoated Au sites retained activity towards Fe(II)/Fe(III) redox (FeRR). Under independent potential control, these microelectrodes showed that selectively deposited TiO2 coatings can suppress the rate of back reactions on neighboring active sites by an order of magnitude compared to uncoated control samples.

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

In-situ Study of Dynamic Phenomena at Metal Nanosolder Interfaces Using Aberration Corrected Scanning Transmission Electron Microcopy.

Controlling metallic nanoparticle (NP) interactions plays a vital role in the development of new joining techniques (nanosolder) that bond at lower processing temperatures but remain viable at higher temperatures. The pr imary objective of this project is t o develop a fundamental understanding of the actual reaction processes, associated atomic mechanisms, and the resulting microstructure that occur during thermally - driven bond formation concerning metal - metal nano - scale (%3C50nm) interfaces. In this LDRD pr oject, we have studied metallic NPs interaction at the elevated temperatures by combining in - situ transmission electron microscopy (TEM ) using an aberration - corrected scanning transmission electron microscope (AC - STEM) and atomic - scale modeling such as m olecular dynamic (MD) simulations. Various metallic NPs such as Ag, Cu and Au are synthesized by chemical routines. Numerous in - situ e xperiments were carried out with focus of the research on study of Ag - Cu system. For the first time, using in - situ STEM he ating experiments , we directly observed t he formation of a 3 - dimensional (3 - D) epitaxial Cu - Ag core - shell nanoparticle during the thermal interaction of Cu and Ag NPs at elevated temperatures (150 - 300 o C). The reaction takes place at temperatures as low as 150 o C and was only observed when care was taken to circumvent the effects of electron beam irradiation during STEM imaging. Atomic - scale modeling verified that the Cu - Ag core - shell structure is energetically favored, and indicated that this phenomenon is a nano - scale effect related to the large surface - to - volume ratio of the NPs. The observation potentially can be used for developing new nanosolder technology that uses Ag shell as the "glue" that stic ks the particles of Cu together. The LDRD has led to several journal publications and numerous conference presentations, and a TA. In addition, we have developed new TEM characterization techniques and phase - field modeling tools that can be used for future materials research at Sandia. Acknowledgeme nts This work was supported by the Laboratory Directed Research and Development (LDRD) program of Sandia National Laboratories. Sandia National Laboratories is a multi - program laboratory managed and operated by Sandia Corporation, a wholly owned subsidia ry of Lockheed Martin Corporation, for the US Department of Energy's National Nuclear Security Administration under contract DE - AC04 - 94AL85000.

Lu, Ping↗
DOE OSTI2018

Solving many-body Schrödinger equations with kinetic energy partition method

We present a general formulation of our previously developed kinetic energy partition (KEP) method for solving many-bodySchrödinger equations. In atomic physics, as well as in general molecular and solid state physics, solving many-electronSchrödinger equations is a very challenging task, often called Dirac’s challenge. The central problem is how to properly handle the electron–electron Coulomb repulsion interactions. Using the KEP solution scheme, in addition to dividing the kinetic energy into partial terms, the electron–electron Coulomb interaction is also separated into parts to be associated with a “negative mass” kinetic energy term. Therefore, the full Hamiltonian can be expressed as a simple sum of subsystem Hamiltonians, each representing an effective one-body problem. Using a Hartree-like product in constructing the wave-function, we achieve fast convergence in the calculations of the ground state energies. First, the model Moshinsky atoms are used to illustrate the solution procedure. We then apply this new KEP method to harmonium atoms and obtain precise energies with an error less than 5% using only two basis functions from each subsystem. It is thus very promising that this methodology, when further extended, can be useful for general many-body systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

The use of plasma ashers and Monte Carlo modeling for the projection of atomic oxygen durability of protected polymers in low Earth orbit

The results of ground laboratory and in-space exposure of polymeric materials to atomic oxygen has enabled the development of a Monte Carlo computational model which simulates the oxidation processes of both environments. The cost effective projection of long-term low-Earth-orbital durability of protected polymeric materials such as SiO(x)-coated polyimide Kapton photovoltaic array blankets will require ground-based testing to assure power system reliability. Although silicon dioxide thin film protective coatings can greatly extend the useful life of polymeric materials in ground-based testing, the projection of in-space durability based on these results can be made more reliable through the use of modeling which simulates the mechanistic properties of atomic oxygen interaction, and replicates test results in both environments. Techniques to project long-term performance of protected materials, such as the Space Station Freedom solar array blankets, are developed based on ground laboratory experiments, in-space experiments, and computational modeling.

Banks, Bruce A.↗
DOE OSTI2017

Coupled Cluster and Density Functional Studies of Atomic Fluorine Chemisorption on Coronene as Model Systems for Graphene Fluorination

Potential energy curves (PECs) of atomic fluorine adsorption on coronene as a model for graphene or nanocarbon surfaces have been computed. The PECs were obtained by scanning the fluorine atom distance to one of the center carbon atoms of the coronene molecule as model system in a “top” position from 4.0 to1.0 Å in intervals of 0.1 Å using a variety of quantum chemical methods. Various density functional theory (DFT) functionals, such as B3LYP, PBE, PBE0, CAM-B3LYP, and LC-ωPBE; approximate DFT methods such as several levels of the density-functional tight-binding (DFTB) method, as well as ab initio wave function theory methods such as MP2, CCSD, CCSD(T), and G2MS extrapolations, were used to evaluate energies for B3LYP/cc-pVDZ PEC geometries. G2MS is an approximation to the highly accurate CCSD(T)/cc-pVTZ level of theory in our work. We found that fluorine is chemically adsorbed on coronene with a binding energy of 22.9 and 23.3 kcal/mol at the B3LYP/cc-pVDZ and G2MS levels of theory, respectively, and 18.3 and 19.3 kcal/mol after counterpoise correction. Here, we found that pure DFT functionals and their DFTB approximations fail to predict the correct dissociation limit due to the DFT-inherent self-interaction error and various limitations in the DFTB approximation itself.

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

Features of the microstructure development under conditions, reproducing the process of friction stir welding. Molecular-dynamics study

Friction stir welding is a recently developed technology which is used in various branches of modern engineering. The basis of this technology is the friction of the rotating cylindrical or specially shaped tool between two metal plates brought together either to meet their ends of one above another with the overlap. When applying the FSW process in various economical sectors, the important task is to study the mechanisms and identify the physical laws and factors leading to formation of structural inhomogeneities and discontinuities in the weld seam. This paper analyzes the basic mechanisms behind the structural state generation in the material subjected to severe plastic deformation and heating. To investigate the atomic mechanisms of structural changes in FSW, the modeling at atomic scale has been carried out. Results of work can be a basis for new knowledge about the microstructure evolution in FSW.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗