Solution of the population balance equation using parallel adaptive cubature on GPUs

• Two GPU implementations of an adaptive cubature method were developed.
• One was parallelized at the integrand level and other at the integral formula level.
• Speedup was dependent on the number of integrands and their complexity.
• Maximum speedup was around 70 for 64k integrands.
• Three bivariate PBM were solved with DuQMoGeM with speedups from 11 to 15.

The Dual Quadrature Method of Generalized Moments (DuQMoGeM) is an accurate moment method for solving the population balance equation (PBE). The drawback of DuQMoGeM is the high computational cost associated with numerical integrations of the PBE integral terms in which each integrand can be integrated independently and, therefore, amenable to parallelization on GPUs. In this work, two parallel adaptive cubature algorithms were implemented on a hybrid architecture (CPU–GPU) to accelerate the DuQMoGeM. The speedup and scalability of these parallel algorithms were studied with different types of Genz's test functions. Then, we applied these parallel numerical integration algorithms in the DuQMoGeM solution of the PBE for three bivariate cases, obtaining speedups between 11 and 15.