Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flow fields”

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

Effect of parallel flow on resonant layer responses in high beta plasmas

Abstract Resonant layers in a tokamak respond to non-axisymmetric magnetic perturbations by amplifying the mode amplitude and balancing the plasma rotation through magnetic reconnection and force balance, respectively. This resonant response can be characterized by local layer parameters and especially by a single quantity in the linear regime, the so-called inner-layer Δ. The computation of Δ under two-fluid drift-MHD formalism has been progressed by reducing the order of the system in the phase space, where the shielding current is approximated as being only carried by electrons, a posteriori . In this study, we relax the approximation and compute Δ accounted for by the parallel flow associated with the ion shielding current. The posteriori is numerically verified in great agreement with the original SLAYER developed in a previous paper (J.-K. Park 2022 Phys. Plasmas 29 072506). Extending the resonant layer response theory to high β plasmas, our research findings answer two important questions: how the parallel flow influences the resonant layer response and why the parallel flow effect appears in high β plasmas. The complicated plasma compression in high β regime allows the parallel flow response to give rise to the ion shielding current, which not only shifts the zero-crossing condition of the ExB flow but also enhances the field penetration threshold. Technically, the Riccati matrix transformation method is adapted to handle the numerical stiffness due to the increased order of the system. The high fidelity of this numerical method makes use of further extension of the model to higher-order systems to take other physical phenomena into account. This work is envisaged to predict the resonant layer response under high β fusion reactor conditions.

Lee, Yeongsun (ORCID:000000034474416X)↗
DOE OSTI2025

REBOUND: Reverse Engineering Bidirectional Outflow Under Non-Equilibrium Diffusion

Rare-earth elements (REEs) are essential for electronics, renewable energy, and defense technologies. However, the current supply of REEs relies on mining concentrated in a few countries and energy-intensive separations. DOE’s Basic Energy Sciences (BES) program has launched a grand challenge which aims to ensure a sustainable supply of critical REEs by developing innovative and environmentally friendly separation methods. As an alternative to costly and harmful traditional methods, the Non-Equilibrium Transport Driven Separations (NETS) initiative has created a microfluidic Y-channel co-flow method that applies external fields to exploit magneto- and electrohydrodynamic effects for separating dilute REE ions from complex feedstocks. Computational fluid dynamics (CFD) studies have identified a few operating conditions with promising ion selectivity and separation efficiency. However, challenges remain regarding Y-channel versatility across feedstocks and accurate incorporation of physical phenomena into CFD models. In this work, we develop a multi-fidelity modelling approach which integrates experimental results with CFD simulation to build a surrogate model for the dependence of separation efficiency to variation of design parameters. The surrogate model enables a reinforcement learning (RL) method to adaptively launch CFD and experimental runs, improving model fidelity around optimal Y-channel parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Potential vorticity conservation for plasma turbulence in an inhomogeneous magnetic field: Theory and implications

The concept and theory of potential vorticity in drift wave turbulence are extended to the case of an inhomogeneous magnetic field. A one-field magnetic potential vorticity conserving equation is derived via the use of conservative gyrokinetics. The similarity between the corresponding systems for drift wave turbulence and shallow water theory is discussed in detail. Zonal flow physics in an inhomogeneous magnetic field is discussed. In particular, a Charney–Drazin type nonacceleration theorem is derived from the novel system, which conserves magnetic potential vorticity. Extensions of the turbulent equipartition theory to the transport of magnetic potential vorticity are proposed.

Physics↗
DOE OSTI2025

Antielectrophoretic Response-Driven Bending–Tilting Deformation of Cationic Polyelectrolyte Brushes Drives Nonlinear Electroosmotic Transport in Brush-Grafted Nanochannels

In this paper, we use all-atom molecular dynamics (MD) simulations to describe a non-linearly enhanced electroosmotic (EOS) flow, where, in a nanochannel grafted with cationic PMETAC ([Poly(2-(Methacryloyloxy)Ethyl) Trimethylammonium Chloride]) brushes, a two-fold increase in the electric field strength leads to a several-fold (more than two-fold) increase in the EOS flow strength and volume flow rate. The electric field enforces the PMETAC brushes to undergo a bending-tilting driven deformation with a significant portion of the brush layer becoming parallel to the grafting surface. In response, a substantial fraction of the counterions leave the brush layer (hence become more mobile), but instead of going into the bulk, accumulate at the brush-bulk interface, i.e., stay in proximity of the brush segments aligned parallel to the grafting surface. This creates an interesting situation, where the counterions are not completely within the brush layer, yet they fully screen the brush charges. Such “freer” conditions enable the counterions to achieve very high velocity, thereby ensuring that the water solvating the counterions themselves move very fast triggering the significantly augmented EOS transport. Probing deeper we can identify that the bending-tilting driven brush deformation, enforcing the brushes to align parallel to the substrate, results from the anti-electrophoretic behavior of the brushes, where despite being positively charged, the brushes move against the electric field direction. Such an anti-electrophoretic behavior of the PE brushes, which has not been reported before, can be associated with the very fast velocities of the negatively charged counterions and the electrostatic and hydrodynamic coupling of the counterions with the positive functional groups of the brushes. Here, we anticipate that the findings of this paper will shed light on strategies for nanochannel flow, the anti-electrophoretic response of charged polymer chains, and the significance of capturing the detailed chemical architecture of polyelectrolytes in nanoscale science and engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On anomalous transport of multi-species plasma associated with the resistive ballooning and resistive drift waves driven turbulence

Anomalous transport of multi-species plasma related to the resistive ballooning and resistive drift wave turbulence is considered in a “cold” ion approximation. It is found that similar to the resistive drift wave turbulence [see A. R. Knyazev and S. I. Krasheninnikov, Phys. Plasmas 31, 012502 (2024); and S. I. Krasheninnikov and R. D. Smirnov, Phys. Plasmas (to be published)] the addition of the ballooning drive does not change the main features of anomalous transport of the multi-species plasma: (i) The transport of all ion species is described as a transport of the passive scalars in the turbulent field of the electrostatic potential and electron density perturbation; (ii) the density of ion species with a larger ratio of the mass to charge has the tendency to the accumulation/depletion in the vortices of plasma flow; and (iii) the cross-field transport of all plasma species (including electrons and ions) is described by the same anomalous transport coefficient.

Physics↗
DOE OSTI2025

Linear Discriminant Analysis-Based Machine Learning and All-Atom Molecular Dynamics Simulations for Probing Electro-Osmotic Transport in Cationic-Polyelectrolyte-Brush-Grafted Nanochannels

Deciphering the correct mechanisms governing certain phenomena in polyelectrolyte (PE) brush grafted systems, revealed through atomistic simulations, is an extremely challenging problem. In a recent study, our all-atom molecular dynamics (MD) simulations revealed a non-linearly large electroosmotic (EOS) flow (in the presence of an applied electric field) in nanochannels grafted with PMETAC [Poly(2-(methacryloyloxy)ethyl trimethylammonium chloride] brushes. Given the lack of any formal procedure that would have directed us to identify the correct factors responsible for such an occurrence, we needed to spend several months and devote significant analyses to unravel the involved mechanisms. In this paper, we propose a Linear Discriminant Analysis (LDA) based Machine Learning (ML) approach to address this gap. At first, we obtain data on certain basic features from the all-atom MD data. These basic features represent the number of atoms of certain species around one atom of another (or same) species. Here, we obtain such data on basic features for a reference case (case of an EOS flow in PMETAC-brush-grafted nanochannels with a smaller electric field) and a perturbed case (case of an EOS flow in PMETAC-brush-grafted nanochannels with a larger electric field) in bins in which the nanochannel half height has been divided into. These datasets are high-dimensional dataset, to which the LDA is applied. This leads to the projection of the data (between the reference and the perturbed states) in a highly separated form on a 1D line. From such LDA calculations, we are able to identify the relative importance of the different basic features in ensuring this separation of the data (between the reference and the perturbed states) on the 1D line. This relative importance of the different basic features is quantified as “importance scores” for the different features, which in turn tell us what to study and where to study. Such knowledge enables us to rapidly identify the key factors responsible for the non-linearly large EOS transport in PMETAC-brush-grafted nanochannels.

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

Direct, simple and efficient computation of all components of the virtual-casing magnetic field in axisymmetric geometries with Kapur–Rokhlin quadrature

In a recent publication (Toler et al., J. Plasma Phys., vol. 89, issue 2, 2023, p. 905890210), we demonstrated that for axisymmetric geometries, the Kapur–Rokhlin quadrature rule provided an efficient and high-order accurate method for computing the normal component, on the plasma surface, of the magnetic field due to the toroidal current flowing in the plasma, via the virtual-casing principle. The calculation was indirect, as it required the prior computation of the magnetic vector potential from the virtual-casing principle, followed by the computation of its tangential derivative by Fourier differentiation, to obtain the normal component of the magnetic field. Our approach did not provide the other components of the virtual-casing magnetic field. In this letter, we show that a more direct and more general approach is available for the computation of the virtual-casing magnetic field. The Kapur–Rokhlin quadrature rule accurately calculates the principal value integrals in the expression for all the components of the magnetic field on the plasma boundary, and the numerical error converges at a rate nearly as high as the indirect method we presented previously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Hot-spot model for inertial confinement fusion implosions with an applied magnetic field

Imposing a magnetic field on inertial confinement fusion implosions magnetizes the electrons in the compressed fuel; this suppresses thermal losses, which increases temperature and fusion yield. Indirect-drive experiments at the National Ignition Facility with 12 and 26 T applied magnetic fields demonstrate up to 40% increase in temperature, 3× increase in fusion yield, and indicate that magnetization alters the radial temperature profile [Moody et al., Phys. Rev. Lett. 129, 195002 (2022); Lahmann et al., APS DPP (2022)]. In this work, we develop a semi-analytic hot-spot model, which accounts for the two-dimensional (2D) Braginskii anisotropic heat flow due to an applied axial magnetic field. First, we show that hot-spot magnetization alters the radial temperature profile, increasing the central peakedness, which is most pronounced for moderately magnetized implosions (with 8–14 T applied field), compared to both unmagnetized (with no applied field) and highly magnetized (with 26 T or higher applied field) implosions. This model explains the trend in the experimental data, which finds a similarly altered temperature profile in the 12 T experiment. Next, we derive the hot-spot model for gas-filled (Symcap) implosions, accounting for the effects of magnetization on the thermal conduction and in changing the radial temperature (and density) profiles. Using this model, we compute predicted central temperature amplification and yield enhancement scaling with the applied magnetic field. The central temperature fits the experimental data accurately, and the discrepancy in the yield suggests a systematic (independent of applied field) degradation, such as mix, and additional degradation in the reference unmagnetized shot, such as reduced laser drive, increased implosion asymmetry, or the magnetic field suppressing ablator mixing into the hot-spot.

Alpha particles↗
DOE OSTI2025

Field Testing of an Affordable Zero-Liquid-Discharge Arsenic-Removal Technology for a Small-Community Drinking Water System in Rural California

Arsenic contamination in groundwater threatens public health, particularly in small, low-income communities lacking affordable treatment solutions. This study investigated the field implementation of novel air cathode assisted iron electrocoagulation (ACAIE) technology for arsenic removal in Allensworth, California, where groundwater arsenic concentrations exceeded 250 µg/L. Over four months, a pilot-scale ACAIE system, operating at 600 L/h, consistently reduced arsenic levels to below the EPA’s maximum contaminant level of 10 µg/L. Laboratory experiments informed the optimization of charge dosage and flow rates, which were validated during field testing of the ACAIE 600 L/h system. The in-situ generation of hydrogen peroxide at the cathode speeded up the reaction kinetics, ensuring high arsenic removal efficiency while allowing high throughput, even with a compact reactor size. An economic analysis demonstrated a treatment cost of USD 0.02/L excluding labor, highlighting the system’s affordability compared to conventional methods. Adding labor costs increased the treatment cost to USD 0.09/L. The regeneration of air cathodes extended their operational life, addressing a key maintenance challenge, thus reducing the costs slightly. Intermittent challenges were encountered with filtration and secondary contaminant removal; these issues highlight opportunities for further operational improvements. Despite these challenges, ACAIE’s low operational complexity, scalability, and cost-effectiveness make it a promising solution for underserved small communities. These findings provide critical insights into deploying sustainable arsenic remediation technologies that are tailored to the needs of rural, low-resource communities.

Bandaru, Siva R S↗
DOE OSTI2026

Implementation of Pfirsch–Schlüter parallel flow in x-ray imaging crystal spectrometer inversion analysis

The x-ray imaging crystal spectrometer (XICS) tomographic inversion code for Wendelstein 7-X (W7-X) has been modified to consider the effects of parallel flows and has been applied to analyze measurements taken during recent experimental campaigns. Previous analysis neglected the effects of parallel flows due to the primarily perpendicular geometry of the sightlines and the small magnitude predicted by neoclassical theory. To reconsider these effects, the incompressibility condition for plasma flows is used to calculate the parallel Pfirsch–Schlüter flow component for the equilibrium configuration. By incorporating this condition along with the geometry of the sightlines—i.e. the fractional contributions of perpendicular and parallel flows—, an updated expression for the measured flow is used for the profile inversion. Application of this modified inversion code to data from W7-X shows that the magnitude of the radial electric field and the flux surface averaged perpendicular flow are reduced by approximately a factor of 2 and brought into better agreement with neoclassical predictions and charge exchange recombination spectroscopy measurements.

Pfirsch–Schlüter flows↗
DOE OSTI2026

Towards a Quintic Ginzburg-Landau Description of the (2,7) Minimal Model

We discuss dimensional continuation of the massless scalar field theory with the 𝑖⁢𝜙 5 interaction term. It preserves the so-called 𝒫⁢𝒯 symmetry, which acts by 𝜙 →−𝜙 accompanied by 𝑖 →−𝑖. Below its upper critical dimension 10/3, this theory has interacting infrared fixed points. We argue that the fixed point in 𝑑 = 2 describes the nonunitary minimal conformal model 𝑀⁡(2,7). We identify the operators 𝜙 and 𝜙 2 with the Virasoro primaries 𝜙 1,2 and 𝜙 1,3 , respectively, and 𝑖⁢𝜙 3 with a quasiprimary operator, which is a Virasoro descendant of 𝜙 1,3 . Our identifications appear to be consistent with the operator product expansions and with considerations based on integrability. Using constrained Padé extrapolations, we provide estimates of the critical exponents in 𝑑 = 3. We also comment on possible lattice descriptions of 𝑀⁡(2,7) and discuss RG flows to and from this conformal field theory (CFT). Finally, we conjecture that the minimal models 𝑀⁡(2,2⁢𝑛 +1) are described by the massless scalar field theories with the 𝑖⁢𝜙 2⁢𝑛−1 interaction terms.

Conformal field theory↗
DOE OSTI2024

Impact of the Out‐Of‐Plane Flow Shear on Magnetic Reconnection at the Flanks of Earth's Magnetopause

Abstract Magnetic reconnection changes the magnetic field topology and facilitates the energy and particle exchange at magnetospheric boundaries such as the Earth's magnetopause. The flow shear perpendicular to the reconnecting plane prevails at the flank magnetopause under southward interplanetary magnetic field conditions. However, the effect of the out‐of‐plane flow shear on asymmetric reconnection is an open question. In this study, we utilize kinetic simulations to investigate the impact of the out‐of‐plane flow shear on asymmetric reconnection. By systematically varying the flow shear strength, we analyze the flow shear effects on the reconnection rate, the diffusion region structure, and the energy conversion rate. We find that the reconnection rate increases with the upstream out‐of‐plane flow shear, and for the same upstream conditions, it is higher at the dusk side than at the dawn side. The diffusion region is squeezed in the outflow direction due to magnetic pressure which is proportional to the square of the Alfvén Mach number of the shear flow. The out‐of‐plane flow shear increases the energy conversion rate , and for the same upstream conditions, the magnitude of is larger at the dusk side than at the dawn side. This study reveals that out‐of‐plane flow shear not only enhances the reconnection rate but also significantly boosts energy conversion, with more pronounced effects on the dusk‐side flank than on the dawn‐side flank. These insights pave the way for better understanding the solar wind‐magnetosphere interactions.

Liang, Haoming↗
DOE OSTI2025

Numerical investigation of magneto-inertial fusion targets magnetized by dynamic enforcement of helical current flow

Magnetized Liner Inertial Fusion (MagLIF) targets require premagnetization to reduce thermal conduction losses from the laser-heated fuel to the liner material during implosion. MagLIF targets are typically magnetized by external magnetic field coils which are technologically limited to providing <20 T axial magnetic fields in the fusion fuel. We present a novel target design, AutoMag-D, which employs dynamic enforcement of helical current in the liner, resulting in >20 T axial magnetic field in the fuel region prior to implosion without the need of external magnetic field coils. These self-magnetizing liners are made of helically oriented alternating electrically conductive material and electrically insulating material surrounded by a thin, conductive outer radial layer. As the liner is pulsed with a ∼ 20 MA, ∼ 100 ns rise time drive current from Z, the outer layer of conductive material is shocked and intensely Joule heated, causing it to melt, vaporize, and turn to plasma. This allows the magnetic drive field to diffuse radially inward which dynamically enforces current flow in the helical conduction paths in the liner and produces axial magnetic field in the fuel region prior to implosion of the inner liner surface. AutoMag-D liner designs do not require dielectric breakdown of electrically insulating material (as in traditional auto-magnetizing helical liners) and do not require helical return current geometries (as in dynamic screw pinches). We present results from three-dimensional radiation-magnetohydrodynamic simulations of MagLIF implosions employing AutoMag-D liner designs. Simulated AutoMag-D targets demonstrate improved fuel conditions and thermonuclear yield compared to traditional MagLIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Selecting Appropriate Model Complexity: An Example of Tracer Inversion for Thermal Prediction in Enhanced Geothermal Systems

Abstract A major challenge in the inversion of subsurface parameters is the ill‐posedness issue caused by the inherent subsurface complexities and the generally spatially sparse data. Appropriate simplifications of inversion models are thus necessary to make the inversion process tractable and meanwhile preserve the predictive ability of the inversion results. In this study, we investigate the effect of model complexity on fracture aperture inversion and thermal performance prediction in a field‐scale EGS model. Principal component analysis was used to map the aperture field to a low‐dimensional latent space. The complexity of the inversion model was quantitatively represented by the percentage of total variance in the original aperture fields preserved by the latent space. Tracer, pressure and flow rate data were used to invert for fracture aperture through an ensemble‐based inversion method, and the inferred aperture field was used to predict thermal performance. With an over‐simplified aperture model, ensemble collapse occurred. The inverted aperture models failed to resolve necessary flow and transport features, leading to a biased thermal performance prediction. A complex aperture model involved excessive features and was prone to overinterpreting the inversion data. Both the tracer/pressure/flow rate data reproduction and thermal prediction showed significant uncertainties, making it difficult to properly estimate long‐term thermal performance. Fortunately, our results indicate that there exists an appropriate model complexity which can simultaneously match inversion data and predict thermal performance with an acceptable uncertainty. The quality of the fit of tracer data appears to be a useful indicator of such an appropriate model complexity.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

A comparative analysis of numerical approaches for the description of gas flow in clay-based repository systems: From a laboratory to a large-scale gas injection test

There is nowadays a consensus among many countries that geological disposal is a favourable solution for the long-term management. Although different host formations and different barrier systems are under consideration around the world, clay-based materials form an important component for waste isolation in most national programmes. Hence, a good comprehension of the effect of gas flow on the hydro-mechanical behaviour of clay-based soils is essential, both at laboratory and field scale. Task B under the international cooperative project DECOVALEX-2023 has recently shown that, after some enhancement, models can be employed to reproduce laboratory scale tests, even with different sample geometries 37 . However, further work is required to understand whether they can be applied to simulate a large-scale experiment. Up-scaling of models for the advective transport of gas through clay-based low permeable material presents a number of problems related to the difficulty in obtaining consistent hydrogeological parameters and constitutive relationships at both laboratory and field scale. Based on a unique dataset from a large-scale gas injection test (Lasgit) performed at the Äspö Hard Rock Laboratory (Sweden), Task B within DECOVALEX-2023 has explored the refinement of these numerical strategies applied to the simulation of gas flow. Work performed within the task reveals that codes do not need to be substantially modified from the laboratory models to reproduce full-scale tests: indeed, model parameters calibrated and validated at laboratory scale have been applied to predict field scale gas flow at Lasgit, including peak gas pressure and injected cumulative gas volume. By means of (1) the introduction of interfaces between blocks to reflect the experimental configuration and the (2) adjustment of some parameters (e.g., higher permeability), the updated models are able to represent most of the key features observed in the experimental data, even at a large scale.

Tamayo-Mas, E↗
DOE OSTI2024

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Utah FORGE: Wells 16A(78)-32 and 16B(78)-32 Extended Circulation Test Data - August and September 2024

The dataset includes data collected during an extended circulation test conducted at the Utah FORGE site between August 8 and September 5, 2024. It provides uncorrected, raw digital data from the test, along with calibration information and a final report detailing the test procedures and corrections. The data is bundled in a single .zip file containing both field calibration and test data. Files include manual calibration data for temperature, pressure, and flow meters, as well as uncorrected, time-stamped measurements recorded in 30-second intervals, such as wellhead pressures, flow rates, and temperatures for wells 16A(78)-32 and 16B(78)-32. Additionally, the final report offers insights into the test methodology and provides guidelines on how to apply field calibrations to correct the raw Pason data.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Reactive Carbide‐Based Synthesis and Microstructure of NASICON Sodium Metal All Solid‐State Electrolyte

Reactive carbide precursor-based synthesis of NASICON-type NZSP (Na 1+x Zr 2 Si x P 3-x O 12 ) solid-state electrolyte (SSE) is demonstrated, in contrast to the established oxide-based approach. Exothermic decomposition of ZrC and SiC in air homogenizes microstructure, yielding 98% compact density after conventional sintering at 1200 °C. Quantitative stereology demonstrates that significant microstructural differences are present. Compacts of carbide-derived Carb-NZSP are 98% dense with a secondary zirconium oxide (ZrO 2 ) volume fraction of 0.2% ± 0.3%, versus 93% dense and 3% ± 1% for oxide-derived baseline. For Carb-NZSP, the secondary glassy phosphate phase is agglomerated, while for baseline, it is dispersed and percolated. Electrochemical testing combined with post-mortem analysis demonstrates how microstructural control of secondary phases is critical for dendrite suppression: Carb-NZSP critical current density (CCD) is 3.1 ± 0.8 mA cm −2 at 0.1 mAh cm −2 , versus 1.0 ± 0.7 mA cm −2 at 0.1 mAh cm −2 . Cryogenic focused ion beam (cryo-FIB) analysis demonstrates that in both materials, the porous 2D sheet-like sodium metal dendrites propagate around and subsume NZSP grains, likely following a path enriched with glassy phase and with porosity. Dendrites also flow around isolated zirconia particles. Phase field simulation reveals deflection of dendrites by mechanically tough zirconia, while brittle glassy phase accelerates dendrite growth, especially when finely distributed.

36 MATERIALS SCIENCE↗