Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interfacial voids”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Development and Validation of Two-Phase Flow Models in MOOSE and Application to Molten Salt Reactors

Two-phase flow in Molten Salt Reactors (MSRs) is important as it impacts reactivity evolution, reactor transient response, and the removal of species dissolved in the molten salt through gas phase transfer. Therefore, accurately predicting the gas distribution and the associated liquid-gas interface area in MSRs is essential for their design and operation. Recently, we integrated two new models into Idaho National Laboratory (INL)’s Multiphysics Object-Oriented Simulation Environment (MOOSE): a multi-D generalization of a mixture drift-flux model and a Euler-Euler model. The Euler-Euler model offers higher fidelity, while the mixture drift-flux model provides greater computational efficiency, which is typically preferred for modeling reactor transients. However, the mixture model's accuracy in capturing void distribution and interfacial area in MSRs still needs to be assessed. This article begins with a description of the mathematical framework for the two-phase models implemented in MOOSE. It then presents validation of these models against relevant experimental data. Finally, both models are applied to the Molten Salt Reactor Experiment case study, analyzing various operational conditions such as different rates of fission product volatilization and diverse cover gas entrainment scenarios at the reactor pump. The article concludes by assessing the suitability of both models for capturing the two-phase flow dynamics critical to MSR operations.

42 - ENGINEERING↗

Triiodide Attacks the Organic Cation in Hybrid Lead Halide Perovskites: Mechanism and Suppression

Molecular I 2 can be produced from iodide-based lead perovskites under thermal stress; triiodide, I 3 − , is formed from this I 2 and I − . Triiodide attacks protic cation MA + - or FA + -based lead halide perovskites (MA + , methylammonium; FA + , formamidinium) as explicated through solution-based nuclear magnetic resonance (NMR) studies: triiodide has strong hydrogen-bonding affinity for MA + or FA + , which leads to their deprotonation and perovskite decomposition. Triiodide is a catalyst for this decomposition that can be obviated through perovskite surface treatment with thiol reducing agents. In contrast to methods using thiol incorporation into perovskite precursor solutions, no penetration of the thiol into the bulk perovskite is observed, yet its surface application stabilizes the perovskite against triiodide-mediated thermal stress. Thiol applied to the interface between FAPbI 3 and Spiro-OMeTAD (“Spiro”) prevents oxidized iodine species penetration into Spiro and thus preserves its hole-transport efficacy. Surface-applied thiol affects the perovskite work function; it ameliorates hole injection into the Spiro overlayer, thus improving device performance. It helps to increase interfacial adhesion (“wetting”): fewer voids are observed at the Spiro/perovskite interface if thiols are applied. Perovskite solar cells (PSCs) incorporating interfacial thiol treatment maintain over 80% of their initial power conversion efficiency (PCE) after 300 h of 85 °C thermal stress.

36 MATERIALS SCIENCE↗

A Morphologically Stable Li/Electrolyte Interface for All-Solid-State Batteries Enabled by 3D-Micropatterned Garnet

Morphological degradation at the Li/solid-state electrolyte (SSE) interface is a prevalent issue causing performance fading of all-solid-state batteries (ASSBs). To maintain the interfacial integrity, most ASSBs are operated under low current density with considerable stack pressure, which significantly limits their widespread usage. In this work, a novel 3D-micropatterned SSE (3D-SSE) that can stabilize the morphology of the Li/SSE interface even under relatively high current density and limited stack pressure is reported. Under the pressure of 1.0 MPa, the Li symmetric cell using a garnet-type 3D-SSE fabricated by laser machining shows a high critical current density of 0.7 mA cm –2 and stable cycling over 500 h under 0.5 mA cm –2 . This excellent performance is attributed to the reduced local current density and amplified mechanical stress at the Li/3D-SSE interface. These two effects can benefit the flux balance between Li stripping and creep at the interface, thereby pre-venting interfacial degradation such as void formation and dendrite growth.

3D-micropatterned garnet↗

Benchmark of interfacial area concentration approaches for the two-fluid model in gas-dispersed condensing flow

Accurate predictions of interfacial area concentration (IAC) are required for closure of the condensation phase change term. The interfacial area transport equation (IATE) can dynamically predict interfacial area concentration and is increasingly being used in the benchmarking of phase change flows over static based correlations which require a flow regime map. For higher flow regimes, the two-group two-fluid IATE model can separately model the transport of spherical bubbles and cap/slug bubbles. Recently, a modified formulation of the one-group and two-group two-fluid model with interfacial area transport equation was proposed for gas-dispersed condensing flows. Three datasets have been used for validation consisting of condensing two-phase flow data in the unheated section of a vertical annulus with conditions spanning a range of pressures, inlet subcoolings and mass fluxes for benchmarking the various IAC approaches. For bubbly flows, the coupled void transport - IATE simulation results in good agreement in the predictions of the void fraction and the interfacial area concentration under nearly all sixty conditions. Considering the proposed two-group two-fluid IATE model formulation, the group-2 Nusselt number correlation is investigated to account for the appropriate heat transfer length scale for group-2 bubbles. Validation of the proposed two-group two-fluid model against the new dataset shows good agreement with experimental data with significant improvement over the one-group IATE model for conditions with substantial group-2 void fraction throughout the simulation domain.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗

Experimental and statistical study on the effect of process parameters on the quality of continuous fiber composites made via additive manufacturing

Ongoing research in additive manufacturing towards structural and industrial application has led to the use of commingled roving as a manufacturing feedstock for printing high fiber volume fraction composites. The prospects of using this technology for high performance applications necessitates the need for a comprehensive experimental investigation into the effects of processing parameters on the quality of an additively manufactured composite printed from commingled roving feedstock. Here, in this work, transverse flexure and void fraction matrix pyrolysis testing are both performed to evaluate composite quality. The transverse flexure test is a testing approach that evaluates the quality of the interfacial fiber-matrix bond while the void fraction test estimates the void content in the printed composite. A full observational study consisting of 27 different test combinations is done to investigate the effects of three different process parameters namely, temperature, pressure, and print speed across three different levels. Composite samples were made from commingled roving of E-glass and amorphous PET using an in-house built continuous fiber composite digital manufacturing system. Least squares regression analysis is performed to study the main, interaction and quadratic effects of process parameters. A statistical regression model having an R2 adjusted value of 80.1% is generated from the transverse flexure study, which is used to explain main and interaction effects and also predict performance. Response surface plots are also generated and are used to optimize process parameters which can subsequently be of help in scaling up composite manufacturing. Results show that all three process parameters are highly statistically significant at the 0.01 level of significance. Pressure * Temperature and Pressure * Printspeed are significant interaction terms. Pressure plays a weightier role when print speed is increased or temperature is decreased as it closes more voids that would ordinarily have been introduced because of drop in polymer melt viscosity. Micrographic analysis is also performed.

36 MATERIALS SCIENCE↗

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗

Improvements to CTF Closure Models for Modelling of Two-Phase Flow

This report documents the efforts to improve the CTF prediction of void fraction and two-phase pressure drop by improving two-phase closure models. Previous validation activities have revealed that CTF tends to overpredict void fraction and two-phase pressure drop. In response to this, two approaches were taken to improve CTF’s predictive capabilities. Based on findings that the interfacial drag and subcooled boiling models significantly impact void prediction, alternative closure models for these physical effects were found and implemented into the code. An extensive assessment was performed by using the existing and newly added two-phase experimental data, which showed that void prediction, wall temperature, and two-phase pressure-drop results are improved by using the newly implemented models. A second approach for improving CTF modeling accuracy involved the use of a Bayesian calibration process that uses CTF validation data to optimize selected modeling coefficients and closure model multipliers to achieve a more accurate prediction of experimental results. A similar assessment was performed with calibrated models that improved the void and pressure-drop prediction for test cases that were both included and not included in the calibration dataset. Results of this study will be used to change the models used in CTF to achieve more accurate BWR analyses moving forward. This study also identified opportunities for improving additional closure models and other opportunities to use calibration techniques to improve CTF.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Probing the Impact of Vacancy Diffusion on Void Dynamics at the Lithium Metal–Solid Electrolyte Interface

Lithium (Li) metal-based solid-state batteries (SSBs) are considered promising candidates for next-generation energy storage due to their superior energy density and enhanced safety compared to conventional Li-ion systems. However, their practical application is limited by challenges such as void formation at the Li-solid electrolyte (SE) interface, which disrupts ion transport and accelerates interfacial degradation. This work investigates how the coupled effects of electro-dissolution kinetics and surface diffusion at the Li metal surface govern the evolution of interfacial morphology during stripping. This work examines the influence of three distinct surface diffusion modes, which are terrace diffusion, step diffusion, and interlayer diffusion, on maintaining interfacial stability. In addition, how the dominant surface diffusion mechanism can overcome the contact loss due to high reaction kinetics is explored. Furthermore, the roughness of the Li metal anode surface is quantified, and the influence of different diffusion mechanisms on the evolution of the dynamic solid–solid interface is examined. The critical role of temperature in enhancing Li surface diffusivity and expanding the regime of stable contact is highlighted. By identifying distinct regimes of interface stability, this study analyzes how non-uniform electrochemical dynamics dictate void morphology evolution and interfacial contact. These insights offer guiding principles for engineering robust Li–SE interfaces in SSBs.

lithium metal anode↗

Adverse Effects of Trace Non-polar Binder on Ion Transport in Free-standing Sulfide Solid Electrolyte Separators

Sulfide solid-state electrolyte (SE) possesses high room-temperature ionic conductivity. However, fabrication of the free-standing, sheet-type thin sulfide SE film electrolyte to enable all-solid-state batteries to deliver high energy and power density remains challenging. Herein we show that argyrodite sulfide (Li 6 PS 5 Cl) SE can be slurry cast to form free-standing films with low (≤5 wt%) loadings of poly(isobutylene) (PIB) binder. Two factors contribute to a lower areal specific resistance (ASR) of the thin film SEs benchmarked to the pristine powder pellet SSE counterparts: i) 1–2 orders reduced thickness and ii) reasonably comparable ionic conductivity at room temperature after the isostatic pressing process. Nevertheless, an increasing polymer binder loading inevitably introduced voids in the thin film SEs, compromising anode/electrolyte interfacial ion transport. Our findings highlight that electrolyte/electrode interfacial stability, as well as the selection of slurry components, including sulfide SE, binder, and solvent, play essential roles in thin film sulfide electrolyte development.

25 ENERGY STORAGE↗

Hierarchical assembly of ZnO nanowire trunks decorated with ZnO nanosheets for lithium ion battery anodes

Hierarchical architectures composed of nanomaterials in different forms are essential to improve the performance of lithium-ion battery (LIB) anodes. Here, we systematically studied the effects of hierarchical ZnO nanostructures on the electrochemical performance of LIBs. ZnO nanowire (NW) trunks were decorated with ZnO NWs or ZnO nanosheets (NSs) by successive hydrothermal synthesis to create hierarchical three-dimensional nanostructures. The branched ZnO NSs on the ZnO NW trunk exhibited a two-fold higher specific gravimetric capacity compared to ZnO NWs and branched ZnO NWs on ZnO NW trunks after 100 cycles of charging–discharging at 0.2C (197.4 mA g –1 ). The improvement in battery anode performance is attributable to the increased interfacial area between the electrodes and electrolyte, and the void space of the branched NSs that facilitates lithium ion transport and volume changes during cycling.

25 ENERGY STORAGE↗

Simulated signatures of ignition

Ignition on the National Ignition Facility (NIF) provides a novel opportunity to evaluate past data to identify signatures of capsule failure mechanisms. We have used new simulations of high-yield implosions as well as some from past studies in order to identify unique signatures of different ignition failure mechanisms: jetting due to the presence of voids or defects, jetting due to the capsule fill tube, interfacial mixing due to instabilities or due to plasma transport, radiative cooling due to the presence of contaminant in the hot spot, long-wavelength drive asymmetry, and preheat. Many of these failure mechanisms exhibit unique trajectories that can be distinguished through variations in experimental observables such as neutron yield, down-scattered ratio (DSR), and burn width. Our simulations include capsules using both plastic and high-density carbon ablators and span all high-yield designs considered since the beginning of the National Ignition Campaign in 2011. We observe that the variability in trajectories through the space of neutron yield, DSR, and burn width varies little across capsule design yet are unique to the failure mechanism. The experimental trajectories are most consistent with simulated preheat and jetting due to voids and defects, which are the only failure mechanisms that are indistinguishable in our analysis. This suggests that improvements to capsule compression due to improved capsule quality or reduced preheat have played a primary role in enabling high yields on NIF. Furthermore, our analysis suggests that further improvements have the potential to increase yields further.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Liquid-phase turbulence measurements in air-water two-phase flows using particle image velocimetry

Liquid-phase turbulence plays a vital role in determining various gas-liquid two-phase flow parameters, such as void fraction distribution, bubble morphology, bubble-bubble interactions, and interfacial area concentration. In this study, a two-phase flow database including both the gas- and liquid-phases measurements was developed focusing on three bubbly flow conditions in an air-water two-phase flow loop with a vertical one-inch diameter circular pipe test section. A particle image velocimetry (PIV) system integrating an optical phase separation method, i.e., the planar laser-induced fluorescence (PLIF) technique using fluorescent particles and optical filtration, was applied to measure the liquid-phase turbulence information, including the time-averaged velocity, Reynolds stress, and turbulent kinetic energy for the liquid phase. The PIV measurements were taken at three ports along the test section at 14.5, 51.5, and 88.5 pipe inner diameters downstream of a bubble injector. In addition, a double-sensor conductivity probe was used to measure radial distributions of the local time-averaged void fraction and gas velocity. The measured liquid-phase turbulence was used to benchmark Sato's turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. Lastly, in the two bubbly flows tested that have low void fractions being less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction.

42 ENGINEERING↗

CTF Improved Drag Model and Flow Regime Transition Criteria

The demand for accurate prediction of two-phase flow behavior in a boiling water reactor (BWR) requires a comprehensive understanding of flow regime, void fraction, heat transfer, and pressure drop. The CTF subchannel code, which is used for the Thermal/Hydraulic (T/H) solution in the Consortium for Advanced Simulation of Light Water Reactors (CASL)-developed Virtual Environment for Reactor Application (VERA) core simulator, is being further developed for BWR applications. In support of this goal, the present work highlights some of the two-phase closure model developments towards improving the CTF void fraction prediction, especially for subcooled boiling. The drift-flux approach has been well-developed for upward dispersed two-phase flows and proven to be accurate in predicting void fraction in bubbly and slug flow regimes. In this work, these kinematic constitutive relations for the drift-flux velocity have been implemented into CTF to describe the interfacial drag of bubbly flow as an alternative to the existing model for better void fraction prediction. The success of these constitutive relations also relies on a good flow regime map that accounts for flow conditions and channel geometry. A more reliable flow regime transition criteria that account for the flow condition has also been implemented in this study for modeling the flow regime transition criteria. The newly implemented models are shown to give improved void fraction predictions in comparison to experimental data.

Hizoum, Belgacem↗

Hybrid additive manufacturing of AISI 316L via asynchronous powder and hot-wire laser directed energy deposition

Hybrid Additive Manufacturing (AM) offers a way to leverage the advantages of different AM technologies, enabling the efficient production of sizeable parts without compromising material properties or geometric complexity capabilities. This study presents an asynchronous hybrid Directed Energy Deposition (DED) strategy employing laser powder DED and laser hot-wire DED. AISI 316L parts comprising multiple powder and wire segments were fabricated with optional machining on AISI 316L substrates to investigate how quality is impacted by (i) alternative process sequences (laser powder DED followed by laser hot-wire DED and vice versa), (ii) machined vs. as-printed interfacial conditions, and (iii) material deposition on top vs. alongside previously built segments. Optical microscopy, X-ray computed tomography, and Vickers hardness were used to characterize the morphology and microstructure of the parts, localized porosity and lack of fusion defects, bulk density, and mechanical properties. Interfacial machining was necessary for dimensional control but promoted lack of fusion voids, resulting in a 99.71 ± 0.01% dense part. As-printed interfaces resulted in a denser part (99.82 ± 0.02%) at the expense of dimensional accuracy. The hardness of the parts with as-printed and machined interfaces was 196 ± 0.37 HV and 192 ± 0.40 HV, respectively, compared to 156 ± 1.4 HV for the substrate. Depositing powder alongside or on top of wire sections resulted in interfaces with a hardness of 217 ± 2.2 HV, compared to 185 ± 3.4 HV for the wire-powder interfaces.

36 MATERIALS SCIENCE↗

A mechanism for reduced compression in indirectly driven layered capsule implosions

High-yield implosions on the National Ignition Facility rely on maintaining low entropy in the deuterium–tritium fuel, quantified by its adiabat, in order to efficiently couple energy to the hot spot through high compression of the fuel layer. We present very-high-resolution xRAGE simulation results that study the impacts of interfacial mixing and the jetting of materials due to surface defects, defects on internal interfaces, voids, and engineering features on fuel layer compression. Defects and voids are typically neglected in implosion simulations due to their small size and three-dimensional geometry. Our results showed that supersonic jets of material arise through weak spots in the shell at peak implosion velocity that prevent uniform compression of the fuel layer even when they do not introduce contaminant into the hot spot. This occurs despite maintaining low fuel entropy, since the formation of the weak spots involves nonradial displacement of fuel mass. In contrast, simulations show that fuel–ablator mixing due to interfacial instabilities has a much smaller impact on compression. We show that defects on interior interfaces of plastic capsules decrease compression by 15% to 25% and interfacial mixing between the ablator and fuel decreases compression by less than 1% for implosions with plastic or high-density carbon (HDC) ablators. For low adiabat implosions, the impact of jetting seeded by the support tent can also decrease the compression by 25%. We demonstrate that the inclusion of interior defects in simulations can explain the inferred compression in two fielded plastic capsule implosions and that the inclusion of voids, for which available characterization has large uncertainties, in simulations of HDC capsule implosions has a qualitatively consistent impact. This mechanism offers a potential explanation for persistently overestimated fuel compression in design simulations of layered implosions on the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Designing carbon wetting layers with inorganic salts as additives for enhanced adhesion and molten sodium wettability in electrochemical application

ß” alumina solid-state electrolyte (BASE) is considered one of the most promising materials for sodium batteries due to its high ionic conductivity and exceptional stability with molten sodium. However, at lower operating temperatures, the poor physical contact at the liquid-solid interface between molten sodium and BASE results in high interfacial resistance. In this study, we explore the optimization of carbon wetting layers with various salt additives to enhance adhesion and improve molten sodium wettability. Using a spray-coating method, carbon layers incorporating sodium hexametaphosphate (Na6(PO3)6, SHMP) as an additive demonstrate superior adhesion, mechanical stability, and enhanced wettability with molten sodium. Scanning Electron Microscopy (SEM) analysis reveals that SHMP preserves the porous carbon network required for efficient sodium transport while minimizing structural defects such as cracks and voids. Electrochemical testing using Na symmetric cells confirms that SHMP-modified layers significantly reduce interfacial resistance and overpotential, outperforming No salts counterparts. These findings highlight that the strategic incorporation of suitable additives, such as SHMP, in the design of microscale carbon layers could significantly enhance the performance of molten sodium batteries at lower operating temperatures.

Han, Henry H.↗