A numerical model based on the finite volume scheme is proposed to approximate the binary breakage problems. Initially, it is considered that the particle fragments are characterized by a single property, i.e. particle’s volume. We then investigate the extension of the proposed model for solving breakage problems considering two properties of particles. The efficiency to estimate the different moments with good accuracy and simple extension for multi-variable problems are the key features of the proposed method. Moreover, the mathematical convergence analysis is performed for one-dimensional problems. All mathematical findings and numerical results are validated over several test problems. For numerical validation, we propose the extension of Bourgade & Filbet (2008
Math. Comput.
77
, 851–882. (
doi:10.1090/S0025-5718-07-02054-6
)) model for solving two-dimensional pure breakage problems. In this aspect, numerical treatment of the two-dimensional binary breakage models using finite volume methods can be treated to be the first instance in the literature.