pith. sign in

arxiv: 2605.26677 · v1 · pith:OPLTCF3Wnew · submitted 2026-05-26 · ⚛️ physics.flu-dyn

A total-Lagrangian vectorial lattice Boltzmann method for finite-strain hyperelastic dynamics

Jingsen Feng , Xu Chu This is my paper

Pith reviewed 2026-06-29 16:06 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords lattice Boltzmann methodhyperelastic dynamicsfinite straintotal Lagrangianvectorial formulationD2Q4 stencilelastodynamicsfirst Piola-Kirchhoff stress
0
0 comments X

The pith

A total-Lagrangian vectorial lattice Boltzmann method solves finite-strain hyperelastic dynamics while preserving the collide-stream structure.

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

The paper constructs a total-Lagrangian vectorial lattice Boltzmann formulation for two-dimensional finite-strain hyperelastic dynamics. The governing equations are rewritten as a conservative first-order system for material velocity and the full deformation gradient. This lets the first Piola-Kirchhoff stress be evaluated locally from the current deformation and enter the lattice only through nonlinear flux moments. A D2Q4 stencil with six-component vector populations matches the state and the two material-coordinate fluxes. The method keeps the local collide-stream operations of standard lattice Boltzmann schemes and supplies second-order initialization, trapezoidal body forcing, velocity quadrature for displacement, and half-way boundary reconstructions.

Core claim

By writing the hyperelastic dynamics as a conservative first-order system for the material velocity and the full deformation gradient, the first Piola-Kirchhoff stress is evaluated locally and incorporated solely through nonlinear flux moments. This representation allows a D2Q4 stencil with six-component vector populations to handle the state and fluxes, adapting the vectorial first-order strategy from linear elastodynamics while preserving the local collide-stream structure of lattice Boltzmann methods.

What carries the argument

The total-Lagrangian vectorial lattice Boltzmann formulation using a D2Q4 stencil with six-component vector populations to match the state and two material-coordinate fluxes, with local first Piola-Kirchhoff stress evaluation via nonlinear flux moments.

Load-bearing premise

The governing equations can be written as a conservative first-order system for the material velocity and the full deformation gradient.

What would settle it

A numerical test on a known finite-strain hyperelastic problem, such as uniaxial tension of a neo-Hookean material, where the simulated deformation gradient or stress deviates from the exact solution by more than the expected truncation error would falsify the method.

Figures

Figures reproduced from arXiv: 2605.26677 by Jingsen Feng, Xu Chu.

Figure 1
Figure 1. Figure 1: Periodic manufactured solution at t = 0.1 on a 1282 grid. The panels compare the exact displacement magnitude, the D2Q4×6 result, the pointwise error, and a centerline trace. 5.03 × 10−4 for first Piola stress, and 5.63 × 10−4 for Cauchy stress; these values are included in the manufactured-solution error summary in [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Grid convergence for the periodic manufactured solution. The left panel reports second-order initialized errors in u, P, and σ. The right panel compares equilibrium initialization with the second-order population initialization. 4.2. Manufactured solution with boundary data The second manufactured-solution test uses the same nonlinear reference field in Eq. (28), logarithmic neo-Hookean material, and body … view at source ↗
Figure 3
Figure 3. Figure 3: Manufactured solution with boundary data. The figure compares all-Dirichlet boundary conditions with mixed Dirichlet– Neumann boundary conditions and reports both full-domain and two-layer-stripped interior errors [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: shows the deformed shape, the reference-error field, and the grid convergence for the two material choices. The raw relative L 2 displacement error shows second-order convergence for the vectorial lattice scheme. On the 1602 grid, the raw relative errors are 5.25 × 10−5 for the SVK case and 8.77 × 10−5 for the neo-Hookean case. The comparison with the published moment-chain data of Müller et al. [22] is re… view at source ↗
Figure 5
Figure 5. Figure 5: shows that the vectorial formulation follows the finite-element displacement field over the full loading interval. The convergence has larger error constants than in uniaxial tension, reflecting the stronger shear-driven boundary gradients in this setup. At n = 160, the raw relative displacement errors are 2.95×10−4 for the SVK case and 5.80 × 10−4 for the neo-Hookean case. As above, the moment-chain compa… view at source ↗
Figure 6
Figure 6. Figure 6: Constitutive stress-response curves for homogeneous affine finite-strain paths. Solid curves are analytical stresses from the constitutive laws in Appendix A; open symbols are D2Q4×6 affine-patch samples on an 802 grid. 1/20 1/40 1/80 1/160 Δx 10 −5 10 −4 EL2(F) (a) uniaxial pure shear simple shear stretch-shear O(Δx 3/2 ) uniaxial pure shear simple shear stretch-shear 10 −5 10 −4 EL2 (b) F P σ [PITH_FULL… view at source ↗
Figure 7
Figure 7. Figure 7: Affine finite-strain patch errors. Panel (a) reports grid refinement of the deformation-gradient error for the four affine paths. Panel (b) summarizes the n = 160 relative errors in F , P , and σ, taking the maximum over the six material laws for each path. J. Feng and X. Chu: Preprint submitted to Elsevier Page 16 of 26 [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Wave speeds about predeformed states. Analytical speeds are obtained from the acoustic tensor Q(n; F0), and open symbols denote phase speeds measured from D2Q4×6 simulations. The two branches correspond to the slow and fast eigenmodes of Eq. (39). Along this path the volumetric contribution vanishes and both constitutive laws give P21(X, t) = µF21(X, t) = µAk cos(kX − cskt). (45) The field in Eq. (42) is t… view at source ↗
Figure 9
Figure 9. Figure 9: Finite-amplitude periodic shear wave for compressible neo-Hookean and logarithmic neo-Hookean materials. Solid curves are the exact travelling-wave solution; open symbols are D2Q4×6 results for shear amplitudes Ak = 0.2 and Ak = 0.5. 5.4. Cantilever bending-wave benchmark The final example considers a bounded-domain wave problem on a cantilever beam. The reference domain is [0, 4] × [0, 1], discretized by … view at source ↗
Figure 10
Figure 10. Figure 10: Transient bending waves in the 160 × 40 neo-Hookean cantilever. The panels show the deformed configuration at t = 2, 5, 8, 12, colored by the velocity magnitude. 6. Conclusions and outlook This work has introduced a total-Lagrangian vectorial lattice Boltzmann formulation for two-dimensional finite-strain hyperelastic dynamics. The governing equations were written as a conservative first-order system for … view at source ↗
Figure 11
Figure 11. Figure 11: Cantilever response compared with a finite-element reference. The present D2Q4×6 errors are measured against an independent Q1 reference using the benchmark definitions of E∞ and E2, and are shown together with the τ = 0.55 results of Müller et al. [22]. benchmarks further demonstrate that the same flux and boundary reconstructions carry over to finite-strain benchmark configurations. The method also repr… view at source ↗
read the original abstract

Inspired by the vectorial lattice Boltzmann method for linear elastodynamics \citep{boolakee2025linear}, we construct a total-Lagrangian vectorial lattice Boltzmann formulation for two-dimensional finite-strain hyperelastic dynamics. The governing equations are first written as a conservative first-order system for the material velocity and the full deformation gradient. This representation separates the kinematic part of the dynamics from the constitutive closure: the first Piola--Kirchhoff stress is evaluated locally from the current deformation gradient and enters the lattice only through nonlinear flux moments. A D2Q4 stencil with six-component vector populations is then used to match the state and the two material-coordinate fluxes. The formulation includes a second-order population initialization, trapezoidally centered body forcing, displacement reconstruction by velocity quadrature, and half-way reconstructions for velocity Dirichlet and Neumann traction boundaries on grid-aligned domains. The resulting method preserves the local collide--stream structure of standard lattice Boltzmann schemes while adapting the vectorial first-order strategy from linear elastodynamics to hyperelastic finite-strain dynamics.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 0 minor

Summary. The paper claims to construct a total-Lagrangian vectorial lattice Boltzmann method for two-dimensional finite-strain hyperelastic dynamics. The governing equations are recast as a conservative first-order system in material velocity and the full deformation gradient F; the first Piola-Kirchhoff stress P(F) is evaluated locally from F and enters the scheme only through nonlinear flux moments. A D2Q4 stencil with six-component vector populations is employed, together with second-order population initialization, trapezoidally centered body forcing, velocity quadrature for displacement reconstruction, and half-way reconstructions for velocity Dirichlet and Neumann traction boundaries on grid-aligned domains. The formulation is asserted to preserve the local collide-stream structure while extending the vectorial first-order strategy previously developed for linear elastodynamics.

Significance. If the method can be shown to be accurate and stable, the separation of kinematic evolution from local constitutive closure would constitute a useful extension of lattice Boltzmann techniques to finite-strain solid mechanics, retaining the locality and parallelism advantages of standard LBM schemes.

major comments (1)
  1. [Abstract and main text] Abstract and main text: the manuscript describes the formulation steps in detail but contains no numerical experiments, convergence studies, error norms, or comparisons against analytical solutions or established finite-element results for hyperelastic test cases. Without such evidence the accuracy, stability, and conservation properties of the scheme for finite-strain problems cannot be assessed.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the detailed review and the recommendation for major revision. We agree that numerical validation is required to substantiate the claims regarding accuracy, stability, and conservation, and we will incorporate such evidence in the revised manuscript.

read point-by-point responses

  1. Referee: [Abstract and main text] Abstract and main text: the manuscript describes the formulation steps in detail but contains no numerical experiments, convergence studies, error norms, or comparisons against analytical solutions or established finite-element results for hyperelastic test cases. Without such evidence the accuracy, stability, and conservation properties of the scheme for finite-strain problems cannot be assessed.

    Authors: We acknowledge that the referee's observation is correct: the submitted manuscript presents the derivation of the total-Lagrangian vectorial LBM formulation but does not contain any numerical experiments, convergence studies, or comparisons. The work was initially submitted as a formulation paper extending the linear elastodynamics approach, with the intention that validation would follow in subsequent work. However, we agree this leaves the practical performance unverified. In the revised version we will add a dedicated numerical results section that includes (i) grid-convergence studies for finite-strain hyperelastic problems (e.g., uniaxial and biaxial tension of a square domain), (ii) L2 error norms against analytical solutions where available, (iii) comparisons with established finite-element results, and (iv) assessments of stability and discrete conservation properties under increasing strain levels. These additions will directly address the referee's concern. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The derivation begins by recasting the hyperelastic equations as a conservative first-order system in material velocity and deformation gradient F, with the first Piola-Kirchhoff stress entering only through nonlinear flux moments; this is a direct mathematical rewriting of the governing PDEs rather than a fit or self-definition. The D2Q4 vectorial populations and collide-stream structure are then adapted from the cited linear-elastodynamics reference, but the adaptation introduces new total-Lagrangian elements, boundary reconstructions, and hyperelastic constitutive evaluation that are not forced by the inputs. No self-citation load-bearing, fitted-input predictions, ansatz smuggling, or renaming of known results occurs within the paper itself. The construction remains self-contained against the stated assumptions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The formulation rests on the standard continuum-mechanics assumption that the first Piola-Kirchhoff stress is a local function of the deformation gradient and on the choice of a conservative first-order system representation; no free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption The governing equations can be written as a conservative first-order system for material velocity and deformation gradient.
    Explicitly stated as the starting point for separating kinematics from constitutive closure.

pith-pipeline@v0.9.1-grok · 5708 in / 1051 out tokens · 37826 ms · 2026-06-29T16:06:18.143348+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Total-Lagrangian vectorial lattice Boltzmann method for finite-strain hyperelasticity with curved boundaries

    physics.comp-ph 2026-06 unverdicted novelty 7.0

    A vectorial lattice Boltzmann scheme for hyperelastic dynamics that resolves geometric mismatch at curved boundaries via level-set embedding, cut-link interpolation, and compatibility projection.

Reference graph

Works this paper leans on

35 extracted references · 2 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    Boolakee, M

    O. Boolakee, M. Geier, L. De Lorenzis, Lattice Boltzmann for linear elastodynamics: Periodic problems and Dirichlet boundary conditions, Comput. Methods Appl. Mech. Engrg. 433 (2025) 117469

    2025

  2. [2]

    McNamara, G

    G.R. McNamara, G. Zanetti, Use of the Boltzmann equation to simulate lattice-gas automata, Phys. Rev. Lett. 61 (20) (1988) 2332–2335

    1988

  3. [3]

    Succi, The Lattice Boltzmann Equation for Fluid Dynamics and Beyond, Clarendon Press, Oxford, 2001

    S. Succi, The Lattice Boltzmann Equation for Fluid Dynamics and Beyond, Clarendon Press, Oxford, 2001

    2001

  4. [4]

    Krüger, H

    T. Krüger, H. Kusumaatmaja, A. Kuzmin, O. Shardt, G. Silva, E.M. Viggen, The Lattice Boltzmann Method: Principles and Practice, Springer, Cham, 2017

    2017

  5. [5]

    Haslam, R.S

    I.W. Haslam, R.S. Crouch, M. Seaïd, Coupled finite element–lattice Boltzmann analysis, Comput. Methods Appl. Mech. Engrg. 197 (51–52) (2008) 4505–4511

    2008

  6. [6]

    H. Zhou, G. Mo, F. Wu, J. Zhao, M. Rui, K. Cen, GPU implementation of lattice Boltzmann method for flows with curved boundaries, Comput. Methods Appl. Mech. Engrg. 225–228 (2012) 65–73

    2012

  7. [7]

    Verdier, P

    W. Verdier, P. Kestener, A. Cartalade, Performance portability of lattice Boltzmann methods for two-phase flows with phase change, Comput. Methods Appl. Mech. Engrg. 370 (2020) 113266

    2020

  8. [8]

    Zhang, S

    P. Zhang, S. Sun, Y. Chen, S.A. Galindo-Torres, W. Cui, Coupled material point lattice Boltzmann method for modeling fluid–structure interactions with large deformations, Comput. Methods Appl. Mech. Engrg. 385 (2021) 114040

    2021

  9. [9]

    Feng, J.-C

    J.-S. Feng, J.-C. Min, Lattice Boltzmann method simulation of two-phase flow in horizontal channel, Acta Phys. Sin. 72 (8) (2023) 084701

    2023

  10. [10]

    Y. Liu, J. Feng, J. Min, X. Zhang, Homogenized color-gradient lattice Boltzmann model for immiscible two-phase flow in multiscale porous media, J. Appl. Phys. 135 (18) (2024) 184701

    2024

  11. [11]

    J. Feng, C. Wang, Y. Guo, Y. Liu, J. Min, M. Wang, K. Xu, Helmholtz-guided suppression of spurious currents in multicomponent pseudopotential model for high-fidelity immiscible flow in porous media, Phys. Fluids 38 (2) (2026) 026611

    2026

  12. [12]

    J. Feng, J. Leng, J. Jiang, X. Chu, Entropic lattice Boltzmann method for general anisotropic advection–diffusion, arXiv:2605.01774, 2026

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  13. [13]

    X. Yin, G. Yan, T. Li, Direct simulations of the linear elastic displacements field based on a lattice Boltzmann model, Internat. J. Numer. Methods Engrg. 107 (3) (2016) 234–251

    2016

  14. [14]

    Boolakee, M

    O. Boolakee, M. Geier, L. De Lorenzis, A new lattice Boltzmann scheme for linear elastic solids: periodic problems, Comput. Methods Appl. Mech. Engrg. 404 (2023) 115756

    2023

  15. [15]

    Boolakee, M

    O. Boolakee, M. Geier, L. De Lorenzis, Dirichlet and Neumann boundary conditions for a lattice Boltzmann scheme for linear elastic solids on arbitrary domains, Comput. Methods Appl. Mech. Engrg. 415 (2023) 116225

    2023

  16. [16]

    Marconi, B

    S. Marconi, B. Chopard, A lattice Boltzmann model for a solid body, Internat. J. Modern Phys. B 17 (1–2) (2003) 153–156

    2003

  17. [17]

    O’Brien, T

    G.S. O’Brien, T. Nissen-Meyer, C.J. Bean, A lattice Boltzmann method for elastic wave propagation in a Poisson solid, Bull. Seismol. Soc. Am. 102 (3) (2012) 1224–1234

    2012

  18. [18]

    Murthy, P.K

    J.S.N. Murthy, P.K. Kolluru, V. Kumaran, S. Ansumali, Lattice Boltzmann method for wave propagation in elastic solids, Commun. Comput. Phys. 23 (4) (2018) 1223–1240

    2018

  19. [19]

    Escande, P.K

    M. Escande, P.K. Kolluru, L.M. Cléon, P. Sagaut, Lattice Boltzmann method for wave propagation in elastic solids with a regular lattice: theoretical analysis and validation, arXiv:2009.06404, 2020

    work page arXiv 2009

  20. [20]

    Schlüter, C

    A. Schlüter, C. Kuhn, R. Müller, Lattice Boltzmann simulation of antiplane shear loading of a stationary crack, Comput. Mech. 62 (5) (2018) 1059–1069

    2018

  21. [21]

    Faust, A

    E. Faust, A. Schlüter, H. Müller, F. Steinmetz, R. Müller, Dirichlet and Neumann boundary conditions in a lattice Boltzmann method for elastodynamics, Comput. Mech. 73 (2) (2024) 317–339

    2024

  22. [22]

    Müller, E

    H. Müller, E. Faust, A. Schlüter, R. Müller, Extending the lattice Boltzmann method to non-linear elastodynamics, Comput. Methods Appl. Mech. Engrg. 443 (2025) 118076. J. Feng and X. Chu:Preprint submitted to ElsevierPage 25 of 26 Total-Lagrangian vectorial LBM for finite-strain hyperelastic dynamics

    2025

  23. [23]

    Dellar, Moment equations for magnetohydrodynamics, J

    P.J. Dellar, Moment equations for magnetohydrodynamics, J. Stat. Mech. Theory Exp. 2009 (6) (2009) P06003

    2009

  24. [24]

    S. Jin, Z. Xin, The relaxation schemes for systems of conservation laws in arbitrary space dimensions, Comm. Pure Appl. Math. 48 (3) (1995) 235–276

    1995

  25. [25]

    Graille, Approximation of mono-dimensional hyperbolic systems: A lattice Boltzmann scheme as a relaxation method, J

    B. Graille, Approximation of mono-dimensional hyperbolic systems: A lattice Boltzmann scheme as a relaxation method, J. Comput. Phys. 266 (2014) 74–88

    2014

  26. [26]

    Dubois, Simulation of strong nonlinear waves with vectorial lattice Boltzmann schemes, Internat

    F. Dubois, Simulation of strong nonlinear waves with vectorial lattice Boltzmann schemes, Internat. J. Modern Phys. C 25 (12) (2014) 1441014

    2014

  27. [27]

    Zhao, W.A

    J. Zhao, W.A. Yong, Vectorial finite-difference-based lattice Boltzmann method: Consistency, boundary schemes and stability analysis, J. Comput. Appl. Math. 441 (2024) 115677

    2024

  28. [28]

    Bhatnagar, E.P

    P.L. Bhatnagar, E.P. Gross, M. Krook, A model for collision processes in gases. I. Small amplitude processes in charged and neutral one-component systems, Phys. Rev. 94 (1954) 511–525

    1954

  29. [29]

    Ogden, Non-Linear Elastic Deformations, Dover Publications, Mineola, 1997

    R.W. Ogden, Non-Linear Elastic Deformations, Dover Publications, Mineola, 1997

    1997

  30. [30]

    Holzapfel, Nonlinear Solid Mechanics: A Continuum Approach for Engineering, Wiley, Chichester, 2000

    G.A. Holzapfel, Nonlinear Solid Mechanics: A Continuum Approach for Engineering, Wiley, Chichester, 2000

    2000

  31. [31]

    Bonet, R.D

    J. Bonet, R.D. Wood, Nonlinear Continuum Mechanics for Finite Element Analysis, second ed., Cambridge University Press, Cambridge, 2008

    2008

  32. [32]

    Hughes, The Finite Element Method: Linear Static and Dynamic Finite Element Analysis, Dover Publications, Mineola, 2000

    T.J.R. Hughes, The Finite Element Method: Linear Static and Dynamic Finite Element Analysis, Dover Publications, Mineola, 2000

    2000

  33. [33]

    M. Junk, A. Klar, L.S. Luo, Asymptotic analysis of the lattice Boltzmann equation, J. Comput. Phys. 210 (2) (2005) 676–704

    2005

  34. [34]

    Banda, W.A

    M.K. Banda, W.A. Yong, A. Klar, A stability notion for lattice Boltzmann equations, SIAM J. Sci. Comput. 27 (6) (2006) 2098–2111

    2006

  35. [35]

    Junk, W.A

    M. Junk, W.A. Yong, WeightedL 2-stability of the lattice Boltzmann method, SIAM J. Numer. Anal. 47 (3) (2009) 1651–1665. J. Feng and X. Chu:Preprint submitted to ElsevierPage 26 of 26

    2009