Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sharp interface”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Exploring the effects of molecular beam epitaxy growth characteristics on the temperature performance of state-of-the-art terahertz quantum cascade lasers

This study conducts a comparative analysis, using non-equilibrium Green’s functions (NEGF), of two state-of-the-art two-well (TW) Terahertz Quantum Cascade Lasers (THz QCLs) supporting clean 3-level systems. The devices have nearly identical parameters and the NEGF calculations with an abrupt-interface roughness height of 0.12 nm predict a maximum operating temperature (T max ) of ~ 250 K for both devices. However, experimentally, one device reaches a T max of ~ 250 K and the other a T max of only ~ 134 K. Both devices were fabricated and measured under identical conditions in the same laboratory, with high quality processes as verified by reference devices. The main difference between the two devices is that they were grown in different MBE reactors. Our NEGF-based analysis considered all parameters related to MBE growth, including the maximum estimated variation in aluminum content, growth rate, doping density, background doping, and abrupt-interface roughness height. From our NEGF calculations it is evident that the sole parameter to which a drastic drop in T max could be attributed is the abrupt-interface roughness height. We can also learn from the simulations that both devices exhibit high-quality interfaces, with one having an abrupt-interface roughness height of approximately an atomic layer and the other approximately a monolayer. However, these small differences in interface sharpness are the cause of the large performance discrepancy. This underscores the sensitivity of device performance to interface roughness and emphasizes its strategic role in achieving higher operating temperatures for THz QCLs. We suggest Atom Probe Tomography (APT) as a path to analyze and measure the (graded)-interfaces roughness (IFR) parameters for THz QCLs, and subsequently as a design tool for higher performance THz QCLs, as was done for mid-IR QCLs. Our study not only addresses challenges faced by other groups in reproducing the record T max of ~ 250 K and ~ 261 K but also proposes a systematic pathway for further improving the temperature performance of THz QCLs beyond the state-of-the-art.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

2-D seismic wave propagation using the distributional finite-difference method: further developments and potential for global seismology

SUMMARY We present a time-domain distributional finite-difference scheme based on the Lebedev staggered grid for the numerical simulation of wave propagation in acoustic and elastic media. The central aspect of the proposed method is the representation of the stresses and displacements with different sets of B-splines functions organized according to the staggered grid. The distributional finite-difference approach allows domain-decomposition, heterogeneity of the medium, curvilinear mesh, anisotropy, non-conformal interfaces, discontinuous grid and fluid–solid interfaces. Numerical examples show that the proposed scheme is suitable to model wave propagation through the Earth, where sharp interfaces separate large, relatively homogeneous layers. A few domains or elements are sufficient to represent the Earth’s internal structure without relying on advanced meshing techniques. We compare seismograms obtained with the proposed scheme and the spectral element method, and we show that our approach offers superior accuracy, reduced memory usage, and comparable efficiency.

Geochemistry & Geophysics↗
DOE OSTI2024

Interface PINNs (I-PINNs): A physics-informed neural networks framework for interface problems

Here, we present a novel physics-informed neural networks (PINNs) framework for modeling interface problems, termed Interface PINNs (I-PINNs). I-PINNs uses different neural networks for any two subdomains separated by a sharp interface such that the neural networks differ only through their activation functions while the other parameters remain identical. The performance of I-PINNs, conventional PINNs, and other existing domain-decomposition PINNs methods such as extended PINNs (XPINNs) and multi-domain PINN (M-PINN) is compared through several one-dimensional, two-dimensional, and three-dimensional benchmark elliptic interface problems. The results demonstrate that I-PINNs provides a root-mean-square-error accuracy, at least two orders of magnitude better than conventional PINNs and XPINNs at approximately one-tenth of the computational cost of conventional PINNs and half the cost of XPINNs. Additionally, while I-PINNs and M-PINN provide comparable accuracies, M-PINN is found to be approximately 50% more expensive.

42 ENGINEERING↗
DOE OSTI2024

A mineral precipitation model based on the volume of fluid method

A novel volume of fluid method is presented for mineral precipitation coupled with fluid flow and reactive transport. The approach describes the fluid-solid interface as a smooth transitional region, which is designed to provide the same precipitation rate and viscous drag force as a sharp interface. Specifically, the governing equation of mineral precipitation is discretized by an upwind scheme, and a rigorous effective viscosity model is derived around the interface. The model is validated against analytical solutions for mineral precipitation in channel and ring-shaped structures. It also compares well with interface tracking simulations of advection-diffusion-reaction problems. Here, the methodology is finally employed to model mineral precipitation in fracture networks, which is challenging due to the low porosity and complex geometry. Compared to other approaches, the proposed model has a concise algorithm and contains no free parameters. In the modeling, only the pore space requires meshing, which improves the computational effciency especially for low-porosity media.

58 GEOSCIENCES↗
DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Method to Probe the Interfaces in La2−xSrxCuO4-LaSrAlO4-La2−xSrxCuO4 Trilayer Junctions

C-axis trilayer cuprate Josephson junctions are essential for basic science and digital circuit applications of high-temperature superconductors. We present a method for probing the interface perfection in La2−xSrxCuO4 (LSCO)-LaSrAlO4 (LSAO)-La2−xSrxCuO4 trilayer junctions. A series of LSCO-LSAO superlattices with atomically smooth surfaces and sharp interfaces were grown by the atomic-layer-by-layer molecular beam epitaxy (ALL-MBE) technique. We have systematically varied the thickness of LSCO and LSAO layers with monolayer precision. By studying the mutual inductance and electrical transport in these superlattices, we detect the non-superconducting (“dead”) layers at the interfaces and quantify their thicknesses. Our results indicate that two optimally doped LSCO monolayers just above and below the one monolayer LSAO barrier are no longer superconducting, rendering the actual barrier thickness of five monolayers. Next, we have shown that introducing a protective highly-overdoped LSCO layer reduces the thickness of dead layers by one or two monolayers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Giant Orbital Anisotropy with Strong Spin–Orbit Coupling Established at the Pseudomorphic Interface of the Co/Pd Superlattice

Orbital anisotropy at interfaces in magnetic heterostructures has been key to pioneering spin–orbit-related phenomena. However, modulating the interface’s electronic structure to make it abnormally asymmetric has been challenging because of lack of appropriate methods. Here, the authors report that low-energy proton irradiation achieves a strong level of inversion asymmetry and unusual strain at interfaces in [Co/Pd] superlattices through nondestructive, selective removal of oxygen from Co3O4/Pd superlattices during irradiation. Structural investigations corroborate that progressive reduction of Co3O4into Co establishes pseudomorphic growth with sharp interfaces and atypically large tensile stress. The normal component of orbital to spin magnetic moment at the interface is the largest among those observed in layered Co systems, which is associated with giant orbital anisotropy theoretically confirmed, and resulting very large interfacial magnetic anisotropy is observed. All results attribute not only to giant orbital anisotropy but to enhanced interfacial spin–orbit coupling owing to the pseudomorphic nature at the interface. They are strongly supported by the observation of reversal of polarity of temperature-dependent Anomalous Hall signal, a signature of Berry phase. This work suggests that establishing both giant orbital anisotropy and strong spin–orbit coupling at the interface is key to exploring spintronic devices with new functionalities.

42 ENGINEERING↗
DOE OSTI2021

Investigating the magnetic and atomic interface configuration for a model Fe/CrN bilayer system

A bilayer of iron on chromium nitride (Fe/CrN) is an interesting system for exchange biasing and sensing applications as the Néel temperature of CrN is 280 K and the Curie temperature of Fe is 1043 K. In this paper, we study the crystal and magnetic structures of the Fe/CrN interface at the atomic level. High quality epitaxial Fe/CrN bilayers prepared by molecular beam epitaxy grow in 001 orientation on MgO(001) substrates with uniform layer thicknesses and sharp interfaces. Our data reveal the epitaxial correlation between Fe and CrN crystals and their magnetic structures at the interface. The magnetic anisotropy directions of Fe and CrN are found parallel to [110] MgO . We studied the electronic and magnetic properties of the interface by performing the first-principles total-energy calculations. Here, we present a model that combines the crystal and magnetic structures of the Fe/CrN bilayer and fully explains all results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1997

Diffuse-Interface Methods in Fluid Mechanics

The authors review the development of diffuse-interface models of hydrodynamics and their application to a wide variety of interfacial phenomena. The authors discuss the issues involved in formulating diffuse-interface models for single-component and binary fluids. Recent applications and computations using these models are discussed in each case. Further, the authors address issues including sharp-interface analyses that relate these models to the classical free-boundary problem, related computational approaches to describe interfacial phenomena, and related approaches describing fully-miscible fluids.

SURFACE DIFFUSION↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs): An Efficient Physics-Informed Neural Networks Framework for Interface Problems

Here, we present an efficient physics-informed neural networks (PINNs) framework, termed Adaptive Interface-PINNs (AdaI-PINNs), to improve the modeling of interface problems with discontinuous coefficients and/or interfacial jumps. This framework is an enhanced version of its predecessor, Interface PINNs or I-PINNs (Sarma et al.; https://doi.org/10.1016/j.cma.2024.117135), which involves domain decomposition and assignment of different predefined activation functions to the neural networks in each subdomain across a sharp interface, while keeping all other parameters of the neural networks identical. In AdaI-PINNs, the activation functions vary solely in their slopes, which are trained along with the other parameters of the neural networks. This makes the AdaI-PINNs framework fully automated without requiring preset activation functions. Comparative studies on one-dimensional, two-dimensional, and three-dimensional benchmark elliptic interface problems reveal that AdaI-PINNs outperform I-PINNs, reducing computational costs by 2-6 times while producing similar or better accuracy.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

A finite element model of conduction, convection, and phase change near a solid/melt interface

Detailed understanding of heat transfer and fluid flow is required for many aerospace thermal systems. These systems often include phase change and operate over a range of accelerations or effective gravitational fields. An approach to analyzing such systems is presented which requires the simultaneous solution of the conservation laws of energy, momentum, and mass, as well as an equation of state. The variable property form of the governing equations are developed in two-dimensional Cartesian coordinates for a Newtonian fluid. A numerical procedure for solving the governing equations is presented and implemented in a computer program. The Galerkin form of the finite element method is used to solve the spatial variation of the field variables, along with the implicit Crank-Nicolson time marching algorithm. Quadratic Langrangian elements are used for the internal energy and the two components of velocity. Linear Lagrangian elements are used for the pressure. The location of the solid/liquid interface as well as the temperatures are determined form the calculated internal energy and pressure. This approach is quite general in that it can describe heat transfer without phase change, phase change with a sharp interface, and phase change without an interface. Analytical results from this model are compared to those of other researchers studying transient conduction, convection, and phase change and are found to be in good agreement. The numerical procedure presented requires significant computer resources, but this is not unusual when compared to similar studies by other researchers. Several methods are suggested to reduce the computational times.

Viterna, Larry A.↗
DOE OSTI2023

Time-reversal symmetry breaking superconductivity between twisted cuprate superconductors

Twisted interfaces between stacked van der Waals (vdW) cuprate crystals present a platform for engineering superconducting order parameters by adjusting stacking angles. Using a cryogenic assembly technique, we construct twisted vdW Josephson junctions (JJs) at atomically sharp interfaces between Bi 2 Sr 2 CaCu 2 O 8+x crystals, with quality approaching the limit set by intrinsic JJs. Near 45° twist angle, we observe fractional Shapiro steps and Fraunhofer patterns, consistent with the existence of two degenerate Josephson ground states related by time-reversal symmetry (TRS). By programming the JJ current bias sequence, we controllably break TRS to place the JJ into either of the two ground states, realizing reversible Josephson diodes without external magnetic fields. Furthermore, our results open a path to engineering topological devices at higher temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Atomic-scale manipulation of buried graphene–silicon carbide interface by local electric field

Abstract Precision of scanning tunneling microscopy (STM) enables control of matter at scales of single atoms. However, transition from atomic-scale manipulation strategies to practical devices encounters fundamental problems in protection of the designer structures formed atop the surface. In this context, STM manipulation of subsurface structures on technologically relevant materials is encouraging. Here, we propose a material platform and protocols for precise manipulation of a buried graphene interface. We show that an electric field from the STM tip reversibly controls breaking and restoring of covalent bonds between the graphene buffer layer and the SiC substrate. The process involves charge redistribution at the atomically sharp interface plane under the epitaxial graphene layer(s). This buried manipulation platform is laterally defined by unit cells from the corresponding (6×6) SiC moiré lattice of the epitaxial graphene. Local and reversible electric-field-induced patterning of graphene heterostructures from the bottom interface creates an alternative architecture concept for their applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Ensemble Kalman inversion of induced polarization data

SUMMARY This paper explores the applicability of ensemble Kalman inversion (EKI) with level-set parametrization for solving geophysical inverse problems. In particular, we focus on its extension to induced polarization (IP) data with uncertainty quantification. IP data may provide rich information on characteristics of geological materials due to its sensitivity to characteristics of the pore–grain interface. In many IP studies, different geological units are juxtaposed and the goal is to delineate these units and obtain estimates of unit properties with uncertainty bounds. Conventional inversion of IP data does not resolve well sharp interfaces and tends to reduce and smooth resistivity variations, while not readily providing uncertainty estimates. Recently, it has been shown for DC resistivity that EKI is an efficient solver for inverse problems which provides uncertainty quantification, and its combination with level set parametrization can delineate arbitrary interfaces well. In this contribution, we demonstrate the extension of EKI to IP data using a sequential approach, where the mean field obtained from DC resistivity inversion is used as input for a separate phase angle inversion. We illustrate our workflow using a series of synthetic and field examples. Variations with uncertainty bounds in both DC resistivity and phase angles are recovered by EKI, which provides useful information for hydrogeological site characterization. Although phase angles are less well-resolved than DC resistivity, partly due to their smaller range and higher percentage data errors, it complements DC resistivity for site characterization. Overall, EKI with level set parametrization provides a practical approach forward for efficient hydrogeophysical imaging under uncertainty.

Geochemistry & Geophysics↗
DOE OSTI2021

Micro-continuum approach for mineral precipitation

Abstract Rates and extents of mineral precipitation in porous media are difficult to predict, in part because laboratory experiments are problematic. It is similarly challenging to implement numerical methods that model this process due to the need to dynamically evolve the interface of solid material. We developed a multiphase solver that implements a micro-continuum simulation approach based on the Darcy–Brinkman–Stokes equation to study mineral precipitation. We used the volume-of-fluid technique in sharp interface implementation to capture the propagation of the solid mineral surface. Additionally, we utilize an adaptive mesh refinement method to improve the resolution of near interface simulation domain dynamically. The developed solver was validated against both analytical solution and Arbitrary Lagrangian–Eulerian approach to ensure its accuracy on simulating the propagation of the solid interface. The precipitation of barite (BaSO 4 ) was chosen as a model system to test the solver using variety of simulation parameters: different geometrical constraints, flow conditions, reaction rate and ion diffusion. The growth of a single barite crystal was simulated to demonstrate the solver’s capability to capture the crystal face specific directional growth.

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

Sequential multidimensional heteroepitaxy of chalcogen-sharing 3D ZnSe and 2D MoSe 2 with quasi van der Waals interface engineering

Two-dimensional (2D) materials are emerging as a promising platform for epitaxial growth, largely free from the constraints of lattice constant and thermal expansion coefficient mismatches. Among them, transition metal dichalcogenides (TMDs), known for their superior electrical properties, are ideal for ultrathin semiconductor applications. Their unique epitaxial characteristics enable seamless integration with 3D materials, facilitating the development of gate stacks and heterojunction devices. In this regard, developing a process for growing high-quality 3D epitaxial materials before and after the growth of 2D TMDs and understanding the 2D/3D interface are crucial. This study demonstrates the sequential growth of fully epitaxial ZnSe/MoSe 2 /ZnSe heterostructures using metal-organic chemical vapor deposition. ZnSe and MoSe 2 , sharing chalcogen elements, enable large-area quasi van der Waals epitaxy with sharp interfaces without intermediate phase. Multiscale analysis involving transmission electron microscopy and density functional theory calculation reveals lattice commensurability, van der Waals gaps, termination, and interfacial reconstruction. Understanding these interactions is crucial for advancing multidimensional integration of 2D and 3D materials.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Turbulence Modeling near Gas-Liquid Interfaces for Cryogenic Fluid Management

In cryogenic fluid management (CFM), gas-liquid interface turbulence transport (ITT) heavily influences crucial design parameters of spaceflight missions. ITT and phase change processes require length-scales and time-scales that are around 4-8 orders of magnitude smaller than the scales of CFM operations in spaceflights. A sharp interface method based parsimonious yet elegant turbulence modeling approach was developed from experimental observations and high-fidelity numerical studies. The technique was validated using flight-scale CFM tests.

cryogenic fluid management↗
DOE OSTI2023

A sharp immersed method for 2D flow-body interactions using the vorticity-velocity Navier-Stokes equations

Immersed methods discretize boundary conditions for complex geometries on background Cartesian grids. Here, this makes such methods especially suitable for two-way coupled flow-body problems, where the body mechanics are partially driven by hydrodynamic forces. However, for the vorticity-velocity form of the Navier-Stokes equations, existing immersed geometry discretizations for two-way coupled problems only achieve first order spatial accuracy near solid boundaries. Here we introduce a sharp-interface approach based on the immersed interface method to handle the one- and two-way coupling between an incompressible flow and one or more rigid bodies using the 2D vorticity-velocity Navier-Stokes equations. Our main contributions are three-fold. First, we develop and analyze a moving boundary treatment for sharp immersed methods that can be applied to PDEs with implicitly defined boundary conditions, such as those commonly imposed on the vorticity field. Second, we develop a two-way coupling methodology for the vorticity-velocity Navier-Stokes equations based on control-volume momentum balance that does not require the pressure field. Third, we show through extensive testing and validation that our resulting flow-body solver reaches second-order accuracy for most practical scenarios, and provides significant efficiency benefits compared to a representative first-order approach.

42 ENGINEERING↗