Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solid modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Synthesis of magnesiowüstite nanocrystallites embedded in an amorphous silicate matrix via low energy multiple ion implantations

The synthesis process is presented for experimentally simulating modifications in cosmic dust grains using sequential ion implantations or irradiations followed by thermal annealing. Cosmic silicate dust analogues were prepared via implantation of 20–80 keV Fe - , Mg - , and O - ions into commercially available p-type silicon (100)wafers. The as-implanted analogues are amorphous with a Mg/(Fe+Mg) ratio of 0.5 tailored to match theoretical abundances in circumstellar dusts. Before the ion implantations were performed, Monte-Carlo-based ion-solid interaction codes were used to model the dynamic redistribution of the implanted atoms in the silicon substrate. 600 keV helium ion irradiation was performed on one of the samples before thermal annealing. Two samples were thermally annealed at a temperature appropriate for an M-class stellar wind, 1000 K, for 8.3 h in a vacuum chamber with a pressure of 1 x 10 -7 torr. The elemental depth profiles were extracted utilizing Rutherford Backscattering Spectrometry (RBS) in the samples before and after thermal annealing. X-ray diffraction (XRD)analysis was employed for the identification of various phases in crystalline minerals in the annealed analogues. Transmission electron microscopy (TEM) analysis was utilized to identify specific crystal structures. RBS analysis shows redistribution of the implanted Fe, Mg, and O after thermal annealing due to incorporation into the crystal structures for each sample type. XRD patterns along with TEM analysis showed nanocrystalline Mg and Fe oxides with possible incorporation of additional silicate minerals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Flexible Zn-MOF with Rare Underlying scu Topology for Effective Separation of C6 Alkane Isomers

Adsorptive separation by porous solids provides an energy-efficient alternative for the purification of important chemical species compared to energy-intensive distillations. Particularly, the separation of linear hexane isomers from its branched counterparts is crucial to produce premium grade gasoline with high research octane number (RON). Herein, we report the synthesis of a new, flexible zinc-based metal-organic framework, [Zn 5 (μ 3 -OH) 2 (adtb) 2 (H 2 O) 5 ·5 DMA] (Zn-adtb), constructed from a butterfly shaped carboxylate linker with underlying (4,8)-connected scu topology capable of separating the C6 isomers nHEX, 3MP, and 23DMB. The sorbate-sorbent interactions and separation mechanisms were investigated and analyzed through in situ FTIR, solid state NMR measurements and computational modeling. Here these studies reveal that Zn-adtb discriminates the nHEX/3MP isomer pair through a kinetic separation mechanism and the nHEX/23DMB isomer pair through a molecular sieving mechanism. Column breakthrough measurements further demonstrate the efficient separation of linear nHEX from the mono- and dibranched isomers.

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

Machine learning potential assisted exploration of complex defect potential energy surfaces

Abstract Atomic-scale defects generated in materials under both equilibrium and irradiation conditions can significantly impact their physical and mechanical properties. Unraveling the energetically most favorable ground-state configurations of these defects is an important step towards the fundamental understanding of their influence on the performance of materials ranging from photovoltaics to advanced nuclear fuels. Here, using fluorite-structured thorium dioxide (ThO 2 ) as an exemplar, we demonstrate how density functional theory and machine learning interatomic potential can be synergistically combined into a powerful tool that enables exhaustive exploration of the large configuration spaces of small point defect clusters. Our study leads to several unexpected discoveries, including defect polymorphism and ground-state structures that defy our physical intuitions. Possible physical origins of these unexpected findings are elucidated using a local cluster expansion model developed in this work.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Systematic softening in universal machine learning interatomic potentials

Machine learning interatomic potentials (MLIPs) have introduced a new paradigm for atomic simulations. Recent advancements have led to universal MLIPs (uMLIPs) that are pre-trained on diverse datasets, providing opportunities for universal force fields and foundational machine learning models. However, their performance in extrapolating to out-of-distribution complex atomic environments remains unclear. In this study, we highlight a consistent potential energy surface (PES) softening effect in three uMLIPs: M3GNet, CHGNet, and MACE-MP-0, which is characterized by energy and force underprediction in atomic-modeling benchmarks including surfaces, defects, solid-solution energetics, ion migration barriers, phonon vibration modes, and general high-energy states. The PES softening behavior originates primarily from the systematically underpredicted PES curvature, which derives from the biased sampling of near-equilibrium atomic arrangements in uMLIP pre-training datasets. Our findings suggest that a considerable fraction of uMLIP errors are highly systematic, and can therefore be efficiently corrected. We argue for the importance of a comprehensive materials dataset with improved PES sampling for next-generation foundational MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A representation-independent electronic charge density database for crystalline materials

Abstract In addition to being the core quantity in density-functional theory, the charge density can be used in many tertiary analyses in materials sciences from bonding to assigning charge to specific atoms. The charge density is data-rich since it contains information about all the electrons in the system. With the increasing prevalence of machine-learning tools in materials sciences, a data-rich object like the charge density can be utilized in a wide range of applications. The database presented here provides a modern and user-friendly interface for a large and continuously updated collection of charge densities as part of the Materials Project. In addition to the charge density data, we provide the theory and code for changing the representation of the charge density which should enable more advanced machine-learning studies for the broader community.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Nanopores with dynamic pore opening diameter

Solid state nanopores have emerged as model systems for understanding transport properties on the nanoscale. They serve as templates for both preparing mimics of biological channels and designing biological sensors. Unlike their biological inspirations, the majority of nanopores prepared thus far, however, have been structurally static devices such that the pore opening diameter is fixed. If we could prepare nanopores whose opening diameter fluctuated in time with controlled amplitude at known locations in the pore, we could create ionic memristors as well as achieve new transport modes. Here we present ∼10 nm diameter single nanopores drilled through a 10 nm thick gold layer positioned on top of a 30 nm thick silicon nitride film. Two types of devices were prepared; one containing single stranded DNA and the other containing hairpin DNA attached to the discrete layer of gold using thiol chemistry. When an external electric field was applied across a nanopore with single stranded DNA, the nanoconfined DNA molecules exhibited steric and electrical constraints that led to memristor-like behavior in the current–voltage curves. The degree of hysteresis was controlled by salt concentration, magnitude of voltage and pore diameter. In contrast, nanopores containing DNA hairpins conducted similar currents in forward and reverse bias in agreement with the rigidity of the hairpin molecule. The experiments are explained by Brownian dynamics simulations that reveal voltage and salt concentration induced changes in DNA extension. The degree of DNA extension was also found to be dependent on the location of the molecules along the pore axis. The nanopores presented here provide the first steps towards preparation of non-equilibrium nanopore systems.

Vlassiouk, Ivan [ORNL] (ORCID:0000000254940386)↗
DOE OSTI2020

Massive formation of early diagenetic dolomite in the Ediacaran ocean: Constraints on the “dolomite problem”

Paleozoic and Precambrian sedimentary successions frequently contain massive dolomicrite [CaMg(CO 3 ) 2 ] units despite kinetic inhibitions to nucleation and precipitation of dolomite at Earth surface temperatures (<60 °C). This paradoxical observation is known as the “dolomite problem.” Accordingly, the genesis of these dolostones is usually attributed to burial–hydrothermal dolomitization of primary limestones (CaCO 3 ) at temperatures of >100 °C, thus raising doubt about the validity of these deposits as archives of Earth surface environments. We present a high-resolution, >63-My-long clumped-isotope temperature (T Δ47 ) record of shallow-marine dolomicrites from two drillcores of the Ediacaran (635 to 541 Ma) Doushantuo Formation in South China. Our T ∆47 record indicates that a majority (87%) of these dolostones formed at temperatures of <100 °C. When considering the regional thermal history, modeling of the influence of solid-state reordering on our T Δ47 record further suggests that most of the studied dolostones formed at temperatures of <60 °C, providing direct evidence of a low-temperature origin of these dolostones. Furthermore, calculated δ 18 O values of diagenetic fluids, rare earth element plus yttrium compositions, and petrographic observations of these dolostones are consistent with an early diagenetic origin in a rock-buffered environment. We thus propose that a precursor precipitate from seawater was subsequently dolomitized during early diagenesis in a near-surface setting to produce the large volume of dolostones in the Doushantuo Formation. Our findings suggest that the preponderance of dolomite in Paleozoic and Precambrian deposits likely reflects oceanic conditions specific to those eras and that dolostones can be faithful recorders of environmental conditions in the early oceans.

58 GEOSCIENCES↗
DOE OSTI2025

Protracted core formation and impact disruptions shaped the earliest outer Solar System planetesimals

The distinct compositions of metallic cores from noncarbonaceous (NC) and carbonaceous (CC) iron meteorite parent bodies (IMPBs) reflect differences in accretion and differentiation histories of earliest inner and outer Solar System planetesimals. Compared to NC IMPBs, CC IMPBs have smaller, sulfur (S)–poor, highly siderophile element (HSE)–enriched cores and younger core formation ages. However, the origins of these differences remain debated. Using equilibrium partitioning models between the S-poor solid and the S-rich liquid metal, we argue that HSE enrichment in IID, IIF, IIIF, and IVB cores resulted from a multistage evolutionary sequence: (i) segregation of S-rich, HSE-depleted protocores during initial planetesimal heating; (ii) collisional disruption before S-poor metal segregation; (iii) reaccretion of mantle fragments into daughter planetesimals; and (iv) further 26 Al-driven heating producing HSE-enriched, S-poor cores. We suggest that iron meteorites from these CC IMPBs originate from such second-generation cores. Accounting for “missing” S-rich protocores helps reconcile several NC-CC IMPB differences and highlights the role of early collisional processing in shaping planetesimal chemical evolution.

Geosciences↗
DOE OSTI2024

Thermodynamic Modeling of Point Defects in Triple Conducting Perovskite Ba 0.95 La 0.05 FeO 3- $_δ$ with Incorporation of the Hydride Defect Formation Reaction for Solid Oxide Cells

Distinct from the proton defect, the hydride defect species may be present in certain perovskite materials in reducing environments such as in fuel electrodes of solid oxide cells (SOCs) or on the reducing side of ceramic membranes. A generalized defect thermodynamic model was developed for the triple-conducting perovskites (La,Ba)Fe 1-x M x O 3-δ (M = Y and Zr) to allow inclusion of the hydride defect formation reaction in addition to the other three main defect reactions, namely, the oxygen vacancy formation, hydration, and charge disproportionation reactions. This comprehensive defect model also allows the incorporation of polynomial functional forms of oxygen nonstoichiometry δ to describe the defect reaction energies and entropies and to enable refinements of the defect reaction equilibrium constants in the defect thermodynamic analysis. As a first step, the developed model is applied to the Ba 0.95 La 0.05 FeO 3-δ material as an illustrative system to obtain its Brouwer diagrams with both the proton and hydride defects in relevant SOC conditions, particularly for more reducing environments. In conclusion, the results provide direct guidance on the influence of electronic and ionic defect concentrations upon thermodynamic properties and ultimately on the performance of Ba 0.95 La 0.05 FeO 3-δ and potentially other (La,Ba)Fe 1-x M x O 3-δ perovskite materials involved in SOC applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Development of Corrosion- and Erosion-Resistant Coatings for Advanced Ultra-Supercritical Materials

This final report summarized the research efforts and major findings of the Phase I Project “Development of Corrosion- and Erosion-Resistant Coatings for Advanced Ultra-Supercritical Materials”, for the period of October 1, 2019 – Sept. 30, 2021. This project is a collaborative endeavor between Tennessee Tech University, Purdue University, Oak Ridge National Laboratory, Siemens Corporation, and Eastern Plating, LLC, aiming at improving the durability and lifetime of high-pressure (HP) steam turbine blades in advanced ultra-supercritical (A-USC) coal-fired power plants through the development of corrosion/erosion-resistant coatings manufactured via a low-cost electrolytic codeposition process. While Tribaloy alloy T-400C was identified by the U.S. A-USC Materials Consortium as a promising coating composition, further composition optimization is needed to enhance its corrosion and erosion resistance for protecting the A-USC Ni-base turbine components. An integrated computational and experimental approach was employed to optimize coating composition/microstructure and processing parameters. In order to identify candidate coating compositions that could offer balanced properties, thermodynamic calculations were performed to explore the γ+Laves composition space in the Co-Ni-Cr-Mo-Si system with different alloying additions at 600-800 °C. Guided by the calculation results, experimental assessment of selected alloys led to the development of a new generation of Tribaloy compositions with the optimal levels of Cr, Mo and Si, reactive element (e.g., 0.4-0.6 wt.% Y) and other alloying additions. The low-cost and non-line-of-sight electro-codeposition process was employed to deposit a Ni(Co)-CrMoSiY composite coating on commercial Haynes 282 (H282) Ni-base alloy. A diffusion treatment was subsequently applied to convert the composite to the Tribaloy-type coating. Both the codeposition parameters and heat treatment conditions were varied to achieve the desired coating composition, microstructure and phase constituents. In addition, since additive manufacturing (AM) may be an alternative cost-saving option for potential A-USC turbine repair, laser direct deposition was explored to fabricate the H282 alloy with minimal defects. The electro-codeposited Tribaloy coating was also applied to the AM H282 substrate to demonstrate the viability of the coating process in improving the surface finish of AM alloys. Both high-temperature oxidation performance and solid particle erosion (SPE) resistance of model alloys and electro-codeposited coatings were evaluated. About 25 model alloys with various Cr/Mo ratios and reactive element levels, as well as partial substitution of Mo with Nb were evaluated with regard to their oxidation resistance in both air and pure steam at 760-800°C. Compositions based on Ni(Co)-20Cr-18Mo-2.6Si-0.6Y (wt.%) showed significantly improved oxidation resistance over the baseline T400-C. Based on the alloy development results, three generations of new Tribaloy coatings (Gen-1, Gen-2, and Gen-3) with various Mo/Cr contents and Y levels were prepared via electro-codeposition and their microstructure/performance were evaluated. Outstanding air and steam oxidation resistance was achieved for the Gen-2 and Gen-3 coatings. Furthermore, while the SPE resistance of the coatings depended on many factors such as temperature, environment, erodent, velocity and impact angle, the coated samples exhibited similar or better SPE resistance compared to the H282 substrate when magnetite was used as erodent (which is a realistic erodent in A-USC steam turbines). The new Tribaloy coatings also had good long-term compatibility with the H282 alloy substrate. The two large-sized rotating barrels were designed, constructed, and employed to coat dummy HP blades. Uniform coating thickness and microstructure were achieved at various blade locations. Also, a preliminary techno-economic analysis of the proposed coating process was conducted to quantify the cost-effectiveness and to assess the commercial viability of the corrosion- and erosion-resistant coatings. It is estimated that a cost reduction of ~30% could be achieved with the electro-codeposition coating process over the state-of-the-art high velocity oxygen fuel (HVOF) thermal spray. Compared to the leading Tribaloy coating technologies such as HVOF and plasma transfer alloying, electro-codeposition based process has advantages such as low-cost process equipment, uniform deposition even for complex shapes, low levels of contaminants/porosities, reduction of powder waste, and potentially better surface finish and longer turbine service life. The Phase 1 study has demonstrated that it is feasible to develop an electro-codeposited Tribaloy coating with balanced corrosion and erosion properties, even though additional research efforts such as further coating process scale-up and longer-term performance evaluation under realistic A-USC conditions are clearly needed.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Modeling and Analysis of Climate Variation Effects on Fixed-Bed Direct Air Capture Systems

This poster was presented at the 2024 FECM/ NETL Carbon Management Research Project Review Meeting held August 5th-9th in Pittsburgh, PA. The poster discusses the modeling, optimization, and analysis of a solid sorbent based direct air capture process. Analysis of the performance of the system when exposed to varying ambient conditions is also presented.

Hughes, Ryan↗
DOE OSTI2024

Designing SOC electrodes for better lifetime performance

Presentation for talk at the 2024 Pacific Rim Meeting of the Electrochemical Society October 6-11, 2024. The talk covers a review of NETL's solid oxide cell (SOC) performance degradation modeling framework and present case studies on how simulations are used to design better SOC electrodes.

Abernathy, Harry↗
DOE OSTI2023

Geochemical Modeling Using OLI Systems Mixed Solvent Electrolyte Model in Support of the Waste Isolation Pilot Plant

Aqueous electrolyte thermodynamic models were developed for synthetic Salado (GWB) and Castile (ERDA-6) brines to inform the experimental design and to facilitate laboratory data interpretation in geochemical studies supporting disposal in the Waste Isolation Pilot Plant (WIPP). The thermodynamic models, which include the relevant vapor-liquid-solid equilibria for the oxidation-reduction conditions of interest, enable computation of system parameters such as activity-based pH, E h [vs. standard hydrogen electrode (SHE)], ionic strength, total dissolved solids, and solution density; the chemical speciation of all phases; saturation tendencies for thermodynamically feasible solids; and changes in brine chemistry because of the addition of H 2 (g), argon, MgO, Fe(0), and other waste form components. The impacts of oxidation-reduction potential (ORP) and pH on aqueous speciation and dominant solid phases were also assessed with the aid of real-solution Pourbaix diagrams for key elements of interest. In addition to the OLI modeling results, limited measurement data for pH and ORP standards are provided, analyzed, and discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Two Dimensional Topology Optimization of Heat Exchangers with the Density and Level-Set Methods

We design heat exchangers using two topology optimization approaches: the density, i.e. volume fraction and level set methods. Our goal is to maximize the heat exchange between two fluids in separate channels while constraining the pressure drop across each channel. The heat exchanger is modeled with a coupled thermal-flow formulation. The flow is governed by an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation with high Peclet number. We solve one set of Stokes-Brinkman equations per fluid. Each Brinkman term in the flow equation serves to model the other phase as a solid, thereby preventing mixing. We first represent the solid and fluid phases using a volume fraction variable and apply a SIMP-like penalization in the Brinkman term to drive the optimization to a discrete design. The cost and constraint function derivatives are automatically calculated with the library pyadjoint and the optimization is performed by the Method of Moving Asymptotes. In a second optimization formulation, we use the level set approach to define the interface that separates the two fluids. Pyadjoint calculates the shape derivatives of the cost and constraint functions and the Hamilton-Jacobi advects the interface, allowing for topological changes. We present results in two dimensions and discuss the advantages and disadvantages of each approach.

42 ENGINEERING↗
DOE OSTI2024

A hybrid finite volume method and smoothed particle hydrodynamics approach for efficient and accurate blast simulations

Modeling strong shock waves in fluids remains a persistent challenge in computational physics. Essential to research efforts in industry and defense, numerous methods have been devised to improve the accuracy and efficiency of shock simulations. A novel, hybrid Finite Volume Method (FVM)-Smoothed Particle Hydrodynamics (SPH) approach is capable of further improving efficiency and retaining accuracy by exploiting the favorable characteristics of each respective method. This hybrid approach is presented for shock capturing in compressible fluids. The Python framework Pyro2 is employed to simulate a coarse FVM mesh, while the Python framework PySPH is utilized to model the fluid in regions with high gradients through SPH particles. The performance of the hybrid FVM-SPH scheme, compared to the individual FVM and SPH methods, is assessed in 1 kt and 10 kt blast simulations. Our results indicate that the hybrid approach offers higher computational efficiency than SPH while preserving its accuracy and characteristics. The hybrid approach had a relative speedup of 11.3x and 22.3x over the FVM and SPH approaches for the 1 kt simulation and a relative speedup of 14.7x and 20.9x over the FVM and SPH approaches for the 10 kt simulation. The hybrid SPH algorithm enables future compressible fluid simulations with more extensive capabilities than grid-based methods alone, presenting potential applications in modeling fluid-structure interactions and solid deformation and fracturing in blast simulations.

Myers, Conner↗
DOE OSTI2025

Angular Correlation Date Measurements with the GeRMAC system

Advanced modeling and simulation efforts have improved at Idaho National Laboratory in recent years with a solid foundation of experimental results. Current computational methods represent significant modeling capabilities but are limited by the accuracy and availability of nuclear data. The creation of pre- and post-processing software tools to address these limitations is fundamental to the improvement of nuclear science modeling capacities. One aspect of predictive modeling tools deals with gamma-rays emitted from radionuclides, including fissile or fissionable material, fission products, or activation products, produced in reactor experiments or other neutron environments. The resulting radionuclides decay in unique ways, providing complications upon measurement as a result of random and cascade, or true, coincidence summing. These effects are not easily quantified during modeling efforts of gamma-ray source terms., The germanium rotational measurements for angular correlation (GeRMAC) system was built to quantify the relative angles for gamma rays emitted by radionuclides of interest to investigate true coincidence, or cascade, summing as well as the nuclear energy levels of decay schemes of interest. Proof of concept studies utilize a series of laboratory check sources to provide validity, and it will soon be used to perform the same measurements for fission products of interest. The resulting data can be used to implement into a Monte Carlo code, such as Geant4, to provide more precise gamma-ray source terms following irradiations of materials.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗