Engineering Papers⌕ Search

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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

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↗
DOE OSTI2023

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.↗
DOE OSTI2024

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↗
DOE OSTI2022

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↗
DOE OSTI2025

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↗
DOE OSTI2025

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↗
DOE OSTI2022

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↗
DOE OSTI2025

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↗
DOE OSTI2024

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↗
DOE OSTI2020

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↗
DOE OSTI2021

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↗
DOE OSTI2023

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↗
DOE OSTI2024

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↗
DOE OSTI2024

Enhanced Ferromagnetism in Atomically Thin Oxides Achieved by Interfacial Reconstruction

Discoveries of ferromagnetic materials with ultrathin thickness are of great importance for both fundamental science and technological applications. Transition metal oxides (TMOs) provide promising candidates in the context of next‐generation spintronics, despite the severe decay of ferromagnetism as the thickness reduces to the nanometer regime. Here, in this work, an efficient strategy to eliminate the magnetic dead layer in atomically thin oxides is presented, by using the epitaxial interface of 3 d and 5 d oxide monolayers that reconciles both strong exchange interaction and large uniaxial magnetic anisotropy. Combining multiple experimental methods, a ferromagnetic transition in an ultrathin oxide heterostructure comprised of only one La 0.2 Sr 0.8 MnO 3 monolayer sandwiched by SrIrO 3 monolayer (total thickness of three unit‐cells) is unambiguously demonstrated. Remarkably, a largely enhanced saturation magnetization (2 µ B Mn −1 ) and Curie temperature (80 K) are observed for the single manganite monolayer, as compared to previously reported ferromagnetic monolayer oxides. The results demonstrate a general strategy for creating robust ferromagnetism in ultrathin TMOs, potentially enabling novel oxide spin‐orbitronic devices.

2D ferromagnetism↗
DOE OSTI2020

Multi-material hydrodynamics with algebraic sharp interface capturing

A finite volume method for Eulerian multi-material hydrodynamics with sharp interface capturing is presented here. The pressure-temperature non-equilibrium multi-material system with finite-rate pressure relaxation in mixed-cells is considered here. This pressure closure facilitates material-property-dependent pressure relaxation, rather than instantaneous pressure equilibration, which in turn allows the use of unsplit high-order time-integrators. A modified tangent of hyperbola for interface capturing (THINC) method is used to reconstruct multi-material (> 2) interfaces, on three-dimensional unstructured meshes. A simple modification which extends the THINC reconstruction to interfaces between more than two materials is proposed. It is demonstrated that the modified THINC can capture multi-material interfaces within 2–4 tetrahedral cells. Since no geometric reconstructions are required by the THINC method, the presented multi-material method is algorithmically simple, and computationally efficient. Consistent reconstructions of conserved quantities at material interfaces ensure that conservation and closure laws are satisfied at the discrete level. Through a suite of test problems solved on unstructured meshes, it is demonstrated that the presented method is a promising candidate for accurate and efficient multi-material hydrodynamics computations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Ferroelectric Self-Polarization Controlled Magnetic Stratification and Magnetic Coupling in Ultrathin La 0.67 Sr 0.33 MnO 3 Films

Multiferroic oxide heterostructures consisting of ferromagnetic and ferroelectric components hold the promise for nonvolatile magnetic control via ferroelectric polarization, advantageous for the low-dissipation spintronics. Modern understanding of the magnetoelectric coupling in these systems involves structural, orbital, and magnetic reconstructions at interfaces. Previous works have long proposed polarization-dependent interfacial magnetic structures; however, direct evidence is still missing, which requires advanced characterization tools with near-atomic-scale spatial resolutions. Here, extensive polarized neutron reflectometry (PNR) studies have determined the magnetic depth profiles of PbZr 0.2 Ti 0.8 O 3 /La 0.67 Sr 0.33 MnO 3 (PZT/LSMO) bilayers with opposite self-polarizations. When the LSMO is 2–3 nm thick, the bilayers show two magnetic transitions on cooling. However, temperature-dependent magnetization is different below the lower-temperature transition for opposite polarizations. PNR finds that the LSMO splits into two magnetic sublayers, but the inter-sublayer magnetic couplings are of opposite signs for the two polarizations. Near-edge X-ray absorption spectroscopy further shows contrasts in both the Mn valences and the Mn–O bond anisotropy between the two polarizations. Finally, this work completes the puzzle for the magnetoelectric coupling model at the PZT/LSMO interface, showing a synergic interplay among multiple degrees of freedom toward emergent functionalities at complex oxide interfaces.

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

Electronic reconstruction at the polar (111)-oriented oxide interface

Atomically flat (111) interfaces between insulating perovskite oxides provide a landscape for new electronic phenomena. For example, the graphene-like coordination between interfacial metallic ion layer pairs can lead to topologically protected states [Xiao et al., Nat. Commun. 2, 596 (2011) and A. Rüegg and G. A. Fiete, Phys. Rev. B 84, 201103 (2011)]. The metallic ion/metal oxide bilayers that comprise the unit cell of the perovskite (111) heterostructures require the interface to be polar, generating an intrinsic polar discontinuity [Chakhalian et al., Nat. Mater. 11, 92 (2012)]. Here, we investigate epitaxial heterostructures of (111)-oriented LaAlO 3 /SrTiO 3 (LAO/STO). We find that during heterostructure growth, the LAO overlayer eliminates the structural reconstruction of the STO (111) surface with an electronic reconstruction, which determines the properties of the resulting two-dimensional conducting gas. This is confirmed by transport measurements, direct determination of the structure and atomic charge from coherent Bragg rod analysis, and theoretical calculations of electronic and structural characteristics. Interfacial behaviors of the kind discussed here may lead to new growth control parameters useful for electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Continuous Family of Conformal Field Theories and Exactly Marginal Operators

Does a conformal manifold imply the existence of exactly marginal operators? We answer this question affirmatively under the assumption that there exists a conformal interface with certain properties connecting nearby conformal field theories. We show that the exactly marginal operator that connects the conformal field theories can be reconstructed from the interface displacement operator. Our construction is model independent and based on the general principles of conformal symmetry.

conformal field theory↗