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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Phase-field approach to cellular blebbing

Bulges in the plasma membrane of cells known as blebs can form spontaneously in a wide range of biological processes, but what controls their shape and stability remains incompletely understood. Here, to address this we introduce a dual phase-field model with coupled order parameters representing the cell cortex and plasma membrane that can quantitatively model blebbing in three dimensions. Simulations and sharp-interface analyses reveal that, depending on whether blebbing occurs by detachment of the plasma membrane or rupture of the actin cortex, blebs can form discontinuously through a saddle-node bifurcation or continuously with increasing cortical tension. The model predictions are in good quantitative agreement with existing experimental data for laser-induced cortex rupture.

Biological and medical sciences↗
DOE OSTI2020

Topological Gaseous Plasmon Polariton in Realistic Plasma

Nontrivial topology in bulk matter has been linked with the existence of topologically protected interfacial states. In this work, we show that a gaseous plasmon polariton (GPP), an electromagnetic surface wave existing at the boundary of magnetized plasma and vacuum, has a topological origin that arises from the nontrivial topology of magnetized plasma. Because a gaseous plasma cannot sustain a sharp interface with discontinuous density, one must consider a gradual density falloff with scale length comparable to or longer than the wavelength of the wave. We show that the GPP may be found within a gapped spectrum in present-day laboratory devices, suggesting that platforms are currently available for experimental investigation of topological wave physics in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fragmentation of faceted crystalline wires

The reduction of surface energy is a well-known driving force for the fragmentation of nanowires. Fragmentation can limit the service life of those small low-dimensional building blocks or can be controlled to induce beneficial shape changes. For isotropic surface energy, fragmentation is analogous to the classic Rayleigh-Plateau instability of liquid jets. However, commonly synthesized crystalline nanowires have strongly anisotropic surface energies and exhibit facets coinciding with cusps in the Wulff plot. Depending on growth orientation, different fragmentation behaviors have been seen in nanowires. Using phase-field simulations, we show that fragmentation of faceted nanowires with cubic crystal symmetry may occur by a finite-amplitude nonlinear instability, as opposed to a Rayleigh-Plateau-like linear instability, depending on nanowire growth orientation. Here, we carry out a weakly nonlinear analysis based on sharp interface theory to characterize the faceted nanowire shape corresponding to a surface-energy saddle point. The analysis predicts that the minimum amplitude of a periodic shape perturbation to trigger fragmentation increases with cusp strength but decreases inversely proportionally to the perturbation wavelength for long wavelengths, in good quantitative agreement with phase-field predictions. The results provide the theoretical foundation to predict nanowire stability as a function of length and surface energy in diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

GeSn alloys with ∼21% Sn grown by effusion cell molecular beam epitaxy

We report the epitaxial growth of high quality GeSn alloys with Sn compositions up to 21.25 ± 1% by effusion-cell molecular beam epitaxy (MBE). Achieving such a high Sn content with effusion cells is particularly significant, since these sources typically impose stronger radiative heating on the substrate, which is a factor long considered a major obstacle to high-Sn GeSn epitaxy. Contrary to the prevailing assumption that high growth rates are good for suppressing Sn segregation, we demonstrate that carefully controlled low fluxes, combined with stable ultra-low substrate temperatures, enable significant Sn incorporation even under effusion-cell conditions. Structural analyses by high-resolution x-ray diffraction, reciprocal space mapping, atomic force microscopy, and transmission electron microscopy confirm sharp interfaces, high crystallinity, and smooth surfaces at high Sn contents. Furthermore, these findings establish the feasibility of MBE-grown high-Sn GeSn even under diffusion cells, thereby broadening opportunities for Si-compatible photonic and quantum devices.

Alloys↗
DOE OSTI2023

Virtual Flow Solver - Geophysics: A 3D Incompressible Navier-Stokes Solver

Virtual Flow Solver - Geophysics (VFS-Geophysics) is a three-dimensional (3D) incompressible Navier-Stokes solver based on the Curvilinear Immersed Boundary (CURVIB) method. The CURVIB is a sharp interface type of immersed boundary (IB) method that enables the simulation of fluid flows in the presence of geometrically complex moving bodies. The CURVIB method can be applied to wind/MHK turbine simulations and energy applications. VFS-Geophysics is the result of many years of research work by several graduate students, post-docs, and research associates that have been involved in the Computational Hydrodynamics and Biofluids Laboratory directed by Professor Fotis Sotiropoulos. The preparation of the present manual has been supported by the U.S. Department of Energy (DE-EE 0005482).

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Phase Field Modeling of Near-Wellbore Hydraulic Fracture Nucleation and Propagation

The goal of Task 2 of the project titled "Closing the loop between in situ stress complexity and EGS fracture complexity" is to model fracture nucleation and propagation in the near wellbore region and, in combination with the experimental study conducted as part of Task 4, investigate the relationship between the in situ stress and the near-wellbore fracture patterns. To this end, We have identified phase field as a promising modeling approach to model the nucleation and propagation process in the near wellbore region. In phase field, fractures are represented by a diffuse variable, i.e the damage. Compared to approaches that model fractures as sharp interfaces, phase field does not require any re-meshing or element insertion algorithm, which makes the treatment of complex fracture geometries straightforward, even with simple meshes.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Ginzburg--Landau functionals in the large-graph limit

Ginzburg–Landau (GL) functionals on graphs, which are relaxations of graph-cut functionals on graphs, have yielded a variety of insights in image segmentation and graph clustering. In this paper, we study large-graph limits of GL functionals by taking a functional-analytic view of graphs as nonlocal kernels. For a graph Wn with n nodes, the corresponding graph GL functional GL W n ϵ is an energy for functions on Wn. We minimize GL functionals on sequences of growing graphs that converge to functions called graphons. For such sequences of graphs, we show that the graph GL functional Γ-converges to a continuous and nonlocal functional that we call the graphon GL functional. We investigate the sharp-interface limits of the graph GL and graphon GL functionals, and we relate these limits to a nonlocal total-variation (TV) functional. We express the limiting GL functional in terms of Young measures and thereby obtain a probabilistic interpretation of the minimization problem in the large-graph limit. Finally, to develop intuition about graphon GL functionals, we determine the GL minimizer for several example families of graphons.

Zhang, Edith↗
DOE OSTI2022

Record Efficiency Multijunction Solar Cells with Strain-Balanced Quantum Well Superlattices

The absorption edge of a GaAs solar cell can be tuned to longer wavelengths with strain-balanced quantum wells, as an alternative to metamorphic epitaxy. For one-sun photovoltaic applications, an optimum bandgap for both single and multijunction solar cells is near the broad water absorption peak at 930 nm and can be accessed by a Ga0.9In0.1As alloy. This alloy is strained with respect to the GaAs substrate, and so only very thin ~10 nm layers can be grown without relaxation via dislocation formation. However, solar cells require significant thickness for complete absorption of the incoming solar spectrum, requiring many such layers. We have designed quantum well superlattices with up to 300 sets of carefully tuned strain-balanced GaInAs/GaAsP pairs by limiting composition modulation and material degradation. We have incorporated quantum well solar cells into multijunction devices to achieve record 32.9% two-junction and 39.5% three-junction solar cells under standard one-sun illumination. We will describe the physics of the quantum well devices and the growth conditions that lead to sharp interfaces and good material quality, as well as the characterization of the efficiencies.

GaAs↗
DOE OSTI2021

Multilayer Lateral Heterostructures of Van Der Waals Crystals with Sharp, Carrier–Transparent Interfaces

Abstract Research on engineered materials that integrate different 2D crystals has largely focused on two prototypical heterostructures: Vertical van der Waals stacks and lateral heterostructures of covalently stitched monolayers. Extending lateral integration to few layer or even multilayer van der Waals crystals could enable architectures that combine the superior light absorption and photonic properties of thicker crystals with close proximity to interfaces and efficient carrier separation within the layers, potentially benefiting applications such as photovoltaics. Here, the realization of multilayer heterstructures of the van der Waals semiconductors SnS and GeS with lateral interfaces spanning up to several hundred individual layers is demonstrated. Structural and chemical imaging identifies {110} interfaces that are perpendicular to the (001) layer plane and are laterally localized and sharp on a 10 nm scale across the entire thickness. Cathodoluminescence spectroscopy provides evidence for a facile transfer of electron‐hole pairs across the lateral interfaces, indicating covalent stitching with high electronic quality and a low density of recombination centers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

3D printed transparent ceramic YAG laser rods: Matching the core-clad refractive index

Yttrium Aluminum Garnet (YAG) solid state laser gain media rods with an active Neodymium-doped core and an optically-clear cladding region were additively manufactured via direct-ink-writing (DIW), followed by sintering and hot isostatic pressing to form fully dense optical ceramics. Lutetium and Gadolinium were chosen as optically-inert co-doping ions in the clad to match the increase in refractive index caused by the Neodymium in the core. Additionally, either 11.6% Lutetium or 3.8% Gadolinium can be used to match the index change from 2% neodymium; however, differences in diffusion distances across the core-clad interface lead to large fluctuations in index in that region. These index fluctuations can be minimized either by matching dopants with similar diffusion distances, or by implementing a gradual gradient in the doping profile, possible through DIW, rather than a sharp compositional interface attainable via more standard fabrication methods. This improvement in index homogeneity resulted in a 40% improvement in lasing performance compared with that of a core-clad rod fabricated with a sharp interface between the doped and undoped regions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sharp front tracking with geometric interface reconstruction

Here, this paper presents a novel sharp front-tracking method designed to address limitations in classical front-tracking approaches, specifically their reliance on smooth interpolation kernels and extended stencils for coupling the front and fluid mesh. In contrast, the proposed method employs exclusively sharp, localized interpolation and spreading kernels, restricting the coupling to the interfacial fluid cells–those containing the interface/front. This localized coupling is achieved by integrating a divergence-preserving velocity interpolation method with a piecewise parabolic interface calculation (PPIC) and a polyhedron intersection algorithm to compute the indicator function and local interface curvature. Surface tension is computed using the Continuum Surface Force (CSF) method, maintaining consistency with the sharp representation. Additionally, we propose an efficient local roughness smoothing implementation to account for surface mesh undulations, which is easily applicable to any triangulated surface mesh. Building on our previous work, the primary innovation of this study lies in the localization of the coupling for both the indicator function and surface tension calculations. By reducing the interface thickness on the fluid mesh to a single cell, as opposed to the 4–5 cell spans typical in classical methods, the proposed sharp front-tracking method achieves a highly localized and accurate representation of the interface. This sharper representation mitigates parasitic currents and improves force balancing, making it particularly suitable for scenarios where the interface plays a critical role, such as microfluidics, fluid-fluid interactions, and fluid-structure interactions. The proposed method is comprehensively validated and tested on canonical interfacial flow problems, including stationary and translating Laplace equilibria, oscillating droplets, and rising bubbles. The presented results demonstrate that the sharp front-tracking method significantly outperforms the classical approach in terms of accuracy, stability, and computational efficiency. Notably, parasitic currents are reduced by approximately two orders of magnitude and stable results are obtained for parameter ranges where classical front tracking fails to converge.

42 ENGINEERING↗
DOE OSTI2022

Temporal evolution of scalar modes in Richtmyer–Meshkov instability of inclined interface using high-speed PIV and PLIF measurements at 60 kHz

Abstract The current work presents simultaneous, high-speed measurements at 60,000 fields per second of velocity and mole fraction using particle image velocimetry (PIV) and planar laser induced acetone-fluorescence in a Richtmyer–Meshkov instability of an inclined interface (Atwood number, At = 0.22). Specifically, around 2 ms of temporal evolution of the vortex structures and their associated scalar modes immediately following the interface-reshock interaction is presented. Two initial interface conditions are discussed—(a) a sharp, inclined ‘single mode’ interface and (b) a ‘multi-mode’ interface where small perturbations are imposed on the single mode case. A 2D wavelet decomposition of the scalar flow field shows a highly intermittent distribution of small-scale variance throughout the interface even at late times. These are correlated strongly with the vortex structures and local turbulence intensity, where each small-scale scalar mode is sandwiched between two co-rotating vortex structures. This indicates that the interstitial regions between the vortices are significant hotspots of entrainment, which is then dispersed by the induced, counter-flow velocity fields. The multimode case demonstrates similar organization at large scales, while the scalar field is much more homogeneous at smaller scales. These observations highlight the importance of capturing the early time vortex evolution to accurately estimate any late time intermittency, especially where deposition of intense vorticity on sharp interfaces is present.

Engineering↗
DOE OSTI2020

Mechanistic and predictive understanding of needle litter decay in semi-arid mountain ecosystems experiencing unprecedented vegetation mortality. Final report

This research compared needle litter decay rates and underlying mechanisms between litter released from insect-killed trees to that of naturally senesced litter from healthy trees across climate variables with elevation and snowmelt rates as variables. Our research was designed to address questions of whether there is a significant impact to C and N cycling at the terrestrial-atmospheric interface in association with pine needle decomposition and the role of changes in soil moisture and temperature as a function of these processes. Field activities were primarily conducted within the Rocky Mountain Biological Laboratory’s (RMBL) East River watershed across a spectrum of ecoregions that spanned approximately 700 vertical meters. Forested catchments at these elevations in southwestern Colorado are threatened by heat-induced drought, altered snowmelt patterns, and invasive insect infestation that have all been associated with climate change. Needle litter from two members of the Pinaceae family experiencing large-scale insect infestations (spruce and lodgepole pine) were collected, redeployed and then sampled over the course of three years. We experimentally tested three tractable parameters that are known to impact litter decay and biogeochemical cycling and are relevant to bark beetle-induced tree mortality: (1) moisture content, (2) temperature by using elevation as a surrogate, and (3) changes to litter compositional chemistry that result from beetle infestation. We examined decay processes and monitored the migration of soluble and gaseous nitrogen and carbon species as well as microbiological signatures of decay such as biodegradable fractions, metal availability and the temporal succession of associated microorganisms at the terrestrial interface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗
DOE OSTI2023

Multipurpose active scanning probe cantilevers for near-field spectroscopy, scanning tunnel imaging, and atomic-resolution lithography

In this work, we report progress on developing a multipurpose scanning probe cantilever applying gallium nitride nanowires as the probe tip. Gallium nitride nanowires possess high potential as probes due to their straight profile, tunable electrical and optical properties, high Young’s Modulus, durability, and high-yield fabrication process. Their wide bandgap enables them to be pumped to emit ultraviolet pulses which can be used for optical imaging and spectroscopy. They can be doped during growth to be electrically conductive, and their sharp tips obtained during epitaxial growth enable confinement of a high electric field at tip–sample interface. Their sharp tips are obtained during fabrication by their epitaxial growth which eliminates the need for postprocess sharpening that is typically required for standard STM tips. We present results of using gallium nitride nanowires for scanning tunnel microscopy applications of atomic-resolution imaging and lithography, and atomic force microscopy applications of imaging and lithography in vacuum and atmospheric environments.

Engineering↗
DOE OSTI2024

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗
DOE OSTI2024

A Thermo‐Flow‐Mechanics‐Fracture Model Coupling a Phase‐Field Interface Approach and Thermo‐Fluid‐Structure Interaction

This work proposes a novel approach for coupling non-isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid-filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface-capturing phase-field fracture method and interface-tracking thermo-fluid-structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase-field approach to represent fractures and reconstruct the geometry to frame a thermo-fluid-structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase-field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid-filled fracture with the phase-field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase-field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.

fracture↗