Pith. sign in

REVIEW 2 major objections 5 minor 21 references

Gas-induced bulging in pouch-cell batteries: a mechanical model

T0 review · 2 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The bulge shape of a gas-filled pouch cell is controlled by two dimensionless numbers, and matching it to one cathode layer predicts internal pressure and gas volume.

desk verdict Clean shape model for pouch-cell bulging with a solid pressure estimate, but the gas-volume equation (38) is internally inconsistent by a factor of 2-3. read the letter →

arxiv 2411.13197 v1 pith:EZYVPRRN submitted 2024-11-20 cond-mat.soft

Andrea Giudici , Colin Please , Jon Chapman This is my paper · ORCID
classification cond-mat.soft
keywords pouchcellbatteriesgas-inducedbulgingmechanicalmodelWinklerfoundationhomogenisedelasticitystateofhealthX-raytomographydimensionlessparameters
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper claims that the bulge shape of a gas-filled lithium-ion pouch cell is set by two dimensionless parameters: one measuring how far bending deformations persist from the clamped edges, and one measuring the through-cell strain produced by internal pressure. Fitting the predicted shape to the top cathode layer in X-ray tomography images from a cycled cell yields quantitative estimates of internal pressure and gas volume. The authors derive a homogenised elastic model in which stiff cathode and current-collector layers bend while soft anode layers act as springs, leading to a closed-form shape formula. If correct, the model lets a single external image of a bulge serve as a non-invasive state-of-health monitor for sealed batteries.

What carries the argument

The central object is the homogenised energy functional in which each anode is a Winkler foundation, a bed of independent springs of stiffness $K$, and each cathode and current collector is an inextensible bending sheet of stiffness $B$. Minimising this energy leads to the fourth-order PDE above, with two dimensionless parameters: $\gamma$, the bending-persistence ratio, and $\epsilon_n$, the pressure-induced strain. The same equation also governs a single bending sheet on an elastic foundation, so the multi-layer result is built on that classical solution; the machinery converts a complex layered deformation into a one-parameter family of shapes.

What would settle it

Instrument a pouch cell with a pressure sensor, image the bulge as gas is injected, and compare the pressure inferred from the fitted shape with the direct sensor reading; systematic disagreement at large strain would falsify the linear-elastic foundation assumption.

Watch

Extended reading notes

Core claim

For a pouch cell with many thin layers, the through-cell displacement field $v(x,y)$ obeys a homogenised plane-strain equation $\partial^4 \bar v/\partial \bar x^4 = 4\gamma^4 \, \partial^2 \bar v/\partial \bar y^2$, whose closed-form solution is a Fourier series with coefficients set by the boundary conditions at the casing. The shape is fully determined by $\gamma = (\hat K W^4/4\hat B T^2)^{1/4}$, the ratio of cell width to a bending-persistence length, and the amplitude is set by $\epsilon_n = P/\hat K$. Fitting this solution to the experimental data at cycles 100, 150, and 200 gives $\gamma=3.21$ and strains $\epsilon_n=0.41,0.62,0.77$; with an assumed cathode bending stiffness this implies an effective anode stiffness $K=92.5$ kPa and pressures $P=37.2,57.2,71.1$ kPa, about one atmosphere at the latest cycle. Although only the top cathode layer was fitted, all other layers are reproduced, and the pressure follows a super-linear relation with gas volume through the ideal gas law.

Load-bearing premise

The anodes behave as linear springs even when they are stretched or filled with gas to roughly 70 percent strain.

Editorial extensions

If this is right

  • Fitting only the outer cathode shape gives estimates of internal pressure and gas amount, enabling state-of-health monitoring without opening the cell.
  • Gas production grows super-linearly with pressure, so tracking pressure alone underestimates chemical degradation.
  • The predicted stress and bending-moment fields show where mechanical damage concentrates: tensile stress in the middle and bending moments near the clamped edges.
  • The homogenised model can be embedded in larger battery simulations to capture mechanical degradation over cycling.
  • The fitted effective anode stiffness is far below pristine values, indicating that tensile softening and structural damage are absorbed into one effective spring constant.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the shape is measured across cycles, the fitted $\gamma$ could track how the effective stiffness $K$ degrades over time, turning a static fit into a dynamic degradation monitor.
  • The linear-spring assumption is least secure at the observed 70% strains; a nonlinear foundation or poroelastic anode model would likely change the pressure-volume relation at high gas content.
  • The same dimensionless reduction may apply to other sealed layered systems that swell and buckle, such as composite panels or packaged food, with different constitutive laws for the soft core.
  • A direct test would be to instrument a pouch cell with a pressure port, image it while injecting gas, and compare the inferred pressure with the sensor reading.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. The paper proposes a homogenised mechanical model for gas-induced bulging of pouch-cell batteries. The cell is represented as alternating bending sheets (cathodes/current collectors) connected by Winkler foundations (anodes). In the many-layer limit the authors derive a PDE whose analytical solution gives the bulge shape in terms of two dimensionless parameters, γ (bending persistence length) and ε_n (pressure-induced through-cell strain). They fit the shape to published X-ray tomography data, obtain γ=3.21 and strains 0.41–0.77 at cycles 100–200, and, using an assumed cathode bending stiffness, estimate internal pressures of 37–71 kPa. They also propose Eq. (38) to infer the amount of gas from pressure and shape. The homogenised solution is compared with a numerical minimisation of the discrete model, and the agreement is good for γ>3. The central claims are that the bulge shape is governed by two parameters and that fitting a single layer yields quantitative pressure and gas-volume estimates for non-invasive state-of-health monitoring.

Significance. If correct, the model provides a simple, analytically solvable framework for non-invasive SOH monitoring, and the shape prediction is a falsifiable, mechanically grounded result. The transparent derivation, the closed-form solution, and the systematic comparison with the discrete numerical model for n=5 are strengths. However, the gas-volume relation contains an internal inconsistency, and the linear-elastic Winkler-foundation assumption is applied far beyond its small-strain range; these issues currently limit the quantitative claims and require revision.

major comments (2)
  1. [Estimate of the pressure and gas production, Eq. (38)] The derivation of Eq. (38) implicitly sets V = V0 + ΔV in the ideal gas law, i.e., it assumes the gas occupies the entire pristine cell volume plus the expansion volume. In a pouch cell the gas is confined to the void and delamination volume, not to the solid electrode volume, so the correct gas volume is approximately ΔV (plus initial porosity). With V ≈ ΔV = V0 ε_n g(γ), the dimensionless gas amount becomes n_g(V0 K̂/RT) = ε_n^2 g(γ), not ε_n(1 + ε_n g(γ)) as written. For the fitted strains 0.41–0.77 and g(γ) of order one, the two expressions differ by a factor of order 2–3, so the reported gas amounts and the SOH-monitoring claim are quantitatively affected. This is an internal inconsistency, not merely a parameter-uncertainty issue, and the equation should be corrected and the experimental gas amounts recomputed.
  2. [Mechanical model and Appendix A] The Winkler-foundation representation is derived under small strains, small gradients, and thin-layer assumptions in Appendix A, but the experimental data cited in the introduction show through-cell strains up to 70% and large gas pockets in the anode layers. The fitted K = 92.5 kPa is three orders of magnitude below the pristine anode modulus, indicating that K is an effective parameter that absorbs the observed degradation and damage. The paper should state the model's intended range of validity and either incorporate a nonlinear foundation or explicitly frame K and the resulting pressures and gas amounts as effective quantities, with a sensitivity analysis over the assumed bending stiffness and the strain range.
minor comments (5)
  1. [Eq. (9)] The expression for b is missing the denominator: it should read b = 2 sin(γ1/2) sinh(γ1/2) / (cos γ1 + cosh γ1).
  2. [Eq. (12)] The notation "1/2t" is ambiguous and should be written as "1/(2t)" to avoid confusion with (1/2)t.
  3. [Single layer and homogenised model, boundary conditions] The boundary condition v = 0 and v'' = 0 at x = ±W/2 is called "clamped ends and no bending moment," but zero bending moment corresponds to a pinned or simply supported end, not a clamped end. The terminology should be corrected, since the physical condition is only zero displacement at the casing.
  4. [Appendix A, Eq. (48)] The derivation of the Winkler energy mixes notation for the length L, width W, and the z-direction; the integration domain should be written explicitly to avoid ambiguity about whether the energy is per unit length or integrated over the full cell.
  5. [Comparison with experiments] The pressure estimates rely on an assumed cathode bending stiffness B_C in the range 1–10 GPa, which introduces at least a factor-of-10 uncertainty in the derived K and P. The paper should state confidence intervals or at least discuss how the reported pressures depend on this assumption.

Circularity Check

1 steps flagged · score 6.0 of 10

Pressure and gas-amount predictions are algebraic rescalings of fitted strain and assumed stiffness; the analytic shape solution itself is independently derived.

  1. fitted input called prediction [Section 'Comparison with experiments and prediction of the pressure', subsection 'Estimate of the pressure and gas production'; Eqs. (19), (21), (35)]
    "We fit our model to the experimental data in [4] to obtain γ and three values of the through-cell strain at the different cycle number ... γ = 3.21, ε100n = 0.41, ε150n = 0.62, ε200n = 0.77. ... We can use equation (21) and γ = 3.21 from our fit to the experimental data to find that the stiffness of the substrate and the pressure in the system for the 100th, 150th and 200th cycles are: K = 92.5 kPa, P100 = 37.2 kPa, P150 = 57.2 kPa, P200 = 71.1 kPa."

    The pressure is not independently measured or predicted. From Eq. (19), v = εn T v̄, so the fitted εn is the through-cell strain amplitude, i.e. the shape-fit parameter. From Eq. (21), γ determines K̂ once the bending stiffness B̂ is assumed, and the paper then uses P = εn K̂ to report pressure. Thus P100, P150 and P200 are exactly the fitted εn values multiplied by a constant derived from the same fitted γ and an assumed bending stiffness; they carry no information beyond those inputs. The gas amount in Eq. (38) is likewise an algebraic function of the same fitted εn, γ and assumed K̂.

full rationale

The core mechanical derivation is not circular: the homogenised energy (13) is minimised to give the PDE (20) and the closed-form shape (34), with no pressure or gas datum used as an input, and the numerical discrete minimisation compared in Fig. 3 provides an independent check. The fitted parameters γ and εn are legitimate calibration variables of the model. The circularity-adjacent step is the interpretation of the pressure estimate: Eq. (35) is simply P = εn K̂, with K̂ obtained from the same fitted γ and an assumed bending stiffness, so the 'predicted' pressure is a re-parameterization of the two fitted numbers rather than an independent estimate. The gas-amount relation (38) is a genuine derived consequence of the ideal-gas law, but it inherits the same fitted/assumed inputs and is not externally validated; additionally, its use of V0 as part of the gas volume is a physical modeling error, since gas occupies the expansion (and void) volume rather than the pristine solid volume, though this is a correctness issue rather than a circularity. The only self-citations ([1], [2]) are contextual and not load-bearing, and no uniqueness theorem is imported from the authors' prior work. Overall, the analytic mechanics is self-contained, but the abstract and conclusion overstate the pressure and gas 'predictions' as if they were independent outputs, justifying a partial circularity score.

Assumptions & free parameters 6 free parameters · 6 assumptions · 0 invented entities

The model introduces no new physical entities. Its quantitative outputs depend on fitted dimensionless parameters (gamma and three strains) and on assumed material inputs (cathode bending stiffness and anode fraction), which are the main free parameters. The physical assumptions are standard for this class of homogenised mechanical models, but the use of linear springs at large strains is the most fragile input.

free parameters (6)
  • gamma = 3.21
    Dimensionless ratio W/l_n fitted to the experimental top-layer shape at all three cycle numbers simultaneously.
  • epsilon_n (cycle 100) = 0.41
    Through-cell strain at cycle 100, fitted to the top cathode layer shape.
  • epsilon_n (cycle 150) = 0.62
    Through-cell strain at cycle 150, fitted to the top cathode layer shape.
  • epsilon_n (cycle 200) = 0.77
    Through-cell strain at cycle 200, fitted to the top cathode layer shape.
  • cathode bending stiffness BC = about 7e-5 Pa m^2
    Assumed from cathode Young's modulus 1 to 10 GPa and thickness up to 100 micrometers (Appendix B). The pressure estimate scales with this value.
  • anode fraction phi = about 0.5
    Assumed fraction of the layer thickness occupied by the anode, used to convert the fitted substrate stiffness to the homogenised stiffness K_hat.
assumptions (6)
  • domain assumption Anode layers behave as linear Winkler foundations with a constant spring stiffness K.
    Appendix A derives this under the assumptions of thin anode layers and small gradients, but experimental strains reach 70%, so linear elasticity may be violated.
  • domain assumption Cathode layers and current collectors are inextensible and do not change thickness.
    Based on the experimental observation that the in-plane dimension does not change and cathode thickness variation is small (Section 'Mechanical model').
  • domain assumption The deformation is two-dimensional plane strain, independent of the z coordinate in the bulk.
    Justified by L > 2W and the boundary layer length ln = W/gamma being much smaller than L (about 7 mm versus 49 mm).
  • domain assumption Displacements and slopes are small enough to use the linearized curvature v'' and small strain approximations.
    Assumed throughout the energy derivation, but the experimental bulge has large local strains in the anodes, potentially violating this for the substrate.
  • domain assumption The gas inside the cell obeys the ideal gas law at constant temperature.
    Equation (36) assumes P V = n_g R T with T constant, which is reasonable for the pressure range but not explicitly verified.
  • standard math The homogenisation limit n -> infinity with t -> 0 is valid, and the alternating bending stiffness term vanishes by the alternating series test.
    Used to pass from the discrete energy (12) to the continuous energy (13); the validity is checked numerically against the discrete model.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Gas-induced bulging in pouch-cell batteries: a mechanical model." pith.science (2026). https://pith.science/paper/EZYVPRRN

@misc{pith2026241113197,
  author       = {Pith},
  title        = {Pith review of: Gas-induced bulging in pouch-cell batteries: a mechanical model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EZYVPRRN}},
  note         = {Machine review of arXiv:2411.13197}
}
read the original abstract

Over the long timescale of many charge/discharge cycles, gas formation can result in large bulging deformations of a Lithium-ion pouch cell, which is a key failure mechanism in batteries. Guided by recent experimental X-ray tomography data of a bulging cell, we propose a homogenised mechanical model to predict the shape of the deformation and the stress distribution analytically. Our model can be included in battery simulation models to capture the effects of mechanical degradation. Furthermore, with knowledge of the bending stiffness of the cathode electrodes and current collectors, and by fitting our model to experimental data, we can predict the internal pressure and the amount of gas in the battery, thus assisting in monitoring the state of health (SOH) of the cell without breaking the sealed case.

Figures

Figures reproduced from arXiv: 2411.13197 by the authors.

Figure 1
Figure 1. FIG. 1. a) A 3D schematic of the deformation in the pouch cell due to gas pressure. b) Side by side of a schematic view of the [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. At the top, we show the shape of a sheet attached to a soft substrate subject to no displacement boundary condition, [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Comparison of the shape of the top layer between the continuous (homogenised) model and the (numerical) discrete [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: FIG. 4. a) Comparison between the experimental data at different charge cycle and the fitted theoretical model, where orange [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. (a) The relationship between the amount of gas [PITH_FULL_IMAGE:figures/full_fig_p011_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Schematics of the discretisation used to minimise the energy in equation (10). [PITH_FULL_IMAGE:figures/full_fig_p013_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

21 extracted references · 9 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archive author booktitle chapter doi edition editor eid eprint howpublished institution journal key keywords month note number organization pages publisher school series title type url volume year archivePrefix primaryClass adsurl adsnote version label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sent...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION add.period duplicate empty 'skip "." * add.blank if FUNCTION if.digit duplicate "0" = swap duplicate "1" = swap duplicate "2" = swap duplicate "3" = swap duplicate "4" = swap duplicate "5" = swap duplicate "6" = swap duplicate "7" = swap duplicate "8" = swap "9" = or or or or or or or or or FUNCTION ...

  3. [3]

    write newline

    " write newline " cite write " FUNCTION editor.postfix editor num.names #1 > "( )" "( )" if FUNCTION editor.trans.postfix editor num.names #1 > "( )" "( )" if FUNCTION trans.postfix translator num.names #1 > "( )" "( )" if FUNCTION authors.editors.reflist.apa5 'field := 'dot := field num.names 'numnames := numnames 'format.num.names := format.num.names na...

  4. [4]

    Available from:

    ENTRY address assignee author booktitle chapter cartographer day edition editor howpublished institution inventor journal key keywords month note number organization pages part publisher school series title type volume word year eprint doi url lastchecked updated archive archivePrefix primaryClass eid adsurl adsnote version label INTEGERS output.state bef...

  5. [5]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize ":" * " " *...

  6. [6]

    write newline

    " write newline "" before.all 'output.state := FUNCTION string.to.integer 't := t text.length 'k := #1 'char.num := t char.num #1 substring 's := s is.num s "." = or char.num k = not and char.num #1 + 'char.num := while char.num #1 - 'char.num := t #1 char.num substring FUNCTION find.integer 't := #0 'int := int not t empty not and t #1 #1 substring 's :=...

  7. [7]

    write newline

    " write newline "" before.all 'output.state := FUNCTION output.doi doi empty skip "doi:" doi * "" * output if FUNCTION format.archive archivePrefix empty "" archivePrefix ":" * if FUNCTION format.primaryClass primaryClass empty "" " [" primaryClass * "] " * if FUNCTION format.eprint eprint empty "" archive empty " https://arxiv.org/abs/" eprint * " " * " ...

  8. [8]

    write newline

    " write newline "" before.all 'output.state := FUNCTION string.to.integer 't := t text.length 'k := #1 'char.num := t char.num #1 substring 's := s is.num s "." = or char.num k = not and char.num #1 + 'char.num := while char.num #1 - 'char.num := t #1 char.num substring FUNCTION find.integer 't := #0 'int := int not t empty not and t #1 #1 substring 's :=...

Show all 21 references
  1. [9]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter doi edition editor eid howpublished institution journal key keywords month note number organization pages publisher school series title type url volume year eprint archive archivePrefix primaryClass adsurl adsnote version label INTEGERS o...

  2. [10]

    write newline

    " write newline "" before.all 'output.state := FUNCTION if.digit duplicate "0" = swap duplicate "1" = swap duplicate "2" = swap duplicate "3" = swap duplicate "4" = swap duplicate "5" = swap duplicate "6" = swap duplicate "7" = swap duplicate "8" = swap "9" = or or or or or or...

  3. [11]

    , " * write output.state after.block = add.period write newline

    ENTRY address archive author booktitle chapter edition editor eprint howpublished institution journal key keywords month note number organization pages publisher school series title type url doi volume year archivePrefix primaryClass eid adsurl adsnote version label INTEGERS o...

  4. [12]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 gl...

  5. [13]

    Timms , author S

    author author R. Timms , author S. Psaltis , author C. P. \ Please ,\ and\ author S. J. \ Chapman ,\ @noop journal journal Journal of the Mechanics and Physics of Solids \ volume 175 ,\ pages 105269 ( year 2023 ) NoStop

    2023

  6. [14]

    Giudici , author J

    author author A. Giudici , author J. Chapman ,\ and\ author C. Please ,\ @noop journal journal arXiv preprint arXiv:2407.00373 \ ( year 2024 ) NoStop

    2024 arXiv

  7. [15]

    author author S. S. \ Zhang ,\ @noop journal journal Frontiers in energy research \ volume 2 ,\ pages 59 ( year 2014 ) NoStop

    2014

  8. [16]

    Du , author R

    author author W. Du , author R. E. \ Owen , author A. Jnawali , author T. P. \ Neville , author F. Iacoviello , author Z. Zhang , author S. Liatard , author D. J. \ Brett ,\ and\ author P. R. \ Shearing ,\ @noop journal journal Journal of Power Sources \ volume 520 ,\ pages 23...

    2022

  9. [17]

    Kim , author B

    author author J. Kim , author B. Gerelt-Od , author E. Shin , author H. Kang , author N. Kim , author C. Jo , author H. Son ,\ and\ author S. Yoon ,\ @noop journal journal Journal of Electroanalytical Chemistry \ volume 907 ,\ pages 115892 ( year 2022 ) NoStop

    2022

  10. [18]

    \"U c el , author P

    author author \.I . \"U c el , author P. Gupta ,\ and\ author P. Gudmundson ,\ @noop journal journal Journal of Energy Storage \ volume 46 ,\ pages 103845 ( year 2022 ) NoStop

    2022

  11. [19]

    Gupta , author I

    author author P. Gupta , author I. \"U c el , author P. Gudmundson ,\ and\ author E. Olsson ,\ @noop journal journal Experimental Mechanics \ volume 60 ,\ pages 847 ( year 2020 ) NoStop

    2020

  12. [20]

    Schiele , author T

    author author A. Schiele , author T. Hatsukade , author B. B. \ Berkes , author P. Hartmann , author T. Brezesinski ,\ and\ author J. Janek ,\ @noop journal journal Analytical chemistry \ volume 89 ,\ pages 8122 ( year 2017 ) NoStop

    2017

  13. [21]

    author author L. D. \ Landau , author E. M. \ Lifshitz , author R. Atkin ,\ and\ author N. Fox ,\ in\ @noop booktitle Physics of Continuous Media \ ( publisher CRC Press ,\ year 2020 )\ pp.\ pages 167--178 NoStop

    2020

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.