Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical mixing”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Potential Vorticity Mixing in a Tangled Magnetic Field

A theory of potential vorticity (PV) mixing in a disordered (tangled) magnetic field is introduced. The analysis is in the context of β -plane MHD, with a special focus on the physics of momentum transport in the stably stratified, quasi-2D solar tachocline. A physical picture of mean PV evolution by vorticity advection and tilting of magnetic fields is proposed. In the case of weak field perturbations, quasi-linear theory predicts that the Reynolds and magnetic stresses balance as turbulence Alfvénizes for a larger mean magnetic field. Jet formation is explored quantitatively in the mean field–resistivity parameter space. However, since even a modest mean magnetic field leads to large magnetic perturbations for large magnetic Reynolds number, the physically relevant case is that of a strong but disordered field. We show that numerical calculations indicate that the Reynolds stress is modified well before Alfvénization—i.e., before fluid and magnetic energies balance. To understand these trends, a double-average model of PV mixing in a stochastic magnetic field is developed. Calculations indicate that mean-square fields strongly modify Reynolds stress phase coherence and also induce a magnetic drag on zonal flows. The physics of transport reduction by tangled fields is elucidated and linked to the related quench of turbulent resistivity. We propose here a physical picture of the system as a resisto-elastic medium threaded by a tangled magnetic network. Applications of the theory to momentum transport in the tachocline and other systems are discussed in detail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Mixed-Flux Techniques for Rational Synthesis and Structural Control in Silver Chalcogenides

The functionality of materials is intrinsically linked to their structures, an axiom encapsulated in the principle of structure-property relationships. The pinnacle of materials design is the tailoring of its structure for a specific function, which requires the ability of rational synthesis and the development of synthesis science. This idea, however, remains elusive for the synthesis of complex extended solids. A major obstacle is the difficulty in using established chemical principles selectively to control reaction paths and favor certain structural patterns over numerous other possible results. In this context, we are developing a synthesis science approach that facilitates the control of the structure and bonding to create new structures. This is achieved by employing a two-component flux consisting of mixed hydroxides and halides as the reaction medium. This enables reaction conditions that allow better control of the structure dimensionality and composition by manipulating the temperature and solvent basicity (via the flux component ratio). Here, we demonstrate the efficacy of this method in controlling their structural motifs to arrive at 23 unreported compositions and 6 unique structure types. These materials are expected to exhibit a broad range of properties, from metallic to semiconducting, with calculations suggesting the potential for emergent phenomena such as Dirac semimetals. The reaction paths afforded by these mixed fluxes establish a direct correlation between the synthetic variables and properties, providing significant insight into a broadly applicable approach for new materials.

Zhou, Xiuquan [Argonne National Laboratory (ANL), ↗
DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗
DOE OSTI2022

Stable and unstable supersonic stagnation of an axisymmetric rotating magnetized plasma

The Naval Research Laboratory ‘Mag Noh problem’, described in this paper, is a self-similar magnetized implosion flow, which contains a fast magnetohydrodynamic (MHD) outward propagating shock of constant velocity. We generalize the classic Noh ( OSTI Tech. Rep. 577058, 1983) problem to include azimuthal and axial magnetic fields as well as rotation. Our family of ideal MHD solutions is five parametric, each solution having its own self-similarity index, gas gamma, magnetization, the ratio of axial to the azimuthal field and rotation. While the classic Noh problem must have a supersonic implosion velocity to create a shock, our solutions have an interesting three-parametric special case with zero initial velocity in which magnetic tension, instead of implosion flow, creates the shock at $t=0+$ . Our self-similar solutions are indeed realized when we solve the initial value MHD problem with the finite volume MHD code Athena. We numerically investigated the stability of these solutions and found both stable and unstable regions in parameter space. Stable solutions can be used to test the accuracy of numerical codes. Unstable solutions have also been widely used to test how codes reproduce linear growth, transition to turbulence and the practically important effects of mixing. Now we offer a family of unstable solutions featuring all three elements relevant to magnetically driven implosions: convergent flow, magnetic field and a shock wave.

Mechanics↗
DOE OSTI2024

On the accuracy of compressibility transformations

This study highlights the importance of satisfying the eddy viscosity equivalence below the logarithmic layer, to deriving accurate compressibility transformations. First, we analyze the ability of known transformations to satisfy the eddy viscosity equivalence and show that the accuracy of these transformations is strongly dependent on this ability. Second, in a step-by-step manner, we devise new transformations that satisfy this hypothesis. An approach based on curve fitting of the incompressible Direct Numerical Simulation data for eddy viscosity profiles below the logarithmic layer provides an extremely accurate transformation, which motivates self-contained methods, making use of mixing length formulas in the inner region. It is shown that the accuracy of existing transformations can be significantly improved by applying these ideas, below the logarithmic layer. Motivated by the effectiveness of the formulations derived from eddy viscosity equivalence, we introduce a new integral transformation based on Reynolds number equivalence between compressible and incompressible flows. This approach is based on defining a new compressible velocity scale, which affects the accuracy of transformations. Several choices for the velocity scale are tested, and in each attempt, it is shown that the eddy viscosity equivalence plays a very important role for the accuracy of compressibility transformations.

42 ENGINEERING↗
DOE OSTI2021

Resilience-Motivated Distribution System Restoration Considering Electricity-Water-Gas Interdependency

A major outage in the electricity distribution system may affect the operation of water and natural gas supply systems, leading to an interruption of multiple services to critical customers. Therefore, enhancing resilience of critical infrastructures requires joint efforts of multiple sectors. In this paper, a distribution system service restoration method considering the electricity-water-gas interdependency is proposed. The objective is maximizing the supply of electricity, water, and gas to critical customers after an extreme event. The operational constraints of electricity, water, and natural gas networks are considered. Additionally, the characteristics of electricity-driven coupling components, including water pumps and gas compressors, are also modeled. Relaxation techniques are applied to non convex constraints posed by physical laws of those networks. Consequently, the restoration problem is formulated as a mixed-integer second-order cone program, which can readily be solved by the off-the-shelf solvers. The proposed method is validated by numerical simulations on an electricity-water-gas integrated system, developed based on benchmark models of the subsystems. The results indicate that considering the interdependency refines the allocation of limited generation resources and demonstrate the exactness of the proposed convex relaxation

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Effects of Testosterone on Mixed-Muscle Protein Synthesis and Proteome Dynamics During Energy Deficit

Effects of testosterone on integrated muscle protein metabolism and muscle mass during energy deficit are undetermined. The objective was to determine the effects of testosterone on mixed-muscle protein synthesis (MPS), proteome-wide fractional synthesis rates (FSR), and skeletal muscle mass during energy deficit. This was a randomized, double-blind, placebo-controlled trial. The study was conducted at Pennington Biomedical Research Center. Fifty healthy men. The study consisted of 14 days of weight maintenance, followed by a 28-day 55% energy deficit with 200 mg testosterone enanthate (TEST, n = 24) or placebo (PLA, n = 26) weekly, and up to 42 days of ad libitum recovery feeding. Mixed-MPS and proteome-wide FSR before (Pre), during (Mid), and after (Post) the energy deficit were determined using heavy water (days 1-42) and muscle biopsies. Muscle mass was determined using the D3-creatine dilution method. Mixed-MPS was lower than Pre at Mid and Post (P < 0.0005), with no difference between TEST and PLA. The proportion of individual proteins with numerically higher FSR in TEST than PLA was significant by 2-tailed binomial test at Post (52/67; P < 0.05), but not Mid (32/67; P > 0.05). Muscle mass was unchanged during energy deficit but was greater in TEST than PLA during recovery (P < 0.05). The high proportion of individual proteins with greater FSR in TEST than PLA at Post suggests exogenous testosterone exerted a delayed but broad stimulatory effect on synthesis rates across the muscle proteome during energy deficit, resulting in muscle mass accretion during subsequent recovery.

Endocrinology & Metabolism↗
DOE OSTI2024

Interface-aware Methods for Diffusion Equations [Slides]

In this talk, we will present interface-aware methods for diffusion equation. The focus is on the necessity of interface-aware, mixed cells, interface reconstruction, discontinuity of material diffusion coefficients, order of accuracy in time, property for large time steps, and simpler interface-aware methods. Numerical examples are provided to demonstrate the features on these aspects of the methods

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Topological Gauge Actions on the Lattice as Overlap Fermion Determinants

Overlap fermion on the lattice has been shown to properly reproduce topological aspects of gauge fields. In this paper, we review the derivation of Overlap fermion formalism in a torus of three space-time dimensions. Using the formalism, we show how to use the Overlap fermion determinants in the massless and infinite mass limits to construct different continuum topological gauge actions, such as the level-k Chern–Simons action, “half-CS” term and the mixed Chern–Simons (BF) coupling, in a gauge-invariant lattice UV regulated manner. Taking special Abelian and non-Abelian background fields, we demonstrate numerically how the lattice formalism beautifully reproduces the continuum expectations, such as the flow of action under large gauge transformations.

2 + 1 dimensions↗
DOE OSTI2024

Kaon mixing beyond the standard model with physical masses

We present nonperturbative results for beyond the standard model kaon mixing matrix elements in the isospin symmetric limit ( m u = m d ) of QCD, including a complete estimate of all dominant sources of systematic error. Our results are obtained from numerical simulations of lattice QCD with N f = 2 + 1 flavors of dynamical domain wall fermions. For the first time, these quantities are simulated directly at the physical pion mass m π ∼ 139 MeV for two different lattice spacings. We include data at three lattice spacings in the range a = 0.11 – 0.07 fm and with pion masses ranging from the physical value up to 450 MeV. Compared to our earlier work, we have added both direct calculations at physical quark masses and a third lattice spacing making the removal of discretization effects significantly more precise and eliminating the need for any significant mass extrapolation beyond the range of simulated data. We renormalize the lattice operators nonperturbatively using RI-SMOM off-shell schemes. These schemes eliminate the need to model and subtract nonperturbative pion poles that arises in the RI-MOM scheme and, since the calculations are performed with domain wall fermions, the unphysical mixing between chirality sectors is suppressed. Our results for the bag parameters in the MS ¯ scheme at 3 GeV are B K ≡ B 1 = 0.5240 ( 17 ) ( 54 ) , B 2 = 0.4794 ( 25 ) ( 35 ) , B 3 = 0.746 ( 13 ) ( 17 ) , B 4 = 0.897 ( 02 ) ( 10 ) and B 5 = 0.6882 ( 78 ) ( 94 ) , where the first error is from lattice uncertainties and the second is the uncertainty due to the perturbative matching to MS ¯ . Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics

The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Predicting Frigid Mixed-Phase Clouds for Pristine Coastal Antarctica

Supercooled water is common in the clouds near coastal Antarctica and occasionally occurs at temperatures at or below -30°C. Yet the ice physics in most regional and global numerical models will glaciate out these clouds. This presents a challenge for the simulation of highly supercooled clouds that were observed at McMurdo, Antarctica during the Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE) project during 2015–2017. The polar optimized version of the Weather Research and Forecasting model (Polar WRF) with the recently developed two-moment P3 microphysics scheme was used to simulate observed supercooled liquid water cases during March and November 2016. Nudging of the simulations to observed rawinsonde profiles and Antarctic automatic weather station observations provided increased realism and much greater cloud water amounts. Furthermore, sensitivity tests that adjust the ice physics for extremely low ice nucleating particle (INP) concentrations decrease cloud ice and increases the cloud liquid water closer to observed amounts. In these tests, a liquid layer near cloud top is simulated, in agreement with observations. Accurate representation of INP concentrations appears to be critical for the simulation of coastal Antarctic clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

I.1 Heavy-Duty Diesel Combustion (Sandia National Laboratories)

Regulatory drivers and market demands for lower pollutant emissions, lower carbon dioxide emissions, and lower fuel consumption motivate the development of cleaner and more fuel-efficient engine operating strategies. Most current production heavy-duty diesel engines use a combination of both in-cylinder and exhaust emissions-control strategies to achieve these goals. The emissions and efficiency performance of in-cylinder strategies depend strongly on flow and mixing processes that can be influenced by using multiple fuel injections. Past work performed under this project showed that adding a second injection can reduce soot to levels below what would have been produced by an unchanged first injection, thereby increasing load while decreasing soot and potentially reducing brake specific fuel consumption. Information characterizing the important in-cylinder processes with multiple injections has been gleaned from ensemble-averaged planar laser-induced incandescence (PLII) imaging visualizing the soot cloud and planar induced fluorescence (PLIF) of OH characterizing the soot oxidation regions. PLII showed a consistent disruption of the first injection soot cloud by the second injection. In conjunction with OH-PLIF, differences in soot oxidation patterns for multiple injections compared to single injections were observed. This understanding was further enhanced in FY20, when high-speed imaging resolving the above-mentioned effects in a single cycle were combined with direct numerical simulations investigating the multiple-injection ignition process on the microscopic level of turbulence and chemistry interaction. In FY21, these findings in conjunction with findings from other researchers published in the scientific literature were composed into a preliminary multiple-injection conceptual model of fuel-mixing, injection and ignition processes. Remaining key research questions were also highlighted. In addition, wall heat flux was investigated experimentally and with numerical simulations to understand the potential of multiple injections to reduce the engine heat losses and further enhance the efficiency.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Refueling infrastructure planning in intercity networks considering route choice and travel time delay for mixed fleet of electric and conventional vehicles

The range anxiety has been a major factor that affects the market acceptance of electric vehicles. Even with the recent development of battery technologies, a lack of charging stations and range anxiety are still significant concerns, specifically for intercity trips. This calls for more investments in building charging stations and advancing battery technologies to increase the market share of electric vehicles and improve sustainability. This study suggests a configuration for plug-in electric vehicle charging infrastructure to support long-distance intercity trips of electric vehicles at the network level. A model is proposed to minimize the total system cost including infrastructure investment (building charging stations/spots) and travel time delays (charging time, waiting time in the queue, and detour time to access charging stations). This study fills existing gaps in the literature by capturing realistic patterns of travel demand and considering flow-dependent charging delays at charging stations. Furthermore, the proposed model, which is formulated as a mixed-integer program with nonlinear constraints, solves the optimization problem at the network level. At the network level, impacts of charging station locations on the traffic assignment problem with a mixed fleet of electric and conventional vehicles need to be considered. To this end, a traffic assignment module is integrated with a simulated annealing algorithm. The numerical experiments show a satisfactory application of the model for a full-scale case study (intercity network in Michigan). The solution quality and efficiency of the proposed solution algorithm are evaluated against those of an enumeration approach for a small case study. The results suggest that even for the current market share and charging stations’ setting, a significant investment is needed to support intercity trips without range anxiety issues and with acceptable delays. Additionally, through sensitivity analyses, the required infrastructure and battery investments to support intercity trips with acceptable delays are established for hypothetical increased market shares and battery size in the future.

42 ENGINEERING↗
DOE OSTI2024

Modeling Flow and Particle Deposition in a Spent Nuclear Fuel Assembly

CFD (Computational Fluid Dynamic) simulation of aerosol-laden natural convective flow and particle deposition in a spent fuel storage canister with 37 assemblies is currently computationally prohibitive. PWR (Pressurized Water Reactor) assemblies have up to 289 pins or tubes with several spacer grids to align the pins. Spacer grids with mixing vanes induce swirling during operation to increase heat transfer. Each spacer grid contains hundreds of small structures such as retaining clips, channel walls, and openings. The largest canisters store 37 PWR assemblies thus, there are numerous pins, tubes, and spacer grids for which the flow region between and around these structures needs to be determined along with the movement and deposition of aerosol particles. Because of the complicated geometry, modeling the intricate flow even for just one assembly is currently impractical. Nonetheless, we are developing techniques for a practical model to assess the natural aerosol particle deposition process in a canister in the event that a release occurs from one or more fuel pins. In the previous work it was demonstrated that CFD can model the flow through a PWR spacer grid with mixing vanes, including particle deposition, in a reasonable amount of time on a personal computer. In this work, the analysis is extended to include the bypass region between an assembly and the canister basket walls. It is shown that the flow velocity in the bypass region is about three times that of the interstitial region between the pins. The lengths before and after the spacer grid are also extended to determine when the flow becomes fully developed. In addition, the approach of computationally “stitching together” segments of an assembly is demonstrated with the plan to ultimately model a full assembly. The fraction of particles that are deposited in a segment with a spacer grid is determined as a function of particle size and flow velocity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Advection algorithms for quantum neutrino moment transport

Neutrino transport in compact objects is an inherently challenging multidimensional problem. Here, this difficulty is compounded if one includes flavor transformation—an intrinsically quantum phenomenon requiring one to follow the coherence between flavors and thus necessitating the introduction of complex numbers. To reduce the computational burden, simulations of compact objects that include neutrino transport often make use of momentum-angle-integrated moments (the lowest order ones being commonly referred to as the energy density and flux) and these quantities can be generalized to include neutrino flavor, i.e., they become quantum moments. Numerous finite-volume approaches to solving the moment evolution equations for classical neutrino transport have been developed based on solving a Riemann problem at cell interfaces. In this paper we describe our generalization of a Riemann solver for quantum moments, specifically decomposing complex numbers in terms of a (signed) magnitude and phase instead of real and imaginary parts. We then test our new algorithm in numerous cases showing a neutrino fast flavor instability, varying from toy models with analytic solutions to snapshots from neutron star merger simulations. Compared to previous algorithms for neutrino transport with flavor mixing, we find uniformly smaller growth rates of the flavor transformation along with concomitantly larger length-scales, and that the results are a better match with the growth rates seen from multiangle codes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Transport and mechanical behavior in PEO-LLZO composite electrolytes

Abstract Composite solid electrolytes (CEs), wherein ion-conducting polymer and ceramic/glass is mixed, are promising candidates for all-solid-state batteries due to their promise of acceptable ionic conductivity and mechanical properties compared to their individual constituents. While numerous studies have focused on improving the performance of CEs, it is still unclear what the material targets are that can result in improved macroscopic performance especially in light of the coupled needs for high transport and high mechanical strength in these materials. In this study, a two-dimensional (2D) mathematical model is developed to investigate electrochemical and mechanical characteristics of CEs. The model is compared to CEs consisting of poly-ethylene-oxide (PEO) polymer and lithium lanthanum zirconium oxide (LLZO) ceramic material with examination of the impact of varying LLZO volume fractions. The potential drop at the PEO-LLZO interface is evaluated using the junction potential theory. Using experimental data from the literature, the model estimates the ionic conductivity, effective transference number, and mechanical stiffness of the CEs. While the mechanical stiffness improves with increasing volume fraction of LLZO, the impact on conductivity and transference number depends on interfacial resistance at the interface. Finally, the study reports CE’s potential to enhance Li-ion transport and mechanical properties to inhibit lithium (Li) dendrite growth.

25 ENERGY STORAGE↗