pith. sign in

arxiv: 2104.06000 · v1 · pith:ZKWQO4RRnew · submitted 2021-04-13 · 🧮 math-ph · math.MP

A Geometric Formulation of Linear Elasticity Based on Discrete Exterior Calculus

Pieter D. Boom , Odysseas Kosmas , Lee Margetts , Andrey Jivkov This is my paper

Pith reviewed 2026-05-24 12:43 UTC · model grok-4.3

classification 🧮 math-ph math.MP
keywords discrete exterior calculuslinear elasticitycell complexeslattice structuresmusical isomorphismsHodge starconstitutive relationsdiscrete differential geometry
0
0 comments X

The pith

Linear elasticity on cell complexes is formulated directly using discrete exterior calculus, with the macroscopic constitutive relation enforced at primal 0-cells via musical isomorphisms.

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

The paper develops a geometric approach to linear elasticity for cell complexes using discrete exterior calculus. Displacements are treated as primal vector-valued 0-cochains, their differences as 1-cochains, and internal forces as dual 2-cochains. The key step is enforcing the macroscopic constitutive relation at each primal 0-cell through musical isomorphisms that convert between discrete cochains and smooth fields. This setup reduces the problem to solving a Laplace equation involving a non-local material Hodge star operator. If successful, it provides a way to derive consistent force-displacement relations for any lattice that reproduces a given macroscopic elastic response.

Core claim

The authors present a direct formulation of linear elasticity based on discrete exterior calculus where the primary unknowns are displacements as primal vector-valued 0-cochains. Displacement differences are primal vector-valued 1-cochains and internal forces are dual vector-valued 2-cochains. The macroscopic constitutive relation is enforced at primal 0-cells using musical isomorphisms to map cochains to smooth fields. Balance of linear momentum is established at the same points, and the system is solved as a Laplace equation with a non-local and non-diagonal material Hodge star. Numerical tests on classical problems confirm agreement with analytic solutions, establishing a method to relate

What carries the argument

Musical isomorphisms that map primal vector-valued cochains to smooth vector fields at 0-cells to enforce the macroscopic constitutive relation, combined with the resulting non-local non-diagonal material Hodge star.

If this is right

  • Relations between displacement differences and internal forces can be calculated for any lattice structure required to follow a prescribed macroscopic elastic behaviour.
  • The governing equations reduce to a Laplace equation with a non-local and non-diagonal material Hodge star.
  • Balance of linear momentum is satisfied at primal 0-cells.
  • Numerical simulations of classical problems with analytic solutions show good agreement with known results.
  • This constitutes the first step toward formulations for dissipative processes in cell complexes.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The approach could support direct design of lattice geometries that achieve target continuum elastic properties without separate homogenization.
  • The non-local Hodge star might naturally incorporate long-range effects when modeling heterogeneous or graded discrete materials.
  • Similar machinery could be adapted for other tensor fields such as thermal conductivity on the same cell complexes.

Load-bearing premise

Musical isomorphisms can be used to map discrete cochains to smooth fields at primal 0-cells in a manner that correctly enforces the macroscopic constitutive relation without introducing uncontrolled approximation error for general cell complexes.

What would settle it

For a chosen lattice and prescribed macroscopic elastic tensor, compute the derived force-displacement relations and test whether the lattice under boundary conditions consistent with a uniform strain reproduces the target macroscopic response to within discretization tolerance.

Figures

Figures reproduced from arXiv: 2104.06000 by Andrey Jivkov, Lee Margetts, Odysseas Kosmas, Pieter D. Boom.

Figure 1
Figure 1. Figure 1: Cylinder and corresponding primal simplicial complex [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Two-dimensional simplices shown with vertex numbering, orientation, and imposed orientation [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Two-dimensional primal simplical complex (red) and dual cell complex (blue). [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: de Rham complex for three dimensional primal and dual complexes [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Relationship between de Rham complex, vector fields, and vector calculus operations. [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Closure for the divergence discrete exterior derivative ( [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Primal complex of circular cylinder of the system matrix. The boundary values are removed from the solution vector, as well as the corresponding columns of the system matrix, and the final system is solved as an overdetermined problem. 5.1. Compression of a circular cylinder The first simulation is the compression of a circular cylinder in order to test the ability of the proposed theory to recover the ful… view at source ↗
Figure 8
Figure 8. Figure 8: Computed Poisson’s ratio for the compression of a circular cylinder [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Computed axial shear stress as a function of radius for the twist of a circular cylinder [PITH_FULL_IMAGE:figures/full_fig_p020_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: The max error in displacements on the coarse grid is on the order of [PITH_FULL_IMAGE:figures/full_fig_p020_10.png] view at source ↗
Figure 10
Figure 10. Figure 10: Computed displacement magnitude divided by height as a function of radius for the twist of a [PITH_FULL_IMAGE:figures/full_fig_p021_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Primal complex of cantilever beam 5.3. Point load on a cantilever beam Next, the deflection of a clamped square section cantilever beam with a single point load at the tip is simulated. In this case, the geometry has both a reasonable amount of displacement and rotation. The beam has a length of 10 as well as both a height and depth of 1. A displacement of −0.19138756 is prescribed to the upper surface at… view at source ↗
Figure 12
Figure 12. Figure 12: Computed deflection of a cantilever beam [PITH_FULL_IMAGE:figures/full_fig_p022_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Primal complex of one-eighth 3D Kirsch problem [PITH_FULL_IMAGE:figures/full_fig_p023_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Computed normal stress for the 3D Kirsch problem on the cube’s mid plane as a function of [PITH_FULL_IMAGE:figures/full_fig_p023_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Relative error in normal stress for the 3D Kirsch problem on the cube’s mid plane as a function [PITH_FULL_IMAGE:figures/full_fig_p024_15.png] view at source ↗
read the original abstract

A direct formulation of linear elasticity of cell complexes based on discrete exterior calculus is presented. The primary unknown are displacements, represented by primal vector-valued 0-cochain. Displacement differences and internal forces are represented by primal vector-valued 1-cochain and dual vector-valued 2-cochain, respectively. The macroscopic constitutive relation is enforced at primal 0-cells with the help of musical isomorphisms mapping cochains to smooth fields and vice versa. The balance of linear momentum is established at primal 0-cells. The governing equations are solved as a Laplace equation with a non-local and non-diagonal material Hodge star. Numerical simulations of several classical problems with analytic solutions are presented to validate the formulation. Good agreement with known solutions is obtained. The formulation provides a method to calculate the relations between displacement differences and internal forces for any lattice structure, when the structure is required to follow a prescribed macroscopic elastic behaviour. This is also the first and critical step in developing formulations for dissipative processes in cell complexes.

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

2 major / 0 minor

Summary. The paper presents a direct DEC-based formulation of linear elasticity on cell complexes. Displacements are represented as primal vector-valued 0-cochains, displacement differences as primal vector-valued 1-cochains, and internal forces as dual vector-valued 2-cochains. The macroscopic constitutive relation is enforced pointwise at primal 0-cells via musical isomorphisms that map cochains to smooth fields (and vice versa), while linear momentum balance is enforced exactly at the same 0-cells. The resulting system is recast as a Laplace equation whose material Hodge star is non-local and non-diagonal. Numerical simulations on several classical problems with known analytic solutions are reported to show good agreement, and the method is claimed to furnish the exact discrete force-displacement relations for any prescribed macroscopic elasticity tensor on an arbitrary lattice.

Significance. If the constitutive mapping via musical isomorphisms can be shown to enforce the macroscopic tensor exactly at the discrete level without uncontrolled interpolation error on general cell complexes, the approach would supply a geometrically intrinsic, parameter-free route to lattice elasticity that respects DEC duality and extends naturally to dissipative models. The numerical examples on standard test problems constitute a positive but preliminary indication of practical utility.

major comments (2)
  1. [Abstract / constitutive enforcement paragraph] The central claim that the formulation exactly reproduces any prescribed macroscopic elasticity on arbitrary lattices rests on the constitutive step at primal 0-cells. The abstract states that musical isomorphisms are used to map discrete cochains to smooth fields so that the macroscopic tensor can be imposed, yet no explicit construction of the induced non-local material Hodge star, no proof that the pointwise mapping preserves the exact linear relation between 1-cochains and 2-cochains, and no error analysis for irregular complexes are supplied. This directly affects the load-bearing assertion that the discrete relations are parameter-free and exact.
  2. [Numerical simulations section] Numerical validation is reported only as “good agreement with known solutions.” No error norms, convergence rates under mesh refinement, comparison against standard finite-element or finite-volume schemes, or details on how the non-local Hodge star is assembled and inverted are provided. Without these quantitative diagnostics the soundness of the discrete constitutive map cannot be assessed.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and outline planned revisions to improve clarity and completeness.

read point-by-point responses

  1. Referee: [Abstract / constitutive enforcement paragraph] The central claim that the formulation exactly reproduces any prescribed macroscopic elasticity on arbitrary lattices rests on the constitutive step at primal 0-cells. The abstract states that musical isomorphisms are used to map discrete cochains to smooth fields so that the macroscopic tensor can be imposed, yet no explicit construction of the induced non-local material Hodge star, no proof that the pointwise mapping preserves the exact linear relation between 1-cochains and 2-cochains, and no error analysis for irregular complexes are supplied. This directly affects the load-bearing assertion that the discrete relations are parameter-free and exact.

    Authors: We agree that the manuscript would benefit from greater explicitness on this point. The constitutive mapping is performed pointwise at each primal 0-cell using the musical isomorphisms to pull back the macroscopic tensor, which induces the non-local material Hodge star appearing in the discrete Laplace equation. However, we did not supply a separate derivation of the resulting operator or a formal proof that the linear relation between 1-cochains and 2-cochains is preserved exactly for arbitrary cell complexes. In the revised version we will add an appendix containing the explicit construction of the material Hodge star together with a short discussion of the conditions under which the mapping remains exact. We will also include a brief remark on possible interpolation effects on highly irregular complexes and note this as an item for future analysis. These additions will clarify rather than alter the central claim. revision: partial

  2. Referee: [Numerical simulations section] Numerical validation is reported only as “good agreement with known solutions.” No error norms, convergence rates under mesh refinement, comparison against standard finite-element or finite-volume schemes, or details on how the non-local Hodge star is assembled and inverted are provided. Without these quantitative diagnostics the soundness of the discrete constitutive map cannot be assessed.

    Authors: We accept that the numerical section would be strengthened by quantitative measures. The present examples were chosen to illustrate qualitative agreement on standard analytic problems; quantitative diagnostics were omitted to keep the focus on the geometric formulation. In the revision we will report L2 error norms for the displacement and stress fields, include convergence studies under successive refinement where the test problems permit, and provide a short description of the assembly and sparse direct solution of the non-local system. A brief comparison with a standard linear finite-element discretization on the same meshes will also be added. revision: yes

Circularity Check

0 steps flagged

Derivation self-contained from DEC principles with external validation

full rationale

The paper derives the discrete formulation directly from discrete exterior calculus: primal vector-valued 0-cochains for displacements, 1-cochains for differences, dual 2-cochains for forces; constitutive law injected at primal 0-cells via musical isomorphisms; momentum balance enforced exactly at 0-cells; system solved as Laplace equation with non-local material Hodge star. Numerical results are compared to independent analytic solutions for classical problems, providing external falsifiability. No quoted step reduces a claimed prediction or result to a fitted parameter, self-citation chain, or definitional equivalence; the central claim (relations for arbitrary lattices matching prescribed macroscopic elasticity) follows from the DEC setup rather than being presupposed by it.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The formulation rests on the existence of musical isomorphisms that faithfully transfer the continuum constitutive law to the discrete setting and on the assumption that balance of linear momentum can be localized at primal 0-cells for any cell complex.

axioms (2)
  • domain assumption Musical isomorphisms map cochains to smooth fields such that the macroscopic constitutive relation holds at primal 0-cells.
    Invoked to enforce the stress-strain law on the discrete structure.
  • domain assumption Balance of linear momentum is established directly at primal 0-cells.
    Used to close the governing equations.

pith-pipeline@v0.9.0 · 5706 in / 1392 out tokens · 29966 ms · 2026-05-24T12:43:46.309006+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages · 1 internal anchor

  1. [1]

    Abraham, J.E

    R. Abraham, J.E. Marsden, R. Ratiu, Manifolds, tensor analysis, and applications: 2nd Edition, Springer-Verlag, Berlin, Heidelberg, 1988

    work page 1988

  2. [2]

    Angoshtari, A

    A. Angoshtari, A. Yavari, A geometric structure-preserving discretization scheme for incompressible linearized elasticity, Computer Methods in Applied Mechanics and Engineering 259 (2013) 130 – 153. doi:https://doi.org/10.1016/j.cma.2013.03.004

    work page doi:10.1016/j.cma.2013.03.004 2013

  3. [3]

    Bell, A.N

    N. Bell, A.N. Hirani, Pydec: software and algorithms for discretization of exterior calculus, ACM Trans. Math. Softw. 39 (2012). doi: 10.1145/2382585.2382588

    work page doi:10.1145/2382585.2382588 2012

  4. [4]

    Chen, W.C

    S. Chen, W.C. Chew, Discrete electromagnetic theory with exterior calculus, in: 2016 Progress in Electromagnetic Research Symposium (PIERS), pp. 896–897

    work page 2016

  5. [5]

    Discrete Exterior Calculus

    M. Desbrun, A.N. Hirani, M. Leok, J.E. Marsden, Discrete exterior calculus, 2005. arXiv:math/0508341

    work page internal anchor Pith review Pith/arXiv arXiv 2005

  6. [6]

    Gillette, Notes on discrete exterior calculus, Technical Report, University of Texas at Austin, 2009

    A. Gillette, Notes on discrete exterior calculus, Technical Report, University of Texas at Austin, 2009

    work page 2009

  7. [7]

    Gillette, C

    A. Gillette, C. Bajaj, Dual formulations of mixed finite element methods with applications, Computer- Aided Design 43 (2011) 1213 – 1221. doi: https://doi.org/10.1016/j.cad.2011.06.017

    work page doi:10.1016/j.cad.2011.06.017 2011

  8. [8]

    Grady, J.R

    L.J. Grady, J.R. Polimeni, Discrete calculus: applied analysis on graphs for computational science, 1st ed., Springer Publishing Company, Incorporated, 2010

    work page 2010

  9. [9]

    Hirani, Discrete exterior calculus, Ph.D

    A.N. Hirani, Discrete exterior calculus, Ph.D. thesis, USA, 2003. AAI3086864. 26

    work page 2003

  10. [10]

    Hirani, K

    A.N. Hirani, K.B. Nakshatrala, J.H. Chaudhry, Numerical method for darcy flow derived using dis- crete exterior calculus, International Journal for Computational Methods in Engineering Science and Mechanics 16 (2015) 151–169. doi: 10.1080/15502287.2014.977500

    work page doi:10.1080/15502287.2014.977500 2015

  11. [11]

    Mohamed, A.N

    M.S. Mohamed, A.N. Hirani, R. Samtaney, Discrete exterior calculus discretization of incompressible navier–stokes equations over surface simplicial meshes, Journal of Computational Physics 312 (2016) 175–191. doi:10.1016/j.jcp.2016.02.028

    work page doi:10.1016/j.jcp.2016.02.028 2016

  12. [12]

    Phillips, Crystals, defects and microstructures: modelling across scales, 1st ed., Oxford University Press, Cambridge, 2009

    R.B. Phillips, Crystals, defects and microstructures: modelling across scales, 1st ed., Oxford University Press, Cambridge, 2009

    work page 2009

  13. [13]

    Phlipot, K.D

    G.P. Phlipot, K.D. M., A quasicontinuum theory for the nonlinear mechanical response of general periodic truss lattices, Journal of the Mechanics and Physics of Solids 124 (2019) 758–780. doi: 10. 1016/j.jmps.2018.11.014

    work page 2019

  14. [14]

    Si, Tetgen, a delaunay-based quality tetrahedral mesh generator, ACM Transactions on Mathematical Software 41 (2) (2015) 11

    H. Si, Tetgen, a delaunay-based quality tetrahedral mesh generator, ACM Trans. Math. Softw. 41 (2015). doi:10.1145/2629697

    work page doi:10.1145/2629697 2015

  15. [15]

    Silling, Reformulation of elasticity theory for discontinuities and long-range forces, Journal of the Mechanics and Physics of Solids 48 (2000) 175–209

    S.A. Silling, Reformulation of elasticity theory for discontinuities and long-range forces, Journal of the Mechanics and Physics of Solids 48 (2000) 175–209. doi: 10.1016/S0022-5096(99)00029-0

    work page doi:10.1016/s0022-5096(99)00029-0 2000

  16. [16]

    Steinhauser, Computational multiscale modelling of liquids and solids: theory and applications, 2nd ed., Springer-Verlag, Berlin, Heidelberg, 2017

    M.O. Steinhauser, Computational multiscale modelling of liquids and solids: theory and applications, 2nd ed., Springer-Verlag, Berlin, Heidelberg, 2017

    work page 2017

  17. [17]

    S.P. Timoshenko, On the correction for shear of the differential equation for transverse vibrations of prismatic bars, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 41 (1921) 744–746. doi: 10.1080/14786442108636264

    work page doi:10.1080/14786442108636264 1921

  18. [18]

    Tonti, The mathematical structure of classical and relativistic physics: a general classification dia- gram, 1st ed., Springer Science+Business Media, New York, 2013

    E. Tonti, The mathematical structure of classical and relativistic physics: a general classification dia- gram, 1st ed., Springer Science+Business Media, New York, 2013

    work page 2013

  19. [19]

    Y. Wang, M. Peter, Macroscopic elastic properties of regular lattices, Journal of the Mechanics and Physics of Solids 56 (2008) 3459–3474. doi: 10.1016/j.jmps.2008.08.011

    work page doi:10.1016/j.jmps.2008.08.011 2008

  20. [20]

    L. Ward, W. Chris, Atomistic calculations and materials informatics: a review, Current Opinion in Solid State and Materials Science 21 (2017) 167–176. doi: 10.1016/j.cossms.2016.07.002

    work page doi:10.1016/j.cossms.2016.07.002 2017

  21. [21]

    Whitney, Geometric integration theory, Princeton University Press, 1957

    H. Whitney, Geometric integration theory, Princeton University Press, 1957

    work page 1957

  22. [22]

    Yavari, On geometric discretization of elasticity, Journal of Mathematical Physics 49 (2008) 022901

    A. Yavari, On geometric discretization of elasticity, Journal of Mathematical Physics 49 (2008) 022901. doi:10.1063/1.2830977

    work page doi:10.1063/1.2830977 2008

  23. [23]

    H. Zhu, Z. Zhou, R. Yang, A. Yu, Discrete particle simulation of particulate systems: theoretical developments, Chemical Engineering Science 62 (2007) 3378–3396. doi: 10.1016/j.ces.2006.12.089. 27

    work page doi:10.1016/j.ces.2006.12.089 2007