Search NASA⌕ Search

SEARCH · Search NASA

Results for “Atomistic simulations”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Theories of homogeneous and electrochemical electron transfer in complex media and interfaces (Final Technical Report)

This project makes the next step in establishing practical theories of charge transfer in complex media. The development of formal models is supported by extensive atomistic simulations, quantum calculations of force-field parameters, and direct measurements of charge-transfer spectra. All theory development is supported by experiment, extensive numerical simulations, and through external collaborations.

14 SOLAR ENERGY↗
DOE OSTI2024

Ultimate compressive strength and severe plastic deformation of equilibrated single-crystalline copper nanoparticles

Mechanical properties and deformation mechanisms of defect-free copper nanoparticles are investigated by combining experiments with atomistic simulations. The compressive strength of the particles increases with decreasing size and tends to saturate near the theoretical strength in the small-size limit. In this limit, the intrinsic size dependence of the strength is governed by the stochastic nature of dislocation nucleation near the particle surface. The particle deformation process evolves from the initial strain softening to strain hardening as the particle accumulates residual damage. The normalized strength-size relation for Cu is compared with those for Au, Ni, and Pt. The lack of universal behavior among the four FCC metals is discussed. Heavily deformed Cu nanoparticles develop polycrystalline structures and change the lattice orientation from [111] to [110]. The experiments and simulations reveal the twinning mechanism of the lattice rotation leading to the new grain formation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Optimizing the Heisenberg Vortex Tube for Hydrogen Cooling (Final Technical Report)

Hydrogen utilization at Plug Power sites ranges between 70-93% with the remainder vented and lost to the atmosphere. The goal of this project was to improve hydrogen utilization at Plug Power fulfillment centers via the patented Heisenberg vortex tube (HVT). The HVT combines conventional counterflow vortex tubes with para-orthohydrogen conversion catalyst to allow the cooling of hydrogen flows with no moving parts or external energy inputs. In year 1 of this project the HVT was evaluated to improve hydrogen utilization via the following concepts: 1) increasing liquid hydrogen pump volumetric efficiency by 20% through vapor separation and subcooling of the liquid, 2) reducing liquid hydrogen storage tank boil-off losses by 20% through thermal vapor shielding (TVS), and 3) increase isentropic efficiency of supercritical hydrogen expansion at 40-50 K, from 31% to over 40%. The TVS concept was selected as the most favorable application for further development. While the conceptual analysis was underway, the cryocatalysis hydrogen experiment facility (CHEF) was retrofitted with a new cryocooler, higher pressure condenser tanks, and in-situ fiber-optic Raman probes for ortho-parahydrogen composition analysis at the inlets and outlets of the HVT. To our knowledge this is the first in-situ implementation of cryogenic Raman probes for ortho-parahydrogen analysis. Subsequent testing of catalyzed and non-catalyzed HVT identified a low flow rate, high-conversion efficacy regime well suited for TVS development. Computational fluid dynamics (CFD), Reduced order modeling (ROM), and quantum Monte Carlo atomistic simulations were applied to optimize the design of the HVT implemented for TVS. The simulations and analysis identified a ruthenium-based catalyst as the most optimal for high conversion with little pressure losses and was matched to experimental measurements. These results indicated the design could achieve the 20% boil-off reduction target. An HVT field-trial was designed and constructed by Plug Power and implemented at a customer site. However, budget limitations reduced the amount of catalyst that could be applied to the HVT. With the reduced amount of catalyst, CFD analysis estimated about 3% reduction in boil-off. A similar amount of boil-off reduction was observed during preliminary measurements during tank commissioning. However, the results have a large margin for error and were operated at relatively low tank temperatures when the HVT has the lowest potential for improvement. This project demonstrated the use of para-orthohydrogen catalysis for reducing liquid hydrogen tank boil-off losses. The field trial system is anticipated to provide long-term experimental measurements on tank operational performance that will allow Plug Power to instrument additional tanks in the future. Recommendations for future work include the development of higher activity, lower cost para-orthohydrogen conversion catalyst for use in the HVT. Tank process optimizations could maximize the efficacy of the HVT and further improve hydrogen utilization.

08 HYDROGEN↗
DOE OSTI2025

Multimetallic Layered Composites (MMLCs) for Rapid, Economical Advanced Reactor Deployment (Final Report)

This project focused on the development of multi-metallic layered composites (MMLCs) for advanced fission reactor technologies. There are many instances where one alloy or material simply cannot meet all the demands thrown at it by a reactor system, or cannot allow it to perform as strongly as one would like. Instead of focusing all our effort on developing one perfect alloy, we seek to leverage the design principle of “separation of functionality,” used in many other arenas in design, to boost performance beyond single alloys alone. One illustrative example shows the power of this approach for molten salt-cooled reactors: A three meter tall, three meter diameter reactor vessel made of Incoloy 800 was quoted at $\$$500k in 2018. A Hastelloy N vessel was quoted at $\$$5M. An MMLC vessel, in which a layer of Hastelloy N would be weld-overlaid onto Incoloy 800, was quoted at $\$$700k, and it would achieve the same performance. The potential economic gains of leveraging this approach are therefore substantial. At a minimum, each MMLC would contain one core structural layer and one coolant-facing corrosion-resistant layer. Sometimes, MMLCs required buffer layers, as the structural and corrosion-resistant layers were metallurgically incompatible. In other words, they didn’t always play nice, thus separating layers compatible with both functioned as intermediaries to keep the composite together. However, in doing so we inevitably produce new interfaces, where new issues can arise. Therefore, this project focused on what happens at these interfaces from a combination of high temperatures, irradiation, corrosion, and time. After all, a reactor makes money when it is operating, and outages of any kind erode its economic viability. First, we set out to experimentally prove that MMLCs for at least two advanced reactor systems can be made, today, in US domestic facilities. In this respect we were successful – one MMLC (a Ni-201/Incoloy 800H composite) was successfully made and drawn into two-inch coolant piping. Others were attempted, though new issues relating to cracking in vanadium layers for one and radiation damage performance of the corrosion-resistant layer in another prevented us from moving further in those specific arenas – these are engineering problems which deserve continued focus after this project. Additional experimental work focused on long-term corrosion testing of the outermost layers of the salt-cooled and liquid lead-cooled MMLC concepts, which would then be fed into predictions of how long the MMLCs could last. Next, computational (thermodynamics and atomistic) simulation studies studied how much we expect the interfaces to “blend,” due to the mixing action of neutron irradiation. This eats into both the margin for the structural layer of each MMLC, as dilution from the corrosion-resistant layer into the structural layer would decrease the total load-bearing capacity of an MMLC of finite size. On the other hand, dilution of the corrosion-resistant layer into the structural layer further reduced the margin of corrodible material, reducing the lifetime of the MMLC or necessitating extra thickness to be imparted to the MMLC to meet its functional requirements. Work here focused on irradiation-induced segregation to predict new phases which may embrittle the MMLCs, as well as quantifying irradiation-induced mixing at each interface. The results showed that mixing is expected, but it is both steady and therefore predictable, and not lifetime-limiting for most MMLC concepts – it simply has to be accounted for in calculations of reactor performance when utilizing an MMLC. Then, full-core simulations using the experimentally-derived corrosion data, the computationally discovered irradiation-induced mixing data (partially validated by experiment), and existing, benchmarked core designs for large and small sized reactor concepts (one salt-cooled, one lead-cooled) were conducted to quantify any expansion of reactor operating envelopes achieved by utilizing these MMLCs. This new framework, called REX (Reactor Envelope Expansion), incorporates a combination of core neutronics, thermal hydraulics, and the material performance data derived from this project to see how using an MMLC expands advanced fission reactor operating envelopes. It was discovered that in some cases, MMLC utilization does indeed increase the maximum operating temperatures and cycle lengths of reactor concepts, while in other cases it does not. Finally, our tech-to-market (T2M) strategy was not necessarily to create specific embodiments of MMLCs for immediate sale (because getting into the nuclear market is incredibly slow and laden with regulation, this is a long-term goal), but rather immediate stimulation of US industry using the design approach of MMLCs derived from this project. In this respect we were successful, as one of the PhD students funded on this project co-founded Allium Engineering, Inc., which created a stainless steel / low-alloy steel MMLC to function as chloride corrosion-resistant rebar for embedding into concrete structures. Allium Engineering continues to be successful, having recently opened their first factory as of this writing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Tailoring composition and deformation modes at the microstructural level for next generation low-cost high-strength austenitic stainless steels

The objective of this project is to enable deliberate development of cost-effective, hydrogen resistant alloys by establishing detailed relationships specific to the effects of alloy composition, short-range order (SRO), and microsegregation in the presence of hydrogen on the transition between homogeneous deformation and localized plasticity in shear bands. In collaboration with the International Institute for Carbon-Neutral Energy Research, I2CNER, at Kyushu University in Japan, we conceptualized, designed, and manufactured four austenitic alloys that maintain corrosion resistance and ensure lower cost relative to baseline commercial alloys. The mechanical properties and deformation modes of the novel alloys (KU alloys) were assessed in the presence of hydrogen (H). Correlations between composition and performance revealed that two of the KU alloys are suitable replacements for 316 steel, while another is a viable replacement for 304 steel at room temperature. We found that, in the presence of other austenite stabilizing elements namely Mn and N, replacing Ni with Cu does not lead to martensite formation as has been previously reported.1–3 Furthermore, we found that the addition of Cu leads to an earlier onset of multiple slip resulting in an relative earlier onset of a higher work hardening rate (WHR). Greater understanding of the relationships between alloy composition and SRO required the development of a novel advanced electron diffraction methodology to characterize SRO in complex FCC alloys. This innovative approach, which combines fluctuation and correlation analyses of diffuse-scattering signals, successfully differentiated between SRO and long-range ordering (LRO). Further investigations into annealed austenitic stainless steels could provide insights into manipulating SRO and its effects on material properties. Atomistic simulations provided understanding of SRO behavior that was difficult to capture experimentally. This project created the first spin cluster expansion model that is able to capture and describe SRO effects in Fe-Ni-Cr FCC alloys, accounting for the non-negligible effects of magnetism. An automated computational workflow was established to provide reliable predictions of SRO in Fe-Ni-Cr austenitic alloys, both with and without the presence of H atoms. Analysis of the propensity for SRO in Fe-Ni-Cr alloys revealed that H tends to cluster with specific, well-defined SRO domains. The computational framework is general purpose and can be extended to realistic stainless steels across diverse composition ranges. With confidence that SRO is possible in austenitic stainless steels, we developed a discrete dislocation finite element code to understand the interaction of dislocations with SRO in the presence of H. By incorporating H effects on the dislocation emission and SRO stress field we show that the critical stress for the dislocation pileup to breakthrough the SRO domain decreases in the presence of H, which directly contributes localized deformation at the macroscale. Through the simulation of a uniaxial tension test, we demonstrated that H-induced weakening of SRO stress field and H-enhanced dislocation emission can lead to the onset of shear localization at lower macroscopic strains. As a whole, this project identified three novel alloys that show improvements in performance and cost efficiency for H-facing applications by studying correlations between alloy chemistry and deformation behavior. We also made significant advancements to experimental and computational methodologies necessary to study the chemistry and distribution of SRO across a range of alloys, which in turn allowed us to demonstrate how deformation mechanisms change due to the contributions of SRO in austenitic alloys in the presence of H. The combined advancements in fundamental understanding with novel alloy development in this project has increased the viability of next generation H-technologies for the broader public through accessible low-cost alloys and accelerated development towards future H-infrastructure.

08 HYDROGEN↗
DOE OSTI2024

Towards engineering the growth morphology of anisotropic bicrystals

Eutectic alloys, such as Pb-free Bi-Sn-Ag or Sn-Ag-Cu, are essential for their low melting points in solder joints but suffer from brittle intermetallic phases, reducing reliability. Traditional solder design relies on equilibrium phase diagrams, often neglecting the kinetic influence of crystalline anisotropy during solidification. This study hypothesizes that the energetics at the solid-solid-liquid trijunction dictate crystalline facet selection. Using the Al-Al 3 Ni alloy as a model, this project combines atomistic simulations and experimental methods to understand interfacial energies and develop predictive models. This approach aims to mitigate the negative effects of intermetallics, enhancing component reliability in electronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Rigorous incorporation of pH effects into ab initio electrochemical models

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0023410. The theme of this project was to update electronically grand canonical atomistic simulations to properly incorporate pH effects, on a rigorous thermodynamic basis, and to compare this to experimental findings. This project was primarily conducted by PI Andrew Peterson and his group at Brown University, with an unfunded collaboration with the Quantum Simulations Group at Lawrence Livermore National Laboratory with Dr. Joel Varley, who hosted the primary graduate student (Alexander Lackey) funded by this project and contributed to all major intellectual pursuits, utilizing the expertise and capabilities of the national laboratory to advance the project goals. Postdoctoral associate Dr. Juye Kim was also a major contributor to this project, and implementation of the final section included cooperation with several electronic-structure open-source software developers, including Colin Baker (Brown), Sandeep Sharma (CalTech), Georg Kastlunger (DTU), and Jens Jørgen Mortensen (DTU).

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

Electronic materials high-T(sub c) superconductivity polymers and composites structural materials surface science and catalysts industry participation

The fifth year of the Center for Advanced Materials was marked primarily by the significant scientific accomplishments of the research programs. The Electronics Materials program continued its work on the growth and characterization of gallium arsenide crystals, and the development of theories to understand the nature and distribution of defects in the crystals. The High Tc Superconductivity Program continued to make significant contributions to the field in theoretical and experimental work on both bulk materials and thin films and devices. The Ceramic Processing group developed a new technique for cladding YBCO superconductors for high current applications in work with the Electric Power Research Institute. The Polymers and Composites program published a number of important studies involving atomistic simulations of polymer surfaces with excellent correlations to experimental results. The new Enzymatic Synthesis of Materials project produced its first fluorinated polymers and successfully began engineering enzymes designed for materials synthesis. The structural Materials Program continued work on novel alloys, development of processing methods for advanced ceramics, and characterization of mechanical properties of these materials, including the newly documented characterization of cyclic fatigue crack propagation behavior in toughened ceramics. Finally, the Surface Science and Catalysis program made significant contributions to the understanding of microporous catalysts and the nature of surface structures and interface compounds.

Source record↗
NASA NTRS1990

Energetics and crystal chemistry of Ruddlesden-Popper type structures in high T(sub c) ceramic superconductors

The formation of Ruddlesden-Popper type layers (alternating slabs of rocksalt and perovskite structures) is seen in these oxides which is similar in many respects to what is seen in the system Sr-Ti-O. However, it was observed that there are some significant differences, for example the rocksalt and perovskite blocks in new superconducting compounds are not necessarily electrically, unlike in Sr-Ti-O systems. This will certainly render an additional coulombic bonding energy between two different types of blocks and may well lead to significant differences in their structural chemistry. In the higher order members of the various homologous series, additional Cu-O planes are inserted in the perovskite blocks. In order for the unit cell to electrically neutral the net positive charge on rocksalt block (which remains constant throughout the homologous series) should be balanced by an equal negative charge on perovskite block. It, thus becomes necessary to create oxygen vacancies in the basic perovskite structure, when width of the perovskite slab changes on addition of extra Cu-O planes. Results of atomistic simulations suggest that these missing oxygen ions allow the Cu-O planes to buckle in these compounds. This is also supported by the absence of buckling in the first member of Bi-containing compounds in which there are no missing oxygen ions and the Sr-Ti-O series of compounds. Additional results are presented on the phase stability of polytypoid structures in these crystal chemically complex systems. The studies will focus on the determination of the location of Cu(3+) in the structures of higher order members of the La-Cu-O system and whether Cu(3+) ions or oxygen vacancies are energetically more favorable charge compensating mechanism.

Dwivedi, Anurag↗
NASA NTRS1993

Diffusion on Cu surfaces

Understanding surface diffusion is essential in understanding surface phenomena, such as crystal growth, thin film growth, corrosion, physisorption, and chemisorption. Because of its importance, various experimental and theoretical efforts have been directed to understand this phenomena. The Field Ion Microscope (FIM) has been the major experimental tool for studying surface diffusion. FIM have been employed by various research groups to study surface diffusion of adatoms. Because of limitations of the FIM, such studies are only limited to a few surfaces: nickel, platinum, aluminum, iridium, tungsten, and rhodium. From the theoretical standpoint, various atomistic simulations are performed to study surface diffusion. In most of these calculations the Embedded Atom Method (EAM) along with the molecular static (MS) simulation are utilized. The EAM is a semi-empirical approach for modeling the interatomic interactions. The MS simulation is a technique for minimizing the total energy of a system of particles with respect to the positions of its particles. One of the objectives of this work is to develop the EAM functions for Cu and use them in conjunction with the molecular static (MS) simulation to study diffusion of a Cu atom on a perfect as well as stepped Cu(100) surfaces. This will provide a test of the validity of the EAM functions on Cu(100) surface and near the stepped environments. In particular, we construct a terrace-ledge-kink (TLK) model and calculate the migration energies of an atom on a terrace, near a ledge site, near a kink site, and going over a descending step. We have also calculated formation energies of an atom on the bare surface, a vacancy in the surface, a stepped surface, and a stepped-kink surface. Our results are compared with the available experimental and theoretical results.

Karimi, Majid↗
NASA NTRS2000

Influence of the Polysilicon Gate on the Random Dopant Induced Threshold Voltage Fluctuations in Sub 100 nm MOSFETS with Thin Gate Oxides

In this paper for the first time we study the influence of the polysilicon gate on the random dopant induced threshold voltage fluctuations in sub 100 nm MOSFETs with tunnelling gate oxides. This is done by using an efficient 3D 'atomistic' simulation technique described elsewhere. Devices with uniform channel doping and with low doped epitaxial channels have been investigated. The simulations reveale that the polysilicon gate is responsible for a substantial fraction of the threshold voltage fluctuations in both devices when the gate oxide is scaled to tunnelling thickness in the range of 1 - 2 nm.

Asenov, Asen↗
NASA NTRS2000

Polysilicon Gate Enhancement of the Random Dopant Induced Threshold Voltage Fluctuations in Sub-100 nm MOSFET's with Ultrathin Gate Oxide

In this paper, we investigate various aspects of the polysilicon gate influence on the random dopant induced threshold voltage fluctuations in sub-100 nm MOSFET's with ultrathin gate oxides. The study is done by using an efficient statistical three-dimensional (3-D) "atomistic" simulation technique described else-where. MOSFET's with uniform channel doping and with low doped epitaxial channels have been investigated. The simulations reveal that even in devices with a single crystal gate the gate depletion and the random dopants in it are responsible for a substantial fraction of the threshold voltage fluctuations when the gate oxide is scaled-in the range of 1-2 nm. Simulation experiments have been used in order to separate the enhancement in the threshold voltage fluctuations due to an effective increase in the oxide thickness associated with the gate depletion from the direct influence of the random dopants in the gate depletion layer. The results of the experiments show that the both factors contribute to the enhancement of the threshold voltage fluctuations, but the effective increase in the oxide-thickness has a dominant effect in the investigated range of devices. Simulations illustrating the effect or the polysilicon grain boundaries on the threshold voltage variation are also presented.

Asenov, Asen↗
NASA NTRS1999

Quantum Mechanical Enhancement of the Random Dopant Induced Threshold Voltage Fluctuations and Lowering in Sub 0.1 Micron MOSFETs

A detailed study of the influence of quantum effects in the inversion layer on the random dopant induced threshold voltage fluctuations and lowering in sub 0.1 micron MOSFETs has been performed. This has been achieved using a full 3D implementation of the density gradient (DG) formalism incorporated in our previously published 3D 'atomistic' simulation approach. This results in a consistent, fully 3D, quantum mechanical picture which implies not only the vertical inversion layer quantisation but also the lateral confinement effects manifested by current filamentation in the 'valleys' of the random potential fluctuations. We have shown that the net result of including quantum mechanical effects, while considering statistical fluctuations, is an increase in both threshold voltage fluctuations and lowering.

Asenov, Asen↗
NASA NTRS2001

Electronic Transport Through Carbon Nanotubes: Effects of Structural Deformation and the Tube Chirality

Atomistic simulations using a combination of classical force field and Density-Functional-Theory (DFT) show that carbon atoms remain essentially sp2 coordinated in either bent tubes or tubes pushed by an atomically sharp AFM tip. Subsequent Green's-function-based transport calculations reveal that for armchair tubes there is no significant drop in conductance, while for zigzag tubes the conductance can drop by several orders of magnitude in AFM-pushed tubes. The effect can be attributed to simple stretching of the tube under tip deformation, which opens up an energy gap at the Fermi surface.

Maiti, Amitesh↗
NASA NTRS2002

Surface Composition of NiPd Alloys

Surface segregation in Ni-Pd alloys has been studied using the BFS method for alloys. Not only does the method predict an oscillatory segregation profile but it also indicates that the number of Pd-enriched surface planes can vary as a function of orientation. The segregation profiles were computed as a function of temperature, crystal face, and composition. Pd enrichment of the first layer is observed in (111) and (100) surfaces, and enrichment of the top two layers occurs for (110) surfaces. In all cases, the segregation profile shows oscillations that are actually related to weak ordering tendencies in the bulk. An atom-by-atom analysis was performed to identify the competing mechanisms leading to the observed surface behaviors. Large-scale atomistic simulations were also performed to investigate the temperature dependence of the segregation profiles as well as for analysis of the bulk structures. Finally, the observed surface behaviors are discussed in relation to the bulk phase structure of Ni-Pd alloys, which exhibit a tendency to weakly order.

Noebe, Ronald D.↗
NASA NTRS2002

Atomistic Modeling of Quaternary Alloys: Ti and Cu in NiAl

The change in site preference in NiAl(Ti,Cu) alloys with concentration is examined experimentally via ALCHEMI and theoretically using the Bozzolo-Ferrante-Smith (BFS) method for alloys. Results for the site occupancy of Ti and Cu additions as a function of concentration are determined experimentally for five alloys. These results are reproduced with large-scale BFS-based Monte Carlo atomistic simulations. The original set of five alloys is extended to 25 concentrations, which are modeled by means of the BFS method for alloys, showing in more detail the compositional range over which major changes in behavior occur. A simple but powerful approach based on the definition of atomic local environments also is introduced to describe energetically the interactions between the various elements and therefore to explain the observed behavior.

Bozzolo, Guillermo↗
NASA NTRS2002

Atomistic Modeling of Co Growth on Cu(111]

The BFS method for alloys is applied to the study of Co growth on Cu(111). The parameterization of the Co-Cu system is obtained from first-principles calculations, and tested against known experimental features for low coverage Co deposition on Cu(100) and Cu(111). Atomistic simulations are performed to investigate the behavior of Co on Cu(111) as a function of coverage.

Khalil, Joseph↗