[论文解读] Numerical simulation of solutions and moments of the smoluchowski coagulation equation

Dustin D. Keck, David M. Bortz|arXiv (Cornell University)|Dec 27, 2013

Advanced Mathematical Modeling in Engineering被引用 4

一句话总结

本文比较了有限元法（FEM）与有限体积法（FLFM）在数值求解Smoluchowski凝聚方程中的表现，重点关注解矩的近似精度与计算成本。研究发现，对于凝聚较慢的系统，FEM在显著更低的计算成本下提供了更优的精度；而对于快速凝聚，FLFM在零阶矩上略为精确，且在第一阶矩上始终更精确，两种方法均表现出二阶收敛性。

ABSTRACT

Researchers have employed variations of the Smoluchowski coagulation equation to model a wide variety of both organic and inorganic phenomena and with relatively few known analytical solutions, numerical solutions play an important role in studying this equation. In this article, we consider numerical approximations, focusing on how different discretization schemes impact the accuracy of approximate solution moments. Pursuing the eventual goal of comparing simulated solutions to experimental data, we must carefully choose the numerical method most appropriate to the type of data we attain. Within this context, we compare and contrast the accuracy and computational cost of a finite element approach and a finite volume-based scheme. Our study provides theoretical and numerical evidence that the finite element approach achieves much more accuracy when the system aggregates slowly, and it does so with much less computation cost. Conversely, the finite volume method is slightly more accurate approximating the zeroth moment when the system aggregates quickly and is much more accurate approximating the first moment in general. Lastly, our study also provides numerical evidence that the finite element method (conventionally considered first order) actually belongs to a class of discontinuous Galerkin methods that exhibit superconvergence, or second order in our case.

研究动机与目标

• 比较有限元法（FEM）与有限体积法（FLFM）在求解Smoluchowski凝聚方程时的精度与计算成本。
• 评估不同离散化格式对解矩（特别是零阶与第一阶矩）近似的影晌，尤其在不同凝聚动力学条件下的表现。
• 研究截断参数 $x_{\text{max}}$ 与网格间距 $\Delta x$ 对两种方法数值精度的影响。
• 根据实验数据类型（如库尔特计数器数据，对应零阶矩；或流式细胞术数据，对应第一阶矩）指导方法选择。
• 提供理论与数值证据，支持FEM的收敛速率与超收敛行为。

提出的方法

• 本研究采用基于Banks与Kappel前期工作的有限元法（FEM），并由Ackleh与Fitzpatrick进一步扩展，用于凝聚方程的空间离散化。
• 采用一种有限体积型格式——Filbet与Laurentc̆ot通量法（FLFM）作为对比基准，该方法在 $L^1$ 范数下已知具有二阶收敛性。
• 对两种凝聚核进行数值求解：常数核 $K_A(x,y) \equiv 1$ 与乘法核 $K_A(x,y) = xy$。
• 通过与精细网格参考解及已知矩的比较，评估FEM与FLFM的精度。
• 通过改变 $\Delta x$ 与 $x_{\text{max}}$ 进行网格细化研究，以评估对敏感性与收敛行为的影响。
• 通过 $M_i(f; x_1,x_2) = \int_{x_1}^{x_2} x^i f(t,x) \, dx$ 计算矩，其中 $f$ 表示粒径分布，$g = x f$ 表示体积分布。

实验结果

研究问题

• RQ1 在近似Smoluchowski凝聚方程的完整解时，FEM与FLFM在精度上如何比较？
• RQ2 在慢速与快速凝聚动力学下，哪种方法对零阶与第一阶矩的近似更准确？
• RQ3 截断参数 $x_{\text{max}}$ 对每种方法的精度有何影响，尤其是在实验数据仅限于小体积范围时？
• RQ4 尽管传统上认为FEM为一阶方法，但其是否如数值证据所示表现出超收敛？
• RQ5 在给定精度水平下，哪种方法具有更高的计算效率？

主要发现

• FEM在 $L^1[\mathbf{X}]$ 范数下表现出二阶收敛性，支持其具有超收敛的猜想，这与其传统上被归类为一阶方法的观念相悖。
• 对于凝聚较慢的系统（$K_A(x,y) \equiv 1$），FEM在与FLFM相近的解精度下，计算成本显著更低。
• 在慢速凝聚中，FEM对零阶矩的近似精度高于FLFM，而FLFM在第一阶矩上略为精确。
• 在快速凝聚（$K_A(x,y) = xy$）情况下，FLFM在所有情形下对第一阶矩的近似更精确，且在零阶矩上也略为优越。
• FEM对 $x_{\text{max}}$ 的选择敏感度低于FLFM，尤其在乘法核下表现更优，因此在实验数据上限检测范围有限时更具鲁棒性。
• 在同等精度下，FEM在800个网格点上的计算成本远低于FLFM在100个网格点上的成本，展现出显著的计算效率优势。