Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic model refinement”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Moiré metrology of energy landscapes in van der Waals heterostructures

The emerging field of twistronics, which harnesses the twist angle between two-dimensional materials, represents a promising route for the design of quantum materials, as the twist-angle-induced superlattices offer means to control topology and strong correlations. At the small twist limit, and particularly under strain, as atomic relaxation prevails, the emergent moiré superlattice encodes elusive insights into the local interlayer interaction. Here we introduce moiré metrology as a combined experiment-theory framework to probe the stacking energy landscape of bilayer structures at the 0.1 meV/atom scale, outperforming the gold-standard of quantum chemistry. Through studying the shapes of moiré domains with numerous nano-imaging techniques, and correlating with multi-scale modelling, we assess and refine first-principle models for the interlayer interaction. We document the prowess of moiré metrology for three representative twisted systems: bilayer graphene, double bilayer graphene and H-stacked MoSe 2 /WSe 2 . Moiré metrology establishes sought after experimental benchmarks for interlayer interaction, thus enabling accurate modelling of twisted multilayers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Systematic discovery of new nano-scale metastable intermetallic eutectic phases in laser rapid solidified Aluminum-Germanium alloy

Laser surface remelting of as-cast Al-Ge eutectic alloy is shown to produce ultrafine lamellar eutectic morphology with interlamellar spacing refined up to ∼60 nm and composed of FCC Al solid solution and unusual Al x Ge y intermetallic phases that do not form during near-equilibrium solidification. The microstructures are characterized and analyzed using a combination of selected area electron diffraction, high-resolution scanning transmission electron microscopy, energy dispersive X-ray spectroscopy to obtain high-resolution elemental maps, and atomistic modeling using density functional theory followed by atomic-scale image simulation. Depending on the local solidification conditions, the crystallography of the Al x Ge y intermetallic phases in the eutectic microstructure is either monoclinic (C 2/c) or monoclinic (P 2 1 ), with high densities of defects in both cases. This is in sharp contrast to the as-cast alloys that showed nominally pure Al and Ge phases with significant solute partitioning and equilibrium FCC and diamond cubic crystal structures, respectively. Corresponding kinetic phase diagrams are proposed to interpret the evolution of nano-lamellar eutectic morphologies with equilibrium Al and metastable Al x Ge y phases, and to explain increased solid solubility in the Al phases manifested by precipitation of ultrafine clusters of Ge. Furthermore, the reasons for the formation of these metastable eutectics under laser rapid solidification are discussed from the perspective of the competitive growth criterion.

Al-Ge eutectic↗
DOE OSTI2021

High Density of Strong yet Deformable Intermetallic Nanorods Leads to an Excellent Room Temperature Strength-Ductility Combination in a High Entropy Alloy

This paper introduces a new microstructural template for high entropy alloys (HEAs), where the face centered cubic (FCC) complex concentrated solid solution is reinforced with a high density of strong, yet deformable, nanorods of an ordered multi-component intermetallic L12 compound. Thermodynamic modeling has been employed to design this HEA with a large L12 volume fraction. Thermo-mechanical processing by isothermal annealing of the conventionally processed bulk cold-rolled alloy directly at precipitation temperatures, has been applied to produce a high density of uniformly distributed L12 nanorods within refined FCC grains, resulting from concomitant recrystallization and discontinuous precipitation processes. The nanorod morphology of the discontinuous L12 product has been established from three-dimensional atom probe tomography. The refined grains result in a complete coverage of the microstructure with discontinuously precipitated intermetallic nanorods. This nanorod strengthened HEA exhibits an exceptionally high room temperature yield strength of ~1630 MPa, good tensile ductility of ~15%, and an ultimate tensile strength of ~1720 MPa. Furthermore, a single L12 phase alloy, melted based on the precipitate composition in the two-phase FCC + L12 HEA, exhibits very high compressive deformability and strain hardenability, unusual for ordered intermetallic compounds. These results open a new strategy for design of fine-grained microstructures strengthened via ordered intermetallic phases, exploiting the beneficial effects of discontinuous precipitation, for achieving very high room temperature tensile strengths while maintaining good ductility.

Gwalani, Bharat↗
DOE OSTI2022

The effect of saturated thermal conduction on clouds in a hot plasma

We numerically investigate the internal evolution of multiphase clouds, which are at rest with respect to an ambient, highly ionized medium (HIM) representing the hot component of the circumgalactic medium. Time-dependent saturated thermal conduction and its implications like condensation rates and mixing efficiency are assessed in multiphase clouds. Our simulations are carried out by using the adaptive mesh refinement code FLASH. The model clouds are initially in both hydrostatic and thermal equilibrium and are in pressure balance with the HIM. Thus, they have steep gradients in both temperature and density at the interface to HIM leading to non-negligible thermal conduction. Several physical processes are considered numerically or semi-analytically: thermal conduction, radiative cooling and external heating of gas, self-gravity, mass diffusion, and dissociation of molecules and ionization of atoms. It turns out that saturated thermal conduction triggers a continuous condensation irrespective of cloud mass. Dynamical interactions with ambient HIM all relate to the radial density gradient in the clouds: (1) mass flux due to condensation is the higher the more homogeneous the clouds are; (2) mixing of condensed gas with cloud gas is easier in low-mass clouds, because of their shallower radial density gradient; and thus (3) accreted gas is distributed more efficiently. A distinct and sub-structured transition zone forms at the interface between cloud and HIM, which starts at smaller radii and is much narrower as deduced from analytical theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗
DOE OSTI2023

Rapid grain refinement and compositional homogenization in a cast binary Cu50Ni alloy achieved by friction stir processing

Friction stir processing (FSP) has been increasingly adopted for joining and processing materials in automotive, aerospace and industrial construction. During FSP a dynamic competition between high-speed shear deformation and deformation-induced heating brings about a complex competition between multiple dynamic microstructural evolution mechanisms making it difficult to predict the microstructural evolution pathway. Hence, improved understanding of microstructural evolution mechanisms during FSP can be beneficial for continued growth in the adoption of FSP for demanding applications of future. Towards this goal, this study uses a model binary Cu – 50 at.% Ni alloy to clarify the effect of single and double pass FSP on the microstructural evolution of a coarse grained and compositionally heterogeneous cast microstructure. High energy synchrotron X-ray diffraction, electron backscatter diffraction, and nanoindentation are used to clarify the microstructural evolution due to FSP. The process of compositional homogenization of as-cast segregations is studied by energy dispersive spectroscopy and atom probe tomography. Our results show that a single fast FSP pass at 30 mm.s -1 produces a 100 µm deep layer of submicrometric and hall-petch hardened CuNi grains. The initial cast compositional heterogeneities in a micrometric scale is rapidly transformed to nano-sized domains, mainly confined at grain boundaries. Double pass FSP increases the penetration depth of the processed layer and leads to a 2.9 times grain growth relative to single pass FSP. Grain fragmentation, discontinuous dynamic recrystallization, grain growth, and twinning mechanisms are discussed. In conclusion, these results highlight the value of FSP for ultrafast grain refinement and compositional homogenization of cast alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid–base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

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

Thermally driven surface phase separation in intermetallic alloys

Intermetallic compounds are widely recognized for their high-temperature phase stability and resistance to composition and structural changes. However, we reveal athermally activated bulk-to-surface mass exchange mechanism that drives surface phase separation, resulting in the formation of surface precipitates with distinct composition andstructure from the bulk matrix. Using the archetypal β-NiAl system, we show that asymmetries in vacancy formation energies between Ni and Al atoms induce preferential Ni segregation to the surface, forming Ni-rich γ'-Ni 3 Al precipitates. By integrating in-situ electron microscopy, synchrotron X-ray absorption spectroscopy and first-principles computational modeling, we establish a direct mechanistic connection between bulk thermal defect dynamics, surface compositional evolution, and phase segregation behavior. This bulk-surface coupling mechanism can be a driver of surface phase separation in multicomponent alloys under thermal stress. In conclusion, these results refine the thermodynamic boundaries of intermetallic stability and provide insights into managingthe performance and durability of intermetallic alloys for demanding high-temperature applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Optimal invariant sets for atomistic machine learning

The representation of atomic configurations for machine learning models has led to numerous sets of descriptors. However, many descriptor sets are incomplete and/or functionally dependent. Incomplete sets cannot faithfully represent atomic environments. Yet complete constructions often suffer from a high degree of functional dependence, where some descriptors are functions of others. These redundant descriptors do not improve discrimination between atomic environments. We employ pattern recognition techniques to remove dependent descriptors to produce the smallest possible set that satisfies completeness. We apply this in two ways: First, we refine an existing description, the atomic cluster expansion. Second, we augment an incomplete construction, yielding a new message-passing neural network architecture that can recognize up to 5-body patterns. This architecture shows strong accuracy on state-of-the-art benchmarks while retaining low computational cost. Our results demonstrate the utility of this strategy to optimize descriptor sets across a range of descriptors and application datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Evaluation of the GaAs displacement damage metric using updated nuclear data

The emerging use of the physics-based athermal recombination-corrected displacement per atom (arc-dpa) model for the displacement damage efficiency has motivated a re-evaluation of the historical empirically-derived GaAs damage response function with the purpose of highlighting needs for future analytical and experimental work. The 1-MeV neutron damage equivalence methodology used in the ASTM E-722 standard for GaAs has been re-evaluated using updated nuclear data. This yielded a higher fidelity representation of the GaAs displacement kerma and, through the use of the refined PKA recoil energy-dependent damage efficiency model, an updated 1-MeV(GaAs) displacement damage function. This re-evaluation included use of the Norgett-Robinson-Torrens (NRT) model for an updated threshold treatment, rather than the sharp-threshold Kinchin-Pease model used in the current ASTM standard. The underlying nuclear data evaluations have been updated to use the ENDF/VIII.0 {sup 75}As and TENDL-2019 {sup 71}Ga/{sup 69}Ga evaluations. The displacement kerma and 1-MeV-equivalent damage responses were calculated using a modified NJOY-2016 code which allowed for refinements in some of the damage models. This paper shows that an updated displacement damage function, based upon the latest nuclear data, is consistent with the experimental data used to develop the current ASTM E-722 GaAs standard. Using a double ratio approach to compare the available experimental data with the calculated response, the average legacy double ratio was found to be 0.97 ± 0.05 and the average updated double ratio was found to be 0.94 ± 0.05. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomistic structure of the SARS-CoV-2 pseudoknot in solution from SAXS-driven molecular dynamics

SARS-CoV-2 depends on -1 programmed ribosomal frameshifting (-1 PRF) to express proteins essential for its replication. The RNA pseudoknot stimulating -1 PRF is thus an attractive drug target. However, the structural models of this pseudoknot obtained from cryo-EM and crystallography differ in some important features, leaving the pseudoknot structure unclear. We measured the solution structure of the pseudoknot using small-angle X-ray scattering (SAXS). The measured profile did not agree with profiles computed from the previously solved structures. Beginning with each of these solved structures, we used the SAXS data to direct all atom molecular dynamics (MD) simulations to improve the agreement in profiles. In all cases, this refinement resulted in a bent conformation that more closely resembled the cryo-EM structures than the crystal structure. Applying the same approach to a point mutant abolishing -1 PRF revealed a notably more bent structure with reoriented helices. This work clarifies the dynamic structures of the SARS-CoV-2 pseudoknot in solution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Large-eddy simulation of non-vaporizing sprays using the spectral-element method

Predictive simulations of high-pressure sprays require accurate representation of the turbulent gaseous flow field generated by liquid jet. Typically, the accuracy that can be obtained with low-order numerical methods (e.g. finite volume, finite element) is limited by stability issues in fine grids and the order of convergence of the method. In this work, we resolve the turbulent flow field in an Eulerian manner using the high-order spectral element method, coupled with a Lagrangian parcels approach to model the atomizing liquid jet. Large eddy simulations of single-hole sprays under non-evaporative conditions were conducted and compared against experimental data from Margot et al. (2008) and Spray A data from the Engine Combustion Network. The sensitivity of liquid penetration and droplet sizes to different breakup model parameters was studied. The effect of different numerical parameters, such as polynomial order of the solution (grid resolution), on liquid penetration was also analyzed. The method achieved grid-independent results using p-refinement, achieving finer resolution (by a factor of x 1.7 - x 3.5) in the gas-phase solution than in state-of-the-art simulations using the finite-volume method. Here the results showed good agreement with experimental data, demonstrating the ability of the current method to accurately capture liquid penetration and the shape of the spray.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Large Angle Rocking Beam Electron Diffraction Utilizing Electron Direct Detector

Electron diffraction of a crystal is fundamentally a function of the potential of that crystal. The intensity of electron diffraction patterns is, as a result, sensitive to the charge density of atoms and bonding inside crystals. Experimentally, the traditional method to probe this information is quantitative Convergent Beam Electron Diffraction (CBED). Quantitative CBED or QCEBD is a method which uses dynamic diffraction theory to quantify CBED intensities and to extract information about crystal structure and bonding. The sensitivity of dynamical scattering is leveraged to measure crystal symmetry and crystal structure factors. At its limit, electron structure factors are measured at high accuracy allowing access to chemical bonding information, which corroborate theoretical calculations through multipole model refinements of the experimental charge density. The primary limit to QCBED is the requirement of a non-overlapping convergent beam. This subsequently limits QCBED to crystals with small unit cells which are stable under the focused probe.

Busch, Robert↗
DOE OSTI2023

Synthesis, crystal and electronic structure of the Zintl phase Ba 16 Sb 11 . A case study uncovering greater structural complexity via monoclinic distortion of the tetragonal Ca 16 Sb 11 structure type

The binary Zintl phase Ba 16 Sb 11 has been synthesized and structurally characterized. Detailed studies via single-crystal X-ray diffraction methods indicate that although Ba 16 Sb 11 appears to crystallize in the tetragonal Ca 16 Sb 11 structure type (space group $P\bar{4}2_1$m with a=13.5647(9) Å, c=12.4124(12)Å, Z=2, R 1 = 3.14%; wR 2 = 4.77%), there exists an extensive structural disorder. Some Ba 16 Sb 11 crystals were found to be monoclinic and the structure was solved and refined in space group P2 1 (a=18.3929(12) Å, b=13.5233(8) Å, c=18.3978(12) Å, β=94.6600(10)°; Z=4, R 1 =5.84 %; wR 2 =9.58 %). The latter corresponds to a 2-fold superstructure of the tetragonal one, which provides a disorder-free structural model. In both descriptions, the disordered tetragonal and the ordered monoclinic superstructure, the basic building units that make up the structure of this Ba-rich compound are pairs of face-shared square antiprisms of Ba atoms, which are centered by Sb atoms. The dimerized antiprisms are linked into parallel chains via square prisms of Ba atoms, which are also centered by Sb atoms. The Zintl concept can be applied in a straightforward manner and as result, the structure of Ba 32 Sb 22 (=2×Ba 16 Sb 11 ) can be rationalized as (Ba 2+ ) 32 (Sb 3– ) 20 [Sb 2 ] 4– . Notably, the partitioning of the valence electrons is done taking into an account the homoatomic Sb–Sb contacts (d=3.01 Å), which can be clearly distinguished in the lower symmetry space group. Electronic structure calculations of Ba 16 Sb 11 are in good accordance with the Zintl rationalization and predict a semiconductor with a band gap of 0.77 eV.

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

Atomic Layer Deposition for Improved Biomass Conversion Catalysts

Heterogeneous catalysts are a key enabler of the transition towards a sustainable, bio-based economy for fuels and chemicals. However, the harsh conditions in many biomass conversion processes lead to nanoparticle sintering, support collapse, and metal leaching in conventional PGM catalysts. Next-generation catalysts must be developed to address these stability challenges. This presentation will discuss ongoing work between the Catalytic Carbon Transformation & Scale-Up Center at NREL and various industrial partners to develop scalable, cost-effective atomic layer deposition (ALD) coatings for improving the performance of biomass conversion catalysts. Substantial focus will be devoted to the benefits of TiO2 ALD coatings on supported Pd hydrogenation catalysts. Ten cycles of TiO2 ALD were found to dramatically improve the activity of a conventional Pd/Al2O3 catalyst towards aromatic hydrogenation, despite partial coverage of the Pd sites by the ALD layer. Subsequent advanced characterization and atomic-scale computational modeling revealed that the ALD coating influenced the adsorption strength of hydrogenation surface intermediates in unexpected ways, leading to higher reaction rates. Reaction testing after exposure to sulfur impurities, high temperature oxidation, and hydrothermal treatment demonstrated the improved stability of the ALD-coated catalyst. Additionally, the ALD synthesis process was found to be scalable over two orders of magnitude with minimal deviation in synthesized catalyst properties. These results were contextualized with cost models of industrial ALD coating and aromatic hydrogenation processes to further refine the value proposition of ALD coatings. Given the demonstrated improvements in hydrogenation, TiO2 ALD coatings have also been applied to supported Pt catalysts for use in other biomass conversion reactions, such as hydrodeoxygenation for production of sustainable aviation fuel (SAF) and the oxidation of glucose for bio-derived gluconic acid. Recent findings from these experimental campaigns will also be shared. ALD technology holds great potential in the development of next-generation catalysts for biofuels and bioproducts and this work constitutes an important examination of the impact of ALD coatings in a variety of reaction environments.

bioenergy↗
DOE OSTI2022

Mechanical properties of zeolite-templated carbons from approximate density functional theory calculations

Zeolite-templated carbon (ZTC) is a unique porous carbonaceous material whose structure is ordered at the nanometre scale, enabling a representative periodic description at the atomistic level. Utilizing an existing, well-defined reference model for ZTCs, a structural library of varying compositions was developed by refinement using density-functional tight-binding (DFTB) potentials parameterized for materials science applications. We first determined the quantum chemical-refined structures of models with CH, CHO, CHON, CHOB, and CHOBN compositions with various degrees of heteroatom substitution. These structural models comprise the characteristic morphological features of highly porous carbon materials, such as open-blade surfaces, edges, saddles, and closed-strut formations, spanning a range of curvatures and characteristic sizes. Second, we carried out alternating compression and expansion of the CHO model unit cell to determine the lowest energy structure as well as to obtain its bulk modulus in order to demonstrate a close connection between macroscopic observations and atomic-scale structures. Further, the agreement between experimental measurements and the computational model is remarkable and demonstrates the power of approximate density functional theory as a cost-effective computational tool with chemical accuracy for the investigation of structure/property relationships in real-world carbon-based solids.

03 NATURAL GAS↗
DOE OSTI2022

MELCOR Model Development of the General Atomics Fast Modular Reactor

General Atomics Electromagnetic Systems (GA-EMS) is currently developing a fast modular reactor (FMR) capable of providing safe, carbon-free electricity by 2035. As part of the development process, researchers at the University of Wisconsin-Madison (UWM) use MELCOR 2.2 to develop a simulation model of the proposed design to perform accident analyses. This paper summarizes the efforts of researchers at UWM in the creation of the FMR model. This model was built based on initial FMR design parameters. The design parameters will change as the design evolves and is refined.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

AToM: Advanced Tokamak Modeling Environment

Stability constraints play an important role in integrated modeling. The global plasma stability and macro-instabilities such as internal kink modes, neoclassical tearing modes, and edge localized perturbations limit the plasma performance and can result in large-scale transient events and plasma disruptions. These macro-instabilities can be also beneficial to the plasma performance. For example, the peeling models often lead to the edge localized modes (ELMs). However, the formation of stationary edge harmonic oscillations (EHOs) due to nonlinear interaction of peeling modes can result in a transition to the Quiescent H-mode (QH-mode) without ELMs. The profiles predicted with transport models need to be refined using the MHD stability constraints. This is important especially for transient and nonlinear stages of discharges such as ramp-up access to hybrid and steady state operation, and L- to H-mode transition as well as transition to the QH-mode. Several stability and MHD codes are already included in the OMFIT framework. These codes include BALOO, ELITE, GATO, M3DC1, MARS, and NIMROD. However, the verifications of stability constraints are currently mostly excluded from transport modeling workflows. The only exception is the EPED module which includes the ELITE predictions to limit the pedestal height. Here, we improved the OMFIT workflow to include stability calculations. Including the stability conditions to the integrated modeling workflow improved the robustness of the predictive modeling discharges. Being implemented in the workflow, the stability conditions can be also used for the stability analysis of experimental data. This improved the physics understanding of various discharge scenarios and can be used in experiment planning.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗