A high-order discontinuous Galerkin method for compressible interfacial flows with consistent and conservative Phase Fields
Not Available
White, William J. (ORCID:0000000322164500)↗
SEARCH · Search NASA
Results for “Interfacial flows”
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.
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↗
Atomic Dynamics of Multi‐Interfacial Migration and Transformations
Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Counter-Current Flow Limitation Studies in Complex Geometries Utilizing Interface Capturing Simulations Coupled with PID Flow Rate Controller
In nuclear thermal-hydraulic studies, counter-current flow limitation (CCFL) typically refers to steam rising at a fast rate such that it prevents coolant from draining down within a confined channel. CCFL is a crucial issue in nuclear reactor safety analysis. This study investigates CCFL in debris bed channels using high-resolution interface-capturing simulations. A novel proportional-integral-derivative flow rate controller is developed to efficiently achieve the CCFL conditions. Verification studies confirm that CCFL occurs under the same conditions with or without the controller, demonstrating that PID control ensures accurate prediction. Three debris bed channel geometries were examined: a cylindrical channel, a channel with small obstacles, and a channel with large obstacles. Results show that obstacles significantly impact flow behavior, interfacial shear, wall shear, and pressure gradients required for CCFL. Furthermore, the comparison with experimental data confirmed that simulations incorporating geometric complexities align more closely with experimental CCFL conditions. A pressure gradient correlation was also developed for CCFL prediction.
Counter-current flow limitation↗
Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction
A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.
Kittell, David E. [Sandia National Laboratories (S↗
(2+δ)-dimensional theory of the electromechanics of lipid membranes. III. Constitutive models
Explore the source record for details and available documents.
bending↗
Asymmetric fluctuations and self-folding of active interfaces
We study the structure and dynamics of the interface separating a passive fluid from a microtubule-based active fluid. Turbulent-like active flows power giant interfacial fluctuations, which exhibit pronounced asymmetry between regions of positive and negative curvature. Experiments, numerical simulations, and theoretical arguments reveal how the interface breaks up the spatial symmetry of the fundamental bend instability to generate local vortical flows that lead to asymmetric interface fluctuations. The magnitude of interface deformations increases with activity: In the high activity limit, the interface self-folds invaginating passive droplets and generating a foam-like phase, where active fluid is perforated with passive droplets. These results demonstrate how active stresses control the structure, dynamics, and break-up of soft, deformable, and reconfigurable liquid–liquid interfaces.
active fluid↗
Mechanistic mass transfer in hollow fiber membrane solvent extraction for bio-based isobutanol
Membrane solvent extraction (MSE) has emerged as a promising method for selectively recovering bioproducts from complex aqueous streams. Bio-isobutanol, a next-generation feedstock for biofuel, remains challenging to recover because of its low concentration and the presence of inhibitory substances. This study explores the potential of hollow fiber (HF) MSE for bio-isobutanol recovery and systematically examines the coupled effects of fiber packing, shell-side flow dynamics, and aqueous chemistry on performance. A resistance-in-series model is applied to understand mass transfer in the HF MSE modules, quantify local resistances, and validate overall performance. The results show that increasing the fiber packing provides a larger interfacial area but induces poor flow distribution and channeling, hindering effective isobutanol transport. Meanwhile, increasing the shell-side velocity improves isobutanol recovery due to reductions in the boundary layer thickness. The presence of salts, added to mimic fermentation broth, increases the partition coefficient through salting-out effects, further improving isobutanol flux. A modified correlation for the shell-side mass transfer coefficient (k s,ϕ+v ), integrating geometric and hydrodynamic effects, was developed and validated. The proposed model achieves highly predictive accuracy (r 2 = 0.9808) across a wide range of conditions, outperforming previous models. The findings provide mechanistic insight into the interaction of geometric packing, hydrodynamics, and chemistry in governing mass transfer in HF MSE. Overall, this work demonstrates the potential of HF MSE for efficient bio-isobutanol recovery and also provides practical guidelines on critical factors (packing fraction, partition coefficient, and shell-side velocity), aiding in the design and scaling of MSE systems for resource recovery.
Aqueous chemistry↗
A Eulerian three-fluid flow framework for the study of fuel dispersal behavior under loss-of-coolant accident conditions
As the nuclear industry explores new fuel designs to accommodate increased burnup, studying fuel behavior during loss-of-coolant accidents is essential for ensuring the safe operation of light-water reactors. Here, this study presents a computational model to simulate the intricate three-phase flow dynamics of fuel dispersal following cladding breach, where solid fuel particles and carrier fission gases interact with the surrounding liquid or gas medium. The model utilizes a Eulerian three-fluid framework to depict the phase interaction within interpenetrating continua, treating fragmented fuel as granular material. Closure relationships for bulk-flow properties of the solid phase and detailed interfacial exchange terms in high-temperature, high-speed gas-solid-liquid flow scenarios are proposed. The model was validated using adiabatic experimental results in the literature, focusing on characterizing fuel dispersal behavior during loss-of-coolant accidents conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched probability density functions reported in experimental studies. Key highlights of this study include the theoretical description of a coupled three-phase flow with phase change, interpretation of the high-pressure boundary conditions and insights into the transient behavior of gas-solid-liquid phase dynamics during a simulated high-speed dispersal event.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Structural Gating Enhances Long-Distance Light-Driven Interfacial Electron Transfer
Structural gating provides a molecular means to transfer electrons preferentially in one desired vectorial direction, a behavior needed for applications in artificial photosynthesis. At the interfaces utilized herein, visible-light absorption by a transition metal complex opens a “structural gate” by planarization of otherwise rotating phenyl rings in p-phenylene ethynylene (PE) bridge units. Planarization provides a conjugated pathway for electron flow toward a conductive oxide surface. Interfacial electron transfer to the oxide restores rotation and closes the gate to the unwanted recombination reaction. This structural gating results in nearly quantitative long-distance (>20 Å) interfacial electron transfer that occurs ~1000 times faster than transfer in the opposite direction. A comparative kinetic study of these complexes with those that contain ionic bridge units, without gating function, as a function of the applied potential and hence –ΔG° provided a physical basis for the structural gating. A small distance-dependent reorganization energy with weak electronic coupling underlies the success of this gate that enables efficient long-distance electron transfer and slow recombination.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Phase Separation, Capillarity, and Odd-Surface Flows in Chiral Active Matter
Active phase separations evade canonical thermodynamic descriptions and have thus challenged our understanding of coexistence and interfacial phenomena. Considerable progress has been made towards a nonequilibrium theoretical description of these traditionally thermodynamic concepts. Spatial parity symmetry is conspicuously assumed in much of this progress, despite the ubiquity of chirality in experimentally realized systems. Here, in this Letter, we derive a theory for the phase coexistence and interfacial fluctuations of a system that microscopically violates spatial parity. We find suppression of the phase separation as chirality is increased as well as the development of steady-state currents tangential to the interface dividing the phases. These odd flows are irrelevant to stationary interfacial properties, with stability, capillary fluctuations, and surface area minimization determined entirely by the capillary surface tension. Using large-scale Brownian dynamics simulations, we find excellent agreement with our theoretical scaling predictions.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Tracking Local pH Dynamics during Water Electrolysis via In-Line Continuous Flow Raman Spectroscopy
The performance of electrochemical devices, which play a critical role in decarbonization efforts, is often governed by proton-coupled electron transfer reactions at the electrode–electrolyte interface. These reactions are highly sensitive to the complex and dynamic microenvironment present at the electrode surface. However, characterizing this environment─particularly monitoring interfacial pH and its evolution under reaction conditions─remains challenging, necessitating the development of advanced analytical tools. Here, in this study, we introduce in-line continuous flow Raman spectroscopy (CFRS) as a spectroelectrochemical platform for quantifying interfacial pH swings generated during water-splitting. By monitoring phosphate ion speciation and controlling the hydrodynamics with a flow cell, we measure pH swings as a function of current density, flow rate, and distance from the electrode. Comparison with theoretical models reveals the impact of bulk pH, boundary layer thickness, and bubble dynamics at high current densities. Collectively, these findings establish CFRS as a platform for quantitatively investigating pH dynamics, offering critical insights for advancing electrochemical energy conversion technologies.
Marquez, Raul A. [Univ. of Texas, Austin, TX (Unit↗
Wetting mechanism and alteration of nano-sized shale pores: Insights from contrast variation small angle neutron scattering
Wettability of tight shale is crucial for fluid flow and mass transport process in energy geosciences. However, understanding the interfacial chemistry and wetting mechanisms at sub-nano-pore scales remains a formidable challenge. Here, in this study, the Contrast Variation technique of Small Angle Neutron Scattering (CV-SANS) is employed to investigate shale’s interfacial chemistry using reagents that possess a range of different polarities, including water, n-decane, toluene, and dimethyl methanamide. Through five different experimental strategies, we have demonstrated a successful modification of shale wettability, ranging from enhancement, weakening, to reversal. Delving into the mechanisms, we illustrated the crucial role of pre-existing liquid films in these changes, where the uniquely co-existing polar and non-polar functional groups in dimethyl methanamide acted as a conduit for interfacial chemistry adjustments. Furthermore, a solvent immersion led to matrix dilation as well as liberation of residual oil-occupied pores, resulting in altered pore size distributions, with hydrogen bonding playing a significant role in the polar groups. Interestingly, despite shale exhibiting a stronger affinity for oil over water, hydrophilic solvents induced more substantial dilation than lipophilic ones. Collectively, this work elucidates the dynamic change of interfacial chemistry via the configuration of polarity using chemical reagents, and the CV-SANS technique underscores its invaluable utilities in decoding the interfacial wettability traits in nanopore space of shale.
58 GEOSCIENCES↗
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↗
Tracking Thermal Transport in Colloidal Quantum Dot Films Using in Situ Time-Resolved X-ray Diffraction
Colloidal quantum dots (QDs) and their thin films are increasingly used in electronic and photonic devices replacing traditional bulk semiconductors. However, thermal properties of QDs remain underexplored relative to device development efforts. This study shows the use of time-resolved X-ray diffraction as a contact-free method to probe the thermal response of QDs in environments representative of the active layer in QD optoelectronic devices, providing in situ insights for future thermal management strategies. Through the extraction of Debye–Waller factors on a subnanosecond time scale, we directly capture the heating and cooling of core/shell CdSe/CdS QDs following pulsed optical excitation. In a QD thin film that actively provides optical gain, the thermal conductivity is found to be as low as 0.13 W m –1 K –1 , because of the poor heat flow between close-packed QD solids. For QDs dispersed in liquids, interfacial thermal conductance dominates thermal relaxation, with a conductance of 15 MW m –2 K –1 .
36 MATERIALS SCIENCE↗
Flow reversal benchmark of a one-sided heated narrow rectangular channel with CATHARE and RELAP5
Flow reversal in narrow coolant channels can be a crucial phenomenon for the safety of research reactors with a downward nominal flow direction. During a loss of forced flow accident, the downward flow stagnates briefly before transitioning into an upward natural circulation flow. The fuel may be damaged if dryout occurs and threshold fuel and/or cladding temperatures are exceeded. A comprehensive study is provided for flow reversal in narrow rectangular channels by examining experimental data and conducting software model analyses. The literature on flow reversal was reviewed, and selected experimental datasets were used to benchmark against CATHARE and RELAP5 models and also compare the code calculations with each other. The experimental data comes from flow reversal tests conducted with a narrow rectangular channel with one-sided heating. The results were compared with experimental data for successful flow reversal tests and predicted dryout power for dryout conditions. Also, the study examined the effects of the pump coastdown period, inlet liquid temperature, system pressure, and localized pressure drops. The experimental results showed that shorter coastdown periods, reduced pressure drops, and lower coolant inlet temperatures increased the dryout power. However, the system pressure did not noticeably affect the results. The simulation results showed that both CATHARE and RELAP5 agreed with experimental data, capturing the trends of the experimental results. Slight differences between each code calculation, as well as the predicted and measured dryout powers, were attributed to experimental uncertainties and the modeling of physical phenomena such as wall nucleation, interfacial heat transfer, drag coefficients, and critical heat flux. Overall, this study provides an understanding of flow reversal and the prediction capabilities of thermal-hydraulics software models. In conclusion, a future study of the flow reversal benchmark of a narrow rectangular channel with two-sided heating may provide additional valuable insights.
CATHARE↗
Interface‐Controlled Redox Chemistry in Aqueous Mn 2 ⁺/MnO₂ Batteries
Manganese dioxide (MnO 2 ) deposition/dissolution (Mn 2+ /MnO 2 ) chemistry, involving a two-electron-transfer process, holds promise for safe and eco-friendly large-scale energy storage. However, challenges like electrode/electrolyte interface environment fluctuations (H + and H 2 O activity), irreversible Mn degradation, and limited understanding of degradation mechanisms hinder the reversibility of the Mn 2+ /MnO 2 conversion. This study demonstrates a vanadyl/pervanadyl (VO 2+ /VO 2 + ) redox-mediated interface designed for high-energy Mn 2+ /MnO 2 batteries. Unlike flow systems, this work uncovers, for the first time, the mechanism of a static redox-mediated interface in regulating interfacial H + and H 2 O activities. Significantly, the VO 2+ /VO 2 + chemical redox mediation targets Mn 3+ intermediates, suppressing their hydrolysis and enabling 100% Mn 2+ /MnO 2 conversion. The redox-mediated interface enhances the Mn redox electron transfer process, achieving a stable ≈95% coulombic efficiency and ultrahigh capacity of 100 mAh cm −2 with an areal energy density of 111 mWh cm −2 , outperforming flow systems. The electrode also exhibits an average specific capacity of 593 mAh g −1 , approaching the theoretical limit of 616 mAh g −1 , and a specific energy density of 721 Wh kg −1 at high MnO 2 loadings (50–150 mg cm −2 ). Furthermore, the findings highlight the critical role of interfacial redox mediation in regulating H + and H 2 O activities and underscore the significance of interface dynamics.
Xue, Xinzhe [University of California, Santa Cruz,↗
Dilute Regeneration-Driven Membrane Capacitive Deionization of Synthetic Seawater Using Nanopatterned Membranes and Prussian Blue Analog Electrodes
Membrane capacitive deionization (MCDI) offers energy-efficient seawater desalination but is limited at high salinity by membrane resistance and incomplete electrode regeneration. Nanopatterned ion-exchange membranes, dilute regeneration protocols, and Prussian blue analog (PBA)-functionalized electrodes are combined in a flow-by-MCDI cell. Nanopatterned ion-exchange membranes (hexagonal, octagonal, double-ring, rectangular) enhance interfacial ion transport, with hexagonal geometry delivering ≈12.5% greater surface area and the best performance. PBA-functionalized electrodes increase salt adsorption and charge-transfer kinetic rates. The integrated system lowers the area-specific resistance by 45 Ω cm2, resulting in a 500 mV reduction in the cell voltage for a current density of 2 mA cm−2 for a 35 000 ppm NaCl feed. This improves the energy-normalized salt adsorption six fold (64–382 mmol J−1). Low salinity (2000 ppm) and mixed-salt regeneration sustains a ≈39% water recovery and stable performance for at least seven cycles. Overall, combining nanopatterned membranes, which promote confinement-enhanced ion mobility, and PBA electrodes, which enhance salt adsorption, improved the energy efficiency of MCDI.
Hasan, Mahmudul↗