pith. sign in

arxiv: 2604.10613 · v1 · submitted 2026-04-12 · 🧮 math.NA · cs.NA

A High-Order Conformal FEM for Multidimensional Nonlinear Collisional Breakage Equations: Analysis and Computation

Arushi Arushi , Naresh Kumar This is my paper

Pith reviewed 2026-05-10 15:39 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords nonlinear collisional breakage equationsconformal finite element methodhigh-order Lagrange elementsBDF2 time discretizationmoment conservationmultidimensional population balanceconvergence analysisparticle breakage modeling
0
0 comments X

The pith

Conformal finite element methods with high-order Lagrange elements and BDF2 time stepping solve nonlinear collisional breakage equations in multiple dimensions while preserving total particle count and hypervolume.

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

The authors introduce a numerical framework based on conformal finite elements to solve the nonlinear collisional breakage equations that describe how particles fragment through collisions. The approach discretizes space with high-order polynomial basis functions and advances time with a second-order backward differentiation formula, extending to one-, two-, and three-dimensional problems. If the method works as claimed, it produces solutions that keep the total number of particles and their combined volume unchanged, matching physical conservation laws. Error estimates are derived for both the spatial discretization and the complete time-dependent scheme. Tests in multiple dimensions confirm that the computed solutions stay accurate and respect the conservation properties without extra stabilization.

Core claim

We present the first conformal finite element discretization of the multidimensional nonlinear collisional breakage equation. High-order Lagrange elements on simplicial meshes are combined with the BDF2 time integrator to approximate the nonlinear integral collision terms. The resulting scheme conserves the zeroth and first moments of the particle distribution function, which correspond to total particle number and hypervolume. Rigorous a priori error estimates are established for the semidiscrete problem and the fully discrete scheme, and numerical experiments in one, two, and three dimensions demonstrate optimal convergence rates together with conservation of the physical quantities.

What carries the argument

High-order Lagrange finite elements on conformal simplicial meshes, paired with BDF2 time stepping, that project the nonlinear breakage operator into the discrete space while enforcing moment conservation through the variational formulation.

If this is right

  • Optimal convergence rates hold for both the semidiscrete and fully discrete formulations.
  • The zeroth and first moments remain exactly conserved at the discrete level.
  • The method applies directly to one-, two-, and three-dimensional problems with nonlinear kernels.
  • High accuracy is obtained for the breakage equation without additional stabilization terms.

Where Pith is reading between the lines

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

  • The same variational structure could be reused for related population-balance models that contain integral collision operators.
  • Exact moment preservation supports reliable long-time integration of breakage processes in applications such as granulation or aerosol modeling.
  • The convergence theory supplies a template for analyzing similar high-order discretizations of integro-differential equations arising in kinetic theory.

Load-bearing premise

The nonlinear integral collision kernels can be evaluated accurately inside the high-order finite element space in two and three dimensions without introducing instabilities or destroying conservation.

What would settle it

A three-dimensional simulation in which the computed total particle count deviates from its initial value by more than the expected truncation error after many time steps would falsify the conservation property.

Figures

Figures reproduced from arXiv: 2604.10613 by Arushi Arushi, Naresh Kumar.

Figure 1
Figure 1. Figure 1: Comparison of particle breakage mechanisms [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Comparison of Test Case 1 for different grid resolutions. [PITH_FULL_IMAGE:figures/full_fig_p022_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Comparison of Test Case 2 for different grid resolutions. [PITH_FULL_IMAGE:figures/full_fig_p023_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Comparison of Test Case 3 for different grid resolutions. [PITH_FULL_IMAGE:figures/full_fig_p025_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Comparison of Test Case 4 for different grid resolutions. [PITH_FULL_IMAGE:figures/full_fig_p026_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Comparison of Test Case 5 for different grid resolutions. [PITH_FULL_IMAGE:figures/full_fig_p027_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Comparison of Test Case 6 for different grid resolutions. [PITH_FULL_IMAGE:figures/full_fig_p028_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Comparison of numerical solutions of Test Case 1 and Test Case 2 with the respective [PITH_FULL_IMAGE:figures/full_fig_p031_8.png] view at source ↗
read the original abstract

Particle breakage due to collisional interactions plays a vital role in the development of several phenomena in science and engineering. The nonlinear collisional breakage equations (NCBEs) are a significant set of equations in this context. Solving the NCBE is computationally challenging due to its nonlinearity, high dimensionality, and complex kernel interactions. Solving NCBE problems is more complex in two- and three-dimensional problems. In these problems, it is more challenging to evaluate multidimensional moments and integrals, maintain solution stability, and achieve computational efficiency. Despite the importance of the NCBE in science and engineering, the development of efficient numerical methods for solving it in two- and three-dimensional problems has not been adequately explored. In this work, we have introduced a new framework for solving the NCBE across multiple dimensions using the conformal finite element method (FEM). To the best of our knowledge, this is the first work to solve the NCBE using the conformal FEM. The new framework employs high-order Lagrange elements in conjunction with the BDF2 scheme for time discretization. The present method preserves the important physical quantities such as the total count and hypervolume of the population particles. Convergence results for error estimates have also been derived for both semidiscrete and fully discrete schemes. Numerical experiments have been carried out for one-, two-, and three-dimensional problems. The numerical experiments have shown that the proposed method achieved high accuracy, optimal convergence rates, and computational efficiency.

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

0 major / 3 minor

Summary. The manuscript develops a conformal finite element method using high-order Lagrange elements together with BDF2 time discretization for the numerical solution of nonlinear collisional breakage equations (NCBEs) in one to three dimensions. It supplies a weak formulation, derives a priori error bounds for the semidiscrete and fully discrete schemes under standard Lipschitz and growth assumptions on the kernel, proves exact discrete conservation of the zeroth and first moments via suitable test functions and quadrature, and reports numerical experiments confirming optimal convergence rates together with invariant preservation.

Significance. If the error estimates and conservation properties hold under the stated assumptions, the work supplies a stable high-order framework for a class of nonlinear integro-PDEs that are difficult to discretize in higher dimensions. The exact moment preservation and the extension to 2D/3D without added stabilization are notable strengths, as is the combination of rigorous analysis with computational validation. This could be useful for modeling particle breakage processes in materials science and engineering.

minor comments (3)
  1. The abstract states that convergence results have been derived but does not indicate the norms in which the estimates hold or the expected orders; this information should appear explicitly in the abstract or introduction.
  2. In the numerical experiments, the specific quadrature rules used for the nonlinear collision integrals in 2D and 3D should be described in more detail to allow readers to assess how quadrature error interacts with the theoretical rates.
  3. A short discussion of the computational scaling with dimension and the cost of kernel evaluations would help substantiate the efficiency claims made for multidimensional problems.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and for recommending minor revision. The recognition of the exact discrete conservation properties, the extension to 2D/3D, and the combination of analysis with numerical validation is appreciated. No specific major comments were provided in the report.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The manuscript introduces a conformal FEM discretization of the nonlinear collisional breakage equation using high-order Lagrange elements and BDF2 time stepping. It derives a priori error bounds for the semidiscrete and fully discrete schemes directly from the weak formulation under standard Lipschitz and growth assumptions on the breakage kernel, without invoking fitted parameters or self-referential definitions. Discrete conservation of the zeroth and first moments follows immediately from the choice of test functions and quadrature rules in the weak form. Numerical experiments confirm the predicted convergence rates in 1D–3D. No load-bearing step reduces by construction to an input, self-citation chain, or ansatz smuggled from prior work; the analysis rests on classical FEM theory applied to the given integral equation.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework rests on standard finite-element theory for error estimates and on the stability properties of BDF2; no new free parameters, ad-hoc axioms, or invented entities are introduced in the abstract.

axioms (2)
  • standard math Standard assumptions of conforming finite element spaces and quadrature rules suffice for convergence of the semidiscrete and fully discrete schemes
    Invoked when deriving error estimates for the nonlinear problem.
  • domain assumption BDF2 time discretization preserves the discrete invariants (total count and hypervolume) for the chosen spatial discretization
    Central to the claimed physical fidelity of the method.

pith-pipeline@v0.9.0 · 5564 in / 1380 out tokens · 37011 ms · 2026-05-10T15:39:11.673996+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

36 extracted references · 36 canonical work pages

  1. [1]

    Arora, S

    G. Arora, S. Hussain, and R. Kumar. Comparison of variational iteration and adomian de- composition methods to solve growth, aggregation and aggregation-breakage equations.J. Comput. Sci., 67:101973, 2023

    work page 2023

  2. [2]

    Arora, S

    G. Arora, S. Hussain, G. Walker, C. Cattani, S. Heinrich, and M. Singh. An accurate approach and its convergence analysis for the multidimensional nonlinear collisional breakage equations. Chemical Engineering Science, page 123400, 2026

    work page 2026

  3. [3]

    M. M. Attarakih, C. Drumm, and H.-J. Bart. Solution of the population balance equation using the sectional quadrature method of moments (sqmom).Chem. Eng. Sci., 64(4):742–752, 2009

    work page 2009

  4. [4]

    P. K. Barik and A. K. Giri. Global classical solutions to the continuous coagulation equation with collisional breakage.Z. Angew. Math. Phys., 71(1):38, 2020

    work page 2020

  5. [5]

    S. K. Bariwal, S. Hussain, and R. Kumar. Non-linear collision-induced breakage equation: Finite volume and semi-analytical methods.Acta Appl. Math., 192(1):6, 2024

    work page 2024

  6. [6]

    Brilliantov, P

    N. Brilliantov, P. Krapivsky, A. Bodrova, F. Spahn, H. Hayakawa, V. Stadnichuk, and J. Schmidt. Size distribution of particles in saturn’s rings from aggregation and fragmen- tation.Proc. Natl. Acad. Sci. USA, 112(31):9536–9541, 2015

    work page 2015

  7. [7]

    Cheng and S

    Z. Cheng and S. Redner. Scaling theory of fragmentation.PRL, 60(24):2450, 1988

    work page 1988

  8. [8]

    Cheng and S

    Z. Cheng and S. Redner. Kinetics of fragmentation.J. Phys. A: Math. Gen., 23(7):1233–1258, 1990

    work page 1990

  9. [9]

    A. Das, J. Kumar, M. Dosta, and S. Heinrich. On the approximate solution and modeling of the kernel of nonlinear breakage population balance equation.SIAM J. Sci. Comput., 42(6):B1570–B1598, 2020

    work page 2020

  10. [10]

    Dubovskii, V

    P. Dubovskii, V. Galkin, and I. Stewart. Exact solutions for the coagulation-fragmentation equation.J. Phys. A: Math. Gen., 25(18):4737–4744, 1992

    work page 1992

  11. [11]

    M. H. Ernst and I. Pagonabarraga. The nonlinear fragmentation equation.J. Phys. A: Math. Theor., 40(17):F331–F337, 2007

    work page 2007

  12. [12]

    Fries, S

    L. Fries, S. Antonyuk, S. Heinrich, D. Dopfer, and S. Palzer. Collision dynamics in fluidised bed granulators: A dem-cfd study.Chem. Eng. Sci., 86:108–123, 2013

    work page 2013

  13. [13]

    Fuerstenau, A

    D. Fuerstenau, A. De, and P. Kapur. Linear and nonlinear particle breakage processes in comminution systems.Int. J. Miner. Process., 74:S317–S327, 2004

    work page 2004

  14. [14]

    A. K. Giri and P. Lauren¸ cot. Existence and nonexistence for the collision-induced breakage equation.SIAM J. Math. Anal., 53(4):4605–4636, 2021. 33

    work page 2021

  15. [15]

    Hussain, R

    S. Hussain, R. Kumar, et al. Hybridized iterative scheme for solving non-linear collisional- breakage equation.J. Comput. Sci., page 102731, 2025

    work page 2025

  16. [16]

    Keshav, M

    S. Keshav, M. Singh, S. Singh, G. Walker, and J. Kumar. A new and unified semi-analytical method with a convergence acceleration parameter for linear and nonlinear fragmentation equations.Proc. R. Soc. A., 481(2312), 2025

    work page 2025

  17. [17]

    Kiefer and H

    J. Kiefer and H. Straaten. A model of ion track structure based on classical collision dynamics (radiobiology application).Phys. Med. Biol., 31(11):1201–1209, 1986

    work page 1986

  18. [18]

    Kostoglou and A

    M. Kostoglou and A. Karabelas. A study of the nonlinear breakage equation: analytical and asymptotic solutions.Phys. A: Math. Gen., 33(6):1221–1232, 2000

    work page 2000

  19. [19]

    Kushwah, A

    P. Kushwah, A. Das, J. Saha, and A. B¨ uck. Population balance equation for collisional break- age: a new numerical solution scheme and its convergence.Commun. Nonlinear Sci. Numer. Simul., 121:107244, 2023

    work page 2023

  20. [20]

    Kushwah, A

    P. Kushwah, A. Ghorai, and J. Saha. Convergence-rate and error analysis of sectional-volume average method for the collisional breakage equation with multi-dimensional modelling.arXiv preprint arXiv:2504.15365, 2025

    work page arXiv 2025

  21. [21]

    Lauren¸ cot and D

    P. Lauren¸ cot and D. Wrzosek. The discrete coagulation equations with collisional breakage. J. Stat. Phys., 104(1):193–220, 2001

    work page 2001

  22. [22]

    S. L. Leong, M. Singh, F. Ahamed, S. Heinrich, S. I. X. Tiong, I. M. L. Chew, and Y. K. Ho. A comparative study of the fixed pivot technique and finite volume schemes for multi-dimensional breakage population balances.Adv. Powder Technol., 34(12):104272, 2023

    work page 2023

  23. [23]

    Lombart, C.-E

    M. Lombart, C.-E. Br´ ehier, M. Hutchison, and Y.-N. Lee. General non-linear fragmentation with discontinuous galerkin methods.MNRAS., 533(4):4410–4434, 2024

    work page 2024

  24. [24]

    Lombart, M

    M. Lombart, M. Hutchison, and Y.-N. Lee. Fragmentation with discontinuous galerkin schemes: non-linear fragmentation.MNRAS, 517(2):2012–2027, 2022

    work page 2012

  25. [25]

    G. M. McFarquhar. A new representation of collision-induced breakup of raindrops and its implications for the shapes of raindrop size distributions.J. Atmos. Sci., 61(7):777–794, 2004

    work page 2004

  26. [26]

    J. Paul, A. Das, and J. Kumar. Moments preserving finite volume approximations for the non-linear collisional fragmentation model.Appl. Math. Comput.., 436:127494, 2023

    work page 2023

  27. [27]

    Paul and J

    J. Paul and J. Kumar. An existence-uniqueness result for the pure binary collisional breakage equation.Math. Meth. Appl. Sci., 41(7):2715–2732, 2018

    work page 2018

  28. [28]

    Piotrowski

    S. Piotrowski. The collisions of asteroids.Acta Astron., 6:115–138, 1953

    work page 1953

  29. [29]

    Planchette, H

    C. Planchette, H. Hinterbichler, M. Liu, D. Bothe, and G. Brenn. Colliding drops as coalescing and fragmenting liquid springs.J. Fluid Mech., 814:277–300, 2017. 34

    work page 2017

  30. [30]

    Saha and M

    J. Saha and M. Singh. Rate of convergence and stability analysis of a modified fixed pivot technique for a fragmentation equation.Numer. Math., 153(2):531–555, 2023

    work page 2023

  31. [31]

    Thom´ ee.Galerkin finite element methods for parabolic problems, volume 25

    V. Thom´ ee.Galerkin finite element methods for parabolic problems, volume 25. Springer Science & Business Media, 2007

    work page 2007

  32. [32]

    Tournus, M

    M. Tournus, M. Escobedo, W.-F. Xue, and M. Doumic. Insights into the dynamic trajectories of protein filament division revealed by numerical investigation into the mathematical model of pure fragmentation.PLoS Comput. Biol., 17(9):e1008964, 2021

    work page 2021

  33. [33]

    R. D. Vigil, I. Vermeersch, and R. O. Fox. Destructive aggregation: Aggregation with collision- induced breakage.J. Colloid Interface Sci., 302(1):149–158, 2006

    work page 2006

  34. [34]

    L. G. Wang, R. Ge, X. Chen, R. Zhou, and H.-M. Chen. Multiscale digital twin for particle breakage in milling: From nanoindentation to population balance model.Powder Technol., 386:247–261, 2021

    work page 2021

  35. [35]

    Yadav, A

    N. Yadav, A. Das, M. Singh, S. Singh, and J. Kumar. Homotopy perturbation method and its convergence analysis for nonlinear collisional fragmentation equations.Proc. R. Soc. A, 479(2279):20230567, 2023

    work page 2023

  36. [36]

    R. M. Ziff and E. D. McGrady. The kinetics of cluster fragmentation and depolymerisation. J. Phys. A: Math. Gen., 18(15):3027–3037, 1985. 35

    work page 1985