Abstract: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 provides 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 10000 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. (C) 2014 Elsevier Inc. All rights reserved.

Optimized Parallelization of Boundary Integral Poisson-Boltzmann Solvers

Optimized Parallelization of Boundary Integral Poisson-Boltzmann Solvers

Accurate, robust and reliable calculations of Poisson-Boltzmann solvation energies

Parallel Higher-order Boundary Integral Electrostatics Computation on Molecular Surfaces with Curved Triangulation

Coupling finite and boundary element methods to solve the Poisson--Boltzmann equation for electrostatics in molecular solvation

A combined steepest descent and genetic algorithm (SD/GA) approach for the optimization of solvation parameters

Poisson-Boltzmann model for protein-surface electrostatic interactions and grid-convergence study using the PyGBe code

Parallel Adaptive Finite Element Algorithms for Solving the Coupled Electro-diffusion Equations

A GPU Poisson-Fermi Solver for Ion Channel Simulations

Quantum Chemistry for Solvated Molecules on Graphical Processing Units (GPUs)using Polarizable Continuum Models

A Massively Parallel Performance Portable Free-space Spectral Poisson Solver

A Direct Simulation Approach for the Poisson-Boltzmann Equation Using the Random Batch Method.

Fast Iterative Method for Local Steric Poisson-Boltzmann Theories in Biomolecular Solvation

Accurate, Robust, And Reliable Calculations Of Poisson-Boltzmann Binding Energies

Automating GPU Scalability for Complex Scientific Models: Phonon Boltzman Transport Equation

Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques

Employing Graphics Processing Unit Technology, Alternating Direction Implicit Method And Domain Decomposition To Speed Up The Numerical Diffusion Solver For The Biomedical Engineering Research

Heuristic Tree-Partition-Based Parallel Method for Biophysically Detailed Neuron Simulation

Accelerating Population Balance-Monte Carlo Simulation for Coagulation Dynamics from the Markov Jump Model, Stochastic Algorithm and GPU Parallel Computing

Poisson-Boltzmann based machine learning (PBML) model for electrostatic analysis

GPU_PBTE: an efficient solver for three and four phonon scattering rates on graphics processing units