skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A moment projection method for population balance dynamics with a shrinkage term

Abstract

A new method of moments for solving the population balance equation is developed and presented. The moment projection method (MPM) is numerically simple and easy to implement and attempts to address the challenge of particle shrinkage due to processes such as oxidation, evaporation or dissolution. It directly solves the moment transport equation for the moments and tracks the number of the smallest particles using the algorithm by Blumstein and Wheeler (1973) . The performance of the new method is measured against the method of moments (MOM) and the hybrid method of moments (HMOM). The results suggest that MPM performs much better than MOM and HMOM where shrinkage is dominant. The new method predicts mean quantities which are almost as accurate as a high-precision stochastic method calculated using the established direct simulation algorithm (DSA).

Authors:
 [1]; ; ;  [2];  [3];  [1];  [2];  [4]
  1. Department of Mechanical Engineering, National University of Singapore, Engineering Block EA, Engineering Drive 1, 117576 (Singapore)
  2. Department of Chemical Engineering and Biotechnology, University of Cambridge, New Museums Site, Pembroke Street, Cambridge, CB2 3RA (United Kingdom)
  3. School of Chemical and Biomedical Engineering, Nanyang Technological University, 62 Nanyang Drive, 637459 (Singapore)
  4. (Singapore)
Publication Date:
OSTI Identifier:
22622252
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 330; Other Information: Copyright (c) 2016 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ACCURACY; ALGORITHMS; BALANCES; COMPUTERIZED SIMULATION; DISSOLUTION; EQUATIONS; EVAPORATION; MOMENTS METHOD; OXIDATION; PARTICLE TRACKS; PARTICLES; PERFORMANCE; SHRINKAGE; STOCHASTIC PROCESSES; TRANSPORT THEORY

Citation Formats

Wu, Shaohua, Yapp, Edward K.Y., Akroyd, Jethro, Mosbach, Sebastian, Xu, Rong, Yang, Wenming, Kraft, Markus, E-mail: mk306@cam.ac.uk, and School of Chemical and Biomedical Engineering, Nanyang Technological University, 62 Nanyang Drive, 637459. A moment projection method for population balance dynamics with a shrinkage term. United States: N. p., 2017. Web. doi:10.1016/J.JCP.2016.10.030.
Wu, Shaohua, Yapp, Edward K.Y., Akroyd, Jethro, Mosbach, Sebastian, Xu, Rong, Yang, Wenming, Kraft, Markus, E-mail: mk306@cam.ac.uk, & School of Chemical and Biomedical Engineering, Nanyang Technological University, 62 Nanyang Drive, 637459. A moment projection method for population balance dynamics with a shrinkage term. United States. doi:10.1016/J.JCP.2016.10.030.
Wu, Shaohua, Yapp, Edward K.Y., Akroyd, Jethro, Mosbach, Sebastian, Xu, Rong, Yang, Wenming, Kraft, Markus, E-mail: mk306@cam.ac.uk, and School of Chemical and Biomedical Engineering, Nanyang Technological University, 62 Nanyang Drive, 637459. Wed . "A moment projection method for population balance dynamics with a shrinkage term". United States. doi:10.1016/J.JCP.2016.10.030.
@article{osti_22622252,
title = {A moment projection method for population balance dynamics with a shrinkage term},
author = {Wu, Shaohua and Yapp, Edward K.Y. and Akroyd, Jethro and Mosbach, Sebastian and Xu, Rong and Yang, Wenming and Kraft, Markus, E-mail: mk306@cam.ac.uk and School of Chemical and Biomedical Engineering, Nanyang Technological University, 62 Nanyang Drive, 637459},
abstractNote = {A new method of moments for solving the population balance equation is developed and presented. The moment projection method (MPM) is numerically simple and easy to implement and attempts to address the challenge of particle shrinkage due to processes such as oxidation, evaporation or dissolution. It directly solves the moment transport equation for the moments and tracks the number of the smallest particles using the algorithm by Blumstein and Wheeler (1973) . The performance of the new method is measured against the method of moments (MOM) and the hybrid method of moments (HMOM). The results suggest that MPM performs much better than MOM and HMOM where shrinkage is dominant. The new method predicts mean quantities which are almost as accurate as a high-precision stochastic method calculated using the established direct simulation algorithm (DSA).},
doi = {10.1016/J.JCP.2016.10.030},
journal = {Journal of Computational Physics},
number = ,
volume = 330,
place = {United States},
year = {Wed Feb 01 00:00:00 EST 2017},
month = {Wed Feb 01 00:00:00 EST 2017}
}
  • This paper proposes a comprehensive framework for accelerating population balance-Monte Carlo (PBMC) simulation of particle coagulation dynamics. By combining Markov jump model, weighted majorant kernel and GPU (graphics processing unit) parallel computing, a significant gain in computational efficiency is achieved. The Markov jump model constructs a coagulation-rule matrix of differentially-weighted simulation particles, so as to capture the time evolution of particle size distribution with low statistical noise over the full size range and as far as possible to reduce the number of time loopings. Here three coagulation rules are highlighted and it is found that constructing appropriate coagulation rule providesmore » a route to attain the compromise between accuracy and cost of PBMC methods. Further, in order to avoid double looping over all simulation particles when considering the two-particle events (typically, particle coagulation), the weighted majorant kernel is introduced to estimate the maximum coagulation rates being used for acceptance–rejection processes by single-looping over all particles, and meanwhile the mean time-step of coagulation event is estimated by summing the coagulation kernels of rejected and accepted particle pairs. The computational load of these fast differentially-weighted PBMC simulations (based on the Markov jump model) is reduced greatly to be proportional to the number of simulation particles in a zero-dimensional system (single cell). Finally, for a spatially inhomogeneous multi-dimensional (multi-cell) simulation, the proposed fast PBMC is performed in each cell, and multiple cells are parallel processed by multi-cores on a GPU that can implement the massively threaded data-parallel tasks to obtain remarkable speedup ratio (comparing with CPU computation, the speedup ratio of GPU parallel computing is as high as 200 in a case of 100 cells with 10 000 simulation particles per cell). These accelerating approaches of PBMC are demonstrated in a physically realistic Brownian coagulation case. The computational accuracy is validated with benchmark solution of discrete-sectional method. The simulation results show that the comprehensive approach can attain very favorable improvement in cost without sacrificing computational accuracy.« less
  • The method of moments is a simple but efficient method of solving the population balance equation which describes particle dynamics. Recently, the moment projection method (MPM) was proposed and validated for particle inception, coagulation, growth and, more importantly, shrinkage; here the method is extended to include the fragmentation process. The performance of MPM is tested for 13 different test cases for different fragmentation kernels, fragment distribution functions and initial conditions. Comparisons are made with the quadrature method of moments (QMOM), hybrid method of moments (HMOM) and a high-precision stochastic solution calculated using the established direct simulation algorithm (DSA) and advantagesmore » of MPM are drawn.« less
  • A mathematical model has been constructed and verified to simulate dynamics of a microbial community in typical tundra. The model contains the following state variables: the population densities of three competing microbial species (exemplified by Arthrobacter, Pseudomonas, and Bacillus), indexes of their physiological state, concentration of available organic substrate, plant litter reserves, the amount of microbiovorous protozoans, and temperature. The mathematical model simulates adequately the qualitative features of microbial seasonal dynamics observed in tundra. The global warming and associated increase in primary productivity, as predicted by simulation, will relieve the pressure of L-selection and thus result in stabilization of themore » tundra microbial community. The model also predicts that aerobic decomposition of dead organic matter in solid will be accelerated compared to its formation. 24 refs., 7 figs., 1 tab.« less