×

Convergence analysis of finite volume scheme for nonlinear aggregation population balance equation. (English) Zbl 1416.65298

Summary: In this work, we introduce the convergence analysis of the recently developed finite volume scheme to solve a pure aggregation population balance equation that is of substantial interest in many areas such as chemical engineering, aerosol physics, astrophysics, polymer science, pharmaceutical sciences, and mathematical biology. The notion of the finite volume scheme is to conserve total mass of the particles in the system by introducing weight in the formulation. The consistency of the finite volume scheme is also analyzed thoroughly as it is an influential factor. The convergence study of the numerical scheme shows second order convergence on uniform, nonuniform smooth (geometric) as well as on locally uniform meshes independent of the aggregation kernel. Moreover, the first-order convergence is shown when the finite volume scheme is implemented on oscillatory and random meshes. In order to check the accuracy, the numerical experimental order of convergence is also computed for the physically relevant as well as analytically tractable kernels and validated against its analytical results.

MSC:

65M08 Finite volume methods for initial value and initial-boundary value problems involving PDEs
35R09 Integro-partial differential equations
65R20 Numerical methods for integral equations
Full Text: DOI

References:

[1] Terrazas‐VelardeK, PeglowM, TsotsasE. Investigation of the kinetics of fluidized bed spray agglomeration based on stochastic methods. AIChE J. 2011;57(11):3012‐3026.
[2] Terrazas‐VelardeK, PeglowM, TsotsasE. Kinetics of fluidized bed spray agglomeration for compact and porous particles. Chem Eng Sci. 2011;66(9):1866‐1878.
[3] HussainM, KumarJ, TsotsasE. A new framework for population balance modeling of spray fluidized bed agglomeration. Particuology. 2015;19:141‐154.
[4] KumarA, VercruysseJ, VanhoorneV, et al. Conceptual framework for model‐based analysis of residence time distribution in twin‐screw granulation. Eur J Pharm Sci. 2015;71:25‐34.
[5] IsmailHY, SinghM, DarwishS, et al. Developing ann‐kriging hybrid model based on process parameters for prediction of mean residence time distribution in twin‐screw wet granulation. Powder Technol. 2019;343:568‐577.
[6] KaurG, SinghM, MatsoukasT, KumarJ, De BeerT, NopensI. Two‐compartment modeling and dynamics of top‐sprayed fluidized bed granulator. Appl Math Model. 2019;68:267‐280.
[7] HoYK, KirseC, BriesenH, SinghM, ChanC‐H, KowK‐W. Towards improved predictions for the enzymatic chain‐end scission of natural polymers by population balances: the need for a non‐classical rate kernel. Chem Eng Sci. 2018;176:329‐342.
[8] RamkrishnaD. Population Balances, Theory and applications to particulate systems in engineering. New York: Academic Press; 2000.
[9] HulburtHM, KatzS. Some problems in particle technology: A statistical mechanical formulation. Chem Eng Sci. 1964;19:555‐574.
[10] ScottWT. Analytic studies of cloud droplet coalescence. J Atmos Sci. 1968;25:54‐65.
[11] AldousDJ. Deterministic and stochastic models for coalescence (aggregation and coagulation): a review of the mean‐field theory for probabilists. Bernoulli. 1999;5(1):3‐48. · Zbl 0930.60096
[12] RigopoulosS, JonesAG. Finite‐element scheme for solution of the dynamic population balance equation. AIChE J. 2003;49:1127‐1139.
[13] KimYP, SeinfeldJH. Simulation of multicomponent aerosol dynamics. J Colloid Interface Sci. 1992;149:425‐449.
[14] LaurenziIJ, BartelsJD, DiamondSL. A general algorithm for exact simulation of multicomponent aggregation processes. J Comput Phys. 2002;177:418‐449. · Zbl 1139.82331
[15] AlexopoulosA, RoussosA, KiparissidesC. Part V: dynamic evolution of the multivariate particle size distribution undergoing combined particle growth and aggregation. Chem Eng Sci. 2009;64(14):3260‐3269.
[16] ChakrabortyJ, KumarJ, SinghM, MahoneyA, RamkrishnaD. Inverse problems in population balances. Determination of aggregation kernel by weighted residuals. Ind Eng Chem Res. 2015;54(42):10530‐10538.
[17] StrumendoM, ArastoopourH. Solution of bivariate population balance equations using the finite size domain complete set of trial functions method of moments. Ind Eng Chem Res. 2008;48(1):262‐273.
[18] RezakhanlouF. Moment bounds for the solutions of the smoluchowski equation with coagulation and fragmentation. Proc R Soc Edinburgh Sec A Math. 2010;140(5):1041‐1059. · Zbl 1206.35055
[19] SinghM, KumarJ, BückA, TsotsasE. A volume‐consistent discrete formulation of aggregation population balance equations. Math Methods Appl Sci. 2015;39(9):2275‐2286. · Zbl 1338.35453
[20] LaurençotP, MischlerS. From the discrete to the continuous coagulation-fragmentation equations. Proc R Soc Edinburgh Sec A Math. 2002;132(5):1219‐1248. · Zbl 1034.35011
[21] FilbetF, LaurençotP. Numerical simulation of the smoluchowski coagulation equation. SIAM J Sci Comput. 2004;25(6):2004‐2028. · Zbl 1086.65123
[22] BourgadeJ‐P, FilbetF. Convergence of a finite volume scheme for coagulation‐fragmentation equations. Math Comput. 2008;77(262):851‐882. · Zbl 1140.82032
[23] Forestier‐CosteL, ManciniS. A finite volume preserving scheme on nonuniform meshes and for multidimensional coalescence. SIAM J Sci Comput. 2012;34(6):B840‐B860. · Zbl 1259.82054
[24] SinghM, KumarJ, BückA. A volume conserving discrete formulation of aggregation population balance equations on non‐uniform meshes. IFAC‐PapersOnLine. 2015;48(1):192‐197.
[25] KumarJ, KaurG, TsotsasE. An accurate and efficient discrete formulation of aggregation population balance equation. Kinetic Related Models. 2016;9(2):373‐391. · Zbl 1333.45008
[26] KumarS, RamkrishnaD. On the solution of population balance equations by discretization‐i. a fixed pivot technique. Chem Eng Sci. 1995;51(8):1311‐1332.
[27] KumarJ, WarneckeG. Convergence analysis of sectional methods for solving breakage population balance equations‐I: the fixed pivot technique. Numerische Mathematik. 2008;111(1):81‐108. · Zbl 1154.65381
[28] KumarJ, PeglowM, WarneckeG, HeinrichS. The cell average technique for solving multi‐dimensional aggregation population balance equations. Comput Chem Eng. 2008;32(8):1810‐1830.
[29] KumarJ, WarneckeG. Convergence analysis of sectional methods for solving breakage population balance equations‐II: the cell average technique. Numerische Mathematik. 2008;110(4):539‐559. · Zbl 1187.65146
[30] SinghM, KaurG, KumarJ, De BeerT, NopensI. A comparative study of numerical approximations for solving the smoluchowski coagulation equation. Braz J Chem Eng. 2018;35(4):1339‐1350.
[31] SahaJ, KumarJ. Development of a mass conserving discretization technique for breakage problems and its convergence analysis. Int J Adv Eng Sci Appl Math. 2015;7(1‐2):51‐61. · Zbl 1342.82128
[32] GiriAK, HausenblasE. Convergence analysis of sectional methods for solving aggregation population balance equations: the fixed pivot technique. Nonlinear Anal Real World Appl. 2013;14(6):2068‐2090. · Zbl 1323.92007
[33] GiriAK, NagarAK. Convergence of the cell average technique for smoluchowski coagulation equation. ESAIM Math Model Numer Anal. 2015;49(2):349‐372. · Zbl 1315.65109
[34] GiriAK. Mathematical and numerical analysis for coagulation‐fragmentation equations. PhD Thesis; 2011.
[35] KumarR, KumarJ, WarneckeG. Convergence analysis of a finite volume scheme for solving non‐linear aggregation‐breakage population balance equations. Kinetic Related Models. 2014;7(4):713‐737. · Zbl 1316.65116
[36] KaurG, SinghM, KumarJ, De BeerT, NopensI. Mathematical modeling and simulation of spray fluidized bed granulator. Processes. 2018;6(10):195.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.