Numerical solution of breakage population balance equations using differential algebraic equations form

In this paper a new numerical approach is formulated to solve the breakage population balance equations along with its moments. In this approach the breakage equations are reformulated to a new class of differential equations known as differential algebraic equations. These equations involve differential equations along with algebraic constraints. Here we consider the moment preserving constraint as the algebraic constraint. The resulting differential algebraic equations are solved numerically using multistep methods of the MATLAB ODE suite.