Modeling Si/SiGe quantum dot variability induced by interface disorder reconstructed from multiperspective microscopy
this is a poster that covers a recent publication that was passed through IR - should look familiar to the DC.
Fabian pena, Luis↗
SEARCH · Engineering Papers
Results for “Interface reconstruction”
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.
Modeling Si/SiGe quantum dot variability induced by interface disorder reconstructed from multiperspective microscopy
10 minute talk on fundamental research of a type in academic literature
Bussmann, Ezra [Sandia National Laboratories (SNL-↗
Revisiting W–ZrC interfaces: A first principles study
We investigate W–ZrC interfaces using first-principles calculations based on the density functional theory. There have been theoretical and experimental studies exploring W–ZrC interfaces, however, the debate regarding the most stable interface continues to persist. In this study, we systematically simulated various W–ZrC interfaces merging W and ZrC surfaces with different orientations. Subsequently, we evaluated their stabilities and explained the corresponding stabilities in terms of the nature of bonding and charge-transfer processes at the interface. We find ZrC(111)–W(110) is the most stable interface with higher adhesive energy than the other interfaces. The additional stability associated with the ZrC(111)–W(110) results from significant interface reconstruction. Three layers of W and ZrC adjacent to the interface are involved in the charge-transfer process leading to stronger ionic bonds in ZrC(111)–W(110) as compared to the other potential candidate: ZrC(100)–W(100). The C and W atoms are found to be displaced from their symmetric position during the reconstruction process at the interface to facilitate stronger bonds with shorter W–C and W–Zr bonds in ZrC(111)–W(110) as compared to ZrC(100)–W(100). This leads to stronger covalent bonds in ZrC(111)–W(110) than that in ZrC(100)–W(100). Therefore, we conclude that the stronger covalent and ionic forces in ZrC(111)–W(110) than those in ZrC(100)–W(100) are responsible for making ZrC(111)–W(110) to be the most stable interface. In conclusion, this study addresses the long-standing question of the most stable W–ZrC interface and derives a number of implications for other W-transition metal carbide interfaces which are potential candidates for improving the mechanical properties of plasma facing materials.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Phosphorus diffusion and deactivation during SiGe oxidation
Dopant profiles near the semiconductor–oxide interface are critical for microelectronic device performance. As the incorporation of Si 1−x Ge x into transistors continues to increase, it is necessary to understand the behavior of dopants in Si 1−x Ge x . In this paper, the diffusion and electrical activation of phosphorus within a strained, single-crystal Si 0.7 Ge 0.3 layer on Si during oxidation are reported. Both layers were uniformly doped, in situ, with an average phosphorus concentration of 4 × 10 19 atoms/cm 3 . After high-temperature oxidation, secondary ion mass spectrometry measurements revealed that the bulk of the phosphorus diffuses out of only the SiGe layer and segregates at the oxidizing SiGe–SiO 2 interface. Hall effect measurements corroborate the observed phosphorus loss and show that the phosphorus diffusing to the oxidizing interface is electrically inactive. Through density functional theory (DFT) calculations, it is shown that phosphorus interstitials prefer sites near the SiGe–SiO 2 interface. Finally, based on a combination of experimental data and DFT calculations, we propose that the phosphorus atoms are displaced from their lattice sites by Ge interstitials that are generated during SiGe oxidation. The phosphorus atoms then migrate toward the SiGe–SiO 2 interface through a novel mechanism of hopping between Ge sites as P–Ge split interstitials. Once they reach the interface, they are electrically inactive, potentially in the form of interstitial clusters or as part of the reconstructed interface or oxide.
Thornton, Chappel S.↗
Interface-aware Methods for Diffusion Equations [Slides]
In this talk, we will present interface-aware methods for diffusion equation. The focus is on the necessity of interface-aware, mixed cells, interface reconstruction, discontinuity of material diffusion coefficients, order of accuracy in time, property for large time steps, and simpler interface-aware methods. Numerical examples are provided to demonstrate the features on these aspects of the methods
97 MATHEMATICS AND COMPUTING↗
Building interface bonding and shield for stable Li-rich Mn-based oxide cathode
Implementation of Li-rich Mn-based oxide cathode with high-energy-density has been restrained by capacity/voltage degradation that results from irreversible lattice oxygen loss and structure rearrangements. To resolve these challenges, in this work, Li 1.2 Mn 0.54 Ni 0.13 Co 0.13 O 2 encapsulated by amorphous Co x B (CB-LRM) is rationally designed via autocatalytic plating for highly reversible cationic/anionic hybrid cathode material. Band coherency is ingeniously evoked by interface-reconstruction between bulk structure and amorphous coating layer, which lower the energy of O $2p$ states. Furthermore, this is associated with strengthened orbital hybridization of O $2p$-Mn $3d$ and increased formation energy of oxygen vacancy, which mitigates the lattice oxygen loss considerably. Additionally, interface shielding effects that protect electrode against electrolyte corrosion and the reduction of surface oxygen are also present with fully coverage of amorphous Co x B coating layer. As a result, the as-designed CBLRM cathode exhibits excellent cycle stability after 100 loops with only 0.154% per cycle capacity fade and improved voltage degradation. Given this, this work provides a potential avenue for rational design of lattice oxygen-based electrode materials with high-energy-density.
25 ENERGY STORAGE↗
Multi-material ALE remap with interface sharpening using high-order matrix-free finite element methods
The arbitrary Lagrangian-Eulerian (ALE) technique involves remapping field quantities from a Lagrangian mesh to an optimized mesh in a conservative, accurate and bounds-preserving manner. For methods based on arbitrary order finite elements, as described in a reference, material volume fractions are advected in pseudo-time using flux-corrected transport (FCT) without any form of interface reconstruction. In practice, this can lead to excessive propagation of small volume fractions throughout the domain. In addition, this method requires assembly of a global advection matrix to compute the bounds-preserving low-order FCT solution. In this work, we introduce a new approach for ALE remap using a high-order matrix-free technique which incorporates a flux modification to sharpen material interfaces in a conservative manner. Our approach begins with computing a bounds-preserving low-order solution to the ALE remap equations at the element level. We then compute a sharp interface solution (not guaranteed to be bounds-preserving) which comes from solving an augmented version of the ALE remap equations with a conservative flux modification which acts to sharpen material volume fractions based on their gradients and transport directions. Using the sharp interface solution, we make global corrections to the bounds-preserving solution while maintaining preservation of bounds. By blending with the sharpened solution at the global level we are able to globally conserve mass without hindering the remap pseudo-time step. This new interface-aware ALE remap method is based entirely on partial assembly techniques where globally assembled matrix operators are no longer needed, resulting in a globally matrix-free FCT method for multi-material, multi-field ALE remap with high performance on GPU architectures. We present results of our new remap method on 1D, 2D and 3D benchmarks and describe the algorithmic tailoring for GPU architectures that was developed.
Vargas, Arturo [Lawrence Livermore National Labora↗
A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer
We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.
97 MATHEMATICS AND COMPUTING↗
Benchmark Solutions for Radiation Transport in Stochastic Media with Inhomogeneous Material Statistics
Accurately solving implicit Monte Carlo (IMC) thermal photon transport problems with mixed material cells is important in realistic applications. The production IMC package at LLNL treats mixed material cells arising from ALE remap and hydrodynamics using the same approximate model. The new Imp IMC thermal photon transport package currently under development has both a material interface reconstruction (MIR) algorithm and a Levermore-Pomraning (LP) stochastic medium algorithm for treating mixed material cells. Existing stochastic medium algorithms for treating mixed material cells in IMC lack a complete theoretical basis. The IMC LP algorithm implementation has been demonstrated to reproduce published deterministic LP solutions for the particular case of spatially homogeneous material statistics. Realistic simulations will include spatially inhomogeneous material statistics (material mean chord lengths). In a previous investigation, the LP-model for transport in binary stochastic media in rod geometry was generalized to accommodate spatially varying material chord lengths, i.e., the mixing statistics were allowed to be nonhomogeneous. Analytical solutions were obtained and used to produce a verifi cation suite for the Imp IMC Levermore-Pomraning implementation for different spatial variations of the chord lengths. However, the accuracy of the LP model when the mixing statistics are nonhomogeneous has not been assessed and leaves open the question of whether local accuracy is improved or further degraded when chord lengths are not uniform. This shortcoming is rectifi ed here by developing benchmark analytic solutions for transport in binary Markovian stochastic mixtures in rod geometry with nonhomogeneous mixing statistics, using spatially varying chord lengths considered in the previous investigation based on the LP model. Methods for sampling a nonhomogeneous Poisson process (NHPP) are first described and used to construct individual realizations of the binary mixtures in rod geometry. Analytic solutions are then obtained for the forward and backward directed fluxes on a given realization, now viewed as a deterministic medium with alternating layers of the two materials with known interface locations. Finally, material averaged scalar fluxes are obtained using these sampling schemes with spatially linear and quadratic chord lengths and used to assess the accuracy of the previously obtained LP-model results.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
A scalable compressible volume of fluid solver using a stratified flow model
Here, in this paper, a Volume of Fluid method which uses a Stratified Flow model for flux calculation is proposed to solve a compressible multiphase flow in a way that has high parallel efficiency. The solver is part of an open-source computing tool to study fuel entrainment and combustion in the application of paraffin-based hybrid rocket motors. This paper focuses on the compressible Volume of Fluid solver for the liquid and gaseous phases in this environment. This solver uses perfect gas and stiffened gas models for equation of state, with the ability to easily add other models as desired. Additionally a comparison of a Riemann solver versus an AUSM+up scheme is conducted, showing improvement with the later in the solution of multiphase flows. This solver successfully reproduces shock tube cases with the expected accuracy, as well as simulating two-dimensional shear and gravity-driven flows. Results also show that the code is able to utilize a Stratified Flow model to evaluate shear flow without the need for interface reconstruction or gradient calculations. Finally, it is shown that the solver has near-ideal scaling under strong scaling tests as well as good performance in static scaling, giving improved performance over current options and implying future advancement of high-performance computing efficiency for multiphase flow solvers.
97 MATHEMATICS AND COMPUTING↗
A Geometric Volume of Fluid-Based Multiphase Flow Solver Extension to the Reacting Flow Solver, PeleLM
A new algorithm is presented to simulate multiphase flows with surface tension in a pathway for spray combustion simulation. The algorithm combines capabilities from two open-source packages, including the interface reconstruction library (IRL), a library of computational geometry routines to enable the volume of fluid (VOF) method, and PeleLM, a solver for the reacting Navier-Stokes equations. Additionally, surface tension is implemented using the continuum surface force (CSF) model with an improved height function technique in the volume fraction field. Spurious errors in volume fraction arising from our combined strategy are corrected through a topology-based method that improves both numerical stability and accuracy. Multiple validation simulations are conducted, including (i) translations and rotations of Zalesak's disk, (ii) a stationary circular droplet with surface tension, (iii) an oscillating elliptical droplet, and (iv) three-dimensional deformation of a spherical droplet. Results indicate that the combined scheme retains the favorable properties of each of the component algorithms.
42 ENGINEERING↗
Tracking the dynamics of catalytic Pt/CeO2 active sites during water-gas-shift reaction
Abstract Understanding the atomistic structure of the active site during catalytic reactions is of paramount importance in both fundamental studies and practical applications, but such studies are challenging due to the complexity of heterogeneous systems. Here, we use Pt/CeO 2 as an example to study the dynamic nature of active sites during the water-gas-shift reaction (WGSR) by combining multiple in situ characterization tools. We show that the different concentrations of interfacial Pt δ+ – O – Ce 4+ moieties at Pt/CeO 2 interfaces are responsible for the rank of catalytic performance of Pt/CeO 2 catalysts: Pt/CeO 2 -rod > Pt/CeO 2 -cube > Pt/CeO 2 -oct. For all the catalysts, metallic Pt is formed during the WGSR, leading to the transformation of the active sites to Pt 0 – O v – Ce 3+ and interface reconstruction. These findings shed light on the nature of the active site for the WGSR on Pt/CeO 2 and highlight the importance of combining complementary in situ techniques for establishing structure-performance relationships.
36 MATERIALS SCIENCE↗
Development and Validation of a Two-Phase Thermal-Hydraulic CFD Code NEK-2P
A project is underway to develop, verify and validate an advanced two-phase flow modeling capability for the highly-scalable, high-performance Computational Fluid Dynamics (CFD) code NEK5000. The goal of this work is to verify and validate the two-phase version of the NEK5000 code, named NEK-2P, to simulate the two-phase flow and heat transfer phenomena that occur in a Boiling Water Reactor (BWR) fuel bundle under various operating conditions. The NEK-2P two-phase flow models follow the approach used for the Extended Boiling Framework (EBF) previously developed at Argonne but include more fundamental physical models of boiling phenomena and advanced numerical algorithms for improved computational accuracy, robustness, and computational speed. The development of the NEK-2P two-phase solver and the implementation of the Extended Boiling Framework two-phase models were initially supported by Argonne National Laboratory (Argonne) through a Laboratory Directed Research and Development (LDRD) project during FY2014-2016. The development and validation of the two-phase models through analyses of selected two-phase boiling flow experiments was supported by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY2017-2020. This report focuses on verification and validation of the water-steam boiling model NEK-2P Two-Phase, CFD code. The NEK-2P was validated with Nuclear Power Engineering Corporation (NUPEC) Pressurized Water Reactor (PWR) Sub-channel and Bundle Test (PSBT) void distribution benchmark. Three different simulations were performed and analyzed for various operating conditions such as wall-heat flux and sub-cooled inlet temperatures. Reasonably good agreement with measured data was obtained in predicting the measured void distributions. Simulations were performed for Virginia Tech. (VT) 3x3 rod bundle geometry with and without spacers. The preliminary results were presented for Simplified Spacer Grid (SSG). In addition, the implementation of interface reconstruction model was tested with one of the Becker benchmark Critical Heat Flux (CHF) experiments.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Numerical Diffusion (Mixing) of Material in Numerical Simulations of Hydrodynamics [Slides]
It is often assumed that a material interface between two materials is spread over a few numerical cells in numerical simulations for hydrodynamics. Also, we have the impression that higher order methods introduce less numerical diffusion (mixing) of material. As we know one of the purposes of adaptive mesh refinement (AMR) is to resolve interfaces between materials, but we would like to know how effective AMR is to reduce numerical diffusion of material. We will present our investigation about numerical diffusion (mixing) of material in xRage. The result of the investigation indicates that the assumptions mentioned above are not always valid. In this talk, we will also demonstrate the effectiveness of numerical techniques to reduce numerical diffusion of material, including contact discontinuity steepening, isotropic interface steepening, max interface steepening, material interface reconstruct.
36 MATERIALS SCIENCE↗
Brazing of Stainless Steels to Yttria Stabilized Zirconia (YSZ) Using Silver -Base Brazes
Three silver-base brazes containing either noble metal palladium (Palcusil-10 and Palcusil-15) or active metal titanium (Ticusil) were evaluated for high-temperature oxidation resistance, and their effectiveness in joining yttria stabilized zirconia (YSZ) to a corrosion-resistant ferritic stainless steel. Thermogravimetric analysis (TGA), and optical- and scanning electron microscopy (SEM) coupled with energy dispersive spectrometry (EDS) were used to evaluate the braze oxidation behavior and the structure and chemistry of the YSZ/braze/steel joints. The effect of the braze type and processing conditions on the interfacial microstructure and composition of the joint regions is discussed with reference to the chemical changes that occur at the interface. It was found that chemical interdiffusion of the constituents of YSZ, steel and the brazes led to compositional changes and/or interface reconstruction, and metallurgically sound joints.
Singh, Mrityunjay↗
Robust Interfacial Effect in Multi-interface Environment through Hybrid Reconstruction Chemistry for Enhanced Energy Storage
Electrochemical-oxidation-driven reconstruction has emerged as an efficient approach for developing advanced materials, but the reconstructed microstructure still faces challenges including inferior conductivity, unsatisfying intrinsic activity, and active-species dissolution. In this report we present hybrid reconstruction chemistry that synergistically couples electrochemical oxidation with electrochemical polymerization (EOEP) to overcome these constraints. During the EOEP process, the metal hydroxides undergo rapid reconstruction and dynamically couple with polypyrrole (PPy), resulting in an interface-enriched microenvironment. We observe that the interaction between PPy and the reconstructed metal center (i.e., Mn > Ni, Co) is strongly correlated. Theoretical calculation results demonstrate that the strong interaction between Mn sites and PPy breaks the intrinsic limitation of MnO 2 , rendering MnO 2 with a metallic property for fast charge transfer and enhancing the ion-adsorption dynamics. Operando Raman measurement confirms the promise of EOEP-treated Mn(OH) 2 (E-MO/PPy) to stably work under a 1.2 V potential window. The tailored E-MO/PPy exhibits a high capacitance of 296 F g –1 at a large current density of 100 A g –1 . Here our strategy presents breakthroughs in upgrading the electrochemical reconstruction technique, which enables both activity and kinetics engineering of electrode materials for better performance in energy-related fields.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Turbulent Mixing of Multiphase Flow
Thus we conduct numerical simulations of multiphase fluids stirred by two-dimensional turbulence to assess the possibility of self-similar drop size distribution in turbulence. In our turbulence simulations, we also explore the non-diffusive limit, where molecular mobility for the interface is vanishing. Special care is needed to transport the non-diffusive interface. Numerically, we use the particle level set method to evolve the interface. Instead of using the usual methods to calculate the surface tension force from the level set function, we reconstruct the interface based on phase- field modeling, and calculate the continuum surface tension forcing from the reconstructed interface.
Young, Y.-N.↗
Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods
A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.
97 MATHEMATICS AND COMPUTING↗