Search NASA⌕ Search

SEARCH · Search NASA

Results for “alignment”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical framework for integrated additive manufacturing-compression molding (AM-CM) of thermoplastic composites

Additive manufacturing-compression molding (AM-CM) has emerged as a transformative technology in advanced composite manufacturing. Additive manufacturing (AM) offers high design flexibility and the ability to produce complex geometries with precisely aligned fibers in the preferred orientation. Compression molding (CM) enhances composite materials by providing excellent dimensional stability, reduced porosity, high production rates, and a smooth surface finish. Despite these advantages, extensive integrated analysis is required to optimize processing conditions for improved fiber orientation distribution (FOD) and porosity control. Here, this study develops a comprehensive numerical model to simulate the AM-CM manufacturing process. The model isolates the effects of both the AM and CM phases while also capturing their integration. Additionally, it accounts for heat transfer, temperature-dependent viscosity, and fiber orientation in the extruded fiber-filled polymer, accurately representing material behavior during processing. This approach enables the analysis of interactions between deposited beads of complex strand shapes and their interface regions after full compression. Moreover, the model predicts key parameters such as polymer flowability, fiber orientation, and temperature evolution in AM-CM parts. By optimizing processing conditions, it facilitates a controlled and predictable microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Initial assessment of alternative carbon fiber geometries for design of cost-effective compressive performance: Size effect studies

Carbon fiber provides opportunity to reduce weight in structural composites, including wind turbine blades, due to the material's superior specific stiffness and specific strength compared to alternatives. Despite these advantages, cost and compressive performance are considered weaknesses for carbon fiber products available today. Studies to produce low-cost carbon fiber alternatives, including the use of textile-derived precursor systems, have shown progress and merit through the DOE/ORNL low-cost carbon fiber initiatives. Here, this work focuses on enabling increases in compressive strength through design of the carbon fiber geometry, applicable to both textile and conventional precursor systems, while also providing opportunities to reduce carbon fiber processing costs. Fiber-resin interface and fiber alignment are among the most frequently cited factors controlling composite compressive performance. However, it is believed that there is opportunity in traditionally unexplored routes to increasing compressive strength through alteration of the carbon fiber geometry by increasing the fiber area moment of inertia and/or the fiber perimeter and interfacial area. This paper presents initial results from manufacturing carbon fiber materials to assess the impacts of carbon fiber size on tested composite compressive performance with projected neutral or even beneficial impact on fiber and composite manufacturing economics. Carbon fiber systems with increasing size illustrate a favorable correlation for compressive performance greater than predicted from a micromechanical failure model. The manufacturing and mechanical test results support the hypothesis of this work that alterations to fiber geometry can be used to produce improvements of the compressive strength of carbon fiber reinforced polymers and provide incentive for related work in designing alternative shapes to further enhance compressive performance.

42 ENGINEERING↗
DOE OSTI2025

Genesis of a novel high-rate composite manufacturing process using large-scale additive manufacturing – compression molding (AM-CM) system: Possibilities and limitations

Oak Ridge National Laboratory (ORNL) has developed a highly automated manufacturing process for thermoplastic composites that combines the benefits of Additive Manufacturing and Compression Molding (AM-CM) to produce high-performance functional composite structures at automotive production rates. Here, the AM-CM process creates highly precise preforms by additively placing extruded fiber-filled polymers (with controlled fiber orientations and multi-material configurations) in the desired mold location before undergoing a secondary compression molding process immediately before the preform cools down. Preforms can be in the form of short, long-chopped, or continuous fiber-filled thermoplastic polymers (e.g., CF/GF-filled ABS, PC, LM-PAEK, etc.). The AM-CM process combines the benefits of controlled fiber alignment, that is only achievable in AM-printed parts with the classical CM process, which eliminates porosity and good surface finish. A preform created using AM-CM can integrate various materials to enable additional architectural functionalities, including over-molding, selective stiffening, and the incorporation of electrically or thermally conductive channels. All these advantages come with a fast part production cycle time. The AM-CM process can manufacture multi-material, multi-functional parts in under 3 min, starting from raw material (pellets) to the final product. The novel AM-CM process offers superior microstructural control and enhanced multi-functionality previously unattainable with any other traditional high-rate thermoplastic composite manufacturing method. This work covers the manufacturing concept, system development, materials and applications of AM-CM process in detail.

Kumar, Vipin [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2024

Examining infrared thermography based approaches to rapid fatigue characterization of additively manufactured compression molded short fiber thermoplastic composites

A novel additive manufacturing (AM) methodology combined with a compression molding (CM) process has been developed to optimize the microstructure of short fiber thermoplastic composites (SFTs)with higher fiber alignment and lower porosity, yielding superior stiffness, strength, and structural integrity. Here, the current work examines the efficacy of the ‘passive’ infrared thermography (IRT) techniques for rapid fatigue characterization of SFTs that use the surface temperature evolution during cyclic loading due to self-heating as a fatigue indicator. A comparison of fatigue limits obtained from traditional stress-life (SN) (≈53.1%σ uts ) and IRT (≈54.1%σ uts ) shows a close match. However, the SN curve required 18 specimens and two weeks of continuous cyclic testing, while IRT used three specimens with 5 hours of testing. Thus, the IRT approach provides an accelerated testing framework for rapidly estimating the fatigue limit. Additionally, existing phenomenological approaches to IRT fatigue characterization have been examined.

42 ENGINEERING↗
DOE OSTI2024

SMATool: Strength of materials analysis toolkit

The study of the strength of materials is a cornerstone in material science and engineering, playing a critical role in shaping the progress and application of materials in diverse industrial sectors. The strength of a material is meticulously examined to understand the behavior of the material under different stress conditions and environments, thereby guiding material selection and structural design. Herein, we introduce the SMATool, a computational toolkit for the efficient calculation and analysis of material strength at both zero and finite temperatures for 3D, 2D, 1D, and tubular 2D-based nanostructures and nanotubes, as well as 1D nanoribbons. The toolkit is capable of calculating tensile, shear, ultimate, yield, and indentation (Vickers' hardness) strengths in various dimensions, as well as the energy storage capacity. We conducted several calculations both at zero and finite temperatures to validate the accuracy and reliability of the developed software. Here, the results show that the SMATool package provides accurate predictions that align with existing data on material strength. SMATool integrates seamlessly with widely used electronic structure codes like VASP and Quantum Espresso, providing a user-friendly interface catering to academic researchers and industry professionals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Directional finite difference method for directly solving 3D gyrokinetic field equations with enhanced accuracy

The gyrokinetic (GK) field equation is a three-dimensional (3D) elliptic equation, but it is often simplified to a set of two-dimensional (2D) equations by assuming that the field does not vary along a specific direction. However, this simplification can introduce inevitable 0th-order numerical errors, as nonlinear mode coupling in toroidal geometry can produce undesirable harmonic modes that violate the assumption. In this work, we propose a novel directional finite difference method (FDM) with a local coordinate transformation to better resolve the target field of interest. The directional FDM can accurately solve 3D GK field equations without simplifications, which can overcome the limitations of conventional methods. The accuracy and efficiency of different FDMs are analyzed in great detail for a variety of geometries, from simple 2D Cartesian coordinates to realistic 3D curvilinear coordinates. The 0th-order numerical errors of simplified 2D GK equations were found to be more problematic for low-harmonic modes and low aspect ratio geometries such as spherical tokamaks. On the other hand, the directional 3D FDM can accurately resolve a much wider range of harmonic modes aligned to the direction of interest, including the low-harmonic modes. In conclusion, we demonstrate that the directional 3D FDM is a highly effective algorithm for solving the 3D GK field equations, achieving accuracy improvements of 10 to 100 times or more, particularly for low-harmonic modes in spherical tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Parallel diffusion operator for magnetized plasmas with improved spectral fidelity

Diffusive transport processes in magnetized plasmas are highly anisotropic, with fast parallel transport along the magnetic field lines sometimes faster than perpendicular transport by orders of magnitude. This constitutes a major challenge for describing non-grid-aligned magnetic structures in Eulerian (grid-based) simulations. Here, the present paper describes and validates a new method for parallel diffusion in magnetized plasmas based on the anti-symmetry representation [Halpern and Waltz, Phys. Plasmas 25, 060703 (2018)]. In the anti-symmetry formalism, diffusion manifests as a flow operator involving the logarithmic derivative of the transported quantity. Qualitative plane wave analysis shows that the new operator naturally yields better discrete spectral resolution compared to its conventional counterpart. Numerical simulations comparing the new method against existing finite difference methods are carried out, showing significant improvement. In particular, we find that combining anti-symmetry with finite differences in diagonally staggered grids essentially eliminates the so-called “artificial numerical diffusion” that affects conventional finite difference and finite volume methods.

Anisotropic diffusion↗
DOE OSTI2026

Fast solvers for tokamak fluid models with PETSc

Multigrid (MG) is widely recognized as a highly effective solver for the model problem, the Laplacian, but textbook MG fails on most problems of interest. MG methods have been applied to complex, real-world applications with careful consideration of the physical model and discretization. In this work we develop the first step in applying MG methods to science and engineering relevant magnetohydrodynamics (MHD) tokamak models in the M3D-C1 (https://m3dc1.pppl.gov) fusion energy science code. The semi-implicit time integrator in M3D-C1 is composed of many linear solves. The implicit advance of the momentum equation is the most challenging and is the focus of this work. The current production solver in M3D-C1 is a block Jacobi (BJ) preconditioner within a Krylov solver, where blocks group degrees of freedom on planes of constant toroidal coordinate. BJ convergence degrades as the number of planes increases due to the spectral properties of the matrix preconditioned with BJ. The partially magnetic field-aligned, regular toroidal grid structure in M3D-C1 is amenable to semi-coarsening geometric MG in the toroidal direction. This paper develops such a solver and demonstrates competitive performance on a runaway electron model of a SPARC (https://cfs.energy/technology/sparc) disruption, and superior robustness on a stellarator model on which the BJ solver fails to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Predicting receptor-ligand pairing preferences in plant-microbe interfaces via molecular dynamics and machine learning

Microbiome assembly, structure, and dynamics significantly influence plant health. Secreted microbial signaling molecules initiate and mediate symbiosis by binding to structurally compatible plant receptors. For example, lipo-chitooligosaccharides (LCOs), produced by nitrogen-fixing rhizobial bacteria and various fungi, are recognized by plant lysin motif receptor-like kinases (LysM-RLKs), which activate the common symbiotic pathway. Accurately predicting these molecular interactions could reveal complementary signatures underlying the initial stages of endosymbiosis. Despite the breakthrough in protein-ligand structure prediction with deep learning-based tools, such as AlphaFold3, the large size and highly flexible nature of signaling compounds like LCOs present major challenges for detailed structural characterization and binding-affinity prediction. Typical structure-/physics-based methods of ligand virtual screening are designed for small, drug-like molecules, often rely on high-resolution, experimentally determined structures of the protein receptors, and rarely achieve sufficient sampling to obtain converged thermodynamic quantities with large ligands. In this study, we developed a hybrid molecular dynamics/machine learning (MD/ML) approach capable of predicting binding affinity rankings with high accuracy in systems involving large, flexible ligands, despite limited experimental structural information. Using coarse initial structural models, the predictions using the MD/ML workflow achieved strong alignment with experimental trends, particularly in the top-affinity tier for four legume LysM-RLKs (LYR3) binding to LCOs and a chitooligosaccharide. Furthermore, the MD-based conformation selection protocol provided critical structural insights into substrate specificity and binding mechanisms. This study demonstrates a powerful method to screen for challenging cognate ligand-receptors and advance our understanding of the molecular basis of microbial colonization in plants.

Lipo-chitooligosaccharides↗
DOE OSTI2025

Probing signals of self-interacting dark matter core collapse in H I -rich galaxies

Here, we analyze rotation curves of five Hi-rich galaxies recently discovered with MeerKAT. These galaxies exhibit sharply rising rotation curves, while their baryonic components are not dynamically dominant, suggesting that their dark matter halos have high inner densities. When fitting the standard Navarro–Frenk–White (NFW) halo model, four galaxies require extremely high halo concentrations, exceeding the cosmological median by 5σ. In contrast, self-interacting dark matter (SIDM) halos in the core-collapse phase naturally account for the high densities in these galaxies. For halos with masses around 10 11 M ⊙ , those in cosmic filaments exhibit concentrations consistent with the cosmological average, while halos in cosmic nodes show relatively higher concentrations that align with the SIDM fits but remain insufficient for the NFW fits. Our analysis indicates that these Hi-rich galaxies may have formed in cosmic nodes of dark matter with significant self-interactions.

dark matter↗
DOE OSTI2025

Characterization of thermally heat-treated polyacrylonitrile carbon fibers

This study investigates the graphitization process of polyacrylonitrile (PAN) carbon fibers by subjecting commercial fibers to thermal heat treatment at temperatures ranging from 1400 to 2100 °C in 100 °C increments, using either argon or nitrogen gas atmospheres. Changes in crystallinity, surface morphology, and lattice parameters were analyzed for two commercial carbon fibers using X-ray diffraction, scanning electron microscopy, and Raman spectroscopy. Results indicated minimal changes in surface morphology with increasing heat-treatment temperature; however, crystallinity significantly increased. Crystallinity changes were more strongly dependent on temperature rather than gas atmosphere or fiber type. At intermediate heat-treatment temperatures (1600–1800 °C), fibers treated in argon showed a slight preference for graphitization. The highest level of graphitization was measured at 2100 °C. Crystallite size increased as the intensity ratio of the D1 to G Raman peaks increased, reaching a peak around ~1800 °C, after which the ratio started to decrease. This behavior aligns with Ferrari's three-stage model of carbon crystallization and is consistent with both the Marie-Mering degree of graphitization and Brubaker's Integrated Absolute Differential models, all of which describe the transformation from an amorphous to a more graphitic structure. At the higher heat-treatment temperatures, the changes between atmospheres and fiber types were measured to converge to similar levels of graphitization. In conclusion, this study evaluates the progressive change in commercial grade carbon fibers when heat-treated.

Characterization↗
DOE OSTI2026

High-Resolution South American Wind Resource Data Downscaled with Generative Machine Learning Conditioned on Near-Surface Observations

High-resolution historical wind data was developed for the entirety of South America using the innovative Super-Resolution for Renewable Resource Data (sup3r) machine learning framework. The publicly available Sup3rWind South America dataset represents a significant advancement in wind resource data generation, leveraging generative machine learning conditioned on near-surface observations from the Meteorological Assimilation Data Ingest System (MADIS) to efficiently and accurately downscale coarse reanalysis data from the European Centre for Medium-Range Weather Forecasts (ERA5). This approach produces fine-scale, spatially and temporally coherent wind and meteorological fields hundreds of times more computationally efficient than traditional numerical weather modeling methods, enabling access to high-fidelity wind information across both continental and offshore regions. Sup3rWind South America builds on the earlier Sup3rWind Ukraine dataset through improvements in model architecture and outputs conditioned on near-surface observation inputs. As with the Ukraine data release, this dataset includes wind speed, wind direction, temperature, relative humidity, and pressure at a horizontal resolution of ~2 km, representing a 15x spatial enhancement relative to the 31 km ERA5 grid. Wind speed and direction are provided at 5-minute resolution, a 12x temporal refinement compared to the hourly ERA5 data, while temperature, relative humidity, and pressure remain at hourly resolution. The data covers all years from 2005 to 2024. Before downscaling, ERA5 inputs were bias-corrected using long-term monthly means and a limited number of quality-controlled observations to align large-scale statistics with regional conditions. The resulting dataset is the first publicly available high-resolution timeseries wind record that provides full spatial coverage of South America. Model validation demonstrates strong agreement with observations across several statistical metrics, consistent with other state-of-the-art high-resolution wind resource datasets. The potential applications of Sup3rWind South America span renewable energy resource assessment, energy system modeling, and grid resilience analysis. The 20-year record and high spatial and temporal resolution support accurate estimation of long-term energy yield and the economic feasibility of potential wind development sites. Continuous coverage across both continental and offshore regions enables comprehensive site prospecting within exclusive economic zones. The 2 km, 5-minute resolution data provide the spatial and temporal variability required for power system simulation, operational planning, and regional risk assessments.

17 WIND ENERGY↗
DOE OSTI2024

Advances in building data management for building performance standards using the SEED platform

Reducing energy consumption and greenhouse gas emissions in the built environment is a critical step in achieving emission goals to mitigate climate change impacts. Local, federal, and international jurisdictions are deploying several methods to reduce energy and emissions such as voluntary and mandatory benchmarking and building performance standards, requiring building owners to reach energy and emission targets. Jurisdictions leveraging benchmarking and building performance standards require knowledge of the buildings covered; which is a large task due to staffing constraints, limited information on building characteristics and tax parcel data, and the need for advanced data management techniques to align datasets. This paper describes an open-source platform's recent advances to create consistent taxonomies, identify erroneous data, enable auditability, and track building performance. The paper concludes with two use cases on how the platform has been used by jurisdictions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Critical review of lithium recovery from geothermal brines with implications for Smackover Formation, USA

The rapidly growing demand for lithium, a critical element for energy storage and national security technologies, has intensified concerns over the long-term availability and environmental impact of conventional lithium sources, such as hard-rock mining. To meet future demand, it is vital to explore unconventional resources that can provide sustainable domestic supplies. Geothermal brines, produced as a byproduct of geothermal energy generation, offer a promising alternative for lithium recovery by leveraging existing infrastructure and renewable energy production. In particular, the Smackover Formation, an extensive reservoir of high-salinity brines spanning Arkansas, Texas, Louisiana, Mississippi, and Alabama in the U.S. Gulf Coast, holds significant untapped lithium reserves. Co-producing geothermal energy and lithium from these brines aligns with sustainable extraction objectives while addressing resource scarcity. This review synthesizes current knowledge of lithium occurrence in the Smackover Formation and geothermal resources in the region, while also exploring how emerging tools such as machine learning can enhance resource targeting and co-production efficiency. Finally, we discuss key technical challenges and outline future research directions needed to advance lithium extraction from geothermal brines and secure a resilient domestic supply chain.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Comparison of structurally diverse simulation models for prediction of epidemic outcomes caused by a long-distance dispersed pathogen

Long-distance dispersal (LDD) pathogens pose substantial challenges for epidemic control due to their ability to generate new infection foci at great distances. While various modeling approaches have been developed to understand and manage such outbreaks, little work has compared how models of different structures behave under shared conditions. Here, in this study, we compare four structurally distinct epidemiological models — EPIMUL, GEMF, PoPS, and Warwick — each adapted to simulate the spread of wheat stripe rust (WSR), a wind-dispersed LDD pathogen, under identical epidemiological parameters and dispersal kernel. Using data from a controlled field experiment, we evaluate the ability of each model to replicate disease prevalence under nine intervention scenarios that vary in timing and culling area. While the models differ substantially in design — ranging from spatial grid-based to network-based and raster-based frameworks — the shared dispersal kernel allowed for close alignment in their predictions. All models accurately captured general epidemic trends, particularly the strong effect of early intervention on disease suppression. We qualitatively compared their behavioral responses across scenarios and also evaluated an ensemble prediction by averaging across model outputs. Our findings highlight how integrating shared epidemiological components into distinct modeling frameworks can improve consistency and accuracy, while reinforcing the importance of early culling in managing LDD pathogen outbreaks.

Dispersal kernel↗
DOE OSTI2026

Intelligent industrial demand response to increase grid flexibility and reliability: A review

The rapid transition toward renewable energy has introduced challenges in grid stability due to the intermittency of non-dispatchable sources like solar and wind. Industrial Demand Response (IDR) offers a promising, cost-effective solution that adjusts energy consumption patterns to align with supply, increases renewable utilization, and reduces costs. This review provides an updated analysis of IDR, sorting technologies into five categories: energy storage, scheduled energy usage, operational flexibility, on-site generation, and intelligent operations. Energy storage solutions, while requiring little flexibility, often have the longest payback periods. While slightly better, on-site generation also has longer payback periods, ranging from 5 to 20 years or more. Scheduled energy usage, operational flexibility, and intelligent operations allow significant peak reduction at lower capital costs but require greater flexibility. While 15–20 % peak reduction is within the range of all five categories, scheduled energy use and on-site energy generation are shown to have reductions of up to 70–80 % in select scenarios. Combining multiple IDR strategies from these five categories maximizes both financial and operational benefits. Synergistic approaches are shown to enhance grid stability while reducing costs. As the grid evolves, IDR will enable a more flexible, renewable-powered future that will benefit industrial facilities and the broader energy system.

Demand flexibility↗
DOE OSTI2024

In-situ polarization modulation IRRAS investigation of ammonia electrooxidation on Pt-Ir and Pt-Ru nanoparticles prepared on engineered catalyst supports

The catalytic activity and surface reactivity of monometallic Pt and bimetallic Pt-Ir and Pt-Ru nanoparticles, supported on two distinct Engineered Catalyst Supports (ECSs), were investigated for the Ammonia Electrooxidation Reaction (AmER) in alkaline media. XRD measurements confirmed alloy formation between Pt-Ir and Pt-Ru nanoparticles, as indicated by the shift of the (111) reflection to higher 2θ values. Cyclic voltammetry, linear sweep voltammetry, and chronoamperometry experiments were conducted to assess the catalytic activity of the Pt, Pt-Ir, and Pt-Ru electrocatalysts. All bimetallic catalysts exhibited lower onset potentials compared to Pt. The differing Tafel slopes between Pt (74 mV dec⁻¹), Pt-Ir (152 mV dec⁻¹), and Pt-Ru (118–197 mV dec⁻¹) suggest that alloying Pt with Ir or Ru alters the reaction mechanisms. Furthermore, the bimetallic Pt-Ir and Pt-Ru catalysts demonstrated greater tolerance for concentrated ammonia solutions relative to Pt. In-situ Polarization Modulation Infrared Reflection Absorption Spectroscopy (PM-IRRAS) provided insights into the formation of N-H species, azide anions (N₃⁻), and N-O compounds. For the Pt-Ru catalyst, an additional peak around ~3600 cm⁻¹ was observed, corresponding to OH⁻ species. The PM-IRRAS results align with the Gerischer–Mauerer mechanism, indicating that partially dehydrogenated ammonia adsorbates act as active intermediates in the oxidation of ammonia over Pt-Ir and Pt-Ru catalysts.

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

Enabling fast-charging of lithium-ion batteries through printed electrodes

It has been well recognized that introducing secondary porous networks (SPNs) into the electrodes can effectively improve the electrochemical performance of lithium-ion batteries (LIBs), especially under fast-charging operations. However, the process complexity and high cost limit the commercial success of advanced electrodes with SPNs. To address this issue, we developed a facile screen-printing process to produce structured graphite electrodes with SPNs. The experimental results demonstrated that, by tuning the diameter and center-to-center (C2C) distance of emulsion dots on the stencil screen, the pore diameters and C2C pore distances of SPNs in screenprinted electrodes can be precisely controlled in the range of 100 mu m to 1 mm and 100 mu m to 3 mm respectively. In addition, the SPNs with hexagonal and square-shape pore alignments have also been imprinted onto the electrode coatings through adjusting the patterns of screen stencils. Used as anodes, the printed graphite electrodes demonstrated significantly reduced overpotential and voltage fluctuation under fast-charging operations from 2C to 6C. Coupled with LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622) cathodes, the full cells with printed graphite anodes exhibited an unprecedently stable performance with almost no capacity decay up to 170 cycles when charged to 80 % SOC at 2C. Observations from electron microscopy showed plated lithium undetectable at the surface of printed graphite electrodes after numerous cycles. The electrochemical analysis on the voltage evolution during the cell rest period indicated the significantly delayed onset of lithium plating in the presence of printed graphite electrodes. In conclusion, all these results suggest that the significantly improved cell performance is associated with the shortened Li-ion diffusion distance, reduced polarization and suppressed Li plating in the printed electrodes with patterned SPNs.

25 ENERGY STORAGE↗