Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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 1,729 records · Page 96

Carbonate Management to Enable Energy- and Carbon-Efficient CO 2 Electrolysis (Final Technical Report)

The rapid growth and plummeting cost of solar energy have spurred growing interest in using CO 2 electrolysis to produce chemicals and fuels as an alternative to conventional petrochemical processes. High-temperature (>800 °C) solid oxide electrolyzers that convert CO 2 into CO and O 2 have recently become commercially available. Low-temperature electrolysis cells offer the prospect of more convenient and flexible operation, which is critical for utilizing intermittent solar energy, and provide access to more valuable C 2+ products such as ethylene, ethanol, and propanol. Over the past 10 years, research in this area has yielded substantial progress in both fundamental understanding of the requisite electrocatalytic reactions and design of prototype devices. Leveraging insights from fuel cells and membrane water electrolyzers, researchers have developed electrolysis cells with gas diffusion electrodes (GDE) that have demonstrated high CO 2 electrolysis current densities (>100 mA cm –2 ) as well as promising selectivity and stability. Despite these advances, the energy efficiency (electrical energy-to-product) and carbon efficiency (CO 2 -to-product) of low-temperature CO 2 electrolysis remain far too low for large-scale deployment. A preponderance of evidence indicates that the principal source of efficiency losses is the rapid and thermodynamically favorable reaction of CO 2 with hydroxide (OH – ) to form carbonate (CO 3 2– ). Carbonate formation imposes steady-state electrolysis conditions that result in large voltage and CO 2 losses for all known (photo)electrochemical CO 2 cells. While much current research remains focused on CO 2 reduction catalyst design, this largely overlooked CO 3 2– problem presents a fundamental scientific barrier to creating a viable electrochemical option for converting solar energy into chemicals and fuels. The project pursues an integrated, multi-PI research effort that establishes a fundamental science of CO 3 2– management. PI Kanan and Co-PI Mani’s contribution to the project is to evaluate strategies to mitigate the CO 3 2– problem by changing the properties of the electrolyte and the environment in which CO 2 reduction catalysis takes place. Experimental studies showed that electrolytes composed of a high concentration of both CO 3 2– and HCO 3 – , which serve as moderately alkaline buffers, improved the cell voltage by compared to all-HCO 3 – electrolytes, but these buffered systems still show substantial CO 2 uptake that reduces pH over time. Computational studies developed a homogenized model of a CO 2 reduction catalyst layer that permits a low-cost exploration of the high-dimensional parameter space associated with catalyst layers on gas diffusion electrodes. The model was validated by accurately reproducing experimental data for the related but simpler reaction of CO reduction and then used to probe the effects of catalyst layer architecture on CO 2 reduction. Minimizing the size of catalyst and hydrophobic domains in the catalyst layer is predicted to mitigate CO 3 2– formation and thereby enable prolonged operation at elevated pH. In support of the CO 2 electrolysis studies, a new method for rapidly prototyping electrochemical cells was developed and validated. The method uses a combination of 3D printing and electroless plating to generate conductive cell components for evaluating new cell designs. The carbonate problem encompasses mass transport processes and acid-base reactions that are relevant to many other electrochemical systems. Investigation of strategies to address the carbonate problem led to an additional line of inquiry into the physicochemical phenomena that determine the efficiency of electrochemical acid-base production, which has numerous applications in the broader field of carbon management. New strategies for using the supporting electrolyte to inhibit H + /OH – recombination in electrochemical acid-base production were evaluated, leading to the development of a novel acid-base producing system that eliminates the need for ion exchange membranes and exhibits promising efficiency and current densities for scalable applications.

25 ENERGY STORAGE↗

Current challenges in the physics of white dwarf stars

White dwarfs are a class of stars with unique physical properties. They present many challenging problems whose solution requires advanced theories of dense matter, state-of-the-art experimental techniques, and extensive computing efforts. New ground- and space-based observatories will soon provide an increasingly detailed view of white dwarf stars and reveal new phenomena that will challenge our models. This review is an introduction to the nature of white dwarfs, the physical processes that determine their structure and evolution, and the physical conditions they span. Here we discuss a wide variety of currently unsolved or partially resolved problems in their constitutive physics that are broadly related to equations of state, transport processes and opacities.

79 ASTRONOMY AND ASTROPHYSICS↗

Features in simulation of crystal growth using the hyperbolic PFC equation and the dependence of the numerical solution on the parameters of the computational grid

We investigate the three-dimensional mathematical model of crystal growth called PFC (Phase Field Crystal) in a hyperbolic modification. This model is also called the modified model PFC (originally PFC model is formulated in parabolic form) and allows to describe both slow and rapid crystallization processes on atomic length scales and on diffusive time scales. Modified PFC model is described by the differential equation in partial derivatives of the sixth order in space and second order in time. The solution of this equation is possible only by numerical methods. Previously, authors created the software package for the solution of the Phase Field Crystal problem, based on the method of isogeometric analysis (IGA) and PetIGA program library. During further investigation it was found that the quality of the solution can strongly depends on the discretization parameters of a numerical method. In this report, we show the features that should be taken into account during constructing the computational grid for the numerical simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Non-Equilibrium Doping of Poly-Si:Ga/SiOx Passivating Contacts Solar Cells

In this contribution, we show that Ga-hyperdoped poly-Si/SiOx contacts are realized by pulsed laser annealing and post-hydrogenation. The motivation behind this work is to improve the inferior performance of boron-doped poly- Si/SiOx contacts when compared to the phosphorus-doped contacts. Here we used Ga as a novel p-type dopant to replace B, as B is known to accumulate at the interface between the oxide and the crystalline silicon and cause recombination. Due to a high diffusivity of Ga in oxide and a much higher segregation coefficient (Si:SiOx), Ga is shown to be immune to the oxide pileup phenomenon. However, Ga has a low solid solubility limit in Si (5x10^19 cm-3), which is not enough to form a low resistivity contact. To overcome this problem, we explored pulsed laser annealing via excimer laser by varying the laser energy densities and pulses. The passivation and electrical properties are studied systematically using photoluminescence (PL), secondary ion mass spectrometry (SIMS), and Van der Pauw-Hall measurements. We show an active Ga doping concentration above 1x10^20 cm^-3, which is above the solid solubility limit of Ga in Si. Good passivation quality is demonstrated with the highest implied Voc (iVoc) of 721 mV on single side polished samples. Furthermore, a low diode contact resistivity of ~22 m..omega.. cm^2 is obtained using Ti/Ag contacts, with minimal metal-to-semiconductor contact resistivity of 0.9 m..omega.. cm^2.

non-equilibrium doping↗

CanTrilBat and Cantera_apps v1.0 beta

This application can determine the performance and chemical behavior of batteries in 1D when they are cycled. With CanTrilBat, we are developing predictive phenomenological models for battery systems to predict operating performance and rate limiting steps in the performance of battery models. Particular attention is paid to primary and secondary chemistry mechanisms, such as the thermal runaway mechanisms experienced in secondary lithium ion batteries or self-discharge reaction mechanism that all batteries experience to one extent or another. The first application of this model has been for modeling the performance of thermal batteries. However, an implementation for secondary ion batteries is next. CanTrilBat applications solves transient problems involving batteries. It is a 1-D application that represents 3-D physical systems that can be reduced using the porous flow approximation for the anode, cathode, and separator. A control volume formulation is used to track conserved quantities. An operator-split approach is used to calculate the chemistry, diffusion and electronic transport that occurs within cathode and anode particles, allowing for the reduction in code complexity. All jacobian operations in CanTrilBat utilize numerical jacobians. A home grown predictor corrector scheme is used for time step control, and a home grown newton solver is used to relax the equations at each time step. Trilinos is used to solve for the resulting linear systems equation A block text input format is used to initialize options for the CanTriBat program. Within this block are the names of several XML input files, which specify the chemistry mechanism. These XML input files are read by low level Cantera routines, and serve to initialize the electrode and electrolyte chemistry mechanisms and transport properties. A GUI implementation has been contracted out to a university professor, but has not been implemented yet. It?s expected that CanTriBat will have the capability to do adaptive grid refinement in later versions. Grids will be halved based on an complicated set of heuristic criteria, based on how well the mesh captures changes to the first and second derivatives of selected solution components. Primary emphasis is given towards converged profiles for IV performance.

Moffat, HarryK.↗

Orbital Modification of the Himalia Family during an Early Solar System Dynamical Instability

Among the irregular satellites orbiting Jupiter, the Himalia family is characterized by a high velocity dispersion δv of several hundred m s{sup −1} among its members, inconsistent with a collisional origin. Efforts to account for this through internecine gravitational interactions do not readily reproduce this feature. Here, we revisit the problem in the context of recent cosmogonical models, where the giant planets migrated significantly through interaction with a planetesimal disk and suffered encounters with planetesimals and planet-sized objects. Our starting assumption is that family formation either predated this phase or occurred soon after its onset. We simulate numerically the diffusive effect of three distinct populations of perturbers on a set of test particles representing the family: Moon-sized (MPT) and Pluto-sized (PPT) planetesimals, and planetary-mass objects (PMO) with masses typical of ice-giant planets. We find that PPT flybys are inefficient, but encounters with MPTs raise the δv of ∼60% of our test particles to >200 m s{sup −1} with respect to Himalia, in agreement with observations. As MPTs may not have been abundant in the disk, we simulate encounters between Jupiter and PMOs. We find that too few encounters generate less dispersion than MPTs while too many essentially destroy the family. For PMO masses in the range 5--20 m{sub ⊕}, the family orbital distribution is reproduced by a few tens of encounters.

79 ASTRONOMY AND ASTROPHYSICS↗

Challenges of accelerated aging techniques for elastomer lifetime predictions

Elastomers are often degraded when exposed to air or high humidity for extended times (years to decades). Lifetime estimates normally involve extrapolating accelerated aging results made at higher than ambient environments. Several potential problems associated with such studies are reviewed, and experimental and theoretical methods to address them are provided. The importance of verifying time–temperature superposition of degradation data is emphasized as evidence that the overall nature of the degradation process remains unchanged versus acceleration temperature. The confounding effects that occur when diffusion-limited oxidation (DLO) contributes under accelerated conditions are described, and it is shown that the DLO magnitude can be modeled by measurements or estimates of the oxygen permeability coefficient (P Ox ) and oxygen consumption rate (Φ). P Ox and Φ measurements can be influenced by DLO, and it is demonstrated how confident values can be derived. In addition, several experimental profiling techniques that screen for DLO effects are discussed. Values of Φ taken from high temperature to temperatures approaching ambient can be used to more confidently extrapolate accelerated aging results for air-aged materials, and many studies now show that Arrhenius extrapolations bend to lower activation energies as aging temperatures are lowered. Furthermore, best approaches for accelerated aging extrapolations of humidity-exposed materials are also offered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical investigation of internal high-speed viscous flows using a parabolic technique

A feasibility study has been conducted to assess the applicability of an existing parabolic analysis (ADD-Axisymmetric Diffuser Duct), developed previously for subsonic viscous internal flows, to mixed supersonic/subsonic flows with heat addition simulating a SCRAMJET combustor. A study was conducted with the ADD code modified to include additional convection effects in the normal momentum equation when supersonic expansion and compression waves are present. A set of test problems with weak shock and expansion waves have been analyzed with this modified ADD method and stable and accurate solutions were demonstrated provided the streamwise step size was maintained at levels larger than the boundary layer displacement thickness. Calculations made with further reductions in step size encountered departure solutions consistent with strong interaction theory. Calculations were also performed for a flow field with a flame front in which a specific heat release was imposed to simulate a SCRAMJET combustor. In this case the flame front generated relatively thick shear layers which aggravated the departure solution problem. Qualitatively correct results were obtained for these cases using a marching technique with the convective terms in the normal momentum equation suppressed. It is concluded from the present study that for the class of problems where strong viscous/inviscid interactions are present a global iteration procedure is required.

Anderson, O. L.↗

Flux-difference split parabolized Navier-Stokes algorithm for non-equilibrium chemically reacting flows

A flux-difference split explicit finite-difference algorithm is presented for solving the parabolized form of the equations governing three-dimensional nonequilibrium chemically reacting flows. The algorithm is based on an explicit noniterative, upwind space-marching scheme developed by Korte, but differs in that the unsteady Riemann problem, rather than the steady Riemann problem, is solved. The algorithm allows either a second or an approximately third-order accurate upwind treatment of the convection terms by employing the unsteady approximate Riemann solver of Roe. The source terms of the species transport equations are treated in either an explicit or implicit manner, and the species diffusion terms are modeled with either a Fickian or a multicomponent model. A validation of the algorithm is performed by comparing computational results with the 2-D Mach 14, 15 degree compression-corner data of Holden. The three-dimensional capability of the algorithm is demonstrated by computing Mach 2.7 flow over a swept wedge scramjet fuel injector, and three-dimensional reacting flow capability is demonstrated by a computing a shock-jet interaction concept for mixing and combustion enhancement.

White, J. A.↗

Dynamics and structure of turbulent premixed flames

In earlier work (Mantel & Bilger, 1994) the structure of the turbulent premixed flame was investigated using statistics based on conditional averaging with the reaction progress variable as the conditioning variable. The DNS data base of Trouve and Poinsot (1994) was used in this investigation. Attention was focused on the conditional dissipation and conditional axial velocity in the flame with a view to modeling these quantities for use in the conditional moment closure (CMC) approach to analysis of kinetics in premixed flames (Bilger, 1993). Two remarkable findings were made: there was almost no acceleration of the axial velocity in the flame front itself; and the conditional scalar dissipation remained as high, or higher, than that found in laminar premixed flames. The first finding was surprising since in laminar flames all the fluid acceleration occurs through the flame front, and this could be expected also for turbulent premixed flames at the flamelet limit. The finding gave hope of inventing a new approach to the dynamics of turbulent premixed flames through use of rapid distortion theory or an unsteady Bernoulli equation. This could lead to a new second order closure for turbulent premixed flames. The second finding was contrary to our measurements with laser diagnostics in lean hydrocarbon flames where it is found that conditional scalar dissipation drops dramatically below that for laminar flamelets when the turbulence intensity becomes high. Such behavior was not explainable with a one-step kinetic model, even at non-unity Lewis number. It could be due to depletion of H2 from the reaction zone by preferential diffusion. The capacity of the flame to generate radicals is critically dependent on the levels of H2 present (Bilger, et al., 1991). It seemed that a DNS computation with a multistep reduced mechanism would be worthwhile if a way could be found to make this feasible. Truly innovative approaches to complex problems often come only when there is the opportunity to work close at hand with the (in this case numerical) experimental data. Not only can one spot patterns and relationships in the data which could be important, but one can also get to know the limitations of the technique being used, so that when the next experiment is being designed it will address resolvable questions. A three-year grant from the Australian Research Council has enabled us to develop a small capability at the University of Sydney to work on DNS of turbulent reacting flow, and to analyze data bases generated at CTR. Collaboration between the University of Sydney and CTR is essential to this project and finding a workable modus operandum for this collaboration, given the constraints involved, has been a major objective of the past year's effort. The overall objectives of the project are: (1) to obtain a quantitative understanding of the dynamics of turbulent premixed flames at high turbulence levels with a view to developing improved second order closure models; and (2) to carry out new DNS experiments on turbulent premixed flames using a carefully chosen multistep reduced mechanism for the chemical kinetics, with a view to elucidating the laser diagnostic findings that are contrary to the findings for DNS using one-step kinetics. In this first year the objectives have been to make the existing CTR data base more accessible to coworkers at the University of Sydney, to make progress on understanding the dynamics of the flame in this existing CTR data base, and to carefully construct a suitable multistep reduced mechanism for use in a new set of DNS experiments on turbulent premixed flames.

Bilger, R. W.↗

Lunar and Planetary Science XXXV: Meteorites

The session "Meteorites" included the following reports:Description of a New Stony Meteorite Find from Bulloch County, Georgia; Meteorite Ablation Derived from Cosmic Ray Track Data Dhofar 732: A Mg-rich Orthopyroxenitic Achondrite Halogens, Carbon and Sulfur in the Tagish Lake Meteorite: Implications for Classification and Terrestrial Alteration; Electromagnetic Scrape of Meteorites and Probably Columbia Tiles; Pre-Atmospheric Sizes and Orbits of Several Chondrites; Research of Shock-Thermal History of the Enstatite Chondrites by Track, Thermoluminescence and Neutron-Activation (NAA) Methods; Radiation and Shock-thermal History of the Kaidun CR2 Chondrite Glass Inclusions; On the Problem of Search for Super-Heavy Element Traces in the Meteorites: Probability of Their Discovery by Three-Prong Tracks due to Nuclear Spontaneous Fission Trace Element Abundances in Separated Phases of Pesyanoe, Enstatite Achondrite; Evaluation of Cooling Rate Calculated by Diffusional Modification of Chemical Zoning: Different Initial Profiles for Diffusion Calculation; Mineralogical Features and REE Distribution in Ortho- and Clinopyroxenes of the HaH 317 Enstatite Chondrite Dhofar 311, 730 and 731: New Lunar Meteorites from Oman; The Deuterium Content of Individual Murchison Amino Acids; Clues to the Formation of PV1, an Enigmatic Carbon-rich Chondritic Clast from the Plainview H-Chondrite Regolith Breccia ;Numerical Simulations of the Production of Extinct Radionuclides and ProtoCAIs by Magnetic Flaring.

Source record↗

Does the Hall Effect Solve the Flux Pileup Saturation Problem?

It is well known that magnetic flux pileup can significantly speed up the rate of magnetic reconnection in high Lundquist number resistive MHD,allowing reconnection to proceed at a rate which is insensitive to the plasma resistivity over a wide range of Lundquist number. Hence, pileup is a possible solution to the Sweet-Parker time scale problem. Unfortunately, pileup tends to saturate above a critical value of the Lundquist number, S_c, where the value ofS_c depends on initial and boundary conditions, with Sweet-Parker scaling returning above S_c. It has been argued (see Dorelli and Bim [2003] and Dorelli [2003]) that the Hall effect can allow flux pileup to saturate (when the scale of the current sheet approaches ion inertial scale, di) before the reconnection rate begins to stall. However, the resulting saturated reconnection rate, while insensitive to the plasma resistivity, was found to depend strongly on the di. In this presentation, we revisit the problem of magnetic island coalescence (which is a well known example of flux pileup reconnection), addressing the dependence of the maximum coalescence rate on the ratio of di in the "large island" limit in which the following inequality is always satisfied: l_eta di lambda, where I_eta is the resistive diffusion length and lambda is the island wavelength.

Dorelli, John C.↗

CanTrilBat and Cantera_apps v1.0 beta

This application can determine the performance and chemical behavior of batteries in 1D when they are cycled. With CanTrilBat, we are developing predictive phenomenological models for battery systems to predict operating performance and rate limiting steps in the performance of battery models. Particular attention is paid to primary and secondary chemistry mechanisms, such as the thermal runaway mechanisms experienced in secondary lithium ion batteries or self-discharge reaction mechanism that all batteries experience to one extent or another. The first application of this model has been for modeling the performance of thermal batteries. However, an implementation for secondary ion batteries is next. CanTrilBat applications solves transient problems involving batteries. It is a 1-D application that represents 3-D physical systems that can be reduced using the porous flow approximation for the anode, cathode, and separator. A control volume formulation is used to track conserved quantities. An operator-split approach is used to calculate the chemistry, diffusion and electronic transport that occurs within cathode and anode particles, allowing for the reduction in code complexity. All jacobian operations in CanTrilBat utilize numerical jacobians. A home grown predictor corrector scheme is used for time step control, and a home grown newton solver is used to relax the equations at each time step. Trilinos is used to solve for the resulting linear systems equation A block text input format is used to initialize options for the CanTriBat program. Within this block are the names of several XML input files, which specify the chemistry mechanism. These XML input files are read by low level Cantera routines, and serve to initialize the electrode and electrolyte chemistry mechanisms and transport properties. A GUI implementation has been contracted out to a university professor, but has not been implemented yet. It?s expected that CanTriBat will have the capability to do adaptive grid refinement in later versions. Grids will be halved based on an complicated set of heuristic criteria, based on how well the mesh captures changes to the first and second derivatives of selected solution components. Primary emphasis is given towards converged profiles for IV performance.

Moffat, Harry K.↗

Unlocking Failure Mechanisms and Improvement of Practical Li–S Pouch Cells through In Operando Pressure Study

For Lithium-sulfur (Li-S) battery commercialization, the study at a pouch cell level is essential, as some problems ignored or deemed minimal at the smaller level could have a greater effect on the performance of the larger pouch cell. Herein, the failure mechanisms of Li-S pouch cells are deeply investigated via in operando pressure analysis. It is found that highly porous structures of cathodes/separators and slow electrolyte diffusion through cathodes/separators can both lead to poor initial wetting. Additionally, Li-metal anode dominates the thickness variation of the whole pouch cell, which is verified by in situ measured pressure variation. Consequently, a real-time approach that combined normalized pressure with dP/dV analysis is proposed and validated to diagnose the morphology evolution of Li-metal anode. Moreover, applied pressure and porosity/tortuosity ratio of the cathode are both identified as independent factors that influence anode performance. In addition to stabilizing anodes, high pressure is proven to improve the cathode connectivity and avoid cathode cracking over cycling, which improves the possibility of developing cathodes with high sulfur mass loading. Here, this work provides insights into Li-S pouch cell design (e.g., cathode and separator) and highlights pathways to improve cell capacity and cycling performance with applied and monitored pressure.

25 ENERGY STORAGE↗

Diffusion Monte Carlo approaches for studying nuclear quantum effects in fluxional molecules

Abstract Diffusion quantum Monte Carlo (DMC) provides a powerful approach for obtaining the ground state energy and wave function of molecules, ions, and molecular clusters. The approach is uniquely well suited for studies of fluxional molecules, which undergo large amplitude vibrational motions even in their ground state. In contrast to the electronic structure problem, where the wave function must be antisymmetric with respect to exchange of any pair of electrons, the wave function for the ground vibrational state is nodeless. This greatly simplifies the application of DMC for vibrational problems. Because there is not a single potential function that can be used to describe the intramolecular and intermolecular interactions in all molecular systems, most methods that are used to describe nuclear quantum effects rely on a carefully chosen zero‐order description of the molecular vibrations. In contrast, DMC calculations can be performed in Cartesian coordinates, making the DMC algorithm easily transferable between different chemical systems. In this contribution, the theory that underlies DMC will be discussed along with important considerations for performing DMC calculations. Extensions for evaluating vibrationally excited states and molecular properties are also discussed. Insights that can be obtained from DMC calculations are illustrated in the context of the protonated water clusters. This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Theoretical and Physical Chemistry > Spectroscopy

Chemistry↗

A novel Mo-W interlayer approach for CVD diamond deposition on steel

Steel is the most widely used material in engineering for its cost/performance ratio and coatings are routinely applied on its surface to further improve its properties. Diamond coated steel parts are an option for many demanding industrial applications through prolonging the lifetime of steel parts, enhancement of tool performance as well as the reduction of wear rates. Direct deposition of diamond on steel using conventional chemical vapour deposition (CVD) processes is known to give poor results due to the preferential formation of amorphous carbon on iron, nickel and other elements as well as stresses induced from the significant difference in the thermal expansion coefficients of those materials. This article reports a novel approach of deposition of nanocrystalline diamond coatings on high-speed steel (M42) substrates using a multi-structured molybdenum (Mo) – tungsten (W) interlayer to form steel/Mo/Mo-W/W/diamond sandwich structures which overcome the adhesion problem related to direct magnetron sputtering deposition of pure tungsten. Surface, interface and tribology properties were evaluated to understand the role of such an interlayer structure. The multi-structured Mo-W interlayer has been proven to improve the adhesion between diamond films and steel substrates by acting as an effective diffusion barrier during the CVD diamond deposition.

36 MATERIALS SCIENCE↗

A Domain Decomposition Model Reduction Method for Linear Convection-Diffusion Equations with Random Coefficients

We developed a domain decomposition model reduction method for linear steady- state convection-diffusion equations with random coefficients. Of particular interest to this effort are the diffusion equation with random diffusivity and the convection-dominated transport equation with random velocity. We investigated two types of random fields, i.e., colored noises and discrete white noises, both of which can lead to high-dimensional parametric dependence in practice. The motivation is to exploit domain decomposition to reduce the parametric dimension of local problems in subdomains, such that an entire parametric map can be approximated with a small number of expensive partial differential equation (PDE) simulations. The new method combines domain decomposition with model reduction and sparse polynomial approximation, so as to simultaneously handle the high dimensionality and irregular behavior of the PDEs under consideration. Here, the advantages of our method lie in three aspects: (i) online-offline decomposition, i.e., the online cost is independent of the size of the triangle mesh; (ii) sparse approximation of operators involving nonaffine high-dimensional random fields; (iii) an effective strategy to capture irregular behaviors, e.g., sharp transitions of the PDE solutions. Two numerical examples are provided to demonstrate the advantageous performance of our method.

97 MATHEMATICS AND COMPUTING↗

Unlocking Failure Mechanisms and Improvement of Practical Li-S Pouch Cells Through in Operando Pressure Study

For Lithium-sulfur (Li-S) battery commercialization, the study at a pouch cell level is essential, as some problems ignored or deemed minimal at the smaller level could have a greater effect on the performance of the larger pouch cell. Herein, the failure mechanisms of Li-S pouch cells are deeply investigated via in operando pressure analysis. It is found that highly porous structures of cathodes/separators and slow electrolyte diffusion through cathodes/separators can both lead to poor initial wetting. Additionally, Li-metal anode dominates the thickness variation of the whole pouch cell, which is verified by in situ measured pressure variation. Consequently, a real-time approach that combined normalized pressure with dP/dV analysis is proposed and validated to diagnose the morphology evolution of Li-metal anode. Moreover, applied pressure and porosity/tortuosity ratio of the cathode are both identified as independent factors that influence anode performance. In addition to stabilizing anodes, high pressure is proven to improve the cathode connectivity and avoid cathode cracking over cycling, which improves the possibility of developing cathodes with high sulfur mass loading. This work provides insights into Li-S pouch cell design (e.g., cathode and separator) and highlights pathways to improve cell capacity and cycling performance with applied and monitored pressure

25 ENERGY STORAGE↗