A Conservative and Positivity-Preserving Discontinuous Galerkin Method for the Population Balance Equation

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Xu, Ziyao, Liu, Guanyang, Zhang, Yong-Tao
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909699790602240
author Xu, Ziyao
Liu, Guanyang
Zhang, Yong-Tao
author_facet Xu, Ziyao
Liu, Guanyang
Zhang, Yong-Tao
contents We develop a conservative, positivity-preserving discontinuous Galerkin (DG) method for the population balance equation (PBE), which models the distribution of particle numbers across particle sizes due to growth, nucleation, aggregation, and breakage. To ensure number conservation in growth and mass conservation in aggregation and breakage, we design a DG scheme that applies standard treatment for growth and nucleation, and introduces a novel discretization for aggregation and breakage. The birth and death terms are discretized in a symmetric double-integral form, evaluated using a common refinement of the integration domain and carefully selected quadrature rules. Beyond conservation, we focus on preserving the positivity of the number density in aggregation-breakage. Since local mass corresponds to the first moment, the classical Zhang-Shu limiter, which preserves the zeroth moment (cell average), is not directly applicable. We address this by proving the positivity of the first moment on each cell and constructing a moment-conserving limiter that enforces nonnegativity across the domain. To our knowledge, this is the first work to develop a positivity-preserving algorithm that conserves a prescribed moment. Numerical results verify the accuracy, conservation, and robustness of the proposed method.
format Preprint
id arxiv_https___arxiv_org_abs_2507_16631
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle A Conservative and Positivity-Preserving Discontinuous Galerkin Method for the Population Balance Equation
Xu, Ziyao
Liu, Guanyang
Zhang, Yong-Tao
Numerical Analysis
We develop a conservative, positivity-preserving discontinuous Galerkin (DG) method for the population balance equation (PBE), which models the distribution of particle numbers across particle sizes due to growth, nucleation, aggregation, and breakage. To ensure number conservation in growth and mass conservation in aggregation and breakage, we design a DG scheme that applies standard treatment for growth and nucleation, and introduces a novel discretization for aggregation and breakage. The birth and death terms are discretized in a symmetric double-integral form, evaluated using a common refinement of the integration domain and carefully selected quadrature rules. Beyond conservation, we focus on preserving the positivity of the number density in aggregation-breakage. Since local mass corresponds to the first moment, the classical Zhang-Shu limiter, which preserves the zeroth moment (cell average), is not directly applicable. We address this by proving the positivity of the first moment on each cell and constructing a moment-conserving limiter that enforces nonnegativity across the domain. To our knowledge, this is the first work to develop a positivity-preserving algorithm that conserves a prescribed moment. Numerical results verify the accuracy, conservation, and robustness of the proposed method.
title A Conservative and Positivity-Preserving Discontinuous Galerkin Method for the Population Balance Equation
topic Numerical Analysis
url https://arxiv.org/abs/2507.16631