[Paper Review] A Conservative Log-Size Master Equation for Fragmentation PBEs: Jump Transport, Drift--Diffusion Asymptotics, and PSD Inference

Juan J. Segura|arXiv (Cornell University)|Jan 10, 2026

Coagulation and Flocculation Studies0 citations

TL;DR

The paper derives an exact conservative master equation for fragmentation PBEs in log-size space, shows how a small-jump limit yields drift–diffusion, presents an optional GKSL/Lindblad factorization for inverse modeling, and validates forward and inverse routes numerically.

ABSTRACT

Fragmentation population-balance equations (PBEs) describe how particle size distributions (PSDs) evolve under breakage and daughter fragment redistribution. From a standard self-similar fragmentation class we derive an \emph{exact conservative transport equation in log-size} for the \emph{normalized mass fraction}: a state-dependent \emph{pure-jump} master equation (nonlocal internal-coordinate mass transfer). We also give an explicit Gorini--Kossakowski--Sudarshan--Lindblad (GKSL) factorization whose diagonal sector reproduces this master equation, used here as an \emph{optional} structure-preserving operator representation and constrained parameterization for inverse modeling (rather than a computational necessity). In a controlled small-jump regime, the nonlocal jump transport reduces to a drift--diffusion (Fokker--Planck) operator in log-size space. Under detailed-balance conditions this operator admits the standard symmetrization to a self-adjoint Schrödinger-type spectral problem, enabling compact parametric hypothesis classes for PSD shapes. We then present two inverse routes: (i) time-resolved parametric fitting of transport/spectral parameters, and (ii) a regularized steady-state inversion that reconstructs an effective potential from a measured steady PSD. To address practical validation, we include numerical benchmarks: forward simulation of the jump transport model (CTMC discretization) and its drift--diffusion reduction, quantitative discrepancy metrics, and inverse parameter recovery on an Airy half-line synthetic benchmark under controlled multiplicative noise.

Motivation & Objective

• Motivate fragmentation PBEs and the need for a conservative, log-space transport formulation.
• Derive an exact conservative master equation for normalized mass fraction in log-size space.
• Connect jump transport to drift–diffusion limits and explore spectral/PSP parametrizations.
• Provide inverse modeling routes for time-resolved fitting and steady-state PSD reconstruction.
• Validate the approach with forward simulations and inverse recovery on synthetic benchmarks.

Proposed method

• Formulate the fragmentation PBE in log-size space and derive the exact one-sided jump transport master equation (Eq. 8).
• Define the normalized mass-fraction density p(ξ,t) and the jump kernel K(u) with the rate λ(ξ) (Eq. 9–10).
• Present an optional GKSL/Lindblad factorization that reproduces the diagonal of the master equation for constrained inverse modeling (Eq. 12–13).
• In the small-jump regime, obtain a drift–diffusion (Fokker–Planck) reduction in log-size (Eq. 16).
• Discuss a detailed-balance-based spectral form and Schrödinger-type reduction for PSD parametrizations (Eq. 17–20).
• Outline two inverse routes: time-resolved parametric fitting and regularized steady-state reconstruction (Sec. 9).

Figure 2: Forward simulation of the log-size jump transport model (CTMC discretization of Eq. ( 8 )) at representative times. The solution remains nonnegative and normalized by construction.

Experimental results

Research questions

• RQ1 How can fragmentation PBEs be written as a conservative transport problem in log-size space?
• RQ2 What is the exact jump-transport master equation governing normalized mass fractions, and how does it relate to drift–diffusion in the small-jump limit?
• RQ3 How can GKSL/Lindblad factorization be used as a constrained parameterization for inverse PSD inference?
• RQ4 What spectral/principal-PSD parametrizations arise under detailed balance, and how can they facilitate inverse problems?
• RQ5 What are practical, validated routes for forward modeling and inverse recovery under noise in synthetic benchmarks?

Key findings

• An exact conservative transport master equation in log-size space is derived for normalized mass fractions in fragmentation PBEs (Eq. 8).
• A GKSL/Lindblad factorization provides a structure-preserving, constrained parameterization for inverse modeling (Eq. 12–13).
• In the small-jump regime, the nonlocal jump transport reduces to a drift–diffusion operator in log-size space (Eq. 16).
• Under detailed balance, the operator admits a self-adjoint Schrödinger-type spectral form enabling compact PSD shape parametrizations (Eq. 17–20).
• Two practical inverse strategies are proposed: time-resolved parametric fitting and regularized steady-state inversion to recover effective potentials from PSDs (Sec. 9).
• Numerical validation includes forward CTMC-based jump transport, drift–diffusion reduction, and inverse parameter recovery on a synthetic Airy-half-line benchmark under multiplicative noise (Sec. 10).

Figure 3: Drift–diffusion reduction residual relative to the jump transport solution at the same times as Fig. 2 : $p_{\mathrm{FP}}(\xi,t)-p_{\mathrm{jump}}(\xi,t)$ . This highlights where the Fokker–Planck reduction deviates from the nonlocal jump transport model.