Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface capturing”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Phase-field modeling of diffusion bonding in 316H stainless steel: Impact of processing conditions on grain morphology and bonding quality

A novel multi-phase, multi-component phase‐field model is presented to study the diffusion bonding of 316H stainless steel. Combined with targeted experimental investigations, this model simulates the bond-growth process and predicts the bonding quality. Unlike previous models, our approach captures the simultaneous evolution of voids and grain structures, while quantifying bonding quality using defined bonding ratio. A comprehensive analysis of bond process control is performed by changing temperature, pressure and surface roughness observing the resulting bond structure, which is consistent with experimental observations and analytical predictions. Temperature is determined to be the dominant factor, with the transition from a flat to a robust bond occurring between 1000 °C and 1050 °C. At the ideal bonding temperature of 1050 °C, a surface roughness exceeding 0.6 μm or an applied stress below 4 MPa results in poor bonding quality. Beyond this, higher pressures and smoother surfaces reduce void size, accelerate void shrinkage, and lead to improved bond integrity. This diffuse-interface model can be extended to other material systems if supplied with appropriate thermodynamic and kinetic data. In conclusion, this makes it an effective modeling platform for optimizing high-temperature diffusion bonding and developing reliable bonded components such as compact heat exchangers.

Diffusion bonding↗
DOE OSTI2022

SAM-ML: Integrating data-driven closure with nuclear system code SAM for improved modeling capability

Advanced reactors often involve complicated thermal-fluid (T-F) phenomena. Modeling such phenomena with the traditional one-dimensional (1-D) system code is a challenging task. The System Analysis Module (SAM), a modern nuclear system code, has developed a coarse mesh multi-dimensional (multi-D) flow model to capture the spatial effect of T-F phenomena in advanced reactors. As a coarse mesh solver, constitutive relations are required for SAM's multi-D model for unresolved fine-scale physics, such as turbulence. Here this work presents a novel approach that integrates neural networks as data-driven closure for SAM's multi-D flow model. The data-driven closure is trained with fine-resolution data to ensure its accuracy while maintaining a coarse mesh setup to ensure its efficiency and consistency with SAM. We demonstrate the applicability of this SAM-ML capability in an open volume thermal stratification problem, where a neural network model serves as the eddy viscosity closure. A customized interface between the neural network and SAM is developed to ensure flexible and efficient data exchange. The SAM-ML results demonstrate superior performance compared to SAM's built-in zero-equation eddy viscosity closure. The case study shows that although the generalization capability of the data-driven closure still needs to be improved for different transient case or different geometric setup, SAM -ML demonstrates good potential for challenging simulation problems with improved accuracy and computational efficiency.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗
DOE OSTI2024

First principles density functional theory study of tritium species adsorption on Ni(111) surface and diffusion in nickel-sublayer for tritium storage

The nickel-plated zircaloy-4 is used as a tritium ( 3 H) getter in the tritium-producing burnable absorber rods (TPBARs) to capture 3 H produced in the 6 Li-riched annular γ-LiAlO 2 pellet under neutron irradiation. The experimental data and our previous theoretical results showed that the 3 H species produced from the γ-LiAlO 2 pellet were mainly 3 H 2 and 3 H 2 O. These 3 H species diffuse from the surface of the LiAlO 2 pellet across vacuum to the nickel-plated zircaloy-4 getter and then further diffuse into the getter to chemically form metal hydrides. While a number of studies show that oxygen binds strongly as compared to 3 H on the nickel (Ni) layer, the detailed mechanism of 3 H species absorption and diffusion across the Ni plate and Ni/Zr interface are still unclear. By employing density functional theory calculations, here we explored the 3 H 2 and 3 H 2 O species adsorption and dissociation on the Ni(111) surface and diffusion into the Ni sublayer. Our results indicated that the 3 H 2 and 3 H 2 O dissociate on the Ni(111) surface. The NiO x and Ni(O 3 H) x could be formed in the Ni layer due to the higher oxygen (O) diffusion energy barrier and formation of Ni vacancy defects. The oxygen was found to be retained in the Ni layer from diffusing across the Ni–Zr interface. This was revealed by comparing the diffusion barriers for 3 H with O. 3 H was found to have nearly three times smaller diffusion barrier than for O, making 3 H comparatively easier to diffuse through the Ni layer. In conclusion, the obtained results provide guidelines for experimental measurements on 3 H retention behavior in TPBARs and may open further avenues to explore the impurity effects on 3 H diffusion and storage at the Ni/zircaloy interfaces.

Tafen, De Nyago [National Energy Technology Lab. (↗
DOE OSTI2024

Methane partial oxidation under periodic reaction conditions on Pt/Al 2 O 3

The increasing interest in utilizing methane, the primary component of natural gas, for chemical production has spurred research into methane partial oxidation (MPO) as an alternative to traditional steam methane reforming (SMR). MPO has lower energy requirements and potential for carbon capture, making it an attractive option for hydrogen production. Challenges remain, however, such as carbon deposition leading to degradation and achieving high hydrogen selectivity. Here, the impact of periodic reactor operation on MPO over a Pt/Al 2 O 3 catalyst was studied, primarily via varying reactor inlet compositions. Experiments were conducted using periodic operation strategies to assess the influence of changing reactant inlet concentrations on hydrogen formation during MPO. The results suggest that cycling between mixtures with low and high oxygen content can lead to transient hydrogen formation rates that surpass those achieved at steady state. Control experiments and density functional theory (DFT) calculations show that enhanced hydrogen formation can be attributed to the reaction between CO with hydroxyl groups at the metal and alumina support interface. Furthermore, this work underscores the critical role of surface coverages at the metal–support interface and suggests avenues for future exploration, including alternative support materials with higher OH mobility and changes in the cycling scheme to enhance catalyst performance under periodic conditions.

03 NATURAL GAS↗
DOE OSTI2022

An electrochemical mesoscale tool for modeling the corrosion of structural alloys by molten salt

Understanding the impact of microstructure on corrosion rates can aid the development of corrosion-resistant alloys for molten salt reactors. Here in this work, we develop an electrochemical phase-field model for capturing the microstructure-dependent corrosion of structural alloys by molten salts. As a demonstration problem, we apply this model to capture the selective depletion of Cr from Ni-Cr grain boundaries during corrosion in molten FLiBe salt. We perform sensitivity analysis and model verification on 1D simulations to confirm that the model predicts diffusion-limited kinetics. The model is validated using 1D, 2D, and 3D simulations against experimental data for Ni-5Cr and Ni-20Cr corrosion in molten FLiBe. The 1D simulations predict the corrosion behavior with reasonable accuracy when using an effective diffusion coefficient that accurately represents the grain boundary diffusion. 2D simulations that represent the grain structure underpredict the corrosion. 3D simulations that represent the grain structure predict the corrosion with reasonable accuracy. The corrosion rate predicted by the 3D simulations is proportional to the average grain size at the alloy/salt interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Reactive capture and electrochemical conversion of CO 2 with ionic liquids and deep eutectic solvents

Ionic liquids (ILs) and deep eutectic solvents (DESs) have tremendous potential for reactive capture and conversion (RCC) of CO 2 due to their wide electrochemical stability window, low volatility, and high CO 2 solubility. There is environmental and economic interest in the direct utilization of the captured CO 2 using electrified and modular processes that forgo the thermal- or pressure-swing regeneration steps to concentrate CO 2 , eliminating the need to compress, transport, or store the gas. The conventional electrochemical conversion of CO 2 with aqueous electrolytes presents limited CO 2 solubility and high energy requirement to achieve industrially relevant products. Additionally, aqueous systems have competitive hydrogen evolution. In the past decade, there has been significant progress toward the design of ILs and DESs, and their composites to separate CO 2 from dilute streams. In parallel, but not necessarily in synergy, there have been studies focused on a few select ILs and DESs for electrochemical reduction of CO 2 , often diluting them with aqueous or non-aqueous solvents. The resulting electrode–electrolyte interfaces present a complex speciation for RCC. In this review, we describe how the ILs and DESs are tuned for RCC and specifically address the CO 2 chemisorption and electroreduction mechanisms. Critical bulk and interfacial properties of ILs and DESs are discussed in the context of RCC, and the potential of these electrolytes are presented through a techno-economic evaluation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental and modeling study of 1/ f noise in multilayer MoS2 and MoSe2 field-effect transistors

In field-effect transistors (FETs) with two-dimensional (2D) transition metal dichalcogenide channels, the dependence of field-effect mobility on atomic layer thickness has been studied and interpreted in terms of interface scattering and interlayer coupling resistance (Rint). However, a model for 1/f noise, such as in MoS2 and in MoSe2 FETs, for various contact metals and layer number thicknesses has not been reported. In this work, we have experimentally studied current–voltage and 1/f noise on MoS2 and MoSe2 FETs with source and drain contacts of high and low work function metals to understand both the mobility and the noise behavior. We have developed a noise model incorporating layer number dependent Hooge parameters and Rint. The noise and mobility models utilize screening lengths for charge, mobility, and Hooge parameter to describe the variation of these quantities with a layer number. Using our single model topology with appropriate fitting parameters for each material and each contact metal, the model captures the experimentally observed layer thickness dependence of the Hooge parameter. Our noise analysis is fully comprehensive and, hence, could be applied to any 2D layered systems.

Kwon, Jiseok↗
DOE OSTI2026

DIVA/DeviceEditor v6.1.2

DIVA is an end-to-end DNA design and construction management platform that streamlines how researchers design, build, and receive sequence-verified DNA constructs. Through a web-based BioCAD interface (DeviceEditor), researchers independently design DNA constructs and submit them to a centralized queue with a single action. Designs progress transparently through standardized states which allow researchers to track status and access finished constructs via a central DNA repository. Submitted designs are reviewed by dedicated staff for feasibility and optimization, reducing costly failures and improving downstream execution. Automated DNA assembly software optimizes construction strategies by reusing existing parts where possible and sourcing synthetic DNA only when needed. Standardized, sequence-agnostic assembly methods enable many independent constructs to be built in parallel using lab automation, dramatically increasing throughput. High-throughput next-generation sequencing is used to verify construct accuracy, with flexible platforms selected based on task requirements. Throughout the process, detailed success and failure data are captured and analyzed, enabling continuous improvement of assembly protocols. Compared to traditional, manual DNA construction workflows, DIVA offers higher scalability, transparency, reproducibility, and data-driven optimization.

Plahar, Hector [Lawrence Berkeley National Laborat↗
DOE OSTI2022

Volumetric energy density impact on mechanical properties of additively manufactured 718 Ni alloy

The 718 Ni-based superalloy has gained enormous attention in the additive manufacturing community for its great weldability, allowing for complex geometries to be formed, and its superb mechanical strength. Here we explore the influence of a wide range of major build parameters on microstructures and mechanical properties of 718 Ni alloys. Volumetric energy density appears to be one of the major parameters that captures the relationship between energy input and measured mechanical properties and microstructures of the alloys. A threshold energy density was identified below which the properties of the additively manufactured specimen degrade substantially. Finally, the influence of energy density on the quality of the built specimen simulated using Flow3D AM correlates well with the experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure, dynamics, and electrochemistry of choline chloride/ethylene glycol eutectic solvents at an electrode surface explored by molecular dynamics simulations

Choline chloride and ethylene glycol mixtures with 1:2, 1:4, and 1:6 molar ratios on the surfaces of graphite and gold electrodes were studied using classical molecular dynamics simulations. Both neutral and charged electrodes were considered. The liquid composition, solvation structure, molecular orientation, and dynamics at the electrode surface are significantly different from those of the bulk liquid. These properties strongly depend on the electrode material and charge density, whereas they are less sensitive to the overall solvent composition. The effect of the electrode on the composition, structure, and orientation of the liquid fades beyond ∼10 Å from the surface of the electrode. This distance corresponds to about two layers of the structured solvent, despite the fact that the layered structure extends to at least five layers or about 25 Å. However, the electrode influences solvent dynamics over a longer distance. The electrochemical properties of the eutectic solvent at both electrode surfaces were also studied. The simulations captured the experimental differential capacitance shapes for both electrode systems, although the magnitudes and exact shapes differ. The simulations further revealed that the solvent in the first solvation layer can both overscreen and underscreen the electrode charges depending on the electrode material and electrode potential.

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

Reactor wall effects in Si–Cl 2 –Ar atomic layer etching

This work complements our previous manuscript [J. Vac. Sci. Technol. A41, 062602 (2023)] where predictions from molecular dynamics (MD) simulations of silicon–chlorine–argon (Si–Cl 2 –Ar) atomic layer etching (ALE) are compared to experiments. When etch product distributions for atomic chlorine (Cl) and silicon chlorides were initially compared to optical emission spectroscopy (OES) signals, it appeared that there was a discrepancy between the MD predictions and experimental results at higher ion fluences. Experiments showed a relatively long period of nearly constant Cl-containing etch products released from the ion-bombarded surface (referred to as the “plateau”) but this effect was not observed in MD simulations. In this report, we demonstrate that the “plateau” observed in the OES signals is most likely due to the desorption of Cl-containing etch products from the walls of the reactor and subsequent adsorption on the Si substrate. Experiments varying the gas residence time in the chamber while keeping incoming gas concentrations and pressure constant support this interpretation. We also conducted experiments with an additional Ar-only flow in the chamber to reduce the concentration of Cl-containing species on the chamber walls. For both sets of flow modification experiments, we observe results consistent with the hypothesis that Cl-containing species desorbing from chamber walls are a significant cause of the observed discrepancy between MD predictions and experimental observations. If the measured OES signals are corrected for this “additional” source of Cl-containing species at the surface, the MD predictions and measured OES signals are in excellent agreement. This further supports the predictive capability of MD simulations to accurately capture the relevant physical and chemical processes in plasma-assisted ALE processes. We provide an order of magnitude estimate of the required density of Cl-containing species that would account for the additional etch products observed. Finally, we discuss the implications of this effect on ALE in plasma nanofabrication.

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

A software package for plasma facing component analysis and design: the Heat flux Engineering Analysis Toolkit (HEAT)

The engineering limits of plasma facing components (PFCs) constrain the allowable operational space of tokamaks. Poorly managed heat fluxes that push the PFCs beyond their limits not only degrade core plasma performance via elevated impurities, but can also result in PFC failure due to thermal stresses or melting. Simple axisymmetric assumptions fail to capture the complex interaction between 3D PFC geometry and 2D or 3D plasmas. This results in fusion systems that must either operate with increased risk or reduce PFC loads, potentially through lower core plasma performance, to maintain a nominal safety factor. High precision 3D heat flux predictions are necessary to accurately ascertain the state of a PFC given the evolution of the magnetic equilibrium. A new code, the Heat flux Engineering Analysis Toolkit (HEAT), has been developed to provide high precision 3D predictions and analysis for PFCs. HEAT couples many otherwise disparate computational tools together into a single open source python package. Magnetic equilibrium, engineering CAD, finite volume solvers, scrape off layer plasma physics, visualization, high performace computing, and more, are connected in a single web-based user interface. Linux users may use HEAT without any software prerequisites via an appImage. This manuscript introduces HEAT, discusses the software architecture, presents first HEAT results, and outlines physics modules in development.

divertor physics↗
DOE OSTI2021

Imaging, understanding, and control of nanoscale materials transformations

The development of liquid cells for transmission electron microscopy has enabled breakthroughs in our ability to follow nanoscale structural, morphological, or chemical changes during materials growth and applications. Time-resolved high-resolution imaging and chemical analysis through liquids opened the opportunity to capture nanoscale dynamic processes of materials, including reaction intermediates and the transformation pathways. In this article, a series of work is highlighted with topics ranging from liquid cell developments to in situ studies of nanocrystal growth and transformations, dendrite formation, and suppression of lithium dendrites through in situ characterization of the solid–electrolyte interphase chemistry. The understanding garnered is expected to accelerate the discovery of novel materials for applications in energy storage, catalysis, sensors, and other functional devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tuning spin excitations in magnetic films by confinement

Spin excitations of magnetic thin films are the founding element for novel transport concepts in spintronics, magnonics, and magnetic devices in general. While spin dynamics have been extensively studied in bulk materials, their behaviour in mesoscopic films is less known due to experimental limitations. Here, we employ Resonant Inelastic X-Ray Scattering to investigate the spin excitation spectrum in mesoscopic Fe films, from bulk-like down to 3 unit cells thick. In bulk-like samples, we find isotropic, dispersive ferromagnons consistent with the dispersion observed by neutron scattering in bulk single crystals. As the thickness is reduced, these ferromagnons survive and evolve anisotropically: renormalising to lower energies along the out-of-plane direction while retaining their dispersion in the in-plane direction. This thickness dependence is captured by simple Heisenberg model calculations accounting for the confinement in the out-of-plane direction through the loss of Fe bonds. Overall, our findings highlight the effects of mesoscopic scaling on spin dynamics and identify thickness as a knob for fine-tuning and controlling magnetic properties in films.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Spin Polarization Enhanced Ethanol Selectivity in Electrocatalytic CO 2 Reduction on the Paramagnetic CuO Surface

We report an electrochemical CO 2 reduction reaction catalyzed by a paramagnetic and conductive CuO/Cu interface with spins polarized by a moderate external magnetic field (MF) of similar to 800 gauss, achieving a similar to 30% increase in CO 2 -to-C 2+ Faradaic efficiency (FE) compared to that in the absence of the MF in a flow cell electrolyzer. At a current density of 400 mA/cm 2 , the CO 2 -to-C 2+ FE reached 86.7 ± 2.7% with 47.9 ± 1.4% cathodic energy efficiency (EE) in contrast to the CO 2 -to-C 2+ FE of 67.6% with 36.4% of EE in the absence of MF. Notably, ethanol production exhibits a much higher response to the MF (similar to 55.6% increase in FE) than ethylene (similar to 6.4% increase in FE) at 400 mA/cm 2 . In situ surface-enhanced Raman spectroscopy (SERS) captured magnetic-field-enhanced *CO coverage and ethanol-forming C 2 intermediates on CuO/Cu, providing direct spectroscopic evidence of spin-modulated pathway selection. Here, computational study suggests that the enhancement of ethanol selectivity is due to the reduced reaction kinetic barrier under MF, while the ethylene selectivity is less affected, mainly due to the insensitivity of the kinetic barriers under MF.

10 SYNTHETIC FUELS↗
DOE OSTI2024

DyG-DPCD: A Distributed Parallel Community Detection Algorithm for Large-Scale Dynamic Graphs

Dynamic (Temporal) graphs capture the valuable evolution of real-world systems, from the continuously evolving patterns of social interactions and genetic pathways to the dynamic fluctuations of economic forces. Detecting communities for such evolving networks poses unique challenges. Detecting and analyzing the evolution of communities within dynamic graphs unlocks valuable insights into the underlying structural and temporal patterns of real-world systems. However, the sheer volume of modern graph data and the inherent complexity of the temporal dimension pose significant challenges to scalable community detection algorithms. Addressing this gap, our work explores the limited landscape of scalable distributed-memory parallel methods specifically designed for dynamic network community detection. We propose a novel parallel algorithm, DyG-DPCD (Dynamic Graph Distributed Parallel Community Detection), to detect communities in dynamic networks using the Message Passing Interface (MPI) framework. We present a vertex-centric approach, allowing us to detect communities through local optimization. Furthermore, we enhance our baseline algorithm by incorporating three heuristics, which improve the algorithm’s performance significantly while maintaining the quality of the solutions. We demonstrate the efficiency of our algorithm by experimenting on several real-world large-scale networks with hundreds of millions of edges spanning diverse domains. Notably, DyG-DPCD achieves speedups between 25× and 30× for large networks that we experimented on using NERSC compute nodes. In conclusion, our algorithm outperforms the STINGER parallel re-agglomeration algorithm by 30×.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗