pith. machine review for the scientific record. sign in

arxiv: 2604.27407 · v2 · submitted 2026-04-30 · 🧮 math.NA · cs.NA· math.AP

Recognition: unknown

A Shifted Cohesive-Zone Method for Non-Interface-Fitted Meshes with Applications to Crystal Plasticity

Cheng-Hau Yang , Mark C. Messner , Tianchen Hu
Authors on Pith no claims yet

Pith reviewed 2026-05-07 10:28 UTC · model grok-4.3

classification 🧮 math.NA cs.NAmath.AP
keywords cohesive zone methodshifted boundary methodcrystal plasticitynon-interface-fitted meshesfinite element methodinterface mechanicssurrogate interfacedamage evolution
0
0 comments X

The pith

The shifted cohesive zone method simulates interface mechanics accurately on non-interface-fitted meshes.

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

This paper develops the Shifted Cohesive Zone Method to enforce traction-separation laws on a surrogate interface rather than the true material interface. The shift permits standard finite-element spaces on meshes that do not conform to the interfaces, removing the need to generate fitted quadrilateral or hexahedral meshes around complex microstructures. Verification shows first-order convergence together with close agreement to reference solutions for reaction forces, energy release, deformation, stress, and damage evolution under crystal plasticity. The approach matters because interface-conforming meshes remain difficult or impossible to construct for many realistic grain geometries. A geometry-aware classification step accelerates construction of the surrogate domains.

Core claim

By shifting the enforcement of traction-separation laws from the true interface to a nearby surrogate interface, SCZM enables the use of standard finite element spaces while avoiding the meshing burden associated with interface-conformal discretizations, yielding first-order convergence and close quantitative agreement with interface-fitted reference solutions in reaction forces, surface energy release, deformation, stress fields, and damage evolution for both linear elasticity and history-dependent crystal plasticity.

What carries the argument

The surrogate interface to which the traction-separation laws are shifted, which carries the weak-form enforcement and produces the nonlinear residual and consistent tangent.

If this is right

  • Standard finite-element spaces suffice for interface problems once the traction-separation laws are applied on the surrogate surface.
  • First-order convergence is recovered for non-interface-fitted discretizations of cohesive-zone problems.
  • Reaction forces, energy release rates, deformation, stress, and damage histories remain close to those obtained with conforming meshes.
  • The same framework accommodates multiple traction-separation laws together with history-dependent crystal plasticity.

Where Pith is reading between the lines

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

  • The surrogate-interface idea could be combined with local mesh refinement near the true interface to recover higher-order accuracy without global conformity.
  • The same shift technique may apply to other interface conditions, such as frictional contact or thermal resistance, provided the constitutive response is not overly sensitive to small geometric offsets.
  • Automated surrogate construction via point classification suggests a route toward fully automated simulation pipelines for polycrystalline aggregates.

Load-bearing premise

The error introduced by moving the traction-separation laws from the true interface to a nearby surrogate interface remains negligible for nonlinear history-dependent constitutive models such as crystal plasticity.

What would settle it

A benchmark in which SCZM on a non-fitted mesh produces stress or damage fields that deviate substantially from a highly refined interface-fitted reference solution under the same crystal-plasticity loading would falsify the accuracy claim.

read the original abstract

The accurate simulation of interface-dominated solid mechanics problems on complex microstructures remains challenging, particularly when interface-fitted quadrilateral or hexahedral meshes are difficult to generate. We extend the shifted boundary method (SBM) to cohesive-zone formulations and introduce the Shifted Cohesive Zone Method (SCZM), with applications to crystal plasticity on non-interface-fitted meshes. By shifting the enforcement of traction-separation laws from the true interface to a nearby surrogate interface, SCZM enables the use of standard finite element spaces while avoiding the meshing burden associated with interface-conformal discretizations. We present a simplified SCZM weak form defined on the surrogate interface, leading to a straightforward implementation of the nonlinear residual and consistent tangent matrix. The method is implemented in the open-source MOOSE framework and coupled with constitutive models from NEML2, enabling simulations with linear elasticity, multiple traction-separation laws, and history-dependent crystal plasticity. We further develop a geometry-aware, PCA-enhanced point classification algorithm to accelerate surrogate-domain construction. Verification and benchmark studies in two and three dimensions demonstrate that SCZM achieves first-order convergence for non-interface-fitted interface problems and closely matches interface-fitted reference solutions in terms of reaction forces, surface energy release, deformation, stress fields, and damage evolution. These results indicate that SCZM provides an accurate and efficient framework for modeling interface mechanics in complex microstructures without requiring interface-fitted meshes.

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 / 2 minor

Summary. The paper introduces the Shifted Cohesive Zone Method (SCZM) as an extension of the shifted boundary method to cohesive-zone formulations. It shifts traction-separation law enforcement to a nearby surrogate interface on non-interface-fitted meshes, presents a simplified weak form, develops a PCA-enhanced point classification algorithm for surrogate-domain construction, implements the method in MOOSE coupled with NEML2 for linear elasticity and history-dependent crystal plasticity, and reports verification studies showing first-order convergence together with close agreement to interface-fitted references in reaction forces, surface energy release, deformation, stress fields, and damage evolution.

Significance. If the numerical results hold, SCZM offers a practical route to interface-mechanics simulations on complex microstructures without the burden of interface-conformal meshing. The open-source implementation, coupling to nonlinear constitutive models, and direct numerical comparison to fitted-mesh references for both linear and crystal-plasticity cases are concrete strengths that could facilitate adoption in solid-mechanics modeling.

major comments (1)
  1. [Verification and benchmark studies] Verification and benchmark studies: the central claim that SCZM achieves first-order convergence for non-interface-fitted problems rests on the reported numerical results; the manuscript should explicitly state the error norms employed (e.g., L2 displacement or traction error) and demonstrate that the observed order remains first-order when the shift distance is comparable to the mesh size h for the history-dependent crystal-plasticity model with evolving damage.
minor comments (2)
  1. [Weak-form derivation] The derivation of the consistent tangent matrix from the simplified weak form on the surrogate interface would benefit from one additional intermediate step showing how the jump operator and traction-separation law are linearized.
  2. [Numerical results] Figure captions comparing SCZM and reference solutions should explicitly note the mesh size and shift distance used so that readers can assess the visual agreement quantitatively.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment of our work and the constructive comment on the verification studies. We address the point below and will incorporate the requested clarifications in a revised manuscript.

read point-by-point responses

  1. Referee: [Verification and benchmark studies] Verification and benchmark studies: the central claim that SCZM achieves first-order convergence for non-interface-fitted problems rests on the reported numerical results; the manuscript should explicitly state the error norms employed (e.g., L2 displacement or traction error) and demonstrate that the observed order remains first-order when the shift distance is comparable to the mesh size h for the history-dependent crystal-plasticity model with evolving damage.

    Authors: We agree that greater explicitness on the error norms and on the shift-distance regime for the crystal-plasticity cases will strengthen the presentation. In the revised manuscript we will state that all convergence studies employ the L2 norm of the displacement field (integrated over the domain) together with the L2 norm of the traction error evaluated on the surrogate interface. Our existing two- and three-dimensional benchmarks already include crystal-plasticity simulations with evolving damage in which the surrogate-interface shift distance is O(h); the observed rates remain first-order. To make this transparent we will add a short dedicated paragraph and a supplementary convergence plot that isolates the crystal-plasticity data, explicitly reports the shift-to-h ratios employed, and confirms that the first-order behavior is retained. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper derives the SCZM weak form by direct extension of the standard shifted boundary method to cohesive-zone traction-separation laws, placing the enforcement on a surrogate interface constructed via PCA classification. This yields a standard nonlinear residual and tangent without any fitted parameter being relabeled as a prediction, without self-citation chains supporting the core variational statement, and without any uniqueness theorem imported from the authors' prior work. All central claims (first-order convergence, agreement with interface-fitted references) rest on explicit numerical benchmarks rather than on definitional equivalence or ansatz smuggling. The method is therefore independent of its own outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

Ledger constructed from abstract only; full text unavailable for exhaustive audit.

axioms (1)
  • standard math Standard finite element weak forms and discretization remain valid when boundary conditions are shifted to a surrogate interface.
    The method relies on this to define the simplified weak form.
invented entities (1)
  • Shifted Cohesive Zone Method (SCZM) no independent evidence
    purpose: Enables cohesive zone modeling on non-interface-fitted meshes
    Newly proposed method extending SBM.

pith-pipeline@v0.9.0 · 5565 in / 1174 out tokens · 92546 ms · 2026-05-07T10:28:04.435617+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

53 extracted references · 4 canonical work pages

  1. [1]

    G. N. Wells, L. Sluys, A new method for modelling cohesive cracks using finite elements, International Journal for numerical methods in engineering 50 (2001) 2667–2682

    2001

  2. [2]

    L. Lee, R. J. LeVeque, An immersed interface method for incompress- ible navier–stokes equations, SIAM Journal on Scientific Computing 25 (2003) 832–856

    2003

  3. [3]

    Bänsch, F

    E. Bänsch, F. Haußer, O. Lakkis, B. Li, A. Voigt, Finite element method for epitaxial growth with attachment–detachment kinetics, Journal of Computational Physics 194 (2004) 409–434

    2004

  4. [4]

    Y. Xiao, J. Xu, F. Wang, High-order extended finite element methods for solving interface problems, Computer Methods in Applied Mechanics and Engineering 364 (2020) 112964

    2020

  5. [5]

    Z. Zhao, J. Yan, Enriched immersed boundary method (eibm) for interface-coupled multi-physics and applications to convective conjugate heat transfer, Computer Methods in Applied Mechanics and Engineer- ing 401 (2022) 115667

    2022

  6. [6]

    T. Chen, B. Barua, T. Hu, M. C. Messner, T. El-Wardany, An ICME Modeling Framework for Titanium/Tungsten-Carbide Metal Matrix Composites, Technical Report, Argonne National Laboratory (ANL), Argonne, IL (United States), 2023

    2023

  7. [7]

    T. Hu, H. H. J. Choi, M. C. Messner, A mechanistic model for creep and thermal aging in Alloy 709, Technical Report, Argonne National Laboratory (ANL), Argonne, IL (United States), 2023. 42

    2023

  8. [8]

    Baweja, T

    S. Baweja, T. Hu, M. C. Messner, Predicting long-term stress relaxation on Alloy 709 using the crystal plasticity finite element method, Techni- cal Report, Argonne National Laboratory (ANL), Argonne, IL (United States), 2025

    2025

  9. [9]

    Baweja, M

    S. Baweja, M. C. Messner, Development of predictive model for accu- rate rupture time from multi-axial creep in alloy 709 with physics-based simulations, Computational Materials Science 263 (2026) 114427

    2026

  10. [10]

    Aduloju, T

    S. Aduloju, T. J. Truster, Modeling of stress and temperature effects on creep of reduced activation ferritic-martensitic steel alloy f82h (2025)

    2025

  11. [11]

    Bhatt, M

    S. Bhatt, M. C. Messner, R. Wang, F. Gao, Microstructural models for the creep strength and ductility of diffusion-bonded 316h steel, in: Advances in Materials Technology for Power Plants, volume 84871, ASM International, 2024, pp. 750–759

    2024

  12. [12]

    Nassif, T

    O. Nassif, T. J. Truster, R. Ma, K. B. Cochran, D. M. Parks, M. C. Messner, T. Sham, Combined crystal plasticity and grain boundary modeling of creep in ferritic-martensitic steels: I. theory and implemen- tation, Modelling and Simulation in Materials Science and Engineering 27 (2019) 075009

    2019

  13. [13]

    Z. Li, T. Lin, X. Wu, New cartesian grid methods for interface problems using the finite element formulation, Numerische Mathematik 96 (2003) 61–98

    2003

  14. [14]

    Adjerid, I

    S. Adjerid, I. Babuška, R. Guo, T. Lin, An enriched immersed finite element method for interface problems with nonhomogeneous jump con- ditions, Computer Methods in Applied Mechanics and Engineering 404 (2023) 115770. 43

    2023

  15. [15]

    N. Moës, J. Dolbow, T. Belytschko, A finite element method for crack growth without remeshing, International journal for numerical methods in engineering 46 (1999) 131–150

    1999

  16. [16]

    K. Li, N. M. Atallah, A. Rodríguez-Ferran, D. M. Valiveti, G. Scov- azzi, The shifted fracture method, International Journal for Numerical Methods in Engineering 122 (2021) 6641–6679

    2021

  17. [17]

    K. Li, A. Rodríguez-Ferran, G. Scovazzi, The simple shifted fracture method, International Journal for Numerical Methods in Engineering 124 (2023) 2837–2875

    2023

  18. [18]

    K. Li, A. Rodríguez-Ferran, G. Scovazzi, Crack branching and merging simulations with the shifted fracture method, Computer Methods in Applied Mechanics and Engineering 433 (2025) 117528

    2025

  19. [19]

    A. Main, G. Scovazzi, The shifted boundary method for embedded domain computations. part i: Poisson and stokes problems, Journal of Computational Physics 372 (2018) 972–995

    2018

  20. [20]

    N. M. Atallah, G. Scovazzi, Nonlinear elasticity with the shifted bound- ary method, Computer Methods in Applied Mechanics and Engineering 426 (2024) 116988

    2024

  21. [21]

    X. Zeng, T. Song, G. Scovazzi, A shifted boundary method for the compressible euler equations, Journal of Computational Physics 520 (2025) 113512

    2025

  22. [22]

    A. Main, G. Scovazzi, The shifted boundary method for embedded do- main computations. part ii: Linear advection-diffusion and incompress- ible navier-stokes equations, J. Comput. Phys. 372 (2018) 996–1026. 44

    2018

  23. [24]

    Atallah, C

    N. Atallah, C. Canuto, G. Scovazzi, The shifted boundary method for solid mechanics, International Journal for Numerical Methods in Engi- neering 122 (2021) 5935–5970

    2021

  24. [25]

    K. Li, A. Rodríguez-Ferran, G. Scovazzi, A blended shifted- fracture/phase-field framework for sharp/diffuse crack modeling, In- ternational Journal for Numerical Methods in Engineering 124 (2023) 998–1030

    2023

  25. [26]

    K. Li, A. Gorgi, R. Rossi, G. Scovazzi, The shifted boundary method for contact problems, Computer Methods in Applied Mechanics and Engineering 440 (2025) 117940

    2025

  26. [27]

    K. Li, N. M. Atallah, G. A. Main, G. Scovazzi, The shifted interface method: a flexible approach to embedded interface computations, In- ternational Journal for Numerical Methods in Engineering 121 (2020) 492–518

    2020

  27. [28]

    Colomés, A

    O. Colomés, A. Main, L. Nouveau, G. Scovazzi, A weighted shifted boundary method for free surface flow problems, Journal of Computa- tional Physics 424 (2021) 109837. 45

    2021

  28. [29]

    C.-H. Yang, G. Scovazzi, A. Krishnamurthy, B. Ganapathysubra- manian, A shifted boundary method for thermal flows, Journal of Computational Physics (2025) 114333. URL:https://www. sciencedirect.com/science/article/pii/S0021999125006151. doi:https://doi.org/10.1016/j.jcp.2025.114333

    work page doi:10.1016/j.jcp.2025.114333 2025

  29. [30]

    D. Xu, O. Colomés, A. Main, K. Li, N. M. Atallah, N. Abboud, G. Scov- azzi, A weighted shifted boundary method for immersed moving bound- ary simulations of stokes’ flow, Journal of Computational Physics 510 (2024) 113095

    2024

  30. [31]

    D. Xu, O. Colomés, A. Main, K. Li, N. M. Atallah, N. Abboud, G. Sco- vazzi, A weighted shifted boundary method for the navier-stokes equa- tions with immersed moving boundaries, Journal of Computational Physics (2025) 114571

    2025

  31. [32]

    C.-H. Yang, K. Saurabh, G. Scovazzi, C. Canuto, A. Krishnamurthy, B. Ganapathysubramanian, Optimal surrogate boundary selection and scalability studies for the shifted boundary method on octree meshes, Computer Methods in Applied Mechanics and Engineering 419 (2024) 116686

    2024

  32. [33]

    J. H. Collins, K. Li, A. Lozinski, G. Scovazzi, Gap-sbm: A new concep- tualization of the shifted boundary method with optimal convergence for the neumann and dirichlet problems, arXiv preprint arXiv:2508.09613 (2025)

    work page arXiv 2025

  33. [34]

    Antonelli, A

    N. Antonelli, A. Gorgi, R. Zorrilla, R. Rossi, Isogeometric multipatch coupling with arbitrary refinement and parametrization using the gap– 46 shifted boundary method, Computer Methods in Applied Mechanics and Engineering 456 (2026) 118913

    2026

  34. [35]

    Colomés, J

    O. Colomés, J. Modderman, G. Scovazzi, The generalized shifted bound- ary method for geometry-parametric pdes and time-dependent domains, Computer Methods in Applied Mechanics and Engineering 452 (2026) 118748

    2026

  35. [36]

    C.-H. Yang, G. Scovazzi, A. Krishnamurthy, B. Ganapathysubrama- nian, Octree-based adaptive mesh refinement and the shifted boundary method for efficient fluid dynamics simulations, Advances in Computa- tional Science and Engineering 4 (2025) 57–84

    2025

  36. [37]

    Matthew Weinberg

    S. Karki, M. Shadkhah, C.-H. Yang, A. Balu, G. Scovazzi, A. Krishnamurthy, B. Ganapathysubramanian, Direct flow sim- ulations with implicit neural representation of complex geome- try, Computer Methods in Applied Mechanics and Engineering 446 (2025) 118248. URL:https://www.sciencedirect.com/science/ article/pii/S0045782525005201. doi:https://doi.org/10.10...

    work page doi:10.1016/j 2025

  37. [38]

    Matouš, A

    K. Matouš, A. M. Maniatty, Finite element formulation for mod- elling large deformations in elasto-viscoplastic polycrystals, Interna- tional Journal for Numerical Methods in Engineering 60 (2004) 2313– 2333

    2004

  38. [39]

    Dohrmann, M

    C. Dohrmann, M. Heinstein, J. Jung, S. Key, W. Witkowski, Node- based uniform strain elements for three-node triangular and four-node tetrahedral meshes, International Journal for Numerical Methods in Engineering 47 (2000) 1549–1568. 47

    2000

  39. [40]

    E. D. S. Neto, F. A. Pires, D. Owen, F-bar-based linear triangles and tetrahedra for finite strain analysis of nearly incompressible solids. part i: formulation and benchmarking, International Journal for Numerical Methods in Engineering 62 (2005) 353–383

    2005

  40. [41]

    E. A. de Souza Neto, D. Peric, D. R. Owen, Computational methods for plasticity: theory and applications, John Wiley & Sons, 2008

    2008

  41. [42]

    Cheng, A

    J. Cheng, A. Shahba, S. Ghosh, Stabilized tetrahedral elements for crystal plasticity finite element analysis overcoming volumetric locking, Computational Mechanics 57 (2016) 733–753

    2016

  42. [43]

    A. F. Bower, Applied Mechanics of Solids, CRC press, 2009

    2009

  43. [44]

    T. J. R. Hughes, The Finite Element Method, Linear Static and Dy- namic Finite Element Analysis, New Jersey:Prentice-Hall, 1987

    1987

  44. [45]

    K. B. Nakshatrala, A. Masud, K. D. Hjelmstad, On finite element for- mulations for nearly incompressible linear elasticity, Computational Me- chanics 41 (2008) 547–561

    2008

  45. [46]

    Needleman, A continuum model for void nucleation by inclusion debonding (1987)

    A. Needleman, A continuum model for void nucleation by inclusion debonding (1987)

    1987

  46. [47]

    Needleman, An analysis of decohesion along an imperfect interface, International Journal of Fracture 42 (1990) 21–40

    A. Needleman, An analysis of decohesion along an imperfect interface, International Journal of Fracture 42 (1990) 21–40

    1990

  47. [48]

    X.-P. Xu, A. Needleman, Numerical simulations of fast crack growth in brittle solids, Journal of the Mechanics and Physics of Solids 42 (1994) 1397–1434

    1994

  48. [49]

    P. P. Camanho, C. G. Dávila, Mixed-Mode Decohesion Finite Elements for the Simulation of Delamination in Composite Materials, Technical 48 Report NASA/TM-2002-211737, National Aeronautics and Space Ad- ministration, 2002

    2002

  49. [50]

    M. K. Salehani, N. Irani, A coupled mixed-mode cohesive zone model: An extension to three-dimensional contact problems, arXiv preprint arXiv:1801.03430 (2018)

    work page arXiv 2018

  50. [51]

    Borazjani, L

    I. Borazjani, L. Ge, F. Sotiropoulos, Curvilinear immersed boundary method for simulating fluid structure interaction with complex 3d rigid bodies, Journal of Computational physics 227 (2008) 7587–7620

    2008

  51. [52]

    Haines, Point in polygon strategies., Graphics Gems 4 (1994) 24–46

    E. Haines, Point in polygon strategies., Graphics Gems 4 (1994) 24–46

    1994

  52. [53]

    M. C. Messner, T. Hu, Fully implicit crystal plasticity models repre- senting orientations with modified rodrigues parameters, Mechanics of Materials 208 (2025) 105388

    2025

  53. [54]

    J. L. Blanco, P. K. Rai, nanoflann: a C++ header-only fork of FLANN, a library for nearest neighbor (NN) with kd-trees,https://github. com/jlblancoc/nanoflann, 2014. 49 Appendix A. Supplementary implementation details Appendix A.1. Detailed algorithms for morphology-aware point classification Thissubsectionpresentsthecompletepseudo-codeunderlyingthemorpho...

    2014