[论文解读] 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 Studies被引用 0

一句话总结

本文推导了 fragmentation PBE 的在对数尺度上的严格守恒主方程，展示了小跃迁极限如何得到漂移-扩散，给出用于反演建模的可选 GKSL/Lindblad 因式分解，并在数值上验证前向与反演路径。

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.

研究动机与目标

• 激发对 fragmentation PBE 的兴趣，以及在对数尺度上实现守恒输运形式的必要性。
• 推导在对数尺度上以归一化质量分数表示的严格守恒主方程。
• 将跳跃输运与漂移-扩散极限联系起来，并探索谱形/PSD 参数化。
• 提供用于时间分辨拟合和稳态 PSD 重构的反演建模路径。
• 通过前向仿真和在合成基准上的反演恢复来验证方法。

提出的方法

• 在对数尺度上建立 fragmentation PBE 的形式，并推导严格单向跳跃输运主方程（Eq. 8）。
• 给出归一化质量分数密度 p(ξ,t) 与跳跃核 K(u)，以及速率 λ(ξ)（Eq. 9–10）。
• 给出一个可选的 GKSL/Lindblad 因式分解，用于在受约束的反演建模中再现主方程对角项（Eq. 12–13）。
• 在小跃迁区域，得到对数尺度上的漂移-扩散（Fokker–Planck）简化（Eq. 16）。
• 讨论基于详细平衡的谱形式和 Schrödinger 型化简，以用于 PSD 参数化（Eq. 17–20）。
• 概述两条反演路线：时分辨的参数拟合与正则化的稳态重建（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.

实验结果

研究问题

• RQ1 碎裂 PBEs 如何在对数尺度上被写成守恒输运问题？
• RQ2 支配归一化质量分数的严格跳跃输运主方程是什么，且它如何在小跃迁极限下与漂移-扩散相关？
• RQ3 如何将 GKSL/Lindblad 因式分解用作受约束的反PSD推断参数化？
• RQ4 在详细平衡下，出现哪些谱/主 PSD 参数化形式，以及它们如何促进反问题？
• RQ5 在合成基准下，包含噪声的前向建模和反演恢复有哪些实际且经过验证的路径？

主要发现

• 对于 fragmentation PBE 中归一化质量分数的对数尺度，已经推导出一个严格守恒输运主方程（Eq. 8）。
• GKSL/Lindblad 因式分解提供了一种结构保持、受约束的反演建模参数化（Eq. 12–13）。
• 在小跃迁极限，非局部的跳跃输运简化为对数尺度上的漂移-扩散算子（Eq. 16）。
• 在详细平衡下，该算子具备自伴的 Schrödinger 型谱形式，便于紧凑的 PSD 形状参数化（Eq. 17–20）。
• 提出两种实际的反演策略：时分辨的参数拟合与正则化稳态反演，以从 PSDs 重构有效势（Sec. 9）。
• 数值验证包括基于前向 CTMC 的跳跃输运、漂移-扩散简化，以及在乘性噪声下对一个合成的 Airy 半线基准的反演参数恢复（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.