Search NASA⌕ Search

SEARCH · Search NASA

Results for “DECOUPLING”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

A comparative study of two numerical approaches for solving Kim–Kim–Suzuki phase-field models

Among the standard multi-phase multi-component phase-field (PF) methods, the Kim–Kim–Suzuki (KKS) method has the advantage of decoupling interfacial energy from bulk energy and solving concentration as the conserved variable. There are two approaches to numerically solving a KKS method: the global solution approach (GSA) solves all variables in a global system simultaneously, and the local solution approach (LSA) solves phase concentrations locally using a Newton solver. This work compares the performance of LSA and GSA for solving four KKS models of increasing complexity with the finite element method using the MOOSE framework. The solution accuracy, degrees of freedom (DOFs), memory usage, and computational efficiency are compared. We find that GSA and LSA generate similar solutions, with a maximum difference of only 0.34%. For each model, LSA has a lower number of DOFs, utilizes less memory, and less wall time. Additionally, the savings of memory and wall time in LSA increase with increasing mesh density of the same model and are more pronounced in models with higher dimensionality and more nodes. However, GSA is easier to implement in existing codes and can better solve highly nonlinear systems by utilizing sophisticated solvers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Application of a temporal multiscale method for efficient simulation of degradation in PEM Water Electrolysis under dynamic operating conditions

Hydrogen is emerging as a vital energy carrier, driven by the need to reduce carbon emissions. Proton Electrolyte Membrane Water Electrolysis (PEMWE) enables hydrogen production under fluctuating renewable power conditions but requires improved understanding and stability of the anode catalyst layer under dynamic operating conditions, especially with low noble metal loadings. Long-term degradation experiments are both time-consuming and costly; therefore, a systematic, model-aided approach is essential. In the present work, a temporal multiscale method is applied to reduce the computational effort of simulating long-term degradation processes in PEMWE, with an exemplary focus on catalyst dissolution. A mechanistic model incorporating the oxygen evolution reaction, catalyst dissolution, and hydrogen permeation from the cathode to the anode was hypothesized and implemented. In this way, the local periodicity of transport and reaction processes in dynamic PEMWE operation, which influence the gradual degradation of the catalyst layer, is captured. The temporal multiscale method significantly reduces the computational effort of simulation, decreasing processing time from hours to mere minutes. This efficiency gain is attributed to the limited evolution of Slow-Scale variables during each period of time P of the Fast-Scale variables. Consequently, simulation is required only until local periodicity is achieved within each Slow-Scale time step. Hence, the fully resolved dynamic problem is decoupled into these two scales, employing a heterogeneous multiscale technique. The developed approach effectively accelerates parameter estimation and predictive simulations, supporting systematic modeling of PEMWE degradation under dynamic conditions.

08 HYDROGEN↗
DOE OSTI2026

Phenomena-based graph representations and applications to chemical process simulation

Rapid and robust simulation of chemical processes is critical to conduct process design, optimization, techno-economic analysis, and sustainability analysis. Yet, efficiently solving simulation models remains a challenge due to the highly coupled and nonlinear nature of the underlying algebraic equations that capture the physical phenomena taking place in the process (e.g., material and energy conservation, phase equilibrium, reactions). In this work, we show that graph-theoretic representations of the physical phenomena within unit operations can help navigate and decompose equations to systematically identify alternative approaches for fast and robust numerical solutions. Specifically, we present a graph-theoretic abstraction that captures the connectivity between the model variables/equations and use this abstraction to group variables/equations into fundamental phenomena. We show that phenomena-based decomposition of the underlying equations can help decouple nonlinearities and enforce material/energy conservation at the process level to accelerate convergence. The proposed decomposition approach differs from the more traditional sequential modular simulation approach, in which equations are grouped and decomposed by unit operations. We implemented the phenomena-based decomposition in BioSTEAM—an open-source process simulation platform in Python—and demonstrated that this approach can converge a variety of separation process models. Compared to sequential modular simulation, the phenomena-based approach can converge idealized systems faster, but it can be slower for (or even fail to converge) highly coupled and nonideal process systems.

Convergence↗
DOE OSTI2024

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Temperature dependent thermoelectric properties of cuprous delafossite oxides

The use of nanostructured delafossite oxides in thermoelectric (TE) applications has attracted a great interest due to their high performance and long-term stability at elevated temperatures. Cuprous delafossites, CuMO 2 (M=Al, Cr, Fe, Ga, Mn), compared to conventional TE materials, such as Bi 2 Te 3 , PbTe and SiGe, are non-toxic and more earth abundant. In particular, CuAlO 2 compound shows a great potential for high performance thermoelectric materials. In this work, a systematic study of temperature dependent TE properties of cuprous delafossite materials, CuAlO 2 , is reported. The optimization of the TE properties has been realized by controlling nanostructure size around 80 nm CuAlO 2 powder was prepared using a solid-state synthesis method at ~1373 K in nitrogen/air atmosphere. The nanostructure size was controlled by a high energy ball milling process. Reducing the particle size of nanostructured bulk materials decouples interdependent electron and phonon transport and results in a lattice thermal conductivity decrease without deteriorating electrical conductivity. As a result, the high effective mass plays a dominant role in the high Seebeck coefficient and low electrical conductivity. The power factor reached ~0.78×10 -5 W/mK 2 at 780 K. Temperature dependent TE properties, including Seebeck coefficient, electrical conductivity, and thermal conductivity are analyzed. The processing-structure-property correlation of these materials are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

An improved methodology for modeling short pulse buried layer x-ray emission spectra

Radiation-hydrodynamic and spectroscopic modeling are important aspects of high energy density experimental design. In this paper, we improve the performance and capabilities over those obtainable with a previous methodology used for simulating x-ray emission spectra from buried layer targets heated by short pulse lasers. The improvement incorporates post-processing HYDRA radiation-hydrodynamic output with a non-local thermodynamic equilibrium atomic-kinetics radiation transport code, Cretin. Each code uses an independent radiation field which allows decoupling HYDRA's radiation group structure from Cretin's spectral output to improve the speed and flexibility of the design methodology. The execution time decreases from 2–3 days to a few hours while the flexibility of the improved methodology allows for performing sensitivity studies including a comparison of steady-state and time-dependent atomic kinetics and differences in the radiation group structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A universal implementation of radiative effects in neutrino event generators

Due to the similarities between electron-nucleus (eA) and neutrino-nucleus scattering (νA), eA data can contribute key information to improve cross-section modeling in eA and hence in νA event generators. However, to compare data and generated events, either the data must be radiatively corrected or radiative effects need to be included in the event generators. We implemented a universal radiative corrections program that can be used with all reaction mechanisms and any eA event generator. Our program includes real photon radiation by the incident and scattered electrons, and virtual photon exchange and photon vacuum polarization diagrams. It uses the “extended peaking” approximation for electron radiation and neglects charged hadron radiation. This method, validated with GENIE, can also be extended to simulate νA radiative effects. This work facilitates data-event-generator comparisons used to improve νA event generators for the next-generation of neutrino experiments. Program Title: emMCRadCorr CPC Library link to program files:https://doi.org/10.17632/hmsxg82vnf.1 Developer's repository link:https://github.com/e4nu/emMCRadCorr Licensing provisions: AGPLv3 Programming language:C++ Nature of problem: Radiative effects can significantly modify the event kinematics and the resulting cross-sections. Such effects must be accounted for when comparing event generators to eA data. Existing radiative correction codes are tailored to specific processes and topologies, and are limited to a restricted phase space defined by the spectrometer acceptance. Therefore, a more general approach is required to apply radiative corrections to semi-inclusive and exclusive eA measurements. Solution method: Our program incorporates real photon radiation from both the incident and scattered electrons, as well as virtual photon exchange and photon vacuum polarization effects. It employs the “extended peaking” approximation for electron radiation while neglecting contributions from charged hadron radiation. The code is fully decoupled from event generator codes and can be used for all event generators in the market.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗
DOE OSTI2022

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗
DOE OSTI2025

Enhancing EnergyPlus capabilities to model dynamic building envelopes using python plugin

Nearly half of the energy consumption in the United States is related to buildings, resulting in an urgent need to develop innovative technologies to improve building energy efficiency. Dynamic building envelopes, comprising switchable insulation and thermal energy storage materials, have been proposed recently as a promising solution to reduce buildings' heating and cooling loads by thermally coupling the indoor environment with the ambient environment when beneficial while decoupling them when outdoor conditions are not favorable. Although various related technologies are still underway, the whole-building energy modeling tools, like EnergyPlus, do not have the capability to simulate the transient and dynamic nature of dynamic envelope materials and components to accurately predict their impact on building energy use. The objective of this study is to formulate a method in EnergyPlus simulation engine to model multilayer envelopes, comprising dynamic building materials with variable thermophysical properties, and discuss the changes made to the program using a Python plugin. Furthermore, the thermal performance of the dynamic envelopes using the proposed method is compared and verified with the results from a well-established commercial code, COMSOL Multiphysics. A parametric assessment is also conducted to evaluate the energy efficiency benefits of dynamic envelopes in a single-family residential building, demonstrating total annual energy savings up to 11.6 %, when a dynamic envelope operates alone, and up to 18.2 % when it is combined with a thin layer of phase change material as a thermal storage medium. Finally, a United States wide energy efficiency assessment is presented to showcase the geographical spread of the energy savings. The method designed and implemented in this study provides the researchers with the ability to implement their dynamic insulation methods in EnergyPlus and evaluate the whole building energy impact.

25 ENERGY STORAGE↗
DOE OSTI2026

Physics vs structure: A systematic benchmark of learning strategies for multi-zone building thermal dynamics

Recent advances in physics-informed and data-driven machine learning promise improved thermal models for advanced building control, yet there is limited quantitative evidence on when added physics structure and architectural complexity are beneficial. Here, this work presents a systematic benchmark of five representative system identification methods for modeling multi-zone building thermal dynamics: linear state-space models, multi-layer perceptrons, neural state-space models, neural ordinary differential equations, and physically-consistent neural networks. The methods are evaluated across multiple data regimes and zone coupling strategies. Using a high-fidelity multi-zone commercial building emulator, we examine short-term and long-term prediction accuracy, computational efficiency, and ease of development. Our results reveal critical trade-offs between prediction performance, model complexity, and physical consistency. We demonstrate that decoupled, nonlinear black-box models consistently outperform coupled physics-constrained architectures in both predictive accuracy and out-of-distribution robustness in majority of the test cases for the building type considered in the study. Our findings quantify the cost of complexity in building thermal modeling and provide concrete, actionable, scenario-based guidelines for selecting model classes for control-oriented applications.

Building thermal modeling↗
DOE OSTI2021

Energy-efficient sorption-based gas clothes dryer systems

Standard electric resistance and fuel-driven dehydration technologies exhibit a maximum coefficient of performance of well below 1 mainly due to enthalpy losses associated with the air leaving the dehydration system. To improve energy efficiency, condensing dryer systems condense the moisture captured from a product in a closed-loop air circulation cycle. Existing condensing dehydration systems including heat pump dryers, however, need to significantly cool the air to achieve dehumidification. The added cooling and subsequent heating to return the air to a desired drying temperature consume substantial energy and thus reduce drying performance. Here, an innovative sorption-based gas dehydration system is proposed to overcome barriers deteriorating energy efficiency in existing gas, electric, or heat pump dryer systems. Decoupling latent and sensible loads, the system employs a liquid-desiccant solution to directly capture air humidity, thereby allowing circulation of the air in a closed loop to achieve high drying energy efficiency. In other words, the system captures waste latent heat from the moisture produced during the dehydration process and reuses it to improve energy efficiency. This study focuses on a comprehensive quasi-steady-state thermodynamic modeling of the proposed sorption-based dehydration concept employed for a gas clothes dryer application to predict transient response and overall drying performance (i.e., time and energy metrics). The analysis indicates the proposed sorption-based gas clothes dryer system can deliver a specific moisture extraction rate of 1.71 kg of water per kWh (i.e., a combined energy factor of 3.167 kg (6.98 lbm) of dry cloth per kWh) with a drying time of 44 min. This is a 112% energy improvement compared with state-of-the-art gas clothes dryers exhibiting a combined energy factor of 1.50 kg (3.3 lbm) of dry cloth per kWh. The technology pursued here can potentially be employed as a platform for many fuel-driven equipment to take advantage of available waste thermal energy in the environment instead of simply burning a fuel.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

The reservoir network: A new network topology for district heating and cooling

Thermal district networks are effective solutions to substitute fossil fuels with renewable energy sources for heating and cooling. Moreover, thermal networking of buildings allows energy efficiency to be further increased. The waste heat from cooling can be reused for heating in thermal district systems. Because of bidirectional energy flows between prosumers, thermal networks require new hydraulic concepts. In this work, we present a novel network topology for simultaneous heating and cooling: the reservoir network. The reservoir network is robust in operation due to hydraulic decoupling of transfer stations, integrates heat sources and heat sinks at various temperature levels and is flexible in terms of network expansion. Furthermore, we used Modelica simulations to compare the new single-pipe reservoir network to a basecase double-pipe network, taking yearly demand profiles of different building types for heating and cooling. The electric energy consumed by the heat pumps and circulations pumps differs between the reservoir and base case networks by less than 1%. However, if the reservoir network is operated with constant instead of variable mass flow rate, the total electrical consumption can increase by 48% compared to the base case. As with any other network topology, the design and control of such networks is crucial to achieving energy efficient operation. Investment costs for piping and trenching depend on the district layout and dimensioning of the network. If a ring layout is applied in a district, the reservoir network with its single-pipe configuration is more economical than other topologies. For a linear layout, the piping costs are slightly higher for the reservoir network than for the base case because of larger pipe diameters.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Experimental study on a novel three-phase absorption thermal battery with high energy density applied to buildings

Buildings contribute to 75% of total electricity consumption and 80% of peak electricity demand in the United States. A significant portion of building electricity consumption goes to space conditioning. Therefore, using thermal energy storage (TES) to decouple the building electricity consumption and the fluctuating space-conditioning load has significant potential to shift peak electricity demand and improve the stability of the grid. To provide effective TES within the spatial constraints in buildings, a three-phase absorption thermal battery (TATB) system with very high energy storage density was designed and studied. The TATB stores low-temperature heat in hydrate crystals of salts and supports space conditioning through absorption systems. In this study, a benchtop TATB using lithium chloride (LiCl) hydrate crystals as the storage material was tested under practical operating conditions for its energy storage density and discharge rate. We report the TATB successfully generated LiCl hydrate crystals during charging mode and then dissolved the crystals during discharging mode. A material-based energy storage density of 300 kWh/m 3 (specific energy of 903 kJ/kg) was achieved with a discharge rate of up to 1.3 kW thermal energy. The test results verified the high performance of the TATB and encouraged further development of this technology.

25 ENERGY STORAGE↗
DOE OSTI2022

A fast and accurate modeling approach for water and steam thermodynamics with practical applications in district heating system simulation

In U.S. district heating (DH) systems, steam is the most common heat transport medium. Industry demand for new advanced modeling capabilities of complete steam DH systems is increasing; however, the existing models for water/steam thermodynamics are too slow for large system simulations because of computationally expensive algebraic loops that require the solution to nonlinear systems of equations. For practical applications, this work presents a novel split-medium approach that implements numerically efficient liquid water models alongside various water/steam models, breaking costly algebraic loops by decoupling mass and energy balance equations. New component models for steam DH systems are also presented. We implemented the models in the equation based Modelica language and evaluated accuracy and computing speed across multiple scales: from fundamental thermodynamic properties to complete districts featuring 10 to 200 buildings. Compared to district models with the IF97 water/steam model and equipment models from the Modelica Standard Library, the new implementation improves the scaling rate for large districts from cubic to quadratic with negligible compromise to accuracy. Additionally, for an annual simulation with 180 buildings, this translates to a computing time reduction from 33 to 1-1.5 h. These results are critically important for industry practitioners to simulate steam DH systems at large scales.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Anomalously high elastic modulus of a poly(ethylene oxide)-based composite electrolyte

The practical use of lithium metal anodes in solid-state batteries requires a polymer membrane with high lithium-ion conductivity, thermal/electrochemical stability, and mechanical strength. The primary challenge is to effectively decouple the ionic conductivity and mechanical strength of the polymer electrolytes. We report a remarkably facile single step synthetic strategy based on in-situ crosslinking of poly(ethylene oxide) (xPEO) in the presence of a woven glass fiber (GF). Such a simple method yields composite polymer electrolytes (CPE) of anomalously high elastic modulus up to 2.5 GPa over a broad temperature range (20 °C – 245 °C) that has never been previously documented. An unsupervised machine learning algorithm, K-mean clustering analysis, was implemented on the hyperspectral Raman mapping at the xPEO/GF interface. Using such a unique means, we show for the first time that the promoted mechanical strength originates from xPEO and GF interactions through dynamic hydrogen and ionic bonding. High ionic conductivity is achieved by the addition plasticizer (e.g. tetraglyme), where trifluoromethanesulfonate anions are tethered to the xPEO matrix and Li + cations are favorably transported through coordination with the plasticizer. Further, stringent galvanostatic cycling tests indicates the CPE can be stably cycled for >3000 h in a Li-metal symmetric cell at a moderate temperature (nearly 1500 Coulombs/cm 2 Li equivalents), outperforming most of the PEO-based electrolytes. The GF reinforced CPE reported here has multifunctional uses, such as solid electrolytes for all solid-state batteries and membranes for redox-flow batteries. Although the focus of this study is on lithium-based batteries, the results are equally promising for other alkali metal based batteries such as sodium and potassium.

25 ENERGY STORAGE↗