Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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

A dynamic 2D Borehole Thermal Energy Storage (BTES) model for enhanced computational efficiency

Progressing toward a future increasingly reliant on renewable energy sources, the development of effective, durable energy storage solutions becomes essential to balance supply and demand fluctuations. Borehole Thermal Energy Storage (BTES) is a long-duration thermal energy storage technology that captures excess heat generated from renewable energy sources and stores it underground for later use, enabling the efficient utilization of sustainable energy. This approach is particularly valuable in district energy networks when integrated with Ground Source Heat Pumps (GSHP) to provide stable heating and cooling. However, traditional three-dimensional (3D) numerical models of BTES systems demand extensive computational resources, limiting their practicality for real-time and large-scale applications. This study introduces a novel two-dimensional (2D) modeling approach that reduces computational costs while maintaining high accuracy. By employing a radial ring-based discretization method, the model simulates heat injection, retention, and retrieval dynamics over seasonal cycles. A new thermal-mass weighted-average temperature parameter is introduced to evaluate the performance of BTES systems. Model validation against FEFLOW simulations demonstrates a 17-fold improvement in computational speed compared to traditional Computational Fluid Dynamics (CFD) models while achieving a mean absolute percentage error (MAPE) of 2 % during charging and 4 % during discharging. Additionally, a trade-off analysis between computational efficiency and accuracy is conducted, ensuring the model's applicability for real-world scenarios. The findings of this research contribute to the development of computationally efficient BTES models, facilitating better optimization, control, and integration into renewable energy systems. This work provides a foundation for further studies in techno-economic analysis, multi-year performance evaluation, and real-time operational strategies for BTES applications, supporting a more sustainable energy future.

2D modeling↗
DOE OSTI2022

AGR-5/6/7 Final Release-to-Birth Ratio Data Analysis

AGR-5/6/7 is the last of a series of Advanced Gas Reactor (AGR) experiments conducted in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) in support of development and qualification of tristructural isotropic (TRISO) low-enriched fuel for use in the high-temperature gas cooled reactor (HTGR). AGR configuration and irradiation conditions are based on prismatic HTGR technology that is distinguished primarily through the use of helium coolant, a low-power-density ceramic core capable of withstanding very high temperatures, and TRISO-coated particle fuel. The AGR tests provide valuable irradiation-performance data to support fuel process development, qualify fuel for normal operation and accident conditions, and support development and validation of fuel performance and fission-product (FP) transport models and codes. Each AGR test consists of multiple independently controlled and monitored capsules containing fuel compacts placed in a graphite cylinder shrouded by a steel shell. Release-to-birth ratios (R/B) for fission-gas isotopes released from each capsule are calculated from release rates, measured by germanium detectors in the Fission Product Monitoring System (FPMS) installed downstream from each capsule, and birth rates calculated using numerical models of FG generation. The R/Bs are a critical measure of the ability of the fuel kernel, the particle coating layers, and the compact matrix to retain fission-gas atoms, preventing their release into the sweep-gas flow, and the impact of initially defective particles and/or particle-coating failures that occur during irradiation. For fission-gas isotopes, particle failure is defined as failure of all coating layers, allowing gaseous fission atoms to escape from a particle. During the first five cycles (162B ? 165A), R/Bs were stable in the 10-8?10-6 range, and no in-pile particle failures were observed, based on the gross gamma counts. The maximum R/B value of around 2 ? 10-6 for Kr-85m resulted from the presence of as-fabricated exposed kernels (based on the high exposed kernel fraction), the dispersed uranium, and high fuel particle temperatures in Capsule 1. Comparison of capsule-measured R/Bs from these early cycles to predictions using the previously developed AGR R/B model demonstrated FG release from the AGR-5/6/7 TRISO fuel was comparable to that of previous experiments. In addition, the Kr-85m R/B per-exposed-kernel values are comparable to R/B values obtained in AGR-3/4 irradiation experiment and four irradiation experiments performed during 1980s: (1) HRB-17/18, (2) COMEDIE-BD1, (3) HFR B1, and (4) HRB-21. In contrast, all measured R/B values are lower than predictions by the commonly used Richards and German models, which are intentionally conservative. A large number of in-pile particle failures occurred in Capsule 1 by the end of Cycle 166A. During the final four cycles (166A ? 168A), apparent damage to the Capsule 1 gas line appeared to cause FG leakage from that capsule into the other four capsules, resulting in an increase in fission gas (FG) detected in the effluent for all capsules. Isolation of the Capsule 1 gas line during the last three cycles also prevented measurement of its FG release. Thus, R/Bs in all capsules after Cycle 166A are highly uncertain because of undefined amount of leakage from Capsule 1, especially for long-lived isotopes. A few hundred in-pile particle failures were estimated for Capsule 1 before the end of Cycle 166A, but the total number of failures is unknown due to the lack of FG release data in the later cycles. Based primarily on evidence from the gross gamma counts during Cycle 168A, approximately 15 particles failed in Capsule 3 and four particles failed in Capsule 2. In-pile failures in Capsule 3 were anticipated because this capsule was designed to operate beyond the HTGR normal operating temperature range. In contrast, no in-pile failures were identified in the top two capsules (4 and 5) based on the absence of the typical spikes in gross gamma counts and low failure estimates using the AGR model, developed in INL/EXT-14-32970, for R/B of the short-lived isotopes (Kr-89 and Xe-137) with minimal leakage from Capsule 1.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Community-Informed Urban Flood Modeling for Impact Mitigation

The intensification of the hydrologic cycle due to climate change poses a threat to aging and under-designed water infrastructure systems which cannot adequately manage intense storm events. Developing a comprehensive plan for managing rain-driven flooding events is challenging due to uncertainties in the magnitude and frequency of future storm events and conflicting stakeholder objectives. In the City of Baltimore, Maryland, stormwater infrastructure is struggling to keep up with rainfall-driven (pluvial) flooding events, which regularly damage housing and disrupt transportation for residents. In this study, a hybrid of community engagement, numerical modeling, and artificial intelligence techniques are employed to explore prospective urban flooding adaptations. Community engagement drives the development of an urban flooding model (EPA Storm Water Management Model) for the Baltimore Harbor watershed. The model integrates complex surface and subsurface stormwater infrastructure data from the City, high-resolution spatial data, insights from local public works experts, and the lived experiences of City residents. This co-developed model simulates adaptations of interest to stakeholders in the city, including green and grey infrastructure and operational management strategies. Stormwater management scenarios focused on inlet cleaning and spatially concentrated green infrastructure are found to be the most effective in reducing flood depths in community priority locations. Together, these adaptations can reduce the duration of intersection inundation by more than twenty minutes, allowing for quicker emergency response and restoration of typical transportation systems. Future work will combine this community engaged flooding model with the Deep Uncertainties Pathways framework to explore tradeoffs between adaptations and develop dynamic adaptations which align with community objectives, enhance climate resilience in Baltimore, and can be adjusted in response to changing future conditions.

Ava, Spangler [Pennsylvania State University]↗
DOE OSTI2022

Measuring binary fluidization of nonspherical and spherical particles using machine learning aided image processing

Abstract The binary fluidization of Geldart D type nonspherical wood particles and spherical low density polyethylene (LDPE) particles was investigated in a laboratory‐scale bed. The experiment was performed for varying static bed height, wood particles count, as well as superficial gas velocity. The LDPE velocity field were quantified using particle image velocimetry (PIV). The wood particles orientation and velocity are measured using particle tracking velocimetry (PTV). A machine learning pixel‐wise classification model was trained and applied to acquire wood and LDPE particle masks for PIV and PTV processing, respectively. The results show significant differences in the fluidization behavior between LDPE only case and binary fluidization case. The effects of wood particles on the slugging frequency, mean, and variation of bed height, and characteristics of the particle velocities/orientations were quantified and compared. This comprehensive experimental dataset serves as a benchmark for validating numerical models.

Li, Cheng↗
DOE OSTI2021

Experimental and modeling investigation of binary liquid mixtures

Binary mixtures of liquids may be encountered in industrial or remote sensing scenarios and present challenges to positive identification compared to neat single-component liquids. Our investigation examines whether one can predict the optical properties of the mixture, i.e. its complex index of refraction, by assuming a linear superposition of the real and imaginary components of the index of refraction in proportion to the ratio of each constituent. To investigate this hypothesis various liquid mixtures were created using mass ratios. The mixtures were then characterized as to their complex index of refraction and used in numerical modeling calculations of thin liquid mixture films on surfaces and compared with composed mixtures using linear n and k synthetic mixtures where the n and k components of the complex index of refraction were combined in similar ratios. The comparison of modeling and experimental results is presented with recommendations for further investigation.

infrared (IR) spectroscopy, Liquid mixtures, compl↗
DOE OSTI2025

Decoupling Carrier Dynamics and Energy Transport in Ultrafast Near-Field Nanoscopy

Ultrafast near-field optical nanoscopy has emerged as a powerful platform to characterize low-dimensional materials. While analytical and numerical models have been established to account for photoexcited carrier dynamics, quantitative evaluation of the associated pulsed laser heating remains elusive. Here, we decouple the photocarrier density and temperature increase in near-field nanoscopy by integrating the two-temperature model (TTM) with finite-difference time-domain (FDTD) simulations. These results reveal that the electron–phonon coupling in a silicon film after femtosecond laser excitation is most pronounced within approximately 3 ps–substantially shorter than the photocarrier decay time scale at tens of picoseconds. Moreover, the coupled TTM-FDTD method indicates that ultrafast laser heating can cause up to a 14% variation in the near-field signal at a 220 μJ/cm 2 pump pulse fluence. Our numerical results are further validated by transient experiments, highlighting the potential of this method for investigations of carrier and thermal phenomena in emerging nanomaterials and nanodevices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Physical science research needed to evaluate the viability and risks of marine cloud brightening

Marine cloud brightening (MCB) is the deliberate injection of aerosol particles into shallow marine clouds to increase their reflection of solar radiation and reduce the amount of energy absorbed by the climate system. From the physical science perspective, the consensus of a broad international group of scientists is that the viability of MCB will ultimately depend on whether observations and models can robustly assess the scale-up of local-to-global brightening in today’s climate and identify strategies that will ensure an equitable geographical distribution of the benefits and risks associated with projected regional changes in temperature and precipitation. To address the physical science knowledge gaps required to assess the societal implications of MCB, we propose a substantial and targeted program of research—field and laboratory experiments, monitoring, and numerical modeling across a range of scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Exploration of lattice Hamiltonians for functional and structural discovery via Gaussian process-based exploration–exploitation

Statistical physics models ranging from simple lattice to complex quantum Hamiltonians are one of the mainstays of modern physics that have allowed both decades of scientific discovery and provided a universal framework to understand a broad range of phenomena from alloying to frustrated and phase separated materials to quantum systems. Traditionally, exploration of the phase diagrams corresponding to multidimensional parameter spaces of Hamiltonians was performed using a combination of basic physical principles, analytical approximations, and extensive numerical modeling. However, exploration of complex multidimensional parameter spaces is subject to the classic dimensionality problem, and the behaviors of interest concentrated on low dimensional manifolds remain undiscovered. Here, we demonstrate that a combination of exploration and exploration–exploitation with Gaussian process modeling and Bayesian optimization allows effective exploration of the parameter space for lattice Hamiltonians and effectively maps the regions at which specific macroscopic functionalities or local structures are maximized. We argue that this approach is general and can be further extended well beyond the lattice Hamiltonians to effectively explore the parameter space of more complex off-lattice and dynamic models.

42 ENGINEERING↗
DOE OSTI2026

Advancing the HERO WEC Through Integrated Modeling, Testing, and Field Deployments: Preprint

The Hydraulic and Electric Reverse Osmosis Wave Energy Converter (HERO WEC) was developed by the National Laboratory of the Rockies as a modular platform for advancing wave-powered desalination technologies. Designed to operate in either a hydraulic or electric configuration, the system enables direct comparison of mechanical-to-water and electrical-to-water conversion pathways using a single hardware architecture. Deployments at the Jennette's Pier test site in 2022 and 2024 demonstrated freshwater production between 60 and 300 gallons per day, including successful operation in wave heights as low as 0.25 m. A structured evaluation approach combining numerical modeling of hydrodynamic and PTO response, controlled laboratory testing, and in-ocean field deployments has been used to characterize and refine system performance. Insights from these efforts are now informing the second-generation HERO WEC (V2), which incorporates improved drivetrain components, refined geometry, enhanced control systems, and design updates aimed at increasing robustness and long-duration survivability. The HERO WEC platform continues to serve as an open-access research asset supporting broader marine energy and desalination development.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Using AI to build a hydrobiogeochemical soil model

Soil water content is a function of inputs from precipitation and outputs via evaporation, transpiration, lateral flow, and vertical percolation, and is sensitive to biogeochemical processes. As such, soils serve as an ideal integrator of atmospheric, hydrological, and biogeochemical processes affecting the water cycle. In addition, soil water retention capacity, infiltration rates, and hydraulic conductivity can buffer or exacerbate the effects of extreme precipitation events (e.g., flooding, runoff, subsurface transport, erosion, greenhouse gas emissions) and mitigate the impact of droughts and heat waves on land systems (e.g., fire, crop failure). However, integrating water cycle measurements spanning different land atmosphere compartments across scales is a fundamental barrier for numerical model predictability. A significant challenge is that each domain (soil, hydrology, biology, and atmosphere) typically collects different sets of data at different temporal and spatial frequencies/scales, and even different dimensionalities (2D vs 3D). To implement soil as an integrator of the water cycle in land models, we suggest that novel machine learning (ML) tools can be developed to effectively simulate complex landscapes across various domains and scales, extended to regions with sparse or no data. The ultimate goals are to improve predictive understanding of land-atmosphere interactions and to extend the predictability of current Earth System Models (ESMs) through better integration of hydrological and biogeochemical data. We envision a framework in which: (1) ML-aided data reconstructions enable the merger of data sources into a unified geospatial product; (2) automated detection techniques are used to improve the knowledge of complex soil processes and interactions; and (3) this knowledge is leveraged and incorporated into models through AI-based emulators to distinctly connect the land and atmospheric compartments of the water cycle in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Advancing the HERO WEC Through Integrated Modeling, Testing, and Field Deployments

The Hydraulic and Electric Reverse Osmosis Wave Energy Converter (HERO WEC) was developed by the National Laboratory of the Rockies as a modular platform for advancing wave-powered desalination technologies. Designed to operate in either a hydraulic or electric configuration, the system enables direct comparison of mechanical-to-water and electrical-to-water conversion pathways using a single hardware architecture. Deployments in 2022 and 2024 at the Jennette's Pier test site in North Carolina demonstrated freshwater production between 60 and 300 gallons per day, including successful operation in wave heights as low as 0.25 m. A structured evaluation approach combining numerical modeling of hydrodynamic and power take-off response, controlled laboratory testing, and in-ocean field deployments has been used to characterize and refine system performance. Specifically, insights from these efforts are now informing the second-generation HERO WEC (V2) to increase robustness and overall survivability. Changes from the original prototype include improved drivetrain components, refined geometry, enhanced control systems, and design. The HERO WEC platform continues to serve as an open-access research asset supporting broader marine energy and desalination development.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Beam breakup instability studies of powerful energy recovery linac for experiments

The maximum achievable beam current in an energy recovery linac (ERL) is often constrained by beam breakup (BBU) instability. Our previous research highlighted that filling patterns have a substantial impact on BBU instabilities in multipass ERLs. In this study, we extend our investigation to the eight-cavity model of the Powerful ERL for Experiment (PERLE). We evaluate its requirements for damping cavity higher order modes (HOMs) and propose optimal filling patterns and bunch timing strategies. Our findings reveal a significant new insight: while filling patterns are crucial, the timing of bunches also plays a critical role in mitigating HOM beam loading and BBU instability. This previously underestimated factor is essential for effective BBU control. We estimated the PERLE threshold current using both analytical and numerical models, incorporating the designed PERLE HOM dampers. During manufacturing, HOM frequencies are expected to vary slightly. Our study found no significant difference in BBU suppression for relative rms frequency jitters of 0.001, 0.002, and 0.005 for the same HOM. Introducing a jitter of 0.001 into our models, we found that the dampers effectively suppressed BBU instability, achieving a threshold current an order of magnitude higher than the design requirement. Our results offer new insights into ERL BBU beam dynamics and have important implications for the design of future ERLs. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Creep anisotropy modeling and uncertainty quantification of an additively manufactured Ni-based superalloy

The advantages offered by additive manufacturing over traditional processes has driven a great deal of industrial and academic interest in recent years. However, the process is relatively new and requires additional investigation to become sufficiently mature for wide scale industrial adoption. Electron beam melting powder bed fusion is one technology that has shown promise for fabricating high temperature resistant materials such as nickel based superalloys. The resulting microstructures typically exhibit a strong fiber texture in the build direction giving rise to anisotropic time-dependent deformation behavior. In order to accelerate the qualification of these materials for industrial adoption accurate numerical models are needed for simulating their behavior. In this work a crystal plasticity model including non-Schmid effects is presented for capturing creep anisotropy observed in additively manufactured IN738LC. The model is calibrated via a probabilistic framework where model parameters are treated as random variables. An iterative sequential design strategy is utilized to efficiently identify the probability density of the unknown model parameters. As a case study the model is utilized to investigate the behavior of randomly oriented equiaxed grain clusters sometimes observed embedded in the additively manufactured columnar structure. A synthetic realization is simulated and uncertainty is propagated through to the full-field response. Results indicate that these features are the source of significant creep relaxation and strain accumulation which partially explains observed grain boundary decohesion at these locations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Improved accuracy and transferability of molecular-orbital-based machine learning: Organics, transition-metal complexes, non-covalent interactions, and transition states

We report molecular-orbital-based machine learning (MOB-ML) provides a general framework for the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. The application of Nesbet’s theorem makes it possible to recast a typical extrapolation task, training on correlation energies for small molecules and predicting correlation energies for large molecules, into an interpolation task based on the properties of orbital pairs. We demonstrate the importance of preserving physical constraints, including invariance conditions and size consistency, when generating the input for the machine learning model. Numerical improvements are demonstrated for different datasets covering total and relative energies for thermally accessible organic and transition-metal containing molecules, non-covalent interactions, and transition-state energies. MOB-ML requires training data from only 1% of the QM7b-T dataset (i.e., only 70 organic molecules with seven and fewer heavy atoms) to predict the total energy of the remaining 99% of this dataset with sub-kcal/mol accuracy. This MOB-ML model is significantly more accurate than other methods when transferred to a dataset comprising of 13 heavy atom molecules, exhibiting no loss of accuracy on a size intensive (i.e., per-electron) basis. It is shown that MOB-ML also works well for extrapolating to transition-state structures, predicting the barrier region for malonaldehyde intramolecular proton-transfer to within 0.35 kcal/mol when only trained on reactant/product-like structures. Finally, the use of the Gaussian process variance enables an active learning strategy for extending the MOB-ML model to new regions of chemical space with minimal effort. We demonstrate this active learning strategy by extending a QM7b-T model to describe non-covalent interactions in the protein backbone–backbone interaction dataset to an accuracy of 0.28 kcal/mol.

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

Deterministic and Statistical Modeling of a New Thermal Breakout Technology for Measuring the Maximum Horizontal In-Situ Stress

The current state-of-the-art technology for in-situ stress measurements involves an integrated approach that combines borehole breakout observations, drilling-induced tensile fractures, and hydraulic fracturing tests (i.e., “mini-fracs”). This methodology has achieved wide application in the oil and gas industry but has several limitations that often prevent successful in-situ stress measurements. One major limitation is that breakouts do not appear in all boreholes and are generally only a natural occurrence that cannot easily be controlled. Because borehole breakouts are used to directly measure the maximum horizontal in-situ stress magnitude, the absence of borehole breakouts presents a major data gap for in-situ stress measurements. In response to this data gap, a new US Department of Energy (US DOE)-sponsored thermal breakout technology that will provide a method for thermally inducing borehole breakouts and allow the consistent measurement of the maximum horizontal stress magnitude is currently in development. This thermal breakout technology involves heating the borehole and increasing the thermoelastic compressive stress in the rock until a breakout develops, which can be directly correlated to the maximum horizontal stress magnitude. The first step in this project was an analytical modeling study of the thermal breakout process. Based on the Kirsch solution (Kirsch 1898), a deterministic and statistical analysis was performed on the pertinent parameters that influence the maximum horizontal stress calculation. As a result of the analysis, the findings indicate that the thermal breakout technology is feasible and provides improved accuracy and/or an enhanced ability to measure the maximum horizontal stress. Future work as part of this US DOE-sponsored project includes additional validation through more detailed numerical modeling, laboratory testing, and field testing of the thermal breakout technology.

58 GEOSCIENCES↗
DOE OSTI2021

Assessment of Turbulent Combustion Models for Simulating Prechamber Ignition in a Natural Gas Engine

Lean combustion is a promising strategy to increase thermal efficiency in an internal combustion engine, by exploiting a favorable specific heat ratio of the fresh mixture while simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition and slow flame propagation at fuel-lean conditions lead to large cycle-to-cycle variability and limit the high-efficiency engine operating range. Prechamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling a faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling standpoint to date science base and available simulation tools are inadequate to properly understand and predict the combustion processes in prechamber ignited engines. In this paper, conceptually different Reynolds-averaged Navier–Stokes (RANS) combustion models widely adopted in the engine modeling community are used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a prechamber spark-ignition system. Here, a flamelet-based turbulent combustion model, i.e., G-equation, and a multizone well-stirred reactor model are employed for this modeling study. Simulation results are compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two evaluated approaches.

42 ENGINEERING↗
DOE OSTI2024

A Phenomenological Thermal Spray Wall Interaction Modeling Framework Applied to a High Temperature Ignition Assistant Device

Airborne compression ignition engines must operate with reliable ignition systems to achieve proper ignition at every cycle, particularly at high altitudes. Glow-plug-based ignition-assistant (IA) devices can provide the necessary energy to preheat the fuel and ensure ignitability of the fuel-air mixture. Ignitability of liquid sprays can be facilitated via direct impingement onto the hot IA surface, however this comes with adverse effects on the IA durability. Therefore, optimizing an IA’s design requires detailed understanding of the physics of fuel spray impingement of superheated surfaces. While spray impingement on relatively low wall temperatures has been extensively studied and appropriate numerical models have been proposed through the years, fundamental understanding of high-speed liquid spray impingement on superheated walls is still elusive. This work aims to formulate a phenomenological thermal spray-wall interaction framework for modeling the film-boiling-induced heat transfer, atomization, and dispersion of fuel spray droplets impinging onto a superheated IA device. A qualitative comparison of the new phenomenological model is performed against optical experiments from the literature of an F-24 fuel spray injected onto an IA device located 12 mm away from the injector tip. The temperature of the IA was set at 1400 K. The fuel injection pressure was 400 bar, while the ambient gas pressure and temperature were 30 bar and 800 K, respectively. The performance of the phenomenological model is evaluated in comparison with two other state-of-art models from the literature. A qualitative analysis of the different spray and fuel-air mixture characteristics is performed to outline the differences in the predictions offered by the new phenomenological model and the two state-of-art spray-wall interaction models.

droplet dispersion↗
DOE OSTI2020

Polysulfide entrapment and retardation in gel electrolyte Li–S batteries: experiments and modeling

Gelling the electrolyte in lithium–sulfur (Li–S) batteries provides advantages in terms of safety and mitigating polysulfide shuttling leading to reduced side reactions and dendrite growth on the anode side. In the current study, organic and ceramic crosslinkers have been used to improve the entrapment of polysulfides and prevent them from shuttling to the anode side. However, the performance of gel electrolytes is still far from being on par with liquid electrolytes in terms of raw capacity. The discharge curves exhibit an extended first plateau for gel electrolyte batteries compared to liquid electrolytes, showing that polysulfides are better trapped in the gel system. Yet, in the conversion of soluble polysulfides into insoluble polysulfides, the second plateau was delayed. Qualitative and quantitative analysis of the polysulfide diffusion in the liquid electrolyte and gel electrolyte systems were also conducted to predict the relative diffusivity ratio between the liquid electrolyte and gel electrolyte for simulation. To probe what happens at the end of the discharge, where the mass transport in the gel electrolyte appears to be slow, a modified discharge procedure was used where the cell is first paused after discharge and then discharged a second time before recharging normally. Coin cell tests for two gel electrolytes in addition to liquid electrolyte show up to 7% capacity recovery after the pause, which indicates that the diffusion of species during the pause is responsible for the failure recovery. Additionally, the recovery is higher for gel electrolytes compared to liquid, for higher C-rates, and for longer pause times. To better understand this behavior, a Li–S numerical model with mass transport limited reactions was used to examine the polysulfide diffusion coefficients expected in different electrolyte systems. The model is able to reproduce the trends seen in experiments and yields a higher recovery for smaller Li 2 S 4 diffusion coefficients, suggesting that insufficient Li 2 S 4 mass transport is responsible for failure at the end of discharge. By identifying the effect that gel electrolytes have on the state of the battery during discharge using both experimental and simulation methods, this work is a step towards leveraging the advantages of gel electrolytes and mitigating their weaknesses.

25 ENERGY STORAGE↗