Results for “Diffuse interface methods” — Page 5

DOE OSTI2026

Probing and Tuning Spatiotemporal pH Evolution in Aqueous Zinc Ion Batteries

Aqueous zinc-ion batteries (AZIBs) offer a combination of safety and low cost, but they suffer from dynamic pH fluctuations that drive parasitic reactions and capacity fading. In this review, we report the current mechanistic understanding of the dynamic pH evolution across the electrolyte-electrode interface and its role in triggering interphase instability. We summarize advanced strategies to regulate pH, including electrolyte engineering, electrode surface modification, and interphase construction. We further discuss current methods for probing pH dynamics and outline potential future approaches for gaining deeper mechanistic and design insights.

Zheng, Xueli [SLAC]↗

DOE OSTI2025

Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

Sitaraman, Hariswaran↗

NASA NTRS2026

Thermal Inspection for the Assessment of Adhesively Bonded Metal Adherents

Certification of adhesively bonded structures is currently a challenge for aircraft manufacturers. The ability to certify bonds in primary structure can reduce dependence on fasteners and this will enable more efficient manufacturing. For example, drilling holes and installing fasteners in bonded joints can be a potentially significant bottleneck in airframe manufacturing. In addition, fasteners increase the airframe weight and can add stress concentration areas. Reducing fastener count can accelerate manufacturing and improve aircraft performance. NASA is currently investigating nondestructive evaluation (NDE) techniques to assess bond integrity. For example, bond thickness influences bond strength and therefore NDE techniques are being investigated to determine bond thickness. For bonded metal structures, there is a large difference in the thermal diffusivity between the 7075 aluminum alloy metal adherends and aircraft grade adhesive. Multi-layered thermal models show a large variation in the thermal response for bondline adhesive thicknesses that vary from 100 to 300 microns for adherend thicknesses of 0.163 cm. Experimental through-transmission thermal measurements reveal promise to quantitatively characterize the bondline thickness. Results were validated with X-ray computed tomography and optical microscopy measurements, and influence of porosity on the thermal model and measurements are investigated.

bondline thickness↗

NASA NTRS2026

Thermal Inspection for the Assessment of Adhesively Bonded Metal Adherends

Certification of adhesively bonded structures is currently a challenge for aircraft manufacturers. The ability to certify bonds in primary structure can reduce dependence on fasteners and this will enable more efficient manufacturing. For example, drilling holes and installing fasteners in bonded joints can be a potentially significant bottleneck in airframe manufacturing. In addition, fasteners increase the airframe weight and can add stress concentration areas. Reducing fastener count can accelerate manufacturing and improve aircraft performance. NASA is currently investigating nondestructive evaluation (NDE) techniques to assess bond integrity. For example, bond thickness influences bond strength and therefore NDE techniques are being investigated to determine bond thickness. For bonded metal structures, there is a large difference in the thermal diffusivity between the 7075 aluminum alloy metal adherends and aircraft grade adhesive. Multi-layered thermal models show a large variation in the thermal response for bondline adhesive thicknesses that vary from 100 to 300 microns for adherend thicknesses of 0.163 cm. Experimental through-transmission thermal measurements reveal promise to quantitatively characterize the bondline thickness. Results were validated with X-ray computed tomography and optical microscopy measurements, and influence of porosity on the thermal model and measurements are investigated.

X-ray computed tomography↗

NASA NTRSDate not supplied

Experimental and Numerical Investigation of the NASA High Efficiency Centrifugal Compressor Vaned Stage Geometry and Aerodynamic Performance

Since its inception in the early 2010s, the NASA High Efficiency Centrifugal Compressor (HECC) has been enigmatic for the propulsion research community: experimental data and numerical simulations of the stage have generally not aligned in their quantifications of performance metrics. Typically, the fault for these disagreements is assigned to the numerical simulations as the simulations are models, and models are inherently incomplete representations of the experiment. This “incompleteness” may manifest in assumptions regarding roughness or heat transfer, simplifications of the flow path (i.e., neglecting bleed flows), or the oft-scapegoated turbulence model. In the case of HECC, recent work showed unexpected discrepancies between the intended impeller geometry defined in the design report (termed Design-Intent) and the manufactured impeller used in the experimental campaigns (termed As-Manufactured). That work used numerical simulations to establish that the geometric differences between the Design-Intent and As-Manufactured impellers were significant enough to result in drastically different performance predictions for the HECC vaneless diffuser configuration. The Design-Intent impeller simulations over predicted the performance relative to the experiment, whereas the As-Manufactured simulations better represented the experimental data, both in terms of one-dimensional performance metrics and spanwise flow profiles. This effort expands on that work by examining in detail the geometry and aerodynamic performance of the HECC vaned diffuser configuration. Further differences between the Design-Intent and As-Manufactured geometries have been discovered in the vaned diffuser and exit guide vanes, and these differences are documented herein. The summations of the geometric differences for all of the components were used to create two numerical models of HECC vaned diffuser configuration: the Design-Intent simulations which are generated from the original geometry definitions given in the design report and the As-Manufactured simulations which are the best available representation of the manufactured compressor hardware used in the experimental test campaigns. In congruence with the earlier vaneless diffuser work, the numerical predictions of the Design-Intent choked mass flow rate, total pressure ratio, and efficiency were notably greater than that of the As-Manufactured simulations. To increase confidence in the experimental dataset, measurements from a recent test campaign conducted in 2024 are used to validate the original experimental data acquired from 2012 to 2014 with good repeatability overall, especially considering the passage of time and differences in the data acquisition systems between the test campaigns. Both numerical simulations were then extensively evaluated against the experimental data. The As-Manufactured simulations provided better estimates of the stage performance than the Design-Intent cases in terms of most performance metrics. Nonetheless, more detailed results still show opportunities for improvement. Despite a more accurate representation of the physical hardware, characterization of the impeller work input remains a challenge even for rigorously developed numerical models.

vaned diffuser↗

NASA NTRSDate not supplied

Experimental and Numerical Investigation of the NASA High Efficiency Centrifugal Compressor Vaned Stage Geometry and Aerodynamic Performance

Since its inception in the early 2010s, the NASA High Efficiency Centrifugal Compressor (HECC) has been enigmatic for the propulsion research community: experimental data and numerical simulations of the stage have generally not aligned in their quantifications of performance metrics. Typically, the fault for these disagreements is assigned to the numerical simulations as the simulations are models, and models are inherently incomplete representations of the experiment. This “incompleteness” may manifest in assumptions regarding roughness or heat transfer, simplifications of the flow path (i.e., neglecting bleed flows), or the oft-scapegoated turbulence model. In the case of HECC, recent work showed unexpected discrepancies between the intended impeller geometry defined in the design report (termed Design-Intent) and the manufactured impeller used in the experimental campaigns (termed As-Manufactured). That work used numerical simulations to establish that the geometric differences between the Design-Intent and As-Manufactured impellers were significant enough to result in drastically different performance predictions for the HECC vaneless diffuser configuration. The Design-Intent impeller simulations over predicted the performance relative to the experiment, whereas the As-Manufactured simulations better represented the experimental data, both in terms of one-dimensional performance metrics and spanwise flow profiles. This effort expands on that work by examining in detail the geometry and aerodynamic performance of the HECC vaned diffuser configuration. Further differences between the Design-Intent and As-Manufactured geometries have been discovered in the vaned diffuser and exit guide vanes, and these differences are documented herein. The summations of the geometric differences for all of the components were used to create two numerical models of HECC vaned diffuser configuration: the Design-Intent simulations which are generated from the original geometry definitions given in the design report and the As-Manufactured simulations which are the best available representation of the manufactured compressor hardware used in the experimental test campaigns. In congruence with the earlier vaneless diffuser work, the numerical predictions of the Design-Intent choked mass flow rate, total pressure ratio, and efficiency were notably greater than that of the As-Manufactured simulations. To increase confidence in the experimental dataset, measurements from a recent test campaign conducted in 2024 are used to validate the original experimental data acquired from 2012 to 2014 with good repeatability overall, especially considering the passage of time and differences in the data acquisition systems between the test campaigns. Both numerical simulations were then extensively evaluated against the experimental data. The As-Manufactured simulations provided better estimates of the stage performance than the Design-Intent cases in terms of most performance metrics. Nonetheless, more detailed results still show opportunities for improvement. Despite a more accurate representation of the physical hardware, characterization of the impeller work input remains a challenge even for rigorously developed numerical models.

centrifugal compressor↗

DOE OSTI2025

Active Fluids Form System-Spanning Filamentary Networks

Recent experimental realizations of liquid-liquid phase separation of active liquid crystals have offered an insight into the interaction between phase separation, ubiquitous in soft matter and biology, and chaotic active flows. Here, in this Letter, we use continuum theory to examine phase separation of an active liquid crystal and a passive fluid and report two new results. First, we provide an analytical derivation of the activity-induced suppression of the phase boundary of the coexistence region—a result first reported in simulations and experiments. We show that the shift in the critical point is a result of the balance between self-stirring active flows and phase-separating diffusive fluxes. Second, we show that this same balance is responsible for dramatically changing the morphology of the phase separated state, resulting in the emergence of a new mixed active phase consisting of a dynamical filamentous active network that invades the entire system area, trapping droplets of passive material. This structure exists even for very low volume fractions of active material. Our work provides an important step towards the goal of understanding how to use activity as a new handle for sculpting interfaces.

active nematics↗

NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

DOE OSTI2026

Visualization of the surface distribution of photocatalytic activity at the anatase/rutile TiO2 interface using X-ray photoelectron emission microscopy

We applied X-ray photoelectron emission microscopy (XPEEM) to visualize the surface distribution of photocatalytic activity at the anatase/rutile (A/R) interface of titanium dioxide. Acetic acid was adsorbed on the surface, and its decomposition and desorption reactions under ultraviolet light irradiation were monitored by observing C 1s spectra, enabling direct observation of the spatial distribution of photocatalytic activity. Complementary crystal structure information was obtained using spatially resolved low-energy electron diffraction, low-energy electron microscopy, X-ray absorption spectroscopy, and micro-focused Raman spectroscopy. Comparison with these structural data allowed for the evaluation of the photocatalytic activity as a function of the distance from the A/R interface. The activity is enhanced in the vicinity of the interface and gradually decreases toward both the anatase and rutile regions. These results demonstrate that XPEEM is an effective technique for spatially resolved mapping of photocatalytic activity, which has\\r\\npreviously been inaccessible using conventional characterization techniques.

25 ENERGY STORAGE↗

DOE OSTI2026

Elucidating the impact of stress states on grain boundary passivation in Ni-5Cr using the Rhines pack method

The oxidation behavior of a Ni–5Cr (at.%) alloy was evaluated at 420 ° C using the Rhines pack method and simultaneous tensile and compressive stress states via a miniature four-point bending fixture. At this moderate temperature, grain boundaries dominate mass transport and the resulting oxidation response. Oxidation produced approximately 1μ⁢m-wide protective Cr 2 O 3 films capping some grain boundaries and penetrative, intergranular Cr-rich oxides at other grain boundaries. Externally applied tensile and compressive stress during oxidation increased the prevalence of Cr 2 O 3 cap formation compared to no applied stress, with tensile stress resulting in more Cr 2 O 3 caps than compressive stress. The observed oxide cap morphology was similar across all test conditions. Regions under compressive stress showed an order of magnitude greater Cr depletion depth along the grain boundary. Chromium nitrides (CrN), likely from N contamination of the Rhines pack cell, were observed both at the oxide–metal interface and intergranularly within the alloy. Collectively these results demonstrate that applied stress promotes localized protective oxide cap formation over grain boundaries, with compressive stress additionally promoting deeper Cr depletion. Furthermore the experimental approach helped separate out the effects of stress on local oxide formation and grain boundary passivation.

4 point bend↗

DOE OSTI2025

Reflection and refraction of directrons at the interface

Reflection and refraction are ubiquitous phenomena with extensive applications, yet minimizing energy loss and information distortion during these processes remains a significant challenge. This study examines the behavior of structurally stable solitons, known as directrons, in nematic liquid crystals interacting with an interface where the director field orientation changes, despite identical physical properties, external potentials, and boundary anchoring in the two regions. During reflection and refraction, the directrons maintain nearly constant structure and velocity, ensuring energy conservation and information integrity. Microscopic analyses of the director field and macroscopic evaluations of effective potential are employed to elucidate the dependence of reflection and refraction probabilities on the directron’s incident angle and the orientation difference across the interface. The findings provide valuable insights into the dynamics of solitary waves in structured liquid crystal systems, offering significant implications for the development of tunable photonic devices, reconfigurable optical systems, and nanoscale material engineering.

Science & Technology - Other Topics↗

DOE OSTI2026

Interface-induced fast Li+ transport in mixed ionic–electronic conductors

Interfacial instability between lithium metal and solid-state electrolytes limits the performance of all-solid-state lithium metal batteries (ASSLBs), leading to parasitic reactions, non-uniform Li+ flux, and dendrite growth. Here, we develop a composite interlayer composed of the anti-perovskite Li2OHCl0.75Br0.25 (AP) and carbon nanotubes (CNTs) to enhance both interfacial stability and ionic transport. The AP–CNT interlayer exhibits enhanced Li+ conductivity arising from interfacial electron transfer from AP to CNTs, which generates a built-in electric field that facilitates Li+ migration. Lithium symmetric cells incorporating this interlayer achieve a high critical current density of 2.4 mA cm−2 at 55 °C. This design integrates chemical robustness with coupled ion–electron transport, offering a generalizable strategy for safe, dendrite-free, and high-performance ASSLBs.

He, Chenche↗

DOE OSTI2026

Generalizable Web User Interface for Scalable and Streamlined Deployment of Building Energy Management Systems in Small and Medium-Sized Commercial Buildings

Small and medium-sized commercial buildings (SMCBs) comprise 94% of US commercial buildings yet face significant barriers to implementing building energy management systems despite advances in smart device technology. Existing solutions present critical limitations: cloud-based API solutions simplify deployment but create vendor lock-in constraints; commercial integrated software solutions ensure compatibility via standardized protocols but require substantial cost and technical expertise; open-source IoT platforms offer cost-effective vendor independence but provide insufficient standardized protocol support for commercial building automation. This research presents a generalizable web user interface framework that bridges the gap between evolving smart device capabilities and lagging software infrastructure for SMCBs. The proposed system integrates VOLTTRON open-source middleware with an automated configuration converter that transforms unified specifications written in YAML, a human-readable data-serialization format, into system-specific files, streamlining manual setup processes. The vendor-agnostic architecture supports industry-standard protocols (BACnet and Modbus) and semantic building models while providing adaptive web interfaces that dynamically adjust to various building configurations. Demonstrations through simulation-based testing and a field deployment show automatic interface adaptation across heterogeneous HVAC systems and multizone monitoring. The automated configuration converter also substantially reduces labor-intensive setup.

Chung, Jihoon [ORNL] (ORCID:0000000184880815)↗

DOE OSTI2025

Improving Ionic Conformality Across Polymer Electrolyte|Electrode Interfaces

Maintaining uniform ionic transport at electrode|electrolyte interfaces, i.e., ionic conformality, remains challenging in polymer electrolyte (PE)-based solid-state batteries. Morphological conformality does not necessarily imply ionic conformality. In PEs, which typically consist of a mechanically supporting component and distinct ionically conductive components, the rearrangement or depletion of mobile ion-conductive domains at interfaces can disrupt ionic transport pathways. Such localized ionic depletion contributes to interfacial instability and capacity degradation in high-voltage lithium-metal batteries. Herein, an electrolyte design approach aimed at minimizing interfacial heterogeneities is demonstrated through compositional adjustments, characterized by spatially resolved structural and chemical X-ray techniques and NMR diffusometry to elucidate ion transport dynamics. This approach improves ionic conformality at electrode interfaces, enhancing cycling stability in Li||LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NMC811) coin and pouch cells cycled at high voltages. These results contribute to understanding interfacial behaviors in multiphase PEs and inform strategies for improving stability across solid-state battery interfaces.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Heterogeneous Solid Electrolyte Interphase Interactions Dictate Interface Instability in Sodium Metal Electrodes

Sodium (Na) metal batteries have attracted recent attention due to their low cost and high abundance of Na. However, the advancement of Na metal batteries is impeded due to key challenges such as dendrite growth, solid electrolyte interphase (SEI) fracture, and low Coulombic efficiency. This study examines the coupled electro‐chemo‐mechanical interactions governing the electrodeposition stability and morphological evolution at the Na/electrolyte interface. The SEI heterogeneities influence transport and reaction kinetics leading to the formation of current and stress hotspots during Na plating. Further, it is demonstrated that the heterogeneity‐induced Na metal evolution and its influence on the stress distribution critically affect the mechanical overpotential, contributing to a faster SEI failure. The analysis reveals three distinct failure mechanisms—mechanical, transport, and kinetic—that govern the onset of instabilities at the interface. Finally, a comprehensive comparative study of SEI failure in Na and lithium (Li) metal anodes illustrates that the electrochemical and mechanical characteristics of the SEI are crucial in tailoring the anode morphology and interface stability. This work delineates mechanistic stability regimes cognizant of the SEI attributes and underlying failure modes and offers important guidelines for the design of artificial SEI layers for stable Na metal electrodes.

25 ENERGY STORAGE↗

DOE OSTI2025

Geometric Delocalization in Two Dimensions

We demonstrate the existence of transient two-dimensional surfaces where a random-walking particle escapes to infinity in contrast to localization in standard flat two-dimensional space. We first prove that any rotationally symmetric two-dimensional membrane embedded in flat three-dimensional space cannot be transient. Then we formulate a criterion for the transience of a general asymmetric two-dimensional membrane. We use it to explicitly construct a class of transient two-dimensional manifolds with a nontrivial metric and height function but “zero average curvature,” which we dub “tablecloth manifolds.” The absence of the logarithmic infrared divergence of the Laplace-Beltrami operator in turn implies the absence of weak localization, nonexistence of bound states in shallow potentials, and breakdown of the Mermin-Wagner theorem and Kosterlitz-Thouless transition on the tablecloth manifolds, which may be realizable in both quantum simulators and corrugated two-dimensional materials.

Anderson localization↗

DOE OSTI2026

Performance of Heterostructural TaC/AlGaN Schottky Diodes Based on First Principles Electronic Structure Properties

Advances in ultra-wide bandgap materials, such as high Al-content AlxGa1-xN (AlGaN), are essential for next generation power electronics, but the requirement for lattice matched substrates is currently a significant obstacle. Recently, conductive TaC has emerged as a promising virtual substrate for AlGaN heteroepitaxy, with wurtzite (0001) AlxGa1-xN lattice-matched to rocksalt (111) TaC at x ~ 0.5. Thus, understanding and controlling the electronic properties of the TaC/AlGaN interface is key for developing technological applications based on TaC/AlGaN devices. Using density functional theory and electronic structure calculations, we here investigate TaC/Al0.5Ga0.5N interfaces, where we include explicit alloy models in the slab calculations. We predict the Schottky barrier height and the electric field discontinuity resulting from interface charges. Considering all possible combinations of (Ta, C) substrate termination, (Al/Ga, N) nucleation, and (Al/Ga, N) polarity, we construct a chemical potential phase diagram to identify the stable interfaces that can be accessed through variation of the synthesis conditions. The predicted interface electronic properties are implemented in device performance simulations to demonstrate a practical design for a strain-free, high-efficiency TaC/AlGaN Schottky diode with a low barrier height and without interface charges, underscoring the potential of TaC as a substrate for ultra-wide bandgap devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Bi‐Continuous W‐Rich Refractory High Entropy Alloy‐Cu Composite: Toward Material Innovation of Nuclear Reactor Coolant System

Abstract Refractory high‐entropy alloys (RHEAs) are considered promising candidate materials for next‐generation nuclear reactors due to their superior mechanical strength, irradiation resistance, and thermal stability at high temperatures. However, the significant positive heat of mixing between refractory alloying elements and Cu, commonly used in cooling systems, poses challenges in forming composite structures. This study addresses the issue using a liquid metal dealloying (LMD) process. A precursor alloy (WTaVTi) with a directional dendrite‐interdendrite structure is fabricated and reacted with molten Cu at 1200 °C for 96 h. This approach produced a RHEA‐Cu composite with a stable interface between RHEA (W 31.5 Ta 30.9 V 21.4 Ti 14.3 ) and Cu, featuring a spontaneously formed W‐rich interlayer that enhances interfacial bonding. The composite showed excellent irradiation resistance, with 30% less swelling under α‐ion irradiation than pure W. It also exhibited low thermal conductivity at room temperature, but reached ≈120 W m −1 ·K −1 at ≈650 °C, surpassing pure W. This temperature‐dependent rise in κ, with a positive gradient of +0.075 W m −1 ·K − 2 , is attributed to decreasing diffuse mismatch at elevated temperatures. The large‐scale reaction and stable microstructure achieved through LMD process highlight its industrial potential. This work offers a strategy for developing high‐performance materials by combining RHEA's radiation resistance with Cu's thermal conductivity for extreme environments.

Chemistry↗